- Grassland Agri-Husbandry Research Center, College of Grassland Science, Qingdao Agricultural University, Qingdao, China
Leaves are the most critical portion of forage crops such as alfalfa (Medicago sativa). Leaf senescence caused by environmental stresses significantly impacts the biomass and quality of forages. To understand the molecular mechanisms and identify the key regulator of the salt stress-induced leaf senescence process, we conducted a simple and effective salt stress-induced leaf senescence assay in Medicago truncatula, which was followed by RNA-Seq analysis coupled with physiological and biochemical characterization. By comparing the observed expression data with that derived from dark-induced leaf senescence at different time points, we identified 3,001, 3,787, and 4,419 senescence-associated genes (SAGs) for salt stress-induced leaf senescence on day 2, 4, and 6, respectively. There were 1546 SAGs shared by dark and salt stress treatment across the three time points. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses showed that the 1546 SAGs were mainly related to protein and amino acids metabolism, photosynthesis, chlorophyll metabolism, and hormone signaling during leaf senescence. Strikingly, many different transcription factors (TFs) families out of the 1546 SAGs, including NAC, bHLH, MYB, and ERF, were associated with salt stress-induced leaf senescence processes. Using the transient expression system in Nicotiana benthamiana, we verified that three functional NAC TF genes from the 1546 SAGs were related to leaf senescence. These results clarify SAGs under salt stress in M. truncatula and provide new insights and additional genetic resources for further forage crop breeding.
Introduction
In nature, leaf senescence is an inevitable process by which plants reallocate nutrients, thus altering their metabolism and energy cycles to redirect more nutrients from senescent or dying leaves to new leaves, flower development, seed maturation, or other tissues and storage organs (Hörtensteiner and Feller, 2002; Yoshimoto et al., 2009; Bar-Dror et al., 2011). Therefore, leaf senescence is a key developmental process that ensures optimal plant performance and directly enhances plant adaptability and viability in the face of stress (Peleg and Blumwald, 2011). However, premature senescence caused by stress can significantly result in yield loss and food quality reduction (Uauy et al., 2006; Yamada et al., 2009; Gregersen et al., 2013). Through delaying the senescence of functional leaves by 1 day, it has been calculated that yields can be increased by 2–10% in major crops such as maize, soybeans, cotton, rice, and wheat (Buchanan-Wollaston et al., 2003; Liu et al., 2015; Schippers et al., 2015). As for forage crops such as alfalfa, most of the nutrients are indeed stored in the leaves. Therefore, leaf senescence has a more direct impact on the production of high-quality forage (Du et al., 2020). As a result, many researchers focused on improving crop yields and agricultural productivity by delaying the leaf senescence process (Breeze, 2020; Bian et al., 2021; Li et al., 2021).
Leaf senescence is a highly complex process that is tightly regulated at multiple levels, including the chromatin level, the mRNA level (via transcriptional and posttranscriptional regulation), and the protein level (via translational and posttranslational regulation). It has been noted that along with leaf senescence, thousands of leaf senescence-associated genes (SAGs) such as transcription factors (TFs) are induced (Buchanan-Wollaston, 1997; Woo et al., 2001). A wide range of TF families (e.g., NAC, WRKY, CDF, bZIP, and MYB) has been shown to participate, often critically, in the regulation of senescence in plants (Zhou et al., 2011b; Rushton et al., 2012; Jaradat et al., 2013; Kim et al., 2013; Mao et al., 2017; Woo et al., 2019; Xu et al., 2020). These regulatory networks are interconnected with multiple crosstalk pathways, thus enabling responses to and regulation on various internal and external environmental senescence signals (Buchanan-Wollaston et al., 2005; Lim et al., 2007; Park et al., 2007; Breeze et al., 2011; Allu et al., 2014). Unfavorable environmental conditions include prolonged drought, darkness, salinity, heat, cold, pests, and diseases (Wingler et al., 2006; Lim et al., 2007; Parlitz et al., 2011; He et al., 2018). Salt stress is one of the common abiotic stresses limiting plant growth and development, and it has become a severe problem restricting crop production as well as ecological environment construction (Lutts et al., 1996; Munné-Bosch and Alegre, 2004). Salt stress can induce senescence in plant leaves, but the regulatory mechanisms are not fully understood (Munns and Tester, 2008; Balazadeh et al., 2010; Allu et al., 2014).
Over the last decade, significant progress has been made in understanding the key molecular principles of leaf senescence in model plants including Arabidopsis thaliana (Lim et al., 2003; Lin and Wu, 2004; Buchanan-Wollaston et al., 2005; Breeze et al., 2011), rice (Lee et al., 2001), tobacco (Pageau et al., 2006) and other crops, such as wheat (Uauy et al., 2006), corn (Zhang et al., 2014), and cotton (Kong et al., 2013). Few such studies have been performed in leguminous forage species. Medicago truncatula, a close relative of alfalfa (Medicago sativa), has been developed as an excellent model legume system, especially for studying forage crop improvements. To elucidate salt-induced leaf senescence in a time-saving and well-targeted way, salt stress-induced detached leaf senescence assays have been applied as useful tools for studying senescence in Arabidopsis and rice (Sakuraba et al., 2018).
To discover the key regulators of leaf senescence in forage crops, in this study, we examined salt- and dark-induced leaf senescence in Medicago truncatula through collecting and analyzing physiological, biochemical, and transcriptional data over the time course of leaf senescence. This work presents detailed expression profiles of leaf senescence induced by salt stress in Medicago truncatula. Many key active regulatory networks and senescence signaling pathways were illustrated. Several important TFs involved in this process were also characterized using the transient expression system in Nicotiana benthamiana. Our work reveals new molecular regulatory mechanisms of leaf senescence in Medicago truncatula and provides a valuable genetic resource to improve the biomass and quality of forage crops under salt stress.
Materials and Methods
Plant Material and Stress Treatments
Medicago truncatula ecotype R108 was used as the experimental plant material. For seed germination, we used sandpaper to lightly scratch the coats of Medicago truncatula seeds that had already been vernalized for 2 days. Seeds were then transferred to dishes with moistened filter paper and incubated in a light incubator for 7 days before being transferred to hydroponic growth cultivation in Hoagland's nutrient solution; during this period the culture medium was changed every 3 days. Plants were placed in an artificial climate incubator with a 16 h photoperiod, day/night temperatures of 25°C/22C, and a relative humidity of 60–70%.
Leaves were sampled when plants were 34 days old, and the sampling site was the third compound leaf of the Medicago truncatula plant (Supplementary Figure 1). The removed leaves were immediately placed in Petri dishes containing different treatment solutions. The Petri dishes were then placed in light or dark conditions according to the experimental design, with the same growth conditions as those generally used for the hydroponic growth of plants. The salt stress treatment fluid was the same as that used for the control treatments (2.4 g of Murashige–Skoog medium, 3 mM MES buffer, pH 5.8), except for its 100, 150, or 200 mM NaCl concentrations achieved by addition of NaCl.
We collected all samples on day 0, 2, 4, and 6 and stored them at −80°C after quickly freezing them in liquid nitrogen. Sampled materials were used to assay physiological indicators and for transcriptomic sequencing.
Physiological Measurements
To determine the chlorophyll content, 0.1 g samples of previously frozen leaves were weighed, frozen, ground, and placed in 96% (v/v) acetone. The 3-h extraction was performed at room temperature under dark conditions, and then, samples were centrifuged at 4°C 6,000 rpm for 30 s, then the supernatant was collected. The chlorophyll content in the supernatant was quantified spectrophotometrically (Zhang and Guo, 2018).
The H2O2 contents were determined spectrophotometrically, while the ABA content was determined by high performance liquid chromatography.
At least three biological replicates were analyzed for each of the above trials. All data were subjected to single-factor ANOVA using SPSS 26 (IBM Corp., Armonk, NY, USA). The differences among the different treatment means were evaluated using Duncan's multiple range test, and P < 0.05 were considered statistically significant. All charts were created with Excel 2019 (Microsoft Corp., Redmond, WA, USA).
RNA Isolation, Quantification, and Quality Controls
Total RNA was extracted from the tissue using TRIzol® Reagent (Plant RNA Purification Reagent for plant tissue; Invitrogen Corp., Carlsbad, CA, USA) according the manufacturer's instructions and genomic DNA was removed using DNase I (TaKaRa, Kusatsu, Japan). Then, RNA quality was determined with a 2100 Bioanalyser (Agilent; Santa Clara, CA, USA) and quantified using the ND-2000 platform (NanoDrop Technologies, Wilmington, DE, USA). Only high-quality RNA samples [OD260/280 = 1.8–2.2, OD260/230 ≥ 2.0, RNA integrity number (RIN) ≥ 6.5, 28S:18S ratio ≥ 1.0, >2 μg samples] were used to construct sequencing libraries.
cDNA Library Construction and RNA Sequencing
RNA-Seq transcriptome libraries were prepared using the TruSeqTM RNA sample preparation Kit from Illumina (San Diego, CA, USA) with 1 μg of total RNA. In brief, messenger RNA was first isolated using the polyA selection method with oligo(dT) beads and then fragmented with fragmentation buffer. Second, double-stranded cDNA was synthesized using a SuperScript double-stranded cDNA synthesis kit (Invitrogen) with random hexamer primers (Illumina). Then, the synthesized cDNAs were subjected to end-repair, phosphorylation, and “A” base addition according to Illumina's library construction protocol. Libraries were size selected for 200–300 bp cDNA target fragments using 2% Low Range Ultra Agarose electrophoresis followed by PCR amplified using Phusion DNA polymerase (New England Biolabs, Ipswich, MA, USA) for 15 cycles. After quantification with a TBS380 fluorometer (Turner BioSystems, Sunnyvale, CA, USA), the paired-end RNA-seq sequencing library was sequenced with the Illumina HiSeq xten/NovaSeq 6000 sequencer (2 × 150 bp read length).
The raw paired end reads were trimmed and quality-filtered by Fastp version: 0.19.5 (https://github.com/OpenGene/fastp) for acquiring clean reads.
All the obtained high-quality and clean reads were mapped against the reference genome of Medicago_truncatula (reference genome version is MedtrA17_4.0; Reference genome source can be reached via http://plants.ensembl.org/Medicago_truncatula/Info/Index) by Mapping tools hisat2 (Version 2.1.0, http://ccb.jhu.edu/software/hisat2/index.shtml), and the matching rate ranged from 6.44 to 90.06%.
To evaluate the similarity of RNA-Seq data from different samples, the mapped RNA-Seq data from different samples were first subjected to a standardized analysis (Z-score) and then to a principal component analysis (PCA) with the analysis software Python.
All information related to project sample information, raw data, quality control data, and comparison data can be found in Supplementary Table 1. Reference genome annotation information are listed in Supplementary Table 2.
Differential Gene Expression Analysis and Functional Enrichment
To identify differentially expressed genes (DEGs) between two different samples, the expression level of each transcript was calculated according to the transcripts per million reads (TPM) method. The expression quantification software tool RSEM (RNA-Seq by Expectation-Maximization, http://deweylab.github.io/RSEM/) was applied to quantify gene abundances for each group and each time point, respectively, using the average transcripts per million reads (TPM) values of biological replicates as a quantitative metric.
R statistical package software EdgeR (Empirical analysis of Digital Gene Expression in R, (http://www.bioconductor.org/packages/2.12/bioc/html/edgeR.html) was adapted to perform differential expression analysis. The threshold values for screening DEGs were |log2FC| ≥ 1 or ≤ −1 & P-adjust < 0.05. The method of multiple test correction is BH (FDR correction with Benjamini and Hochberg). In order to control the probability or frequency of error in the overall inference result, the p-value obtained from the statistical test is corrected, i.e., multiple test correction is performed, and the corrected p-value is called P-adjust.
In addition, Enrichment analysis using Fisher's precision test, functional-enrichment analysis including GO and KEGG enrichment analyses, were performed to identify which DEGs were significantly enriched within specific GO terms and metabolic pathways at a BH-corrected P-value (p-adjust) < 0.05 threshold compared with the whole-transcriptome background. GO functional enrichment and KEGG pathway analysis were conducted using Goatools (https://github.com/tanghaibao/Goatools) and KOBAS (http://kobas.cbi.pku.edu.cn/home.do), respectively. The GO and KEGG enrichment analysis of senescence-associated genes (SAGs) is detailed in Supplementary Tables 3, 4, respectively. Specific information on the software names, version details, and URL links used in this study are shown in Supplementary Table 5.
Quantitative Real-Time PCR Analysis
In order to verify the expression patterns of genes observed in the RNA-seq analysis, we used the total RNA of all samples subjected to RNA-Seq analysis for qRT-PCR. Sixteen DEGs were selected for qRT-PCR detection, and the MtACTIN gene was used as a reference gene. RNAs were reverse transcribed into cDNAs using M5 Super Plus qPCR RT kit with gRNA remover (mei5 Biotech Co., Ltd; Beijing, China). ChamQ SYBR color qPCR Master Mix (Vazyme Biotech Co., Ltd; Nanjing, China) was used for qRT-PCR, and three technical repeats were performed for each reaction. All primers were designed using Primer 5.0 software and are listed in Supplementary Table 6.
Functional Verification of Genes by Agrobacterium-Mediated Transient Expression
We selected 9 genes for transient expression in tobacco (Nicotiana benthamiana), among which SGR was selected as a known positive control. Using the cDNAs from the cDNA library as templates, nine target genes were cloned and sequenced to verify the sequence accuracy. Then, using seamless DNA cloning (ClonExpress® Ultra One Step Cloning Kit, Vazyme Biotech Co., Ltd), the target genes were ligated to the pBI121 vector. The ligation products were transformed into E. coli strain DH5α (Vazyme Biotech Co., Ltd), and a positive E. coli colony was selected and verified by sequencing. The plasmid was extracted and then transformed into Agrobacterium tumefaciens, and a single positive Agrobacterium tumefaciens colony was chosen for Nicotiana benthamiana injection using the method described by Li et al. (2008).
The software tool CE Design (http://www.vazyme.com) was used to design specific amplification primers (Supplementary Table 6).
Transcription Factor Family Analysis
TFs usually consist of a DNA-binding domain (DBD), a transcription-activating domain (TAD), and other TF-binding domains (e.g., signal sensing domain). Depending on their functional domains, TFs can be grouped into different families. TFs in the same family generally contain the same DBD and different TADs. Analysis of the domain information contained in gene transcription products was used for TF prediction and family analysis of genes. We used Plant TFDB 4.0 (http://planttfdb.cbi.pku.edu.cn/) for TF analysis. For blastp, an E < 1e-5 threshold was used. For a detailed table of TF statistics, please see Supplementary Table 7. We used the software FIMO to predict the target genes of TF MTR_3g088110 based on the information of transcription factor binding sites included in the Plant TFDB database. We screened target genes with high correlation (correlation> 0.5, P < 0.05) with TF MTR_3g088110 from 3001,3787,4419 SAGs (Figures 4A–C). The information related to the target genes obtained at 2, 4, and 6 days of prediction is shown Supplementary Table 8.
Short Time-Series Expression Miner and Functional Similarity Analyses
Time-series analysis (time sequence expression analysis) involves the comparison of gene expression on the same (or similar) samples at different time points, in order to observe changes in gene expression between each time point, and to elucidate the patterns of interdependence between genes. Short Time-series Expression Miner (STEM, Version 1.3.11) with a p < 0.05 threshold was used for time-series expression trend analysis (Supplementary Table 9).
Microarray Significant Profiles Analysis
The maSigPro was applied to perform the analysis of gene expression data at different time point. MaSigPro (Version 1.56.0, http://www.bioconductor.org/packages/release/bioc/html/maSigPro.html) is an R package for analyzing time series data based on multiple linear regression models to fit the time, experimental conditions, and gene expression, then apply a stepwise regression method to find the best combination of independent variables. Time-series difference analysis based on maSigPro was used to obtain the genes that have changed differentially throughout the sample of time series nodes (Supplementary Table 10).
Results
Physiological and Biochemical Changes During Salt-Induced Leaf Senescence
To illustrate the effects of salt stress on leaf senescence, we collected fully expanded single leaves without senescence signs from 5-week-old Medicago truncatula plants (Supplementary Figure 1A) as experimental materials. Different concentrations of salt (NaCl) treatments were investigated in this experiment. As the salt concentration increased, the leaf senescence developed more dramatically (Supplementary Figure 1B). Based on our preliminary experimental results, the 150 mM NaCl treatment was observed to induce leaf senescence similar to that observed under the dark treatment. Dark-induced leaf senescence has been reported to be similar to age-dependent leaf senescence when plants develop normally (Sobieszczuk-Nowicka et al., 2018).
To remove the background effect of the development process, we established a light control group (control-light). To better screen the senescence-related genes, we set up a dark treatment group as a positive control in the salt-induced senescence experiment performed on isolated leaves. Thus, this experimental design could exclude the negative control conditions and be compared with the positive control senescence conditions. The salt-induced senescence experiments were performed under conditions with 150 mM NaCl and light. It was observed that the detached leaves showed different senescence phenotypes across the different groups (Figure 1). In the salt-induced group, the leaves started to turn yellow on day 4 and gradually turned 50–90% yellow with eroded, partially blackened, or transparent leaves on day 6. In the dark-induced group, leaves showed signs of yellowing on day 4 and increased yellowing on day 6; however, in the light control group, entire leaves stayed green and showed no signs of senescence until day 6.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/666660/fpls-12-666660-HTML/image_m/fpls-12-666660-g001.jpg)
Figure 1. Salt stress-induced leaf senescence in Medicago truncatula. The leaf color and senescence progression of detached M. truncatula leaves under dark (Dark), normal light (Control-light), and 150 mM NaCl (150 mM NaCl-light) conditions for 0, 2, 4, and 6 days. DT, day(s) of treatment. Scale bar, 1 cm.
We investigated the physiological and biochemical responses of in vitro leaves treated with 150 mM NaCl through measuring chlorophyll, H2O2, and abscisic acid (ABA) contents. Chlorophyll is commonly used as a marker and representative indicator of plant senescence to evaluate the leaf senescence process. The chlorophyll content decreased significantly in both the salt-induced group and the dark-induced group and was only slightly degraded in the light control group (Figure 2A). H2O2 is an important biochemical indicator that reflects the emergency response of oxidation under stress. H2O2 content showed a complex profile throughout leaf senescence. In the salt-induced group, H2O2 content increased initially up to its maximum on day 4 and then decreased. The levels of H2O2 in the dark-induced group and salt-induced group followed similar trends (Figure 2B). ABA content increased during leaf senescence in both the salt or dark-induced groups and the control group. In the salt-induced group, ABA content peaked on day 2 and then decreased, but remained above the initial value (Figure 2C).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/666660/fpls-12-666660-HTML/image_m/fpls-12-666660-g002.jpg)
Figure 2. Physiological analysis of salt-induced leaf senescence in Medicago truncatula. Levels of (A) chlorophyll (a+b), (B) hydrogen peroxide (H2O2), and (C) abscisic acid (ABA) in detached Medicago truncatula leaves exposed to different conditions during leaf senescence. The values are presented as the mean ± SE values of three independent biological replicates per time point. Different letters indicate significant differences among treatments according to an analysis of variance (ANOVA, P < 0.05). Error bars correspond to standard errors. FW indicates fresh weight.
Transcriptomic Analysis Identified Thousands of Differentially Expressed Senescence-Associated Genes at Different Time Points
We collected leaves sampled on day 0, 2, 4, and 6 of treatment for sequencing, with three biological replicates for each group at each time point, thus obtaining a total of 222.57 Gb of clean data (Supplementary Table 1) from the 10 groups (30 samples in total) of leaf samples. The clean data for each sample exceeded 5.91 Gb. Approximately 85% of clean reads matched the reference genomic sequence, and the Q30 base percentage was above 91.87%.
Principal component analysis (PCA) showed that replicate samples from the same group were tightly clustered together, suggesting that each biological replicate sample had low variability. Among different groups, the transcriptomic data differed substantially across time points, suggesting that the salt stress treatment caused significantly differential expression of genes (Figure 3A). Validation of RNA-seq data by quantitative real-Time PCR (qRT-PCR) analysis revealed similar expression profiles (Supplementary Figure 2). The observed high inter-sample correlation also validated the experimental design (Figure 3B). Using the software edgeR with screening thresholds of |log2FC| ≥ 1 or ≤ −1 together with P < 0.05, the differentially expressed genes (including up- and down-regulated genes) obtained from each group compared to the control group were analyzed (Figure 3C, Supplementary Table 11).Compared with the 0 day time point, both salt- and dark-induced groups had large amounts of up- and down-regulated genes, ranging from 3,700 to 6,700, while in the light control group differentially expressed genes numbered ~2,000 for both up- and down-regulated genes (Figure 3C).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/666660/fpls-12-666660-HTML/image_m/fpls-12-666660-g003.jpg)
Figure 3. Transcriptomic overview of salt stress-induced detached leaf senescence in Medicago truncatula. (A) Principal component analysis (PCA) plot of transcriptome profiles from different conditions. (B) Correlation analysis showed the direct association between samples in the same treatment. (C) The number of differentially expressed genes that were up- or down-regulated during leaf senescence compared with 0 day control prior to treatment, as assessed using the difference analysis software edgeR at thresholds of |log2FC| ≥ 1 or ≤ −1 & P-adjust < 0.05.
To obtain genes associated with leaf senescence, we generated Venn diagrams of the differentially expressed genes at each time point under different treatments to identify the genes that were induced in both the salt- and dark-induced groups but were not identified in the light control group. This analysis identified 3,001, 3,787, and 4,419 senescence-associated genes (SAGs) on day 2, 4, and 6, respectively (Figures 4A–C, Supplementary Table 12). Then, we assessed the Venn diagrams of the differentially expressed SAGs shared by salt- and dark-induced groups at three time points, thus obtaining 1546 SAGs shared by the three time points (Figure 4D).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/666660/fpls-12-666660-HTML/image_m/fpls-12-666660-g004.jpg)
Figure 4. Differentially expressed senescence-associated genes (SAGs) analysis. Differentially expressed SAGs in detached leaves under darkness and 150 mM NaCl treatment on day 2 (A), 4 (B), and 6 (C). (D) Venn diagram of SAGs shared by the darkness and 150 mM NaCl treatments on day 2 (A, 3001), 4 (B, 3787), and 6 (C, 4419). There were 1546 SAGs shared by darkness and 150 mM NaCl treatments across three time points.
Gene Ontology Enrichment and Kyoto Encyclopedia of Genes and Genomes Pathway Enrichment Analysis of SAGs Throughout Leaf Senescence
We investigated 1546 differentially expressed SAGs by using their GO classification. According to the GO enrichment analysis results, differentially expressed genes were mainly enriched within the cellular carbohydrate biosynthetic process, plastid membrane organization, sulfur compound biosynthetic process, and small molecule catabolic process categories (Table 1). The directed acyclic plots of GO terms revealed that senescing leaves subjected to salt stress in vitro were significantly affected in terms of the genes associated with the structure and function of intracellular chloroplasts and thylakoids, resulting in a dramatic decrease in leaf photosynthesis (Supplementary Figure 3).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/666660/fpls-12-666660-HTML/image_m/fpls-12-666660-t001.jpg)
Table 1. Gene ontology enrichment and Kyoto encyclopedia of genes and genomes pathway enrichment analysis of SAGs (Figure 4D, 1546 genes) shared by darkness and salt-induced leaf senescence at three time points.
KEGG pathway enrichment analysis of the 1546 differentially expressed SAGs indicated that there were 633 SAGs associated with 121 KEGG pathways, with 3 pathways showing significant expression differences throughout this process (P-adjust < 0.05). The most significantly enriched pathways were “Porphyrin and chlorophyll metabolism” (13 genes), “Glyoxylate and dicarboxylate metabolism” (16 genes), and “Photosynthesis” (18 genes). Throughout the salt-stress period, the involved metabolic processes were mostly associated with nutrient recycling and photosynthesis, which correspond with the progression of the leaf senescence process (Table 1).
In addition, we also performed separate GO enrichment and KEGG pathway enrichment analyses for differentially expressed SAGs shared by salt- and dark-induced groups at three time points (3,001, 3,787, and 4,419 SAGs on day 2, 4, and 6, respectively), to assess these genes in more detail (Supplementary Figures 4, 5). The analysis showed that the changes in photosynthesis and nutrient recycling began to occur on day 2 of senescence in salt-stressed leaves. The metabolism of secondary substances was enhanced, with increased expression of many antioxidants and antioxidant enzymes involved in the senescence process associated with changes in plant hormone signal transduction. The genes involved in photosynthesis, light reactions, and carbon fixation in non-senesced leaves were down-regulated as leaves senesced, along with an enhancement of genes associated with amino acid metabolism, lipid metabolism, nutrient recycling, and secondary material metabolism over time.
Dynamic SAG Expression During Leaf Senescence
In order to analyze the dynamic expression of 1546 SAGs at different time points (day 0, 2, 4, and 6) during leaf senescence, we used Short Time-series Expression Miner (STEM) for temporal pattern analysis. There were three main gene expression patterns observed (P < 0.05), including two up-regulated patterns (blue, profiles 21 and 24), four down-regulated patterns (orange-red, profiles 0, 1, 3, and 4), and one other pattern (green profiles 6) (Figure 5). Each profile contains genes with the same expression trend.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/666660/fpls-12-666660-HTML/image_m/fpls-12-666660-g005.jpg)
Figure 5. Short Time-series Expression Miner (STEM) analysis of darkness and salt treatment-enriched 1546 SAGs across three time points. Each square box indicates a type of expression profile, with the profile order on the upper left and the P-value on the bottom left. Only the significantly enriched cluster profile with a P-adjust < 0.05 threshold is shown.
The KEGG pathway enrichment analysis of these three clusters revealed that the metabolic pathways of these genes differed among clusters. The genes from the up-regulated patterns of profiles 21 and 24 were associated with the organic substance catabolic process and amino acid metabolism, including cysteine & methionine metabolism and arginine & proline metabolism; in contrast, SAGs with down-regulated patterns corresponding to profiles 0, 1, 3, and 4 were enriched mainly in energy metabolism and carbohydrate metabolism (Supplementary Figure 6).
Cluster Analysis
Based on MaSigPro (Version 1.56.0) analysis, we obtained the transcriptome profile of salt stress-induced leaf senescence for 2, 4, and 6 days in Medicago truncatula, and the results revealed 601 (9 clusters) SAGs which from 1546 SAGs (Figure 4D). The actual gene expression differences between the control-light groups and 150 mM NaCl-light are showed in Figure 6. Compared to the control-light group, the SAGs expression in each cluster of the 150 mM NaCl-light group varied within a specific pattern. These up-regulated genes were classified into six groups (clusters 1, 2, 5, 6, 8, and 9), while the down-regulated genes were classified into three groups (clusters 3, 4, and 7). It is noteworthy that 358 genes (59.6%) were enriched into these three clusters. Clustering of expression patterns of genes in all samples is shown in Supplementary Figure 7.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/666660/fpls-12-666660-HTML/image_m/fpls-12-666660-g006.jpg)
Figure 6. Data visualization with maSigPro (https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4155246/). Each line chart shows the average expression profile of the gene clusters. Horizontal axis represents different time points, vertical axis represents average expression values. Different colored lines represent different groups (red and green lines represent the control-light and 150 mM NaCl-light groups, respectively). Dotted lines represent Fitted curves of the gene expression at each time point.
Transcription Factor Analysis
We used Plant TFDB 4.0 to perform TF prediction and statistical analysis of the 1546 SAGs identified. There were 146 TFs among the 1546 SAGs annotated by Plant TFDB 4.0, and the top three representative transcription families with the most members were bHLH (eight members), MYB (seven members), and ERF (seven members). The transcription factor families NAC, B3, CO-like, Dof, MYB-related, WRKY, and bZIP each had three members included (Figure 7A). Most of these TFs were down-regulated throughout the senescence process, but the NAC TF family was up-regulated (Supplementary Table 13). GO enrichment analysis revealed that the TFs were involved in biological regulation, organelle, cell part and transcription regulator activity categories (Figure 7B). Furthermore, we identified 351, 440, and 545 TFs among the significantly differentially expressed SAGs shared by darkness and salt treatment on day 2, 4, and 6, respectively (Figures 7C–E, Supplementary Table 7). The most differentially expressed senescence-associated TFs belonged to the bHLH, ERF, MYB, NAC, and WRKY TF families, which are the major families known to be involved in the regulation of senescence.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/666660/fpls-12-666660-HTML/image_m/fpls-12-666660-g007.jpg)
Figure 7. Transcription factor families of the senescence-associated genes (SAGs) on day 2, 4, and 6 shared by the darkness and salt treatments. Significantly enriched transcription factor families were identified by Hmmscan E-value (1.0E-5) and Blast E-value (1.0E-5). (A) The transcription factor families from 1546 darkness and salt treatment SAGs across three time points. (B) Scatter diagram of transcription factor families in (A) based on the Gene Ontology enrichment (GO) enrichment analysis. (C) Transcription factor families in the 3,001 darkness and salt treatment SAGs on day 2. (D) Transcription factor families in the 3,787 darkness and salt treatment SAGs on day 4. (E) Transcription factor families in the 4,419 darkness and salt treatment SAGs on day 6.
Gene Function Verification by Agrobacterium-Mediated Transient Expression
To identify the functions of differentially expressed SAGs in our results, we performed a transient overexpression assay in Nicotiana benthamiana. Based on the expression patterns of the 146 TFs among the different clusters (Figure 5), eight senescence-associated TFs were chosen, and their constructs were then agro-infiltrated into Nicotiana benthamiana leaves (Table 2).
We found that the green leaf spots infiltrated by the number 5, 7, 8, and SGR constructs (MTR_4g081870, MTR_6g012670, MTR_3g088110, and MTR_5g011120, respectively), began to turn yellow 24 h after being injected with Agrobacterium tumefaciens (Figure 8), but no change occurred in the leaves injected with H2O or blank vectors. These results indicated that the three TFs may play important roles in the leaf senescence process.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/666660/fpls-12-666660-HTML/image_m/fpls-12-666660-g008.jpg)
Figure 8. Functional validation of selected senescence-associated genes (SAGs) using an Agrobacterium-mediated transient expression assay. (A) The relative expression fold change of selected SAGs under dark (Dark), normal light (Control-light), and 150 mM NaCl (150 mM NaCl-light) conditions on day 2, 4, and 6 compared with day 0 by RNA-seq and qRT-PCR. (B) Representative Nicotiana benthamiana leaf senescence symptoms 24 h after infiltration with different constructs encoding MTR_5g011120 (SGR), MTR_4g081870 (5), MTR_6g012670 (7), MTR_3g088110 (8), and empty vector with GUS, respectively.
We predicted that the TF MTR_3g088110 co-regulates 56, 69, and 94 target genes on day 2, 4, and 6, respectively (Supplementary Figure 8). Moreover, these genes were associated with the response to oxygen-containing compounds, carbohydrate biosynthetic process, vacuole organization, sphingolipid metabolism, and plant hormone signal transduction.
Data Integration Analysis
To gain a comprehensive understanding of the leaf senescence process, we integrated transcriptomic data with physiological and biochemical data for a combined analysis. We found that the expression levels of the differentially expressed SAGs were closely related to plant metabolite activity. Many genes related to lipid oxidation and antioxidant defense were differentially expressed, which is consistent with the increased H2O2 contents observed. We also found differentially expressed SAGs involved in plant hormone signal transduction along with an increase in ABA content at the start of the stress treatment and then a decrease as leaf senescence proceeded (Supplementary Figure 9A, Supplementary Table 14). In addition, during leaf senescence, the genes related to nutrient cycling were abundantly expressed, indicating the nutrients' movement and redistribution. To elucidate the relationship between leaf senescence and leaf nutrients, we focused on the expression of SAGs related to protein and amino acid metabolism. We found 247 protein-metabolism-related genes among the 1546 SAGs, out of which 135 were related to amino acid metabolism, and most of them were down-regulated. We performed a STEM analysis of these genes and found that compared with the control-light group, most of them in the salt-induced group were down-regulated across the three time points, suggesting protein metabolism increased as leaves senesced (Supplementary Figure 9B).
Discussion
In this study, we examined salt stress-induced senescence in Medicago truncatula leaves using a time-course analysis of the expression patterns of SAGs. To obtain accurate SAG expression data during the complex process of senescence, we sampled highly controlled repeatable detached leaves and measured them at many time points. As a positive control for our data analysis, we also included a dark-induced senescence treatment. This experimental design enabled the superimposition of the data from both treatments and helped us to obtain time-course differential expression network data for those genes specifically associated with salt-induced leaf senescence.
Our data showed a significant downward trend in chlorophyll content during leaf senescence from 0 to 6 days under salt stress, which corresponds to the observed leaf senescence symptoms. We also found that 252 SAGs are related to chlorophyll and chloroplast metabolism (Supplementary Figure 9). Over a long period of time, many genes encoding proteins involved in the chlorophyll catabolism pathway have been discovered and cloned, including STAY-GREEN (Zhou et al., 2011a), SAG12 (Noh and Amasino, 1999), and WRKY53 (Zentgraf et al., 2010). Similarly, many enzymes related to chlorophyll degradation have been identified, such as chlorophyllase (AtCLH), pheophorbide a oxygenase (AtPaO), and red chlorophyll catabolite reductase (RCCR) (Mach et al., 2001; Benedetti and Arruda, 2002; PruŽinská et al., 2003).
H2O2 is commonly used as an indicator of salt stress in plants; thus, we examined the levels of this substances in our assay (Figure 2B). In our experiment, the H2O2 content accumulated throughout leaf senescence, although with some variation. H2O2 is a hub for ROS interconversion, which can directly or indirectly oxidize biomolecules, such as intracellular nucleic acids and proteins, and damage cell membranes, thereby accelerating cellular senescence and disintegration (Bieker et al., 2012; Chen et al., 2012). The structural integrity and physiological functions of the cell membrane were disrupted when the leaf cells were subjected to oxidative damage (Halliwell, 1984).
Leaf senescence is associated with many kinds of plant hormones, and ABA is an important stress hormone involved in leaf senescence (Hung and Kao, 2004). We assessed ABA levels throughout the salt-induced leaf senescence process and found that ABA increased during leaf senescence (Figure 2C). ABA can promote elevated protein hydrolase and RNase activities and accelerate protein and nucleic acid degradation, thereby reducing photosynthetic rates and promoting leaf senescence (Marschner, 2011). Liang et al. (2014) identified OsNAP, an essential component of the ABA-mediated plant senescence signaling pathway in rice, and elucidated the molecular regulatory mechanism of ABA-mediated senescence in rice leaves for the first time (Liang et al., 2014). Similarly, Mao et al. (2017) revealed the regulatory network underlying ABA synthesis using the ABA-NAC-SAGs model, providing a new theoretical basis for the study of the molecular mechanisms of ABA regulation of leaf senescence (Mao et al., 2017).
During plant senescence, gene expression levels are closely related to metabolite content, with many metabolites involved in amino acid and lipid metabolism increasing and those involved in chloroplast metabolism decreasing (Breeze et al., 2011; Moschen et al., 2016; Sekhon et al., 2019). In our assay, we identified 3,001, 3,787, and 4,419 SAGs in total on day 2, 4, and 6, respectively, through transcriptomic analysis of salt-treated detached leaves, which appear to be involved in multiple signaling pathways and regulatory networks in leaf senescence. GO and KEGG enrichment analyses indicated that differentially expressed genes were mainly enriched in lipid or carbohydrate metabolism, amino acid metabolism and synthesis of secondary metabolites. In addition, genes involved in nutrient recycling and photosynthesis pathways were significantly enriched. Changes in the expression of these genes may impact the final forage quality in agricultural production.
TFs have been reported to play central roles in leaf senescence (Balazadeh et al., 2008). In recent years, many senescence-associated TFs have been discovered in model plants. For example, Guo et al. (2004) identified at least 96 TFs with 3-fold changes in gene expression during leaf senescence (Guo et al., 2004). In our study, we predicted TFs from among significantly differentially expressed SAGs and found that many of these predicted TFs belonged to the bHLH, MYB, ERF, and NAC families, which are the major families known to be involved in regulating senescence. We selected nine TFs with significant expression differences, including members of the NAC, MYB, WRKY, and NAP families, and transiently expressed them in Nicotiana benthamiana to verify their functions. Three of the selected TFs caused senescence phenotypes in Nicotiana benthamiana leaves: MTR_4g081870, MTR_6g012670, MTR_3g088110. We found the high sequence identities of MTR_4g081870, MTR_6g012670, and MTR_3g088110 from Medicago truncatula with AT1G69490 (NAP), AT5G39610 (ANAC092), and AT1G01720 (NAC002) from Arabidopsis thaliana, respectively. Therefore, MTR_4g081870, MTR_6g012670, and MTR_3g088110 may be involved in the regulation of leaf senescence. Further bioinformatic analysis and laboratory experiments may help to characterize their functions in more details.
Members of the NAC TF family in Arabidopsis thaliana play an important role in natural senescence or the biotic and abiotic stress-induced leaf senescence process. NAC1 is regulated by miR164, which regulates leaf senescence and causes cell death. However, Arabidopsis thaliana mutants that lack NAC1 did not show significant leaf senescence, and thus, NAC1 may not play a key role in leaf senescence (Guo et al., 2005; D'haeseleer et al., 2011). AT1G69490 (AtNAP) has been confirmed to be directly related to leaf senescence (Guo and Gan, 2006). AT5G39610 (ANAC092) plays a major role in the leaf senescence-related gene regulatory network and is an important network element that regulates salt stress and promotes senescence (Balazadeh et al., 2010). Garapati et al. (2015) suggested that AT1G01720 (NAC002) executes its physiological role by affecting both key chloroplast maintenance and senescence-promoting TFs (Garapati et al., 2015). Therefore, NAC TF family may also play a critical role in salt induced leaf senescence in Medicago truncatula, although the details of molecular mechanism is not well-known.
We also found that MTR_2g091190 (TF bHLH30), MTR_1g052470 (TF bHLH122-like protein), and MTR_1g080890 (TF UNE12) were highly up-regulated during leaf senescence, but these genes have not been reported to be directly related to plant senescence until now. MTR_3g103580, which is a member of the MYB family, showed high expression levels during leaf senescence in Medicago truncatula. This TF has high sequence similarity to AT1G10070 (BCAT2) from Arabidopsis thaliana. BCAT2 encodes a chloroplast branched-chain amino acid aminotransferase and is involved in cell wall development (Angelovici et al., 2013). Therefore, we have reason to believe that this TF of MTR_3g103580 regulates the metabolic activity of chloroplasts.
Many studies have shown that ethylene is involved in leaf abscission and may thus be associated with the initiation of leaf senescence (Aharoni and Lieberman, 1979). We also identified one ERF TF that showed high expression levels. MTR_1g074370 (ethylene-responsive TF 1B) had high sequence similarity to AT3G23240 (ERF1B) from Arabidopsis thaliana. ERF1 regulates ethylene overproduction and is involved in the ethylene signaling cascade. ERF1 was highly induced by high salinity stress in Arabidopsis and plays a positive role in salt stress tolerance (Lorenzo et al., 2003; Cheng et al., 2013, 2017). More laboratory experiments may help us understand their functions in more details.
Plant senescence is an important part of plant growth and development, and the timing of crop senescence profoundly affects the stability of crop yield, quality, and regional distribution. Breeders consider appropriate senescence time to be a key agronomic trait and thus aim to obtain superior germplasms with ideal senescence times in order to improve agricultural production. Alfalfa (Medicago sativa L.) is one of the world's most important forage crops owing to its high yield, high nutritional quality and wide adaptability. Additionally, most of the nutrients in alfalfa are stored in its leaves, and leaf senescence may therefore greatly affect the overall nutritional quality of alfalfa, especially when induced by external factors, such as salinity, drought, or cold stress, which can lead to premature aging of leaves. Therefore, preventing premature aging or senescence to increase biomass accumulation is profoundly important for improving alfalfa production. Furthermore, our previous research has shown that silencing the endogenous SGR gene (MsSGR) in alfalfa can result in the retention of more than 50% of their chlorophyll content during senescence, thus significantly increasing the crude protein contents compared to control lines (Zhou et al., 2011a).
In summary, M. truncatula, a close relative of alfalfa (M. sativa), has been widely used as a model in legume biology and is therefore a valuable experimental system for studying forage senescence. We have presented and analyzed a highly informative gene expression dataset over the course of senescence under salt stress and darkness. Our knowledge about the mechanisms regulating leaf senescence can be applied to improving the biomass and quality of forage under salt stress through future forage breeding.
Data Availability Statement
The original contributions presented in the study are publicly available. This data can be found at: NCBI repository, accession number: PRJNA701852 (https://www.ncbi.nlm.nih.gov/bioproject/701852).
Author Contributions
Z-YW, MC, and SD designed and planned the experiments. SD and LS performed the experiments. SD analyzed the data and wrote the manuscript. MC, HX, and Z-YW revised the manuscript. All authors read, revised, and approved the final manuscript.
Funding
This research was supported by the National Natural Science Foundation of China (U1906201) and the First Class Grassland Science Discipline Program of Shandong Province, China.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We are grateful to our colleagues in our laboratory for their advice in the execution of this study, and we also thank Peng Yu for help with sample collection.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2021.666660/full#supplementary-material
Supplementary Figure 1. The sampled leaf tissues and the relationship between leaf senescence and different NaCl concentration treatments. (A) The leaves detached for treatment-induced senescence. (B) The relationship between leaf senescence and different concentrations of NaCl.
Supplementary Figure 2. The validation of RNA-seq data via qRT-PCR. The relative expression fold change of the representative senescence-associated genes (SAGs) from different clusters in Figure 5 in different conditions compared with day 0.
Supplementary Figure 3. Significant Gene Ontology (GO) distribution for molecular function of 1546 senescence-associated genes (SAGs) in Figure 4D.
Supplementary Figure 4. Gene Ontology (GO) enrichment analysis of senescence-associated genes (SAGs) shared by darkness and salt-induced leaf senescence on day 2 (A for Figure 4A, 3011 enrichment analysis of genes enrichment analysis), 4 (B for Figure 4B, 3787 enrichment analysis of genes enrichment analysis), and 6 (C for Figure 4C, enrichment analysis of 4419 genes) treatment in detached Medicago truncatula leaves. Only the GO term corresponding with the top 20 enrichment levels are shown. GO term are shown along the vertical axis. In contrast, the horizontal axis represents the enrichment factor, i.e., the ratio of the number of genes/transcripts (sample number) enriched in this pathway relative to the number of annotated genes/transcripts (background number). The greater the enrichment factor, the greater the degree of enrichment. The size of the point indicates the number of genes in this pathway, while the colors of the points correspond to different P-value ranges.
Supplementary Figure 5. Scatter diagram based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of the senescence-associated genes (SAGs) shared by darkness and salt-induced leaf senescence on day 2 (A for 3,011 genes), 4 (B for 3,787 genes), and 6 (C for 4,419 genes) in detached Medicago truncatula leaves. Only the KEGG pathways within the top 20 enrichment levels are displayed. Pathway names are shown along the vertical axis, while the horizontal axis represents the enrichment factor, i.e., the ratio of the number of genes/transcripts (sample number) enriched in this pathway relative to the number of annotated genes/transcripts (background number). The greater the enrichment factor, the greater the degree of enrichment. The size of the point indicates the number of genes in this pathway, while the colors of the points correspond to different P-adjust value ranges.
Supplementary Figure 6. Scatter diagram based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of the senescence-associated genes (SAGs) in Figure 5 for up-regulated (A) and down-regulated genes (B).
Supplementary Figure 7. The maSigPro analysis of Differentially expressed senescence-associated genes (SAGs) in all samples. The horizontal axis represents specific samples, the vertical axis represents the expression values after treatment log10(x+1), and the different colored lines represent the trend of different genes in all samples.
Supplementary Figure 8. A hypothetical network of transcription factor MTR_3g088110. It showed the predicted potential target genes of transcription factors MTR_3g088110 across three time points: MTR_3g088110-2, MTR_3g088110-4, and MTR_3g088110-6 display the genes from day 2 Figure 4A 3011genes, day 4 Figure 4B 3787genes, and day 6 Figure 4C 4419genes co-regulated with MTR_3g088110, respectively.
Supplementary Figure 9. Differentially expressed senescence-associated genes (SAGs) related to physiological and biochemical indicators. (A) The number of genes from darkness and salt treatment enriched 1546 SAGs was associated with physiological and biochemical changes during salt induced leaf senescence. (B) Short Time-series Expression Miner (STEM) analysis of the genes in (A).
References
Aharoni, N., and Lieberman, M. (1979). Ethylene as a regulator of senescence in tobacco leaf discs. Plant Physiol. 64, 801–804. doi: 10.1104/pp.64.5.801
Allu, A. D., Soja, A. M., Wu, A., Szymanski, J., and Balazadeh, S. (2014). Salt stress and senescence: identification of cross-talk regulatory components. J. Exp. Bot. 65, 3993–4008. doi: 10.1093/jxb/eru173
Angelovici, R., Lipka, A. E., Deason, N., Gonzalez-Jorge, S., Lin, H., Cepela, J., et al. (2013). Genome-wide analysis of branched-chain amino acid levels in Arabidopsis seeds. Plant Cell. 25, 4827–4843. doi: 10.1105/tpc.113.119370
Balazadeh, S., Riaño-Pachón, D., and Mueller-Roeber, B. (2008). Transcription factors regulating leaf senescence in Arabidopsis thaliana. Plant Biol. 10, 63–75. doi: 10.1111/j.1438-8677.2008.00088.x
Balazadeh, S., Wu, A., and Mueller-Roeber, B. (2010). Salt-triggered expression of the ANAC092-dependent senescence regulon in Arabidopsis thaliana. Plant Signal. Behav. 5, 733–735. doi: 10.4161/psb.5.6.11694
Bar-Dror, T., Dermastia, M., Kladnik, A., Žnidarič, M. T., Novak, M. P., Meir, S., et al. (2011). Programmed cell death occurs asymmetrically during abscission in tomato. Plant Cell. 23, 4146–4163. doi: 10.1105/tpc.111.092494
Benedetti, C. E., and Arruda, P. (2002). Altering the expression of the chlorophyllase geneATHCOR1 in transgenic Arabidopsis caused changes in the chlorophyll-to-chlorophyllide ratio. Plant Physiol. 128, 1255–1263. doi: 10.1104/pp.010813
Bian, Z., Gao, H., and Wang, C. (2021). NAC transcription factors as positive or negative regulators during ongoing battle between pathogens and our food crops. Int. J. Mol. Sci. 22:81. doi: 10.3390/ijms22010081
Bieker, S., Riester, L., Stahl, M., Franzaring, J., and Zentgraf, U. (2012). Senescence-specific alteration of hydrogen peroxide levels in Arabidopsis thaliana and oilseed rape spring variety Brassica napus L. cv. Mozart. J. Integrat. Plant Biol. 54, 540–554. doi: 10.1111/j.1744-7909.2012.01147.x
Breeze, E. (2020). Make, modify, move: multilayered regulation of ONAC054 during ABA-induced leaf senescence. Plant Cell. 32, 532–533. doi: 10.1105/tpc.20.00003
Breeze, E., Harrison, E., Mchattie, S., Hughes, L., Hickman, R., Hill, C., et al. (2011). High-resolution temporal profiling of transcripts during Arabidopsis leaf senescence reveals a distinct chronology of processes and regulation. Plant Cell. 23, 873–894. doi: 10.1105/tpc.111.083345
Buchanan-Wollaston, V. (1997). The molecular biology of leaf senescence. J. Exp. Bot. 48, 181–199. doi: 10.1093/jxb/48.2.181
Buchanan-Wollaston, V., Earl, S., Harrison, E., Mathas, E., Navabpour, S., Page, T., et al. (2003). The molecular analysis of leaf senescence–a genomics approach. Plant Biotechnol. J. 1, 3–22. doi: 10.1046/j.1467-7652.2003.00004.x
Buchanan-Wollaston, V., Page, T., Harrison, E., Breeze, E., Lim, P. O., Nam, H. G., et al. (2005). Comparative transcriptome analysis reveals significant differences in gene expression and signalling pathways between developmental and dark/starvation-induced senescence in Arabidopsis. Plant J. 42, 567–585. doi: 10.1111/j.1365-313X.2005.02399.x
Chen, H.-J., Lin, Z.-W., Huang, G.-J., and Lin, Y.-H. (2012). Sweet potato calmodulin SPCAM is involved in salt stress-mediated leaf senescence, H2O2 elevation and senescence-associated gene expression. J. Plant Physiol. 169, 1892–1902. doi: 10.1016/j.jplph.2012.08.004
Cheng, M.-C., Liao, P.-M., Kuo, W.-W., and Lin, T.-P. (2013). The Arabidopsis ethylene response factor1 regulates abiotic stress-responsive gene expression by binding to different cis-acting elements in response to different stress signals. Plant Physiol. 162, 1566–1582. doi: 10.1104/pp.113.221911
Cheng, M. C., Kuo, W. C., Wang, Y. M., Chen, H. Y., and Lin, T. P. (2017). UBC 18 mediates ERF 1 degradation under light-dark cycles. New Phytologist. 213, 1156–1167. doi: 10.1111/nph.14272
D'haeseleer, K., Den Herder, G., Laffont, C., Plet, J., Mortier, V., Lelandais-Brière, C., et al. (2011). Transcriptional and post-transcriptional regulation of a NAC1 transcription factor in Medicago truncatula roots. New Phytologist. 191, 647–661. doi: 10.1111/j.1469-8137.2011.03719.x
Du, J., Lu, S., Chai, M., Zhou, C., Sun, L., Tang, Y., et al. (2020). Functional characterization of Petiolule-Like Pulvinus (PLP) gene in abscission zone development in Medicago truncatula and its application to genetic improvement of alfalfa. Plant Biotechnol. J. 19, 351–364. doi: 10.1111/pbi.13469
Garapati, P., Xue, G.-P., Munné-Bosch, S., and Balazadeh, S. (2015). Transcription factor ATAF1 in Arabidopsis promotes senescence by direct regulation of key chloroplast maintenance and senescence transcriptional cascades. Plant Physiol. 168, 1122–1139. doi: 10.1104/pp.15.00567
Gregersen, P. L., Culetic, A., Boschian, L., and Krupinska, K. (2013). Plant senescence and crop productivity. Plant Mol. Biol. 82, 603–622. doi: 10.1007/s11103-013-0013-8
Guo, H.-S., Xie, Q., Fei, J.-F., and Chua, N.-H. (2005). MicroRNA directs mRNA cleavage of the transcription factor NAC1 to downregulate auxin signals for Arabidopsis lateral root development. Plant Cell. 17, 1376–1386. doi: 10.1105/tpc.105.030841
Guo, Y., Cai, Z., and Gan, S. (2004). Transcriptome of Arabidopsis leaf senescence. Plant Cell Environ. 27, 521–549. doi: 10.1111/j.1365-3040.2003.01158.x
Guo, Y., and Gan, S. (2006). AtNAP, a NAC family transcription factor, has an important role in leaf senescence. Plant J. 46, 601–612. doi: 10.1111/j.1365-313X.2006.02723.x
Halliwell, B. (1984). Chloroplast metabolism: the structure and function of chloroplasts in green leaf cells. Oxford University Press. 184:259. doi: 10.1016/0014-5793(85)80685-2
He, L., Wu, W., Zinta, G., Yang, L., Wang, D., Liu, R., et al. (2018). A naturally occurring epiallele associates with leaf senescence and local climate adaptation in Arabidopsis accessions. Nat. Commun. 9, 1–11. doi: 10.1038/s41467-018-02839-3
Hörtensteiner, S., and Feller, U. (2002). Nitrogen metabolism and remobilization during senescence. J. Exp. Bot. 53, 927–937. doi: 10.1093/jexbot/53.370.927
Hung, K. T., and Kao, C. H. (2004). Hydrogen peroxide is necessary for abscisic acid-induced senescence of rice leaves. J. Plant Physiol. 161, 1347–1357. doi: 10.1016/j.jplph.2004.05.011
Jaradat, M. R., Feurtado, J. A., Huang, D., Lu, Y., and Cutler, A. J. (2013). Multiple roles of the transcription factor AtMYBR1/AtMYB44 in ABA signaling, stress responses, and leaf senescence. BMC Plant Biol. 13, 1–19. doi: 10.1186/1471-2229-13-192
Kim, Y.-S., Sakuraba, Y., Han, S.-H., Yoo, S.-C., and Paek, N.-C. (2013). Mutation of the Arabidopsis NAC016 transcription factor delays leaf senescence. Plant Cell Physiol. 54, 1660–1672. doi: 10.1093/pcp/pct113
Kong, X., Luo, Z., Dong, H., Eneji, A. E., Li, W., and Lu, H. (2013). Gene expression profiles deciphering leaf senescence variation between early-and late-senescence cotton lines. PLoS ONE 8:e69847. doi: 10.1371/journal.pone.0069847
Lee, R. H., Wang, C. H., Huang, L. T., and Chen, S. C. G. (2001). Leaf senescence in rice plants: cloning and characterization of senescence up-regulated genes. J. Exp. Bot. 52, 1117–1121. doi: 10.1093/jexbot/52.358.1117
Li, C., Liu, C.-Q., Zhang, H.-S., Chen, C.-P., Yang, X.-R., Chen, L.-F., et al. (2021). LPS1, encoding iron-sulfur subunit SDH2-1 of succinate dehydrogenase, affects leaf senescence and grain yield in rice. Int. J. Mol. Sci. 22:157. doi: 10.3390/ijms22010157
Li, X., Chanroj, S., Wu, Z., Romanowsky, S. M., Harper, J. F., and Sze, H. (2008). A distinct endosomal Ca2+/Mn2+ pump affects root growth through the secretory process. Plant Physiol. 147, 1675–1689. doi: 10.1104/pp.108.119909
Liang, C., Wang, Y., Zhu, Y., Tang, J., Hu, B., Liu, L., et al. (2014). OsNAP connects abscisic acid and leaf senescence by fine-tuning abscisic acid biosynthesis and directly targeting senescence-associated genes in rice. Proc. Nat. Acad. Sci. U.S.A. 111, 10013–10018. doi: 10.1073/pnas.1321568111
Lim, P. O., Kim, H. J., and Gil Nam, H. (2007). Leaf senescence. Annu. Rev. Plant Biol. 58, 115–136. doi: 10.1146/annurev.arplant.57.032905.105316
Lim, P. O., Woo, H. R., and Nam, H. G. (2003). Molecular genetics of leaf senescence in Arabidopsis. Trends Plant Sci. 8, 272–278. doi: 10.1016/S1360-1385(03)00103-1
Lin, J. F., and Wu, S. H. (2004). Molecular events in senescing Arabidopsis leaves. Plant J. 39, 612–628. doi: 10.1111/j.1365-313X.2004.02160.x
Liu, J. D., Goodspeed, D., Sheng, Z., Li, B., Yang, Y., Kliebenstein, D. J., et al. (2015). Keeping the rhythm: light/dark cycles during postharvest storage preserve the tissue integrity and nutritional content of leafy plants. BMC Plant Biol. 15, 1–9. doi: 10.1186/s12870-015-0474-9
Lorenzo, O., Piqueras, R., Sánchez-Serrano, J. J., and Solano, R. (2003). Ethylene Response Factor1 integrates signals from ethylene and jasmonate pathways in plant defense. Plant Cell. 15, 165–178. doi: 10.1105/tpc.007468
Lutts, S., Kinet, J., and Bouharmont, J. (1996). NaCl-induced senescence in leaves of rice (Oryza sativa L.) cultivars differing in salinity resistance. Ann. Bot. 78, 389–398. doi: 10.1006/anbo.1996.0134
Mach, J. M., Castillo, A. R., Hoogstraten, R., and Greenberg, J. T. (2001). The Arabidopsis-accelerated cell death gene ACD2 encodes red chlorophyll catabolite reductase and suppresses the spread of disease symptoms. Proc. Nat. Acad. Sci. U.S.A. 98, 771–776. doi: 10.1073/pnas.98.2.771
Mao, C., Lu, S., Lv, B., Zhang, B., Shen, J., He, J., et al. (2017). A rice NAC transcription factor promotes leaf senescence via ABA biosynthesis. Plant Physiol. 174, 1747–1763. doi: 10.1104/pp.17.00542
Marschner, H. (2011). Marschner's Mineral Nutrition of Higher Plants. Academic press. doi: 10.1016/C2009-0-63043-9
Moschen, S., Bengoa Luoni, S., Di Rienzo, J. A., Caro, M. D. P., Tohge, T., Watanabe, M., et al. (2016). Integrating transcriptomic and metabolomic analysis to understand natural leaf senescence in sunflower. Plant Biotechnol. J. 14, 719–734. doi: 10.1111/pbi.12422
Munné-Bosch, S., and Alegre, L. (2004). Die and let live: leaf senescence contributes to plant survival under drought stress. Functional Plant Biol. 31, 203–216. doi: 10.1071/FP03236
Munns, R., and Tester, M. (2008). Mechanisms of salinity tolerance. Annu. Rev. Plant Biol. 59, 651–681. doi: 10.1146/annurev.arplant.59.032607.092911
Noh, Y.-S., and Amasino, R. M. (1999). Identification of a promoter region responsible for the senescence-specific expression of SAG12. Plant Mol. Biol. 41, 181–194. doi: 10.1023/A:1006342412688
Pageau, K., Reisdorf-Cren, M., Morot-Gaudry, J.-F., and Masclaux-Daubresse, C. (2006). The two senescence-related markers, GS1 (cytosolic glutamine synthetase) and GDH (glutamate dehydrogenase), involved in nitrogen mobilization, are differentially regulated during pathogen attack and by stress hormones and reactive oxygen species in Nicotiana tabacum L. leaves. J. Experi. Botany 57, 547–557. doi: 10.1093/jxb/erj035
Park, S.-Y., Yu, J.-W., Park, J.-S., Li, J., Yoo, S.-C., Lee, N.-Y., et al. (2007). The senescence-induced staygreen protein regulates chlorophyll degradation. Plant Cell. 19, 1649–1664. doi: 10.1105/tpc.106.044891
Parlitz, S., Kunze, R., Mueller-Roeber, B., and Balazadeh, S. (2011). Regulation of photosynthesis and transcription factor expression by leaf shading and re-illumination in Arabidopsis thaliana leaves. J. Plant Physiol. 168, 1311–1319. doi: 10.1016/j.jplph.2011.02.001
Peleg, Z., and Blumwald, E. (2011). Hormone balance and abiotic stress tolerance in crop plants. Curr. Opin. Plant Biol. 14, 290–295. doi: 10.1016/j.pbi.2011.02.001
PruŽinská, A., Tanner, G., Anders, I., Roca, M., and Hörtensteiner, S. (2003). Chlorophyll breakdown: pheophorbide a oxygenase is a Rieske-type iron-sulfur protein, encoded by the accelerated cell death 1 gene. Proc. Nat. Acad. Sci. U.S.A. 100, 15259–15264. doi: 10.1073/pnas.2036571100
Rushton, D. L., Tripathi, P., Rabara, R. C., Lin, J., Ringler, P., Boken, A. K., et al. (2012). WRKY transcription factors: key components in abscisic acid signalling. Plant Biotechnol. J. 10, 2–11. doi: 10.1111/j.1467-7652.2011.00634.x
Sakuraba, Y., Kim, D., and Paek, N.-C. (2018). Salt treatments and induction of senescence. Plant Senescence. 1744, 141–149. doi: 10.1007/978-1-4939-7672-0_13
Schippers, J. H., Schmidt, R., Wagstaff, C., and Jing, H.-C. (2015). Living to die and dying to live: the survival strategy behind leaf senescence. Plant Physiol. 169, 914–930. doi: 10.1104/pp.15.00498
Sekhon, R. S., Saski, C., Kumar, R., Flinn, B. S., Luo, F., Beissinger, T. M., et al. (2019). Integrated genome-scale analysis identifies novel genes and networks underlying senescence in maize. Plant Cell. 31, 1968–1989. doi: 10.1105/tpc.18.00930
Sobieszczuk-Nowicka, E., Wrzesiński, T., Bagniewska-Zadworna, A., Kubala, S., Rucińska-Sobkowiak, R., Polcyn, W., et al. (2018). Physio-genetic dissection of dark-induced leaf senescence and timing its reversal in barley. Plant Physiol. 178, 654–671. doi: 10.1104/pp.18.00516
Uauy, C., Distelfeld, A., Fahima, T., Blechl, A., and Dubcovsky, J. (2006). A NAC gene regulating senescence improves grain protein, zinc, and iron content in wheat. Science 314, 1298–1301. doi: 10.1126/science.1133649
Wingler, A., Purdy, S., Maclean, J. A., and Pourtau, N. (2006). The role of sugars in integrating environmental signals during the regulation of leaf senescence. J. Exp. Bot. 57, 391–399. doi: 10.1093/jxb/eri279
Woo, H. R., Chung, K. M., Park, J.-H., Oh, S. A., Ahn, T., Hong, S. H., et al. (2001). ORE9, an F-box protein that regulates leaf senescence in Arabidopsis. Plant Cell. 13, 1779–1790. doi: 10.1105/TPC.010061
Woo, H. R., Kim, H. J., Lim, P. O., and Nam, H. G. (2019). Leaf senescence: systems and dynamics aspects. Annu. Rev. Plant Biol. 70, 347–376. doi: 10.1146/annurev-arplant-050718-095859
Xu, P., Chen, H., and Cai, W. (2020). Transcription factor CDF4 promotes leaf senescence and floral organ abscission by regulating abscisic acid and reactive oxygen species pathways in Arabidopsis. EMBO Rep. 21:e48967. doi: 10.15252/embr.201948967
Yamada, T., Ichimura, K., Kanekatsu, M., and Van Doorn, W. G. (2009). Homologs of genes associated with programmed cell death in animal cells are differentially expressed during senescence of Ipomoea nil petals. Plant Cell Physiol. 50, 610–625. doi: 10.1093/pcp/pcp019
Yoshimoto, K., Jikumaru, Y., Kamiya, Y., Kusano, M., Consonni, C., Panstruga, R., et al. (2009). Autophagy negatively regulates cell death by controlling NPR1-dependent salicylic acid signaling during senescence and the innate immune response in Arabidopsis. Plant Cell. 21, 2914–2927. doi: 10.1105/tpc.109.068635
Zentgraf, U., Laun, T., and Miao, Y. (2010). The complex regulation of WRKY53 during leaf senescence of Arabidopsis thaliana. Eur. J. Cell Biol. 89, 133–137. doi: 10.1016/j.ejcb.2009.10.014
Zhang, W. Y., Xu, Y. C., Li, W. L., Yang, L., Yue, X., Zhang, X. S., et al. (2014). Transcriptional analyses of natural leaf senescence in maize. PLoS ONE 9:e115617. doi: 10.1371/journal.pone.0115617
Zhang, Z., and Guo, Y. (2018). Hormone treatments in studying leaf senescence. Plant Senescence. 1744, 125–132. doi: 10.1007/978-1-4939-7672-0_11
Zhou, C., Han, L., Pislariu, C., Nakashima, J., Fu, C., Jiang, Q., et al. (2011a). From model to crop: functional analysis of a STAY-GREEN gene in the model legume Medicago truncatula and effective use of the gene for alfalfa improvement. Plant Physiol. 157, 1483–1496. doi: 10.1104/pp.111.185140
Keywords: salt stress, leaf senescence, Medicago truncatula, comparative transcriptome analysis, RNA-Seq
Citation: Dong S, Sang L, Xie H, Chai M and Wang Z-Y (2021) Comparative Transcriptome Analysis of Salt Stress-Induced Leaf Senescence in Medicago truncatula. Front. Plant Sci. 12:666660. doi: 10.3389/fpls.2021.666660
Received: 10 February 2021; Accepted: 14 June 2021;
Published: 09 July 2021.
Edited by:
László Szabados, Hungarian Academy of Sciences (MTA), HungaryCopyright © 2021 Dong, Sang, Xie, Chai and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Maofeng Chai, mfchai@qau.edu.cn orcid.org/0000-0001-9915-0321