Abstract

Drought is a major abiotic stress that negatively impacts plant growth and crop production. Among various techniques used to alleviate drought stress in plants, nanoparticle application is considered to be effective and promising. In this study, the responses of plants treated with iron, copper, cobalt, and zinc oxide nanoparticles (NPs) were analyzed in soybean under drought-induced conditions. The obtained results indicated that these metal-based NPs supported the drought tolerance of NP-treated plants. The desired physiological traits, viz., relative water content, drought tolerance index, and biomass reduction rate, were significantly improved, especially in iron NP-treated plants. At the molecular level, quantitative PCR analysis of several drought-responsive genes revealed a gene-, tissue-, and NP-dependent upregulation of gene expression. Iron NP treatment promoted the expression of all tested genes in roots; additionally, the expression of three drought-responsive genes increased in leaves of all NP-treated plants, while the expression of GmERD1 (Early Responsive to Dehydration 1) was induced in both roots and shoots under the four NP treatments tested. Our findings suggest that NP application can improve drought tolerance of soybean plants by triggering drought-associated gene expression.

1. Introduction

Among the top 10 most widely grown crop species, soybean (Glycine max (L.) Merrill) is a globally important crop used for both food and animal feed [1]. Unfortunately, it is also among the most drought-susceptible crop plants [2], with up to 40% of annual yield losses due to water shortage. The problem is even more serious in tropical countries, such as Vietnam, where the prevalence and duration of drought have increased over the years. To overcome this challenge, alternative breeding and biotechnological approaches, such as micronutrient modification, have been extensively studied [36]. The supplementation of mineral microelements, such as iron, copper, cobalt, magnesium, manganese, nickel, and zinc, is known to increase the crop yield, especially under adverse environmental conditions [7, 8].

The rapid development of nanotechnology has allowed metal-based nanoparticles (NPs) of 1–100 nm in size and in extremely low quantity to be tested as substitute plant mineral nutrients and stimulants. The application of NPs enhances plant responses to drought stress compared to conventional bulk fertilization forms [9]. For example, the application of ZnO NPs increased soybean seed germination percentage under water stress [10]. Similarly, sorghum productivity and nitrogen acquisition improved remarkably with ZnO NP application to the soil [11]. Furthermore, the addition of ZnO NPs to the culture medium promoted somatic embryo formation and plant regeneration and increased stress tolerance in in vitro-grown banana plants [12]. In turn, Cu and Zn NPs effectively attenuated drought effects on wheat plants by increasing antioxidant enzyme activities and relative water content, reducing thiobarbituric acid reactive substance (TBARS) accumulation, and stabilizing photosynthetic pigment content in leaves [13]. The use of zerovalent Fe NPs was found to support normal drought sensitivity maintenance, plasma membrane Hþ-ATPase activation, and stomatal opening, as well as chlorophyll and plant biomass content, while CO2 assimilation of Arabidopsis thaliana increased during water stress [14].

Although the possible roles of NPs in boosting plant defense responses to abiotic stress, including water stress, have been extensively studied, many aspects remain to be elucidated. When analyzing related researches and their own work, Khan et al. [15] hypothesized that NPs may mediate plant adaption and tolerance to stress through the activation of defense systems and stress-related gene expression. Plants respond to environment stress factors by activating various transcellular membrane sensors, particularly Ca2+ channels and Ca2+-binding proteins [16]. Subsequent downstream events lead to changes in gene expression and ultimately to plant adaptation to stress. Among stress-responsive genes, RD20A (Response to Desiccation 20A) and ERD1 (Early Responsive to Dehydration 1) are among the most responsive stress marker genes in plants, including soybean [17, 18]. Regarding transcription factor (TF) genes, families such as bZIP (basic leucine zipper), DREB (dehydration-responsive element-binding factor), NAC (NAM (no apical meristem), ATAF (Arabidopsis transcription activation factor), and CUC (cup-shaped cotyledon)), WRKY (composed of a conserved WRKYGQK motif), MYB (myeloblastosis), ERF (ethylene response factor), and ABF/AREB (abscisic acid-responsive element-binding factor) are associated with plant adaptation to environmental stress [19, 20]. The overexpression of GmFDL19, of the bZIP family, significantly improved transgenic soybean tolerance to water shortage [21]. Similarly, GmDREB2 overexpression led to enhanced drought and high-salt stress tolerance of transgenic Arabidopsis plants [22]. Recently, GmWRKY27, GmMYB118, and GmMYB174 were found to be strongly upregulated in soybean during abiotic stress, including drought [23, 24]. Furthermore, GmNAC11 regulated DREB1A and other stress-related genes (COR15A, ERF5, ERD11, etc.) in transgenic soybean [25].

The contribution of metal-based NPs to the improvement of plant adaptation to abiotic stress has been extensively demonstrated[2632]; however, their role in inducing stress-related gene expression has not been evaluated yet. Therefore, in the present study, our objective was to evaluate physiological and molecular responses of soybean plants treated with Fe, Cu, Co, and ZnO NPs under water stress. Physiological indicators studied included relative water content (RWC), drought tolerance index (DTI), and biomass reduction rate. As for the expression of drought tolerance marker genes, our analysis included GmRD20A, GmDREB2, GmERD1, GmFDL19, GmNAC11, GmWRKY27, GmMYB118, and GmMYB174. A model explaining the likely mode of action of NPs in triggering plant tolerance to water stress was proposed based on the results reported herein and the literature available.

2. Materials and Methods

2.1. Sample Preparation and Drought Treatments

Seeds of the local soybean (Glycine max (L.) Merrill) cultivar DT26 harvested in the winter of 2017 were provided by the Legumes Research and Development Center of the Food Crops Research Institute at the Vietnam Academy of Agriculture Sciences. Healthy, uniformly sized seeds were selected for NP seed treatment. NPs, including Fe, Cu, Co, and ZnO with and size of 20–60 nm, provided by the Institute of Environmental Technology, VAST [33], were suspended in distilled deionized water separately to generate 50 mg/L of Fe, ZnO, and Cu NPs and 0.05 mg/L Co NP; subsequently, each suspension was sonicated with Elmasonic S 100/H (Elma Schmidbauer GmbH, Singen, Germany) at 200 W and 37 kHz for 30 min and immediately used for seed treatment. One hundred and fifty seeds (approximately 30 g) were soaked in 10 mL NP suspensions for 30 min. The control batches of seeds were soaked in distilled deionized water. The concentrations of each NP type were selected based on our previous experimental results [34].

Drought treatment was carried out as described by Thu et al. [35]. Briefly, in all experiments, plants were grown in a premixed, standard clean black soil (Namix, Ho Chi Minh City, Vietnam) rich in nutrients and germ-free. Control (untreated) and NP-treated seeds were sown in PVC sheets rolled into tubular shape (80 cm in height, 10 cm in diameter, 1 seed/tube) for maximum axial root growth. Plants were grown in a net house under natural light, 28–30°C ambient temperature, and 60–70% relative humidity; they were watered every day until the V2-V3 (vegetative second and third nodes) stage (12 days) and then divided into well-watered control and induced drought stress groups for which watering was completely stopped for the following 15 days. Soil moisture content (SMC) was monitored at 5-day intervals () using the HydroSense II system (Campbell Scientific, Inc., Logan, UT, USA) to maintain approximately 65–70% SMC for the control group while for the induced drought stress treatment groups, SWC was allowed to gradually decrease to 30–40% by the end of the experimental period.

After drought treatment, plants ( for each treatment) were gently removed and the lengths of shoots and roots of each individual plant were recorded. Roots and shoots were then dried separately at 65°C for 48 h to determine dry weight (DW). The specific root length (SRL) index was calculated based on root length and dry mass.

2.2. Physiological Parameter Assays
2.2.1. Relative Water Content

Determination of relative water content (RWC) was carried out 27 days after planting as described by Hossain et al. [36]. The aerial parts of the plants () from NP-treated and nontreated groups were collected and weighed immediately for calculating fresh weight (FW). Turgid weight (TW) of samples was estimated by soaking samples in deionized water for 24 h at , placing them on absorbent paper to remove the surface water, and weighing them. Finally, samples were dried at 65°C for 48 h for dry weight (DW) measurements. RWC was calculated as per the following formula:

2.2.2. Drought Tolerance Index

Examination of the drought tolerance index (DTI) was performed as described by Thu et al. [35]. Seeds of each NP treatment and control were germinated in plastic pots (25 cm in height and 30 cm in diameter, five seeds/pot; ) and watered regularly for 12 days. Drought treatment was carried out as described above. During the 15 days of drought treatment duration, the number of nonwithered plants was recorded every two days. When SMC decreased to 30% or lower and the number of nonwithered plants reached more than 50% of all plants in the group, plants were watered again for the following 15 days. The number of plants showing recovery was determined at 2-day intervals. DTI of each NP treatment was calculated as per the following formula: where is the percentage of nonwilted plants after day(s) of drought treatment, is the percentage of plants showing recovery after day(s) of rewatering, , and the number of equal inner angles () is 16.

2.2.3. Biomass Reduction Rate

Plant biomass reduction under drought was calculated according to Nguyen et al. [37]. Each plant was grown in a plastic pot containing 600 g soil and 250 mL water. Experiments involving drought treatment were performed as described above. After 15 days, plants were allowed to recover by rewatering for the following seven days. The entire shoot of each plant was collected and dried at 65°C for 48 h to determine plant dry weight ( per treatment). The biomass reduction rate (%) was calculated as per the following formula: where WW is the dry weight of plants under well-watered conditions and DR is the dry weight of plants showing recovery from drought treatment.

2.3. Gene Expression Analyses

To analyze the effects of NP treatment on the expression of target genes under drought conditions, roots and leaves of plants were collected for RNA isolation. Total RNA was isolated and purified using the GeneJET Plant RNA Purification Mini Kit (Thermo Fisher Scientific, Waltham, MA, USA), and DNA was removed by DNaseI using the RapidOut DNA removal kit (Thermo Fisher Scientific, Waltham, MA, USA). RNA concentrations were determined using the NanoDrop™ 2000/2000c Spectrophotometers (Thermo Fisher Scientific, Waltham, MA, USA). cDNA synthesis was performed using 1 μg total RNA and the Maxima First Strand cDNA Synthesis Kit (Thermo Fisher Scientific, Waltham, MA, USA).

Primer pairs for drought-responsive genes GmDREB2, GmNAC11, GmWRKY27, GmMYB118, and GmMYB174 were designed using the Primer3Plus software (http://frodo.wi.mit.edu/primer3/), while primer sequences for the other genes (GmRD20A, GmERD1, and GmFDL19) were designed as reported previously [18, 21]. Primer specificity was confirmed by BLAST (http://blast.ncbi.nlm.nih.gov/Blast.cgi), followed by analyzing the melting curves and amplicon fragments. The sequences of primers used are listed in Table 1. The GmFbox gene was selected as a reference gene [38].

Reverse transcription- (RT-) qPCR was performed in 96-well plates on the Eppendorf Realplex 4 Mastercycler (Eppendorf AG, Hamburg, Germany) using the SYBR Green PCR Master Mix (Thermo Fisher Scientific, Waltham, MA, USA). RT-qPCR reaction was carried out as described by Thao et al. [39]. Relative expression levels of each gene were calculated using the 2∆∆Ct method [40] with expression normalized against the internal reference gene, GmFbox. Three independent biological samples were used from each experimental group.

2.4. Statistical Analysis

Data were statistically analyzed according to Student’s -test (significant differences with are indicated by asterisks) or one-way analysis of variance (ANOVA) with significant differences among treatments at determined according to Tukey’s test.

3. Results

3.1. Effect of NP Treatments on the Development of the Soybean DT26 Cultivar under Drought Conditions

The performance of NP-treated plants was analyzed under induced drought conditions to understand the possible benefits of NP application to the enhancement of soybean tolerance to a water deficit. The experiments were conducted under strictly monitored conditions to maintain soil moisture content (% SMC) between 60 and 70% in the control treatment groups and to gradually reduce SMC to 30–40% by the last day of drought treatment. The development of shoots and roots of these plants was analyzed after 15 days of induced drought treatment.

Shoot morphology analysis (Figures 1(a) and 1(b)) indicated that NP treatment promoted shoot development in both well-watered and drought-stressed plants at the vegetative stage. Shoot length reached approximately 50 cm in NP-treated plants and only 40 cm in the NP-untreated control plants (Figure 1(c)), under well-watered conditions. NP treatment effectiveness was ranked as follows: Fe NP > Co NP > Cu NP > ZnO NP. Furthermore, although plant height increased with NP treatment, shoot dry weight decreased from 1.2 g for controls to 0.9 g for Fe NP-treated plants under well-watered conditions, with the rest of NP treatments lying intermediate between these two extremes (Figure 1(d)). Conversely, although under drought stress conditions the differences in shoot length among Co, Cu, and ZnO NP plants and controls were not statistically significant (Figure 1(c)), shoot dry weight of NP-treated plants was higher (; Tukey’s test) than that of untreated plants (Figure 1(d)).

Root growth was improved by NP treatment only under induced drought conditions. The roots of NP-treated plants were significantly longer (; Tukey’s test) than those of control plants; particularly, Fe NP-treated plants produced the longest roots (average 80.1 cm), followed by Cu, Co, and ZnO NP plants with lengths between 75.4 and 78.1 cm, compared to 67.6 cm in root length for control plants. Although root length increased, root dry weight was not significantly different between NP treatments and controls, neither in well-watered nor in induced drought conditions (Figures 2(a) and 2(b)). Such root response resulted in greater specific root length (SRL) of NP-treated plants under induced drought conditions (Figure 2(c)).

3.2. Assessment of NP Treatment Effects on the Drought Tolerance Index

To assess the physical responses of NP-treated soybean plants under induced drought conditions, we measured relative water content (RWC), drought tolerance index (DTI), and biomass reduction rate, which are common physiological traits used to estimate plant drought tolerance. While RWC values provide information on plant tolerance level under drought conditions, DTI reveals more accurately plant capacity for drought tolerance. Thus, the higher the RWC and DTI values, the more drought tolerant the cultivar is.

Almost all NP-treated plant groups displayed RWC values similar to that of the untreated groups under well-watered conditions (Figure 3(a)). Conversely, all plant groups showed a similar trend of considerable reduction of RWC after the experimental drought period. However, under such induced drought conditions, Fe and Cu NP-treated plants maintained RWC at 71%, which was significantly higher (; Tukey’s test) than the RWC of control plants (64%) (Figure 3(a)).

With respect to DTI, results shown in Figure 3(b) indicate that NP-treated plants exhibited higher DTI values compared to controls. Among NP-treated plants, Fe NP-treated plants showed the highest drought tolerance capacity, with a DTI value of 0.97, which was sevenfold higher than the value for untreated plants, followed by Co, Cu, and ZnO NPs (Figure 3(b)). Significantly, although Fe and Cu NP-treated plants had a similar RWC, the DTI value of Fe NP-treated plants was twofold higher and significantly different from the DTI value of Cu NP-treated plants (0.41).

Plant biomass of Cu and Zn NP-treated plants (0.63 g) was higher than that of Fe or Co NP-treated plants (0.53 and 0.59 g, respectively) under well-watered conditions. However, after the experimental drought period and subsequent rewatering, plant biomass of the former were 0.04–0.06 g lower than those of Fe or Co NP-treated plants or that of the control (Figure 4(a)). Therefore, the lowest biomass reduction rate was observed for Fe and Co NP-treated plants (36% and 40%, respectively), whereas it was 49% for Cu NP, 48% for Zn NP, and 46% for the control treatment (Figure 4(b)). The calculated biomass reduction rates indicated that Fe and Co NP-treated plants maintained higher biomass under drought than untreated or Cu and Zn NP-treated plants.

3.3. Expression Analysis of Drought-Responsive Genes in NP-Treated Soybean Plants

The expression of drought tolerance marker genes in leaves and roots of NP-treated and nontreated plants (controls) was analyzed using a qPCR method. The results shown in Figure 5 reveal that under induced drought conditions, selected drought tolerance marker genes (GmRD20A, GmDREB2, GmERD1, GmFDL19, GmNAC11, GmWRKY27, GmMYB118, and GmMYB174) were upregulated in roots or shoots of NP-treated plants, depending on the specific gene, the plant tissue, and the NP used. In Fe NP-treated plants, the expression of all genes analyzed was significantly higher in the roots ( or ; -test), while GmERD1, GmWRKY27, GmMYB118, and GmMYB174 showed significantly elevated expression in the leaves ( or ; -test). Additionally, under drought-induced conditions, Cu NP treatment raised the transcript level for GmDREB2, GmERD1, GmNAC11, and GmMYB174 in roots and for GmDREB2, GmERD1, GmNAC11, GmWRKY27, GmMYB118, and GmMYB174 in leaves. In turn, Co NP treatment increased the expression level of GmRD20A, GmERD1, GmNAC11, and GmWRKY27 in roots and of GmERD1, GmWRKY27, GmMYB118, and GmMYB174 in leaves. Lastly, fewer marker genes were upregulated in both roots and leaves of ZnO NP-treated plants; these were GmDREB2, GmERD1, and GmNAC11 in roots and GmERD1, GmMYB118, and GmMYB174 in leaves.

The increase in the expression level of GmWRKY27, GmMYB118, and GmMYB174 in the leaves of all four NP-treated plant groups was the most remarkable change in gene expression observed, namely, 8- to 10-fold, 6- to 20-fold, and 2- to 8-fold, respectively. Among all analyzed drought tolerance marker genes, GmERD1 was the only one that was upregulated in both roots and leaves, regardless of the NP used. Further, the expression of GmERD1 in leaves was highest (2.16-fold) in Fe NP-treated plants relative to untreated plants, followed by Co NP (2.03-fold), Cu NP (1.73-fold), and ZnO NP (1.4-fold). Thus, our results clearly demonstrated that the expression of drought tolerance marker genes was effectively induced by Fe, Cu, Co, and ZnO treatments.

4. Discussion

4.1. Metal-Based Nanoparticles Promote Plant Growth Responses to Water Stress Adaptation and Tolerance

The results described above clearly indicate that the NPs tested herein contributed to the increase in the tolerance of the soybean cultivar DT26 to drought, as reflected by our measurements of plant growth, RWC, DTI, and biomass reduction rate. NP-treated plants, especially with Fe NP, grew to a larger extent than control plants under both well-watered and induced drought conditions. Although the shoot length of NP-treated plants was similar to that of control plants, the dry weight of NP-treated plants markedly increased above the control level by the end of the drought period. This observation suggests continued photosynthetic activity in the leaves of NP-treated plants. It is also important to observe that root length and specific root length index (SRL, i.e., root length per unit dry mass) (Figure 2) of NP-treated plants increased compared to that of controls. They are among several root morphological characteristics associated with sustained productivity under drought conditions [41]. The higher SRL observed here suggests that NP-treated soybean plants might have supported root elongation under water stress, thereby permitting soybean plants to adapt/tolerate drought by accessing deeper soil water. Consistently, according to Comas et al. [42], small roots with larger SRL allow the plant to easily increase hydraulic conductivity by increasing the root surface area in contact with soil water and the soil volume that can be exploited for water.

The positive effects of metal-based NPs have been previously examined for drought tolerance in several other plant species. Thus, for example, Mozafari et al. [43] reported that the addition of iron NPs in an in vitro culture medium increased drought stress tolerance of Fragaria × ananassa Duch. Similarly, the drought stress endurance of Ocimum basilicum improved upon application of TiO2 NPs combined with gibberellin (GA3) [44]. Furthermore, wheat productivity improved under water deficit conditions with Cu and Fe NP application through the increase in superoxide dismutase (SOD) activity and sugar content [45].

In addition to improving morphological characteristics, NP treatment-induced increase in drought tolerance was also reflected by RWC, which is the most appropriate measurement for plant water status in case of a cellular water deficit. In this study, Fe and Cu NPs seemingly supported cell water retention under water deficit conditions. Consistently with our results, Cu and Zn NPs induced an increase in RWC by 8%–10% in leaves of seedlings of two different wheat varieties exposed to drought [13]. In addition, such physiological parameters have been frequently used to screen for drought-tolerant genotypes in wheat [46], soybean [47], and finger millet [48].

Although the effects of Fe and Cu NPs on water retention were equal, Fe NP was found to be most effective in increasing drought tolerance in soybean (7-fold increase in DTI compared to the control), followed by Co NP (with 4-fold increase in DTI; Figure 3(b)). DTI is considered a more accurate index for drought tolerance assessment because it combines all tolerance parameters, i.e., recovery factor, tolerance index, and recovery index. In practice, DTI has been calculated as a single factor or a combination with others to estimate drought adaptability in maize [49], rice [50], sugarcane [51], and soybean [35]. Therefore, Fe and Co NP treatments might be applied to field-grown soybean under water deficit conditions. In plants, Fe and Co are indispensable elements and especially important for legume species such as soybean, to sustain plant development, nodule production, and nitrogen fixation [6, 52]. The availability of these elements was necessary for enhancing or preparing plants to better tolerate drought stress [7, 5356]. The effect of Fe NP has also been investigated in cowpea [57], but this is the first report related to NP-induced drought tolerance in soybean.

Shoot dry weight of NP-treated and untreated plants declined considerably under drought conditions compared with that of well-watered controls; however, it was higher in NP-treated plants than in untreated plants after rewatering. Consequently, the biomass reduction rate was lower in NP-treated plants than in controls, and Fe and Co NP-treated plants showed the lowest biomass reduction rate. Overall, our observations indicated that NPs, especially Fe NPs, had positive effects on soybean plants that enabled them to better adapt to drought conditions.

4.2. Nanoparticles Induce Drought-Related Gene Expression

In this study, the combination of physiological responses and expression analysis of representative genes may help understand how soybean plants and NPs interact at a molecular level to mitigate drought stress. The expression of the selected drought tolerance marker genes, GmRD20A, GmDREB2, GmERD1, GmFDL19, GmNAC11, GmWRKY27, GmMYB118, and GmMYB174, was found to be upregulated in roots or shoots (or both) of NP-treated plants under drought. The above genes have been employed as water deficit tolerance markers in several studies on Arabidopsis [17] and other species, including soybean [18, 58]. Although these are drought-induced genes, according to Neves-Borges et al. [18] their expression pattern and induction levels were dependent on each specific gene and on drought experimental design. Among the tested NPs in our study, Fe NP was the most interesting because it was the only one to induce all marker genes in roots under drought stress conditions.

The increased expression of GmWRKY27, GmMYB118, and GmMYB174 in leaves proves that NP treatments do indeed trigger drought tolerance responses in plants. GmWRKY27 belongs to a superfamily of plant transcription factor genes, WRKY, and is involved in various physiological and developmental processes in plants, such as hormone signaling, catabolism, secondary biosynthesis, phosphate acquisition, lignin biosynthesis, seed germination, and stress responses [59]. The most important role of the WRKY transcription factor family is the regulation of plant defense and stress responses. WRKY proteins are known to regulate plant responses to a wide range of abiotic stress conditions including drought, flooding, heat, cold, heavy metal toxicity, low humidity, and osmotic, oxidative, salt, and UV stress [59]. According to Zhou et al. [60], GmWRKY27 was expressed at a low level in soybean roots and leaves under normal conditions but was highly induced by drought stress, salt stress, and abscisic acid (ABA) treatment. Notably, the expression of GmWRKY27 was found to increase strongly in the very early phase (0.5–1 h) of drought, salinity, or ABA treatment and decrease in the late phase (12 h) of treatment [60]. The high expression level of GmWRKY27 in the leaves of all NP-treated plants under drought conditions (Figure 5), but not in roots, suggests that NPs may be linked to the regulation of ABA biosynthesis and stomatal function under drought conditions. This is because plant cells are known to trigger ABA biosynthesis and accumulate ABA in the cytosol under stress conditions, leading to the regulation of stomatal opening and closure. This event is important to control water loss through transpiration under water deficit conditions [61]. On the contrary, WRKY27 might interact with other transcriptional activators to induce the expression of genes related to the ABA-dependent pathway, respond to water deficit stress, and trigger drought tolerance in soybean. The combined expression of GmWRKY27 and GmMYB174 is believed to suppress the expression of GmNAC29, a negative effector, resulting in plants becoming tolerant to abiotic stress [23]. Overexpression of GmWRKY27 also produced transgenic soybean plants tolerant to salinity and water deficit stress [23].

GmMYB118 is one of the 156 GmMYB genes identified in soybean by Liao et al. [62]. In this study, all NPs promoted increased expression of GmMYB118 in leaves, but only Fe NP treatment enhanced its expression in both roots and leaves under drought conditions. Du et al. [24] postulated that the GmMYB118 transcription factor might improve tolerance to drought and salt stress by promoting the expression of stress-associated genes. The expression of GmMYB118 (synonym: GmMYB173) is also known to induce the salinity tolerance process by regulating flavonoid biosynthesis to reduce reactive oxygen species (ROS) content [63]. Therefore, GmMYB118 might be considered one of the most important transcription factors regulating soybean adaptation to drought conditions.

The upregulation of GmERD1 in both leaves and roots makes it the most sensitive gene under induced drought conditions in NP-treated plants. ERD1 functions in the ABA-independent pathway; it is a functional gene involved in a cascade of reactions acting directly in response to abiotic stress [58, 64]. It was shown to be strongly expressed in a drought-sensitive soybean cultivar under water deficit conditions, but not in the tolerant counterpart [18, 58]. According to Neves-Borges et al. [18], ERD1 plays an important role at an earlier stage in the drought response pathway. Therefore, the increased expression of GmERD1 in leaves and roots of NP-treated plants suggests that ERD1 has a function in common responses at an early stage of an adverse signal response pathway in soybean. Neves-Borges et al. [18] considered that the upregulation of GmERD1 might be important for plants lacking traits needed for plant adaptation to drought, which might explain the drought stress responses of the non-drought-tolerant DT26 soybean cultivar, through the enhancement of GmERD1 expression by NP treatment.

In contrast with ERD1, RD20A is considered part of the ABA-dependent pathway [64]. The study of Neves-Borges et al. [18] on soybean indicated that the expression of RD20A dramatically increased in roots of drought-stressed plants growing in a potted soil system (Psy), but not in plants abruptly subjected to drought, as compared to plants growing in a hydroponic system (Hsy) that can be made to change suddenly. In this study, water stress treatment was carried out in potted soil in which soil moisture content was controlled by a method similar to those in the Psy system, which allowed plants to gradually adapt to water shortage. In this system, the expression of RD20A in our experiments was remarkably enhanced only in roots of Fe or Co NP-treated plants. The results indicated that Fe and Co NPs had a strong effect on the expression of drought marker genes in roots, which might explain the morphological changes of roots and DTI values in Fe and Co NP-treated plants.

In addition to RD20A, the expression of GmNAC11 increased in roots of all NP treatments, while that of DREB2 increased in roots of the Fe, Cu, and ZnO NP-treated plants. NAC transcription factors are also known to activate gene expression via stimulating other transcriptional activators such as DREB2 or ERD1 [64]. According to Hao et al. [25], GmNAC11 expression was induced by adverse environmental factors such as salinity and dehydration; transgenic plants were resistant to salinity following overexpression of GmNAC11. Hao et al. [25] also indicated that GmNAC11 had a role in transcription initiation in combination with DREB1A. The increased expression of GmNAC11 and GmDREB2 in the roots of most NP-treated plants strengthens the hypothesis that NPs in our study relate to the expression of genes involved in the ABA-independent pathway. Overall, metal-based NPs can affect the expression of drought-inducible genes through both ABA-dependent and ABA-independent pathways under drought conditions, subject to particular NP type and specific plant tissues.

It is well known that the response of plant cells to environmental stress is a complex process involving a network of signaling pathways [64]. Although precise modes of action of NPs in such conditions have been investigated, several hypotheses have emerged. The hypothesis that best matched our findings is that of Khan et al. [15], in which NPs are considered to act as Ca2+ or signaling molecules in the cytoplasm, identified by the calcium-binding protein (CaBP) complex or nanoparticle-specific proteins once inside plant cells. Consequently, a series of intracellular signaling pathways are induced and the expression of associated genes is activated, leading to enhanced plant tolerance responses to adverse conditions. In Arabidopsis, the expression of the RD20 gene was induced by drought and salinity conditions or by ABA treatment [65]. Further analysis indicated that the RD20 protein harbors a conservative and specific region for calcium ion binding (EF-hand) [65]. The enhancement of GmRD20A expression in this study, especially in Fe and Co NP-treated plants, clearly supports the proposed hypothesis that nanoparticles may participate in inducing Ca2+-binding protein expression. In addition, NPs are also believed to play a role in ROS scavenging by simulating antioxidative enzyme activities. Recently, Sun et al. [66] provided strong evidence for the above theory, in which the expression of Fe/Mn SOD, Cu/Zn SOD, APX, and CAT was significantly enhanced in ZnO NP-treated plants under drought. Based on our own findings and other reported studies, a possible mechanism of NP-mediated induction of drought tolerance of plants is illustrated in Figure 6.

As shown, NPs enter the cells by penetrating or by transport through specific channels in the cell membrane. They may act as stress signaling molecules and induce expression of drought-responsive genes (including regulatory factors) leading to defense system activation and tolerance to stress. NPs are proposed to mimic Ca2+ ions and bind with calcium-binding proteins (CaBP, CDPK, and CBL) [15]. The activated CaBP directly binds to the promoters of drought-responsive genes in the Ca2+ signaling pathway and then triggers the expression of downstream drought-related genes. Additionally, a suitable amount of NPs can maintain an appropriate level of ROS to stimulate the ROS signaling network for activation of the plant defense system under drought conditions. In this study, seven regulatory genes—GmWRKY27, GmMYB174, GmMYB118, GmERD1, GmDREB2, GmRD20A, and GmNAC11 (shown in orange boxes)—were highly expressed in NP-treated soybean plants under induced drought conditions. GmWRKY27 and GmMYB174 factors cooperatively inhibit the expression of GmNAC29, a transcriptional repressor, that leads to reduced ROS levels [23]. The high expression of GmWRKY27 in leaves may regulate ABA biosynthesis and the ABA signaling pathway. As previously reported, ERD1 and DREB2 belong to the ABA-independent pathways, while RD20A plays roles in the ABA-dependent pathway [64]. The NAC11 factor may activate gene expression by stimulating other transcriptional activators such as ERD1 or DREB2 in ABA-independent pathways [64]. In Arabidopsis, the RD20 protein contains a conserved EF-hand Ca2+-binding domain that has been shown to bind to Ca2+ [65], suggesting that GmRD20A may relate to the calcium signaling pathway and interact with Ca2+ ions or NPs. The GmMYB118 factor might improve drought tolerance by promoting the expression of drought-related genes, reducing the ROS level, and regulating osmotic substances and flavonoid biosynthesis [24, 63].

5. Conclusions

The results of the present study demonstrated that Fe, Cu, Co, and ZnO NP treatments effectively helped soybean plants at an early vegetative stage to adapt to drought stress. In addition to the improvement of shoot and root morphology and drought tolerance indices, the expressions of the tested drought tolerance marker genes were significantly enhanced under water deficit conditions in NP-treated plants, compared with untreated control plants. Therefore, metal-based NPs might promote plant tolerance to drought stress through the induction of drought-related gene expression. The usefulness of NP treatment for coping with drought stress in soybean and other plant species of commercial interest at other plant development stages warrants further research.

Data Availability

The data used to support the findings reported herein are available from the corresponding author upon request.

Conflicts of Interest

The authors declare no conflict of interest.

Acknowledgments

We gratefully thank A/Prof. Dr. N.P. Thao and her research group (School of Biotechnology, International University, Vietnam National University HCMC) for the technical support in drought experiments. We would also like to thank Prof. Rony Swennen, Dr. Jassmine Zorrilla, and Mrs. Hien Do (Laboratory of Tropical Crop Improvement, Department of Biosystems, Faculty of Bioscience Engineering, KU Leuven, Belgium) for their support in the qPCR experiments and analysis. Finally, we thank Editage (https://www.editage.com) for English language editing. This research was funded by the key project “Application of Nanotechnology in Agriculture” (code: VAST.TĐ.NANO-NN/15-18) financed by the Vietnam Academy of Science and Technology.