- Research
- Open access
- Published:
Dynamic physiological and transcriptomic changes reveal memory effects of salt stress in maize
BMC Genomics volume 24, Article number: 726 (2023)
Abstract
Background
Pre-exposing plants to abiotic stresses can induce stress memory, which is crucial for adapting to subsequent stress exposure. Although numerous genes involved in salt stress response have been identified, the understanding of memory responses to salt stress remains limited.
Results
In this study, we conducted physiological and transcriptional assays on maize plants subjected to recurrent salt stress to characterize salt stress memory. During the second exposure to salt stress, the plants exhibited enhanced salt resistance, as evidenced by increased proline content and higher POD and SOD activity, along with decreased MDA content, indicative of physiological memory behavior. Transcriptional analysis revealed fewer differentially expressed genes and variations in response processes during the second exposure compared to the first, indicative of transcriptional memory behavior. A total of 2,213 salt stress memory genes (SMGs) were identified and categorized into four response patterns. The most prominent group of SMGs consisted of genes with elevated expression during the first exposure to salt stress but reduced expression after recurrent exposure to salt stress, or vice versa ([+ / −] or [− / +]), indicating that a revised response is a crucial process in plant stress memory. Furthermore, nine transcription factors (TFs) (WRKY40, WRKY46, WRKY53, WRKY18, WRKY33, WRKY70, MYB15, KNAT7, and WRKY54) were identified as crucial factors related to salt stress memory. These TFs regulate over 53% of SMGs, underscoring their potential significance in salt stress memory.
Conclusions
Our study demonstrates that maize can develop salt stress memory, and the genes identified here will aid in the genetic improvement of maize and other crops.
Background
Soil salinization, typically caused by the accumulation of NaCl, is widespread abiotic stress that has deleterious impacts on plant growth and distribution. The adverse effects of salt stress can be attributed to two primary causes as follows [1,2,3]: On one hand, there is osmotic stress, which results in reduced root water absorption, stomatal closure, and inhibition of cell expansion. On the other hand, there is ionic toxicity, which causes oxidative damage and excessive generation of reactive oxygen species (ROS). Excessive ROS can damage membranes, disturb cellular metabolism, and accelerate plant senescence and death. Therefore, comprehending the underlying mechanisms governing plant responses to salt stress and the formulation of efficacious strategies for enhancing salt tolerance emerge as pivotal approaches to mitigating the loss of arable land due to soil salinization.
Plants have evolved sophisticated mechanisms to respond and adapt to abiotic stresses, which involve complex networks of numerous genes. A pivotal aspect of stress response is the regulation of functional genes by transcription factors (TFs), including various members of the NAC, WRKY, MYB, and bHLH families [2, 4,5,6], which activate protective mechanisms enabling plants to cope with abiotic stress. Notably, recent studies have demonstrated that plants can memorize their previous exposure to abiotic stress at physiological and transcriptional levels [7,8,9], leading to enhanced tolerance or resistance upon subsequent stress encounters, which is called stress memory [10]. For example, dehydration stress memory has been observed in several plant species, including Arabidopsis [11, 12], maize [13, 14], and switchgrass [7]. Plants that have experienced dehydration/recovery cycles exhibit improved performance, including an increased ability to maintain leaf relative water content when subsequently exposed to dehydration stress. By using transcriptome deep sequencing (RNA-seq), genes associated with dehydration stress memory have been identified [7, 13]. These genes exhibited significantly divergent transcript responses between the initial and subsequent stress events, and they can be classified into four categories as follows: (1) genes that are up-regulated during the initial stress event and exhibit elevated expression during the subsequent stress event (designated as [+ / +]), (2) genes that are down-regulated during the initial stress event and display reduced expression during the subsequent stress event ([-/-]), (3) genes that are up-regulated during the initial stress event but experience reduced expression during the subsequent stress event ([+/-]) and (4) genes that exhibit the opposite pattern ([-/ +]). Stress memory results in fewer changes in gene expression during subsequent exposure to the same type of stress [15, 16], suggesting a reduction in the disturbance of gene expression. Interestingly, studies comparing Arabidopsis, maize, and switchgrass have revealed the existence of both conserved and species-specific mechanisms that regulate plant responses to recurring dehydration stress [7], highlighting the intricate regulatory landscape of stress memory across plant species. Meanwhile, there is emerging evidence to suggest that epigenetic mechanisms, including DNA methylation and histone modifications such as glycosylation, acetylation, ubiquitination, and methylation alterations, serve as pivotal contributors to the establishment of environmental stress memory [17,18,19,20].
Maize (Zea mays), a globally significant cereal crop, plays a pivotal role in providing food, feed and bioenergy. However, its susceptibility to salt stress, particularly during the seedling stage, imposes considerable constraints on its growth, development, yield, and quality [21, 22]. To mitigate this challenge, it becomes imperious to: (1) elucidate the molecular mechanisms underlying salt tolerance in maize, (2) identify novel stress-resistant genes, and (3) breed new cultivars with enhanced resilience to salt stress. Stress memory enables plants to mount vigorous and rapid responses upon subsequent encounters with such stresses, consequently facilitating their recuperation. Therefore, the salt stress resistance and productivity of maize can be improved by inducing stress memory. However, our current understanding of salt stress memory and the underlying molecular mechanisms in maize remains limited.
In this study, we investigated the phenomenon of salt stress memory in maize using a precisely controlled liquid culture-based experimental framework. Maize plants were subjected to recurrent salt stress (200 mM NaCl for 0 h, 3 h, and 48 h), and their phenotypic, physiological, and transcriptomic responses to salt stress were evaluated. The results revealed that salt stress memory was established in maize. Protective components, such as proline, and processes, notably photosynthesis, were found to play an important role in the development of salt stress memory. Furthermore, we revealed the significance of TFs in initiating the salt stress memory in maize and identified nine pivotal TFs associated with stress memory. These results provide valuable insights into the mechanisms underlying salt stress memory in maize and may help to design genetic improvement strategies aimed at enhancing the salt stress resistance and productivity of maize.
Results
Phenotypic and physiological analyses indicated that maize displays salt stress memory
A salt stress experiment was conducted to investigate the presence of salt stress memory in maize plants (Fig. 1A). Briefly, maize seedlings were cultivated in pots for one week and subsequently transferred to a hydroponic culture using Hoagland’s liquid medium for five days. Subsequently, the plants were exposed to 200 mM NaCl for 3 h, and then divided into two distinct experimental groups: Group A and Group B. Plants in Group A were subjected to an additional 45 h exposure to 200 mM NaCl (a total of 48 h). Plants in Group B were transferred to Hoagland’s liquid medium for three days of recovery culture before the second exposure to salt stress. Plants subjected to salt stress for 0 h, 3 h, and 48 h during both the first and second stress events were photographed and observed phenotypically (Fig. 1B). After 48 h of salt stress, it was evident that plants in Group A exhibited a more severe stress phenotype compared to those in Group B, characterized by pronounced leaf wilting, suggesting the development of salt stress memory following the initial salt stress treatment.
To validate these findings, root samples were collected for the measurement of physiological parameters, including proline and malondialdehyde (MDA) content, and the activity of superoxide dismutase (SOD) and peroxidase (POD) (Fig. 1C). Compared with non-saline conditions, exposure to salt stress for 3 h and 48 h induced a significant increase in proline content and POD activity, while decreasing MDA content. Notably, these changes were more pronounced in plants subjected to pretreatment. Moreover, salt stress led to a substantial reduction in SOD activity in both un-pretreated and pretreated plants; however, SOD activity trended to be higher in pretreated plants than in un-pretreated plants. Overall, these results indicated that a 3 h salt pretreatment significantly enhances the salt tolerance ability of maize plants, affirming the presence of salt stress memory in salt-pretreated plants.
RNA-seq analysis revealed distinct transcriptional responses to repeated salt stress events
To examine the impact of recurrent exposure to salt stress at the transcriptomic level, RNA-seq was performed using root samples collected at 0 h, 3 h, and 48 h of salt treatment (designated as T1_0h, T1_3h, and T1_48h, respectively), and 3 h and 48 h of a subsequent exposure to salt stress (designated as T2_3h and T2_48h, respectively) (Fig. 1A). A total of 579.5 million paired-end reads were obtained, with an average of 38.6 million reads per sample (Table S1). These reads were mapped to the maize reference genome, with an average mapping rate of 72.9%, consistent with what was previously reported [23]. Pearson correlation coefficients and Principal component analysis (PCA) demonstrated a high similarity among biological replicates (Fig. 2A and B), suggesting the high quality and reproducibility of the RNA-seq data.
A total of 11,074 differentially expressed genes (DEGs) (fold change ≥ 2 and adjusted P-value ≤ 0.05) were identified in the T1_3h, T1_48h, T2_3h, and T2_48h groups when compared to the T1_0h control group (Fig. 2C, Table S2). Notably, the number of DEGs increased with the increasing duration of stress exposure during both the first and second stress events (T1_3h: 6,267 DEGs; T1_48h: 6,981 DEGs; T2_3h: 4,899 DEGs; T2_48h: 5,241 DEGs), indicating the activation or suppression of a greater number of genes as stress duration prolonged. Intriguingly, we observed that the number of DEGs in the T2_3h and T2_48h plants decreased by 21.8% and 24.9% when compared with the T1_3h and T1_48h plants, respectively. The greater number of DEGs identified at the first exposure to salt stress may be responsible for the more severe phenotypic changes (Fig. 1B). To further illustrate the functional roles of these DEGs, we performed Gene Ontology (GO) enrichment analysis for each of the four sets of DEGs (Fig. 2D, Table S3). This analysis revealed that they shared several enriched GO terms involved in the oxidation–reduction process, response to oxidative stress, trehalose biosynthetic process, and cellular ion homeostasis (Fig. 2D), which are critical processes for plant stress adaptation. However, there were also substantial differences in enriched GO terms between the first and second salt stress events. GO terms related to photosynthesis, water response, anion transport, and transcription regulation were predominantly enriched in the T1_3h samples, whereas GO terms related to ion transport, G-protein coupled receptor signalling pathway, and nitrate assimilation were enriched in the T2_3h samples (Fig. 2D). DEGs in the T1_48h samples enriched GO terms related to translation, amide biosynthetic and metabolic process, and organonitrogen compound biosynthetic and metabolic process, whereas DEGs in the T2_48h samples enriched GO terms related to drug response and transport, and nitrogen compound transport. Taken together, despite a partial overlap in stress response, there were distinct transcriptional responses to recurrent salt stress in maize.
Identification of salt stress memory genes in maize
To explore the mechanisms underlying salt memory in maize, we identified salt stress memory genes (SMGs). Stress memory genes refer to genes whose transcriptional responses significantly differ between the first and subsequent stress events, which can be categorized into four types [7, 11, 16]. Accordingly, SMGs were identified by conducting a comparative analysis of DEGs between T1_3h vs. T1_0h and T2_3h vs. T1_3h, as previously done in switchgrass and rice [7, 16]. Genes displaying an up-regulation in both T1_3h vs T1_0h and T2_3h vs T1_3h ([+ / +]), down-regulation in both comparisons ([-/-]), or up-regulation in T1_3h vs T1_0h but down-regulation in T2_3h vs T1_3h or vice versa ([+/-] or [-/ +]) were considered to be SMGs. A total of 2,213 SMGs were identified, consisting of 94 [+ / +] genes, 75 [-/-] genes, 1,392 [+/-] genes, and 652 [-/ +] genes (Fig. 3A, Table S4). Notably, the number of [+ / −] genes was higher than the other three types of genes, accounting for 62.9% of all SMGs.
Then, we performed GO enrichment analysis on these SMGs, unveiling their significant enrichment in 30 pathways (adjusted P-value ≤ 0.05) (Table S5). The representative enriched pathways are involved in oxidation–reduction process, photosynthesis, regulation of gene expression, and trehalose biosynthetic process (Fig. 3B). For example, the Zm00001d024852 (ACO5) gene, which plays a role in the oxidation of 1-aminocyclopropane-1-carboxylic acid (ACC) to ethylene [24], exhibited an up-regulation in both T1_3h vs T1_0h and T2_3h vs T1_3h ([+ / +]) (Fig. 3C). In contrast, the Zm00001d052227 (TPPD) gene, which converts trehalose-6-phosphate (T6P) to trehalose [25], displayed a down-regulation in both comparisons ([-/-]) (Fig. 3D). The gene related to the light reactions of photosynthesis, Zm00001d032197 (CAB4) [26, 27], was up-regulated in T1_3h vs T1_0h but down-regulated in T2_3h vs T1_3h ([+/-]) (Fig. 3E). However, the Zm00001d031164 (ANAC042) gene, which encodes a transcription factor [28], was down-regulated in T1_3h vs T1_0h but up-regulated in T2_3h vs T1_3h ([-/ +]) (Fig. 3F).
Transcription factors network analysis for identification of crucial regulators
Given the important role of TFs in the formation and maintenance of memory behavior in both animals and plants [29,30,31,32,33], we wonder whether TFs participated in salt stress memory in maize. By comparing sequence homologies with entries in the PlantTFDB database, 11.5% (255/2,213) of SMGs were identified as TF genes distributed across 38 families. Notably, the NAC (33 genes) and MYB (33 genes) families exhibited the highest abundance, closely followed by the WRKY (31 genes) and bHLH (29 genes) families (Table S6). Interestingly, the proportion of TF genes in SMGs was significantly higher than that in the maize genome (6.0% of all genes [2,361/39,498] were TF genes; P < 0.01, Fisher’s exact test) (Fig. 4A). Additionally, we assessed the enrichment of each TF family within SMGs and found that seven TF families, including NAC, WRKY, MYB, bHLH, HD-ZIP, G2-like, and TALE, exhibited marked enrichment (P < 0.01, Fisher’s exact test) (Fig. 4B). Altogether, these results indicate that specific TFs are involved in the regulation of salt stress memory in maize.
To further reveal potential relationships among the enriched TFs, we established a protein–protein interaction (PPI) network using the STRING database. As shown in Fig. 4C, the PPI network consisted of 92 nodes (TFs) and 219 edges (interactions). Subsequently, we identified hub TFs within the network using the Molecular Complex Detection (MCODE) and Maximal Clique Centrality (MCC) algorithms, yielding 18 and 10 TFs, respectively. Notably, nine TFs, namely WRKY40, WRKY46, WRKY53, WRKY18, WRKY33, WRKY70, MYB15, KNAT7, and WRKY54, were considered crucial hub TFs due to their identification by both methods (Fig. 4D-E). Interestingly, except for WRKY18 and WRKY46, the remaining seven TFs belong to [+ / −] memory genes. Several of these TFs have been demonstrated to play an important role in the regulation of salt stress responses. For example, WRKY46 regulates lateral root development under salt stress conditions by regulating ABA signaling and auxin homeostasis [34, 35]. Furthermore, WRKY33 [36, 37], WRKY40 [38], WRKY53 [39], and MYB15 [40] act as transcriptional activators enhancing salt tolerance, whereas WRKY18 heightens sensitivity to salt and osmotic stress in plants [41].
Validation of expression pattern of TF genes using qRT-PCR
To validate the results of RNA-seq, we conducted quantitative real-time PCR (qRT-PCR) experiments to quantify the expression levels of the nine crucial TF genes (Fig. 5, Figs. S1 and S2, Table S7). The measured expression profiles of these genes were consistent with the results obtained from the RNA-seq analysis, with the exception of Zm00001d005056. Zm00001d005056 was identified as a [+ / +] gene via RNA-seq but as a [+ / −] gene via qRT-PCR. This discrepancy may be attributed to technical or biological reasons [42]. Overall, the RNA-seq data demonstrated high reliability in assessing gene expression in response to salt stress in maize. The nine TFs identified could potentially contribute to the improvement in salt tolerance during recurrent exposure to salt stress.
Predicted crucial TFs directly regulated SMGs
We want to know whether the nine crucial TFs identified have the capacity to directly regulate SMGs. To this end, we scanned the promoter regions (defined as 1 kb upstream of the transcription start site [TSS] of a gene) of SMGs using the corresponding TF binding motifs to infer potential TF binding sites, and thus identified the target genes for each TF. Since only six of the nine TFs had accessible motifs in public databases, we concentrated our analysis on these six TFs, namely, WRKY40, WRKY46, WRKY18, WRKY33, WRKY70, and MYB15. We discovered 393 predicted target genes for WRKY40, 376 for WRKY46, 332 for WRKY18, 396 for WRKY33, 304 for WRKY70, and 701 for MYB15, respectively (Fig. 6A). Representative examples are shown in Fig. 6B. In total, 1,192 target genes were detected, accounting for more than half (53.9%) of all SMGs. Thus, these TFs tend to directly regulate the expression of SMGs. Interestingly, our analysis revealed that out of the 1,192 target genes, 36 were predicted to be the common targets for all six TFs (Fig. 6A), consisting of one [+ / +] gene, two [-/-] genes, 21 [+/-] genes, and 12 [-/ +] genes. We performed GO analysis on this group of genes, revealing their distribution across diverse categories (Fig. 6C). Within these categories, cellular process (27.8%, 10/36) and metabolic process (27.8%, 10/36), membrane (13.9%, 5/36), and binding (44.4%, 16/36) were prominently represented in biological process, cellular component and molecular function, respectively.
Discussion
Plants subjected to abiotic stress can acquire stress memory, enhancing their responsiveness to subsequent stress events. Despite its importance, the molecular mechanisms driving this process remain unclear. To gain insights into this, we conducted experiments involving cycles of salt stress and restoration treatments in maize.
First, to assess salt stress memory in maize, we examined proline and MDA contents and SOD and POD activities, which are common indicators of plant stress tolerance [43,44,45]. Elevated SOD and POD activities, along with higher proline content and lower MDA content, are indicative of enhanced stress tolerance. Within the same duration of salt treatment, we observed higher proline content, SOD, and POD activities in T2 plants (those that were exposed to a second round of salt stress) than in T1 plants (those that were first exposed to salt stress). Conversely, MDA content was lower in T2 plants than in T1 plants (Fig. 1C), suggesting an enhanced adaptive response to salt stress during the second salt stress event. These memory-associated changes are similar to what has been observed in other plants [46, 47], indicating the presence of a physiological memory effect in maize following the first stress event.
In addition to physiological changes, we also observed differences in gene expression profiles between plants with recurrent exposure to salt stress and those encountering salt stress for the first time (Fig. 2). Consistent with results reported in coffee and rice [15, 16], we observed a reduced number of DEGs in maize plants previously exposed to salt stress. Further functional analysis revealed several GO terms enriched in both the first and second salt stress responses, including the oxidation − reduction process, response to oxidative stress, trehalose biosynthetic process, and cellular ion homeostasis. These enriched GO terms have been demonstrated to play a crucial role in enhancing stress tolerance in plants [48,49,50], indicating their importance across different stress events. In addition, the stress resistance response might reduce photosynthesis efficiency, impacting short-term plant productivity, while heightening tolerance to subsequent stress events and promoting long-term productivity gains [47, 51]. This phenomenon was verified in the present study as we observed a perturbation in photosynthesis only during the first salt stress event (3h) that recovered during subsequent stress events (Fig. 2D). Furthermore, we observed significant differences in the response processes during the second salt stress compared to the first. These contrasting responses demonstrate the presence of transcriptional memory of salt stress exposure in maize.
Genes associated with stress memory have been identified through transcriptomic analysis in several plants [7, 11, 13, 16]. These genes are distinguished by their capacity to respond to the initial stress event and subsequently adjust their expression during repeated exposure to the same stress, enabling the plant to refine its responses to recurrent stress events [16]. For example, in Arabidopsis, over 2,000 genes have exhibited memory responses [11, 47, 52]. In this study, we employed RNA-seq to identify genes related to salt stress memory in maize, resulting in the identification of 2,213 SMGs that displayed distinct expression profiles in subsequent stress events compared to the initial event (Fig. 3A). Among these SMGs, the most prominent group comprised genes with high expression levels during the first stress event but reduced expression during the second stress event or vice versa ([+ / −] or [− / +]). This finding is consistent with that of previous studies on Arabidopsis, switchgrass, and rice [7, 16]. Genes displaying [+ / −] or [− / +] response patterns have been classified as “revised response” genes because they “revise” their initial behavior following subsequent exposure to the same stress [11, 13, 16]. This underscores the pivotal role of revised response in the plant stress memory mechanism.
TFs have been proposed as pivotal factors in determining memory behavior [53,54,55]. In this study, we observed that 11.5% of the identified SMGs in maize belong to the TF category (Fig. 4A), indicating the involvement of TFs in the establishment and orchestration of salt stress memory in maize. Analysis of the PPI network suggested that a complex regulatory network was responsible for inducing stress memory. Notably, we obtained nine hub TFs within the network (Fig. 4D), which act as robust candidates for salt memory factors that activate protective genes in response to subsequent salt stress. Interestingly, seven out of these nine TFs were from the WRKY family, which are known to play important roles in plant responses to various biotic and abiotic stresses [56,57,58]. Despite the relatively limited knowledge about their roles in salt stress memory, six of these TFs, namely, WRKY18 [41], WRKY33 [36, 37], WRKY40 [38], WRKY46 [34, 35], WRKY53 [39], and MYB15 [40], have been demonstrated to be involved in the response to salt stress. Therefore, further functional studies, particularly with regard to transcriptional memory, are warranted. Furthermore, we found that over 53% of all identified SMGs can be directly regulated by the nine identified TFs (Fig. 6A), providing evidence of their active involvement in the processes of salt stress memory in maize. However, there were differences in response patterns, such as MYB15 exhibited [+/-] expression, whereas WRKY18 exhibited [-/ +] expression. Together, our results indicate that salt stress memory can be established in maize and plants pre-exposed to salt stress exhibit better resistance to subsequent salt stresses.
Conclusions
Our study unveiled salt stress memory in maize, identifying 2,213 salt stress memory genes (SMGs) grouped into four distinct subcategories with altered responses to repeated stress. The most notable SMG groups exhibited revised responses to repeated stress, highlighting the significance of revised responses in plant stress memory. Moreover, we identified nine pivotal transcription factors (TFs) that govern these SMGs, offering potential genetic targets for enhancing salt stress resilience and crop productivity.
Materials and methods
Plant materials and growth conditions
Seeds of the maize inbred line B73 were germinated and grown in a greenhouse (13/11 h light/dark, 28/26 °C light/dark, and 50% relative humidity). Zhiying Zhu and Jinlei Han undertook the formal identification of the plant material used in the study, and no voucher specimens were collected. After the maize seedlings grew to the two leaves and one heart stage, they were transplanted and hydroponically cultivated in 1/2X Hoagland nutrient solution for five days. Thereafter, the plants were treated with 200 mM NaCl for 3 h (this corresponded to the first salt stress event, T1_3h). Non-stressed plants were simultaneously sampled as controls (T1_0h). Subsequently, half of the salt-stressed plants were incubated with 200 mM NaCl for another 45 h (a total exposure of 48 h, T1_48h). The remaining plants were transferred to 1/2X Hoagland nutrient solution for three days of recovery culture and subsequently exposed to salt stress (200 mM NaCl) for 3 h and 48 h (the second salt stress event, T2_3h and T2_48h) (Fig. 1A). Root tissues from both control and salt-treated plants were collected for further analysis.
Measurement of physiological parameters
Physiological parameters were measured using corresponding assay kits according to the manufacturer’s protocols. Specifically, the Proline Content Assay Kit (Solarbio, Cat#BC0290) was employed to quantify proline content using an acidic ninhydrin-based method. The absorbance was detected at a wavelength of 520 nm with a spectrophotometer. The Malondialdehyde (MDA) Content Assay Kit (Solarbio, Cat#BC0020) was utilized to measure MDA content, employing the thiobarbituric acid (TBA)-based colorimetric method, with absorbance measured at 532 nm. Superoxide Dismutase (SOD) activity was determined using a nitroblue tetrazolium colorimetric SOD Activity Detection Kit (Solarbio, Cat#BC0170), with absorbance measured at 560 nm. Peroxidase (POD) activity was assessed through the guaiacol method with the colorimetric POD Activity Detection Kit (Solarbio, Cat#BC0090), and absorbance measurements were conducted at 470 nm. All experiments were performed in triplicate.
RNA-seq and data analysis
Total RNA was extracted from root samples using the Omega Plant RNA kit (Omega Bio-tek, Cat. no. R6827-01) and used to make RNA-seq libraries with the Illumina TruSeq RNA Kit (NEB, Cat. no. E7530) following the manufacturer's instructions. The libraries were sequenced with the Illumina NovaSeq 6000 platform in the 150 bp paired-end mode. Three biological replicates were performed for the RNA-seq experiment. Sequencing raw reads were quality filtered and trimmed by using the Trim_Galore v.0.6.7 package (https://www.bioinformatics.babraham.ac.uk/projects/trim_galore/). Generated clean reads were then mapped to the maize reference genome using TopHat2 v.2.1.1 [59] with default settings. The genome sequence and annotation files for maize (RefGen_V4) were downloaded from Phytozome (https://phytozome-next.jgi.doe.gov). Cufflinks v.2.2.1 [60] was employed to calculate the gene expression level represented by FPKM (fragments per kilobase of transcript per million mapped reads) values. To identify differentially expressed genes (DEGs), featureCounts v.2.0.1 [61] was used to calculate the read counts for each gene, and differential expression analysis was performed using DESeq2 [62]. Genes with an adjusted p-value < 0.05 and |log2(fold change)|> 1 were considered differentially expressed. Gene Ontology (GO) enrichment analysis was performed using the Omicshare online tools (www.omicshare.com/tools).
To visualize the mapped reads, the BAM format files were converted to the bigwig format using the bamCoverage function in deepTools v.3.1.3 [63] with a bin size of 10 bp and RPKM normalization. Genome browser images were made using the Integrative Genomics Viewer (IGV) (https://software.broadinstitute.org/software/igv/) with bigwig files processed as described above.
Identification of salt stress memory genes
Salt stress memory genes (SMGs) were identified as described previously [7]. Briefly, DEGs between the T1_3h vs T1_0h and T2_3h vs T1_3h were compared, and four types of genes were considered to be SMGs: (1) those up-regulated in both T1_3h vs T1_0h and T2_3h vs T1_3h ([+ / +]), (2) down-regulated in both comparisons ([-/-]), and (3) up-regulated in T1_3h vs T1_0h and down-regulated in T2_3h vs T1_3h ([+/-]) or (4) vice versa ([-/ +]).
Protein–protein interaction network analysis
A protein–protein interaction (PPI) network was constructed based on the STRING database (https://string-db.org/) and visualized using Cytoscape v.3.9.1 [64]. To identify hub genes, Cytoscape cytohHubba [65] plugin with the Maximal Clique Centrality (MCC) method or the Molecular Complex Detection (MCODE) [66] plugin was utilized.
Defining predicted binding sites for TFs
The TF DNA binding motifs were downloaded from the PlantTFDB database (http://planttfdb.gao-lab.org/). Potential TF binding sites were determined by scanning motif occurrences in the region of interest by using the FIMO v.5.4.1 [67] with default settings. Candidate target genes for a given TF were defined as those with the TF binding motif in the promoter region, which was defined as 1 kb upstream of the transcription start site (TSS) of the gene.
Quantitative RT-PCR analysis
To generate a qRT-PCR template, the extracted RNA was reverse transcribed using the StarScript II First-strand cDNA Synthesis Mix With gDNA Remover (Genstar, Cat#A224-10). We performed qRT-PCR using 2X RealStar Green Fast Mixture (Genstar) on a CFX Connect Real-Time PCR Detection System (Bio-Rad). The 2−ΔΔCT method [68] was utilized for normalizing and calculating relative expression levels. Each qRT-PCR assay was repeated three times with three independent RNA preparations, and the maize EF1-α gene (NM_001112117) was used as an internal control for normalization [69, 70]. The nucleotide sequences of primers are listed in Table S7.
Availability of data and materials
The RNA-seq data described in this work have been deposited to the Genome Sequence Archive (GSA) database (http://gsa.big.ac.cn/) under the accession number CRA010351.
References
Chen HC, Cheng WH, Hong CY, Chang YS, Chang MC. The transcription factor OsbHLH035 mediates seed germination and enables seedling recovery from salt stress through ABA-dependent and ABA-independent pathways, respectively. Rice. 2018;11(1):50.
Zhang H, Guo J, Chen X, Zhou Y, Pei Y, Chen L, UlHaq S, Lu M, Gong H, Chen R. Pepper bHLH transcription factor CabHLH035 contributes to salt tolerance by modulating ion homeostasis and proline biosynthesis. Hortic Res. 2022;9:uhac203.
Zhu JK. Abiotic stress signaling and responses in plants. Cell. 2016;167(2):313–24.
Hoang XLT, Nhi DNH, Thu NBA, Thao NP, Tran LP. Transcription factors and their roles in signal transduction in plants under abiotic stresses. Curr Genomics. 2017;18(6):483–97.
Lim C, Kang K, Shim Y, Yoo SC, Paek NC. Inactivating transcription factor OsWRKY5 enhances drought tolerance through abscisic acid signaling pathways. Plant Physiol. 2022;188(4):1900–16.
Chen L, Yang H, Fang Y, Guo W, Chen H, Zhang X, Dai W, Chen S, Hao Q, Yuan S, et al. Overexpression of GmMYB14 improves high-density yield and drought tolerance of soybean through regulating plant architecture mediated by the brassinosteroid pathway. Plant Biotechnol J. 2021;19(4):702–16.
Zhang C, Peng X, Guo X, Tang G, Sun F, Liu S, Xi Y. Transcriptional and physiological data reveal the dehydration memory behavior in switchgrass (Panicum virgatum L.). Biotechnol Biofuels. 2018;11:91.
Ling Y, Serrano N, Gao G, Atia M, Mokhtar M, Woo YH, Bazin J, Veluchamy A, Benhamed M, Crespi M, et al. Thermopriming triggers splicing memory in Arabidopsis. J Exp Bot. 2018;69(10):2659–75.
Charng YY, Mitra S, Yu SJ. Maintenance of abiotic stress memory in plants: lessons learned from heat acclimation. Plant Cell. 2023;35(1):187–200.
Sharma M, Kumar P, Verma V, Sharma R, Bhargava B, Irfan M. Understanding plant stress memory response for abiotic stress resilience: molecular insights and prospects. Plant Physiol Biochem. 2022;179:10–24.
Ding Y, Liu N, Virlouvet L, Riethoven JJ, Fromm M, Avramova Z. Four distinct types of dehydration stress memory genes in Arabidopsis thaliana. BMC Plant Biol. 2013;13:229.
Virlouvet L, Fromm M. Physiological and transcriptional memory in guard cells during repetitive dehydration stress. New Phytol. 2015;205(2):596–607.
Ding Y, Virlouvet L, Liu N, Riethoven JJ, Fromm M, Avramova Z. Dehydration stress memory genes of Zea mays; comparison with Arabidopsis thaliana. BMC Plant Biol. 2014;14:141.
Virlouvet L, Avenson TJ, Du Q, Zhang C, Liu N, Fromm M, Avramova Z, Russo SE. Dehydration stress memory: gene networks linked to physiological responses during repeated stresses of zea mays. Front Plant Sci. 2018;9:1058.
Guedes FAdF, Nobres P, Rodrigues Ferreira DC, Menezes-Silva PE, Ribeiro-Alves M, Correa RL, DaMatta FM, Alves-Ferreira M. Transcriptional memory contributes to drought tolerance in coffee (Coffea canephora) plants. Environ Exp Bot. 2018;147:220–33.
do Amaral MN, Arge LWP, Auler PA, Rossatto T, Milech C, Magalhães AM Jr, Braga EJB. Long-term transcriptional memory in rice plants submitted to salt shock. Planta. 2020;251(6):111.
Chang YN, Zhu C, Jiang J, Zhang H, Zhu JK, Duan CG. Epigenetic regulation in plant abiotic stress responses. J Integr Plant Biol. 2020;62(5):563–80.
Friedrich T, Faivre L, Bäurle I, Schubert D. Chromatin-based mechanisms of temperature memory in plants. Plant Cell Environ. 2019;42(3):762–70.
Halder K, Chaudhuri A, Abdin MZ, Majee M, Datta A. Chromatin-based transcriptional reprogramming in plants under abiotic stresses. Plants. 2022;11(11):1449.
Nguyen NH, Vu NT, Cheong JJ. Transcriptional stress memory and transgenerational inheritance of drought tolerance in plants. Int J Mol Sci. 2022;23(21):12918.
Chen F, Ji X, Bai M, Zhuang Z, Peng Y. Network analysis of different exogenous hormones on the regulation of deep sowing tolerance in maize seedlings. Front Plant Sci. 2021;12:739101.
Luo M, Zhang Y, Li J, Zhang P, Chen K, Song W, Wang X, Yang J, Lu X, Lu B, et al. Molecular dissection of maize seedling salt tolerance using a genome-wide association analysis method. Plant Biotechnol J. 2021;19(10):1937–51.
Cao S, Wang L, Han T, Ye W, Liu Y, Sun Y, Moose SP, Song Q, Chen ZJ. Small RNAs mediate transgenerational inheritance of genome-wide trans-acting epialleles in maize. Genome Biol. 2022;23(1):53.
Rauf M, Arif M, Fisahn J, Xue GP, Balazadeh S, Mueller-Roeber B. NAC transcription factor speedy hyponastic growth regulates flooding-induced leaf movement in Arabidopsis. Plant Cell. 2013;25(12):4941–55.
Yadav UP, Ivakov A, Feil R, Duan GY, Walther D, Giavalisco P, Piques M, Carillo P, Hubberten HM, Stitt M, et al. The sucrose-trehalose 6-phosphate (Tre6P) nexus: specificity and mechanisms of sucrose signalling by Tre6P. J Exp Bot. 2014;65(4):1051–68.
Azari R, Reuveni M, Evenor D, Nahon S, Shlomo H, Chen L, Levin I. Overexpression of UV-DAMAGED DNA BINDING PROTEIN 1 links plant development and phytonutrient accumulation in high pigment-1 tomato. J Exp Bot. 2010;61(13):3627–37.
Thum KE, Shin MJ, Palenchar PM, Kouranov A, Coruzzi GM. Genome-wide investigation of light and carbon signaling interactions in Arabidopsis. Genome Biol. 2004;5(2):R10.
Zhu F, Luo T, Liu C, Wang Y, Zheng L, Xiao X, Zhang M, Yang H, Yang W, Xu R, et al. A NAC transcription factor and its interaction protein hinder abscisic acid biosynthesis by synergistically repressing NCED5 in Citrus reticulata. J Exp Bot. 2020;71(12):3613–25.
Kambona CM, Koua PA, Léon J, Ballvora A. Stress memory and its regulation in plants experiencing recurrent drought conditions. Theor Appl Genet. 2023;136(2):26.
Friedrich T, Oberkofler V, Trindade I, Altmann S, Brzezinka K, Lämke J, Gorka M, Kappel C, Sokolowska E, Skirycz A, et al. Heteromeric HSFA2/HSFA3 complexes drive transcriptional memory after heat stress in Arabidopsis. Nat Commun. 2021;12(1):3426.
Hu G, Chen J. A genome-wide regulatory network identifies key transcription factors for memory CD8+ T-cell development. Nat Commun. 2013;4:2830.
Bolton JL, Schulmann A, Garcia-Curran MM, Regev L, Chen Y, Kamei N, Shao M, Singh-Taylor A, Jiang S, Noam Y, et al. Unexpected transcriptional programs contribute to hippocampal memory deficits and neuronal stunting after early-life adversity. Cell Rep. 2020;33(11):108511.
Hirano Y, Ihara K, Masuda T, Yamamoto T, Iwata I, Takahashi A, Awata H, Nakamura N, Takakura M, Suzuki Y, et al. Shifting transcriptional machinery is required for long-term memory maintenance and modification in Drosophila mushroom bodies. Nat Commun. 2016;7:13471.
Ding ZJ, Yan JY, Xu XY, Yu DQ, Li GX, Zhang SQ, Zheng SJ. Transcription factor WRKY46 regulates osmotic stress responses and stomatal movement independently in Arabidopsis. Plant J. 2014;79(1):13–27.
Ding ZJ, Yan JY, Li CX, Li GX, Wu YR, Zheng SJ. Transcription factor WRKY46 modulates the development of Arabidopsis lateral roots in osmotic/salt stress conditions via regulation of ABA signaling and auxin homeostasis. Plant J. 2015;84(1):56–69.
Wang H, Li Z, Ren H, Zhang C, Xiao D, Li Y, Hou X, Liu T. Regulatory interaction of BcWRKY33A and BcHSFA4A promotes salt tolerance in non-heading Chinese cabbage [Brassica campestris (syn. Brassica rapa) ssp. chinensis]. Hortic Res. 2022;9:113.
Jiang Y, Deyholos MK. Functional characterization of Arabidopsis NaCl-inducible WRKY25 and WRKY33 transcription factors in abiotic stresses. Plant Mol Biol. 2009;69(1–2):91–105.
Huang Z, Wang J, Li Y, Song L, Chen D, Liu L, Jiang CZ. A WRKY protein, MfWRKY40, of resurrection plant myrothamnus flabellifolia plays a positive role in regulating tolerance to drought and salinity stresses of Arabidopsis. Int J Mol Sci. 2022;23(15):8145.
Zheng Y, Ge J, Bao C, Chang W, Liu J, Shao J, Liu X, Su L, Pan L, Zhou DX. Histone deacetylase HDA9 and WRKY53 transcription factor are mutual antagonists in regulation of plant stress response. Mol Plant. 2020;13(4):598–611.
Ding Z, Li S, An X, Liu X, Qin H, Wang D. Transgenic expression of MYB15 confers enhanced sensitivity to abscisic acid and improved drought tolerance in Arabidopsis thaliana. J Genet Genomics. 2009;36(1):17–29.
Chen H, Lai Z, Shi J, Xiao Y, Chen Z, Xu X. Roles of arabidopsis WRKY18, WRKY40 and WRKY60 transcription factors in plant responses to abscisic acid and abiotic stress. BMC Plant Biol. 2010;10:281.
Stone JD, Koloušková P, Sloan DB, Štorchová H. Non-coding RNA may be associated with cytoplasmic male sterility in Silene vulgaris. J Exp Bot. 2017;68(7):1599–612.
Hu L, Yao Y, Cai R, Pan L, Liu K, Bai L. Effects of fenclorim on rice physiology, gene transcription and pretilachlor detoxification ability. BMC Plant Biol. 2020;20(1):100.
Li Q, Jiang W, Jiang Z, Du W, Song J, Qiang Z, Zhang B, Pang Y, Wang Y. Transcriptome and functional analyses reveal ERF053 from Medicago falcata as key regulator in drought resistances. Front Plant Sci. 2022;13:995754.
Wang X, Zeng J, Li Y, Rong X, Sun J, Sun T, Li M, Wang L, Feng Y, Chai R, et al. Expression of TaWRKY44, a wheat WRKY gene, in transgenic tobacco confers multiple abiotic stress tolerances. Front Plant Sci. 2015;6:615.
Feng XJ, Li JR, Qi SL, Lin QF, Jin JB, Hua XJ. Light affects salt stress-induced transcriptional memory of P5CS1 in Arabidopsis. Proc Natl Acad Sci USA. 2016;113(51):E8335–43.
Li P, Yang H, Wang L, Liu H, Huo H, Zhang C, Liu A, Zhu A, Hu J, Lin Y, et al. Physiological and transcriptome analyses reveal short-term responses and formation of memory under drought stress in rice. Front Genet. 2019;10:55.
Zhou GA, Chang RZ, Qiu LJ. Overexpression of soybean ubiquitin-conjugating enzyme gene GmUBC2 confers enhanced drought and salt tolerance through modulating abiotic stress-responsive gene expression in Arabidopsis. Plant Mol Biol. 2010;72(4–5):357–67.
Li M, Chen R, Jiang Q, Sun X, Zhang H, Hu Z. GmNAC06, a NAC domain transcription factor enhances salt stress tolerance in soybean. Plant Mol Biol. 2021;105(3):333–45.
Rathor P, Borza T, Liu Y, Qin Y, Stone S, Zhang J, Hui JPM, Berrue F, Groisillier A, Tonon T, et al. Low mannitol concentrations in Arabidopsis thaliana expressing ectocarpus genes improve salt tolerance. Plants. 2020;9(11):1508.
Bruce TJA, Matthes MC, Napier JA, Pickett JA. Stressful “memories” of plants: evidence and possible mechanisms. Plant Sci. 2007;173(6):603–8.
Ding Y, Fromm M, Avramova Z. Multiple exposures to drought “train” transcriptional responses in Arabidopsis. Nat Commun. 2012;3:740.
Liu N, Ding Y, Fromm M, Avramova Z. Different gene-specific mechanisms determine the “revised-response” memory transcription patterns of a subset of A. thaliana dehydration stress responding genes. Nucleic Acids Res. 2014;42(9):5556–66.
Khan A, Khan V, Pandey K, Sopory SK, Sanan-Mishra N. Thermo-priming mediated cellular networks for abiotic stress management in plants. Front Plant Sci. 2022;13:866409.
Lakhina V, Arey RN, Kaletsky R, Kauffman A, Stein G, Keyes W, Xu D, Murphy CT. Genome-wide functional analysis of CREB/long-term memory-dependent transcription reveals distinct basal and memory gene expression programs. Neuron. 2015;85(2):330–45.
Guo C, Guo R, Xu X, Gao M, Li X, Song J, Zheng Y, Wang X. Evolution and expression analysis of the grape (Vitis vinifera L.) WRKY gene family. J Exp Bot. 2014;65(6):1513–28.
Su T, Xu Q, Zhang FC, Chen Y, Li LQ, Wu WH, Chen YF. WRKY42 modulates phosphate homeostasis through regulating phosphate translocation and acquisition in Arabidopsis. Plant Physiol. 2015;167(4):1579–91.
Chen C, Chen X, Han J, Lu W, Ren Z. Genome-wide analysis of the WRKY gene family in the cucumber genome and transcriptome-wide identification of WRKY transcription factors that respond to biotic and abiotic stresses. BMC Plant Biol. 2020;20(1):443.
Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009;25(9):1105–11.
Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, Salzberg SL, Wold BJ, Pachter L. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010;28(5):511–5.
Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014;30(7):923–30.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.
Ramírez F, Ryan DP, Grüning B, Bhardwaj V, Kilpert F, Richter AS, Heyne S, Dündar F, Manke T. deepTools2: a next generation web server for deep-sequencing data analysis. Nucleic Acids Res. 2016;44(W1):W160–5.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.
Chin CH, Chen SH, Wu HH, Ho CW, Ko MT, Lin CY. cytoHubba: identifying hub objects and sub-networks from complex interactome. BMC Syst Biol. 2014;8(Suppl 4):S11.
Bader GD, Hogue CW. An automated method for finding molecular complexes in large protein interaction networks. BMC Bioinformatics. 2003;4:2.
Grant CE, Bailey TL, Noble WS. FIMO: scanning for occurrences of a given motif. Bioinformatics. 2011;27(7):1017–8.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method. Methods. 2001;25(4):402–8.
Chen Z, Wang HC, Shen J, Sun F, Wang M, Xu C, Tan BC. PPR-SMR1 is required for the splicing of multiple mitochondrial introns, interacts with Zm-mCSF1, and is essential for sesed development in maize. J Exp Bot. 2019;70(19):5245–58.
Wei S, Xia R, Chen C, Shang X, Ge F, Wei H, Chen H, Wu Y, Xie Q. ZmbHLH124 identified in maize recombinant inbred lines contributes to drought tolerance in crops. Plant Biotechnol J. 2021;19(10):2069–81.
Acknowledgements
We are grateful for the funding support.
Funding
This work was supported by the National Natural Science Foundation of China (32070544 and 32200543), the Natural Science Foundation of Jiangsu Province (BK20220604), and the Science and Technology Program of Nantong City (JC12022064).
Author information
Authors and Affiliations
Contributions
JH, KW, and BW planned and designed the research. ZZ, YD, GY, XZ, and HZ performed the experiments. JH, ZZ, QC, and XK analyzed the data. JH, KW, BW, BZ, GR and EM wrote and revised the manuscript. All authors read and approved the final version for publication.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
The seeds of B73 plants were provided by Prof. Kai Wang at Nantong University. The appropriate permissions were obtained for all materials used in this study. We complied with all relevant institutional, national, and international guidelines and legislation.
Consent for publication
Not applicable.
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
Additional file 1: Table S1.
Sequencing and mapping statistics for the RNA-seq data.
Additional file 2: Table S2.
Differentially expressed genes (DEGs) between the control and salt stressed samples.
Additional file 3: Table S3.
GO enrichment in biological process of DEGs.
Additional file 4: Table S4.
Salt stress memory genes (SMGs) identified in maize.
Additional file 5: Table S5.
GO enrichment in biological process of SMGs.
Additional file 6: Table S6.
Transcription factor involved in salt stress memory.
Additional file 7: Table S7.
List of the primers for qRT-PCR analysis of the genes.
Additional file 8: Figure S1.
PCR Efficiency.
Additional file 9: Figure S2.
Melting curves of 10 genes were obtained using CFX Manager (Bio-Rad).
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, 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 changes were made. 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/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Zhu, Z., Dai, Y., Yu, G. et al. Dynamic physiological and transcriptomic changes reveal memory effects of salt stress in maize. BMC Genomics 24, 726 (2023). https://doi.org/10.1186/s12864-023-09845-w
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12864-023-09845-w