Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Exogenous melatonin reduces the inhibitory effect of osmotic stress on photosynthesis in soybean

  • Mingcong Zhang,

    Roles Writing – original draft

    Affiliation College of Agronomy, Heilongjiang Bayi Agricultural University, Daqing, P.R. China

  • Songyu He,

    Roles Data curation

    Affiliation College of Agronomy, Heilongjiang Bayi Agricultural University, Daqing, P.R. China

  • Yingce Zhan,

    Roles Data curation

    Affiliation College of Agronomy, Heilongjiang Bayi Agricultural University, Daqing, P.R. China

  • Bin Qin,

    Roles Formal analysis

    Affiliation College of Agronomy, Heilongjiang Bayi Agricultural University, Daqing, P.R. China

  • Xijun Jin,

    Roles Software

  • Mengxue Wang,

    Roles Writing – original draft

    Affiliations College of Agronomy, Heilongjiang Bayi Agricultural University, Daqing, P.R. China, College of Agronomy, Heilongjiang Bayi Agricultural University, Daqing, P.R. China

  • Yuxian Zhang ,

    Roles Writing – review & editing

    zyx_lxy@126.com

    Affiliation College of Agronomy, Heilongjiang Bayi Agricultural University, Daqing, P.R. China

  • Guohua Hu,

    Roles Writing – review & editing

    Affiliation College of Agronomy, Heilongjiang Bayi Agricultural University, Daqing, P.R. China

  • Zhanlin Teng,

    Roles Data curation

    Affiliation Huanan Agrotechnical Extension Center, Jiamusi, P.R. China

  • Yaokun Wu

    Roles Investigation

    Affiliation Daqing Branch of Heilongjiang Academy of Sciences, Daqing, P.R. China

Abstract

Understanding the relationship between exogenous melatonin and water deficit stress is crucial for achieving high yields and alleviating the effects of water deficit stress on soybean (Glycine max (L.) Merrill) plants in agriculture. This study investigated the effects of exogenous melatonin on soybean photosynthetic capacity under water deficit stress induced by polyethylene glycol (PEG) 6000. We conducted a potting experiment in 2018 using the soybean (Glycine max L. Merrill) cultivar Suinong 26. We identified the impacts of a concentration of PEG 6000 simulating drought (15%, w/v) and an appropriate melatonin concentration (100 μmol/L) on the growth of soybean seedlings and flowering stages in a preliminary test. We applied exogenous melatonin by foliar spraying and root application to determine the effects on leaf photosynthesis during water deficit stress. Our results indicated that 15% PEG 6000 had an obvious inhibitory effect on the growth of soybean seedlings and flowering stages, causing oxidative stress and damage due to reactive oxygen species (ROS) (H2O2 and O2·-) accumulation and potentially reducing air exchange parameters and photosystem II (PSII) efficiency. The application of exogenous melatonin significantly relieved the inhibitory effects of PEG 6000 stress on seedlings and flowering growth, and gas exchange parameters, potentially improved PSII efficiency, improved the leaf area index (LAI) and the accumulation of dry matter, slowed down oxidative stress and damage to leaves by increasing the activity of antioxidant enzymes (SOD, POD, and CAT), reduced the content of malondialdehyde (MDA), and ultimately improved soybean yield. Overall, the results of this study demonstrated that application of exogenous melatonin at the seedlings and flowering stages of soybean is effective in alleviating plant damage caused by water deficit stress and improving the drought resistance of soybean plants. In addition, the results showed that application of exogenous melatonin by root is superior to foliar spraying.

Introduction

Drought is a major threat to plant growth and crop yields [1]. In the soybean (Glycine max L. Merrill) production process, water shortage directly reduces photosynthesis in leaves by affecting plant growth and development [2]. For example, the absorption of CO2 is limited by stomata and mesophyll organs [3]. In addition, under conditions of water deficit stress, the expansion of leaf cells is reduced, and the absorption of CO2 is limited to inhibit the growth of plants [4]. Excessive reactive oxygen species (ROS) are formed and rapidly accumulate in plants [5]. ROS (H2O2 and O2·-) inhibit leaf cell wall enzymes to affect the enzymatic activity of leaf photosynthesis [6]. The tolerance of plants to water deficit varies among different crops and varieties [7].

Water deficit stress reduces net photosynthetic rates in leaves by reducing nonstomatal factors such as ribulose bisphosphate carboxylase (Rubisco) and ribulose-1,5-bisphosphate (RuBP) [8]. Some crops close their stomata to reduce plant damage under drought conditions [9], which reduces the amount of CO2 absorption to decrease the photosynthetic capacity of mesophyll organs [10]. However, plants can adapt their photosynthetic metabolism upon application of exogenous hormones to survive under stress [11]. It has been reported that under water deficit stress, compared with no treatment, the application of exogenous hormones can significantly alleviate growth and developmental impairments in soybean and improve the activity of antioxidant enzymes [12] to enhance the ability of endogenous antioxidants to remove excessive ROS [8, 13].

Chlorophyll fluorescence (ChlF) is used to evaluate the integrity and efficiency of photosynthetic organs as well as the overall health of plant tissues [14]. Photosynthetic organs, especially PSII, are sensitive to water stress [1516]. ChlF emission changes, which are mainly caused by PSII, provide information concerning almost all aspects of photosynthesis [1718] and reflect a plant’s tolerance to environmental stresses, including drought. Fluorescence parameters have been used to identify leaf damage in the absence of obvious symptoms [18], and ChlF is often used as a potential index of environmental stress and a screening method for disease-resistant plants [19]. Fluorescence efficiency can be described by the optimal/maximal quantum yield (Fv/Fm), nonphotochemical quenching (NPQ), electron transport rate (ETR) and actual photochemical quantum yield (ΦPSII) [20]. Studies have shown that the ChlF parameter can be used to screen drought-tolerant germplasm resources of soybean and to evaluate the degree of drought stress present in soybean [21].

It is difficult to control watering to reduce soil moisture at the vegetative or reproductive growth stage of soybean with drought stress. Li [22] suggested that vermiculite and perlite be used to replace agricultural soil, and polyethylene glycol (PEG) with a molecular weight of 6,000 was used to simulate drought stress in plants at any stage of soybean growth. PEG 6000 is an inert substance that cannot enter cell wall tissues, and it has been demonstrated that crops have difficulty absorbing PEG with a molecular weight exceeding 3000 [23]. As reported, PEG 6000 causes obvious water stress to plants without any toxic effects [24]. Our preliminary experiment [25] confirmed that a PEG 6000 concentration of 15% (w/v) had a strong water deficit stress effect on Suinong 26 at V3 stage; however, application of the appropriate melatonin concentration (100 μmol/L) was effective for alleviating plant damage caused by water deficit stress (Fig 1).

thumbnail
Fig 1. Effect of melatonin addition on the growth of soybean plant in our preliminary experiment.

WW: half-strength Hoagland’s nutrient solution alone; PEG: half-strength Hoagland’s nutrient solution plus 15% PEG 6000 treatment; PEG+ML: half-strength Hoagland’s nutrient solution plus 15% PEG 6000 treatment and 100 μmol/L melatonin by leaf; PEG+MR: half-strength Hoagland’s nutrient solution plus 15% PEG 6000 treatment and 100 μmol/L melatonin by root.

https://doi.org/10.1371/journal.pone.0226542.g001

Melatonin (N-acetyl-5-methoxytryptamine), also known as the pineal hormone [26], was first detected in plants by Hattori [2728]. Subsequent studies have detected melatonin in seeds, roots, fruits and leaves of plants [2930]. One important function of melatonin is as an antioxidant that can alleviate plant damage under adverse conditions. Such antioxidant effects have been reported in several plants (e.g., apple [31], maize [32], wheat [33] and cucumber [34]). Research has shown that under water stress, melatonin can directly scavenge ROS and also improve the activity of antioxidant enzymes to alleviate the oxidative damage caused by water stress in cells [35]. Recent research results have suggested that melatonin can improve plant photosynthetic rates under water stress [36], inhibit chlorophyll degradation, and increase leaf photosynthetic rates and the accumulation of dry matter [37]. Previous studies have mainly focused on the beneficial effects of melatonin on plant stress resistance in maize [38], rice [39] and horticultural [3] crops, while few reports have investigated the effect of exogenous melatonin on photosynthetic reactions of soybean seedlings and flowering stages under water stress or on the adaptation of photosynthetic organs to water stress.

Therefore, the purpose of this study was to clarify the adaptive mechanism of soybean seedlings and flowering stages to water stress under exogenous melatonin to provide a theoretical basis for improving the drought resistance of soybean seedlings and flowering stages in actual production. The mechanism was clarified by analyzing the photosynthetic rates and fluorescence parameters of soybean at the V1, V3 and flowering stages under PEG 6000 stress.

Materials and methods

Test material

The soybean variety Suinong 26 was used as the test material in this experiment. This variety features infinite podding and a growth period of approximately 120 d. Melatonin was purchased from Sigma-Aldrich (St. Louis, MO, USA); the molecular formula of melatonin is C13H16N2O2, and the relative molecular mass is 232.28.

Test method

The following experiments were conducted between May and October 2018 at the rainproof greenhouse of the National Engineering and Technology Research Center for Grains in the High-tech Zone, Daqing City, Heilongjiang Province, China (latitude 45°46' to 46°55', longitude 124°19' to 125°12', elevation 146 m). The experiment adopted a potting method. The pots were 36 cm tall with a diameter of 29 cm. On the bottoms of the pots, 5 holes (0.5 cm diameter) were drilled, and gauze was applied. Vermiculite and perlite were mixed 1:1 to serve as the culture medium. Full seeds with uniform colors and sizes were selected and disinfected with 5% sodium hypochlorite before they were sown on May 18. Nine seeds were sown in each pot. During the seedling period, 3 seedlings were consistently maintained evenly in each pot. The ripe seeds were harvested on October 2.

Experimental treatments

The water deficit stress treatment was carried out at the V1 (first unrolled trifoliate leaf and June 4), V3(third unrolled trifoliate leaf and June 17) and R1 (beginning of the flowering stage, fifth unrolled trifoliate leaf and June 30) stages of soybean in this test, and the experimental plants were divided into six groups as follows: (i) the well-watered (WW) group, which was watered with 1 L of half-strength Hoagland’s nutrient solution every 3 d during the whole growth period (water leaked out of the bottom); (ii) the well-watered + melatonin by leaf (WW+ML) group, which received half-strength Hoagland’s nutrient solution and was treated with 100 μmol/L melatonin by leaf spraying (on the top and bottom sides of all leaves until water dripped); (iii) the well-watered + melatonin by root irrigation (WW+MR) group, which received half-strength Hoagland’s nutrient solution and was treated with 100 μmol/L melatonin by root; (iv) the drought stress (PEG) group, which received water deficit stress treatment with half-strength Hoagland’s nutrient solution plus 150 g of PEG 6000 (15%, w/v) (drought; water potential of -0.30 MPa); (v) the drought stress + melatonin by leaf (PEG+ML) group, which received half-strength Hoagland’s nutrient solution plus 15% PEG 6000 and was treated with 100 μmol/L melatonin by leaf spraying (on the top and bottom sides of all leaves until water dripped); and (vi) the drought stress + melatonin by root irrigation (PEG+MR) group, which received half-strength Hoagland’s nutrient solution plus 15% PEG 6000 and was treated with 100 μmol/L melatonin by root (in the preliminary test, it had already been demonstrated that 100 μmol/L melatonin had an ameliorative effect on soybean under water deficit stress). The water deficit stress treatment was conducted on day 0 and day 3 of the V1, V3 and R1 stages of soybean growth by applying half-strength Hoagland’s nutrient solution with 15% PEG 6000. Sampling and photosynthetic parameter determination were conducted at 9:00 am on the 6th day (nutrient solution not applied) after treatment of the V1, V3 and R1 stages. On the evening of the sampling day of sampling, clean water was used to flush the residual PEG 6000 until the solution was consistent with the WW solution, and 1/2 Hoagland’s nutrient solution was normally applied in the later stage to ensure normal growth and development.

Testing items and methods

Leaf area index (LAI) analysis.

Altogether, 5 pots of soybean were selected at random. A LI-3100C area meter (LiCor, Lincoln, NE) was used to measure the areas of all the leaves.

Leaf area duration (LAD) analysis.

LAD (m2·d·m-2) = [(LA1+LA2)/2] ×(t2-t1), LA1 and LA2 represent the former and latter measurements of leaf area, and t1 and t2 represent the times at the former and latter points.

Dry matter accumulation analysis.

Different soybean organs (e.g., leaves, petioles, stems, and roots) were subjected to 105 °C treatment for 30 min and dried at 70 °C to a constant weight, and the biomasses of the different parts of the plants were weighed. Each treatment had 5 replicates.

Determination of antioxidant enzyme activity and membrane lipid peroxidation.

Five pots of soybean were randomly selected, and fully grown 2nd leaves were taken from each plant, wrapped in foil, and stored at -80°C. The activity of superoxide dismutase (SOD) was determined by using nitroblue tetrazole (NBT) method described by Beauchamp et al.[40], the activity of peroxidase (POD) was determined by using the guaiacol method described by Rao et al.[41], the activity of catalase (CAT) was determined by using the hydrogen peroxide method described by Hamurcu et al.[42], and the levels of malondialdehyde (MDA) were determined by the thiobarbituric acid method described by Draper et al.[43].

The superoxide anion production rate was determined. Briefly, 0.5 g of leaf tissue was ground into 10 ml of prechilled PBS (pH 7.8), and each homogenate was centrifuged at 13,000 x g for 20 minutes at 4°C. For the following analysis, 0.5 ml of the supernatant was added to 0.5 ml of PBS and 1 ml of 10 mmol/L hydroxylamine hydrochloride. After 1 hour, 2 ml of ether and 1 ml of 7 mmol α-naphthylamine were added at 4°C, and the mixtures were centrifuged at 3,000 x g for 20 minutes. The absorbance was measured at 530 nm using 0.5 ml of PBS as a reference blank according to the methods of Ke [44]. Each treatment was repeated 5 times.

Histochemical detection of H2O2 and O2·- in soybean leaves.

Detection of hydrogen peroxide (H2O2) in soybean leaves was performed by 3,3-diaminobenzidine (DAB) staining, and detection of superoxide radical (O2·-) was performed by NBT staining [45]. Soybean leaves were immersed and infiltrated under vacuum with 1.25 mg/mL DAB staining solution (pH 7.8) dissolved in H2O for 6 h and with 3 mg/mL nitroblue tetrazolium (NBT) staining solution in 10 mM potassium phosphate buffer (pH 7.0) containing 10 mM NaN3 for 30 min at room temperature. The stained leaves were bleached in an acetic acid: glycerol: ethanol (1:1:3 v/v) solution at 100°C for 5 min and stored in a glycerol: ethanol (1:4 v/v) solution until photographed.

Measurement of gas exchange parameters.

The net photosynthetic rate (Pn), stomatal conductance (Gs), transpiration rate (Tr), and intercellular carbon dioxide concentration (Ci) of the middle part of the middle lobule of the 2nd fully unfolded functional compound leaf from the main stem of each soybean plant were measured using a Li-6400 portable photosynthetic assay system (LICOR Inc., USA) between 09:00 and 11:30. The light source was an internally installed red-blue light, the plants were exposed to a constant photosynthetically active radiation of 1,200 μmol/m2/s, the air velocity was 500 mol/s, the CO2 supply concentration was 400 μmol/mol, the leaf chamber temperature was 25°C, the relative humidity was approximately 25% [46], and each treatment had 5 replicates.

Chlorophyll fluorescence.

Between 09:00 and 11:30 on the date of sampling, an portable chlorophyll fluorometer (FMS-2, Hansatech, Englang) was used to detect the chlorophyll fluorescence parameters of the middle part of the middle lobule of the 2nd fully unfolded functional compound leaf from the main stem of each soybean plant. Before measurement, the leaves were dark-adapted for 20 min; the minimum fluorescence under dark adaptation (Fo) was measured with low light, and then the maximum fluorescence (Fm) was measured with a saturating light pulse. When the fluorescence dropped from Fm to nearly Fo, the steady-state fluorescence (Fs) was measured with continuous actinic light. Later, a beam of saturating light was overlaid to detect the maximum fluorescence (Fm') under light. Finally, the actinic light was stopped, and far-red light was immediately turned on to measure the minimum fluorescence (Fo') under light and the electron transfer rate (ETR).

According to calculations based on Baker’s [46] chlorophyll fluorescence kinetic curve, variable fluorescence (Fv) = Fm–Fo, PSII primary light energy conversion efficiency (Fv/Fm) = (Fm-Fo)/Fm, the PSII photochemical quantum yield (ΦPSII) = (Fm'–Fs)/Fm', and the nonphotochemical quenching coefficient (NPQ) = Fm/Fm'-1. Five plants were measured under each treatment method, and each leaf was recorded three times. The average value was taken as the value of the photosynthetic parameter.

Yield and yield composition.

Sampling were conducted at harvest time. Five pots of soybean under each treatment were randomly selected for assessment of plant height, number of pods per plant, total number of pods, and hundred-grain weight.

Statistical analysis

Analysis of variance was performed by using Microsoft Office Excel 2019. The data are reported as the mean and standard error (SE) values of triplicate experiments and were analyzed by using SPSS 17.0. Function plots were constructed using Origin 2018. One-way analysis of variance (ANOVA) and Duncan’s multiple range tests were used to determine the significance of the differences among samples with a significance level of 0.05 or 0.01.

Results

Effect of exogenous melatonin on individual LAI and LAD values of soybean under osmotic stress

According to Table 1, LAI and LAD in stages V1, V3 and R1 showed no significant difference between the well-watered with melatonin (WW+ML and WW+MR) and the simple well-watered (WW) application; however, the PEG-simulated water stress treatments (PEG, PEG+ML, and PEG+MR) significantly reduced the LAI per plant in the seedling and flowering stages, and the LAI per plant values under V1, V3 and R1 water stress were significantly lower than under the WW condition. The amplitude was decreased by 19.7%~22.8% (P < 0.01) at the V1 stage, 26.0%~37.8% (P < 0.01) at the V3 stage and 36.3%~47.3% (P < 0.01) at the R1 stage. Application of exogenous melatonin could increase the LAI per plant in the soybean seedling stage. The LAI per plant values were higher in the groups treated with melatonin (PEG+ML and PEG+MR) under water stress at the V1 stage than in the groups without melatonin treatment (PEG), but the difference did not reach significance (5%). In stage V3, compared with PEG treatment, melatonin treatment by root irrigation (PEG+MR) increased the LAI per plant by 9.37% (P<0.05), while melatonin treatment by foliar spraying (PEG+ML) increased the value by 3.45% under water stress; this difference, however, did not reach significance. Additionally, the LAI per plant values were higher with the application of exogenous melatonin treatment (PEG+ML and PEG+MR) under water stress at the R1 stage than without melatonin treatment (PEG), with increases by 5.24% and 8.07% (P<0.05), respectively. Water stress markedly reduced the LAD of leaves at the V1-V3 and R1-V3 stages. Compared with the WW condition, PEG treatment decreased the LAD by 56.3% and 69.0% (P<0.05), respectively, and the application of exogenous melatonin increased the LAD of soybean leaves at the V1-V3 and R1-V3 stages under water stress. Compared with PEG, PEG+ML and PEG+MR increased the LAD by 4.68% and 4.90% (P<0.05) in the V1-V3 stage, and PEG+MR increased the LAD by 5.04% (P<0.05) in the R1-V3 stage.

thumbnail
Table 1. Effect of exogenous melatonin on individual LAI and LAD values of soybean under osmotic stress.

https://doi.org/10.1371/journal.pone.0226542.t001

Effect of exogenous melatonin on dry matter accumulation of soybean under osmotic stress

The accumulation of dry matter in various organs showed no significant difference between the well-watered with melatonin application (WW+ML and WW+MR) and simple well-watered (WW) groups (Table 2); however, compared with the WW treatment at the V1, V3 and R1 stages, the accumulation of dry matter in various organs was significantly reduced under water stress (PEG, PEG+ML, and PEG+MR). Compared with the WW treatment, PEG treatment in the V1 stage decreased dry matter accumulation in leaves, petioles, stems and roots by 14.3%, 60.0%, 42.9% and 32.0% (P<0.01), respectively; however, there was no significant difference between water deficit stress with melatonin application (PEG+ML and PEG+MR) and simple water deficit stress (PEG). The trend of water stress (PEG, PEG+ML, and PEG+MR) and WW treatment in stage V3 was equivalent to stage V1, but the trend between plants subjected to water stress with melatonin application (PEG+ML and PEG+MR) and to water stress alone (PEG) in stage V3 differed from those in stage V1. The accumulation of dry matter in leaves under water deficit stress treated with melatonin by leaf spraying (PEG+ML) was 9.09% (P<0.05) higher than in leaves under PEG treatment alone, and the accumulation of dry matter in roots under water deficit stress and treatment with melatonin by root irrigation (PEG+MR) was 17.3% (P<0.05) higher than in roots under PEG alone in stage V3. Compared with PEG treatment in stage R1, the accumulation of dry matter in leaves under water deficit stress and treated with melatonin by root irrigation (PEG+MR) was 9.48% (P<0.05) higher than in leaves under PEG treatment alone.

thumbnail
Table 2. Effect of exogenous melatonin on dry matter accumulation under osmotic stress.

https://doi.org/10.1371/journal.pone.0226542.t002

Effect of exogenous melatonin on antioxidant systems and MDA levels of soybean leaves under osmotic stress

The SOD activity of soybean leaves under water deficit stress (PEG, PEG+ML, and PEG+MR) was higher than in soybean leaves in the WW condition in stage V1 (Fig 2A), while the SOD activity under melatonin treatment and water deficit stress (PEG+ML and PEG+MR) was higher than under water deficit stress alone (PEG). The activity was 6.95% higher (P<0.05) in plants under PEG treatment with melatonin by root irrigation (PEG+MR) than in plants under PEG treatment alone. Stages V3 and R1 showed the same trend as Stage V1, but the activity under PEG+MR was 10.22% and 7.96% higher (P<0.05) than under PEG alone respectively, and thus, the rate of increase in stage V3 was greater compared with the other stages.

thumbnail
Fig 2. Effect of melatonin addition on the activity level of superoxide dismutase (SOD) (A), peroxidase (POD) (B), and catalase (CAT) (C) and the content of malondialdehyde (MDA) (D) under osmotic stress.

WW: half-strength Hoagland’s nutrient solution alone; WW+ML: half-strength Hoagland’s nutrient solution and 100 μmol/L melatonin by leaf; WW+MR: half-strength Hoagland’s nutrient solution and 100 μmol/L melatonin by root; PEG: half-strength Hoagland’s nutrient solution plus 15% PEG 6000 treatment; PEG+ML: half-strength Hoagland’s nutrient solution plus 15% PEG 6000 treatment and 100 μmol/L melatonin by leaf; PEG+MR: half-strength Hoagland’s nutrient solution plus 15% PEG 6000 treatment and 100 μmol/L melatonin by root. The data represent the mean ± SE of five replicate samples. * indicate significant differences according to Duncan’s multiple range test (P<0.05).

https://doi.org/10.1371/journal.pone.0226542.g002

POD and CAT jointly eliminate H2O2 produced by SOD disproportionation reactions and maintain reactive oxygen at normal levels [47]. As shown in Fig 2B and 2C, POD and CAT activities were higher in soybean leaves in stage R1 than stages V1 and V3, and the results showed essentially the same trend as reported for stages V1, V3 and R1. POD and CAT activities were higher in soybean leaves under water stress (PEG, PEG+ML, and PEG+MR) than leaves under the WW condition, while they were higher under PEG and melatonin treatment than under PEG alone. The activities of both POD and CAT showed greater increases with melatonin application by root irrigation than by leaf spraying. The POD activity levels were increased with PEG+MR treatment in V1, V3 and R1 compared with PEG alone by 12.8%, 6.21% and 5.21%(P<0.05), and CAT was increased by 9.88%, 5.49% and 6.20%(P<0.05), respectively.

To some extent, the MDA content can reflect the degree of drought stress in plants. The higher the content, the more severe is the cell membrane lipid peroxidation [48]. As shown in Fig 2D, the MDA content of soybean leaves increased rapidly after water deficit stress. The content increased by 42.1%, 9.93% and 14.5% (P<0.05), respectively, in stages V1, V3 and R1 under PEG treatment compared with WW treatment; however, melatonin treatment (PEG+ML and PEG+MR) under water deficit stress reduced the MDA content. Compared with PEG, PEG+ML reduced the content by 21.9% in stage V1 and by 9.40% in stage R1 (P<0.05), and PEG+MR reduced the content by 42.1% in stage V1, by 9.93% in stage V3 and by 14.5% in stage R1 (P<0.05). Furthermore, PEG+MR induced a decrease by 16.5% in stage V1, 16.3% in stage V3 and 4.63% in stage R1 (P<0.05) compared with PEG+ML. These findings suggest that melatonin administration by root irrigation could better regulate the activity of antioxidant enzymes in soybean leaves and reduce the degree of damage as compared to administration by leaf spraying.

Effect of exogenous melatonin on H2O2 contents and O2·- production rates of soybean leaves under osmotic stress

As shown in Fig 3(A) and 3(B), compared with the WW condition, water stress treatment significantly increased the ROS (H2O2 and O2·-) content in leaves, indicating that water deficit stress promoted the accumulation of anions in leaves, while melatonin treatment under water deficit stress significantly reduced the ROS content in leaves compared with PEG treatment alone. In stage V1, compared with PEG, PEG+ML and PEG+MR reduced H2O2 by 16.6% and 15.1% (P<0.05) (Fig 3A), and reduced O2·- by 21.6% and 23.9% (P<0.05) (Fig 3B); compared with PEG in stage V3, PEG+ML and PEG+MR reduced H2O2 by 13.5% and 28.5%, and reduced O2·- by 20.6% and 28.6% (P<0.05); compared with PEG in stage R1, PEG+ML and PEG+MR reduced H2O2 by 9.10% and 17.1%, and reduced O2·- by 9.50% and 22.0% (P<0.05); respectively. This result indicates that melatonin can reduce ROS (H2O2 and O2·-) content in leaves under water deficit stress, alleviate ROS damage to plant cells in leaves, and maintain the balance of reactive oxygen in plant leaf tissues.

thumbnail
Fig 3. Effect of melatonin addition on hydrogen peroxide (H2O2) content (A) and superoxide anion (O2·-) production rates (B) of soybean leaves under osmotic stress.

WW: half-strength Hoagland’s nutrient solution alone; WW+ML: half-strength Hoagland’s nutrient solution and 100 μmol/L melatonin by leaf; WW+MR: half-strength Hoagland’s nutrient solution and 100 μmol/L melatonin by root; PEG: half-strength Hoagland’s nutrient solution plus 15% PEG 6000 treatment; PEG+ML: half-strength Hoagland’s nutrient solution plus 15% PEG 6000 treatment and 100 μmol/L melatonin by leaf; PEG+MR: half-strength Hoagland’s nutrient solution plus 15% PEG 6000 treatment and 100 μmol/L melatonin by root. The data represent the mean ± SE of five replicate samples. * indicate significant differences according to Duncan’s multiple range test (P<0.05).

https://doi.org/10.1371/journal.pone.0226542.g003

According to the method of Awasthi [45], H2O2 and O2·- production in water deficit-stressed (PEG 6000; water potential of -0.30 MPa) and unstressed soybean seedling leaves was investigated qualitatively using DAB and NBT histochemical staining, respectively (Fig 4A and 4B). Under normal physiological water conditions, soybean seedling leaves showed low production of H2O2 and O2·-, and there was no significant different between the well-watered with melatonin application (WW+ML and WW+MR) and well-water alone (WW). However, soybean seedling leaves showed higher production under water deficit stress (PEG), while the leaves under PEG+ML and PEG+MR treatment exhibited markedly lower DAB and NBT staining, particularly compared with those under PEG treatment, which is an indication of less ROS production and less oxidative damage (Fig 4A and 4B).

thumbnail
Fig 4. Effect of melatonin addition on histochemical localization of H2O2 (A) by the DAB uptake method and of O2·- (B) by NBT staining at V3 stage during water deficit stress treatment.

WW: half-strength Hoagland’s nutrient solution alone; WW+ML: half-strength Hoagland’s nutrient solution and 100 μmol/L melatonin by leaf; WW+MR: half-strength Hoagland’s nutrient solution and 100 μmol/L melatonin by root; PEG: half-strength Hoagland’s nutrient solution plus 15% PEG 6000 treatment; PEG+ML: half-strength Hoagland’s nutrient solution plus 15% PEG 6000 treatment and 100 μmol/L melatonin by leaf; PEG+MR: half-strength Hoagland’s nutrient solution plus 15% PEG 6000 treatment and 100 μmol/L melatonin by root.

https://doi.org/10.1371/journal.pone.0226542.g004

Effect of exogenous melatonin on fluorescence characteristics of soybean leaves under osmotic stress

Fig 5 shows that water deficit stress had an impact on the chlorophyll fluorescence parameters of soybean leaves in seedlings and flowering stages. Overall, the Fv/Fm, ΦPS(II) and ETR of soybean seedlings and flowering stages demonstrated considerable downward trends under water deficit stress, while NPQ showed an upward trend. Compared with those under the WW condition, PSII under water deficit stress exhibited maximum decreases in Fv/Fm, ΦPS (II) and ETR of 47.6%, 45.5% and 50.7% (P<0.05) in stage V1, 39.2%, 44.9% and 61.2% (P<0.05) in stage V3, and 32.0%, 55.8 and 49.4% (P<0.05) in stage R1, respectively. NPQ increased by 31.1%, 32.2% and 36.0% (P<0.05) in V1, V3 and R1, respectively. The application of exogenous melatonin under PEG+ML and PEG+MR treatment increased, Fv/Fm, ΦPS(II) and ETR to obviously higher levels than those under PEG treatment; the levels in the PEG+ML and PEG+MR groups increased by 7.17% and 11.2%, 7.93% and 9.39%, 9.76% and 6.03% (P<0.05) in V1, by 9.80% and 17.4%, 8.54% and 13.8%, 10.3% and 17.1% (P<0.05) in V3, and by 6.67% and 10.7%, 6.73% and 12.0%, 5.91% and 9.17 (P<0.05) in stage R1, respectively. There was no significant difference between the PEG+ML and PEG+MR groups under water deficit stress in stage V1, while PEG+MR treatment in stage V3 resulted in significantly higher values than PEG+ML, increasing the Fv/Fm, ΦPS(II) and ETR by 6.91%, 4.82% and 6.18% (P<0.05), respectively, indicating that root application under water deficit stress has a better alleviating effect than foliar spray in stage V3.

thumbnail
Fig 5. Effect of melatonin addition on PSII primary light energy conversion efficiency under osmotic stress.

The PSII primary light energy conversion efficiency (Fv/Fm) (A), photochemical quenching coefficient (NPQ) (B), PSII actual photochemical quantum yield (ΦPSII) (C), and electron transfer rate (ETR) (D) of soybean leaves in seedlings and flowering stages during water deficit stress treatment. WW: half-strength Hoagland’s nutrient solution alone; WW+ML: half-strength Hoagland’s nutrient solution and 100 μmol/L melatonin by leaf; WW+MR: half-strength Hoagland’s nutrient solution and 100 μmol/L melatonin by root; PEG: half-strength Hoagland’s nutrient solution plus 15% PEG 6000 treatment; PEG+ML: half-strength Hoagland’s nutrient solution plus 15% PEG 6000treatment and 100 μmol/L melatonin by leaf; PEG+MR: half-strength Hoagland’s nutrient solution plus 15% PEG 6000 treatment and 100 μmol/L melatonin by root. The data represent the mean ± SE of five replicate samples. * indicate significant differences according to Duncan’s multiple range test (P<0.05).

https://doi.org/10.1371/journal.pone.0226542.g005

NPQ reflects the heat dissipation of PSII-absorbed light energy in the process of photosynthetic electron transport [49]. The NPQ trend was contrary to the trend in Fv/Fm. PEG treatment resulted in a significantly higher NPQ than the WW condition. PEG+ML and PEG+MR reduced NPQ in soybean leaves under water deficit stress. Compared with PEG alone, the NPQ of PEG+ML and PEG+MR decreased by 9.35% and 12.2% in stage V1, by 9.21% and 5.13% in stage V3, and by 5.02% and 7.31% in stage R1, respectively, and the differences reached significance (5%).

Effect of exogenous melatonin on gas exchange parameters of soybean leaves under osmotic stress

Based on Fig 6, water deficit stress significantly reduced the net photosynthetic rate (Pn) (Fig 6A); however, application of exogenous melatonin increased the Pn. PEG+ML and PEG+MR in stages V3 and R1 increased the Pn by 9.56% and 19.1%, 5.00% and 8.42% (P<0.05), respectively, compared to PEG treatment alone. Under the influence of water deficit stress, the intercellular carbon dioxide concentration (Ci) (Fig 6B) in leaves also decreased significantly. However, the content of Ci under PEG+MR treatment in stage V3 and stage R1 compared with PEG treatment alone was increased by 6.91% and 5.86% (P<0.05), respectively.

thumbnail
Fig 6. Effect of melatonin addition on the net photosynthetic rate under osmotic stress.

The net photosynthetic rate (Pn) (A), intercellular CO2 concentration (Ci) (B), transpiration rate (Tr) (C), and stomatal conductance (Gs) (D) of soybean leaves in seedlings and flowering stages during water deficit stress treatment. WW: half-strength Hoagland’s nutrient solution alone; WW+ML: half-strength Hoagland’s nutrient solution and 100 μmol/L melatonin by leaf; WW+MR: half-strength Hoagland’s nutrient solution and 100 μmol/L melatonin by root; PEG: half-strength Hoagland’s nutrient solution plus 15% PEG 6000 treatment; PEG+ML: half-strength Hoagland’s nutrient solution plus 15% PEG 6000 treatment and 100 μmol/L melatonin by leaf; PEG+MR: half-strength Hoagland’s nutrient solution plus 15% PEG 6000 treatment and 100 μmol/L melatonin by root. The data represent the mean ± SE of five replicate samples. * indicate significant differences according to Duncan’s multiple range test (P<0.05).

https://doi.org/10.1371/journal.pone.0226542.g006

Water deficit stress also significantly reduced stomatal conductance (Gs) (Fig 6D) in soybean leaves, and the decrease in Gs reduced stomatal activity, thus slowing transpiration and reducing the transpiration rate (Tr). The application of exogenous melatonin increased the Tr and Gs of soybean leaves under water deficit stress; the increases under PEG+MR treatment in stage V3 and stage R1 compared with PEG treatment were 10.7% (Fig 6C) and 17.6% (Fig 6D), and 7.64% and 26.1% (P<0.05), respectively. Thus, exogenous melatonin could alleviate the influence of water deficit stress on photosynthesis in soybean leaves by improving the Pn, Gs, Tr and Ci in soybean leaves and by maintaining strong photosynthesis to ensure the normal growth and development of soybean leaves.

Effect of exogenous melatonin on soybean yield and yield composition under osmotic stress

As shown in Table 3, soybean yield and yield components in stages V1, V3 and R1 showed no significant difference between the well-watered with melatonin application (WW+ML and WW+MR) and simple well-watered (WW) groups. However, the soybean plant height was significantly reduced under water deficit stress. Compared with the WW condition, PEG treatment reduced the height in stages V1, V3 and R1 by 27.1%, 24.1% and 19.4% (P<0.05), respectively, indicating that water deficit stress had a great impact on plant height. However, plants under PEG+MR treatment in stages V1, V3 and R1 were significantly taller than those under PEG treatment alone, while PEG+MR treatment in stages V1 and R1 resulted in significantly higher values than PEG+ML, increasing the plant height by 6.18% and 8.59% (P<0.05), respectively. The numbers of soybean pods were also affected by water deficit stress; the numbers of 2-, 3- and 4-grain pods were significantly reduced under water stress but increased after concurrent application of exogenous melatonin, such that the total grain numbers of the plants were increased. Table 3 shows that compared with PEG, PEG+ML and PEG+MR increased the total numbers of grains by 29.1% and 39.0%, 27.4% and 26.4% (P<0.05), and 9.75% and 22.3% (P<0.05) in stages V1, V3 and R1, respectively. According to the yield results (Table 3), compared with PEG, PEG+ML and PEG+MR increased the yield by 31.6% and 43.3%, 29.8% and 31.8%, and 12.6% and 25.0% (P<0.05) in the V1, V3 and R1 stages, respectively, which indicated that PEG+MR was superior to PEG+ML.

thumbnail
Table 3. Effects of external melatonin on soybean yield and yield components.

https://doi.org/10.1371/journal.pone.0226542.t003

Principal component analysis

According to the information from the GGEbiplot (Fig 7) in stages V1, V2 and R1, PC1 was explained 77.61%, 82.26% and 89.29% of the variation, and PC2 was 16.29%, 10.95% and 7.89%, showing that 93.89%, 93.21% and 97.18% of all true information was processed and that the analysis results were highly reliable.

thumbnail
Fig 7. Analysis of all indexes with different treatment by GGE biplot.

WW: half-strength Hoagland’s nutrient solution alone; WW+ML: half-strength Hoagland’s nutrient solution and 100 μmol/L melatonin by leaf; WW+MR: half-strength Hoagland’s nutrient solution and 100 μmol/L melatonin by root; PEG: half-strength Hoagland’s nutrient solution plus 15% PEG 6000 treatment; PEG+ML: half-strength Hoagland’s nutrient solution plus 15% PEG 6000 treatment and 100 μmol/L melatonin by leaf; PEG+MR: half-strength Hoagland’s nutrient solution plus 15% PEG 6000 treatment and 100 μmol/L melatonin by root. PH: plant height; SP: seed pod; SPP: seeds per plant; SOD: superoxide dismutase; POD: peroxidase; CAT: catalase; MDA: the content of malondialdehyde; H2O2: hydrogen peroxide content; O2·-: superoxide anion production raes; Fv/Fm: energy conversion efficiency; NPQ: photochemical quenching coefficient; ΦPSII: PSII actual photochemical quantum yield; ETR: electron transfer rate; Pn: net photosynthetic rate; Ci: intercellular CO2 concentration; Tr: transpiration rate; Gs: stomatal conductance.

https://doi.org/10.1371/journal.pone.0226542.g007

A polygon was formed by connecting the outermost treatment points (PEG, PEG+MR, WW+MR, WW+ML, WW) with straight lines. The three dotted lines perpendicular to each side from the center point divide the entire figure into three areas. All indicators are distributed among the three regions. We can see that four indicators NPQ, O2·−, MDA, and H2O2 are divided into PEG treatment areas in stages V1, V2 and R1. Therefore, the four indicators had the greatest influence under water deficit stress, which is also reflected in the previous analysis. In addition, the angle between these four indicators was less than 90°, indicating a significant positive correlation between them. Moreover, the angle between H2O2 and O2·− was relatively small, showing that the correlation between the two was relatively strong. Therefore, it can be preliminarily concluded that the content of NPQ, O2·−, MDA, and H2O2 of soybean leaves under PEG treatment were relatively high; thus, the degree of damage was greater than that in the other treatments.

CAT, POD and SOD were distributed in the region in which the PEG+MR treatment was the apex, indicating that the three indicators were highest under PEG+MR treatment. Those indicators were a positive correlation, because the angle is less than 90°. In addition, the angle between antioxidant enzyme activity (CAT, POD and SOD) and ROS (H2O2 and O2·−) is also less than 90°, indicating that ROS accumulation under PEG+MR treatment induced antioxidant enzyme activity and increased osmotic adjustment substances. PEG stress led to an increase in the activity of antioxidant enzymes and osmotic adjustment substances, and the activity and content of the two substances further increased in response to melatonin by root irrigation to remove excess H2O2 and O2·−, regulate cell osmotic potential, and alleviate cell damage, all of which are essentially consistent with previous findings (Figs 24). These results reflect the regulatory effect of melatonin on the antioxidant enzyme activity and osmosis-regulating contents under PEG stress.

Other indicators such as fluorescence characteristics, gas exchange parameters, two sp, three sp and yield were distributed within the well-watered region. These results showed that, compared with the other treatments, the well-watered treatment had a significant effect on the fluorescence characteristics, gas exchange parameters and yield of soybean, all of which are essentially consistent with previous findings (Figs 5 and 6, Table 3). Moreover, the angle between the fluorescence characteristics (Fv/Fm, ΦPS (II) and ETR) and gas exchange parameters (Pn, Gs, Ci and Tr) was less than 90°, indicating a positive correlation.

Discussion

Water-deficit stress is an important environmental factor inhibiting plant growth and development [50], and it is closely related to plant growth and development. PEG 6000 at a concentration of 15% had a particularly prominent inhibitory effect on soybean growth (Table 1). PEG-induced water deficit has an adverse effect on the physiological characteristics of soybean [51]. However, the application of exogenous substances can alleviate PEG-induced water deficit stress [52], and the effect varies for some physiological characteristics of soybean [53].

Drought results in a decrease in the leaf area index of soybean [54] and inhibits increases in the stem weight and root weight of soybean [55]. In our experiment, compared with the WW condition, PEG 6000 treatment simulated water deficit stress and significantly reduced the LAI, LAD and accumulation of dry matter in soybean seedlings and flowering stages. Compared with PEG alone, the application of exogenous melatonin under water deficit stress, improved the LAI and LAD. In addition, the accumulation of dry matter in roots, stems and leaves were also increased, possibly because melatonin maintained a high photosynthetic rate in the plant leaves, perhaps in particular by improving the efficiency of photosystem II (Fv/Fm). Ren [56] described the soaking of soybean seeds with melatonin that were sowed in an experimental field and then harvested at the end of a three-week growing period. The application of exogenous melatonin significantly increased the number of nodules and shoot biomass. Our research indicated that in stage V1, dry matter accumulation in different organs did not differ under water deficit stress between plants that were treated with melatonin and plants that were not treated, while the dry matter accumulation in leaves and roots in stage V3 and leaves in stage R1 was significantly higher in treated than in nontreated plants. This finding suggests that melatonin under water deficit stress has no alleviating effect on soybean morphogenesis in stage V1, but a significant effect in stage V3 and R1. This difference is possibly because the nutrients needed for organ building in stage V1 mainly came from seeds, and therefore, exogenous additives had little influence on the dry matter accumulation of soybean in this stage; in contrast, in stages V3 and R1, exogenous additives had a significant influence on the dry matter accumulation in soybean leaves and roots, consistent with the results of a study [57] on soybean nutrient requirements.

Inhibition of photosynthesis is one of the main hazards of water deficit stress [58]. Water deficit may limit photosynthesis by inducing stomatal closure (stomatal restriction), damaging the metabolic process (metabolic restriction) or both [59]. The consistent trends among intercellular carbon dioxide concentrations (Ci), net photosynthetic rates (Pn) and stomatal conductance (Gs) indicated that the number of CO2 molecules in plant cells decreased to influence the effectiveness of CO2 in the plant body. The photosynthetic efficiency of the crop plant was affected by changes in carbon metabolic approaches in the leaves, suggesting that the prevailing influencing factor in photosynthesis in this study was stomatal. In contrast, the prevailing factors in other studies have been nonstomatal [6061]. The results of the present study show that during water deficit stress, the Gs of soybean leaves was significantly reduced, and the stomatal conductance decreased such that the leaf Tr and Ci decreased. CO2 is an important raw material for photosynthesis. Reductions of photosynthetic raw materials lead to declines in the Pn in soybean leaves, partly due to of water stress-induced stomatal closure. Stomatal restriction is generally considered a short-term response to water stress that is mediated by regulation or closure of stomata, while nonstomatal effects are generally considered to involve damage to photosynthetic cells or organs due to long-term water stress [6263]. In this experiment, the Gs, Tr and Ci were improved in soybean leaves under exogenous melatonin treatment with water deficit stress compared with water deficit stress alone, as was Pn (Fig 2); these findings indicated that exogenous melatonin could regulate Gs in soybean leaves under water deficit stress to improve the Ci, increase the availability of raw materials for photosynthesis to increase the Pn, alleviate the negative effect of stomatal closure on the net photosynthesis rate under water deficit stress, and prevent irreversible damage to photosynthetic organs. Therefore, it is inferred that the factors influencing the decrease in the net photosynthetic rate under water deficit stress in this study were stomatal. Under water stress, plants accumulate the signaling hormone abscisic acid, causing rapid stomatal closure to reduce water loss, which is a short-term adaptation to water stress. In addition, the results of the gas exchange parameter test in this study also proved that PEG 6000 at a concentration of 15% did not cause serious damage to photosynthetic organs in a short period of time; furthermore, short-term damage caused by water deficit stress to the photosynthetic organs of plant leaves could be alleviated by exogenous melatonin.

Water stress directly affects the structure and activity of photosynthetic organs in plant leaves as well as the photochemical reactions and dark reactions of photosynthesis [64]. Therefore, chlorophyll fluorescence parameters are very sensitive to water deficit stress and can reflect the absorption, transmission and conversion of light energy in leaves [65]. Carcia et al. [66] suggested that water stress has no effect on Fv/Fm. In contrast, in our experiment, the Fv/Fm of soybean under water stress decreased to different degrees (Fig 3). Zhou et al. [67] indicated that the Fv/Fm reflects the health of the PSII system, suggesting that water stress affects the initial PSII light energy conversion rate and potentially affects PSII activity in soybean leaves to reduce photosynthesis. ΦPS(II) is the actual photosynthetic efficiency of PSII [68]. In this study, comparison with gas exchange parameters suggested a positive correlation between ΦPS(II) and net photosynthetic efficiency (Pn) to a certain extent. ΦPS(II) decreases mainly because of the photoelectric transfer rate (ETR) [69]. In this experiment, a significant reduction in ETR occurred under water stress, reducing the energy distribution and electron transfer between PSI and PSII in soybean leaves and further reducing the effectiveness of the plants to absorb light energy to influence soybean photosynthesis. NPQ reflects the heat consumption capacity and apparent light protection function of plants and can protect the photosynthetic organs of soybean seedlings [70]. In this study, water stress treatment significantly increased NPQ in soybean leaves, showing that under adverse conditions, soybean seedling consume light energy in the form of heat to reduce irreversible damage to soybean leaves’ PSII systems. At present, there are many research reports about exogenous hormone alleviation of plant damage under adverse conditions, such as high temperatures in rice [71] and drought in corn [72] and apple [73]. Studies have shown that under water stress conditions, exogenous hormones improve the Fv/Fm, ΦPS(II) and ETR and reduce NPQ, consistent with the results of our research. In this experiment, the application of exogenous melatonin under water stress improved the Fv/Fm, ETR and ΦPS(II) in soybean leaves, enhanced the photosynthetic electron conduction efficiency of soybean leaves under water stress, and improved the capability of leaves to capture light energy. NPQ was significantly reduced. These findings suggest that exogenous melatonin enhances the electron capture efficiency in the reaction centers of soybean leaves and maintains a high electron transfer rate. With these changes, absorbed light energy can be better used for photosynthetic electron transfer, reducing excess light energy dissipation in the form of heat, enhancing photorespiration, and protecting photosynthetic mechanisms.

Many studies have shown that water stress can lead to the disturbance of antioxidant enzyme activity in plants [74]. Excessive ROS generated under water stress significantly increase MDA levels and osmotic conductivity in plant cells, and ROS levels are considered markers of the degree of membrane lipid peroxidation [75]. Membrane lipid peroxidation in biological membranes impairs membrane function, reduces fluidity, and reduces membrane binding receptor and enzyme activity [76]. In this study, antioxidant enzyme activity and MDA levels in soybean leaves were increased under water stress, indicating that excessive ROS production was induced under water stress, and the accumulation of MDA increased the degree of membrane lipid peroxidation (Fig 2), consistent with the results of previous studies on soybean and corn [77]. Under water stress, plant cells initiate various stress responses to cope with the effects of water shortage. Melatonin is a high-intensity free radical scavenger with spectral properties that can easily pass through the cell membrane and enter the cell body to fundamentally maintain the antioxidant activity of a cell [77]. Therefore, melatonin can directly scavenge excess ROS [78], and it can also promote the accumulation of other antioxidants and antioxidant enzymes to indirectly remove ROS [79]. In our study, the application of exogenous melatonin under water stress significantly increased the activity of antioxidant enzymes, reduced the accumulation of MDA, and decreased the content of ROS (H2O2 and O2·-) in leaves (Fig 3), thus regulating the oxidative balance of soybean, maintaining intracellular ROS concentrations, protecting against ROS accumulation, and reducing the damage caused by water stress on membrane lipid peroxidation. The application of melatonin by root irrigation generated the best effect.

In addition, the application of exogenous melatonin increased the plant height, pod numbers and total grain numbers of soybean under water stress to increase the soybean yield under water stress, consistent with the results of a study in which the seeds were coated with melatonin [80]. The yield was slightly higher in stage V1 than stages V3 and R1, indicating that the effect of water stress on soybean yield was greater in stages V3 and R1 than stage V1; the reason for this difference may be that nutrients in vegetative organs in stage V1 derive from seeds, while exogenous substances have no effect on the morphogenesis of vegetative organs. The applied melatonin entered the soybean plant to maintain the oxidative balance in cells and removed excessive ROS to better alleviate the impact of external water stress on soybean in the V1 stage compared with the V3 and R1 stages. Nutrients in the V3 and R1 stages are absorbed from the external environment and obtained through nitrogen fixation during root nodulation; these processes are greatly affected by external water stress, thus affecting soybean morphogenesis and significantly reducing the soybean yield in V3 and R1 compared with V1, consistent with the results of previous studies with Ren [56] in which the seeds were soaked in melatonin. Based on the results, it can conclude that the application of melatonin by root irrigation contributes more to the improvement of soybean yield than application by foliar spraying.

The results show that melatonin increases plant growth, seed production, and abiotic stress tolerance in soybean plants, possibly through the enhancement of photosynthesis, carbohydrate metabolism, and antioxidative actions. These discoveries provide a basis for the potential application of exogenous melatonin in increasing soybean production and warrant further studies to dissect the molecular mechanism(s) that regulate melatonin-mediated water-deficient stress in soybean plants.

Conclusion

In summary, the application of exogenous melatonin under water stress in the V1 and V3 stages of soybean increases the leaf area index and leaf area duration, increases the dry matter accumulation, enhances antioxidant enzyme (SOD, POD, and CAT) activity, reduces ROS (H2O2 and O2·-) levels, reduces MDA content, improves the Fv/Fm, ETR, ΦPS (II), NPQ and other chlorophyll fluorescence characteristics, and improves the photosynthetic rate in soybean leaves to enhance photosynthesis and ultimately improve the soybean yield under water stress. Application of melatonin by root irrigation is better than application by foliar spraying. This is the first report of the application of exogenous melatonin in soybean seedlings and flowering stages during induced drought stress and the first analysis to use an experimental approach investigating leaf chloroplast morphology and photosynthesis.

Supporting information

S1 File. Supporting information data.

This file contains date including antioxidant activity, MDA levels, H2O2 contents, O2·- production rates, fluorescence characteristics, gas exchange parameters, principal component analysis data.

https://doi.org/10.1371/journal.pone.0226542.s001

(XLSX)

Acknowledgments

The sequencing and analysis were conducted at the Key Laboratory of National Coarse Cereals Engineering Research Center, Ministry of Education (at Heilongjiang Bayi Agricultural University) and were financially supported by the National Natural Science Foundation of China (No. 31801303), the China Postdoctoral Science Foundation (No. 2016M591561), the Natural Science Foundation of Heilongjiang Province of China (No. C2016042), the Doctoral Scientific Research Start-up Foundation of Heilongjiang Bayi Agricultural University of China (No. XYB2014-04), and Key Scientific Research Projects of Heilongjiang Farms and Land Reclamation administration (No. HKKY190206-1).

References

  1. 1. Jin ZN, Ainsworth EA, Leakey ADB. Lobell DB. Increasing drought and diminishing benefits of elevated carbon dioxide for soybean yields across the US Midwest. Global change biology. 2018; 24: e522–e533. https://doi.org/10.1111/gcb.13946 pmid:29110424
  2. 2. Gray SB, Dermody O, Klein SP, Locke AM, McGrath JM, Paul RE. Intensifying drought eliminates the expected benefits of elevated carbon dioxide for soybean. Nature Plants. 2016; 2:16132. https://doi.org/10.1038/NPLANTS.2016.132 pmid:27595230
  3. 3. Meng JF, Xu TF, Wang ZZ, Fang YL, Xi ZM, Zhang ZW. The ameliorative effects of exogenous melatonin on grape cuttings under water-deficient stress: antioxidant metabolites, leaf anatomy, and chloroplast morphology. Journal of Pineal Research. 2014; 57: 200–212. https://doi.org/10.1111/jpi.12159 pmid:25039750
  4. 4. Golldack D, Li C, Mohan H, Probst N. Tolerance to drought and salt stress in plants: Unra- veling the signaling networks. Front in Plant Sci. 2014; 5: 151. https://doi.org/10.3389/fpls.2014.00151
  5. 5. Graham N, Amna M, Christine HF. The Roles of Reactive Oxygen Metabolism in Drought: Not So Cut and Dried. Plant Physiology. 2014; 164: 1636–1648. https://doi.org/10.1104/pp.113.233478 pmid:24715539
  6. 6. Xu LX, Yu JJ, Han LB, Huang B. Photosynthetic enzyme activities and gene expression associated with drought tolerance and post-drought recovery in Kentucky bluegrass. Environmental and Experimental Botany. 2013; 89: 28–35. https://doi.org/10.1016/j.envexpbot.2012.12.001
  7. 7. Cominelli E, Conti L, Tonelli C, Galbiati M. Challenges and perspectives to improve crop drought and salinity tolerance. New Biotechnology. 2013; 4: 355–361. https://doi.org/10.1016/j.nbt.2012.11.001
  8. 8. Zhang C, Zhan DX, Luo HH, Zhang YL, Zhang WF. Photorespiration and photoinhibition in the bracts of cotton under water stress. Photosynthetica. 2016; 1: 12–18. https://doi.org/10.1007/s11099-015-0139-9
  9. 9. Bollig C, Feller U. Impacts of drought stress on water relations and carbon assimilation in grassland species at different altitudes. Agric. Ecosyst. Environ. 2014; 188: 212–220. https://doi.org/10.1016/j.agee.2014.02.034
  10. 10. Osakabe Y, Osakabe K, Shinozaki K, Tran LSP. Response of plants to water stress. Frontiers in Plant Science. 2014; 5: 86. https://doi.org/10.3389/fpls.2014.00086 pmid:24659993
  11. 11. De Souza TC, Magalhães PC, de Castro EM, de Albuquerque PEP, Marabesi MA. The influence of ABA on water relation, photosynthesis parameters, and chlorophyll fluorescence under drought conditions in two maize hybrids with contrasting drought resistance. Acta Physiol Plant. 2013; 35:515–527. https://doi.org/10.1007/s11738-012-1093-9
  12. 12. Wu YY, Tian YD, Liu LX, Xing XH, Jiang HQ, Xing H, et al. Effects of Triadimefon on physiological characteristics and yield of soybean under drought and rewatering at flowering stage. Journal of Nuclear Agricultural Sciences. 2013; 27: 1749–1755.
  13. 13. Xing XH, Jiang HQ, Zhou Q, Xing H, Jiang H, Wang S. Improved drought tolerance by early IAA- and ABA-dependent H2O2 accumulation induced by a-naphthaleneacetic acid in soybean plants. Plant Growth Regulation. 2016; 3: 303–314. https://doi.org/10.1007/s10725-016-0167-x
  14. 14. Wang W, Wang C, Pan D, Zhang Y, Luo B, Ji J. Effects of drought stress on photosynthesis and chlorophyll fluorescence images of soybean (Glycine max) seedlings. International Journal of Agricultural and Biological Engineering. 2018; 11: 196–201.
  15. 15. Wang X, Khodadadi E, Fakheri B, Komatsu S. Organ-specific proteomics of soybean seedlings under flooding and drought stresses. Journal of Proteomics. 2017; 162: 62–72. https://doi.org/10.1016/j.jprot.2017.04.012 pmid:28435105
  16. 16. Chartzoulakis K, Patakas A, Kofidis G, Bosabalidis A, Nastou A. Water stress affects leaf anatomy, gas exchange, water relations and growth of two avocado cultivars. Scientia horticulturae. 2002; 95: 39–50. https://doi.org/10.1016/S0304-4238(02)00016-X
  17. 17. Mathobo R, Marais D, Steyn JM. The effect of drought stress on yield, leaf gaseous exchange and chlorophyll fluorescence of dry beans (Phaseolus vulgaris L.). Agricultural Water Management. 2017; 180: 118–125. https://doi.org/10.1016/j.agwat.2016.11.005
  18. 18. Mo C, Kim MS, Kim G, Cheong EJ, Yang J, Lim J. Detecting Drought Stress in Soybean Plants Using Hyperspectral Fluorescence Imaging. Journal of Biosystems Engineering. 2015; 40: 335–344.
  19. 19. Yao CJ, Guo SM, Ma YC, Lai XL, Yang XH. Effect of drought stress on characteristics of photosynthesis and chlorophyll fluorescence of four species of Cassia. Pratacultural Science. 2017; 34: 1880–1888. https://doi.org/10.11829/j.issn.1001-0629.2016-0588
  20. 20. Li YB, Song H, Zhou L, Xu Z, Zhou G. Tracking chlorophyll fluorescence as an indicator of drought and rewatering across the entire leaf lifespan in a maize field. Agricultural Water Management. 2019; 211: 190–201. https://doi.org/10.1016/j.agwat.2018.09.050
  21. 21. Wu X L, Liang H Y, Yang F, Liu WG, She Y.H, Yang WY. Comprehensive evaluation and screening identification indexes of shade tolerance at seedling in soybean. Scientia Agricultura Sinica. 2015; 48: 2497–2507. https://doi.org/10.3864/j.issn.0578-1752.2015.13.002
  22. 22. Gao W, Zhang Y, Feng Z, Bai Q, He J, Wang Y. Effects of Melatonin on Antioxidant Capacity in Naked Oat Seedlings under Drought Stress. Molecules. 2018; 23, 1580. https://doi.org/10.3390/molecules23071580
  23. 23. Mether , Shivakrishna P, Reddy KA, Rao DM. Effect of PEG-6000 imposed drought stress on RNA content, relative water content (RWC), and chlorophyll content in peanut leaves and roots. Saudi Journal of Biological Sciences. 2018; 25:285–289. https://doi.org/10.1016/j.sjbs.2017.04.008 pmid:29472779
  24. 24. Emmerich WE, Hardegree SP. Polyethylene glycol solution contact effects on seed germinat- ion. Agronomy Journal. 1990; 6: 1103–1107. https://doi.org/10.2134/agronj1990.00021962008200060015x
  25. 25. He SY, Qin B, Zhang MC, Jin XJ, Wang MX, Ren CY, Zhang YX. Effects of Exogenous Melatonin on Antioxidant Properties and Yield of Soybean Seedling Under Water Stress. Soybean Science, 2019, 38: 407–412. https://doi.org/10.11861/j.issn.1000-9841.2019.03.0407
  26. 26. Arnao MB, Hernández RJ. Functions of melatonin in plants: a review. Journal of Pineal Research. 2015; 2: 133–150. https://doi.org/10.1111/jpi.12253
  27. 27. Dubbels R, Reiter RJ, Klenke E, Goebel A, Schnakenberg E, Ehlers C, et al. Melatonin in edible plants identified by radioim-munoassay and by high performance liquid chromatography-mass spectrometry. Journal of pineal research, 1995; 18, 28–31. https://doi.org/10.1111/j.1600-079X.1995.tb00136.x pmid:7776176
  28. 28. Arnao MB. Phytomelatonin: discovery, content, and role in plants. Advances in Botany. 2014; 2: 1–12. http://doi.org/10.1155/2014/815769
  29. 29. Murch SJ, KrishnaRaj S, Saxena P. Tryptophan is a precursor for melatonin and serotonin biosynthesis in vitro regenerated St. John’s wort (Hypericum perforatumL. cv. Anthos) plants. Plant Cell Reports. 2000; 19: 698–704. https://doi.org/10.1007/s002990000206 pmid:30754808
  30. 30. Manchester LC, Tan DX, Reiter RJ, Park W, Monis K, Qi W. High levels of melatonin in the seeds of edible plants: Possible function in germ tissue protection. Life Sciences. 2000; 25: 3023–3029. https://doi.org/10.1016/S0024-3205(00)00896-1
  31. 31. Wang P, Sun X, Li C, Wei ZW, Liang D, Ma FW. Long-term exogenous application of melatonin delays drought-induced leaf senescence in apple. Journal of Pineal Research. 2013; 54: 292–302. https://doi.org/10.1111/jpi.12017 pmid:23106234
  32. 32. Chen YE, Mao JJ, Sun LQ, Huang B, Ding CB, Gu Y, et al. Exogenous melatonin enhances salt stress tolerance in maize seedlings by improving antioxidant and photosynthetic capacity. Physiologia Plantarum. 2018; 164: 349–363. https://doi.org/10.1111/ppl.12737 pmid:29633289
  33. 33. Turk H, Erdal S, Genisel M, Atici O, Demir Y, Yanmis D. The regulatory effect of melatonin on physiological, biochemical and molecular parameters in cold-stressed wheat seedling. Plant Growth Regulation. 2014; 74: 139–152. https://doi.org/10.1007/s10725-014-9905-0
  34. 34. Wang LY, Liu JL, Wang WX, Sun Y. Exogenous melatonin improves growth and photosynthetic capacity of cucumber under salinity-induced stress. Photosynthetica. 2016; 54: 19–27. https://doi.org/10.1007/s11099-015-0140-3
  35. 35. Zhang N, Zhao B, Zhang HJ, Weeda S, Yang C, Yang ZC, et al. Melatonin promotes water-stress tolerance, lateral root formation, and seed germination in cucumber. Journal of Pineal Research. 2013; 54:15–23. https://doi.org/10.1111/j.1600-079X.2012.01015.x pmid:22747917
  36. 36. Yang XL, Xu H, Li D, Gao X, Li TL, Wang R. Effect of melatonin priming on photosynthetic capacity of tomato leaves under low-temperature stress. Photosynthetica. 2018; 56: 884–892. https://doi.org/10.1007/s11099-017-0748-6
  37. 37. Arnao M B, Hernàndez RJ. Protective effect of melatonin against chlorophyll degradation during the senescence of barley leaves. Journal of Pineal Research. 2009; 46: 58–63. https://doi.org/10.1111/j.1600-079X.2008.00625.x pmid:18691358
  38. 38. Jiang XW, Li HQ, Song XY. Seed priming with melatonin effects on seed germination and seedling growth in maize under salinity stress. Pak. J. Bot. 2006; 48: 1345–1352.
  39. 39. Liang C, Zheng G, Li W, Wang Y, Hu B, Wang H, et al. Melatonin delays leaf senescence and enhances salt stress tolerance in rice. Journal of Pineal Research. 2015; 59: 91–101. https://doi.org/10.1111/jpi.12243 pmid:25912474
  40. 40. Beauchamp CO, Fridovich I. Isozymes of superoxide dismutase from wheat germ. Biochim-ica et Biophysica Acta (BBA)-Protein Structure. 1973; 317: 50–64. https://doi.org/10.1016/0005-2795(73)90198-0
  41. 41. Rao M V, Paliyath G, Ormrod DP. Ultraviolet-B- and Ozone-Induced Biochemical Changes in Antioxidant Enzymes of Arabidopsis thaliana. Plant Physiol. 1996; 10: 125–136. https://doi.org/10.1104/pp.110.1.125
  42. 42. Hamurcu M, Sekmen AH, Turkan I, Gezgin S, Demiral T, Bell RW. Induced anti-oxidant activity in soybean alleviates oxidative stress under moderate boron toxicity. Plant Growth Regulation. 2013; 70: 217–226. https://doi.org/10.1007/s10725-013-9793-8
  43. 43. Draper HH, Squires EJ, Mahmoodi H, Wu J, Agarwal S, Hadley M. A comparative evaluation of thiobarbituric acid methods for the determination of malondialdehyde in biological materials. Free Radical Biology and Medicine. 1993; 4: 353–363. https://doi.org/10.1016/0891-5849(93)90035-S
  44. 44. Ke D, Sun G, Wang Z. Effects of superoxide radicals on ACC synthase activity in chilling-stressed etiolated mungbean seedlings. Plant Growth Regulation. 2007; 51: 83–91. https://doi.org/10.1007/s10725-006-9150-2
  45. 45. Awasthi JP, Saha B, Regon P, Sahoo S, Chowra U, Pradhan A, et al. Morphophysiological analysis of tolerance to aluminum toxicity in rice varieties of North East India. PLoS ONE. 2017, 12: e0176357. https://doi.org/10.1371/journal.pone.0176357 pmid:28448589
  46. 46. Baker NR. Chlorophyll fluorescence: A probe of photosynthesis in vivo. Annual Review of Plant Biology. 2008; 59: 89–113. https://doi.org/10.1146/annurev.arplant.59.032607.092759 pmid:18444897
  47. 47. Liu Z, Cai JS, Li JJ, Lu GY, Li CS, Fu GP, et al. Exogenous application of a low concentration of melatonin enhances salt tolerance in rapeseed (Brassica napus L.) seedlings. Journal of Integrative Agriculture. 2018; 17: 328–335. https://doi.org/10.3390/ijms20030709
  48. 48. Taïbi K, Taïbi F, Abderrahim LA, Ennajah A, Belkhodja M, Mulet JM. Effect of salt stress on growth, chlorophyll content, lipid peroxidation and antioxidant defence systems in Phaseolus vulgaris L. South African Journal of Botany. 2016; 105: 306–312. https://doi.org/10.1016/j.sajb.2016.03.011
  49. 49. Hu W, Tian SB, Di Q, Duan SH, Dai K. Effects of exogenous calcium on mesophyll cell ultrastructure, gas exchange, and photosystem II in tobacco (Nicotiana tabacum Linn.) under drought stress. Photosynthetica. 2018; 56: 1204–1211. https://doi.org/10.1007/s11099-018-0822-8
  50. 50. Harb A, Krishnan A, Ambavaram MM, Pereira A. Molecular and Physiological Analysis of Drought Stress in Arabidopsis Reveals Early Responses Leading to Acclimation in Plant Growth. Plant Physiology. 2010; 154: 1254–1271. https://doi.org/10.1104/pp.110.161752 pmid:20807999
  51. 51. Hamayun M, Khan SA, Shinwari ZK, Khan AL, Ahmad N, Lee IJ. Effect of polyethylene glycol induced drought stress on physio-hormonal attributes of soybean. Pakistan Journal of Botany. 2010; 42: 977–986.
  52. 52. Sergiev I, Todorova D, Shopova E, Jankauskiene J, Jankovska-Bortkevič E, Jurkonienė S. Exogenous auxin type compounds amend PEG-induced physiological responses of pea plants. Scientia Horticulturae. 2019; 248: 200–205. https://doi.org/10.1016/j.scienta.2019.01.015
  53. 53. Tian F, Jia T, Yu B. Physiological regulation of seed soaking with soybean isoflavones on drought tolerance of Glycine max and Glycine soja. Plant Growth Regulation. 2014; 74: 229–237. https://doi.org/10.1007/s10725-014-9914-z
  54. 54. Salvatori E, Fusaro L, Manes F. Chlorophyll fluorescence for phenotyping drought-stressed trees in a mixed deciduous forest. Annali di Botanica. 2016; 6: 39–49. https://doi.org/10.4462/annbotrm-13263
  55. 55. Guo S J, Yang K M, Huo J, Zhou YH, Wang YP, Li GQ. Influence of drought on leaf photosynthetic capacity and root growth of soybeans at grain filling stage. Chinese Journal of Applied Ecology. 2015; 26: 1419–1425 pmid:26571660
  56. 56. Ren SX, Jiang GL, Rutto L. Melatonin priming enhances symbiotic nitrogen fixation in soybean, Glycine max L. Journal of Biotech Research. 2019; 10: 136–144. http://www.btsjournals.com/assets/2019v10p136-144.pdf
  57. 57. Dong Z. Soybean yield physiology. Beijing: China Agriculture Press. 2012; 75–82.
  58. 58. Kaiser W. M. Effects of water deficit on photosynthetic capacity. Physiologia Plantarum. 1987; 71: 142–149. https://doi.org/10.1111/j.1399-3054.1987.tb04631.x
  59. 59. Wallin G, Skarby L. The influence of ozone on the stomatal and non-stomatal limitation of photosynthesis in Norway spruce, Picea abies (L.) Karst. Exposed to soil moisture deficit. Trees. 1992; 6: 128–136. https://doi.org/10.1007/bf00202428
  60. 60. Escalona JM, Flexas J, Medrano H. Stomatal and non-stomatal limitations of photosynthesis under water stress in field-grown grapevines. Australian Journal of Plant Physiology. 1999; 26: 421–433. https://doi.org/10.1071/PP99019
  61. 61. Flexas J, Ribas-Carbó M, Bota J, Galmés J, Henkle M, Martínez-Cañellas S, et al. Decreased Rubisco activity during water stress is not induced by decreased relative water content but related to conditions of low stomatal conductance and chloroplast CO2 concentration. New Phytologist. 2006; 172: 73–82. https://doi.org/10.1111/j.1469-8137.2006.01794.x pmid:16945090
  62. 62. Farquhar G D, Sharkey T D. Stomatal conductance and photosynthesis. Annual Review of Plant Physiology. 1982; 33: 317–345. https://doi.org/10.1146/annurev.pp.33.060182.001533
  63. 63. Sharkey TD. Photosynthesis in intact leaves of C3 plants: Physics, physiology and rate limitations. The Botanical Review. 1985; 51: 53–105. https://doi.org/10.2307/4354049
  64. 64. Fracheboud Y, & Leipner J. The Application of Chlorophyll Fluorescence to Study Light, Temperature, and Drought Stress. Practical Applications of Chlorophyll Fluorescence in Plant Biology. 2003; 125–150. https://doi.org/10.1007/978-1-4615-0415-3-4
  65. 65. Li RH, Guo PG, Michael B, Stefania G, Salvatore C. Evaluation of Chlorophyll Content and Fluorescence Parameters as Indicators of Drought Tolerance in Barley. Agricultural Sciences in China. 2006; 5: 751–757. https://doi.org/10.1016/S1671-2927(06)60120-X
  66. 66. García-Sánchez F, Syvertsen JP, Gimeno V, Botía P, Perez-Perez JG. Responses to flooding and drought stress by two citrus rootstock seedlings with different water-use efficiency. Physiologia Plantarum, 2007; 130, 532–542. https://doi.org/10.1111/j.1399-3054.2007.00925.x
  67. 67. Zhou S X, Medlyn BE and Prentice IC. Long-term water stress leads to acclimation of drought sensitivity of photosynthetic capacity in xeric but not riparian Eucalyptus species. Annals of Botany. 2015; 117: 133–44. https://doi.org/10.1093/aob/mcv161 pmid:26493470
  68. 68. Oukarroum A, El Madidi S, Schansker G, Strasser RJ. Probing the responses of barley cultivars (Hordeum vulgare L.) by chlorophyll a fluorescence OLKJIP under drought stress and re-watering. Environmental and Experimental Botany. 2007; 60: 438–446. https://doi.org/10.1016/j.envexpbot.2007.01.002
  69. 69. Ivanov B, & Edwards G. Influence of ascorbate and the Mehler peroxidase reaction on non-photochemical quenching of chlorophyll fluorescence in maize mesophyll chloroplasts. Planta. 2000; 210: 765–774. https://doi.org/10.2307/23385907 pmid:10805448
  70. 70. Fahad Shah, Hussain Saddam, Saud Shah, Hassan S, Ihsan Z, Shah AN, et al. Exogenously Applied Plant Growth Regulators Enhance the Morpho-Physiological Growth and Yield of Rice under High Temperature. Frontiers in Plant Science. 2016; 7: 1–13. https://doi.org/10.3389/fpls.2016.01250
  71. 71. Akter N, Islam MR, Karim MA, Hossain T. Alleviation of drought stress in maize by exogenous application of gibberellic acid and cytokinin. Journal of Crop Science and Biotechnology. 2014; 17: 41–48. https://doi.org/10.1007/s12892-013-0117-3
  72. 72. Zhang L, Gao M, Zhang L, Li B, Han M, Alva AK, et al. Role of exogenous glycine betaine and humic acid in mitigating drought stress-induced adverse effects in Malus robust a seedling. Turkish Journal of Botany. 2013; 37: 920–929. https://doi.org/10.3906/bot-1212-21
  73. 73. Bai LP, Sui FG, Ge TD, Sun ZH, Lu YY, Zhou GS. Effect of Soil Drought Stress on Leaf Water Status, Membrane Permeability and Enzymatic Antioxidant System of Maize. Pedosphere. 2006; 16: 326–332. https://doi.org/10.1016/S1002-0160(06)60059-3
  74. 74. Gharibi S, Tabatabaei BES, Saeidi G, Goli SAH. Effect of Drought Stress on Total Phenolic, Lipid Peroxidation, and Antioxidant Activity of Achillea Species. Applied Biochemistry and Biotechnology. 2016; 178: 796–809. https://doi.org/10.1007/s12010-015-1909-3 pmid:26541161
  75. 75. Hou Q, Ufer G, & Bartels D. Lipid signalling in plant responses to abiotic stress. Plant, Cell & Environment. 2016; 39: 1029–1048. https://doi.org/10.1111/pce.12666
  76. 76. Devi MA, & Giridhar P. Variations in Physiological Response, Lipid Peroxidation, Antioxidant Enzyme Activities, Proline and Isoflavones Content in Soybean Varieties Subjected to Drought Stress. Proceedings of the National Academy of Sciences, India Section B: Biological Sciences. 2015; 85: 35–44. https://doi.org/10.1007/s40011-013-0244-0
  77. 77. Campos CN, Ávila RG, de Souza KRD, Azevedo LM, Alves JD. Melatonin reduces oxidative stress and promotes drought tolerance in young Coffea arabica L. plants. Agricultural Water Management. 2019; 211: 37–47. https://doi.org/10.1016/j.agwat.2018.09.025
  78. 78. Reiter RJ, Tan DX, Burkhardt S, Manchester LC. Melatonin in plants. Nutrition Reviews. 2001; 59: 286–290. https://doi.org/10.1111/j.1753-4887.2001.tb07018.x pmid:11570431
  79. 79. Zhang N, Zhang H, Yang R, Huang Y, Guo Y. Advances in Melatonin and Its Functions in Plants. Agricul- tural Science & Technology. 2012; 13: 1833–1837. https://doi.org/10.3969/j.issn.1009-4229-B.2012.09.006
  80. 80. Wei W, Li QT, Chu YN, Reiter RJ, Yu XM, Zhu DH, et al. Melatonin enhances plant growth and abiotic stress tolerance in soybean plants. Journal of Experimental Botany. 2015; 66: 695–707. https://doi.org/10.1093/jxb/eru392 pmid:25297548