Shanshan Zhu1,2
Junzhen Mi1,2,3
Baoping Zhao1,2,3Zhaoming Wang4Zhixue Yang1,2Mengxin Wang1,2Jinghui Liu1,2,3*- 1Coarse Cereals Industry Collaborative Innovation Center, Inner Mongolia Agricultural University, Hohhot, China
- 2National agricultural scientific research outstanding talents and their innovation team, Inner Mongolia grassland talents innovation team, Hohhot, China
- 3Oat Engineering Research Center of Inner Mongolia Agricultural University, Oat Engineering Laboratory of Inner Mongolia Autonomous Region, Hohhot, China
- 4National Center of Pratacultural Technology Innovation (under way)/M-Grass Ecology And Environment (Group) Co., Ltd., Hohhot, China
Drought stress inhibits oat growth and yield. The application of fulvic acid (FA) can improve the drought resistance of oats, but the corresponding molecular mechanism of FA-mediated drought resistance remains unclear. Here, we studied the effects of FA on the drought tolerance of oat leaves through physiological, transcriptomic, and metabolomics analyses, and identified FA-induced genes and metabolites related to drought tolerance. Physiological analysis showed that under drought stress, FA increased the relative water and chlorophyll contents of oat leaves, enhanced the activity of antioxidant enzymes (SOD, POD, PAL, CAT and 4CL), inhibited the accumulation of malondialdehyde (MDA), hydrogen peroxide (H2O2) and dehydroascorbic acid (DHA), reduced the degree of oxidative damage in oat leaves, improved the drought resistance of oats, and promoted the growth of oat plants. Transcriptome and metabolite analyses revealed 652 differentially expressed genes (DEGs) and 571 differentially expressed metabolites (DEMs) in FA-treated oat leaves under drought stress. These DEGs and DEMs are involved in a variety of biological processes, such as phenylspropanoid biosynthesis and glutathione metabolism pathways. Additionally, FA may be involved in regulating the role of DEGs and DEMs in phenylpropanoid biosynthesis and glutathione metabolism under drought stress. In conclusion, our results suggest that FA promotes oat growth under drought stress by attenuating membrane lipid peroxidation and regulating the antioxidant system, phenylpropanoid biosynthesis, and glutathione metabolism pathways in oat leaves. This study provides new insights into the complex mechanisms by which FA improves drought tolerance in crops.
1 Introduction
Oat (Avena nuda L.) is a globally significant food and feed crop. Its grains are abundant in nutrients, including high levels of β-glucan, protein, fat, and soluble dietary fiber. It exhibits remarkable effects on blood sugar regulation and cholesterol reduction, making it an important medicinal crop (Alemayehu et al., 2023). With the increasing focus on nutrition and dietary health, there has been a gradual rise in market demand for high-quality oat raw materials (Singh et al., 2022). Therefore, the development of the oat industry contributes to agricultural growth and farmers’ income and holds great significance in addressing issues related to residents’ imbalanced nutritional intake. China is the primary origin of large-grain naked oats, with a relatively concentrated planting area primarily located in arid and semi-arid regions such as North and Northwest China. These areas are characterized by harsh natural conditions, and as a result of uneven precipitation distribution and the escalating occurrence, frequency, and intensity of droughts attributed to global warming, drought has emerged as the predominant limiting factor for oat production (Tian et al., 2022; Zhang et al., 2022; Salih et al., 2023). Drought and water scarcity triggers a cascade of physiological and biochemical responses in plant morphology (Wang et al., 2022a). For example, drought stress induces the excessive accumulation of reactive oxygen species (ROS), resulting in oxidative damage to cell membranes and chlorophyll degradation, reduced photosynthetic efficiency, and inhibition of crop growth (Muhammad et al., 2021; Peng et al., 2022; Hu et al., 2023), all of which exert a negative impact on plants. Therefore, exploring appropriate methods to enhance the drought tolerance of oat plants is crucial for optimizing their growth and yield (Gao et al., 2018; Bai et al., 2021). The utilization of plant growth regulators is regarded as a proficient and eco-friendly approach to address this issue (Che et al., 2017; Jesmin et al., 2023).
Fulvic acid (FA) is a low-molecular-weight, highly polar active humic substance that contains active organic functional groups and hydrophilic free radicals (Gong et al., 2020). As a growth regulator, it has multiple regulatory effects on plant growth and development (Liang et al., 2024), can stimulate plant growth, enhance leaf photosynthesis, and regulate antioxidant systems (Justi et al., 2019; Chen et al., 2022b). Most importantly, FA is a key factor affecting plant tolerance to abiotic stresses. For example, its strong antioxidant activity helps maintain the homeostasis of ROS to protect plants from the damage caused by environmental stress (Wang et al., 2019). FA can also activate various antioxidant enzymes, such as superoxide dismutase (SOD), catalase (CAT), peroxidase (POD), and polyphenol oxidase (PPO), to reduce oxidative stress in plants (Yu et al., 2024; Hareem et al., 2024). Plant photosynthesis is highly sensitive to drought stress. FA can improve net photosynthesis, increase chlorophyll content and electron transport, and maintain chloroplast ultrastructure, thus alleviating the adverse effects of abiotic stress on photosynthesis (Liu et al., 2022). In addition, studies on its molecular mechanism reveal that FA can resist drought stress through the gene expression involved in the primary metabolism of maize leaves, particularly photosynthesis, carbon fixation, hormone, and osmotic metabolism pathways (Chen et al., 2022a). FA can also enhance the antioxidant defense ability of tea plants under drought stress by regulating ascorbic acid metabolism, glutathione metabolism, and flavonoid biosynthesis (Sun et al., 2020). At present, most research focuses on the drought resistance of FA-induced crops at the physiological level, while the molecular mechanism of fulvic acid-mediated drought resistance in oats remains unclear. Transcriptome analysis is commonly employed to investigate alterations in gene expression levels under diverse biotic and abiotic stresses, thereby enhancing our understanding of the functionality of the gene regulatory networks and signaling pathways in plants (Zhao et al., 2021; Jiang et al., 2023). Similarly, metabolomics has emerged as a novel tool facilitating comprehensive component analysis, proving valuable for identifying changes in metabolites induced by various environmental fluctuations (Zhang et al., 2021). Integrative transcriptomic and metabolomic analyses establish a robust connection between gene regulation and metabolite generation, facilitating the comprehension of plant responses to dynamic environmental changes (Elakhdar et al., 2023; Cui et al., 2023). Scholars have successfully applied these methods to study plant physiology (Li et al., 2023a; Wang et al., 2023a; Dwivedi et al., 2023). Sun et al. (2020) revealed a crosstalk regulation between FA and ascorbate metabolism and flavonoids biosynthesis, contributing to understanding plant drought tolerance. In our previous experiments, we demonstrated the ability of foliar spraying to mitigate the detrimental effects of drought stress on oat seedlings (Zhu et al., 2022). However, the combined transcriptomic and metabolomic analyses of drought tolerance in oats have rarely been reported. In this study, we employed integrative metabolomic and transcriptomic approaches to investigate the effects of FA on overall leaf responses (including morphology, physiology, gene transcription, and metabolite generation) in oat leaves under drought conditions. The objective of the study is to investigate the physiological and molecular mechanisms by which FA regulates oats under drought stress. Our findings enhance the understanding of FA regulation mechanisms in oat plants and provide valuable strategies for improving crop drought tolerance.
2 Materials and methods
2.1 Plant material and experimental treatments
The drought-sensitive oat variety ‘Bayou 9’ was used in this investigation, characterized by a fertility period of approximately 80 to 90 days. Thirty seeds were sown into plastic pots of 25 centimeters in diameter and 18 centimeters in depth, each containing roughly 3 kilograms of mixed soil (soil, vermiculite, peat soil in a volume ratio of 1:1:1). By the third leaf stage, each pot was thinned to 20 seedlings. The seedlings were grown in a controlled greenhouse at Inner Mongolia Agricultural University (Hohhot, China) at a temperature of 20 ± 2°C, under a photoperiod of 16 hours of light and 8 hours of darkness, with a light intensity of 2000 lux. All pots were at 75% water prior to the jointing stage; water content was maintained by weighing every two days and applying water as needed. Drought stress (water content maintained at 45%) was applied to the seedlings at the jointing stage, and the stress persisted until maturity.
Three treatments were implemented, each with three replicates (pots): (1) spraying distilled water and watering to 75% field capacity (CK); (2) spraying distilled water and exposed to drought stress (DS); and (3) spraying FA and exposed to drought stress (DF). Sprinkle FA once during the plant’s early jointing, heading, and filling stages, at a rate of 15 ml per pot. The FA concentration (600 mg/L) used in this investigation was based on our previous research (Zhu et al., 2022). After 7 days of FA treatment (7 days based on the previous research of the research group) (Li et al., 2019), flag leaf samples were collected, quickly frozen in liquid nitrogen, and stored at -80°C for future measurements.
2.2 Growth analysis
The flag leaf, penultimate leaf, and third to last leaf of the oat plant were selected as the study objects. The length and width of the leaves were measured with a ruler. The leaf area was calculated based on the length–width coefficient method:
where 0.73 is an empirical coefficient. The chlorophyll content of the flag leaves was determined using a chlorophyll meter (SPAD-502, Konica Minolta, Tokyo, Japan). Dry matter quality was estimated by taking the flag leaf, penultimate leaf, and third of 20 oat plants were put into paper bags, and then put into oven at 105°C for 20 min and subsequently dried at 80°C for 48 hours to constant weight. The relative water content (RWC) of leaves was determined by the drying and weighing method (Farooq et al., 2009). Briefly, the fresh weight (W1) of oat flag leaves was weighed and the leaves were then placed on a water surface for 12 h to obtain the saturation weight (WS). Finally, the leaves were dried at 85°C to a constant weight (W2). The RWC of the oat leaves was calculated as follows:
2.3 Measurements of physiological characterization
Malondialdehyde (MDA) content was determined by thiobarbituric acid (TBA) colorimetric method. superoxide dismutase (SOD), peroxidase (POD), catalase (CAT), phenylalanine aminolyase (PAL) activities, hydrogen peroxide (H2O2), 4-coumaroyl coenzyme A ligase (4CL), and dehydroascorbic acid (DHA) were all measured by enzyme-linked immunosorbent assay kit instruction (Keshun Biotechnology Co., Ltd., Shanghai, China).
2.4 RNA extraction, cDNA library construction, and transcriptome sequencing
Total RNA was isolated using the mirVanaTM miRNA isolation kit, following to the manufacturer ‘s instructions. The Tru Seq Stranded mRNA LT Sample Prep Kit was used to create the sequencing library, following to the manufacturer ‘s instructions.
The Agilent 2100 Bioanalyzer (Agilent, Santa Clara, CA, USA) was used to determine the library size and purity of the sample (1 μL). Transcriptome sequencing was carried out by OE Biotechnology Co., Ltd. (Shanghai, China). Trimmomatic software was used to handle raw data, perform quality control, remove the adaptor, and filter to obtain clean reads (Bolger et al., 2014). To get the reference genes, clean reads were aligned to the reference genome using hisat2 (Kim et al., 2015). All sequencing readings were uploaded to the National Center for Biotechnology Information (SRA accession number: PRJNA1111847).
2.5 Gene annotation, differential expression, and enrichment analyses
The gene function was annotated based on the following databases: Gene Ontology (GO) (http://www.geneontology.org/) and Kyoto Encyclopedia of Genes and Genomes (KEGG) (http://www.genome.jp/kegg/). Cufflinks software was utilized to calculate the fragments per kilobase of transcript per million mapped reads (FPKM) value for each gene, while htseq-count was employed to obtain the read counts of each gene. Differentially expressed genes (DEGs) were identified using the DESeq R package by estimating size factors and performing a negative binomial test. Significantly differential expressions were determined by applying the thresholds p-value ≤ 0.05 and |log2FC| ≥ 1. Hierarchical cluster analysis was conducted on DEGs to explore patterns in gene expression. GO and KEGG pathway enrichment analyses were performed on the DEGs using R software version 3.6.2 with employing the hypergeometric distribution.
2.6 Metabolite extraction
The samples (80 mg) were meticulously transferred to a 1.5 mL Eppendorf tube, followed by the addition of two small steel balls. An internal standard solution of L-2-chlorophenylalanine (0.3 mg/mL) dissolved in methanol was added at a volume of 20 μL, and each sample received a mixture of methanol and water (7/3, v/v) at a volume of 1 mL. Subsequently, the samples were placed at −20°C for 2 min. The entire sample was then subjected to ultrasonication in an ice-water bath for 30 min and subsequently cooled in a refrigerator at −20°C for 20 min. After centrifugation at 4°C (13,000 rpm) for 10 min, supernatant aliquots of 150 μL were collected from each tube using a crystal syringe and filtered through a microfilter with a pore size of 0.22 μm before being transferred to LC vials. These vials were stored at −80°C until liquid chromatography-mass spectrometry (LC-MS) analysis was performed. A quality control (QC) sample was prepared by pooling all individual samples into one mixed sample.
LC-MS analysis was conducted using a UHPLC system (Vanquish, Thermo Fisher Scientific, Waltham, MA, USA) coupled to a Q Exactive HFX mass spectrometer (Orbitrap MS, Thermo Fisher Scientific) with a UPLC BEH Amide column (2.1 mm × 100 mm, 1.7 μm). The mobile phase consisted of 25 mmol/L ammonium acetate and 25 mmol/L ammonia hydroxide in water (pH = 9.75) (A) and acetonitrile (B). Gradient elution was employed for the analysis as follows: 0–0.5 min, 95% B; 0.5–7.0 min, 95%–65% B; 7.0–8.0 min, 65%–40% B; 8.0–9.0 min, 40% B; 9.0–9.1 min, 40%–95% B; 9.1–12.0 min, 95% B. The column temperature was maintained at a constant value of 30°C throughout the analysis period, while the auto-sampler temperature was set at −4°C and the injection volume used was fixed at exactly 2 μL. MS spectra were obtained in information-dependent acquisition mode using Xcalibur (Thermo Fisher Scientific) and a QE HFX mass spectrometer equipped with an electrospray ionization source. Measurements were made under the following conditions: sheath gas flow rate of 50 Arb, auxiliary gas flow rate of 10 Arb, capillary temperature of 320°C, full MS resolution set at 6000 resolving power, MS/MS resolution set at 7,500 resolving power with collision energy values in the range 10–60 under normalized collision energy mode depending on specific requirements during the data acquisition process, and positive or negative spray voltage (3.5 kV or −3.2 kV, respectively).
2.7 Data preprocessing and statistical analysis
The raw data were preprocessed and normalized using the metabolomics processing software Progenesis QI v3.0 (Nonlinear Dynamics, Newcastle, UK). The Human Metabolome Database (HMDB), Lipidmaps (v2.3), METLIN, and PMDB database were utilized for the qualitative analysis. The selected data were scored according to the qualitative results of the compound with a standard score of 36 points, and data with a score below 36 points were deleted. The resulting positive and negative ion data were combined into a data matrix table. First, unsupervised principal component analysis (PCA) was used for multivariate statistical analysis to observe the overall distribution among the oat samples and the stability of the whole analysis process. Supervised partial least squares-discriminant analysis (PLS-DA) and orthogonal analysis. Partial least squares analysis (OPLS-DA) was used to distinguish overall differences in metabolic profiles between groups and to find inter-group differences. Differential metabolites. Student’s t-test and fold change analysis were employed. To compare metabolites between the two groups. The method of combining multidimensional analysis and one-dimensional analysis was adopted. The following thresholds were used for the analysis: VIP value of the first principal component of the OPLS-DA model > 1; FC > 1; and t-test p-value < 0.05. Quasi-screening was adopted for differential metabolites between groups. Functional annotation and pathway analyses of differential metabolites were performed using KEGG database enrichment analysis.
2.8 Real-time quantitative reverse transcription PCR analysis
To verify the reliability of the DEGs from the sequencing results, eight DEGs were selected for internal q-PCR. The specific primers are reported in Supplementary Table 1, with the oat actin gene employed as the internal reference gene. The RNA was reverse transcribed into cDNA using TransScript All-in-One First-Strand cDNA Synthesis SuperMIX for q-PCR (All-in-One Gold). q-PCR was performed on an ABI9700 system using the dye method. The reaction system, described in Supplementary Table 2, was executed under the following conditions: 94°C for 30 s; 40 cycles of 94°C for 5 s; and 60°C for 30 s. The temperature was slowly increased from 60°C to 97°C, and the fluorescence signal was collected five times per °C. Each sample was replicated three times, with sterile water used as the control. The relative expression of the genes was calculated with the 2-ΔΔCt algorithm.
3 Results
3.1 Morphological and physiological responses of the oat plant induced by FA under drought stress
We investigated the morphological and physiological changes in FA-treated leaves under drought conditions to confirm the positive effects of FA in alleviating drought stress in oat leaves. Under drought stress, the oat plants exhibited a significant reduction in height, accompanied by a decrease in leaf count. The lower leaves displayed green and yellow discoloration and curling, while the upper leaves showed scorching at the edges. Additionally, there was a notable decline in spikelet numbers. The number of dead leaves, the length of dead leaves, and the degree of yellowing all exhibited a significant reduction in FA-treated leaves under drought stress in comparison to the control group. There was also an observed increase in the number of spikelets (Figure 1A). The leaf area and dry matter quality of the oat leaves were significantly reduced under drought stress, while the application of FA mitigated the detrimental effects of drought stress on these parameters (Figures 1B, C). Under drought stress, the chlorophyll content and RWC declined by 23.29% and 47.48%, respectively, compared to the control, while the application of FA resulted in a significant recovery of chlorophyll content and RWC by 24.5% and 23.33%, respectively (Figures 1D, E). These results indicate that FA alleviates the negative effects of drought on oat leaves and promotes plant growth.
Figure 1. Effects of fulvic acid on the morphology and physiology of oat (Avena sativa L.) under drought stress. CK was the control, DS was exposed to drought stress, and DF was exposed to drought stress and was treated with fulvic acid. (A) growth condition. (B) chlorophyll content. (C) leaf area. (D) dry weight. (E) relative water content. The columns represent the means and the error bars represent ± SD of three replicates. Columns in the same panel and capped with the same lower case letter are not significantly different P > 0.05 according to Duncan’s multiple range tests.
MDA and H2O2 levels were selected as markers of oxidative damage and ROS accumulation. Drought stress resulted in a significant enhancement in the levels of H2O2 and MDA in the oat leaves, while the application of FA effectively mitigated the accumulation of H2O2 and MDA (Figures 2A, B). We estimated the levels of antioxidative enzymes in the oat leaves under different treatments. The changes in these enzymatic antioxidants exhibited a similar variation trend (Figures 2C-G). Under water-deficit conditions, the activities of the antioxidative enzymes (SOD, POD, PAL, CAT, and 4CL) significantly increased in oat leaves, and FA application was able to further improve the levels of SOD, POD, PAL, CAT, and 4CL in the oat leaves. Ascorbic acid is one of the most abundant water-soluble antioxidants and can be used as a ROS scavenger directly against oxidative stress. Drought stress led to an increase in dehydroascorbic acid content in oat leaves, while fulvic acid application decreased its content (Figures 2H). These results suggest that the application of FA can enhance antioxidant defense activity and decrease excessive ROS.
Figure 2. Effects of Fulvic acid on (A) MDA content, (B) H2O2 content, (C) SOD activity, (D) POD activity, (E) CAT activity, (F) PAL activity in oat,(G) 4CL activity, and (H) DHA content (Avena sativa L.) leaf under drought stress. CK was the control, DS was exposed to drought stress, and DF was exposed to drought stress and was treated with fulvic acid. The bars represent the means ± SD of three replicates. Different letters indicate significant differences at P < 0.05 according to Duncan’s multiple range tests.
3.2 Leaf transcriptome profiles in response to different treatments
3.2.1 Transcriptome sequencing analysis
RNA-seq analysis was performed on three biological replicates of leaf samples under CK, DS, and DF treatments. Nine cDNA libraries were prepared and then paired-end sequenced. The raw reads ranged from 44.26 to 49.06 million. After removing low-quality reads, clean reads ranged between 44.16 and 48.96 million. In our libraries, the percentage of clean reads was greater than 99%, and the Q30 values were greater than 90.77%. The mean GC content was determined as 55.14% (Supplementary Table 3). The transcriptome sequencing results exhibit high quality and are suitable for the subsequent analysis.
Pearson’s correlation coefficient analysis of the three replicates indicated that the sequencing data were of high quality and met the requirements for further analyses (Figure 3A). Based on the thresholds of adjusted |log2Fold Change| > 1 and false discovery rate (FDR) of <0.05, a total of 2126, 2786, and 652 DEGs were screened in DS-vs-CK (1138 upregulated and 988 downregulated), DF-vs-CK (1532 upregulated and 1254 downregulated), and DF-vs-DS groups (475 upregulated and 177 downregulated)groups, respectively (Figure 3B). There were more upregulated genes than downregulated genes in the seedlings following the DS and DF treatments, and more genes were differentially regulated by drought alone than by the DF treatment. Among the DEGs, 150 transcripts were commonly induced by the CK, DS, and DF treatments. Moreover, we identified 225 common DEGs between the DF-vs-DS and DS-vs-CK groups. These genes are likely to be involved in the alleviation of drought stress in oats through the action of FA (Figure 3C). The gene expression trend verified by q-PCR was similar to that determined by RNA-seq (Supplementary Figure 1), indicating that the data obtained by transcriptome analysis were highly reliable.
Figure 3. The correlation among samples and the contrast between groups of differentially expressed genes. (A) Heat map of correlation coefficient between samples. (B) The total number of differentially expressed genes (DEGs) and upregulated and downregulated DEGs under different treatments. (C) Venn diagram of DEGs. CK, Normal moisture treatment; DS, drought treatment; DF, Fulvic acid + drought treatment.
3.2.2 GO and KEGG enrichments of DEGs
To further understand the functions and the related biological processes of the DEGs, GO enrichment analyses were conducted. The results are presented in Figures 4A, B. The DEGs were classified into biological processes, cellular components, and molecular functions. In the biological process of DS-vs-CK, the differential genes were mainly enriched in the defense response and flavonoid biosynthetic process. In the cellular component, the most enriched DEGs were the integral component of the membrane and the chloroplast thylakoid membrane. The molecular function mainly included DNA-binding transcription factor activity and peroxidase activity. (Figure 4A). These terms have been reported vital for drought stress tolerance in plants. The most differentially enriched genes of DF-vs-DS in biological processes were the defense response, cell wall organization, and the abscisic acid-activated signaling pathway. The DEGs enriched in cellular components were the nucleus and extracellular region. DNA-binding transcription factor activity and sequence-specific DNA binding were the main types of genes with molecular function enrichment (Figure 4B). Thus, genes encoding functionally diverse proteins contribute to the FA-mediated response of oat plants to drought stress, and DF treatment uniquely and differentially regulates genes for multiple functions.
Figure 4. Top 20 enriched GO terms and KEGG pathways from DEGs. (A) GO terms of DEGs in DS-vs-CK. (B) GO terms of DEGs in DF-vs-DS. (C) KEGG pathway analysis of DEGs in DS-vs-CK. (D) KEGG pathway analysis of DEGs in DF-vs-DS.
To identify the main pathways by which FA contributes to drought tolerance, the DEGs were subjected to KEGG pathway enrichment analysis. Figures 4C, D depicts the scatter plot of enrichment analysis combined with the number of DEGs-identified significant pathways. In total, 785 (DS-vs-CK) and 150 (DF-vs-DS) DEGs were mapped to 108 and 67 KEGG pathways, respectively (Supplementary Table 4). The photosynthesis-antenna proteins, phenylpropanoid biosynthesis, photosynthesis, ascorbate and aldarate metabolism, and glutathione metabolism pathways were enriched in the DS-vs-CK group (Figure 4C). There were more upregulated genes in the pathway of photosynthesis-antenna protein and photosynthesis than downregulated genes, while the genes in the pathway of phenylpropanoid biosynthesis, ascorbate and aldarate metabolism, and glutathione metabolism showed an opposite trend. This indicates that the photosynthesis of oat leaves in the drought treatment group was weakened, and some secondary metabolites and their biosynthetic pathways were activated to enhance the drought tolerance of oat (Supplementary Table 5). Phenylpropanoid biosynthesis and glutathione metabolism pathways were significantly enriched in the DF-vs-DS group (Figure 4D), and the number of upregulated genes was greater than that of downregulated genes (Supplementary Table 5). This suggests that FA improved the antioxidant defense system and drought resistance of plants, mainly by promoting the expression of phenylpropanoid biosynthesis and glutathione metabolism.
Photosynthesis significantly affected the expression of genes encoding key enzymes in the photosynthesis pathway, including photosystem II, photosystem I, photosynthetic electron transport, and F-type ATPase (Supplementary Figure 2). Compared with CK, the genes related to photosynthesis were downregulated under drought stress, with the exception of PetH(ferredoxin–NADP+reductases). Similarly, 25 expressed genes related to the light-harvesting chlorophyll (LHC) protein complex in photosynthetic antenna proteins were downregulated. However, these same genes—encoding PsbS and Lhcb1—were upregulated by FA under drought stress.
3.3 Leaf metabolome profiles in response to different treatments
The dataset employed in this study is complex as the difference between groups is relatively small while the difference within groups is relatively large. Thus, to obtain more reliable differential metabolites between the control groups, OPLS-DA was used to analyze the data collected from the oat leaf samples. In the OPLS-DA score diagram, the separation of two different treatment samples was obvious, and the intra-group clustering was also evident, indicating that there were different metabolites between groups (Figure 5A). The heatmap visualized distinct hierarchical clustering of metabolites obtained in all leaf samples, suggesting that these metabolites have different expression patterns in response to the CK, DS, and DF treatments (Figure 5B). These findings confirmed that the reproducibility among the biological samples was sufficient for further analyses.
Figure 5. Overview of widely targeted metabolome analysis of leaves responsive to CK, DS, and DF treatments. (A) OPLS-DA score map. (B) Heatmap visualization of the metabolites. The transition from blue to red indicates a gradient of expression abundance for metabolites, ranging from low to high; the redder the color, the higher the expression abundance of differential metabolites.
The significant differences in the relative metabolite content were screened using VIP >1 and p-value <0.05 (t-test results). In total, 8054 DEMs were detected in all leaf samples and were mainly classified into benzene and substituted derivatives (289), carboxylic acids and derivatives (633), fatty acyls (181), flavonoids (224), glycerophospholipids (322), organooxygen compounds (660), polyketides (261), prenol lipids (806), steroids and steroid derivatives (264), and others (3414) (Figure 6A; Supplementary Table 6). There were 724, 571, and 700 DEGs screened in DS-vs-CK (433 upregulated and 291 downregulated), DF-vs-DS (329 upregulated and 242 downregulated), and DF-vs-CK groups (448 upregulated and 252 downregulated), respectively (Figure 6B). Therefore, the differences in metabolite expression induced by FA in oats changed significantly under drought stress. Overlapping analysis showed that 132 metabolites commonly responded to DS-vs-CK, DF-vs-DS, and DF-vs-CK (Figure 6C). Interestingly, in the 132 DEMs, fatty acyls, flavonoids, isoflavonoids, and glycerophospholipids were highly accumulated, indicating that FA could improve the drought resistance of oats by regulating the accumulation of lipids and inducing the change of antioxidants (Figure 6D; Supplementary Table 7).
Figure 6. (A) Category of the DAMs. (B) Number of differential metabolites between different control groups (C) Venn diagram illustrating the number of DAMs in DS-vs-CK, DF-vs-DS, and DF-vs-CK groups. (D) Heat map analysis of common DEMs in DS-vs-CK, DF-vs-DS, and DF-vs-CK groups.
3.4 Integrated analysis of the RNA-seq and metabolome data of oat leaves
The DEGs and DAMs from the same group were mapped onto the KEGG pathway database to elucidate the associations between genes and metabolites. A total of 50 and 29 pathways were enriched in DS-vs-CK and DF-vs-DS, respectively, in which genes related to phenylpropanoid biosynthesis and glutathione metabolism pathways were significantly enriched in both control groups, and two metabolic pathways were also enriched in the metabolome (Figure 7; Supplementary Table 8). Moreover, these two pathways are important for the antioxidant system. Therefore, we focused on the phenylpropanoid biosynthesis and glutathione metabolism pathways.
Figure 7. Differential enrichment pathways for metabolites and genes. (A), DS-vs-CK, (B), DF-vs-DS, The horizontal axis represents the number of differentially enriched metabolites and genes in the pathway, and the vertical axis represents the KEGG pathway name. The meta represents the metabolome, and the gene represents the transcriptome.
In the phenylpropanoid biosynthesis pathway, two genes encoding phenylalanine ammonia-lyase (PAL) and 4-coumarate-CoA ligase (4CL) were upregulated, four genes encoding β-glucosidase were upregulated and one gene was downregulated, nine genes encoding peroxidase (POD) were upregulated and seven genes were downregulated, and the metabolites p-coumaric acid and ferulic acid were upregulated under drought stress. These results indicate that the phenylpropanoid biosynthesis pathway experienced significant changes under drought stress, and the aforementioned genes and metabolites played an important role in protecting oats from drought stress. FA also significantly affected the phenylpropanoid biosynthesis pathway. Most genes were still upregulated under FA, such as PAL, 4CL, β-glucosidase, and POD. The expression patterns of metabolites p-coumaric acid and ferulic acid were opposite to those under the drought treatment, both of which were downregulated while sinapitol was upregulated. These results indicated that FA could improve the drought resistance of oats by promoting the expression of genes related to the phenylpropanoid biosynthesis pathway and regulating metabolites (Figure 8A). In the glutathione metabolism pathway, thirteen genes encoding GST were significantly upregulated, the metabolite dehydroascorbic acid was upregulated, and L-glutamate was downregulated under drought stress. FA application increased the expression of six genes encoding GST and the level of L-glutamate and decreased G6PD and dehydroascorbic acid. These results indicate that FA treatment induced relatively higher expression of glutathione metabolism-related genes and L-glutamate biosynthesis than drought (Figure 8B).
Figure 8. The expression of important pathway genes and metabolites in oat leaves under CK, DS and FA treatments. (A) Expression profiles of the genes and metabolites involved in phenylpropanoid biosynthesis. (B) Expression profiles of the genes and metabolites involved in glutathione metabolism. The rectangular patterns represent the genes or metabolites, and the heatmap at the corresponding place depicts the differential expression of each identified gene or metabolite, which ranges from green or blue (low) to red (high). (C) The related networks of DEGs and DAMs in Phenylpropanoid biosynthesis. (D) The related networks of DEGs and DAMs participating in glutathione metabolism and dam. In the related network diagram, blue squares represent metabolites, yellow squares represent genes, red lines indicate positive correlation, blue lines indicate negative correlation, and the thickness of the lines represents the strength of the correlation.
Network correlation analysis was conducted on the DEGs and DAMs in these two pathways. In phenylpropanoid biosynthesis, p-Coumaric acid (C00811), and ferulic acid (C01494) exhibited a negative correlation with the genes encoding 4CL, PAL, and β-glucosidase, a positive correlation with four genes encoding POD, and a negative correlation with three genes. Sinapyl alcohol (C02325) showed an opposite expression pattern with p-Coumaric acid (C00811) and ferulic acid (C01494) (Figure 8C). In glutathione metabolism, L-glutamate (C00025) exhibited a positive correlation with genes encoding GST and a negative correlation with genes encoding G6PD. The expression pattern of dehydroascorbic acid (C05422) was opposite to that of L-Glutamate (C00025) (Figure 8D). These findings further confirmed that fulvic acid induced DEGs and DAMs were mainly concentrated in phenylpropanoid biosynthesis and glutathione metabolism pathways, and the differential genes were highly correlated with corresponding metabolites, indicating their importance in FA mediated drought stress response in plant.
4 Discussion
Under adverse conditions, FA can promote crop growth and increase yield (Chen et al., 2022b; Jesmin et al., 2023). When plants are subjected to drought stress, the water content in the body decreases and growth is delayed, resulting in a decrease in yield. In contrast, spraying FA can significantly improve plant growth (Shokhmgar et al., 2023). In the present study, under drought stress, spraying with FA was observed to reduce the number of dead oat leaves, deepen the green color of the leaves, and increase the number of spikelets (Figures 1A-C). The low molecular weight of FA enables it to penetrate through the pores of the membrane and promote nutrient entry into the cells by forming complexes with cations to increase the leaf area of oat leaves and promote the accumulation of dry matter (Bocanegra et al., 2006; Lalas et al., 2018). As the most important and effective pigment for photosynthesis, chlorophyll can reflect the plant’s ability to assimilate substances to a certain extent, and its content is an important factor in determining the plant’s photosynthetic capacity and substance production (Walker et al., 2018). Leaf RWC reflects the degree of soil aridity as well as the plant’s ability to retain water, and is most relevant to crop drought tolerance (Utsumi et al., 2019). (Figures 1D, E). Consistent with this, the expression levels of genes involved in chlorophyll synthesis were changed by FA (Supplementary Figure 3). FA decreased the expression levels of genes involved in chlorophyll degradation (SGR, NOL, and PAO) and increased the expression levels of genes involved in chlorophyll biogenesis (CHLH and CAO). This indicates that the occurrence of severe cellular water deficit caused by drought-affected chlorophyll synthesis, Which accelerated the decomposition rate of the original chlorophyll and consequently decreased in chlorophyll content (Javadi et al., 2017; Li et al., 2024). In contrast, FA increases chlorophyll content by improving the water retention capacity of oat leaves, allowing the leaves to absorb more light energy, mitigating chlorophyll degradation, and promoting chlorophyll synthesis (Nikoogoftar-Sedghi et al., 2024; Li et al., 2024). This result further supports the expression of antenna proteins, particularly LHCB1 (Supplementary Figure 2), whose genes are involved in the regulation of light capture in photosystem II (Bielczynski et al., 2020). FA may alleviate stress injury by regulating the light-capture protein of oat plants under drought stress. In wheat, the overexpression of LHC in chloroplast tissue enhances stress tolerance (Chen et al., 2023). The up-regulation of the photosynthesis-related protein PsbR further confirmed the role of FA in increasing photosynthetic efficiency of oat plants under drought stress (Wang et al., 2019) (Supplementary Figure 2). In conclusion, FA can enhance the drought tolerance and growth of oats by increasing chlorophyll, antenna protein, and photosynthetic pathway-related gene expression.
Elevated levels of H2O2 in plant tissues under drought stress may lead to oxidative stress, and the excessive accumulation of ROS leads to an increase in MDA, a key indicator of oxidative damage to cell membrane integrity (Begović et al., 2018). In this study, we found that H2O2 and MDA increased significantly under drought stress, but spraying FA inhibited this phenomenon. This indicates that drought stress led to the dehydration and large accumulation of H2O2 and MDA in the oat leaves, which ultimately led to the oxidative damage and growth retardation of leaves, while FA mitigated the oxidative damage of cell membranes by direct scavenging of H2O2 (Figures 2A, B) (Irani et al., 2021). The antioxidant enzyme defense system is another important mechanism for plants to manage drought stress, protecting the cell membrane from oxidative damage (Wang et al., 2023b). FA is a broad-spectrum plant growth regulator with pleiotropic effects on a wide range of unfavorable environmental factors (Sun et al., 2020; Liang et al., 2024). In our study, a significant increase in antioxidant enzyme activities was observed under water deficit conditions, and antioxidant enzyme activities were further increased by FA treatment (Figures 2C-G). This indicates that FA exerted free radical detoxification and cell membrane structure repair and protection effects on plant seedlings, This could maintain the stability of the intracellular environment, effectively alleviating the oxidative damage caused by drought stress, and enhancing the drought resistance of plant seedlings. Similar results were observed in cabbage seedlings under calcium nitrate stress (Wu et al., 2023). The phenylpropanoid biosynthesis pathway not only affects plant growth and development but also stress responses. Dong and Lin (2021) found that the plant phenylpropanoid metabolic pathway, especially lignin synthesis, has an important regulatory function in plant responses to biotic and abiotic stresses. Phenylpropanes are a large class of plant secondary metabolites. In the phenylpropanoid biosynthesis pathway, phenylalanine is deaminated by PAL to produce cinnamic acid, which is then converted to p-coumaric acid. Finally, p-coumaric acid is converted to p-coumaroyl coenzyme A catalyzed by 4-coumaroyl coenzyme A ligase (4CL), and p-coumarin and p-coumarol coenzyme A are ultimately formed into lignin through various metabolic pathways (Deng and Lu, 2017). In this study, we found that PAL genes were upregulated following the spraying of FA under drought stress, indicating that the biosynthesis of coumaric acid and ferulic acid should be higher (Chen et al., 2021). Interestingly, this study found that the levels of metabolites on coumaric acid and ferulic acid decreased under the action of FA, and more lignin synthase genes were also upregulated, such as 4CL. Therefore, it is speculated that more p-coumaric acid and ferulic acid may be converted into lignin monomer sinapyl alcohol under the catalysis of 4CL, which promotes the synthesis of lignin (Zhang et al., 2023). The increase in sinapyl alcohol in this study also confirms this phenomenon. Peroxidase (POD) is a multi-functional enzyme. It is involved in several different plant physiological processes, including stress resistance, oxidation, and the polymerization of lignin monomers after transportation to the cell wall (Lee et al., 2021). Under stress, plants enhance stress resistance by up-regulating the expression of genes encoding POD (Fei et al., 2020; Li et al., 2021). Research has shown that POD is a hub gene related to lignin biosynthesis (Yang et al., 2023). In this study, the upregulated expression of the POD gene further indicates that FA can promote lignin synthesis. Lignin is an important cell-wall component, and the increase in the lignin level contributes to the fixation of the cell wall (Shah et al., 2023). Plant lignification has been reported to improve stress resistance (Chun et al., 2019) and drought (Gu et al., 2020). Therefore, it is speculated that FA can help plants increase the thickness and strength of cell walls by promoting lignin synthesis, which controls water losses, reduces the wilting degree of plants, and enhances the drought resistance of oats (Ibrahim et al., 2019). The changes in PAL and POD activities in the physiological indexes further confirmed the possibility of lignin accumulation in plants.
Glutathione metabolism is known to be involved in the maintenance of cellular redox homeostasis under drought stress. This adaptive mechanism repairs oxidative damage by utilizing a variety of reactive oxygen scavengers (antioxidants) and redox reactions (Raihan et al., 2023; Wang et al., 2023c). Glutathione metabolism-related genes and metabolites regulate cellular redox homeostasis during abiotic stress (Yang et al., 2022; Horváth et al., 2023; Ali et al., 2024). For example, the overexpression of glutathione metabolism-related genes and proteins, including GST and G6PDH, improved drought tolerance in transgenic plants (Wang et al., 2016; Wei et al., 2020). The overexpression of GST activates its antioxidant-related transcripts, reduces ROS accumulation, and enhances plant drought resistance (Nerva et al., 2022; Niu et al., 2024). In this study, the genes encoding GST were all upregulated after spraying FA under drought stress, suggesting that FA may improve the scavenging capacity of ROS in oats under drought stress by increasing the expression of GST. This is similar to the results of a previous study on tea plants (Sun et al., 2020). The pentose phosphate pathway produces of NADPH, which promotes the production of reduced glutathione and thus maintains cellular redox homeostasis (Sharkey, 2021). G6PDH, as the rate-limiting enzyme in the pentose phosphate pathway (PPP), plays a key role in maintaining redox homeostasis (Li et al., 2023b). Under drought and cold stress, G6PDH maintains cellular redox homeostasis by increasing NADPH/NADP+ levels (Landi et al., 2016; Zhang et al., 2020). The expression of G6PDH can control the rate of increase of ROS and improve stress tolerance (Jiang et al., 2022). In this study, we found that the expression of the gene encoding G6PDH was downregulated after spraying with FA under drought stress. This suggests that glutathione reduction in FA-treated oat plants does not strictly depend on the positive PPP response to replenish NADPH for drought resistance. In contrast, Sun et al. (2020) showed that FA may enhance the ROS scavenging capacity of tea plants under moderate to severe drought stress by increasing the expression of G6PDH and the biosynthesis of the reduced form of glutathione, This differs from our results, possibly due to the differences in the species or the degree of stress employed in the experiments. Ascorbic acid is a powerful antioxidant that neutralizes free radicals and reduces oxidants, but in this process, ascorbic acid itself is oxidized to dehydroascorbic acid and loses its antioxidant function (Njus et al., 2020). Glutathione is a tripeptide that is able to reduce dehydroascorbic acid and turn it back into ascorbic acid, allowing it to continues to play an antioxidant role (Foyer and Kunert, 2024). In this study, the level of dehydroascorbic acid decreased after spraying FA under drought stress, and the results were consistent with the physiological level. This may be because the FA treatment accelerated the reduction process of glutathione and converted it back to ascorbic acid, which enhanced the antioxidant effect and alleviated the oxidative damage caused by drought stress to oats (Kim et al., 2020). As a ubiquitous amino acid, glutamate is involved in the synthesis of a wide range of precursors and plant resistance-related proteins or non-protein amino acids (Forde and Roberts, 2014). The exogenous application of glutamate is reported to enhances salt tolerance in lettuce (Franzoni et al., 2022). Glutamate treatment under drought stress significantly increases the proline content in the phloem and xylem of oilseed rape (La et al., 2020). Proline synthesis occurs in plants under two pathways, namely, through the glutamate and ornithine pathways. Among them, glutamate is directly or indirectly involved in altering certain physiological metabolisms (e.g., carbon and nitrogen metabolism) in plants by regulating the expression of related genes and the activities of key enzymes, which helps plants to grow under stress and enhances their adaptability to various adversities (Goto et al., 2020; Teixeira et al., 2020). In this study, FA treatment under drought stress increased the level of the metabolite glutamate. It is hypothesized that FA may promote proline synthesis and enhance drought resistance primarily by increasing the glutamate pathway.
5 Conclusion
This study indicates that FA improves drought tolerance and promotes the growth of oats through physiology and gene and metabolite expression under water-deficit conditions. FA reduces cell membrane damage and ROS accumulation, alleviates drought stress, and promotes oat growth by increasing chlorophyll content and enhancing antioxidant defense activity. It also improves the drought resistance of oats by inducing the expression of genes in phenylpropanoid biosynthesis, glutathione metabolism, and metabolite accumulation (Figure 9). Our findings revealed a new direction and theoretical basis for the future research on the use of FA technology to improve crop drought resistance.
Figure 9. Mechanism of fulvic acid regulating drought resistance of oats.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author contributions
SZ: Conceptualization, Data curation, Formal Analysis, Investigation, Methodology, Software, Visualization, Writing – original draft. JM: Data curation, Formal Analysis, Methodology, Visualization, Writing – review & editing. BZ: Formal Analysis, Methodology, Writing – review & editing. ZW: Formal Analysis, Methodology, Writing – review & editing. ZY: Investigation, Methodology, Writing – review & editing. MW: Investigation, Methodology, Writing – review & editing. JL: Conceptualization, Data curation, Funding acquisition, Methodology, Project administration, Supervision, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by project of the National Natural Science Foundation (32160508); Oats whole industry chain technology innovation team (BR22-12-05); Evaluation of oat resistance germplasm resources, breeding and utilization of new varieties (BR22-11-15); Breeding of new forage oat varieties and demonstration application of improved variety breeding (CCPTZX2023B05); The National Modern Agricultural Industry Technology System (CARS-07); Oat engineering laboratory capacity building project in Inner Mongolia Autonomous Region (BR221023).
Conflict of interest
Author ZW was employed by National Center of Pratacultural Technology Innovation/M-Grass Ecology And Environment Group Co., Ltd.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2024.1439747/full#supplementary-material
References
Alemayehu, G. F., Forsido, S. F., Tola, Y. B., Amare, E., Capurso, C. (2023). Nutritional and phytochemical composition and associated health benefits of oat (Avena sativa) grains and oat-based fermented food products. Sci. World J. 2023, 1–16. doi: 10.1155/2023/2730175
Ali, Q., Sami, A., Haider, M. Z., Ashfaq, M., Javed, M. A. (2024). Antioxidant production promotes defense mechanism and different gene expression level in Zea mays under abiotic stress. Sci. Rep. 14, 7114. doi: 10.1038/s41598-024-57939-6
Bai, J., Jin, K., Qin, W., Wang, Y., Yin, Q. (2021). Proteomic responses to alkali stress in oats and the alleviatory effects of exogenous spermine application. Front. Plant Sci. 12. doi: 10.3389/fpls.2021.627129
Begović, L., Mlinarić, S., Dunić, J. A., Cesar, V. (2018). Induction of oxidative stress in barley by drought and chilling. Free Radical Biol. Med. 124, 566. doi: 10.1016/j.freeradbiomed.2018.05.033
Bielczynski, L. W., Schansker, G., Croce, R. (2020). Consequences of the reduction of the Photosystem II antenna size on the light acclimation capacity of Arabidopsis thaliana. Plant Cell Environ. 43, 866–879. doi: 10.1111/pce.13701
Bocanegra, M. P., Lobartini, J. C., Orioli, G. A. (2006). Plant uptake of iron chelated by humic acids of different molecular weights. Commun. Soil Sci. Plant Anal. 37, 239–248. doi: 10.1080/00103620500408779
Bolger, A. M., Lohse, M., Usadel, B. (2014). Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120. doi: 10.1093/bioinformatics/btu170
Che, R., Huang, L., Xu, J.-W., Zhao, P., Li, T., Ma, H., et al. (2017). Effect of fulvic acid induction on the physiology, metabolism, and lipid biosynthesis-related gene transcription of Monoraphidium sp. FXY-10. Bioresour. Technol. 227, 324–334. doi: 10.1016/j.biortech.2016.12.017
Chen, L., Yang, W., Liu, S., Meng, Y., Zhu, Z., Liang, R., et al. (2023). Genome-wide analysis and identification of light-harvesting chlorophyll a/b binding (LHC) gene family and BSMV-VIGS silencing TaLHC86 reduced salt tolerance in wheat. Int. J. Biol. Macromol. 242, 124930. doi: 10.1016/j.ijbiomac.2023.124930
Chen, O., Deng, L., Ruan, C., Yi, L., Zeng, K. (2021). Pichia galeiformis induces resistance in postharvest citrus by activating the phenylpropanoid biosynthesis pathway. J. Agr. Food Chem. 69, 2619–2631. doi: 10.1021/acs.jafc.0c06283
Chen, Q., Qu, Z., Ma, G., Wang, W., Dai, J., Zhang, M., et al. (2022a). Humic acid modulates growth, photosynthesis, hormone and osmolytes system of maize under drought conditions. Agr. Water Manage. 263, 107447. doi: 10.1016/j.agwat.2021.107447
Chen, X., Zhang, X., Chen, H., Xu, X. (2022b). Physiology and proteomics reveal Fulvic acid mitigates Cadmium adverse effects on growth and photosynthetic properties of lettuce. Plant Sci. 323, 111418. doi: 10.1016/j.plantsci.2022.111418
Chun, H. J., Baek, D., Cho, H. M., Lee, S. H., Jin, B. J., Yun, D. J., et al. (2019). Lignin biosynthesis genes play critical roles in the adaptation of Arabidopsis plants to high-salt stress. Plant Signaling Behav. 14, 1625697. doi: 10.1080/15592324.2019.1625697
Cui, X., Wang, B., Chen, Z., Guo, J., Zhang, T., Zhang, W., et al. (2023). Comprehensive physiological, transcriptomic, and metabolomic analysis of the key metabolic pathways in millet seedling adaptation to drought stress. Physiol. Plant 175, e14122. doi: 10.1111/ppl.14122
Deng, Y., Lu, S. (2017). Biosynthesis and regulation of phenylpropanoids in plants. Crit. Rev. Plant Sci. 36, 257–290. doi: 10.1080/07352689.2017.1402852
Dong, N. Q., Lin, H. X. (2021). Contribution of phenylpropanoid metabolism to plant development and plant–environment interactions. J. Integr. Plant Biol. 63, 180–209. doi: 10.1111/jipb.13054
Dwivedi, A. K., Singh, V., Anwar, K., Pareek, A., Jain, M. (2023). Integrated transcriptome, proteome and metabolome analyses revealed secondary metabolites and auxiliary carbohydrate metabolism augmenting drought tolerance in rice. Plant Physiol. Biochem. 201, 107849. doi: 10.1016/j.plaphy.2023.107849
Elakhdar, A., El-Naggar, A. A., Kubo, T., Kumamaru, T. (2023). Genome-wide transcriptomic and functional analyses provide new insights into the response of spring barley to drought stress. Physiol. Plant 175, e14089. doi: 10.1111/ppl.14089
Farooq, M., Wahid, A., Lee, D. J. (2009). Exogenously applied polyamines increase drought tolerance of rice by improving leaf water status, photosynthesis and membrane properties. Acta Physiol. Plant 31, 937–945. doi: 10.1007/s11738-009-0307-2
Fei, X., Li, J., Kong, L., Hu, H., Tian, J., Liu, Y., et al. (2020). miRNAs and their target genes regulate the antioxidant system of Zanthoxylum bungeanum under drought stress. Plant Physiol. Biochem. 150, 196–203. doi: 10.1016/j.plaphy.2020.01.040
Forde, B. G., Roberts, M. R. (2014). Glutamate receptor-like channels in plants: a role as amino acid sensors in plant defence? F1000Prime Rep. 6, 37. doi: 10.12703/p6-37
Foyer, C. H., Kunert, K. (2024). The ascorbate/glutathione cycle coming of age. J. Exp. Bot. 75, 2682–2699. doi: 10.1093/jxb/erae023
Franzoni, G., Cocetta, G., Trivellini, A., Garabello, C., Contartese, V., Ferrante, A. (2022). Effect of exogenous application of salt stress and glutamic acid on lettuce (Lactuca sativa L.). Sci. Hortic. 299, 111027. doi: 10.1016/j.scienta.2022.111027
Gao, W., Zhang, Y., Feng, Z., Bai, Q., He, J., Wang, Y. (2018). Effects of melatonin on antioxidant capacity in naked oat seedlings under drought stress. Molecules 23, 1580. doi: 10.3390/molecules23071580
Gong, G.-q., Yuan, X., Zhang, Y. j., Li, Y. j., Liu, W. x., Wang, M., et al. (2020). Characterization of coal-based fulvic acid and the construction of a fulvic acid molecular model. RSC Adv. 10, 5468–5477. doi: 10.1039/c9ra09907g
Goto, Y., Maki, N., Ichihashi, Y., Kitazawa, D., Igarashi, D., Kadota, Y., et al. (2020). Exogenous treatment with glutamate induces immune responses in Arabidopsis. Mol. Plant-Microbe Interact. 33, 474–487. doi: 10.1094/mpmi-09-19-0262-r
Gu, H., Wang, Y., Xie, H., Qiu, C., Zhang, S., Xiao, J., et al. (2020). Drought stress triggers proteomic changes involving lignin, flavonoids and fatty acids in tea plants. Sci. Rep. 10, 15504. doi: 10.1038/s41598-020-72596-1
Hareem, M., Danish, S., Obaid, S. A., Ansari, M. J., Datta, R. (2024). Mitigation of drought stress in chili plants (Capsicum annuum L.) using mango fruit waste biochar, fulvic acid and cobalt. Sci. Rep. 14, 14270. doi: 10.1038/s41598-024-65082-5
Horváth, E., Kulman, K., Tompa, B., Hajnal, Á.B., Pelsőczi, A., Bela, K., et al. (2023). Glutathione transferases are involved in the genotype-specific salt-stress response of tomato plants. Antioxidants 12, 1682. doi: 10.3390/antiox12091682
Hu, J., Zhao, X., Gu, L., Liu, P., Zhao, B., Zhang, J., et al. (2023). The effects of high temperature, drought, and their combined stresses on the photosynthesis and senescence of summer maize. Agr. Water Manage. 289, 108525. doi: 10.1016/j.agwat.2023.108525
Ibrahim, W., Zhu, Y. M., Chen, Y., Qiu, C. W., Zhu, S., Wu, F. (2019). Genotypic differences in leaf secondary metabolism, plant hormones and yield under alone and combined stress of drought and salinity in cotton genotypes. Physiol. Plant 165, 343–355. doi: 10.1111/ppl.12862
Irani, H., ValizadehKaji, B., Naeini, M. R. (2021). Biostimulant-induced drought tolerance in grapevine is associated with physiological and biochemical changes. Chem. Biol. Technol. Agric. 8, 1–13. doi: 10.1186/s40538-020-00200-9
Javadi, T., Rohollahi, D., Ghaderi, N., Nazari, F. (2017). Mitigating the adverse effects of drought stress on the morpho-physiological traits and anti-oxidative enzyme activities of Prunus avium through β-amino butyric acid drenching. Sci. Hortic. 218, 156–163. doi: 10.1016/j.scienta.2017.02.019
Jesmin, A., Anh, L. H., Mai, N. P., Khanh, T. D., Xuan, T. D. (2023). Fulvic acid improves salinity tolerance of rice seedlings: evidence from phenotypic performance, relevant phenolic acids, and momilactones. Plants 12, 2359. doi: 10.3390/plants12122359
Jiang, Y., Su, S., Chen, H., Li, S., Shan, X., Li, H., et al. (2023). Transcriptome analysis of drought-responsive and drought-tolerant mechanisms in maize leaves under drought stress. Physiol. Plant 175, e13875. doi: 10.1111/ppl.13875
Jiang, Z., Wang, M., Nicolas, M., Ogé, L., Pérez-Garcia, M. D., Crespel, L., et al. (2022). Glucose-6-phosphate dehydrogenases: the hidden players of plant physiology. Int. J. Mol. Sci. 23, 16128. doi: 10.3390/ijms232416128
Justi, M., Morais, E. G., Silva, C. A. (2019). Fulvic acid in foliar spray is more effective than humic acid via soil in improving coffee seedlings growth. Arch. Agron. Soil Sci. 65, 1969–1983. doi: 10.1080/03650340.2019.1584396
Kim, D., Langmead, B., Salzberg, S. L. (2015). HISAT: a fast spliced aligner with low memory requirements. Nat. Methods 12, 357–360. doi: 10.1038/nmeth.3317
Kim, Y.-S., Park, S.-I., Kim, J.-J., Boyd, J. S., Beld, J., Taton, A., et al. (2020). Expression of heterologous OsDHAR gene improves glutathione (GSH)-dependent antioxidant system and maintenance of cellular redox status in Synechococcus elongatus PCC 7942. Front. Plant Sci. 11. doi: 10.3389/fpls.2020.00231
La, V. H., Lee, B.-R., Islam, M. T., Mamun, M. A., Park, S.-H., Bae, D.-W., et al. (2020). Characterization of glutamate-mediated hormonal regulatory pathway of the drought responses in relation to proline metabolism in Brassica napus L. Plants 9, 512. doi: 10.3390/plants9040512
Lalas, S., Athanasiadis, V., Dourtoglou, V. G. (2018). Humic and fulvic acids as potentially toxic metal reducing agents in water. Clean Soil Air Water 46, 1700608. doi: 10.1002/clen.201700608
Landi, S., Nurcato, R., De Lillo, A., Lentini, M., Grillo, S., Esposito, S. (2016). Glucose-6-phosphate dehydrogenase plays a central role in the response of tomato (Solanum lycopersicum) plants to short and long-term drought. Plant Physiol. Biochem. 105, 79–89. doi: 10.1016/j.plaphy.2016.04.013
Lee, C.-J., Park, S.-U., Kim, S.-E., Lim, Y.-H., Ji, C. Y., Kim, Y.-H., et al. (2021). Overexpression of IbLfp in sweetpotato enhances the low-temperature storage ability of tuberous roots. Plant Physiol. Biochem. 167, 577–585. doi: 10.1016/j.plaphy.2021.08.041
Li, X., Cai, Q., Yu, T., Li, S., Li, S., Li, Y., et al. (2023b). ZmG6PDH1 in glucose-6-phosphate dehydrogenase family enhances cold stress tolerance in maize. Front. Plant Sci. 14. doi: 10.3389/fpls.2023.1116237
Li, Y. H., Liu, J. H., Zhao, B. P., Tian, L., Wang, Q., Zhu, S. S., et al. (2019). Regulatory effect of humic acid on chlorophyll fluorescence characteristics of oat under drought stress. J. Irrig. Drain. 39, 1385–1391. doi: 10.13522/j.cnki.ggps.2019275
Li, C., Qi, Y., Zhao, C., Wang, X., Zhang, Q. (2021). Transcriptome profiling of the salt stress response in the leaves and roots of halophytic Eutrema salsugineum. Front. Genet. 12. doi: 10.3389/fgene.2021.770742
Li, C., Wan, Y., Shang, X., Fang, S. (2023a). Integration of transcriptomic and metabolomic analysis unveils the response mechanism of sugar metabolism in Cyclocarya paliurus seedlings subjected to PEG-induced drought stress. Plant Physiol. Biochem. 201, 107856. doi: 10.1016/j.plaphy.2023.107856
Li, X., Zhang, W., Niu, D., Liu, X. (2024). Effects of abiotic stress on chlorophyll metabolism. Plant Sci. 342, 112030. doi: 10.1016/j.plantsci.2024.112030
Liang, Y., Wang, J., Wang, Z., Hu, D., Jiang, Y., Han, Y., et al. (2024). Fulvic acid alleviates the stress of low nitrogen on maize by promoting root development and nitrogen metabolism. Physiol. Plant 176, e14249. doi: 10.1111/ppl.14249
Liu, C., Lv, C., Ai, X., Bi, H. (2022). Effects of fulvic acid on photosynthetic characteristics, yield and quality of cucumber under drought stress. Chin. J. Appl. Ecol. 33, 1300–1310. doi: 10.13287/j.1001-9332.202205.014
Muhammad, I., Shalmani, A., Ali, M., Yang, Q.-H., Ahmad, H., Li, F. B. (2021). Mechanisms regulating the dynamics of photosynthesis under abiotic stresses. Front. Plant Sci. 11. doi: 10.3389/fpls.2020.615942
Nerva, L., Guaschino, M., Pagliarani, C., De Rosso, M., Lovisolo, C., Chitarra, W. (2022). Spray-induced gene silencing targeting a glutathione S-transferase gene improves resilience to drought in grapevine. Plant Cell Environ. 45, 347–361. doi: 10.1111/pce.14228
Nikoogoftar-Sedghi, M., Rabiei, V., Razavi, F., Molaei, S., Khadivi, A. (2024). Fulvic acid foliar application: a novel approach enhancing antioxidant capacity and nutritional quality of pistachio (Pistacia vera L.). BMC Plant Biol. 24, 241. doi: 10.1186/s12870-024-04974-0
Niu, M. X., Feng, C. H., He, F., Zhang, H., Bao, Y., Liu, S. J., et al. (2024). The miR6445-NAC029 module regulates drought tolerance by regulating the expression of glutathione S-transferase U23 and reactive oxygen species scavenging in Populus. New Phytol. 242, 2043–2058. doi: 10.1111/nph.19703
Njus, D., Kelley, P. M., Tu, Y.-J., Schlegel, H. B. (2020). Ascorbic acid: The chemistry underlying its antioxidant properties. Free Radical Biol. Med. 159, 37–43. doi: 10.1016/j.freeradbiomed.2020.07.013
Peng, X., Feng, C., Wang, Y. T., Zhang, X., Wang, Y. Y., Sun, Y. T., et al. (2022). miR164g-MsNAC022 acts as a novel module mediating drought response by transcriptional regulation of reactive oxygen species scavenging systems in apple. Hortic. Res. 9, uhac192. doi: 10.1093/hr/uhac192
Raihan, M. R. H., Rahman, M., Rastogi, A., Fujita, M., Hasanuzzaman, M. (2023). Exogenous allantoin confers rapeseed (Brassica campestris) tolerance to simulated drought by improving antioxidant metabolism and physiology. Antioxidants 12, 1508. doi: 10.3390/antiox12081508
Salih, W., Epule, T. E., EL Khalki, E. M., Ouatiki, H., Erraki, S., Achli, S., et al. (2023). A comprehensive assessment of satellite precipitation products over a semi-arid region: focus on extreme events. Nat. Hazards 120, 3037–3065. doi: 10.1007/s11069-023-06317-y
Shah, R. S., Senanayake, M., Zhang, H.-H., Pu, Y., Biswal, A. K., Pingali, S. V., et al. (2023). Evidence for lignin–carbohydrate complexes from studies of transgenic switchgrass and a model lignin–pectin composite. ACS Sustain. Chem. Eng. 11, 15941–15950. doi: 10.1021/acssuschemeng.3c04322
Sharkey, T. D. (2021). Pentose phosphate pathway reactions in photosynthesizing cells. Cells 10, 1547. doi: 10.3390/cells10061547
Shokhmgar, M., Seghatoleslami, M., Moosavi, S. G., Baradaran, R. (2023). Plant growth regulators can affect the physiological traits of foxtail millet (Setaria italica) under deficit irrigation. Gesunde Pflanz. 75, 2505–2514. doi: 10.3390/cells10061547
Singh, A., Kumari, A., Chauhan, A. K. (2022). Formulation and evaluation of novel functional snack bar with amaranth, rolled oat, and unripened banana peel powder. J. Food Sci. Tech. 59, 3511–3521. doi: 10.1007/s13197-021-05344-6
Sun, J., Qiu, C., Ding, Y., Wang, Y., Sun, L., Fan, K., et al. (2020). Fulvic acid ameliorates drought stress-induced damage in tea plants by regulating the ascorbate metabolism and flavonoids biosynthesis. BMC Genomics 21, 411. doi: 10.1186/s12864-020-06815-4
Teixeira, W. F., Soares, L. H., Fagan, E. B., da Costa Mello, S., Reichardt, K., Dourado-Neto, D. (2020). Amino acids as stress reducers in soybean plant growth under different water-deficit conditions. J. Plant Growth Regul. 39, 905–919. doi: 10.1007/s00344-019-10032-z
Tian, H., Zhou, Q., Liu, W., Zhang, J., Chen, Y., Jia, Z., et al. (2022). Responses of photosynthetic characteristics of oat flag leaf and spike to drought stress. Front. Plant Sci. 13. doi: 10.3389/fpls.2022.917528
Utsumi, Y., Utsumi, C., Tanaka, M., Ha, C. V., Takahashi, S., Matsui, A., et al. (2019). Acetic acid treatment enhances drought avoidance in cassava (Manihot esculenta Crantz). Front. Plant Sci. 10. doi: 10.3389/fpls.2019.00521
Walker, B. J., Drewry, D. T., Slattery, R. A., VanLoocke, A., Cho, Y. B., Ort, D. R. (2018). Chlorophyll can be reduced in crop canopies with little penalty to photosynthesis. Plant Physiol. 176, 1215–1232. doi: 10.1104/pp.17.01401
Wang, W., Shi, S., Kang, W., He, L. (2023b). Enriched endogenous free Spd and Spm in alfalfa (Medicago sativa L.) under drought stress enhance drought tolerance by inhibiting H2O2 production to increase antioxidant enzyme activity. J. Plant Physiol. 291, 154139. doi: 10.1016/j.jplph.2023.154139
Wang, X., Wang, M., Yan, G., Yang, H., Wei, G., Shen, T., et al. (2023c). Comparative analysis of drought stress-induced physiological and transcriptional changes of two black sesame cultivars during anthesis. Front. Plant Sci. 14. doi: 10.3389/fpls.2023.1117507
Wang, H., Yang, L., Li, Y., Hou, J., Huang, J., Liang, W. (2016). Involvement of ABA-and H2O2-dependent cytosolic glucose-6-phosphate dehydrogenase in maintaining redox homeostasis in soybean roots under drought stress. Plant Physiol. Biochem. 107, 126–136. doi: 10.1016/j.plaphy.2016.05.040
Wang, C., Yang, H., Yue, L., Sun, W., Chen, F., Cao, X., et al. (2022a). Physiological and molecular level understanding of advanced carbon dots to enhance maize drought tolerance: modulation of photosynthesis and signaling molecules. Environ. Sci.: Nano 9, 3821–3832. doi: 10.1039/d2en00176d
Wang, Y., Yang, R., Zheng, J., Shen, Z., Xu, X. (2019). Exogenous foliar application of fulvic acid alleviate cadmium toxicity in lettuce (Lactuca sativa L.). Ecotox. Environ. Safe. 167, 10–19. doi: 10.1016/j.ecoenv.2018.08.064
Wang, J., Yu, Y., Jiang, C., Sun, Z., Wang, X., Wang, Z., et al. (2023a). Comparative analysis of physiology-anatomy and transcriptome-metabolome involving acute drought stress response of root between two distinct peanut cultivars at seedling stage. Environ. Exp. Bot. 214, 105442. doi: 10.1016/j.envexpbot.2023.105442
Wei, B., Hou, K., Zhang, H., Wang, X., Wu, W. (2020). Integrating transcriptomics and metabolomics to studies key metabolism, pathways and candidate genes associated with drought-tolerance in Carthamus tinctorius L. under drought stress. Ind. Crops Prod. 151, 112465. doi: 10.1016/j.indcrop.2020.112465
Wu, X., Zhang, Y., Chu, Y., Yan, Y., Wu, C., Cao, K., et al. (2023). Fulvic acid alleviates the toxicity induced by calcium nitrate stress by regulating antioxidant and photosynthetic capacities and nitrate accumulation in Chinese flowering cabbage seedlings. Appl. Sci. 13, 12373. doi: 10.3390/app132212373
Yang, S., Liu, M., Chu, N., Chen, G., Wang, P., Mo, J., et al. (2022). Combined transcriptome and metabolome reveal glutathione metabolism plays a critical role in resistance to salinity in rice landraces HD961. Front. Plant Sci. 13. doi: 10.3389/fpls.2022.952595
Yang, J., Song, J., Feng, Y., Cao, Y., Fu, B., Zhang, Z., et al. (2023). Osmotic stress-induced lignin synthesis is regulated at multiple levels in alfalfa (Medicago sativa L.). Int. J. Biol. Macromol. 246, 125501. doi: 10.1016/j.ijbiomac.2023.125501
Yu, B., Xue, X., Nie, P., Lu, N., Wang, L. (2024). Fulvic acid alleviates cadmium-induced root growth inhibition by regulating antioxidant enzyme activity and carbon–nitrogen metabolism in apple seedlings. Front. Plant Sci. 15. doi: 10.3389/fpls.2024.1370637
Zhang, X., Liu, W., Lv, Y., Li, T., Tang, J., Yang, X., et al. (2022). Effects of drought stress during critical periods on the photosynthetic characteristics and production performance of Naked oat (Avena nuda L.). Sci. Rep. 12, 11199. doi: 10.1038/s41598-022-15322-3
Zhang, Y., Luo, M., Cheng, L., Lin, Y., Chen, Q., Sun, B., et al. (2020). Identification of the cytosolic glucose-6-phosphate dehydrogenase gene from strawberry involved in cold stress response. Int. J. Mol. Sci. 21, 7322. doi: 10.3390/ijms21197322
Zhang, F., Wu, J., Sade, N., Wu, S., Egbaria, A., Fernie, A. R., et al. (2021). Genomic basis underlying the metabolome-mediated drought adaptation of maize. Genome Biol. 22, 260. doi: 10.1186/s13059-021-02481-1
Zhang, X., Yang, R., Zheng, X., Wang, Q., Silvy, E. M., Li, Y., et al. (2023). UV-C irradiation accelerated the deposition of suberin and lignin at wounds by activating phenylpropanoid metabolism in potato tubers. Sci. Hortic. 309, 111634. doi: 10.1016/j.scienta.2022.111634
Zhao, C., Yang, M., Wu, X., Wang, Y., Zhang, R. (2021). Physiological and transcriptomic analyses of the effects of exogenous melatonin on drought tolerance in maize (Zea mays L.). Plant Physiol. Biochem. 168, 128–142. doi: 10.1016/j.plaphy.2021.09.044
Keywords: oat, drought stress, fulvic acid, phenylpropanoid biosynthesis, glutathione metabolism
Citation: Zhu S, Mi J, Zhao B, Wang Z, Yang Z, Wang M and Liu J (2024) Integrative transcriptome and metabolome analysis reveals the mechanism of fulvic acid alleviating drought stress in oat. Front. Plant Sci. 15:1439747. doi: 10.3389/fpls.2024.1439747
Received: 28 May 2024; Accepted: 30 August 2024;
Published: 19 September 2024.
Edited by:
PingWu Liu, Hainan University, ChinaReviewed by:
Zhen Yang, Qilu University of Technology, ChinaShoukun Dong, Northeast Agricultural University, China
Copyright © 2024 Zhu, Mi, Zhao, Wang, Yang, Wang and Liu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jinghui Liu, Y2F1bGloQDE2My5jb20=