- 1Yunnan Tobacco Company Dali State Branch, Dali, Yunnan, China
- 2College of Agronomy and Biotechnology, Southwest University/Engineering Research Center of South Upland Agriculture, Ministry of Education, Chongqing, China
Background: Intercropping tobacco with soybean is a sustainable approach to improve resource use efficiency and crop resilience. However, the optimal soybean planting density for maximizing the physiological and protective benefits to tobacco remains unclear.
Methods: A field experiment was conducted in Yunnan Province, China, including five treatments: tobacco monoculture and four tobacco–soybean intercropping densities. Photosynthetic parameters, carbon and nitrogen metabolism, defense-related physiology, and leaf chemical composition were measured across key growth stages.
Results: Intercropping density significantly affected photosynthetic and metabolic performance in both species. The medium density with four soybean holes achieved the best results, increasing the net photosynthetic rate of tobacco by 30.8% compared with monoculture during the vigorous growth stage. This treatment also enhanced PSII photochemical efficiency, with Fv/Fm and ΦPSII values both higher than other treatments, and chlorophyll a content increased by 32.9% compared with high-density intercropping. The activities of Rubisco and nitrate reductase rose by 18.8% and 49.2%, respectively. At the same time, this density reduced the incidence of tobacco black shank disease and increased salicylic acid and jasmonic acid contents by 38.9% and 33.7%. Peroxidase, superoxide dismutase, and phenylalanine ammonia-lyase activities were also elevated. Tobacco leaves under this treatment showed a balanced chemical composition with high sugar, high potassium, and low chlorine contents, resulting in superior flue-cured quality and the highest economic return.
Conclusion: The four-hole soybean density optimized photosynthesis, nitrogen metabolism, and defense responses, improving tobacco quality and yield. These findings provide a physiological and agronomic basis for developing efficient and sustainable tobacco–soybean intercropping systems.
1 Introduction
Tobacco (Nicotiana tabacum L.) is one of the most important economic crops in China, and its yield and quality are closely related to the sustainable development of the tobacco industry (Zhu et al., 2024). However, long-term monocropping has led to soil nutrient imbalance, frequent soil-borne diseases, and a decline in leaf quality, all of which have restricted the high-quality development of the tobacco sector (Wang et al., 2024). In recent years, intercropping and relay cropping systems have gained widespread attention as green and efficient cultivation practices, owing to their ability to improve soil conditions, enhance crop stress resistance, and increase overall productivity (Akchaya et al., 2025). Within tobacco production systems, how to enhance photosynthetic efficiency and optimize yield and quality through rational intercropping design has become a major research focus.
Soybean (Glycine max L.) - a typical leguminous crop - possesses the unique advantage of biological nitrogen fixation and can improve nitrogen availability in the soil (Kebede, 2021).Moreover, the differences in growth duration and ecological niches between soybean and tobacco promote a more efficient utilization of light and nutrient resources, and help to alleviate the obstacles caused by continuous tobacco cropping (Zhang et al., 2025). Recent studies have shown that tobacco-legume intercropping can improve soil fertility, promote nutrient accumulation in tobacco, and enhance economic returns. Tobacco - soybean intercropping significantly increased soil nutrient contents, with total nitrogen, available phosphorus, and available potassium all rising by more than 15% compared with monoculture (Gu et al., 2025). In addition, intercropping with soybean increased the populations of nitrifying and denitrifying bacteria in the tobacco rhizosphere, which promoted nitrogen cycling and improved the sugar-to-alkaloid and potassium-to-chlorine ratios in tobacco leaves, ultimately enhancing leaf quality (Tu, 2015).
Nevertheless, current research on tobacco-soybean intercropping has mainly focused on soil fertility improvement and nutrient interactions, whereas studies on the photosynthetic physiology and carbon–nitrogen metabolic regulation of tobacco under different intercropping conditions remain limited (Gu et al., 2025). The spatial configuration of an intercropping system directly influences canopy light distribution and energy balance. Inappropriate planting densities can intensify light competition, reduce light interception by tobacco, and suppress photosynthetic activity and carbon assimilation, ultimately disrupting the source–sink balance and lowering yield and quality (Wang et al., 2023). In contrast, a suitable intercropping density improves light capture and utilization, maintains the activity of photosynthetic organs, and promotes coordination between carbon and nitrogen metabolism, thereby enhancing the conversion and assimilation efficiency of photosynthetic products (Zhou et al., 2023a). Therefore, optimizing intercropping density is crucial for developing an efficient tobacco–soybean system. Understanding the photosynthetic responses of tobacco under different density configurations is essential for revealing interspecific light utilization mechanisms and improving resource-use efficiency and sustainability in tobacco production.
In addition, soil-borne diseases such as tobacco black shank frequently occur in continuous tobacco monocropping systems (Chen et al., 2023). Previous studies have indicated that intercropping can mitigate disease incidence by improving microclimate conditions and modifying rhizosphere microbial community structures (Boudreau, 2013; Yang et al., 2014). Moreover, intercropping can induce plant defense responses, including activation of the salicylic acid (SA) signaling pathway, enhancement of antioxidant enzyme activities, and stimulation of phenolic metabolism, thereby improving crop resistance (Barna et al., 2003; Gullner et al., 2017; Zheng et al., 2022b). However, the physiological and biochemical mechanisms underlying disease resistance regulation in tobacco under different intercropping densities remain poorly understood.
Yunnan Province is one of the major flue-cured tobacco-producing regions in China, and tobacco–soybean intercropping has already been adopted by local farmers (Gu et al., 2025). Nevertheless, due to the lack of systematic studies on intercropping density and spatial configuration, inappropriate soybean densities–either too high or too low–are commonly observed in production, leading to insufficient light interception by tobacco, decreased photosynthetic capacity, and reduced leaf quality. Previous studies have demonstrated that intercropping can enhance canopy photosynthetic efficiency and stress tolerance by optimizing light distribution, promoting carbon–nitrogen coordination, and activating defense-related physiological processes (Yao et al., 2017; Harman et al., 2021).
Therefore, in this study, the main flue-cured tobacco cultivar ‘Hongda’ commonly grown in Yunnan was used to establish a tobacco-soybean intercropping system with varying soybean densities. This study systematically analyzed the changes in photosynthetic physiology, carbon-nitrogen metabolism, and defense-related physiological traits of tobacco under different intercropping densities, aiming to elucidate the synergistic mechanisms of light energy utilization and disease resistance. The findings provide a theoretical basis for optimizing the design of tobacco-soybean intercropping systems, improving leaf quality, and promoting efficient resource utilization in Yunnan’s tobacco production.
2 Materials and methods
2.1 Experimental site
The field experiment was conducted from April to August 2025 in Midu County, Dali Bai Autonomous Prefecture, Yunnan Province, China (25.38°N, 100.41°E).The experimental area is characterized by a mid-subtropical monsoon climate, with an average annual temperature of 17.2 °C and an average annual precipitation of 642.5 mm. The monthly rainfall during the experimental period is shown in Figure 1, and the physicochemical properties of the experimental soil are presented in Table 1.
2.2 Experimental design
The field experiment was conducted on flat land with relatively uniform soil fertility using a randomized complete block design (RCBD) with five treatments: tobacco monoculture (S0) and four tobacco–soybean intercropping densities (S2, S4, S6, and S8). Soybean was sown 15 days after tobacco transplanting, with two plants per hill. Intercropping density was adjusted by varying the number of soybean hills: two hills (S2, 6.6 × 104 plants ha-1), four hills (S4, 1.32 × 105 plants ha-1), six hills (S6, 1.98 × 105 plants ha-1), and eight hills (S8, 2.64 × 105 plants ha-1). Each treatment was replicated five times.
Tobacco seedlings were transplanted on April 26, 2025, and soybean was sown 15 days later. Soybean was planted on the ridge surface, 20 cm away from the tobacco stem base on both sides. In treatments S2, S4, S6, and S8, one, two, three, and four soybean hills were planted on each side of the tobacco plant, respectively. Each plot covered an area of 42.00 m², with a row spacing of 100 cm and a plant spacing of 50 cm. A 1 m-wide buffer zone was set between adjacent plots.
The experimental field had been continuously cultivated with flue-cured tobacco for three consecutive years prior to this study. After tobacco harvest in August of the previous year, no other crops were planted. Before transplanting, the soil was tilled and ridged following standard local tobacco field management practices. Other agronomic measures, including irrigation, fertilization, and pest control, were conducted in accordance with local production guidelines to ensure uniform management across all treatments and replicates. The main flue-cured tobacco cultivar ‘Hongda’ and the local soybean cultivar ‘Yunhuang No. 13’ were used in this study. These cultivars are widely planted in the Dali tobacco-growing region and represent typical varieties for local production systems.
2.3 Photosynthetic characteristics
Photosynthetic parameters of both tobacco and soybean were measured simultaneously at corresponding growth stages: the rosette stage (RS) of tobacco and the flowering stage (R2) of soybean; the fast-growing stage (FGS) of tobacco and the pod-filling stage (R4) of soybean; and the maturity stage (MS) of tobacco and the grain-filling stage (R6) of soybean. Measurements were conducted on clear days between 09:00 and 11:30 a.m. using a portable photosynthesis system (LI-6400XT, LI-COR, USA). The net photosynthetic rate (Pn), stomatal conductance (Gs), intercellular CO2 concentration (Ci), and transpiration rate (Tr) were determined on the third fully expanded leaf from the top of tobacco and the uppermost fully expanded trifoliolate leaf of soybean.
During measurement, the photosynthetic photon flux density (PPFD) was set at 1200 μmol·m-2·s-1. The reference CO2 concentration in the leaf chamber was maintained at 400 μmol·mol-1 using a CO2 cylinder, and the airflow rate was fixed at 500 μmol·s-1. The leaf chamber temperature was set at 25 °C, and relative humidity was maintained between 60% and 70% (Ding et al., 2025). Prior to measurement, leaves were illuminated for approximately 10 minutes to reach a steady photosynthetic state. Each leaf was measured at three different positions, and the average value was calculated. Five plants were measured per treatment.
On the same day as photosynthetic measurements, leaves were dark-adapted for 30 minutes, and the maximal photochemical efficiency of photosystem II (Fv/Fm) was determined using a JUNIOR-PAM chlorophyll fluorometer (Heinz Walz, Germany). Under light-adapted conditions (identical to those used for photosynthetic measurements), the actual photochemical quantum yield (ΦPSII) was recorded.
For pigment analysis, leaf discs were sampled from the same leaves used for photosynthetic gas exchange measurements. Pigments were extracted with 80% acetone, and the concentrations of chlorophyll a, chlorophyll b, and carotenoids were determined spectrophotometrically according to the method of Lichtenthaler and Wellburn (1983).
2.4 Activities of key enzymes involved in carbon and nitrogen metabolism
Leaves previously used for photosynthetic measurements were sampled to determine the activities of key enzymes involved in carbon and nitrogen metabolism. Ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco) activity was measured using the NADH-coupled spectrophotometric method at 340 nm (Li et al., 2025). Nitrate reductase (NR) activity was determined by the Griess colorimetric method at 540 nm (Imran et al., 2019). Glutamine synthetase (GS) activity was assayed using the hydroxylamine colorimetric method at 540 nm (Peng et al., 2016).
2.5 Disease assessment
The incidence of tobacco black shank disease (Phytophthora nicotianae) was investigated every 15 days starting 30 days after transplanting and continued until 90 days. In each plot, disease occurrence was assessed using the five-point sampling method on 20 randomly selected plants. Disease severity was recorded according to the Grading and Investigation Method for Tobacco Diseases and Insect Pests (GB/T 23222–2008). The disease incidence and disease index were calculated as follows:
Based on the disease index obtained from five consecutive surveys, the area under the disease progress curve (AUDPC) was calculated using the following formula:
i = the order of the survey (from the first to the (n-1)th observation); = the disease index recorded at the i and surveys, respectively; and = the number of days after transplanting (DAT) at the i and surveys; n = the total number of surveys.
2.6 Determination of defense-related physiological parameters
Defense-related physiological parameters of tobacco were measured during the fast-growing stage. The concentrations of salicylic acid (SA) and jasmonic acid (JA) were quantified using enzyme-linked immunosorbent assay (ELISA) kits (Shanghai Youxuan Biotechnology Co., Ltd., China) following the manufacturer’s protocol. Peroxidase (POD) activity was determined by the guaiacol oxidation method (Li and Gong, 2008). Superoxide dismutase (SOD) activity was assayed using the nitroblue tetrazolium (NBT) photochemical reduction method (Qu et al., 2014). Phenylalanine ammonia-lyase (PAL) activity was measured spectrophotometrically at 290 nm using L-phenylalanine as the substrate, and expressed as the rate of trans-cinnamic acid formation (Jin et al., 2009). Total polyphenol (TP) content was determined by the Folin–Ciocalteu colorimetric method (Cheng et al., 2022).
2.7 Tobacco leaf quality and yield assessment
After flue-curing, 2 kg of middle tobacco leaves (C3F grade) were collected from each treatment to determine the chemical composition. Total and reducing sugar contents were measured according to YC/T 159–2002 (State Tobacco Monopoly Administration, 2002a), total nitrogen according to YC/T 161–2002 (State Tobacco Monopoly Administration, 2002b), total alkaloids according to YC/T 468–2013 (State Tobacco Monopoly Administration, 2013), potassium content according to YC/T 217–2007 (State Tobacco Monopoly Administration, 2007), and chlorine content according to YC/T 162–2011 (State Tobacco Monopoly Administration, 2011).
In early September 2025, the flue-cured tobacco leaves from each treatment were graded according to GB 2635–1992. The yield, output value, proportion of high- and middle-grade leaves, and average price were recorded. Tobacco prices were based on the local purchase prices set by the Tobacco Company, while soybean income was calculated according to market purchase prices and measured yields.
2.8 Data processing and statistical analysis
Experimental data were initially organized and calculated using Microsoft Excel 2019. Statistical analyses were performed using SPSS 26.0 software. Two-way analysis of variance (ANOVA) was applied to evaluate the effects of intercropping density, growth stage, and their interactions. For data involving a single factor, one-way ANOVA was used. Significant differences among treatments were determined by Duncan’s multiple range test at P < 0.05. All figures and visualizations were prepared using Origin 2022 software.
3 Results
3.1 Photosynthetic parameters of tobacco and soybean
The gas exchange parameters of tobacco are presented in Figures 2a–d. During the RS stage, Pn values were generally low across treatments, but S4 and S6 were significantly higher than S0 and S8 (Figure 2a). At the FGS stage, Pn increased markedly and reached its maximum across the growth period, with S4 and S6 attaining 26.41 and 25.80 μmol CO2·m-2·s-1, respectively, both significantly higher than the other treatments. At the MS stage, Pn decreased overall, yet S4 and S6 maintained higher values, being 30.77% and 25.00% greater than S0, respectively. The variation of Gs followed a pattern similar to Pn (Figure 2b). During the RS stage, S2, S4, and S6 were significantly higher than S0 and S8. At the FGS stage, Gs further increased, reaching 0.56, 0.57, and 0.58 mol·m-2·s-1 in S2, S4, and S6, respectively, all significantly higher than S0 and S8. During the MS stage, Gs declined, but S4 remained significantly higher than S0, S2, and S8. Ci at the RS stage was 287.5 and 291 μmol·mol-1 for S0 and S8, respectively, both significantly higher than S2, S4, and S6 (Figure 2c). At the FGS stage, Ci in S8 rose to 335.2 μmol·mol-1, significantly higher than in S4 and S6, while no significant differences were found among treatments at the MS stage. Tr at the RS stage was significantly higher in S2, S4, and S6 than in S0 and S8 (Figure 2d). During the FGS stage, Tr increased, with S4 showing the highest value of 4.60 mmol·m-2·s-1, 9.52% higher than S0, while S8 showed the lowest value of 4.01 mmol·m-2·s-1. At the MS stage, Tr in S8 remained significantly lower than the other treatments.
Figure 2. Photosynthetic parameters of tobacco and soybean under different intercropping densities. (a) net photosynthetic rate of tobacco; (b) stomatal conductance of tobacco; (c) intercellular CO2 concentration of tobacco; (d) transpiration rate of tobacco; (e) net photosynthetic rate of soybean; (f) stomatal conductance of soybean; (g) intercellular CO2 concentration of soybean; (h) transpiration rate of soybean. RS, rosette stage of tobacco; FGS, fast-growing stage of tobacco; MS, maturity stage of tobacco. Different lowercase letters indicate significant differences at P < 0.05. *, **, and *** represent significance at P < 0.05, P < 0.01, and P < 0.001, respectively.
The gas exchange parameters of soybean are shown in Figures 2e–h. At the R2 stage, Pn was highest in S4 and S6, reaching 31.51 and 30.82 μmol CO2·m-2·s-1, respectively, both significantly higher than S2 and S8 (Figure 2e). At the R4 stage, Pn increased slightly, with S4 reaching 33.04 μmol CO2·m-2·s-1, significantly higher than S2 and S8. At the R6 stage, Pn decreased overall, but S4 and S6 remained significantly higher than S8 by 37.5% and 31.25%, respectively. Gs at the R2 stage was highest in S4, reaching 0.35 mol·m-2·s-1, significantly higher than S2 and S8 (Figure 2f), while S6 did not differ significantly from S4. At the R4 stage, Gs in S4 and S6 remained significantly higher than in S2 and S8, and though Gs declined at the R6 stage, it remained highest in S4. For Ci, S8 exhibited the highest value of 286.3 μmol·mol-1 at the R2 stage, significantly higher than other treatments (Figure 2g). Ci in S4 and S6 was lower, at 261.3 and 268.4 μmol·mol-1, respectively, showing the same pattern at the R4 stage, with S8 remaining significantly highest. At the R6 stage, Ci in S8 increased further to 335.2 μmol·mol-1, significantly higher than in S4 and S6. At the R2 stage, Tr was highest in S4 and S6, reaching 9.40 and 9.21 mmol·m-2·s-1, respectively, both significantly higher than S2 and S8 (Figure 2h). At the R4 stage, Tr further increased, with S4 and S6 being 25.81% and 22.58% higher than S8, respectively. During the R6 stage, Tr decreased, but S4 and S6 still maintained significantly higher values than S2 and S8.
In summary, S4 and S6 consistently maintained higher Pn, Gs, and Tr values, along with lower Ci, throughout the growth period of both tobacco and soybean. This indicates that appropriate intercropping density effectively enhanced photosynthetic carbon assimilation and gas exchange coordination in the system.
3.2 Chlorophyll fluorescence parameters of tobacco and soybean
Different intercropping densities, growth stages, and their interactions had significant effects on the chlorophyll fluorescence parameters of both tobacco and soybean (Tables 2, 3). The chlorophyll fluorescence parameters of tobacco are presented in Table 2. During the RS stage, Fv/Fm values ranged from 0.77 to 0.80, with S4 significantly higher than S0 and S8. After entering the FGS stage, Fv/Fm slightly decreased across treatments, with S4 at 0.74, still significantly higher than the others. By the MS stage, Fv/Fm declined further, but S4 maintained the highest value at 0.69, representing an 11.29% increase compared with S8. ΦPSII showed that S4 and S6 were significantly higher than S0, S2, and S8 during all three stages (RS, FGS, and MS). Across all stages, ΦPSII first increased and then decreased, reaching its maximum at the FGS stage, where S4 and S6 were 19.44% and 16.67% higher than S0, respectively. The chlorophyll fluorescence parameters of soybean are shown in Table 3. At the R2 and R4 stages, Fv/Fm in S4 was significantly higher than in S2 and S8, and at the R6 stage, S4 still maintained the highest value, significantly exceeding S8. ΦPSII in S4 was 7.69% higher than in S8 at the R2 stage and increased slightly at the R4 stage. During the R4 and R6 stages, S4 remained significantly higher than S2 and S8. Overall, S4 and S6 treatments consistently maintained higher Fv/Fm and ΦPSII values throughout the growth period of both tobacco and soybean, indicating enhanced photochemical efficiency and more stable PSII function under moderate intercropping density.
Table 2. Chlorophyll fluorescence parameters of tobacco under different soybean intercropping densities.
Table 3. Chlorophyll fluorescence parameters of soybean under different soybean intercropping densities.
3.3 Photosynthetic pigment contents of tobacco and soybean
The pigment contents of tobacco leaves are shown in Figures 3a–c. Chlorophyll a levels were generally higher during the RS and FGS stages, with S4 and S6 significantly exceeding S8. At the FGS stage, chlorophyll a in S4 was 32.94% higher than in S8. By the MS stage, pigment contents decreased across all treatments, but chlorophyll a in S4 remained 1.83 and 2.06 times higher than in S0 and S8, respectively (Figure 3a). Chlorophyll b showed a similar trend to chlorophyll a during the RS and FGS stages, with S4 and S6 significantly higher than S0 and S8. At the FGS stage, chlorophyll b in S4 was 55.56% higher than in S8, followed by a pronounced decline at the MS stage (Figure 3b). Carotenoid content was highest in S4 during the RS stage (0.24 mg·g-1), significantly higher than S0, S6, and S8. Although carotenoid content declined during the FGS and MS stages, S4 consistently remained the highest among treatments (Figure 3c).
Figure 3. Pigment contents of tobacco and soybean at different growth stages under various intercropping densities. (a) chlorophyll a content of tobacco; (b) chlorophyll b content of tobacco; (c) carotenoid content of tobacco; (d) chlorophyll a content of soybean; (e) chlorophyll b content of soybean; (f) carotenoid content of soybean. Different lowercase letters indicate significant differences at P < 0.05. *, **, and *** represent significance at P < 0.05, P < 0.01, and P < 0.001, respectively.
The pigment contents of soybean leaves are presented in Figures 3d–f. At the R2 stage, chlorophyll a in S2, S4, and S6 was significantly higher than in S8, increasing by 12.14%, 15.71%, and 14.29%, respectively. At the R4 stage, only S2 was significantly higher than S8. By the R6 stage, S2 and S4 reached the highest levels, at 1.46 and 1.49 mg·g-1, respectively (Figure 3d). Chlorophyll b showed a similar trend across the R2, R4, and R6 stages, with S2, S4, and S6 significantly higher than S8. At the R2 stage, these treatments were 30.95%, 40.48%, and 35.71% higher than S8, respectively (Figure 3e). Carotenoid content in S4 and S6 reached 0.32 and 0.33 mg·g-1 at the R2 stage, both significantly higher than S8. During the R4 stage, carotenoid levels increased significantly across all treatments, with S6 reaching the highest value of 0.59 mg·g-1. At the R6 stage, carotenoid content declined markedly, but S4 and S6 still remained significantly higher than S2 and S8 (Figure 3f).
Overall, S4 and S6 treatments maintained higher levels of chlorophylls and carotenoids throughout the growth period of both tobacco and soybean, providing a stable pigment foundation for sustained and efficient photosynthetic performance.
3.4 Activities of key enzymes involved in carbon and nitrogen metabolism in tobacco
The activities of Rubisco, nitrate reductase (NR), and glutamine synthetase (GS) were significantly influenced by intercropping density and growth stage (Figures 4a–c). Overall, the three enzymes exhibited a consistent trend of “promotion under moderate density and inhibition under excessive density.” During the RS stage, Rubisco activity in S4 and S6 was significantly higher than in S0 and S8. At the FGS stage, Rubisco activity reached its maximum, with S4 and S6 showing 18.75% and 16.25% increases over S0, respectively. By the MS stage, Rubisco activity declined in all treatments (Figure 4a). NR activity also varied significantly across growth stages, with S4 and S6 maintaining the highest levels during both the FGS and MS stages. At the FGS stage, NR activity in S4 was 49.17% higher than in S0 (Figure 4b). GS activity followed a similar pattern to Rubisco and NR, peaking at the FGS stage. S4 and S6 were significantly higher than S0, S2, and S8. Although GS activity decreased during the MS stage, S4 and S6 still maintained relatively higher levels (Figure 4c). Overall, S4 and S6 treatments substantially enhanced the activities of Rubisco, NR, and GS, thereby promoting carbon assimilation and nitrogen metabolism in tobacco under optimal intercropping density.
Figure 4. Activities of key enzymes involved in carbon and nitrogen metabolism in tobacco under different intercropping densities. (a) ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco) activity; (b) nitrate reductase (NR) activity; (c) glutamine synthetase (GS) activity. Different lowercase letters indicate significant differences at P < 0.05. *, **, and *** represent significance at P < 0.05, P < 0.01, and P < 0.001, respectively.
3.5 Investigation and evaluation of tobacco black shank disease
The incidence rate and disease index of tobacco black shank increased significantly with the number of days after transplanting, but clear differences were observed among intercropping densities (Table 4). Throughout the investigation period, S4 exhibited the lowest disease occurrence, with both incidence rate and disease index remaining significantly lower than those of other treatments. At 90 days after transplanting, the incidence rate and disease index of S4 were 19.0% and 27.3, respectively, representing reductions of 47.22% and 47.30% compared with monocropped tobacco. Disease severity in S2 and S6 was also significantly reduced, although still higher than in S4. In contrast, S8 showed much greater susceptibility to black shank infection.
Table 4. Incidence rate and disease index of tobacco black shank under different soybean intercropping densities.
After integrating the disease index across the entire infection period, the area under the disease progress curve (AUDPC) differed significantly among treatments (Figure 5). The AUDPC value of S4 was the lowest, significantly smaller than that of all other treatments. No significant difference was found between S2 and S6, but both recorded markedly lower values than the monocropping treatment. In contrast, S8 exhibited the highest AUDPC value. Overall, S2, S4, and S6 effectively suppressed the incidence and progression of tobacco black shank disease, whereas S8 intensified disease development.
Figure 5. Area under the disease progress curve (AUDPC) for tobacco black shank under different intercropping densities. Different lowercase letters indicate significant differences at P < 0.05. *, **, and *** represent significance at P < 0.05, P < 0.01, and P < 0.001, respectively.
3.6 Defense-related enzyme activities and signal molecule contents in tobacco
Significant differences were observed among treatments in the defense-related physiological parameters (Figure 6). The contents of salicylic acid (SA) and jasmonic acid (JA) were significantly influenced by intercropping density, showing a trend of “enhancement under moderate density and reduction under excessive density” (Figures 6a, b). Compared with monocropping, S2, S4, and S6 markedly increased the SA and JA contents in leaves, with the most pronounced effects in S4 and S6. The SA contents in S4 and S6 were 38.93% and 33.67% higher than those in S0, respectively. However, SA and JA levels in S8 were significantly lower than in monocropping.
Figure 6. Defense-related enzyme activities and signal molecule contents in tobacco leaves under different intercropping densities. (a) salicylic acid (SA) content; (b) jasmonic acid (JA) content; (c) total polyphenol (TP) content; (d) peroxidase (POD) activity; (e) superoxide dismutase (SOD) activity; (f) phenylalanine ammonia-lyase (PAL) activity. Different lowercase letters indicate significant differences at P < 0.05.
Total polyphenol (TP) content showed a pattern similar to that of SA and JA, exhibiting a significant increase under the S6 treatment, which was 21.43% higher than S0 (Figure 6c). Meanwhile, peroxidase (POD) and superoxide dismutase (SOD) activities were significantly higher in S4 and S6 than in S0, whereas POD activity in S8 was significantly lower than in S0 (Figures 6d, e). Phenylalanine ammonia-lyase (PAL) activity was also highest in S4 and S6, being 33.33% and 36.36% higher than in monocropping, respectively. S2 showed a moderate increase, while S8 exhibited a marked decline, consistent with the trend observed for TP content (Figure 6f).
In summary, S4 and S6 treatments substantially enhanced SA, JA, and TP contents as well as POD, SOD, and PAL activities in tobacco, thereby strengthening both systemic and non-specific defense responses. In contrast, excessive intercropping density (S8) suppressed these physiological processes, weakening the overall defense capacity of tobacco.
3.7 Chemical quality and economic benefits of flue-cured tobacco
Significant differences were observed in the chemical composition of flue-cured tobacco leaves among treatments (Table 5). For sugar content, variations in total and reducing sugars followed a similar pattern. S4 and S6 exhibited the highest total sugar contents, reaching 26.84% and 26.91%, respectively, both significantly higher than S0 and S2. The reducing sugar content in S4 reached 17.26%, significantly higher than in all other treatments, whereas S8 showed no significant difference compared with S0 and S2. Regarding nitrogenous compounds, total alkaloid content did not differ significantly among treatments. However, total nitrogen content was significantly higher in S4, S6, and S8 compared with S0. Potassium content was highest in S4 (3.02%), significantly higher than in S0 and S8, while chlorine content was greatest in S0 and S8 and lowest in S4. Overall, S4 displayed the most favorable balance among key chemical quality indicators, characterized by high sugar, high potassium, and moderate nitrogen–alkaloid levels, which correspond to superior flue-cured tobacco quality.
Table 5. Chemical quality of flue-cured tobacco leaves under different soybean intercropping densities.
Different intercropping densities had significant effects on the economic performance of the tobacco–soybean system (Table 6). In terms of tobacco yield, there were no significant differences among S0, S2, S4, and S6, whereas S8 showed a marked reduction. The proportion of high-grade leaves was highest in S0 and lowest in S8. The trend in tobacco output value followed a pattern similar to yield, with S0 and S4 significantly higher than S8. In contrast, soybean output value showed the opposite trend: S8 achieved the highest soybean income, significantly exceeding the other intercropping treatments, while S2 was the lowest. From the perspective of total system productivity, S4 and S6 achieved the highest overall economic returns, with total output values significantly higher than S0. Although S8 had the lowest tobacco output, its combined system value did not differ significantly from the other treatments. These results indicate that S4 and S6 provided the optimal balance between tobacco yield and soybean profitability, resulting in the greatest overall system benefit.
Table 6. Yield and economic benefits of the tobacco–soybean intercropping system under different soybean intercropping densities.
4 Discussion
4.1 Effects of soybean intercropping density on the photosynthetic characteristics of tobacco and soybean
This study demonstrated that soybean intercropping density significantly affected the photosynthetic characteristics of tobacco and simultaneously exerted a feedback influence on the photosynthetic performance of soybean itself. Overall, tobacco intercropped with four or six soybean hills exhibited significantly higher photosynthetic rates, chlorophyll fluorescence parameters, and pigment contents in both crops. In contrast, excessive density with eight soybean hills caused photosynthetic inhibition and restricted metabolism, indicating that an appropriate intercropping density can enhance light energy utilization efficiency by regulating canopy structure and resource allocation (Li et al., 2024b).
During the rosette stage of tobacco, the overall photosynthetic rate was low. However, in the treatments with four or six intercropped soybean hills, the net photosynthetic rate (Pn), stomatal conductance (Gs), and transpiration rate (Tr) were significantly higher than in monocropping, while the intercellular CO2 concentration (Ci) was lower, suggesting that photosynthetic limitation at this stage was mainly stomatal. Moderate intercropping improved inter-row ventilation and light conditions, promoting CO2 diffusion and stomatal regulation, thereby enhancing carbon assimilation efficiency (Fan et al., 2018; Dong et al., 2024). During the fast-growing stage, Pn increased markedly, with S4 and S6 treatments showing higher photosynthetic capacity, whereas the S8 treatment exhibited a sharp decline due to excessive canopy shading and insufficient light. At the maturity stage, Pn decreased overall, but moderate intercropping still maintained relatively high levels. In contrast, high-density intercropping showed decreased Gs but increased Ci, indicating that the limiting factor of photosynthesis shifted from stomatal to non-stomatal (Sakoda et al., 2021).
The photosynthetic responses of soybean followed a similar trend to those of tobacco. At the R2 and R4 stages, under the tobacco–soybean (four-hole) intercropping system, Pn, Gs, and Tr in soybean leaves were significantly higher than those under the eight-hole treatment, suggesting that moderate intercropping improved the light energy utilization and CO2 assimilation efficiency of soybean. Because the tobacco canopy exhibited relatively good light transmittance, moderate intercropping provided soybean plants with a more uniform light environment, helping to alleviate photoinhibition and delay leaf senescence (Wang et al., 2024). Under high-density intercropping, however, the lower soybean leaves were continuously exposed to low-light conditions, resulting in decreased photosynthetic rate and elevated intercellular CO2 concentration, indicating restricted carbon assimilation.
Chlorophyll fluorescence parameters reflect the activity of photosystem II (PSII) and energy conversion efficiency (Semer et al., 2019). In this study, tobacco intercropped with four and six soybean holes maintained higher maximum photochemical efficiency (Fv/Fm) and actual photochemical quantum yield (ΦPSII) at all growth stages, suggesting stronger PSII electron transport activity and more efficient light energy conversion (Li et al., 2009). Moderate intercropping can optimize the light distribution within the canopy, where proper shading reduces excessive light stress and the lower leaves receive more scattered light, thereby improving overall light-use efficiency (Chen et al., 2024). For soybean, the S4 treatment exhibited significantly higher Fv/Fm and ΦPSII at the R2, R4, and R6 stages compared with the eight-hole treatment, as the tobacco canopy served as a light buffer, mitigating photodamage and demonstrating a complementary and mutual light benefit effect (Zhang et al., 2021). Thus, a light environment characterized by both moderate shading and adequate transmittance can simultaneously enhance the photosynthetic efficiency of both main and intercrop species in intercropping systems (Brooker et al., 2015).
This study also found that moderate intercropping significantly increased the contents of chlorophyll a, chlorophyll b, and carotenoids in both tobacco and soybean leaves, especially in the four-hole intercropping treatment, where pigment levels remained high during the middle and late growth stages. The maintenance of pigment content not only enhances light-harvesting capacity but also helps delay chloroplast degradation, thereby slowing leaf senescence (Griffiths et al., 2014). Carotenoids play an important role in the antioxidant system by dissipating excess excitation energy through non-photochemical quenching (NPQ) and reducing the accumulation of reactive oxygen species (ROS), thus maintaining the structural stability of PSII (Latowski et al., 2011). In the eight-hole intercropping treatment, however, pigment contents in both tobacco and soybean decreased significantly, likely due to excessive shading and intensified nutrient competition, which accelerated chloroplast degradation and the aging of photosynthetic apparatus (Niinemets, 2023). In summary, a moderate intercropping density is essential for maintaining pigment homeostasis and sustaining efficient photosynthetic function in tobacco and soybean.
4.2 Effects of soybean intercropping density on carbon and nitrogen metabolism in tobacco
Rubisco, nitrate reductase (NR), and glutamine synthetase (GS) are key enzymes that connect photosynthetic carbon fixation with nitrogen assimilation. In the treatments with four and six intercropped soybean hills, the activities of these three enzymes were significantly higher than in monocropped tobacco, indicating that appropriate intercropping density enhanced the carbon and nitrogen metabolism of tobacco (Cong et al., 2015). The increase in Rubisco activity not only represents an enhancement of carboxylation capacity but also reflects a higher level of enzyme protein synthesis and activation (Iñiguez et al., 2021). This suggests that tobacco under moderate intercropping has a stronger capacity for CO2 fixation and energy supply, providing sufficient ATP and NADPH for subsequent nitrogen assimilation reactions (Kramer and Evans, 2011).
Meanwhile, the elevated activities of NR and GS indicate that nitrogen uptake, reduction, and assimilation processes were simultaneously enhanced, resulting in greater accumulation of nitrogen metabolites such as glutamine and glutamate in tobacco leaves. These compounds serve as essential substrates for amino acid and protein synthesis (Akhtar et al., 2024). Because tobacco leaves under moderate intercropping exhibit higher photosynthetic rates, the increased carbon skeletons produced can serve as substrates for NR and GS catalysis (Yoneyama and Suzuki, 2019). Moreover, the enhanced photosynthetic electron transport ensures sufficient NADPH supply from the light reactions, which strengthens NR reduction activity (Yamori et al., 2015). In addition to the stimulation of nitrogen assimilation by improved carbon metabolism, the availability and transformation efficiency of nitrogen itself were also enhanced by intercropping.
Soybean biological nitrogen fixation improves soil nitrogen availability. Through the secretion of organic acids, amino compounds, and flavonoids, soybean roots stimulate rhizosphere microbial activity, increasing the concentrations of ammonium nitrogen and soluble organic nitrogen in the tobacco root zone. These effects, in turn, promote the synthesis and activation of NR and GS (Lin et al., 2024; Molefe et al., 2023). Photosynthetic carbon metabolites also provide carbon skeletons for nitrogen assimilation, while enhanced nitrogen assimilation further facilitates the synthesis of key photosynthetic components such as Rubisco and photosystem proteins (Paul and Pellny, 2003). Under these positive feedback conditions, both photosynthetic productivity and nitrogen metabolism were simultaneously enhanced.
However, under the eight-hole intercropping treatment, excessive canopy density restricted light penetration and intensified belowground competition. Consequently, the synthesis and activation of Rubisco proteins declined, resulting in suppressed CO2 carboxylation rates (Zheng et al., 2022a). Moreover, the excessive root density of soybean increased soil oxygen consumption and led to ammonium accumulation, where high concentrations of NH4+ directly inhibited NR activity. In addition, hypoxic conditions in the rhizosphere impaired root metabolic function and induced rhizosphere acidification, further reducing NR and GS activities (Churchland and Grayston, 2014). These factors disrupted the coupling between carbon and nitrogen metabolism, causing metabolic imbalance. Therefore, a moderate intercropping density not only achieves complementarity in light energy and nutrient use at the canopy level but also maintains coordination between carbon–nitrogen flux and energy utilization at the cellular level. This coordination represents the fundamental physiological basis for the efficient operation of composite agricultural systems.
One limitation of this study is that only enzyme activities were analyzed, while gene expression and metabolite profiling were not included. Integrating transcriptomic and metabolomic analyses in future work would help clarify the molecular mechanisms underlying the coordination of carbon and nitrogen metabolism under different intercropping densities.
4.3 Effects of soybean intercropping density on the disease incidence
This study showed that in the tobacco–soybean intercropping system, moderate soybean intercropping effectively suppressed disease occurrence, whereas excessively high intercropping density produced the opposite effect. This indicates that there exists an optimal density range within the intercropping system that balances interspecific competition and mutualistic interactions.
Crop population density directly determines the transmission opportunities of pathogens. Moderate intercropping can reduce the infection probability by introducing non-host crops, thereby diluting the host density of the pathogen (Douma and Noordhoek, 2025). However, when the overall density becomes too high, the distance between plants decreases, effectively increasing the contact frequency among hosts and accelerating the spread of disease (Luo et al., 2021). Intercropping also influences disease development by altering the field microclimate. The inclusion of soybean modifies canopy ventilation and light penetration, which affects spore dispersal pathways and their survival in the air (Boudreau, 2013). Similar canopy-mediated effects have been reported in potato strip-planting systems, where intercropping significantly reduced the severity of late blight (Homulle et al., 2025). Conversely, excessive planting density leads to elevated canopy humidity and restricted airflow, providing favorable conditions for pathogen infection (Singh et al., 2023). This “density threshold effect” is consistent with findings by Wu et al. (2013) and Phillips (2005) on sheath blight in rice and late blight in potato.
Non-host crops can also indirectly suppress disease through rhizosphere microecological regulation. Isoflavonoid compounds secreted by soybean roots possess antifungal activity and can induce systemic resistance in tobacco (Yang et al., 2021; Sugiyama, 2019). In addition, root exudates from intercrop plants can reshape the rhizosphere microbial community, forming a disease-suppressive microecological network (Zhou et al., 2023b). Antagonistic microbial groups such as Actinobacteria and Bacillus species, enriched in the soybean rhizosphere, can inhibit soil-borne pathogens through competitive exclusion and secretion of antimicrobial metabolites (Lee et al., 2021; Li et al., 2024a). However, when soybean density is excessively high, intensified root competition may disrupt microbial community stability, potentially creating a favorable environment for pathogens (Yang et al., 2014).
In conclusion, moderate tobacco–soybean intercropping density suppresses tobacco black shank disease through a combination of host dilution, canopy regulation, and rhizosphere microecological optimization. This provides a theoretical basis for developing green and sustainable disease management strategies in tobacco cultivation systems.
4.4 Effects of soybean intercropping density on defense-related enzyme activities and signaling molecules in tobacco
The results of this study demonstrated that tobacco–soybean intercropping regulated the defense-related signaling molecules and resistance enzyme activities in tobacco, indicating that the system induced systemic resistance at the physiological and molecular levels. In the treatments with four or six intercropped soybean hills, the contents of salicylic acid (SA) and jasmonic acid (JA) increased significantly, accompanied by enhanced activities of phenylalanine ammonia-lyase (PAL) and antioxidant enzymes (POD and SOD). These findings suggest that intercropping enhanced the defensive potential of tobacco through the synergistic activation of the SA and JA signaling pathways (Qiu et al., 2014; Asghari and Hasanlooe, 2015).
SA and JA constitute two core signaling pathways of the plant immune system. SA primarily mediates systemic acquired resistance (SAR), whereas JA participates in induced systemic resistance (ISR) and defense against necrotrophic pathogens and herbivores (Ali et al., 2018; Heil and Bostock, 2002). Under moderate intercropping density, the simultaneous increase in SA and JA levels indicates that tobacco defense may not rely on a single pathway but rather on the coordinated activation of both, enabling simultaneous responses from SAR and ISR. Previous studies have shown that the SAR pathway induces the expression of pathogenesis-related (PR) proteins such as PR-1 and β-1,3-glucanase, while the ISR pathway upregulates transcription factors such as WRKY and MYC2 via JA and ethylene signaling, thereby activating the antioxidant system (Darby et al., 2000; Javed and Gao, 2023). In this study, the significant enhancement of PAL, POD, and SOD activities was consistent with the changes in SA and JA levels, indicating that intercropping density may have strengthened systemic defense by activating these two signaling cascades. Root exudates from soybean, including flavonoids and phenolic acids, may serve as exogenous signaling molecules that activate the ISR pathway in tobacco roots (Mandal et al., 2010).
At the leaf level, SA accumulation enhances PAL pathway activity, promoting the phenylpropanoid metabolism toward the synthesis of lignin and phenolic compounds, thereby reinforcing structural resistance (Dixon et al., 2002). JA, on the other hand, enhances reactive oxygen signaling and induces the expression of antioxidant enzymes (SOD and POD), contributing to non-specific defense enhancement (Qiu et al., 2014). Therefore, under moderate intercropping density, positive crosstalk between SA and JA signaling likely results in synergistic enhancement of both systemic and local defense responses. In contrast, in the high-density treatment with eight intercropped soybean hills, the contents of SA and JA decreased significantly, and the activities of PAL and POD were suppressed, indicating that excessive density impaired the tobacco defense system.
In summary, the tobacco–soybean intercropping system activates the SA and JA signaling pathways, thereby inducing the coordinated expression of key enzymes in phenylpropanoid metabolism (such as PAL), PR protein-related defense genes, and antioxidant enzymes (POD and SOD). Together, these responses establish an integrated and efficient systemic defense network in tobacco.
4.5 Effects of soybean intercropping density on system yield and economic benefits
The results of this study showed that soybean intercropping density significantly affected the yield and economic benefits of the tobacco–soybean composite system. The treatment with four intercropped soybean hills per tobacco plant achieved the highest economic output per unit land area by increasing soybean income while maintaining stable tobacco quality. In this treatment, both tobacco and soybean contributed to the highest overall system benefit, as the increase in soybean yield complemented the stable economic return from tobacco.
The balance of chemical quality in flue-cured tobacco depends on the coordinated proportions of key components such as sugars, nitrogen, potassium, and chlorine. The treatment with four intercropped soybean hills significantly increased the total sugar and reducing sugar contents of tobacco leaves while maintaining low chlorine and appropriate nitrogen–alkaloid levels. This formed a “high-sugar, high-potassium, low-chlorine” compositional pattern, characteristic of high-quality flue-cured tobacco (Chen et al., 2021). The increase in sugar improves aroma and combustibility, whereas potassium accumulation buffers leaf acidity and enhances burning quality (Jiang et al., 2024; Banožić et al., 2020). These results are consistent with previous findings that intercropping can improve tobacco leaf quality (Gu et al., 2025), likely due to the enhancement of rhizosphere nutrient cycling and regulation of carbon–nitrogen metabolism by soybean roots (Nelson and Sadowsky, 2015; Liu et al., 2018). However, intercropping eight soybean hills per tobacco plant led to decreased sugar and potassium contents and increased chlorine content, likely because excessive planting density restricted photosynthesis and intensified nutrient competition, thereby reducing the accumulation and transformation efficiency of photosynthetic products in the leaves (Chen et al., 2019; Carrier et al., 2019).
The dynamic balance between competition and complementarity among crops is crucial for determining yield and economic efficiency in intercropping systems (Li et al., 2014). In this study, although monocropped tobacco had a higher individual yield value, the additional soybean income in the intercropping system significantly improved land use efficiency. Moderate intercropping achieved a “stable tobacco with enhanced soybean” effect, while excessive density reduced total system benefit due to inhibited tobacco growth and quality decline. Overall, intercropping tobacco with four soybean hills optimized both tobacco chemical composition and economic benefit, achieving the best balance between quality and efficiency. This configuration provides a theoretical basis for developing a green and complementary cropping system characterized by “tobacco supporting soybean and soybean enriching tobacco.”
5 Conclusion
This study demonstrated that an appropriate soybean intercropping density can markedly improve the photosynthetic performance of tobacco while simultaneously enhancing defense capacity and overall system productivity. Moderate intercropping significantly increased the net photosynthetic rate, chlorophyll content, and photosystem II efficiency of tobacco, and promoted carbon and nitrogen metabolism by elevating the activities of key enzymes such as Rubisco, nitrate reductase, and glutamine synthetase. Meanwhile, higher levels of salicylic acid (SA) and jasmonic acid (JA) collectively strengthened tobacco resistance and alleviated black shank disease. Among all treatments, intercropping four soybean hills per tobacco plant showed the best comprehensive performance. This configuration optimized canopy light distribution and resource utilization, improved photosynthetic efficiency and tobacco chemical quality, and formed a desirable “high-sugar, high-potassium, low-chlorine” composition. Consequently, it achieved optimal coordination among tobacco yield, quality, and system economic benefit. In summary, moderate soybean intercropping enhances photosynthetic productivity and physiological resilience in tobacco while improving the sustainability of the cropping system. Future studies should focus on multi-year and multi-regional validation of the optimized model and explore its underlying physiological and ecological mechanisms to support wider application.
Data availability statement
The original contributions presented in the study are included in the article/supplementary material. Further inquiries can be directed to the corresponding author.
Author contributions
XL: Writing – original draft, Conceptualization, Data curation, Methodology, Visualization, Validation, Investigation. YaH: Validation, Methodology, Data curation, Formal Analysis, Software, Writing – original draft, Investigation. CY: Resources, Writing – original draft, Methodology, Validation, Investigation. JL: Writing – original draft, Methodology, Data curation, Software, Investigation. CL: Writing – original draft, Formal Analysis, Data curation, Software, Validation. NT: Supervision, Investigation, Project administration, Writing – original draft. HZ: Writing – original draft, Project administration, Supervision, Investigation. SJ: Methodology, Investigation, Writing – review & editing, Resources. JS: Resources, Investigation, Writing – review & editing, Validation. DW: Methodology, Validation, Supervision, Writing – review & editing. CX: Methodology, Writing – review & editing, Validation, Software. YuH: Writing – review & editing, Methodology, Investigation. ML: Writing – review & editing, Conceptualization, Supervision, Funding acquisition, Project administration.
Funding
The author(s) declare financial support was received for the research and/or publication of this article. This work was supported by the Key Scientific and Technological Projects of Yunnan Tobacco Monopoly Administration (Project No. 2024530000241012); and the Science and Technology Plan Project of Dali Prefecture Branch of Yunnan Tobacco Company (DLYC2023001, DLYC2025001, DLYC2025002). Yunnan Tobacco Company was not involved in the study design, collection, analysis, interpretation of data, the writing of this article, or the decision to submit it for publication.
Conflict of interest
All authors were employed by Yunnan Tobacco Company Dali State Branch.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Akchaya, K., Parasuraman, P., Pandian, K., Vijayakumar, S., Thirukumaran, K., Mustaffa, M. R. A. F., et al. (2025). Boosting resource use efficiency, soil fertility, food security, ecosystem services, and climate resilience with legume intercropping: A review. Front. Sustain. Food Syst. 9. doi: 10.3389/fsufs.2025.1527256
Akhtar, K., Ain, N. U., Prasad, P. V., Naz, M., Aslam, M. M., Djalovic, I., et al. (2024). Physiological, molecular, and environmental insights into plant nitrogen uptake and metabolism under abiotic stresses. Plant Genome 17, e20461. doi: 10.1002/tpg2.20461
Ali, A., Shah, L., Rahman, S., Riaz, M. W., Yahya, M., Xu, Y. J., et al. (2018). Plant defense mechanism and current understanding of salicylic acid and NPRs in activating SAR. Physiol. Mol. Plant Pathol. 104, 15–22. doi: 10.1016/j.pmpp.2018.09.002
Asghari, M. and Hasanlooe, A. R. (2015). Interaction effects of salicylic acid and methyl jasmonate on total antioxidant content, catalase and peroxidase enzyme activity in “Sabrosa” strawberry fruit during storage. Sci. Hortic. 197, 490–495. doi: 10.1016/j.scienta.2015.10.034
Banožić, M., Jokić, S., Ačkar, Đ., Blažić, M., and Šubarić, D. (2020). Carbohydrates–key players in tobacco aroma formation and quality determination. Molecules 25, 1734. doi: 10.3390/molecules25071734
Barna, B., Fodor, J., Pogány, M., and Király, Z. (2003). Role of reactive oxygen species and antioxidants in plant disease resistance. Pest Manage. Sci. 59, 459–464. doi: 10.1002/ps.682
Boudreau, M. A. (2013). Diseases in intercropping systems. Annu. Rev. Phytopathol. 51, 499–519. doi: 10.1146/annurev-phyto-082712-102246
Brooker, R. W., Bennett, A. E., Cong, W. F., Daniell, T. J., George, T. S., Hallett, P. D., et al. (2015). Improving intercropping: a synthesis of research in agronomy, plant physiology and ecology. New Phytol. 206, 107–117. doi: 10.1111/nph.13132
Carrier, M., Gonzalez, F. A. R., Cogliastro, A., Olivier, A., Vanasse, A., and Rivest, D. (2019). Light availability, weed cover and crop yields in second generation of temperate tree-based intercropping systems. Field Crops Res. 239, 30–37. doi: 10.1016/j.fcr.2019.05.011
Chen, J., Li, Y., He, X., Jiao, F., Xu, M., Hu, B., et al. (2021). Influences of different curing methods on chemical compositions in different types of tobaccos. Ind. Crops Prod 167, 113534. doi: 10.1016/j.indcrop.2021.113534
Chen, G., Liu, M., Zhao, X., Bawa, G., Liang, B., Feng, L., et al. (2024). Improved photosynthetic performance under unilateral weak light conditions in a wide–narrow-row intercropping system is associated with altered sugar transport. J. Exp. Bot. 75, 258–273. doi: 10.1093/jxb/erad379
Chen, Y., Ren, K., He, X., Gong, J., Hu, X., Su, J., et al. (2019). Dynamic changes in physiological and biochemical properties of flue-cured tobacco of different leaf ages during flue-curing and their effects on yield and quality. BMC Plant Biol. 19, 555. doi: 10.1186/s12870-019-2189-8
Chen, Y., Yang, L., Zhang, L., Li, J., Zheng, Y., Yang, W., et al. (2023). Autotoxins in continuous tobacco cropping soils and their management. Front. Plant Sci. 14. doi: 10.3389/fpls.2023.1106033
Cheng, J., Wang, F., Li, X., Du, Y., and Xie, Y. (2022). Progress in the study of phenolic allelochemicals: species, extraction, separation, and detection. Jiangsu Agric. Sci. 50, 8–15. doi: 10.15889/j.issn.1002-1302.2022.06.002
Churchland, C. and Grayston, S. J. (2014). Specificity of plant–microbe interactions in the tree mycorrhizosphere biome and consequences for soil C cycling. Front. Microbiol. 5. doi: 10.3389/fmicb.2014.00261
Cong, W. F., Hoffland, E., Li, L., Six, J., Sun, J. H., Bao, X. G., et al. (2015). Intercropping enhances soil carbon and nitrogen. Glob. Change Biol. 21, 1715–1726. doi: 10.1111/gcb.12738
Darby, R. M., Maddison, A., Mur, L. A., Bi, Y. M., and Draper, J. (2000). Cell-specific expression of salicylate hydroxylase in an attempt to separate localized HR and systemic signaling establishing SAR in tobacco. Mol. Plant Pathol. 1, 115–123. doi: 10.1046/j.1364-3703.2000.00014.x
Ding, S., Cheng, T., Wang, B., Yu, D., Rao, D., Meng, F., et al. (2025). Effects of dense planting on canopy photosynthetic production and yield formation of soybean cultivars released in different eras. Acta Agron. Sin. 51, 161–173. doi: 10.3724/SP.J.1006.2025.00161
Dixon, R. A., Achnine, L., Kota, P., Liu, C. J., Reddy, M. S., and Wang, L. (2002). The phenylpropanoid pathway and plant defence–a genomics perspective. Mol. Plant Pathol. 3, 371–390. doi: 10.1046/j.1364-3703.2002.00131.x
Dong, B., Wang, Z., Evers, J. B., Stomph, T. J., Van Der Putten, P. E., Yin, X., et al. (2024). Competition for light and nitrogen with an earlier-sown species negatively affects leaf traits and leaf photosynthetic capacity of maize in relay intercropping. Eur. J. Agron. 155, 127119. doi: 10.1016/j.eja.2024.127119
Douma, J. C. and Noordhoek, R. (2025). The effect of plant host density on disease incidence–A meta-analysis. Plant Cell Environ. 48, 6632–6644. doi: 10.1111/pce.XXXXXX
Fan, Y., Chen, J., Cheng, Y., Raza, M. A., Wu, X., Wang, Z., et al. (2018). Effect of shading and light recovery on the growth, leaf structure, and photosynthetic performance of soybean in a maize–soybean relay-strip intercropping system. PloS One 13, e0198159. doi: 10.1371/journal.pone.0198159
Griffiths, C. A., Gaff, D. F., and Neale, A. D. (2014). Drying without senescence in resurrection plants. Front. Plant Sci. 5. doi: 10.3389/fpls.2014.00036
Gu, K., Liu, X., Liu, M., Wei, X., Li, J., Hu, Y., et al. (2025). Tobacco intercropping enhances soil fertility by improving synergic interactions between soil physicochemical and microbial properties. Front. Microbiol. 16. doi: 10.3389/fmicb.2025.1647493
Gullner, G., Juhász, C., Németh, A., and Barna, B. (2017). Reactions of tobacco genotypes with different antioxidant capacities to powdery mildew and Tobacco mosaic virus infections. Plant Physiol. Biochem. 119, 232–239. doi: 10.1016/j.plaphy.2017.09.001
Harman, G., Khadka, R., Doni, F., and Uphoff, N. (2021). Benefits to plant health and productivity from enhancing plant microbial symbionts. Front. Plant Sci. 11. doi: 10.3389/fpls.2020.610065
Heil, M. and Bostock, R. M. (2002). Induced systemic resistance (ISR) against pathogens in the context of induced plant defences. Ann. Bot. 89, 503–512. doi: 10.1093/aob/mcf076
Homulle, Z., Cassiano, P., Shevchuk, S., Anten, N. P. R., Stomph, T. J., Van Der Werf, W., et al. (2025). Disease-suppressive mechanisms in contrasting potato-based strip-cropping systems. Eur. J. Plant Pathol. 173, 1–21. doi: 10.1007/s10658-025-02872-4
Imran, M., Sun, X., Hussain, S., Ali, U., Rana, M. S., Rasul, F., et al. (2019). Molybdenum-induced effects on nitrogen metabolism enzymes and elemental profile of winter wheat (Triticum aestivum L.) under different nitrogen sources. Int. J. Mol. Sci. 20, 3009. doi: 10.3390/ijms20123009
Iñiguez, C., Aguiló-Nicolau, P., and Galmés, J. (2021). Improving photosynthesis through the enhancement of Rubisco carboxylation capacity. Biochem. Soc Trans. 49, 2007–2019. doi: 10.1042/BST20210080
Javed, T. and Gao, S. J. (2023). WRKY transcription factors in plant defense. Trends Genet. 39, 787–801. doi: 10.1016/j.tig.2023.05.006
Jiang, F., Xiao, X., Tang, L., Kuang, S., Xu, Y., Song, W., et al. (2024). The response of tobacco mineral composition and absorption to application of potassium sulfate fertilizer. Ind. Crops Prod 210, 118155. doi: 10.1016/j.indcrop.2024.118155
Jin, L., Cui, S., Du, J., Jin, C., Wu, Y., and Qirige (2009). Effects of drought stress on PAL and C4H activities in leaves of Amygdalus mongolica under different ecological conditions. Acta Agric. Boreali-Sin 24, 118–122.
Kebede, E. (2021). Contribution, utilization, and improvement of legumes-driven biological nitrogen fixation in agricultural systems. Front. Sustain. Food Syst. 5. doi: 10.3389/fsufs.2021.767998
Kramer, D. M. and Evans, J. R. (2011). The importance of energy balance in improving photosynthetic productivity. Plant Physiol. 155, 70–78. doi: 10.1104/pp.110.165076
Latowski, D., Kuczyńska, P., and Strzałka, K. (2011). Xanthophyll cycle–a mechanism protecting plants against oxidative stress. Redox Rep. 16, 78–90. doi: 10.1179/174329211X13020951739938
Lee, S. M., Kong, H. G., Song, G. C., and Ryu, C. M. (2021). Disruption of Firmicutes and Actinobacteria abundance in tomato rhizosphere causes the incidence of bacterial wilt disease. ISME J. 15, 330–347. doi: 10.1038/s41396-020-00776-x
Li, P., Cheng, L., Gao, H., Jiang, C., and Peng, T. (2009). Heterogeneous behavior of PSII in soybean (Glycine max) leaves with identical PSII photochemistry efficiency under different high temperature treatments. J. Plant Physiol. 166, 1607–1615. doi: 10.1016/j.jplph.2009.04.012
Li, Z. and Gong, M. (2008). Improvement of guaiacol method for determining peroxidase activity in plants. Plant Physiol. Commun. 44, 323–324. doi: 10.13592/j.cnki.ppj.2008.02.018
Li, Y., Gu, X., Yong, T., and Yang, W. (2024b). A global synthesis reveals additive density design drives intercropping effects on soil N-cycling variables. Soil Biol. Biochem. 191, 109318. doi: 10.1016/j.soilbio.2024.109318
Li, X., Huang, R., Wang, Y., Jiang, H., Luo, Y., Wang, C., et al. (2025). Straw application promotes soil carbon storage by affecting aggregate-associated bacterial community structure and RuBisCO activity: a 35-year field experiment. Microbiol. Spectr. 13, e00088–e00025. doi: 10.1128/spectrum.00088-25
Li, L., Tilman, D., Lambers, H., and Zhang, F. S. (2014). Plant diversity and overyielding: insights from belowground facilitation of intercropping in agriculture. New Phytol. 203, 63–69. doi: 10.1111/nph.12778
Li, B., Yang, P., Feng, Y., Du, C., Qi, G., and Zhao, X. (2024a). Rhizospheric microbiota of suppressive soil protect plants against Fusarium solani infection. Pest Manage. Sci. 80, 4186–4198. doi: 10.1002/ps.8161
Lichtenthaler, H. K. and Wellburn, A. R. (1983). Determinations of total carotenoids and chlorophylls a and b of leaf extracts in different solvents. Biochem. Soc Trans. 11, 591–592. doi: 10.1042/bst0110591
Lin, P., Wang, J., Chen, P., Fu, Z., Luo, K., Li, Y., et al. (2024). Relay intercropped soybean promotes nodule development and nitrogen fixation by root exudate deposition. Front. Plant Sci. 15. doi: 10.3389/fpls.2024.1447447
Liu, A., Contador, C. A., Fan, K., and Lam, H. M. (2018). Interaction and regulation of carbon, nitrogen, and phosphorus metabolisms in root nodules of legumes. Front. Plant Sci. 9. doi: 10.3389/fpls.2018.01860
Luo, C., Ma, L., Zhu, J., Guo, Z., Dong, K., and Dong, Y. (2021). Effects of nitrogen and intercropping on the occurrence of wheat powdery mildew and stripe rust and the relationship with crop yield. Front. Plant Sci. 12. doi: 10.3389/fpls.2021.637393
Mandal, S. M., Chakraborty, D., and Dey, S. (2010). Phenolic acids act as signaling molecules in plant–microbe symbioses. Plant Signal. Behav. 5, 359–368. doi: 10.4161/psb.5.4.10871
Molefe, R. R., Amoo, A. E., and Babalola, O. O. (2023). Communication between plant roots and the soil microbiome: involvement in plant growth and development. Symbiosis 90, 231–239. doi: 10.1007/s13199-023-00936-1
Nelson, M. S. and Sadowsky, M. J. (2015). Secretion systems and signal exchange between nitrogen-fixing rhizobia and legumes. Front. Plant Sci. 6. doi: 10.3389/fpls.2015.00491
Niinemets, Ü. (2023). Variation in leaf photosynthetic capacity within plant canopies: optimization, structural, and physiological constraints and inefficiencies. Photosynth. Res. 158, 131–149. doi: 10.1007/s11120-023-01048-2
Paul, M. J. and Pellny, T. K. (2003). Carbon metabolite feedback regulation of leaf photosynthesis and development. J. Exp. Bot. 54, 539–547. doi: 10.1093/jxb/erg052
Peng, I. C., Bott, A. J., and Zong, W. X. (2016). Spectrophotometric determination of glutamine synthetase activity in cultured cells. Bio-Protocol 6, e1959. doi: 10.21769/BioProtoc.1959
Phillips, S. L., Shaw, M. W., and Wolfe, M. S. (2005). The effect of potato variety mixtures on epidemics of late blight in relation to plot size and level of resistance. Annals of Applied Biology 147 (3), 245–252. doi: 10.1111/j.1744-7348.2005.00027.x
Qiu, Z., Guo, J., Zhu, A., Zhang, L., and Zhang, M. (2014). Exogenous jasmonic acid can enhance tolerance of wheat seedlings to salt stress. Ecotoxicol. Environ. Saf. 104, 202–208. doi: 10.1016/j.ecoenv.2014.02.002
Qu, M., Qin, L., Liu, Y., Fan, H., Zhu, S., and Wang, J. (2014). Comparison of two methods for determination of SOD enzyme activity. J. Food Saf. Qual 5, 3318–3323. doi: 10.19812/j.cnki.jfsq11-5956/ts.2014.10.063
Sakoda, K., Yamori, W., Groszmann, M., and Evans, J. R. (2021). Stomatal, mesophyll conductance, and biochemical limitations to photosynthesis during induction. Plant Physiol. 185, 146–160. doi: 10.1093/plphys/kiaa014
Semer, J., Navrátil, M., Špunda, V., and Štroch, M. (2019). Chlorophyll fluorescence parameters to assess utilization of excitation energy in photosystem II independently of changes in leaf absorption. J. Photochem. Photobiol. B Biol. 197, 111535. doi: 10.1016/j.jphotobiol.2019.111535
Singh, B. K., Delgado-Baquerizo, M., Egidi, E., Guirado, E., Leach, J. E., Liu, H., et al. (2023). Climate change impacts on plant pathogens, food security and paths forward. Nat. Rev. Microbiol. 21, 640–656. doi: 10.1038/s41579-023-00936-5
State Tobacco Monopoly Administration (2002a). YC/T 159–2002. Tobacco and tobacco products–Determination of water soluble sugars–Continuous flow method (Beijing: State Tobacco Monopoly Administration).
State Tobacco Monopoly Administration (2002b). YC/T 161–2002. Tobacco and tobacco products–Determination of total nitrogen–Continuous flow method (Beijing: State Tobacco Monopoly Administration).
State Tobacco Monopoly Administration (2007). YC/T 217–2007. Tobacco and tobacco products–Determination of potassium–Continuous flow method (Beijing: State Tobacco Monopoly Administration).
State Tobacco Monopoly Administration (2011). YC/T 162–2011. Tobacco and tobacco products–Determination of chloride–Continuous flow method (Beijing: State Tobacco Monopoly Administration).
State Tobacco Monopoly Administration (2013). YC/T 468–2013. Tobacco and tobacco products–Determination of total alkaloids–Continuous flow (potassium thiocyanate) method (Beijing: State Tobacco Monopoly Administration).
Sugiyama, A. (2019). The soybean rhizosphere: Metabolites, microbes, and beyond–A review. J. Adv. Res. 19, 67–73. doi: 10.1016/j.jare.2019.03.003
Tu, Y. (2015) Field configuration techniques of soybean intercropped with flue-cured tobacco and their effects on yield, quality, and occurrence of diseases and pests. China: Sichuan Agricultural University).
Wang, L., Cheng, B., Zhou, T., Jing, S., Liu, R., Gao, Y., et al. (2023). Quantifying the effects of plant density on soybean lodging resistance and growth dynamics in maize–soybean strip intercropping. Front. Plant Sci. 14. doi: 10.3389/fpls.2023.1264378
Wang, Z., Guo, X., Cao, S., Yang, M., Gao, Q., Zong, H., et al. (2024). Tobacco/Isatis intercropping system improves soil quality and increase total production value. Front. Plant Sci. 15. doi: 10.3389/fpls.2024.1458342
Wu, W., Liao, Y., Shah, F., Nie, L., Peng, S., and Cui, K.. (2013). Plant growth suppression due to sheath blight and the associated yield reduction under double rice-cropping system in central China. Field Crops Research 144, 268–280. doi: 10.1016/j.fcr.2013.01.012
Yamori, W., Shikanai, T., and Makino, A. (2015). Photosystem I cyclic electron flow via chloroplast NADH dehydrogenase-like complex performs a physiological role for photosynthesis at low light. Sci. Rep. 5, 13908. doi: 10.1038/srep13908
Yang, Y., Zhang, H., Fang, Y., Li, Y., Mei, X., Huang, H., et al. (2021). Interference by non-host plant roots and root exudates in the infection processes of Phytophthora nicotianae. Front. Agric. Sci. Eng. 8, 474–486. doi: 10.15302/J-FASE-2021438
Yang, M., Zhang, Y., Qi, L., Mei, X., Liao, J., Ding, X., et al. (2014). Plant–plant–microbe mechanisms involved in soil-borne disease suppression in a maize and pepper intercropping system. PloS One 9, e115052. doi: 10.1371/journal.pone.0115052
Yao, X., Zhou, H., Zhu, Q., Li, C., Zhang, H., Wu, J. J., et al. (2017). Photosynthetic response of soybean leaf to wide light fluctuation in a maize-soybean intercropping system. Front. Plant Sci. 8. doi: 10.3389/fpls.2017.01695
Yoneyama, T. and Suzuki, A. (2019). Exploration of nitrate-to-glutamate assimilation in non-photosynthetic roots of higher plants by studies of 15N-tracing, enzymes involved, reductant supply, and nitrate signaling: a review and synthesis. Plant Physiol. Biochem. 136, 245–254. doi: 10.1016/j.plaphy.2018.12.030
Zhang, G., Li, Z., Zhu, Q., Yang, C., Shu, H., Gao, Z., et al. (2025). Cropping patterns and plant population density alter nitrogen partitioning among photosynthetic components, leaf photosynthetic capacity and photosynthetic nitrogen use efficiency in field-grown soybean. Ind. Crops Prod 226, 120680. doi: 10.1016/j.indcrop.2024.120680
Zhang, J., Shuang, S., Zhang, L., Xie, S., and Chen, J. (2021). Photosynthetic and photoprotective responses to steady-state and fluctuating light in the shade-demanding crop Amorphophallus xiei grown in intercropping and monoculture systems. Front. Plant Sci. 12. doi: 10.3389/fpls.2021.663473
Zheng, Y., Guo, Y., Lv, J., Dong, K., and Dong, Y. (2022b). Mitigation of vanillic acid–promoted faba bean Fusarium wilt by faba bean–wheat intercropping. Plant Pathol. 71, 830–842. doi: 10.1111/ppa.13568
Zheng, B., Zhao, W., Ren, T., Zhang, X., Ning, T., Liu, P., et al. (2022a). Low light increases the abundance of light reaction proteins: proteomics analysis of maize (Zea mays L.) grown at high planting density. Int. J. Mol. Sci. 23, 3015. doi: 10.3390/ijms23063015
Zhou, M., Sun, C., Dai, B., He, Y., and Zhong, J. (2023a). Intercropping system modulated soil–microbe interactions that enhanced the growth and quality of flue-cured tobacco by improving rhizospheric soil nutrients, microbial structure, and enzymatic activities. Front. Plant Sci. 14. doi: 10.3389/fpls.2023.1233464
Zhou, X., Zhang, J., Rahman, M. K. U., Gao, D., Wei, Z., Wu, F., et al. (2023b). Interspecific plant interaction via root exudates structures the disease suppressiveness of rhizosphere microbiomes. Mol. Plant 16, 849–864. doi: 10.1016/j.molp.2023.03.007
Keywords: tobacco–soybean intercropping, intercropping density, photosynthetic efficiency, carbon and nitrogen metabolism, disease resistance, flue-cured tobacco quality
Citation: Liu X, Hu Y, Yang C, Li J, Lu C, Tian N, Zhou H, Jin S, Su J, Wang D, Xu C, Huang Y and Liu M (2025) Effects of soybean intercropping density on photosynthetic characteristics and disease resistance in tobacco. Front. Plant Sci. 16:1724956. doi: 10.3389/fpls.2025.1724956
Received: 14 October 2025; Accepted: 12 November 2025; Revised: 06 November 2025;
Published: 03 December 2025.
Edited by:
Bibi Rafeiza Khan, University of Scranton, United StatesReviewed by:
Bartholomew Saanu Adeleke, Olusegun Agagu University of Science and Technology, NigeriaFaisal Mehdi, Yunnan Academy of Agricultural Sciences, China
Copyright © 2025 Liu, Hu, Yang, Li, Lu, Tian, Zhou, Jin, Su, Wang, Xu, Huang and Liu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ming Liu, bGl1bTA2MTVAMTYzLmNvbQ==
†These authors have contributed equally to this work
Yanxia Hu1†