- Research
- Open access
- Published:
Transcriptomic analysis of melatonin-mediated drought stress response genes in alfalfa during germination period
BMC Plant Biology volume 25, Article number: 637 (2025)
Abstract
Background
Drought stress is a predominant abiotic factor contributing to reduced crop yields globally. Therefore, exploring the molecular mechanism of drought control is of great significance to improve drought resistance and ultimately achieve crop yield increase. As a plant endogenous hormone, melatonin plays a key role in the regulation of abiotic stress, but the key genes and metabolic pathways of melatonin mediated drought resistance regulation in alfalfa have not been fully revealed. Based on transcriptomics and physiological index detection, this study aimed to explore the regulatory mechanism of melatonin in alleviating drought stress during alfalfa germination.
Results
The findings revealed that alfalfa seedlings treated with melatonin exhibited higher germination rates, increased shoot length, and greater fresh weight compared to those exposed solely to drought stress. Additionally, there was a reduction in the levels of malondialdehyde (MDA) and superoxide anion (O2−), while the activity and concentration of superoxide dismutase (SOD), peroxidase (POD), and glutathione (GSH) were enhanced to varying extents. To investigate the molecular mechanism underlying melatonin-mediated drought resistance in alfalfa, we performed a transcriptomic analysis on the seedlings. In the drought treatment group, we identified a total of 1,991 differentially expressed genes (DEGs), comprising 778 up-regulated and 1,213 down-regulated genes. Conversely, in the melatonin-treated group, we discovered 2,336 DEGs, including 882 up-regulated and 1,454 down-regulated genes.
Conclusions
Through the application of GO functional annotation and KEGG pathway enrichment analysis, we discovered that DEGs were predominantly enriched in pathways related to flavonoid and isoflavone biosynthesis, plant hormone biosynthesis and signal transduction, glutathione metabolism, and MAPK signaling, and the ABC transporter signaling. Notably, the DEGs added to the MT group showed greater enrichment in these pathways. This suggests that MT mitigates drought stress by modulating the expression of genes associated with energy supply and antioxidant capacity. These findings hold significant reference value for breeding drought-tolerant alfalfa and other crops.
Introduction
Currently, soil drought is intensifying due to human activities and climate change. Drought stress is one of the most prevalent abiotic challenges in crop production, restricting plant growth and yield by limiting water absorption and utilization [1,2,3]. The impact of drought on crops is multifaceted, affecting photosynthesis, respiration, nutrient absorption, and transportation. The disruption of these physiological processes can result in delayed growth and development, and even plant death, severely impacting agricultural production and hindering local economic development [4,5,6].
Alfalfa (Medicago sativa L.) is a high-yield forage known for its substantial protein content and balanced amino acid composition. In China, it is predominantly cultivated in arid and semi-arid regions, including the northwest and northern areas. However, alfalfa in these regions frequently encounters drought, a condition that not only hampers its growth but also adversely impacts critical physiological processes such as photosynthesis. Research indicates that drought stress is a significant factor curtailing alfalfa productivity [7, 8]. Thus, exploring the molecular mechanisms underpinning drought resistance in alfalfa and enhancing its resilience is vital for increasing yield. Current research on drought stress in alfalfa primarily focuses on identifying specific drought-resistant genes, such as MsCYP71 and MsTHI1 [9, 10]. Nevertheless, these studies often limit themselves to the effects of one or a few genes, indicating the need for more comprehensive exploration. In combating drought, plant hormones play an essential role by regulating growth, development, and metabolism, thereby aiding adaptation to adverse environments [11, 12]. Melatonin, a multifunctional biomolecule, has recently been identified as playing a pivotal role in regulating growth, enhancing antioxidant capacity, and facilitating mineral absorption in plants [13,14,15]. As an antioxidant, melatonin stimulates the expression of antioxidant enzyme genes and increases enzyme activity, thereby bolstering plant stress resistance [16]. Although melatonin's role in enhancing drought resilience has been widely acknowledged across various crops, the precise transcriptomic mechanisms underlying its effects on drought tolerance in alfalfa remain largely unexplored [17, 18]. This study aims to address this critical knowledge gap by examining the gene expression changes triggered by melatonin treatment under drought stress. Previous research has demonstrated that, in Zea mays melatonin promotes the upregulation of antioxidant genes such as SOD and POD, thereby mitigating oxidative damage [19]. Similarly, in Oryza sativa, melatonin interacts with ABA signaling pathways to induce stomatal closure, improving drought tolerance [20]. These findings suggest that comparable regulatory mechanisms may also exist in alfalfa; however, they have yet to be thoroughly investigated. To bridge this gap, this study explores the gene expression changes associated with melatonin treatment during drought stress. To achieve a more comprehensive analysis, Weighted Gene Co-expression Network Analysis (WGCNA) was employed to construct gene expression-related networks. This method allows for the identification of gene modules strongly associated with specific traits, such as drought response, and facilitates the pinpointing of hub genes within the regulatory network. These insights provide a robust framework for systematically analyzing the multi-gene collaborative mechanisms involved in drought tolerance. High-throughput techniques, leveraging omics analysis, can elucidate the regulatory mechanisms of these complex plant traits. In our study, we compared the growth performance and physiological indicators of alfalfa under PEG-simulated drought stress, both with and without melatonin (MT) treatment, and performed whole-gene transcriptome sequencing to identify differentially expressed genes under the influence of MT and PEG. Furthermore, we utilized WGCNA to pinpoint key signaling pathways and hub genes regulated by MT, aiming to clarify their cellular mechanisms of action. Under drought conditions, the synthesis of vital amino acids (such as proline, tryptophan, and arginine), flavonoids (including anthocyanins, quercetin, and licorice), and plant hormones (such as abscisic acid, gibberellin, and melatonin) increases. These findings offer new insights into cultivating drought-tolerant alfalfa.
Result
Effects of melatonin on alfalfa plant growth under drought stress
Various concentrations of melatonin were tested on alfalfa seeds subjected to drought stress. In comparison to the control (CK), a 15% PEG solution significantly reduced the germination rate, seedling fresh weight, and root length of alfalfa seeds (P < 0.05). However, applying low concentrations of melatonin led to an increase in these metrics to varying degrees (Fig. 1 A). Conversely, a higher concentration of melatonin (> 200 μM) significantly diminished these improvements, sometimes even leading to a decline. Under treatment with a 10 μM melatonin solution, alfalfa seeds achieved a germination rate of 78.89% and a root length of 6.21 cm, showing no significant difference from the CK group. However, the fresh weight, recorded at 0.201 g, exceeded that of other concentrations (Fig. 1 B-D). Thus, 10 μM melatonin was determined to be the optimal concentration for mitigating drought stress during the germination of alfalfa seeds.
A Phenotypic display of alfalfa seedlings, with a scale of 2 cm. B Statistics of seed germination rate after 7 days of treatment. The line graph displays the mean and standard deviation of six biological replicates, each containing 30 seeds. C Root length measurement results; The columns with different lowercase letters in the figure represent the significant differences between treatments at the p < 0.05 level according to Duncan's test. The data is presented as the average and standard deviation of three biological replicates, each replicate consisting of 10 plants. D Fresh weight of alfalfa seedlings; This column of data represents the average of three biological replicates. Different lowercase letters indicate significant differences between the treatments identified through Duncan's test. CK represents the control group
Effects of exogenous melatonin pretreatment on antioxidant properties of alfalfa seedlings under drought stress
To further investigate the physiological changes associated with melatonin-mediated drought resistance in alfalfa seedlings, we examined various antioxidant indicators, including superoxide anion, malondialdehyde (MDA), glutathione (GSH), peroxidase (POD), and superoxide dismutase (SOD). Our findings indicate that osmotic stress led to an increase in both superoxide anion and MDA accumulation (Fig. 2 A-E). However, when treated with exogenous melatonin, these levels were reduced by 27.8% and 26.6%, respectively (Fig. 2 A-B). Additionally, while GSH content decreased by 21.8% under osmotic stress, the introduction of exogenous melatonin not only restored but also increased GSH content by 18.2% compared to the control group (Fig. 2 E). Furthermore, under osmotic stress, POD and SOD activities exhibited slight increases, which were further enhanced with exogenous melatonin treatment (Fig. 2 C-D).
Exogenous melatonin can alleviate oxidative stress caused by drought. Specific indicators include O2− level (A), MDA concentration (B), SOD activity (C), POD activity (D), and GSH content (E). Each set of data underwent three biological replicates, with each replicate containing ten independent plant samples. Duncan's test was used for data analysis, with different lowercase letter markers indicating significant differences between treatments
Identification of drought-responsive genes in alfalfa seedlings
A transcriptome analysis was performed on alfalfa seedlings using deionized water, 15% PEG, and a combination of 10 µM MT with 15% PEG as treatments. Through comparison of these treatment methods, a total of 3,134 genes were identified. Only six genes (0.19% of total) were common across all treatments. Additionally, 28 genes were shared between ZMPMT vs. ZMCK and ZMPMT vs. ZMP, comprising only 0.83% of the overall gene count. Eighteen genes were common among ZMPMT, ZMP, and ZMCK, amounting to just 0.57% of the total gene count. In contrast, a significant 1,215 genes were found in all three treatments (ZMPMT, ZMCK, and ZMCK), accounting for approximately 38.77% of the total. Furthermore, 764 genes were unique to the comparison between ZMP and ZMCK, making up about 24.38% of the total gene count (Fig. 3 A). Analyze differentially expressed genes across three treatments, categorize genes with similar or identical expression patterns, and create volcano plots. (Fig. 3 B). The Venn diagram analysis revealed that 778 genes were upregulated in the ZMP vs. ZMCK comparison. In the ZMPMT vs. ZMCK comparison, 882 genes were upregulated, while the ZMPMT vs. ZMP comparison showed a minimal difference, with only 28 genes upregulated (Fig. 3 C). Additionally, variance analysis identified ten differentially expressed genes potentially related to drought tolerance mediated by the melatonin hormone in alfalfa (Fig. 3 D).
Control the comprehensive visualization of PEG, CK, PEG, and MT combined with NaCl (MT + PEG). A The Venn diagram reveals the uniqueness and sharing of DEGs. B The scatter plot illustrates the upregulated and downregulated characteristic genes across each differential group. C Multiple inter group analysis of variance was conducted* Indicates P-value < 0.05, * * indicates P-value < 0.01. D Hierarchical clustering analysis of log2 RPKM expression trends of DEGs treated with PEG or melatonin alone
Enrichment analysis of differentially expressed genes
In the primary classification of GO, DEGs are assigned to three major categories: molecular function, cellular component, and biological process. In the secondary classification, DEGs are more frequently annotated with terms such as binding activity, catalytic activity, cell-binding entity, metabolic process, and cellular process. Notably, the number of DEGs in the ZMP and ZMPMT treatment groups exceeds that in the ZMCK control group (Fig. 4 A). KEGG pathway analysis of DEGs reveals 15 significantly enriched pathways (P < 0.05) in the ZMP versus ZMCK treatment group and 14 significantly enriched pathways in the ZMPMT versus ZMCK treatment group (Fig. 4 B-C). Shared pathways include phenylpropanoid biosynthesis, flavonoid and isoflavone biosynthesis, starch and sucrose metabolism, glutathione metabolism, and MAPK signaling (Fig. 4 D). It is noteworthy that the arginine and proline metabolism pathways are unique to the ZMP versus ZMCK treatment group, whereas the ABC transporter protein pathway is exclusively observed in the ZMPMT versus ZMCK treatment group.
Functional annotation and KEGG enrichment analysis of DEGs were conducted. A The GO annotation analysis categorized the differentially expressed genes into biological processes, cellular components, and molecular functions, represented by distinct colors for group differences. B Enrichment analysis revealed functional classification differences in the comparison between ZMP and ZMCK. C The functional classification and enrichment analysis described the differentially expressed genes between ZMPMT and ZMCK. D Functional classification enrichment analysis was presented for the differentially expressed genes between ZMPMT and ZMP
Gene co-expression network
The response of alfalfa to osmotic stress is regulated by multiple genes. In this study, we performed a transcriptome analysis on Zhongmu-3 under various treatments, identifying numerous differentially expressed genes (DEGs). We utilized Weighted Gene Co-expression Network Analysis (WGCNA) to explore the correlation between DEGs and physiological traits linked to drought stress (Fig. 5 A). We identified twelve co-expressed modules and found that correlation analysis of module traits indicated a significant positive relationship between gene expression levels in the green module and phenotypic characteristics such as germination rate, root length, and fresh weight, with correlation coefficients ranging from 0.767 to 0.867. Conversely, in the blue module, we noted a significant negative correlation, with coefficients ranging from 0.683 to 0.85, among gene expression levels related to similar features. This indicates that genes within the green module may contribute to MT-mediated drought resistance, while those in the blue module may be involved in regulating reactive oxygen species clearance and maintaining oxygen free radical balance, given their significant positive correlation with MDA content and O2− levels, respectively. Furthermore, indicators of antioxidant capacity, such as SOD and POD activity, showed positive correlations with gene expression levels in the yellow module. The gene expression level in the black module is negatively correlated with GSH content. The green module includes 564 DEGs, the blue module comprises 2,748 DEGs, the yellow module contains 970 DEGs, and the black module includes 132 DEGs (Fig. 5 B). In the blue, green, yellow, and black modules, we identified thirty key genes representing corresponding functions (Fig. 5 C-F). These central genes were visualized using CytoScape software version 3.9.1, highlighting their critical roles in flavonoid biosynthesis pathways, plant hormone signaling pathways, soluble sugar metabolism pathways, glutathione metabolism pathways, and ABA transporter regulation.
The WGCNA analysis revealed the correlation between DEGs and physiological traits related to drought stress regulated by melatonin. A The heatmap clearly shows the correlation between gene expression modules and physiological traits related to drought stress regulated by melatonin. B A module classification tree was constructed through hierarchical clustering analysis of WGCNA, demonstrating the similarity and difference between gene expression modules. C In the MEblue module, a detailed co-expression network of DEGs is depicted. D In the MEgreen module, a detailed co-expression network of DEGs is depicted. E In the MEblack module, a detailed co-expression network of DEGs is depicted. F In the Meyellow module, a detailed co-expression network of DEGs is depicted
Transcriptional analysis of several differentially expressed genes
To further confirm the trend of changes in differentially expressed genes (DEGs), qRT-PCR was performed. The eight selected genes were primarily involved in flavonoid synthesis (MsG0180005358.01), ABA biosynthesis and signal transduction (MsG0180001827.01), GA biosynthesis (MsG0180005519.01), and melatonin synthesis (MsG04800200076.01). As illustrated in Fig. 6, the trend of transcript changes in these eight DEGs was consistent with the RNA-seq data, indicating that the data is reproducible and suitable for further analysis. To ensure the reliability of the transcriptome data, we validated it by calculating the linear regression coefficient between the RNA-seq data (log2 (FC)) and the qRT-PCR results (ΔΔCt method). The results showed that for the ZMP vs. ZMCK group, R2 = 0.61, and for the ZMPMT vs. ZMCK group, R2 = 0.70, indicating the consistency between our transcriptome data and the fluorescence quantitative data. The correlation has been illustrated in the attached Figure S1.
Transcriptional analysis of 8 differentially expressed genes. To ensure reproducibility and reliability, qRT-PCR analysis was performed on three independent biological replicates. Calculate the relative expression level using the 2−∆∆Ct method. These columns represent the average of three biological replicates. Using Duncan's test, different lowercase letters indicate significant differences between treatments (p < 0.05)
Discussion
Removal of ROS
In plants, resistance to drought stress is mediated through a complex network of genes, metabolic pathways, and other molecules [21,22,23]. When plants encounter drought stress, they produce a substantial amount of reactive oxygen species (ROS) and harmful substances like malondialdehyde (MDA), leading to increased membrane lipid peroxidation [24, 25]. The plant's elimination of ROS is primarily accomplished through intricate non-enzymatic and enzymatic antioxidant systems [26,27,28]. Non-enzymatic antioxidants mainly include reduced glutathione and osmolytes such as proline, arginine, and soluble sugars, while enzymatic antioxidants include superoxide dismutase, ascorbate peroxidase, glutathione peroxidase, and peroxidase [29,30,31]. Glutathione S-transferase (GST) plays a pivotal role in glutathione metabolism and in maintaining cellular redox homeostasis, contributing significantly to plant resistance against abiotic stress [32]. In our study, 14 differentially expressed genes (DEGs) were identified as enriched in GST. Under drought conditions, these genes showed higher expression levels in the melatonin (MT) treatment group compared to the polyethylene glycol (PEG) treatment group. We hypothesize that melatonin promotes GST expression, accelerates glutathione (GSH) synthesis, and neutralizes ROS. Research indicates that oxidized glutathione (GSSG) can be reverted to GSH through the ascorbic acid-glutathione cycle, maintaining a reduced cellular state and protecting cell membranes from oxidative harm [33, 34]. The upregulation of DHAR, a pivotal gene in the ASA-GSH cycle, observed under melatonin treatment aligns with the increased GSH content (Fig. 2E). This suggests that melatonin enhances the regeneration of glutathione, a critical component for neutralizing ROS and preserving cellular redox balance. Such a mechanism plays a significant role in mitigating oxidative damage caused by drought stress.. Additionally, amino acids and soluble sugars act as accumulated osmolytes and free radical scavengers during stress responses [35, 36]. In our study, the expression levels of crucial genes for proline synthesis, such as ProA and ProB, as well as argF for arginine synthesis, were significantly upregulated in the MT group. Conversely, the main gene for proline degradation, PRODH, was notably downregulated, indicating that melatonin aids in proline accumulation to counter drought stress. In galactose metabolism, we noticed a downregulation of important genes GLA and INV during drought stress. However, their expression was significantly upregulated upon melatonin addition, suggesting that exogenous melatonin may mitigate drought stress by regulating osmotic balance. Furthermore, the expression of genes encoding ROS scavenging enzymes like superoxide dismutase, glutathione peroxidase, ascorbate peroxidase, peroxidase, and catalase, was markedly upregulated under drought conditions. This has significant implications for maintaining ROS homeostasis and enhancing plant drought resistance [37, 38]. Flavonoids, widely present as secondary metabolites in plants, play a crucial role in antioxidant capacity by scavenging free radicals and inhibiting oxidase activity. Studies show that flavonoids can influence the activity of SOD and POD [39,40,41], with upregulation of the anthocyanin reductase (ANR) gene in tobacco enhancing catechin and epicatechin accumulation, thereby boosting the plant's antioxidant capacity [42]. Our transcriptome analysis revealed that under melatonin treatment compared to drought treatment, DEGs associated with flavonoid biosynthesis, including HCT (responsible for catalyzing the synthesis of lignin precursors) and ANR (involved in producing antioxidant catechins), were upregulated in response to melatonin treatment. These metabolites play a pivotal role in scavenging ROS and stabilizing cellular membranes, thereby synergistically enhancing the activities of SOD and POD (Fig. 2 C-D). This combined effect significantly contributes to protecting against oxidative stress caused by drought conditions. This may contribute to increased SOD and POD activity, consistent with our enzyme activity assay results.
Plant hormone regulatory network
Plant hormones are subtle signaling molecules produced within plants, crucial for adapting to adverse environments [43, 44]. Under drought stress, plants synthesize abscisic acid (ABA), which performs multiple functions such as regulating leaf stomatal size, restricting lateral root growth, aiding seed germination, and expediting fruit ripening [45, 46]. In our experiment, we noted a significant upregulation of the key gene NCED (MsG0180001827.01), which is involved in ABA synthesis when subjected to drought stress. This finding suggests that the ABA biosynthesis pathway is activated, enabling plants to manage drought stress by restricting lateral root elongation and regulating stomatal size. Notably, NCED expression diminished following the addition of melatonin. This may occur because melatonin modulates ABA synthesis by inhibiting ABA synthesis gene expression, thus counteracting drought stress [47, 48]. The SnRK2 kinase can induce the expression of the transcription factor ABF2 in alfalfa's ABA signaling pathway. Protein phosphatase PP2C (MsG0080048230.01) inhibits SnRK2 kinase activity via dephosphorylation, while the receptor protein PYL (MsG0480018569.01) counters PP2C activity through ABA. Research indicates that overexpression of the transcription factor ABF2, responsible for synthesizing ABA in alfalfa, can elevate the expression of stress response genes, thereby enhancing the plant's drought resistance [49]. In this study, under PEG stress conditions, PYL and SnRK2 gene expressions were downregulated, whereas ABF2 and PP2C gene expressions increased. This suggests that alfalfa can boost the expression of downstream resistance genes in the root system by enhancing the expression of the transcription factor ABF2, thereby improving the plant's response to drought stress.
Gibberellin (GA) is a trace endogenous signaling molecule synthesized by plants, known for promoting plant cell and stem elongation, leaf enlargement, accelerating growth and development, advancing crop maturity, and increasing yield [50, 51]. Under drought stress, the expression of GA2ox, a key gene responsible for GA degradation, is upregulated [52]. Concurrently, GA modulates the interaction between the DELLA protein and the PIF3 and PIF4 genes, influencing hypocotyl cell elongation. Excessive GA results in the degradation of DELLA protein and inhibition of PIF4 activity [53, 54]. In this study, under PEG stress, the mRNA abundance of GA2ox increased while GA20ox expression decreased under osmotic stress, leading to a reduction in the biosynthesis of biologically active gibberellin and ultimately inhibiting hypocotyl elongation in alfalfa seedlings, which contradicted the observations in Fig. 1C. Furthermore, the expression of PIF4 was significantly upregulated, indicating a direct interaction between GA and PIF4. Upon the addition of melatonin, GA2ox expression was notably downregulated, promoting GA accumulation. This explains why the root length in the melatonin-treated group exceeded that observed under drought stress, highlighting its significance in drought response.
Melatonin is a significant class of indole compounds prevalent in plants. It serves multiple physiological roles, such as protecting chlorophyll, removing excess reactive oxygen species, and bolstering stress resistance in plants [55,56,57]. Numerous studies suggest that the increase in melatonin content in plant tissues following the application of exogenous melatonin may be due to the enhancement of endogenous melatonin synthesis. This involves enzymes like tyrosine decarboxylase (TDC) and acetylserotonin O-methyltransferase (ASMT), which promote the accumulation of endogenous melatonin [58,59,60]. In our experiment, transcriptomic analysis revealed that the key genes responsible for synthesizing ASMT (MsG0480020076.01) were significantly upregulated under drought stress and further intensified with the addition of melatonin. This aligns with findings from studies on other plant species, such as rice, where exogenous melatonin application improved salt tolerance by enhancing the antioxidant defense system [61]. Additionally, in wheat, similar treatments were shown to bolster drought resistance by increasing leaf water content [62].Research indicates that exogenous melatonin application can mitigate plant damage under adverse conditions and improve their tolerance to various abiotic stresses [63, 64]. Our findings showed that alfalfa seedlings treated with exogenous melatonin exhibited significantly reduced wilting and increased fresh weight under drought conditions compared to those not treated with melatonin. This effect is likely due to melatonin enhancing water absorption and organic matter accumulation [65, 66], which is crucial for alfalfa's ability to withstand drought stress. Previous studies indicate that melatonin plays a role in the regulatory network of plant hormone synthesis signals [67]. Furthermore, it can increase the level of GA in plants and significantly boost the expression of ABF genes during drought stress [68, 69]. Therefore, it can be inferred that GA, ABA, and melatonin are all involved in alfalfa's defense response to drought stress, with melatonin potentially being a key hormone in this regulatory network.
Interactions between flavonoids, plant hormones, and ABC transporters
Plant hormones and growth regulators play a crucial role in modulating several physiological processes, including development, metabolism, and aging, by interacting with nucleic acids, proteins, and enzymes in plants [70]. These elements also influence the synthesis of secondary metabolites such as flavonoids and terpenoids. Research has demonstrated that plant hormones can regulate the synthesis of flavonoid metabolites like quercetin, anthocyanins, and glycyrrhizin, key components in flavonoid biosynthesis [68]. In this study, transcriptomic analysis disclosed that several DEGs, primarily involved in flavonoid biosynthesis (such as CHS, CYP73A, and HCT), were induced under osmotic stress, with the transcriptional level of HCT notably enhanced with additional melatonin (MT). Furthermore, ATP-binding cassette (ABC) transporters are known to transport entities including plant hormones, terpenes, and alkaloids [71]. This research, under drought conditions, revealed through transcriptome analysis that the levels of gibberellins (GAs, MsG0180005519.01), MT, and certain flavonoids increased in alfalfa roots. The enrichment of ABC transporter genes (ABCB and ABCC) suggests that melatonin may contribute to regulating the transport of stress-related metabolites, including flavonoids and abscisic acid (ABA), across cellular membranes. This regulatory mechanism potentially strengthens root-to-shoot signaling and reinforces systemic antioxidant defenses, forming a coordinated strategy for adapting to drought conditions., particularly involving the ABCB and ABCC transporter proteins. Additionally, the genes responsible for anthocyanin synthesis within the flavonoid synthesis pathway showed significant upregulation. The presence of ABF2, PP2C, and PYL genes in the ABA signaling pathway also saw substantial enrichment in DEGs. ABA regulates stress responses by modulating stomatal closure and scavenging ROS. Its interaction with flavonoid biosynthesis genes likely enhances antioxidant defenses under drought stress. We propose that anthocyanin synthesis is closely linked to the ABF2, PP2C, and PYL genes [72]. In this process, ABA is detected by the PYL receptor, which inhibits the protein phosphatase PP2C and activates the kinase SnRK2. This activation triggers the phosphorylation of the transcription factor ABF2, facilitating its translocation to the nucleus, where it binds to ABA-responsive elements in the promoters of target genes—including those involved in flavonoid biosynthesis, such as HCT and ANR. Flavonoids, such as epicatechin, are crucial for neutralizing ROS during drought stress, while ABA minimizes water loss by reducing stomatal apertures. The combined action of flavonoid-mediated ROS scavenging and ABA-driven water conservation significantly boosts drought tolerance in alfalfa. Moreover, these findings indicate that ABC transporters may contribute to the movement of plant hormones, further regulating flavonoid synthesis and underpinning a complex, interconnected stress response system (Fig. 7).
Conclusions
Our findings demonstrate that melatonin (MT) significantly enhances the drought tolerance of alfalfa by harmonizing physiological responses and transcriptomic adjustments. From a physiological standpoint, the treatment group exposed to exogenous MT exhibited superior performance in root length and fresh weight compared to the control group, effectively mitigating the adverse effects of drought stress on alfalfa. This improvement is largely attributable to MT's ability to boost the activity of antioxidant enzymes, such as superoxide dismutase (SOD), thereby minimizing oxidative damage. At the molecular level, transcriptome sequencing combined with Weighted Gene Co-expression Network Analysis (WGCNA) revealed that MT regulates the expression of genes involved in flavonoid biosynthesis. Specifically, MT upregulated key genes, including PAL, CHS, and FLS, which facilitated the accumulation of flavonoids that play a critical role in scavenging reactive oxygen species (ROS). In addition, MT plays a pivotal role in the abscisic acid (ABA) and gibberellin (GA) hormone signaling pathways. MT enhances the expression of ABA-responsive genes, such as ABF2 and PYL4, while suppressing the activity of GA synthesis genes, like GA20ox. This dual regulation achieves an optimal balance between stomatal closure and growth recovery, ultimately contributing to alfalfa's improved drought tolerance. In conclusion, our study revealed a melatonin mediated regulatory network in alfalfa, which involves flavonoid biosynthesis, ABA/GA signaling pathway and ABC transporters. These findings provide key insights into the molecular mechanism of melatonin in alleviating drought stress in alfalfa, and provide operational target genes such as ABF2 and ASMT for improving crop drought resistance.
Materials and methods
Experimental materials and pretreatment
The Medicago sativa L. alfalfa cultivar was utilized in this study, provided by Professor Qingchuan Yang from the Institute of Animal Sciences, Chinese Academy of Agricultural Sciences. Alfalfa seeds were soaked in 75% ethanol for 10 min and then rinsed four times with double-distilled water. We germinated Thirty seeds in a Petri dish (with a diameter of 9 cm) containing 4 mL of melatonin-polyethylene glycol (PEG) solution at concentrations of 0, 10, 50, 100, 200, and 300 μM, based on previous research [73, 74]. The seeds were cultivated in a growth chamber at 25 °C with a relative humidity of 70% under a 16-h light and 8-h darkness cycle. The germination rate, fresh weight, and root length were measured for six biological replicates after seven days. For subsequent physiological and transcriptomic analyses, seedlings treated with 0 (designated as CK), 15% PEG (designated as ZMP), or 15% PEG + 10 μM melatonin (designated as ZMPMT) were collected, immediately frozen in liquid nitrogen, and stored at −80 °C.
Measurement of physiological indicators
The levels of GSH and O2− were determined using Solarbio kits (BC1175 and BC1295) from Beijing, China. SOD activity was assessed with the Abbkine kit (KTB1030) from Wuhan, China, whereas POD activity and MDA content were evaluated using Abbkine kits (KTB1150 and KTB1050), respectively. Each physiological index measurement was conducted in triplicate to ensure accuracy and reliability.
Transcriptomic analysis
Total RNA was extracted using the MJZol total RNA extraction kit (Shanghai Meiji Biomedicine Technology Co., Ltd., Shanghai, China). Complementary DNA (cDNA) was synthesized from the short RNA fragments, utilizing random hexamer primers and reverse transcriptase. The library fragments were subsequently purified with the Biowest agarose kit (Biowest, Logroño, Spain) and the RNA purification kit (Shanghai Meizhao Biomedicine Technology Co., Ltd.). Bridge-PCR clusters were generated on cBot. High-throughput sequencing was conducted on the Illumina NovaSeq 6000 platform (Illumina, USA). Upon completion of sequencing, the FASTp tool filtered the sequence quality. This involved removing the sequencing link, low-quality reads (trimming bases at the sequence's 3'end with quality scores below 20). If remaining sequences contained bases with a quality score below 10, they were discarded entirely. Reads with a high N-content ratio (exceeding 10%, where N represents uncertain bases) and sequences shorter than 20 base pairs, post adapter removal and quality pruning, were also eliminated. Mapping these data (reads) to a reference genome (Alfalfa Zhongmu No.1) enabled subsequent transcript assembly and expression calculations [75]. Differential gene expression was functionally annotated and classified using the Blast2GO program (https://www.blast2go.com/) [76, 77]. The KOBAS 3.0 online tool (http://kobas.cbi.pku.edu.cn/) identified drought stress response genes, assessing their enrichment in the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways (http://www.genome.jp/kegg/.) [78,79,80]. The raw sequencing reads underwent quality assessment using FastQC, followed by genome alignment carried out with HISAT2 (v2.2.1). Gene concentration was enriched and analyzed through software, with Fisher's exact test ensuring accuracy. The Benjamini–Hochberg method adjusted the p-values, and functions were deemed significantly enriched at an adjusted p-value (P-adjust) of less than 0.05. Differential expression analysis was conducted using DESeq2 R4.1.2(v1.30.1), a normalization approach grounded in the negative binomial distribution. The significance criteria were defined as |log2 FC|> 1 and FDR < 0.05.
Weighted gene co-correlation network analysis and protein interaction network analysis
The co-expression network was constructed using Weighted Gene Co-expression Network Analysis (WGCNA), identifying modules of closely related genes. Among the frequently used expression gene modules, those associated with phenotypic information are highlighted to explore gene networks and the correlation between phenotype and core network genes. The relevant functions of WGCNA are provided by the WGCNA package (https://cran.r-project.org/web/packages/WGCNA/index.html).
Protein interaction network analysis utilized the STRING (v11.5; https://string-db.org) database to construct a network of differential genes, illustrating the relationships among these differential genes [81].
qRT-PCR analysis
Eight differentially expressed genes were chosen for qRT-PCR analysis. Initially, full-length complementary DNA (cDNA) was synthesized using the PrimeScript RT kit (Cat# RR047A, TaKaRa, Tokyo, Japan), incorporating gDNA Eraser. The one-step qRT-PCR kit (Cat# RR420A; TaKaRa) was utilized following the manufacturer's guidelines. To ensure reliability, three independent RNA preparations served as biological replicates. The housekeeping gene MsACTIN2 was employed as an internal control to normalize gene expression. The primers were designed using Primer-BLAST (NCBI) and Oligo 7 software. The design parameters included a primer length of 18–24 bp, GC content of 40–60%, an annealing temperature (Tm) of 58–62 °C, and a product length of 80–200 bp. Primer specificity was thoroughly validated through melting curve analysis and sequencing. The detailed primer sequences are presented in Table 1.
Data processing
The data were organized and computed using Microsoft Excel 2021, while statistical analyses were conducted with SPSS 26.0. Data visualization was performed with GraphPad Prism 8.
Data availability
The data in this study can be requested from the corresponding author. The raw RNA-seq data has been deposited in the NCBI SRA database under the accession number PRJNA1113790.
References
Anjum SA, Ashraf U, Tanveer M, Khan I, Hussain S, Shahzad B, et al. Drought induced changes in growth, osmolyte accumulation and antioxidant metabolism of three maize hybrids. Front Plant Sci. 2017;8:69. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fpls.2017.00069.
Zhu Z, Liang Z, Han R. Saikosaponin accumulation and antioxidative protection in drought-stressed Bupleurum chinense DC. plants. Environ Exp Bot. 2009;66:326–33. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.envexpbot.2009.03.017.
Grzesiak MT, Marcińska I, Janowiak F, Hura T. Erratum to: the relationship between seedling growth and grain yield under drought conditions in maize and triticale genotypes. Acta Physiol Plant. 2012. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s11738-012-1039-2.
Jaleel CA, Gopi R, Sankar B, Gomathinayagam M, Panneerselvam R. Differential responses in water use efficiency in two varieties of Catharanthus roseus under drought stress. C R Biol. 2008;331:42–7. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.crvi.2007.11.003.
Jaleel CA, Manivannan P, Wahid A, Farooq M, Al-Juburi HJ, et al. Drought stress in plants: a review on morphological characteristics and pigments composition. Int J Agric Biol. 2009;11(1):100–5. https://doiorg.publicaciones.saludcastillayleon.es/10.3763/ijas.2009.0459.
Liu M, Li M, Liu K, Na S. Effects of drought stress on seed germination and seedling growth of different maize varieties. J Agric Sci. 2015;7(5). https://doiorg.publicaciones.saludcastillayleon.es/10.5539/jas.v7n5p231.
Liu Y, Wu Q, Ge G, et al. Influence of drought stress on afalfa yields and nutritional composition. BMC Plant Biol. 2018;18:13. https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12870-017-1226-9.
Zhao Y, Xiong L, Yin J, Zha X, Li W, Han Y. Understanding the effects of flash drought on vegetation photosynthesis and potential drivers over China. Sci Total Environ. 2024;931:172926. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.scitotenv.2024.172926.
Liu J, Shi K, Wang S, Zhu J, Wang X, Hong J, et al. MsCYP71 is a positive regulator for drought resistance in alfalfa. Plant Physiol Biochem. 2023;203:107999. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.plaphy.2023.107999.
Yin H, Wang Z, Li H, et al. MsTHI1 overexpression improves drought tolerance in transgenic alfalfa (Medicago sativa L.). Front Plant Sci. 2022;13:992024. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fpls.2022.992024. Published 2022 Sep 8.
Ku YS, Sintaha M, Cheung MY, Lam HM. Plant hormone signaling crosstalks between biotic and abiotic stress responses. Int J Mol Sci. 2018;19:3206. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/ijms19103206.
Arnao MB, Hernández-Ruiz J. Melatonin and its relationship to plant hormones. Ann Bot. 2018;121:195–207. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/aob/mcx114.
Simlat M, Ptak A, Skrzypek E, Warchoł M, Morańska E, Piórkowska E. Melatonin significantly influences seed germination and seedling growth of Stevia rebaudiana Bertoni. PeerJ. 2018;6:e5009. https://doiorg.publicaciones.saludcastillayleon.es/10.7717/peerj.5009.
Posmyk MM, Bałabusta M, Wieczorek M, Sliwinska E, Janas KM. Melatonin applied to cucumber (Cucumis sativus L.) seeds improves germination during chilling stress. J Pineal Res. 2009;46:214–23. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/j.1600-079X.2008.00652.x.
Arnao MB, Hernández-Ruiz J. Protective effect of melatonin against chlorophyll degradation during the senescence of barley leaves. J Pineal Res. 2009;46:58–63. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/j.1600-079X.2008.00625.x.
Li G, Li Y, Zhu Y, Zheng W, Li M, Hu J, et al. Exogenous application of melatonin to mitigate drought stress-induced oxidative damage in Phoebe sheareri seedlings. PeerJ. 2023;11:e15159. https://doiorg.publicaciones.saludcastillayleon.es/10.7717/peerj.15159.
Tiwari RK, Lal MK, Kumar R, Chourasia KN, Naga KC, Kumar D, et al. Mechanistic insights on melatonin-mediated drought stress mitigation in plants. Physiol Plant. 2021;172:1212–26. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/ppl.13307.
Altaf MA, Shahid R, Ren MX, Naz S, Altaf MM, Khan LU, et al. Melatonin improves drought stress tolerance of tomato by modulating plant growth, root architecture, photosynthesis, and antioxidant defense system. Antioxid Basel Switz. 2022;11:309. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/antiox11020309.
Ahmad S, Cui W, Kamran M, et al. Exogenous application of melatonin induces tolerance to salt stress by improving the photosynthetic efficiency and antioxidant defense system of maize seedling. J Plant Growth Regul. 2021;40:1270–83. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s00344-020-10187-0.
Xie Z, Jin L, Sun Y, et al. OsNAC120 balances plant growth and drought tolerance by integrating GA and ABA signaling in rice. Plant Commun. 2024;5(3):100782. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.xplc.2023.100782.
Wang K, Nan L-L, Xia J, Wu S-W, Yang L-L. Metabolomics reveal root differential metabolites of different root-type alfalfa under drought stress. Front Plant Sci. 2024;15:1341826. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fpls.2024.1341826.
Conti V, Parrotta L, Romi M, Del Duca S, Cai G. Tomato biodiversity and drought tolerance: a multilevel review. Int J Mol Sci. 2023;24:10044. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/ijms241210044.
Meng L-S. Compound synthesis or growth and development of roots/stomata regulate plant drought tolerance or water use efficiency/water uptake efficiency. J Agric Food Chem. 2018;66:3595–604. https://doiorg.publicaciones.saludcastillayleon.es/10.1021/acs.jafc.7b05990.
Miller G, Suzuki N, Ciftci-Yilmaz S, Mittler R. Reactive oxygen species homeostasis and signalling during drought and salinity stresses. Plant Cell Environ. 2010;33:453–67. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/j.1365-3040.2009.02041.x.
Qi J, Song C-P, Wang B, Zhou J, Kangasjärvi J, Zhu J-K, et al. Reactive oxygen species signaling and stomatal movement in plant responses to drought stress and pathogen attack. J Integr Plant Biol. 2018;60:805–26. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/jipb.12654.
Mittler R. Oxidative stress, antioxidants and stress tolerance. Trends Plant Sci. 2002;7:405–10. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/s1360-1385(02)02312-9.
Kumar A, Majeti NVP. Proteomic responses to lead-induced oxidative stress in Talinum triangulare Jacq. (Willd.) roots: identification of key biomarkers related to glutathione metabolisms. Environ Sci Pollut Res Int. 2014;21:8750–64. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s11356-014-2808-9.
Apel K, Hirt H. Reactive oxygen species: metabolism, oxidative stress, and signal transduction. Annu Rev Plant Biol. 2004;55:373–99. https://doiorg.publicaciones.saludcastillayleon.es/10.1146/annurev.arplant.55.031903.141701.
Yun Z, Gao H, Chen X, Chen Z, Zhang Z, Li T, et al. Effects of hydrogen water treatment on antioxidant system of litchi fruit during the pericarp browning. Food Chem. 2021;336:127618. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.foodchem.2020.127618.
Wang L, Chen S, Shao J, Zhang C, Mei L, Wang K, et al. Hydrogen sulfide alleviates chilling injury in peach fruit by maintaining cell structure integrity via regulating endogenous H2S, antioxidant and cell wall metabolisms. Food Chem. 2022;391:133283. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.foodchem.2022.133283.
Gill SS, Tuteja N. Reactive oxygen species and antioxidant machinery in abiotic stress tolerance in crop plants. Plant Physiol Biochem PPB. 2010;48:909–30. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.plaphy.2010.08.016.
Ji W, Zhu Y, Li Y, Yang L, Zhao X, Cai H, et al. Over-expression of a glutathione S-transferase gene, GsGST, from wild soybean (Glycine soja) enhances drought and salt tolerance in transgenic tobacco. Biotechnol Lett. 2010;32:1173–9. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s10529-010-0269-x.
Huang C, Guo T, Zheng SC, Feng QL, Liang JH, et al. Increased cold tolerance in Arabidopsis thaliana transformed with Choristoneura fumiferana glutathione S-transferase gene. Biol Planta. 2009;53:183–7.
Gordillo GM, Biswas A, Khanna S, Spieldenner JM, Pan X, Sen CK. Multidrug resistance-associated protein-1 (MRP-1)-dependent Glutathione Disulfide (GSSG) efflux as a critical survival factor for oxidant-enriched tumorigenic endothelial cells. J Biol Chem. 2016;291:10089–103. https://doiorg.publicaciones.saludcastillayleon.es/10.1074/jbc.M115.688879.
Irfan M, Ahmad A, Hayat S. Effect of cadmium on the growth and antioxidant enzymes in two varieties of Brassica juncea. Saudi J Biol Sci. 2014;21:125–31. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.sjbs.2013.08.001.
Hayat S, Hayat Q, Alyemeni MN, Wani AS, Pichtel J, Ahmad A. Role of proline under changing environments: a review. Plant Signal Behav. 2012;7:1456–66. https://doiorg.publicaciones.saludcastillayleon.es/10.4161/psb.21949.
Yang PM, Huang QC, Qin GY, Zhao SP, Zhou JY. Different drought-stress responses in photosynthesis and reactive oxygen metabolism between autotetraploid and diploid rice. Photosynthetica. 2014;52(2):193–202. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s11099-014-0020-2.
Zhang Z, Liu L, Li H, Zhang S, Fu X, Zhai X, et al. Exogenous melatonin promotes the salt tolerance by removing active oxygen and maintaining ion balance in wheat (Triticum aestivum L.). Front Plant Sci. 2021;12:787062. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fpls.2021.787062.
Treutter D. Significance of flavonoids in plant resistance and enhancement of their biosynthesis. Plant Biol Stuttg Ger. 2005;7:581–91. https://doiorg.publicaciones.saludcastillayleon.es/10.1055/s-2005-873009.
Ammar RB, Bhouri W, Sghaier MB, Boubaker J, Skandrani I, Neffati A, et al. Antioxidant and free radical-scavenging properties of three flavonoids isolated from the leaves of Rhamnus alaternus L. (Rhamnaceae): a structure-activity relationship study. Food Chem. 2009;116:258–64. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.foodchem.2009.02.043.
Sun M, Li L, Wang C, Wang L, Lu D, Shen D, et al. Naringenin confers defence against Phytophthora nicotianae through antimicrobial activity and induction of pathogen resistance in tobacco. Mol Plant Pathol. 2022;23(12):1737–50. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/mpp.13255.
Luo P, Shen Y, Jin S, Huang S, Cheng X, Wang Z, et al. Overexpression of Rosa rugosa anthocyanidin reductase enhances tobacco tolerance to abiotic stress through increased ROS scavenging and modulation of ABA signaling. Plant Sci Int J Exp Plant Biol. 2016;245:35–49. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.plantsci.2016.01.007.
Verma V, Ravindran P, Kumar PP. Plant hormone-mediated regulation of stress responses. BMC Plant Biol. 2016;16:86. https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12870-016-0771-y.
Sabagh AE, Mbarki S, Hossain A, Iqbal MA, Islam MS, Raza A, et al. Potential role of plant growth regulators in administering crucial processes against abiotic stresses. Front Agron. 2021;3:648694. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fagro.2021.648694.
Mukherjee A, Dwivedi S, Bhagavatula L, Datta S. Integration of light and ABA signaling pathways to combat drought stress in plants. Plant Cell Rep. 2023;42:829–41. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s00299-023-02999-7.
Marusig D, Tombesi S. Abscisic acid mediates drought and salt stress responses in Vitis vinifera: a review. Int J Mol Sci. 2020;21:8648. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/ijms21228648.
Li C, Tan DX, Liang D, Chang C, Jia D, Ma F. Melatonin mediates the regulation of ABA metabolism, free-radical scavenging, and stomatal behaviour in two Malus species under drought stress. J Exp Bot. 2015;66(3):669–80. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/jxb/eru476.
Lee SU, Mun BG, Bae EK, Kim JY, Kim HH, Shahid M, et al. Drought stress-mediated transcriptome profile reveals NCED as a key player modulating drought tolerance in Populus davidiana. Front Plant Sci. 2021;12:755539. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fpls.2021.755539.
Fujita Y, Fujita M, Shinozaki K, Yamaguchi-Shinozaki K. ABA-mediated transcriptional regulation in response to osmotic stress in plants. J Plant Res. 2011;124:509–25. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s10265-011-0412-3.
Park EJ, Kim HT, Choi YI, Lee C, Nguyen VP, Jeon HW, et al. Overexpression of gibberellin 20-oxidase1 from Pinus densiflora results in enhanced wood formation with gelatinous fiber development in a transgenic hybrid poplar. Tree Physiol. 2015;35:1264–77. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/treephys/tpv099.
Sakamoto T, Miura K, Itoh H, Tatsumi T, Ueguchi-Tanaka M, Ishiyama K, et al. An overview of gibberellin metabolism enzyme genes and their related mutants in rice. Plant Physiol. 2004;134(4):1642–53. https://doiorg.publicaciones.saludcastillayleon.es/10.1104/pp.103.033696.
Band LR, Nelissen H, Preston SP, Rymen B, Prinsen E, AbdElgawad H, et al. Modeling reveals posttranscriptional regulation of GA metabolism enzymes in response to drought and cold. Proc Natl Acad Sci U S A. 2022;119:e2121288119. https://doiorg.publicaciones.saludcastillayleon.es/10.1073/pnas.2121288119.
Li K, Yu R, Fan LM, Wei N, Chen H, Deng XW. DELLA-mediated PIF degradation contributes to coordination of light and gibberellin signalling in Arabidopsis. Nat Commun. 2016;7:11868. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/ncomms11868.
Wang F, Chen X, Dong S, Jiang X, Wang L, Yu J, Zhou Y. Crosstalk of PIF4 and DELLA modulates CBF transcript and hormone homeostasis in cold response in tomato. Plant Biotechnol J. 2020;18(4):1041–55. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/pbi.13272.
Zhang J, Shi Y, Zhang X, Du H, Xu B, Huang B. Melatonin suppression of heat-induced leaf senescence involves changes in abscisic acid and cytokinin biosynthesis and signaling pathways in perennial ryegrass (Lolium perenne L.). Environ Exp Bot. 2017;138:36–45. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.envexpbot.2017.02.012.
Arnao MB, Hernández-Ruiz J. Melatonin in flowering, fruit set and fruit ripening. Plant Reprod. 2020;33(2):77–87. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s00497-020-00388-8.
Onik JC, Wai SC, Li A, Lin Q, Sun Q, Wang Z, et al. Melatonin treatment reduces ethylene production and maintains fruit quality in apple during postharvest storage. Food Chem. 2021;337:127753. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.foodchem.2020.127753.
Xu F, Liu W, Wang H, Alam P, Zheng W, Faizan M. Genome identification of the tea plant (Camellia sinensis) ASMT gene family and its expression analysis under abiotic stress. Genes. 2023;14:409. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/genes14020409.
Shamloo-Dashtpagerdi R, Lindlöf A, Nouripour-Sisakht J. Unraveling the regulatory role of MYC2 on ASMT gene expression in wheat: implications for melatonin biosynthesis and drought tolerance. Physiol Plant. 2023;175:e14015. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/ppl.14015.
Xie Z, Wang J, Wang W, et al. Integrated analysis of the transcriptome and metabolome revealed the molecular mechanisms underlying the enhanced salt tolerance of rice due to the application of exogenous melatonin. Front Plant Sci. 2021;11:618680. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fpls.2020.618680.
Cui G, Zhao X, Liu S, Sun F, Zhang C, Xi Y. Beneficial effects of melatonin in overcoming drought stress in wheat seedlings. Plant Physiol Biochem. 2017;118:138–49. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.plaphy.2017.06.014.
Yang WJ, Du YT, Zhou YB, Chen J, Xu ZS, Ma YZ, et al. Overexpression of TaCOMT improves melatonin production and enhances drought tolerance in transgenic Arabidopsis. Int J Mol Sci. 2019;20:652. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/ijms20030652.
Zhao D, Wang H, Chen S, Yu D, Reiter RJ. Phytomelatonin: an emerging regulator of plant biotic stress resistance. Trends Plant Sci. 2021;26:70–82. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.tplants.2020.08.009.
Sadak MS, Abdalla AM, Elhamid EMA, Ezzo M. Role of melatonin in improving growth, yield quantity and quality of Moringa oleifera L. plant under drought stress. Bull Natl Res Cent. 2020;44(1). https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s42269-020-0275-7.
Bai Y, Xiao S, Zhang Z, Zhang Y, Sun H, Zhang K, et al. Melatonin improves the germination rate of cotton seeds under drought stress by opening pores in the seed coat. PeerJ. 2020;8:e9450. https://doiorg.publicaciones.saludcastillayleon.es/10.7717/peerj.9450.
Hu W, Yang H, Tie W, Yan Y, Ding Z, Liu Y, et al. Natural variation in banana varieties highlights the role of melatonin in postharvest ripening and quality. J Agric Food Chem. 2017;65:9987–94. https://doiorg.publicaciones.saludcastillayleon.es/10.1021/acs.jafc.7b03354.
Khan M, Ali S, Manghwar H, Saqib S, Ullah F, Ayaz A, et al. Melatonin function and crosstalk with other phytohormones under normal and stressful conditions. Genes (Basel). 2022;13(10):1699. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/genes13101699.
Li J, Xie J, Yu J, Lyv J, Zhang J, Ding D, et al. Melatonin enhanced low-temperature combined with low-light tolerance of pepper (Capsicum annuum L.) seedlings by regulating root growth, antioxidant defense system, and osmotic adjustment. Front Plant Sci. 2022;13:998293. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fpls.2022.998293.
Banerjee A, Roychoudhury A. Melatonin application reduces fluoride uptake and toxicity in rice seedlings by altering abscisic acid, gibberellin, auxin and antioxidant homeostasis. Plant Physiol Biochem PPB. 2019;145:164–73. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.plaphy.2019.10.033.
Bari R, Jones JD. Role of plant hormones in plant defence responses. Plant Mol Biol. 2009;69(4):473–88. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s11103-008-9435-0.
Kang J, Hwang JU, Lee M, Kim YY, Assmann SM, Martinoia E, et al. PDR-type ABC transporter mediates cellular uptake of the phytohormone abscisic acid. Proc Natl Acad Sci U S A. 2010;107(5):2355–60. https://doiorg.publicaciones.saludcastillayleon.es/10.1073/pnas.0909222107.
Wang X, Yin J, Wang J, Li J. Integrative analysis of transcriptome and metabolome revealed the mechanisms by which flavonoids and phytohormones regulated the adaptation of alfalfa roots to NaCl stress. Front Plant Sci. 2023;14:1117868. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fpls.2023.1117868.
Roy M, Niu J, Irshad A, Kareem HA, Hassan MU, et al. Exogenous melatonin protects alfalfa (Medicago sativa L.) seedlings from drought-induced damage by modulating reactive oxygen species metabolism, mineral balance and photosynthetic efficiency. Plant Stress. 2021;2:100044. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.stress.2021.100044.
Aguilera Y, Herrera T, Liébana R, Rebollo-Hernanz M, Sanchez-Puelles C, MartÃn-Cabrejas MA. Impact of melatonin enrichment during germination of legumes on bioactive compounds and antioxidant activity. J Agric Food Chem. 2015;63(36):7967–74. https://doiorg.publicaciones.saludcastillayleon.es/10.1021/acs.jafc.5b03128.
Shen C, Du H, Chen Z, Lu H, Zhu F, Chen H, et al. The chromosome-level genome sequence of the autotetraploid alfalfa and resequencing of core germplasms provide genomic resources for alfalfa research. Mol Plant. 2020;13:1250–61. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.molp.2020.07.003.
Gene Ontology Consortium, Aleksander SA, Balhoff J, Carbon S, Cherry JM, Drabkin HJ, et al. The gene ontology knowledgebase in 2023. Genetics. 2023;224:iyad031. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/genetics/iyad031.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene ontology: tool for the unification of biology. the Gene Ontology Consortium. Nat Genet. 2000;25:25–9. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/75556.
Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28:27–30. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/nar/28.1.27.
Kanehisa M. Toward understanding the origin and evolution of cellular organisms. Protein Sci Publ Protein Soc. 2019;28:1947–51. https://doiorg.publicaciones.saludcastillayleon.es/10.1002/pro.3715.
Kanehisa M, Furumichi M, Sato Y, Kawashima M, Ishiguro-Watanabe M. KEGG for taxonomy-based analysis of pathways and genomes. Nucleic Acids Res. 2023;51:D587–92. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/nar/gkac963.
Szklarczyk D, Gable AL, Lyon D, et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47(D1):D607–13. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/nar/gky1131.
Acknowledgements
We thank Shanghai Majorbio Biotechnology Co., Ltd. for their help with sequencing.
Funding
This work was supported by the earmarked fund for the Biological Breeding-National Science and Technology Major Project (2022ZD04011), the National Key R&D Program of China (2024YFF1001300), the China Postdoctoral Science Foundation (2022M721043), the China Agriculture Research System (CARS-34). The funding body played no role in the design of the study, the collection, analysis, and interpretation of the data, or the writing of the manuscript.
Author information
Authors and Affiliations
Contributions
XLR, XYZ, YHS, CZW, DFL, HS: Conceptualization and methodology; XXZ, WXZ, ZRL, GML, YAL: Resources; XXZ, WXZ, ZRL, GML, JZW, YAL: Investigation; XXZ, WXZ, ZRL: Data curation, Formal analysis, and writing – original draft; HS: Supervision and funding acquisition. All authors have read and agreed to the published version of the manuscript.
Corresponding authors
Ethics declarations
Consent for publication
We declare that all experimental consumables and plants are collected with permission from the local agricultural department, and all plant experimental research complies with the requirements of relevant institutions, national and international standards, and legislation.
Competing interests
The authors declare no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
About this article
Cite this article
Zhang, X., Zhu, W., Liu, Z. et al. Transcriptomic analysis of melatonin-mediated drought stress response genes in alfalfa during germination period. BMC Plant Biol 25, 637 (2025). https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12870-025-06665-w
Received:
Accepted:
Published:
DOI: https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12870-025-06665-w