- College of Grassland Science, Nanjing Agricultural University, Nanjing, China
Salinity is a major limiting factor that affects crop production. Understanding of the mechanisms of plant salt tolerance is critical for improving crop yield on saline land. Alfalfa (Medicago sativa L.) is the most important forage crop, while its salt tolerance mechanisms are largely unknown. The physiological and transcriptomic responses in two contrasting salt tolerant cultivars to salinity stress were investigated in the present study. “Magnum Salt” showed higher salt tolerance than “Adrenalin,” with higher relative germination rate, survival rate, biomass and K+/Na+ ratio after salt treatment. Activities of antioxidant enzymes SOD, CAT and GR, and proline concentrations were upregulated to higher levels in roots and shoots in Magnum Salt than in Adrenalin after salinity stress, except for no difference in GR activity in shoots, and lower levels of O2⋅– and H2O2 were accumulated in leaves. It was interesting to find that salinity caused a decrease in total unsaturated fatty acid in Adrenalin other than Magnum Salt, C18:2 was increased significantly after salinity in Magnum Salt, while it was unaltered in Adrenalin. High quality RNA sequencing (RNA-seq) data was obtained from samples of Magnum Salt and Adrenalin at different time points (0, 2, and 26 h). Generally, “phagosome,” “TCA cycle” and “oxidative phosphorylation” pathways were inhibited by salinity stress. Upregulated DEGs in Magnum Salt were specifically enriched in “fatty acid metabolism,” “MAPK signaling” and “hormone signal transduction” pathways. The DEGs involved in ionic homeostasis, reactive oxygen species (ROS) scavenging and fatty acid metabolism could partially explain the difference in salt tolerance between two cultivars. It is suggested that salt tolerance in alfalfa is associated with regulation of ionic homeostasis, antioxidative enzymes and fatty acid metabolism at both transcriptional and physiological level.
Introduction
Salinity is one of the most important abiotic stresses affecting agricultural productions. Over 800 million hectares of arable lands are facing the problem of salinity in the world (Liu and Wang, 2021). Salt stress can cause growth inhibition even senescence and death in plants. Hyperosmosis is the primary transient effect of salinization to plants that lead to water absorption inhibition and dehydration. As exposure prolong, ionic equilibrium is interrupted with Na+ hyperaccumulation and K+ and Ca2+ influx hinder, followed by secondary injury such as photosynthesis inhibition, oxidative injury and metabolic disorder. Cultivation of salt-tolerant plants is a sustainable strategy to make better use of saline land, contributing to carbon neutral at mean time. The identification of salt tolerant plant germplasm resources and understanding their salt tolerance mechanisms are significant for new salt tolerant plant breeding and saline soil utilization.
Alfalfa (Medicago sativa L.) is the world’s most extensively cultivated forage legume. Improving salt tolerance is an important goal in alfalfa breeding. Different genotypes of alfalfa exhibit wide variations in salt tolerance. Yu et al. (2021) evaluated salt tolerance of 20 alfalfa cultivars at seedling stage and divided them into highly salt-tolerant, salt-tolerant, moderately salt-tolerant and salt-sensitive groups. Cornacchione and Suarez (2017) identified “SISA14” and “SW 8421S” populations as the most salt tolerant out of 15 alfalfa populations with high production in salt condition. In our research, salt-tolerance of 8 alfalfa cultivars was studied and one salt tolerant cultivar “Magnum Salt” and one salt sensitive cultivar “Adrenalin” were identified for further study.
Plants utilize multiple strategies to alleviate salt stress damage. Thus, salt tolerance is a complex trait controlled by multiple genes involving in salt perception, signaling transduction, osmotic regulation, ion transport, hormones synthesis, photosynthesis and metabolization (Deinlein et al., 2014). A few genes involved in salt resistance have been reported in alfalfa. For example, MsCBL4 enhances calcium metabolism but not sodium transport in transgenic tobacco under salt stress (An et al., 2020); MsPP2CA1 gene confers enhanced salt tolerance in Arabidopsis, which may due to osmolytes accumulation other than protecting the cell membrane from the damage by superoxide anion radicals (Guo et al., 2018). Besides, gene function of MsGRP (Long et al., 2013), MsWRKY11 (Wang et al., 2018), MsERF11 (Chen et al., 2012), MsGME (Ma et al., 2014), MsRCI2A/B/C (Li et al., 2021) and MsZEP (Zhang et al., 2016) in salt resistance have been reported.
RNA sequencing allows for analysis of the transcript levels as well as alternative splicing and novel transcript identification, showing great potential in tolerance mechanism studying and salt tolerant genes identification. Gruber et al. (2017) compared the transcriptome profiles of two salinity-tolerant alfalfa breeding populations and a more salinity-sensitive population, found that both saline-tolerant populations showed more substantial up-regulation in redox-related genes and B-ZIP transcripts. Time course full-length transcriptome database for alfalfa root tips under NaCl was constructed by Luo et al. (2019). Total of 8861 NaCl-regulated differentially expressed genes were identified, in which continuously up-regulated genes were mainly involved in “ion homeostasis,” “antiporter activity,” “trehalose biosynthetic process,” “thiamine pyrophosphate binding,” and “ethylene-activated signaling pathway.” By analyzing the transcriptome of salt tolerant “Halo” and salt intolerant “Vernal,” Bhattarai et al. (2021) found out that ion binding was the key molecular activity for salt tolerance in alfalfa. Furthermore, they identified 28 (leaf) and 31 (root) salt responsive candidate genes involved in transmembrane protein function, photosynthesis, carbohydrate metabolism, defense against oxidative damage, cell wall modification and protection against lipid peroxidation.
In our research, time course transcriptome database of roots from salt tolerant cultivar “Magnum Salt” and salt sensitive “Adrenalin” was conducted to understand the salt resistance mechanism of alfalfa. It was found that regulation of genes involved in K+, Na+, and Ca2+ homeostasis, reactive oxygen species (ROS) scavenging, as well as fatty acid metabolization, especially unsaturated fatty acid linoleic acid (C18:2) metabolization, might play important roles in salt resistance in alfalfa.
Materials and Methods
Evaluation of Germination Rate Under Salt Solution
A germination test was carried out using eight alfalfa cultivars, and 150 mM NaCl was used as salt treatment according to pre-experiments. Fifty plump seeds per petri dish were treated with 30 mL ddH2O or 150 mM NaCl solution with three replicates. The criterion of germination is the appearance of a 2 mm radicle protrusion. The relative germination was calculated as 100 × germination rate within 7 days in salt treatment group divided by germination rate in control group.
Evaluation of Survival Rate Under Salt Treatment
Survival rates were evaluated according to Sun et al. (2020) with modification. The four-week old seedings planted in 50-pots tray were treated with 1 L Hoagland nutrient solution with or without 500 mM NaCl (pH 5.8) weekly for 3 weeks. The surviving plants with turgid and green leaves and dead plants were counted for calculation of survival rate after 1 week of recovery by washing NaCl solution with water.
Evaluation of Biomass and K+/Na+/Ca2+ Concentration by Hydroponic Culture
Two-week old Magnum Salt and Adrenalin seedlings grown in soil with a mixture of pearlite, vermiculite and perlite (1:1:2) were transferred to a plate filled with 1/2 Hoagland nutrient solution (pH 5.8) for 1 week recovery growth. After that, seedlings were treated with 1/2 Hoagland nutrient solution containing 0 or 150 mM NaCl for 7 days, which was determined according to pre-experiments. Biomass were analyzed with 10 replicates, and K+/Na+/Ca2+ concentrations were measured as described by Dai et al. (2022).
ROS, Antioxidant Enzyme Activities and Proline Content Analysis
Four-week old Magnum Salt and Adrenalin seedlings grown in soil (pearlite: vermiculite: perlite 1:1:2) were treated with ddH2O or 500 mM NaCl every 3 days, which was modified with reference to method of Dai et al. (2022). Six days later, DAB (3,3′-Diaminobenzidine) and NBT (ρ-nitro blue tetrazolium chloride) staining for H2O2 and O2⋅– were performed as below. Leaves of the same position were collected followed by submerging into 0.1% DAB or NBT solution for 24 h. After decolorizing with 95% alcohol, photographs of the leaves were taken under an optical microscope. Antioxidant enzyme activity and free proline content were analyzed according to Zhuo et al. (2018).
Fatty Acids Content Analysis
Roots from seedlings treated with 500 mM NaCl solution for 3 and 6 days were collected for fatty acid content analysis. Gas chromatographic analysis was used and the procedure was as described previous (Huang et al., 2017). DBI (Double-bond index) was used to measure the degree of lipid total unsaturation. DBI = [1 × (% 18:1) + 2 × (% 18:2) + 3 × (% 18:3)]/100, as described by De Vos et al. (1993).
RNA Extraction, Library Construction, and Sequencing
Two-week old seedlings grown in soil with a mixture of pearlite, vermiculite and perlite (1:1:2) were transferred to a plate filled with 1/2 Hoagland nutrient solution (pH 5.8) for 1 week recovery growth. After that, seedlings were treated with 1/2 Hoagland nutrient solution containing 150 mM NaCl, roots were sampled at 0, 2, and 26 h. Total RNA of roots were extracted using a RNAprep pure Plant Kit (Tiangen Inc., Beijing, China). mRNAs were purified using Dynabeads Oligo(dT)25 (Life Technologies). Ribo-Zero™ Magnetic Kit (Epicentre) was used to remove rRNA. The derived mRNAs were fragmented and reverse transcribed into cDNAs with random hexamer. The cDNAs were purified and ligated to adaptors for Illumina paired-end sequencing. The cDNA library was sequenced using the Illumina HiSeq™4000 system by Gene Denovo Biotechnology Co. (Guangzhou, China). The sequencing data has all been archived in the NCBI Sequence Read Archive (SRA) database under accession number PRJNA830977.
Unigene de novo Assembly and Annotation
Raw reads were initially cleaned by removing adaptor sequences and low quality reads to get clean reads. The Trinity program was carried out to assemble the clean reads to obtain nonredundant unigenes. The assembly unigenes functional annotations were obtained using NCBI nucleotide sequence database (Nt), nonredundant protein database, Swiss-Prot database, Clusters of Orthologous Groups of proteins (COG) database, Gene Ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) annotation database.
Identification of Differentially Expressed Genes
Differentially expressed genes (DEGs) was analyzed by a DESeq R package1 with the thresholds of false discovery rate (FDR) < 0.05 and | Log2Fold-Change| > 4. Then DEGs GO enrichment and KEGG analysis were performed with agriGO 2.02 and KOBAS 3.03, respectively.
Differentially Expressed Genes Temporal Expression Patterns Analysis
To achieve the temporal expression patterns, the DEGs of at 2 and 26 h after salt treatment were analyzed with the Short Time-series Expression Miner (STEM) software (v1.3.11; Ernst and Bar-Joseph, 2006). The potential function annotations and enrichments of DEGs identified in STEM analysis were explored with KEGG analysis.
Heatmap Analysis
FPKM was used for representing the expression abundance of genes. The clustered heatmap was portrayed after normalized with zero-to-one scale method using the TBtools software (Chen et al., 2020).
Quantitative Real-Time PCR Analysis
cDNA synthesis was performed with HiScript III RT SuperMix for qPCR (+gDNA wiper) reagent kit with gDNA Eraser (Vazyme, China). qRT-PCR reaction was performed following the instructions of ChamQ Universal SYBR qPCR Master Mix (Vazyme, China). All primer sequences are shown in Supplementary Table 5. The relative gene expression level was calculated with 2–ΔΔCt method.
Statistical Analysis
All data were subjected to analysis of variances according to the model for completely randomized design using an SPSS program (SPSS Inc., Chicago, United States). Significant differences among means of cultivars in Figure 1 were evaluated by one-way ANOVA at 0.05 probability level. Significant differences among means of both treatments and cultivars were evaluated by LSD test at 0.05 probability level.
Figure 1. Evaluation of salt tolerance of 8 alfalfa cultivars by analyzing relative germination rate and survival rate under salt condition. Alfalfa seeds were treated with 0 or 150 mM NaCl solution, and relative germination rate (A) were calculated within 7 days. Four-week old seedlings treated with Hoagland nutrient solution with 500 mM NaCl (pH 5.8) weekly. Three weeks later, photographs (B) were taken and survival rates (C) were analyzed. Values are mean ± SE (n = 3 replicates). Different letters (P < 0.05, one-way ANOVA) compare different cultivars.
Results
Identification of Salt-Tolerant and Salt-Sensitive Cultivars
Salt tolerance of eight commonly planted cultivars was evaluated based on relative germination rate and survival rate in response to salt treatment. The relative germination rate decreased in “Magnum Salt,” “Magnum 601,” “Magnum II,” “Magnum 5010,” “Zhongmu No.1,” “Sanditi,” and “Adrenalin” sequentially (Figure 1A). “Magnum Salt,” “Magnum 601,” “Magnum II,” “Magnum 5010,” and “Zhongmu No.1” had significantly higher survival rate than “Sanditi” and “Adrenalin” after the seedling were irrigated with 500 mM NaCl solution for 3 weeks (Figures 1B,C). The data indicated that “Magnum Salt” and “Adrenalin” showed the highest and lowest salt tolerance, respectively.
Growth performance of Magnum Salt and Adrenalin was further evaluated by hydroponic culture experiment (Figure 2). Treatment with 150 mM NaCl caused obvious inhibition of root growth and yellowing of leaves in both cultivars (Figure 2A). Dry weight of roots or shoots showed no significant difference between Magnum Salt and Adrenalin under control condition (Figures 2B,C). It was decreased after salt treatment in both cultivars, while that was higher in Magnum Salt than that in Adrenalin (Figures 2B,C).
Figure 2. Effect of salt treatment on the biomass of Magnum Salt and Adrenalin. Three-week old seedlings were treated with 1/2 Hoagland solution containing 0 or 150 mM NaCl for 1 week. Photographs (A) were taken and biomass of shoot (B) and root (C) were analyzed. Values are mean ± SE (n = 10 replicates). Different letters indicate the significant differences (P < 0.05, LSD) between two cultivars or treatments.
Analysis of Na+, K+ and Ca2+ Concentrations
The Na+ concentration was increased in both roots and shoots after 7 d of salt treatment, with higher level in Adrenalin than in Magnum Salt (Figures 3A,B). K+ and Ca2+ concentrations were decreased in both cultivars, with higher levels in Magnum Salt than in Adrenalin after salt treatment (Figures 3C–F). As a consequence, the Na+/K+ ratio in roots and shoots were higher in Adrenalin than in Magnum Salt (Figures 3G,H).
Figure 3. Effect of salt treatment on the Na+, K+ and Ca2+ accumulation in Magnum Salt and Adrenalin. Three-week old seedlings were treated with 1/2 Hoagland solution containing 0 or 150 mM NaCl for 1 week. Na+, K+, Ca2+ content and Na+/K+ ratio in roots (A,C,E,G) and shoots (B,D,F,H) were analyzed. Values are mean ± SE (n = 3 replicates). Different letters indicate the significant differences (P < 0.05, LSD) between two cultivars or treatments.
Analysis of Antioxidants, Proline and ROS Accumulation
Antioxidant enzyme activities and proline concentrations were measured after plants were treated with NaCl. SOD, CAT, GR, and APX activities showed no difference between two cultivars under control condition. SOD, CAT, and GR were increased after salt treatment, and higher levels were observed in roots and shoots in Magnum Salt than that in Adrenalin, except for no difference in GR activity in shoots (Figures 4A–F). APX activity was increased in both cultivars after salt treatment and showed no difference in roots, but higher level was maintained in shoots in Adrenalin than in Magnum Salt (Figures 4G,H). No significant difference in proline concentration was observed between two cultivars under control condition. It was increased greatly in roots and shoots in both cultivars, while higher levels were maintained in both roots and shoots in Magnum Salt than in Adrenalin after salt treatment (Figures 4I,J).
Figure 4. Effect of salt treatment on the antioxidant enzyme activities and proline content in Magnum Salt and Adrenalin. Four-week old seedlings grown in soil with a mixture of pearlite, vermiculite and perlite (1:1:2) were treated with 0 or 500 mM NaCl for 6 days. SOD, CAT, GR, and APX activities and proline contents were analyzed in both root (A,C,E,G,I) and shoot (B,D,F,H,J). Values are mean ± SE (n = 3 replicates). Different letters indicate the significant differences (P < 0.05, LSD) between two cultivars or treatments.
Leaves were sampled for detection of ROS accumulation using NBT and DAB staining. Compared to the control leaves, O2⋅– and H2O2 accumulation was observed in leaves in both cultivar after salt treatment, while lower levels were maintained in Magnum Salt than in Adrenalin (Figure 5).
Figure 5. Effect of salt treatment on the ROS content in Magnum Salt and Adrenalin. Four-week old seedlings grown in soil with a mixture of pearlite, vermiculite and perlite (1:1:2) were treated with 0 or 500 mM NaCl for 6 days. Leaves were sampled for NBT and DAB staining.
Analysis of Fatty Acids Content
Nearly 90% of total fatty acids were C16 and C18 fatty acids in alfalfa roots under control condition, among them were major C18:2 (40%), C18:3 (20%) and C16:0 (20%). The portion of C18:2 increased significantly after 6 d of salt stress in Magnum Salt, but it was slightly decreased after 3 d of salt treatment, followed by a recovery at 6 d in Adrenalin. C18:3 content continuously decreased after salinity treatment, with no significant difference between two cultivars. C16:0 fatty acid was higher in Magnum Salt than in Adrenalin under control condition, and it showed no significant difference between two cultivars under salt stress (Figure 6A). Both saturated fatty acid (SFA) and unsaturated fatty acid (UFA) abundance were not altered in Magnum Salt after salinity treatment, while SFA was slightly increased but UFA was decreased in Adrenalin (Figure 6B). The monounsaturated fatty acids (MUFA) abundance was increased under salt stress, while polyunsaturated fatty acid (PUFA) was decreased, with lower level in Adrenalin. Double bond index (DBI; number of double bonds per mole) showed no difference between two cultivars under control condition, but it was decreased after salinity treatment with higher level in Magnum Salt than in Adrenalin (Figure 6C).
Figure 6. Effect of salt treatment to the Fatty acids content in Magnum Salt and Adrenalin. Four-week old seedlings treated with 0 or 500 mM NaCl. At 0, 3 and 6 days, roots were sampled for fatty acids content analysis (A). The upper and lower edges of the box in the boxplot represent the maximum and minimum values, respectively, and the horizontal line in the middle represents the median values. Relative abundance of saturated fatty acids (SFA), unsaturated fatty acid (UFA), monounsaturated fatty acids (MUFA), polyunsaturated fatty acid (PUFA); (B) and double bond index (DBI; number of double bonds per mole) (C) were calculated. Values are mean ± SE (n = 3 replicates). Different letters indicate the significant differences (P < 0.05, LSD) between two cultivars or treatments.
Identification and Functional Annotation of Differentially Expressed Genes
In order to study the expression profiles of genes in response to salinity, cDNA libraries representing samples of two cultivars (S-Adrenalin; T-Magnum Salt) at different time points (0, 2, and 26 h) were constructed for high-throughput RNA-Seq. A total of 740,274,514 raw reads were obtained (Supplementary Table 1).
A total of 735,954,286 clean reads remained after removing the adaptor sequences and low-quality sequences. Then, 79,712 unigenes with average GC content of 40.71% and N50 value of 1260 bp were obtained (Supplementary Table 1). By searching against four databases (Nr, Swiss-Prot, KOG and KEGG), 53171 unigenes were annotated (Supplementary Table 2). A total of 47365 DEGs matching the thresholds of false discovery rate (FDR) < 0.05 and | Log2Fold-Change| > 4 were identified (Figure 7). At 2 h, 12558 and 7399 DEGs were detected in Adrenalin and Magnum Salt, with 1684 DEGs in common. After 26 h treatment, 12227 and 15181 DEGs were identified in Adrenalin and Magnum Salt, with 6664 DEGs in common. There were 1684 and 6664 genes that differently expressed in both cultivars at 2 and 26 h.
Figure 7. Differentially expressed genes (DEGs) in Magnum Salt (T) and Adrenalin (S) under salt stress. (A) Number of upregulated and downregulated DEGs under salt stress. (B) Venn diagram analysis of the DEGs in different comparison groups of S0-vs-S2, S0-vs-S26, T0-vs-T2 and T0-vs-T26.
Differentially Expressed Genes Temporal Expression Patterns Analysis
Temporal expression patterns of the DEGs at 2 and 26 h after salt treatment were analyzed with the Short Time-series Expression Miner (STEM) software (v1.3.11). Eight different profiles of DEGs identified in STEM analysis were explored and enriched with KEGG analysis (Supplementary Table 3). Profile 0, exhibited a consistent decreasing trend, contains 6464 DEGs enriched in “phagosome,” “endocytosis,” “TCA cycle,” “carbon metabolism,” “glyoxylate and dicarboxylate metabolism,” “steroid biosynthesis,” “tryptophan metabolism,” “oxidative phosphorylation,” “arginine and proline metabolism” in Magnum Salt, and 2845 DEGs enriched in “phagosome,” “oxidative phosphorylation,” “isoquinoline alkaloid biosynthesis,” “cysteine and methionine metabolism,” “ribosome biogenesis,” “TCA cycle,” “steroid biosynthesis,” “cutin, suberine and wax biosynthesis,” “sesquiterpenoid and triterpenoid biosynthesis” in Adrenalin. Profile 7 exhibited an opposite trend. In Magnum Salt, 912 continuously upregulated DEGs were contained in Profile 7, and were enriched mainly in “linoleic acid metabolism,” “MAPK signaling pathway” and “hormone signal transduction” and so on. However, 393 DEGs in Adrenalin were contained in Profile 7, and were enriched mainly in “pyruvate metabolism,” “endocytosis,” “glycolysis.” Notably, “phagosome,” “TCA cycle” and “oxidative phosphorylation” pathways were frequently present in enrichment results of profile 0, 1, and 3, indicating that vesicle transport and aerobic metabolism were inhibited under salinity stress. The DEGs in profile 4, 6, and 7 in Magnum Salt were frequently enriched in “fatty acid metabolism,” “MAPK signaling pathway” and “hormone signal transduction” pathways, implying that fatty acid metabolization and signal transduction played an important role in coping with salt in alfalfa.
Differentially Expressed Genes Involved in Na+, K+ and Ca2+ Transport, ROS Scavenging, and Fatty Acid Metabolization
DEGs involved in K+/Na+/Ca2+ transport were investigated in our research (Figure 8A and Supplementary Table 4). The levels of AKT1, NHX2, HKT1 and HKT6 that involved in sodium/potassium homeostasis were higher in Magnum Salt after 26 h salinity, nevertheless, NHX7, AKT2, HAK23 and KAT3 always had higher expression in Adrenalin. Six unigenes annotated as CNGC3 and CNGC20 were expressed at higher levels in Magnum Salt than in Adrenalin at 2 h. At 26 h, ACA2, ACA8 and ACA12 involved in Ca2+ influx or efflux pathways were upregulated and both ACA2 and ACA12 were expressed at higher levels in Magnum Salt.
Figure 8. Heatmap of the relative expression of the DEGs involved in Na+, K+ and Ca2+ transport (A), ROS-scavenging (B) and fatty acid metabolization (C) in Magnum Salt and Adrenalin under salt stress. FPKM was used for representing the expression abundance of genes. The clustered heatmap was portrayed after normalized with zero-to-one scale method using the TBtools software.
Differentially expressed genes encoding SOD, CAT, APX and GR were analyzed. Most of them were upregulated after 2 or 26 h of salt treatment in Magnum Salt, while less genes were upregulated in Adrenalin (Figure 8B and Supplementary Table 5).
DEGs involved in fatty acid metabolization pathway were highly enriched in profile 7 in Magnum Salt (Supplementary Table 3). KCS2, ACAA1, ACAA1B, LACS4, FAD2 and FAD7 were upregulated to higher levels in Magnum Salt by NaCl. Nevertheless, the expression levels of GPAT5 and GPAT6 increased more strongly in Adrenalin following salt stress (Figure 8C and Supplementary Table 6). The transcript levels of KCS1, LOX1.1 and LOX1.4 were induced by salinity stress, but they were expressed at higher levels in Adrenalin.
RNA-Seq Results Validation With qRT-PCR Analysis
The relative folder change of eleven randomly selected DEGs were analyzed by qRT-PCR to verify the RNA-Seq results. Primers were listed in Supplementary Table 7. A high consistency was detected for these genes between qRT-PCR and RNA-seq results (Figures 9A–K). As shown in Figure 9L, a strong positive correlation (R2 = 0.8917) was obtained by a linear regression analysis, suggesting that our transcriptome data were reliable.
Figure 9. RNA-seq data validation. Random selected genes with different expression patterns were analyzed by qRT-PCR (A–K), correlation between qRT-PCR and RNA-seq data are calculated (L). Values are mean ± SE (n = 3 replicates).
Discussion
By analyzing the germination rates, survival rates and biomass under salt treatment, “Magnum Salt” and “Adrenalin” were identified as salt tolerant and salt sensitive cultivar, respectively (Figures 1, 2). Both physiological and transcriptional responses of these two cultivars to salinity stress were investigated. It was found that regulation of ionic homeostasis, fatty acid metabolization and antioxidative system played important roles in salt tolerance in alfalfa.
Ion Homeostasis
Maintaining ion homeostasis, especially a low cytosolic Na+/K+ ratio, is an important adaptive trait of salt-tolerant plants, which was also observed in salt-tolerant cultivar Magnum Salt in our experiments (Figures 3G,H). The transporters or channels involved in Na+ efflux and K+ uptake play important roles in salt resistance. NHX-type cation/H+ transporters in plants have been shown to mediate Na+(K+)/H+ exchange for salinity tolerance and K+ homoeostasis, providing an important strategy for ionic homeostasis in plants under saline conditions (Apse et al., 1999; Zeng et al., 2018). There were two salt-upregulated NHX genes identified in DEGs, NHX2 and NHX7. NHX2 expression abundance was about three times higher than NHX7 in all samples collected in our research (Supplementary Table 4). Notably, NHX2 increased greater in Magnum Salt than in Adrenalin. High-affinity K+ channel (HKT) and a low-affinity K+ channel Arabidopsis K+ Transporter1 (AKT1) were reported to mediate K+/Na+ influx (Ali et al., 2021, Wang et al., 2021). KAT3 encoded a member of the shaker family of voltage-gated potassium channel subunits and responded to low potassium stress in Arabidopsis, however, its function in salinity response had not been reported (Wang et al., 2016). The expression of AKT1, HKT1 and HKT6 were significantly higher in Magnum Salt after 26 h salinity, while genes encoding HKT23 or KAT3 were highly expressed in Adrenalin (Figure 8A).
Plant response to salt is closely linked to the calcium (Ca2+) channels as well as Ca2+ sensing and signaling (Seifikalhor et al., 2019, Feng et al., 2018). A rapid rise in cytosolic Ca2+ level can be detected within seconds of exposure to salt stress (Jiang et al., 2019). In another hand, salinity inhibited nutrient uptake, including Ca2+, just as showed in this paper (Figures 3E,F). However, Ca2+ concentrations were inhibited in a higher degree in Adrenalin than Magnum Salt by NaCl. The cyclic nucleotide-gated ion channel (CNGC) family have been implicated in the uptake of cations such as Na+, K+ and Ca2+, so that mediate numerous biological processes ranging from plant development to stress tolerance. Among CNGC members in Arabidopsis, AtCNGC3 helps in seed germination and cation transport, whereas AtCNGC19 and AtCNGC20 are associated with salt tolerance (Gobert et al., 2006, Kugler et al., 2009. Both CNGC3 and CNGC20 were expressed at higher levels in Magnum Salt than in Adrenalin at 2 h. Autoinhibited Ca2+-ATPases (ACAs) can mediate Ca2+ efflux to help restore [Ca2+]cyt to resting basal levels, for example, AtACA1, 2, and 7 mediate Ca2+ efflux from cytoplasm to endoplasmic reticulum and contribute to growth and pollen fitness in a redundant way (Ishka et al., 2021). There were four ACAs that were upregulated by salt. Among them, ACA2 and ACA12 involved in Ca2+ homeostasis were expressed at higher levels in Magnum Salt at 26 h. It’s interesting to find out whether they are involved in salt tolerance at future.
Antioxidative System
Oxidative injury happened in about all biotic and abiotic stress as a secondary injury. Thus, the ability to scavenge ROS efficiently is vital for tolerance (Nadarajah, 2020). Salt-sensitive XJ had high levels of ROS and elevation of ROS-related enzyme activity in response to salinity stress, whereas the salt-tolerant ZM showed relatively lower levels of ROS production and unaltered activity of ROS-related enzymes (Lei et al., 2018). Bhattarai et al. (2021) found out that salt tolerance in alfalfa appears to be associated with consistent expression of genes for enhancing defense against oxidative damage and protection against lipid peroxidation. In our research, salinity stress caused accumulation of O2⋅– and H2O2 in alfalfa leaves and upregulation of SOD, CAT, GR and APX activities in both roots and shoots. Magnum Salt accumulated less ROS in leaves than Adrenalin, corresponding with the higher SOD and CAT activities in Magnum Salt than in Adrenalin under salt treatment. Higher levels of SOD, CAT and GR were observed in roots in Magnum Salt than that in Adrenalin (Figures 4A–H and 5). Transcriptional levels of genes encoding antioxidative enzymes were quickly upregulated after salt treatment in Magnum Salt, indicating its fast response to apply antioxidative enzymes to scavenge ROS (Figure 8B and Supplementary Table 5).
Fatty Acids Metabolization
Membrane fluidity and integrity are largely affected by lipid composition and the degree of fatty acid desaturation in plants (Mikami and Murata, 2003). Alterations in membrane lipids in response to salinity have been observed in a number of plant species including both halophytes and glycophytes (Guo et al., 2019). Membrane fatty acids unsaturation was suggested to affect cell salinity tolerance by affect membrane fluidity so that corelated to Na+ and Cl– transport (Lin and Wu, 1996). In peanut, unsaturated fatty acids improved salt tolerance by alleviating photoinhibition (Liu et al., 2017). It has been reported in buffalograss (Buchloe dactyloides; Lin and Wu, 1996), corn (Zea mays L.; Hajlaoui et al., 2009) and Fabaceae (Bejaoui et al., 2016) that salt-tolerant clones or species had higher unsaturated fatty acids ratio than salt-sensitive ones after salinity treatment. However, no research has investigated the fatty acids contents under salt stress in alfalfa. In our research, it was found that palmitic (16:0), linoleic (18:2) and linolenic acid (18:3) accounted for the largest proportion of the fatty acid in alfalfa (Figure 6A), which was conserved across plant species (Guo et al., 2019). Salinity stress caused a decrease in total unsaturated fatty acids in salt-sensitive Adrenalin, which has been observed in many other plants. However, no significant difference existed in unsaturated fatty acids abundance in Magnum Salt between control or treated group (Figure 6B). Among unsaturated fatty acids, linoleic acid in Magnum Salt increased significantly after 6 days’ salinity. However, no significant difference existed in Adrenalin (Figure 6A).
There are four types of key enzyme in fatty acids synthesis, including fatty acid synthases (FAS), 3-ketoacyl-CoA thiolase (ACAA), β-ketoacyl-coenzyme A (CoA) synthases (KCS) and fatty acid desaturase (FAD). ACAA family members are responsible for the thiolytic cleavage of straight chain 3-oxoacyl-CoAs. KCS mediates the synthesis of very-long-chain fatty acids (VLCFAs) from 22 to 26 carbons in length. Research showed that KCS1 was salt-inducible, and might play a role in adapting intracellular membrane compartments to function in balancing the external osmotic pressure in extreme halotolerant alga Dunaliella salina (Azachi et al., 2002). FAD family members are responsible for fatty acid unsaturation. Endoplasmic reticulum localized FAD2 and FAD6 in plastids convert oleic acid (18:1) to linoleic acid (18:2) by inserting a double bond at the ω-6 position. Whereas FAD3, FAD7, and FAD8 convert linoleic acid (18:2) to linolenic acid (18:3). By investigating the phenotype of mutation materials, FAD2 and FAD6 were proved to be required for salt tolerance probably by maintaining proper function of membrane attached Na+/H+ exchangers in Arabidopsis (Zhang et al., 2012). Lipoxygenase (LOX) is also involved in fatty acid metabolization by catalyzing the hydroperoxidation of linolenic acid. High lipoxygenase activity and its further upregulation by salt stress are the unique features of salt-sensitive sunflower seedlings (Gogna and Bhatla, 2020). In our result, fatty acid metabolization pathway were highly enriched in profile 0 DEGs in Magnum Salt (Supplementary Table 3). Expression levels of ACAA1, ACAA1b, KCS1, FAD2, and FAD7 were upregulated by salinity to higher levels in Magnum Salt than that in Adrenalin, except for KCS2, which was highly expressed in Adrenalin. The transcript levels of LOX1.1 and LOX1.4 were higher in Adrenalin than Magnum Salt (Figure 8C and Supplementary Table 6). These results might partially explain the higher content of unsaturated fatty acids in Magnum Salt under salt stress.
Conclusion
In conclusion, compared to Adrenalin, salt-tolerant Magnum Salt applied different strategies to cope with salinity condition by regulation of ionic homeostasis, antioxidative enzyme activities and fatty acid metabolization at both transcriptional and physiological level (Figure 10).
Figure 10. A proposed model showing salt tolerance mechanisms in alfalfa. Salt tolerance in alfalfa is associated with regulation of ionic homeostasis, antioxidative enzymes and fatty acid metabolism. Red and blue genes are proposed positive and negtive regulatory genes in salt tolerance.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author Contributions
ZG and HZ designed the experiments. JL, MM, YS, and PL performed the experiments. JL, HZ, and HS contributed to the data analysis. HZ and ZG wrote and revised this manuscript. All authors approved the submitted version.
Funding
This work was supported by the National Natural Science Foundation of China (Grant numbers 32101439) and the Fundamental Research Funds for the Central Universities (KYQN2022015).
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.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.931619/full#supplementary-material
Footnotes
- ^ http://www.r-project.org/
- ^ http://systemsbiology.cau.edu.cn/agriGOv2/
- ^ http://kobas.cbi.pku.edu.cn/
References
Ali, A., Raddatz, N., Pardo, J. M., and Yun, D. J. (2021). HKT sodium and potassium transporters in Arabidopsis thaliana and related halophyte species. Physiol. Plant. 171, 546–558. doi: 10.1111/ppl.13166
An, Y. M., Yang, X. X., Zhang, L. S., Zhang, J., Du, B. H., Yao, L., et al. (2020). Alfalfa MsCBL4 enhances calcium metabolism but not sodium transport in transgenic tobacco under salt and saline-alkali stress. Plant Cell Rep. 39, 997–1011. doi: 10.1007/s00299-020-02543-x
Apse, M. P., Aharon, G. S., Snedden, W. A., and Blumwald, E. (1999). Salt tolerance conferred by overexpression of a vacuolar Na+/H+ antiport in Arabidopsis. Science 285, 1256–1258. doi: 10.1126/science.285.5431.1256
Azachi, M., Sadka, A., Fisher, M., Goldshlag, P., Gokhman, I., and Zamir, A. (2002). Salt induction of fatty acid elongase and membrane lipid modifications in the extreme halotolerant alga Dunaliella salina. Plant Physiol. 129, 1320–1329. doi: 10.1104/pp.001909
Bejaoui, F., Salas, J. J., Nouairi, I., Smaoui, A., Abdelly, C., Martinez-Force, E., et al. (2016). Changes in chloroplast lipid contents and chloroplast ultrastructure in Sulla carnosa and Sulla coronaria leaves under salt stress. J. Plant Physiol. 198, 32–38. doi: 10.1016/j.jplph.2016.03.018
Bhattarai, S., Fu, Y. B., Coulman, B., Tanino, K., Karunakaran, C., and Biligetu, B. (2021). Transcriptomic analysis of differentially expressed genes in leaves and roots of two alfalfa (Medicago sativa L.) cultivars with different salt tolerance. BMC Plant Biol. 21:446. doi: 10.1186/s12870-021-03201-4
Chen, C., Chen, H., Zhang, Y., Thomas, H. R., Frank, M. H., He, Y., et al. (2020). TBtools: An integrative toolkit developed for interactive analyses of big biological data. Mol. Plant. 13, 1194–1202. doi: 10.1016/j.molp.2020.06.009
Chen, T. T., Yang, Q. C., Zhang, X. Q., Ding, W., and Gruber, M. (2012). An alfalfa (Medicago sativa L.) ethylene response factor gene, MsERF11, enhances salt tolerance in transgenic Arabidopsis. Plant Cell Reports 31, 1737–1746. doi: 10.1007/s00299-012-1287-z
Cornacchione, M. V., and Suarez, D. L. (2017). Evaluation of alfalfa (Medicago sativa L.) populations’ response to salinity stress. Crop Sci. 57, 137–150. doi: 10.2135/cropsci2016.05.0371
Dai, M. T., Huang, R. S., Han, Y. Y., Zhang, Z. Y., Chen, Y. Y., Shi, H. F., et al. (2022). A novel salt responsive PvHAK16 negatively regulates salt tolerance in transgenic Arabidopsis thaliana. Environ. Exp. Bot. 194:104689. doi: 10.1016/j.envexpbot.2021.104689
De Vos, C. H. R., TenBookum, W. M., Vooijs, R., Schat, H., and De Kok, L. J. (1993). Effect of copper on fatty acid composition and peroxidation of lipids in the roots of copper tolerant and sensitive Silene-cucubalus. Plant Physiol. Biochem. 31, 151–158. doi: 10.1111/j.1399-3054.1986.tb05073.x
Deinlein, U., Stephan, A. B., Horie, T., Luo, W., Xu, G., and Schroeder, J. I. (2014). Plant salt-tolerance mechanisms. Trends Plant Sci. 19, 371–379. doi: 10.1016/j.tplants.2014.02.001
Ernst, J., and Bar-Joseph, Z. (2006). STEM: a tool for the analysis of short time series gene expression data. BMC Bioinformatics 7:191. doi: 10.1186/1471-2105-7-191
Feng, W., Kita, D., Peaucelle, A., Cartwright, H. N., Doan, V., Duan, Q. H., et al. (2018). The FERONIA receptor kinase maintains cell-wall integrity during salt stress through Ca2+ signaling. Curr. Biol. 28, 666–675 e5. doi: 10.1016/j.cub.2018.01.023
Gobert, A., Park, G., Amtmann, A., Sanders, D., and Maathuis, F. J. M. (2006). Arabidopsis thaliana Cyclic Nucleotide Gated Channel 3 forms a non-selective ion transporter involved in germination and cation transport. J. Exp. Bot. 57, 791–800. doi: 10.1093/jxb/erj064
Gogna, M., and Bhatla, S. C. (2020). Salt-tolerant and -sensitive seedlings exhibit noteworthy differences in lipolytic events in response to salt stress. Plant Signal. Behav. 15:1737451.
Gruber, M. Y., Xia, J., Yu, M., Steppuhn, H., Wall, K., Messer, D., et al. (2017). Transcript analysis in two alfalfa salt tolerance selected breeding populations relative to a non-tolerant population. Genome 60, 104–127. doi: 10.1139/gen-2016-0111
Guo, P., Shi, W. C., Li, L. L., and Bao, Y. J. (2018). MsPP2C, a protein phosphatase 2C gene of alfalfa, confers enhanced salt tolerance in Arabidopsis. Int. J. Agric. Biol. 20, 62–70. doi: 10.17957/Ijab/15.0362
Guo, Q., Liu, L., and Barkla, B. J. (2019). Membrane lipid remodeling in response to salinity. Int. J. Mol. Sci. 20:4264. doi: 10.3390/ijms20174264
Hajlaoui, H., Denden, M., and Ayeb, N. E. (2009). Changes in fatty acids composition, hydrogen peroxide generation and lipid peroxidation of salt-stressed corn (Zea mays L.) roots. Acta Physiol. Plant. 31, 787–796. doi: 10.1007/s11738-009-0293-4
Huang, J. X., Xue, C. W., Wang, H., Wang, L. S., Schmidt, W., Shen, R. F., et al. (2017). Genes of acyl carrier protein family show different expression profiles and overexpression of acyl carrier protein 5 modulates fatty acid composition and enhances salt stress tolerance in Arabidopsis. Front. Plant Sci. 8:987. doi: 10.3389/fpls.2017.00987
Ishka, M. R., Brown, E., Rosenberg, A., Romanowsky, S., Davis, J. A., Choi, W. G., et al. (2021). Arabidopsis Ca2+-ATPases 1, 2, and 7 in the endoplasmic reticulum contribute to growth and pollen fitness. Plant Physiol. 185, 1966–1985. doi: 10.1093/plphys/kiab021
Jiang, Z. H., Zhou, X. P., Tao, M., Yuan, F., Liu, L. L., Wu, F. H., et al. (2019). Plant cell-surface GIPC sphingolipids sense salt to trigger Ca2+ influx. Nature 572, 341–346. doi: 10.1038/s41586-019-1449-z
Kugler, A., Kohler, B., Palme, K., Wolff, P., and Dietrich, P. (2009). Salt-dependent regulation of a CNG channel subfamily in Arabidopsis. BMC Plant Biol. 9:140. doi: 10.1186/1471-2229-9-140
Lei, Y., Xu, Y., Hettenhausen, C., Lu, C., Shen, G., Zhang, C., et al. (2018). Comparative analysis of alfalfa (Medicago sativa L.) leaf transcriptomes reveals genotype-specific salt tolerance mechanisms. BMC Plant Biol. 18:35. doi: 10.1186/s12870-018-1250-4
Li, C. X., Song, T. T., Zhan, L. F., Cong, C. L., Xu, H. H., Dong, L., et al. (2021). Overexpression of MsRCI2A, MsRCI2B, and MsRCI2C in alfalfa (Medicago sativa L.) provides different extents of enhanced alkali and salt tolerance due to functional specialization of MsRCI2s. Front. Plant Sci. 12:1836. doi: 10.3389/fpls.2021.702195
Lin, H., and Wu, L. (1996). Effects of salt stress on root plasma membrane characteristics of salt-tolerant and salt-sensitive buffalograss clones. Environ. Exp. Bot. 36, 239–254. doi: 10.1016/0098-8472(96)01025-8
Liu, L. L., and Wang, B. S. (2021). Protection of halophytes and their uses for cultivation of saline-alkali soil in China. Biology 10:353. doi: 10.3390/biology10050353
Liu, S. S., Wang, W. Q., Li, M., Wan, S. B., and Sui, N. (2017). Antioxidants and unsaturated fatty acids are involved in salt tolerance in peanut. Acta Physiol. Plant. 39, 1–10. doi: 10.1007/s11738-017-2501-y
Long, R. C., Yang, Q. C., Kang, J. M., Zhang, T. J., Wang, H. M., Li, M. N., et al. (2013). Overexpression of a novel salt stress-induced glycine-rich protein gene from alfalfa causes salt and ABA sensitivity in Arabidopsis. Plant Cell Rep. 32, 1289–1298. doi: 10.1007/s00299-013-1443-0
Luo, D., Zhou, Q., Wu, Y. G., Chai, X. T., Liu, W. X., Wang, Y. R., et al. (2019). Full-length transcript sequencing and comparative transcriptomic analysis to evaluate the contribution of osmotic and ionic stress components towards salinity tolerance in the roots of cultivated alfalfa (Medicago sativa L.). BMC Plant Biol. 19:32. doi: 10.1186/s12870-019-1630-4
Ma, L. C., Wang, Y. R., Liu, W. X., and Liu, Z. P. (2014). Overexpression of an alfalfa GDP-mannose 3, 5-epimerase gene enhances acid, drought and salt tolerance in transgenic Arabidopsis by increasing ascorbate accumulation. Biotechnol. Lett. 36, 2331–2341. doi: 10.1007/s10529-014-1598-y
Mikami, K., and Murata, N. (2003). Membrane fluidity and the perception of environmental signals in cyanobacteria and plants. Prog. Lipid Res. 42, 527–543. doi: 10.1016/s0163-7827(03)00036-5
Nadarajah, K. K. (2020). ROS homeostasis in abiotic stress tolerance in plants. Int. J. Mol. Sci. 21:5208. doi: 10.3390/ijms21155208
Seifikalhor, M., Aliniaeifard, S., Shomali, A., Azad, N., Hassani, B., Lastochkina, O., et al. (2019). Calcium signaling and salt tolerance are diversely entwined in plants. Plant Signal. Behav. 14:1665455. doi: 10.1080/15592324.2019.1665455
Sun, G. L., Zhu, H. F., Wen, S. L., Liu, L. S., Gou, L. M., and Guo, Z. F. (2020). Citrate synthesis and exudation confer Al resistance in alfalfa (Medicago sativa L.). Plant Soil 449, 319–329. doi: 10.1007/s11104-020-04490-8
Wang, X. P., Chen, L. M., Liu, W. X., Shen, L. K., Wang, F. L., Zhou, Y., et al. (2016). AtKC1 and CIPK23 synergistically modulate AKT1-mediated low-potassium stress responses in Arabidopsis. Plant Physiol. 170, 2264–2277. doi: 10.1104/pp.15.01493
Wang, X. S., Zhao, J. L., Fang, Q. W., Chang, X. C., Sun, M. Y., Li, W. B., et al. (2021). GmAKT1 is involved in K+ uptake and Na+/K+ homeostasis in Arabidopsis and soybean plants. Plant Sci. 304:110736. doi: 10.1016/j.plantsci.2020.110736
Wang, Y. J., Jiang, L., Chen, J. Q., Tao, L., An, Y. M., Cai, H. S., et al. (2018). Overexpression of the alfalfa WRKY11 gene enhances salt tolerance in soybean. PLoS One 13:e0192382. doi: 10.1371/journal.pone.0192382
Yu, R. G., Wang, G. L., Yu, X. Y., Li, L. Y., Li, C. X., Song, Y. X., et al. (2021). Assessing alfalfa (Medicago sativa L.) tolerance to salinity at seedling stage and screening of the salinity tolerance traits. Plant Biol. 23, 664–674. doi: 10.1111/plb.13271
Zeng, Y., Li, Q., Wang, H. Y., Zhang, J. L., Du, J., Feng, H. M., et al. (2018). Two NHX-type transporters from Helianthus tuberosus improve the tolerance of rice to salinity and nutrient deficiency stress. Plant Biotechnol. J. 16, 310–321. doi: 10.1111/pbi.12773
Zhang, J. T., Liu, H., Sun, J., Li, B., Zhu, Q., Chen, S. L., et al. (2012). Arabidopsis fatty acid desaturase FAD2 is required for salt tolerance during seed germination and early seedling growth. PLoS One 7:e30355. doi: 10.1371/journal.pone.0030355
Zhang, Z. Q., Wang, Y. F., Chang, L. Q., Zhang, T., An, J., Liu, Y. S., et al. (2016). MsZEP, a novel zeaxanthin epoxidase gene from alfalfa (Medicago sativa), confers drought and salt tolerance in transgenic tobacco. Plant Cell Reports 35, 439–453.
Zhuo, C. L., Liang, L., Zhao, Y. Q., Guo, Z. F., and Lu, S. Y. (2018). A cold responsive ethylene responsive factor from Medicago falcata confers cold tolerance by up-regulation of polyamine turnover, antioxidant protection, and proline accumulation. Plant Cell Environ. 41, 2021–2032. doi: 10.1111/pce.13114
Keywords: alfalfa, salt tolerance, transcriptome, ionic homeostasis, antioxidant defense system, unsaturated fatty acids
Citation: Li J, Ma M, Sun Y, Lu P, Shi H, Guo Z and Zhu H (2022) Comparative Physiological and Transcriptome Profiles Uncover Salt Tolerance Mechanisms in Alfalfa. Front. Plant Sci. 13:931619. doi: 10.3389/fpls.2022.931619
Received: 29 April 2022; Accepted: 19 May 2022;
Published: 09 June 2022.
Edited by:
Maofeng Chai, Qingdao Agricultural University, ChinaCopyright © 2022 Li, Ma, Sun, Lu, Shi, Guo and Zhu. 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: Haifeng Zhu, emh1aGFpZmVuZ0BuamF1LmVkdS5jbg==