Plant materials and growth conditions
This study was conducted in 2024 at the experimental station of the National Coarse Cereals Engineering Technology Research Center, located in the High-Tech Industrial Development Zone of Daqing City, Heilongjiang Province, China. The soybean (Glycine max L.) cultivar ‘Heihe 43’, a predominant cultivar in Heilongjiang Province characterized by semi-determinate growth habit, was used as plant material. Prior to sowing, uniformly sized seeds free from disease spots and physical damage were selected. Surface sterilization was performed by treating seeds with 5% sodium hypochlorite solution for 10 min, followed by three rinses with sterile distilled water. Plants were grown in plastic pots (height: 33 cm; diameter: 30 cm). To prevent waterlogging, five drainage holes (1 cm diameter) were drilled at the bottom of each pot, which was then lined with a mesh screen to contain the root system within the pot. The pots were filled with quartz sand that had been pre-washed with tap water to remove impurities and subsequently rinsed twice with distilled water. The sand was filled to a level 7 cm below the rim of the pot. Before sowing, each pot was irrigated with sufficient distilled water to achieve complete saturation of the sand substrate. Nine seeds were evenly placed on the sand surface and covered with a 2 cm layer of sand. To avoid potential effects of natural precipitation on experimental conditions, all pots were maintained under a movable rain-out shelter throughout the experiment.
Experimental Design
From sowing until the emergence of the first true leaf, each pot received daily irrigation with 500 mL of distilled water. At the full expansion of the first true leaf, three uniformly growing seedlings per pot were retained, and the remaining ones were carefully removed. Upon reaching the V1 stage (first trifoliate fully unfolded), the cotyledons of all seedlings were excised to eliminate potential confounding effects of residual nitrogen reserves within these organs. The seedlings were then randomly assigned to three experimental groups. Each group received daily irrigation with 500 mL of a modified nutrient solution based on half-strength Hoagland's formulation, differing primarily in nitrogen concentration:
CK (Control): Plants were irrigated with a nutrient solution containing a standard nitrate nitrogen concentration (91 mg∙L⁻¹ NO₃⁻-N).
LN (Low Nitrogen): Plants were subjected to nitrogen deficiency stress by irrigation with a nutrient solution containing one-fifth of the standard nitrate nitrogen concentration (1/5 of CK level).
LN + GABA: Plants initially received the low nitrogen solution (identical to LN). At the V2 stage, this group was treated for three consecutive days with the low nitrogen solution supplemented with 5 mmol∙L⁻¹ GABA. Following this 3-day period, irrigation reverted to the standard low nitrogen solution (without GABA) for the remainder of the experiment.
To prevent salt accumulation in the substrate, all pots were leached every 5 days with 3 L of distilled water. The day of initial GABA application (coinciding with the V2 stage) was designated as day 0 of the treatment period. The experiment was terminated 30 days after the initiation of treatments (Day 30).
Composition of nutrient solutions
The nitrogen sources for the standard nitrogen concentration solution were Ca(NO₃)₂ and KNO₃, applied at concentrations of 328 mg∙L⁻¹ and 252 mg∙L⁻¹, respectively. For the low nitrogen stress solution, the concentrations of Ca(NO₃)₂ and KNO₃ were reduced to 65.64 mg∙L⁻¹ and 50.55 mg∙L⁻¹, respectively. To maintain potassium ion balance, KCl was supplemented at 33.55 mg∙L⁻¹ in the low nitrogen solution. All other macro- and micronutrient components remained identical between the two solutions, consisting of: 53.49 mg∙L⁻¹ (NH₄)₂SO₄, 120.37 mg∙L⁻¹ MgSO₄, 1 mL∙L⁻¹ Fe–EDTA stock solution (prepared by dissolving 5.57 g FeSO₄·7H₂O and 7.45 g Na₂EDTA per liter), 8.6 mg∙L⁻¹ ZnSO₄·7H₂O, 6.2 mg∙L⁻¹ H₃BO₃, 0.08 mg∙L⁻¹ CuSO₄·5H₂O, 22.3 mg∙L⁻¹ MnSO₄, and 0.025 mg∙L⁻¹ Na₂MoO₄·H₂O.
Sampling schedule
Plant samplings were conducted at 0, 20, and 30 days after treatment initiation. Whole plants were harvested at these time points for the determination of morphological parameters, biomass, and nitrogen accumulation. Additionally, leaf and root samples were collected at 5, 10, 20, and 30 days. A subset of these samples was immediately used for the assay of nitrate reductase activity, chlorophyll content, and root activity. Another subset was rapidly frozen in liquid nitrogen and subsequently stored at − 80 ℃ for later analysis of enzymatic activities and other physiological indicators. Photosynthetic gas exchange parameters were measured directly using portable instruments without destructive sampling.
Measurement Indices and Methods
Measurement of morphological and dry matter-related parameters
Plant height and root length were measured using a standard ruler. Stem diameter was determined as the mid-internode diameter of the first fully expanded internode (counting the cotyledonary node as node 0) using a digital vernier caliper. Leaf area was quantified with a Yaxin-1241 leaf area meter (Beijing Yaxin Liyi Technology Co., Ltd., China). For root system analysis, roots were scanned using an EPSON Perfection V800 flatbed scanner (Seiko Epson Corporation, Japan) and parameters including root volume, root surface area, and number of root tips were analyzed with the WinRHIZO Pro 2016a image analysis system (Regent Instruments Inc., Canada).
Plants were separated into leaves, stems, and roots, and each component was placed in individual paper envelopes. Samples were first oven-dried at 105 ℃ for 30 min to deactivate enzymes, followed by drying at 80 ℃ until a constant weight was achieved. Dry matter mass was then measured using an analytical balance. The growth rate from day 0 to day 20 was calculated based on dry weight accumulation using the following formula:
Growth Rate (g∙plant− 1∙day− 1\(\:=\frac{\text{D}{\text{M}}_{20}\:\left({\text{g}\bullet\:\text{p}\text{l}\text{a}\text{n}\text{t}}^{-1}\right)-\text{D}{\text{M}}_{0}\:\left({\text{g}\bullet\:\text{p}\text{l}\text{a}\text{n}\text{t}}^{-1}\right)}{21\:\text{d}\text{a}\text{y}\text{s}}\)
\(\:\text{D}{\text{M}}_{20}\) represents the dry weight of soybean plants on the 20th day, respectively, while \(\:\text{D}{\text{M}}_{0}\)indicates the dry weight of soybean plants on day 0.
Determination of root activity
Root activity was assessed using the triphenyltetrazolium chloride (TTC) reduction assay following the established method (Zhang et al. 2022). Briefly, root systems were carefully washed with deionized water to remove adhering soil particles and gently blotted dry. The cleaned roots were then immersed in a 0.4% (w/v) TTC solution and incubated in the dark at 37℃ for 2–4 hours to allow for adequate TTC penetration and enzymatic reduction. After the incubation period, roots were removed from the TTC solution. The formed formazan (TTF) was extracted from the root tissues using ethyl acetate as the solvent. The extract was transferred to centrifuge tubes and subjected to centrifugation at an appropriate speed to obtain a clear supernatant. The absorbance of the supernatant was measured at a wavelength of 485 nm using a spectrophotometer. Root activity, expressed as TTC reduction intensity in units of µg TTC reduced per gram of root fresh weight per hour (µg TTC·g⁻¹·h⁻¹), was calculated based on the measured absorbance values and a pre-established standard curve.
Photosynthetic Indices
Determination of photosynthetic pigment content
The concentrations of chlorophyll a (Chl a), chlorophyll b (Chl b), total chlorophyll (Chl a + b), and carotenoids (Car) were determined according to the method described by Lichtenthaler (Lichtenthaler 1987). Fully expanded functional leaves (100 mg fresh weight) were cut into segments and immersed in 10 mL of absolute ethanol for 24 hours in darkness until the tissues became completely bleached. The optical density (OD) of the extracts was measured at wavelengths of 470, 649, and 665 nm using a Jenway 6850 UV-Vis spectrophotometer (Cole-Parmer Ltd., UK). Pigment concentrations were calculated using the following equations:
Chl a (µg·mL⁻¹) = 13.95 × OD₆₆₅ – 6.88 × OD₆₄₉
Chl b (µg·mL⁻¹) = 24.96 × OD₆₄₉ – 7.32 × OD₆₆₅
Total chlorophyll (µg·mL⁻¹) = Chl a + Chl b
Car (µg·mL⁻¹) = (1000 × OD₄₇₀ – 2.05 × Chl a – 111.48 × Chl b) / 245
Final pigment contents were expressed per gram fresh weight of the leaf tissue.
Measurement of photosynthetic gas exchange parameters
Gas exchange parameters were measured between 9:00 and 11:00 AM following the methodology described by Wang (Wang et al. 2021). Measurements were conducted on the second fully expanded leaf from the apex of the main stem using a Li-6400 portable photosynthesis system (LI-COR Biosciences, Lincoln, NE, USA) equipped with a red-blue LED light source chamber. The parameters recorded included net photosynthetic rate (Pₙ), transpiration rate (Tr), stomatal conductance (Gₛ), and intercellular CO₂ concentration (Ci). During measurements, the following environmental conditions were maintained within the leaf chamber: photosynthetic photon flux density (PPFD) of 1000 µmol·m⁻²·s⁻¹, CO₂ concentration of 400 µmol·mol⁻¹, leaf temperature of 25 ℃, and relative humidity of 25%.
Determination of key nitrogen metabolism enzyme activities and nitrogen-containing compounds
Nitrate reductase (NR) activity was determined according to Mancuso (Mancuso and Caviness
et al. 1991). Fresh leaf samples (0.5 g) were homogenized in 5 mL of phosphate buffer (0.1 mol·L⁻¹, pH 7.5) and 5 mL of potassium nitrate solution (0.2 mol·L⁻¹). The reaction mixture was incubated in darkness at 25 ℃ for 1 h and terminated by adding 1 mL of 30% (w/v) trichloroacetic acid (TCA). Then, 2 mL of the reaction mixture was combined with 8 mL of nitration reagent, incubated at 20 ℃, and the absorbance was measured at 540 nm.
Glutamine synthetase (GS) activity was assayed following the method of Ribeiro (Ribeiro et al. 2000). Fresh leaf or root samples (0.1 g) were ground into powder with liquid nitrogen and extracted with 8 mL of extraction buffer (100 mmol·L⁻¹ Tris-HCl, 0.5 mmol·L⁻¹ EDTA, 5 mmol·L⁻¹ β-mercaptoethanol, pH 7.5). The homogenate was centrifuged at 15,000 × g for 20 min at 4 ℃, and the supernatant was used for enzyme activity determination. A mixture of 1.6 mL reaction buffer and 0.6 mL enzyme extract was pre-incubated at 25 ℃ for 5 min. The reaction was initiated by adding 0.2 mL of hydroxylamine reagent, continued for 15 min at 25 ℃, and stopped with 1 mL of FeCl₃ reagent.
Glutamate synthase (GOGAT) activity was measured according to An (An et al. 2023), using the same extraction method as for GS. The reaction mixture (3 mL total volume) contained 0.4 mL of 20 mmol·L⁻¹ L-glutamine, 0.05 mL of 0.1 mol·L⁻¹ α-ketoglutarate, 0.1 mL of 10 mmol·L⁻¹ KCl, 0.1 mL of 3 mmol·L⁻¹ NADH, and 0.3 mL enzyme extract, with the volume made up with 25 mmol·L⁻¹ Tris-HCl (pH 7.6). The reaction was initiated by adding L-glutamine, and the decrease in absorbance at 340 nm was recorded every 30 s. Enzyme activity was calculated from the linear decrease in optical density.
Glutamate dehydrogenase (GDH) activity was determined as described by Lin and Kao (Lin and Kao 1996), with extraction identical to GS. The assay mixture contained 2.6 mL of assay stock solution (23.1 mmol·L⁻¹ α-ketoglutarate, 231 mmol·L⁻¹ NH₄Cl, and 115.4 mmol·L⁻¹ Tris-HCl buffer, pH 8.0), 0.1 mL of 6 mmol·L⁻¹ NADH, and 0.1 mL of 30 mmol·L⁻¹ CaCl₂. The reaction was started by adding 0.1 mL enzyme extract, and the decrease in absorbance at 340 nm was monitored for 3 min.
Aspartate aminotransferase (GOT) and alanine aminotransferase (GPT) activities were measured according to Zhu (Zhu et al.2021). Samples (0.1 g) were dried at 70 ℃ to constant weight, ground into powder, and extracted with pre-chilled extraction buffer. After 30 min of ice-bath incubation, the homogenate was centrifuged at 10,000 × g for 20 min at 4 ℃. The supernatant was collected for assay. For GOT activity, 0.5 mL enzyme extract was mixed with 4.5 mL reaction solution containing 0.1 mol·L⁻¹ Tris-HCl buffer (pH 7.8), 0.25 mmol·L⁻¹ EDTA-Na₂, and 0.5 mmol·L⁻¹ DL-α-ketoglutarate. After incubation at 37 ℃ for 1 h, the reaction was stopped with 1 mL of 80 g·L⁻¹ TCA. Then, 0.5 mL of 0.6 mol·L⁻¹ NaOH and 3 mL distilled water were added to 1 mL of the reaction mixture. Subsequently, 1 mL of 0.1 mol·L⁻¹ 2,4-dinitrophenylhydrazine reagent was added, followed by incubation at 37 ℃ for 30 min. Then, 5 mL distilled water and 2.5 mL of 4 mol·L⁻¹ NaOH were added. After centrifugation at 3000 × g for 5 min, the absorbance of the supernatant was measured at 520 nm using a 1-cm cuvette. GPT activity was assayed similarly, except that DL-α-ketoglutarate was replaced with pyruvate in the reaction mixture.
Nitrate nitrogen (NO₃⁻–N) content was determined according to Oliveira (Oliveira et al.2013). Leaf samples (1 g) were extracted in 10 mL of ethanol:chloroform:water (12:5:3, v/v/v) for 24 h. After centrifugation at 2000 × g for 30 min, one-fourth volume of chloroform was added. The aqueous phase was collected after 24 h of phase separation, concentrated at 37 ℃ for 15 h, and used for analysis. Then, 0.1 mL of sample was mixed with 0.4 mL of 5% (w/v) salicylic acid–sulfuric acid reagent, incubated at room temperature for 20 min, followed by slow addition of 9.5 mL of 8% (w/v) NaOH. After cooling to room temperature, absorbance was measured at 410 nm against a blank.
Ammonium nitrogen (NH₄⁺–N) content was measured as described by Chien and Kao (Chien and Kao et al.2000). Fresh plant material (0.1 g) was homogenized in 2 mL of 0.3 mmol·L⁻¹ H₂SO₄ (pH 3.5). The homogenate was centrifuged at 10,000 × g for 10 min. Then, 200 µL of supernatant was diluted to 4 mL with 0.3 mmol·L⁻¹ H₂SO₄ (pH 3.5). Subsequently, 0.5 mL of color reagent A (containing 0.5 g·L⁻¹ phenol and 0.25 g·L⁻¹ sodium nitroprusside) and 0.5 mL of color reagent B (containing 400 mL·L⁻¹ of 5% sodium hypochlorite and 0.25 g·L⁻¹ NaOH) were added. The mixture was incubated at 37 ℃ for 20 min with gentle shaking. Absorbance was measured at 625 nm against a blank prepared with distilled water.
Soluble protein concentration was determined using the Bradford method (Bradford 1976). Samples (0.1 g) were homogenized in 5 mL of phosphate buffer (pH 7.0), kept on ice for 30 min, and centrifuged at 12,000 × g for 15 min at 4 ℃. The supernatant was reacted with Coomassie Brilliant Blue G-250 reagent for 5 min, and absorbance was measured at 595 nm. Soluble protein content was calculated based on a standard curve.
Free amino acid content was measured according to Liu (Liu et al.2019). Samples (0.1 g) were extracted in 5 mL distilled water at 80 ℃ for 30 min. After centrifugation at 4000 × g for 10 min, 1 mL of supernatant was mixed with 2 mL ninhydrin reagent, heated in a boiling water bath for 15 min, cooled, and absorbance was measured at 570 nm.
Determination of key sugar contents
Soluble sugar content was determined using the anthrone-sulfuric acid method (Hu et al.2009). A 0.1 g sample was accurately weighed and extracted with 5 mL distilled water at 80 ℃ for 30 min. After cooling, the homogenate was centrifuged at 4000 rpm for 10 min. One milliliter of supernatant was mixed with 4 mL anthrone reagent and incubated in boiling water for 10 min. Following cooling, absorbance was measured at 620 nm.
Sucrose and fructose contents were determined according to the method of Cao (Cao et al.2019). A 0.1 g sample was weighed into a centrifuge tube, homogenized with 5 mL distilled water, and extracted at 80 ℃ for 30 min. After cooling to room temperature, the mixture was centrifuged at 4000 rpm for 10 min to collect the supernatant. One milliliter of supernatant was transferred into two separate tubes: one was mixed with 2 mL resorcinol reagent and 5 mL concentrated hydrochloric acid (for fructose determination), and the other with 2 mL resorcinol reagent and 5 mL concentrated sulfuric acid (for sucrose determination). The mixtures were incubated in boiling water for 10 min, cooled, and absorbance was measured at 480 nm (fructose) and 620 nm (sucrose), respectively. Sucrose and fructose concentrations were calculated based on standard curves.
Determination of Key Enzyme Activities in Sucrose Metabolism
The activities of sucrose synthase (SS) and sucrose phosphate synthase (SPS) were determined according to the method described by Chopra (Chopra et al.2000). Briefly, 0.1 g of sample was homogenized and incubated in the respective reaction buffers for SS (pH 7.0) or SPS (pH 7.5) at 37°C for 30 min. The reaction was terminated by adding DNS reagent, followed by heating in a boiling water bath for 5 min. After cooling, the absorbance was measured at 540 nm. Enzyme activity was calculated based on a standard curve and expressed as the amount of reducing sugars produced per unit time per gram of fresh weight (µmol·min⁻¹·g⁻¹ FW).
The activities of acid invertase (AI) and neutral invertase (NI) were determined according to the method of Tsai (Tsai et al.1970). Briefly, 0.1 g of fresh sample was homogenized and incubated in phosphate buffer (pH 4.6 for AI or pH 7.0 for NI) at 37°C for 30 min. The reaction was terminated by adding DNS reagent, followed by heating in a boiling water bath for 5 min. After cooling, the absorbance was measured at 540 nm. Enzyme activity was calculated based on a standard curve and expressed as the amount of reducing sugar produced per minute per gram of fresh weight (µmol·min⁻¹·g⁻¹ FW).
Determination of plant nitrogen concentration and calculation of nitrogen accumulation
Plant tissue nitrogen concentration was determined using the Zheng method (Zheng et al.2020). Mature samples, dried at 105 ℃, were ground to a fine powder. Approximately 0.2 g of sample was accurately weighed, digested with H₂SO₄–H₂O₂, and analyzed for nitrogen concentration using an automatic Kjeldahl apparatus.
Nitrogen accumulation (g) was calculated as:
Nitrogen accumulation (g) = Dry matter accumulation (g) × Nitrogen concentration (%)
Determination of plant carbon content
Total carbon content in plant tissues was determined using the potassium dichromate-sulfuric acid oxidation method (Starr et al.1964). Briefly, 0.005 g of dried soybean sample was weighed into a test tube, and a known volume of potassium dichromate-sulfuric acid solution was added. The mixture was allowed to stand for 30 min and then boiled for 10 min to ensure complete oxidation of organic matter. After cooling to room temperature, the solution was transferred to a conical flask, and 2–3 drops of sodium diphenylamine sulfonate indicator were added. The solution was then titrated with a standardized ferrous sulfate solution (of known concentration) until the color changed from purple to green, and the titration volume was recorded. A blank titration was performed simultaneously. The total carbon content in the plant sample was calculated based on the difference in titration volumes between the blank and the sample, the concentration of the ferrous sulfate standard solution (C, mol/L), and the mass of the dried sample (m, g), according to the established formula.
Data Processing and Plotting
Microsoft Excel 2020 was used for data entry and organization, and RStudio was used for data analysis and plotting. Data statistics were performed using one-way analysis of variance (ANOVA) and Duncan's multiple range test (p < 0.05). In Mantel analysis, Pearson's r characterizes the linear dependence between two variables, with values close to 0 indicating no correlation; values between 0.1 and 0.3 or -0.1 and − 0.3 indicate weak correlation; values between 0.3 and 0.5 or -0.3 and − 0.5 indicate moderate correlation; and values greater than 0.5 or less than − 0.5 indicate strong correlation. Mantel's r represents the similarity or dissimilarity between two distance matrices, with threshold ranges similar to those of Pearson's r. Mantel's p-value assists in determining the significance of Mantel's r, with p-values less than 0.05 typically considered statistically significant, indicating that the correlation between the two matrices is unlikely to be due to random chance. A p-value less than 0.01 indicates an even stronger level of significance. Conversely, a p-value greater than 0.05 is generally regarded as not significant, suggesting insufficient evidence to indicate a correlation between the two matrices.