De Novo Transcriptional Analysis of Alfalfa in Response to Saline-Alkaline Stress

Saline-alkaline stress, caused by high levels of harmful carbonate salts and high soil pH, is a major abiotic stress that affects crop productivity. Alfalfa is a widely cultivated perennial forage legume with some tolerance to biotic and abiotic stresses, especially to saline-alkaline stress. To elucidate the mechanism underlying plant saline-alkaline tolerance, we conducted transcriptome analysis of whole alfalfa seedlings treated with saline-alkaline solutions for 0 day (control), 1 day (short-term treatment), and 7 days (long-term treatment) using ion torrent sequencing technology. A transcriptome database dataset of 53,853 unigenes was generated, and 2,286 and 2,233 genes were differentially expressed in the short-term and long-term treatment, respectively. Gene ontology analysis revealed 14 highly enriched pathways and demonstrated the differential response of metabolic pathways between the short-term and long-term treatment. The expression levels of 109 and 96 transcription factors were significantly altered significantly after 1 day and 7 days of treatment, respectively. Specific responses of peroxidase, flavonoids, and the light pathway component indicated that the antioxidant capacity was one of the central mechanisms of saline-alkaline stress tolerance response in alfalfa. Among the 18 differentially expressed genes examined by real time PCR, the expression levels of eight genes, including inositol transporter, DNA binding protein, raffinose synthase, ferritin, aldo/keto reductase, glutathione S-transferase, xyloglucan endotrans glucosylase, and a NAC transcription factor, exhibited different patterns in response to saline and alkaline stress. The expression levels of the NAC transcription factor and glutathione S-transferase were altered significantly under saline stress and saline-alkaline stress; they were upregulated under saline-alkaline stress and downregulated under salt stress. Physiology assays showed an increased concentration of reactive oxygen species and malondialdehyde and a decreased content of chlorophyll, indicating that anti-oxidation and detoxification play an important role in response to saline-alkaline stress. Overall, the transcriptome analysis provided novel insights into the saline-alkaline stress tolerance response mechanisms in alfalfa.

Introduction

Abiotic stresses, including drought, salinity, and extreme temperatures, adversely affect crop production worldwide. Soil salinity and alkalinity are among the major abiotic factors that limit crop yield. More than 954 million hectares of the global land surface consist of sodic soil (Zhang et al., 2006). In the Song-nen Plain, the main crop production area of Northeast China, approximately 3.78 million hectares of land are threatened by soil salinization and alkalization (Yang and Wang, 2015). Sodic soils with high pH and salt concentration restrict plant growth and limit agricultural production in this area.

Plant damage caused by saline stress is usually divided into three main categories: osmotic shock, toxic cation stress, and high pH damage. High salt concentrations reduce the soil water potential, resulting in osmotic shock (Wong et al., 2006), which is a state when an excess of sodium ions inhibits the metabolic processes in the cells and generates secondary metabolites such as peroxides (Zhu, 2001). Numerous studies have focused on plant salt tolerance; however, information on the combined saline and alkaline stress (hereafter referred to as saline-alkaline stress) is still limited. Thus, the gene expression regulation and molecular mechanisms involved in plant saline-alkaline stress tolerance are yet to be discovered. Saline-alkaline stress shares many factors with saline stress, such as excess of Na⁺ and water deficit, that trigger similar responses in both stress types (Gong et al., 2014). However, alkaline stress is unique for its high pH levels that initiate the synthesis of reactive oxygen species (ROS) and malondialdehyde that damage the plant cell membrane and intracellular components. High pH due to soil alkalization is a crucial and critical stress condition compared with salinity stress, because increased ion concentrations have more adverse effects on plant growth than soil salinization alone. Therefore, alkaline stress is a more diversified form of stress, and plant tolerance to alkalinity reflects a more complex tolerance mechanism (Peng et al., 2008).

Plant saline-alkaline stress tolerance mechanisms involve the expression of a cluster of genes and the interaction among gene products rather than the expression of individual genes (Xu and Tuyen do, 2012; Gong et al., 2013). Gene expression is affected by many internal and external factors (Yang et al., 2009); hence, a more comprehensive understanding of saline-alkaline stress tolerance mechanisms in plants should be inferred at the gene expression level. Transcriptome analysis is a novel tool for the study of plant saline-alkaline stress tolerance mechanisms (Sun et al., 2013; Wang C. et al., 2014) at the RNA level and uncovers the differences in common stress response pathways as a response to a variety of stresses such as saline, peroxide, and saline-alkaline stress.

In recent years, high throughput technologies have been rapidly developed and used in the genome-wide gene expression profiling of Arabidopsis (Jiang and Deyholos, 2006), rice (Kawasaki et al., 2001), wheat (Kawaura et al., 2006), and soybean (Irsigler et al., 2007; Ge et al., 2010). Next generation sequencing is a large-scale analysis method of differential gene expression using various platforms (Genome Analyzer [Solexa/Illumina], 454 [Roche], ABI-SOLiD [Applied Biosystems], and HeliScope [SeqLL LLC, Woburn, MA, USA]) that facilitate RNA-seq analysis and unravel a variety of stress responses on a transcriptome-wide scale in non-model plant species. Ion torrent sequencing detects hydrogen ions that are released during DNA polymerization, and it differs from other sequencing technologies, since it does not employ any modified nucleotides or optics. This technology is a rapid, compact, and economical sequencing system that can be used for plant transcriptome analysis.

Alfalfa (Medicago sativa) is a major forage legume cultivated worldwide, including the Song-nen Plain, Northeast China (Li et al., 1996). ‘Zhaodong’ is an alfalfa cultivar with resistance to various abiotic stresses such as freezing, drought, salinity, and saline-alkaline stress (Peng et al., 2008). Despite its economic importance, the physiological and molecular mechanisms underlying its response to abiotic stresses are poorly understood. The completed sequencing project of M. truncatula, a close relative of M. sativa, provides a reference genome to identify tentative consensus sequences in alfalfa in response to abiotic stress and generate an expression pattern of stress response.

In this paper, we conducted transcriptional profiling of the whole seedling alfalfa treated with saline-alkaline solution using ion torrent to assess changes in gene expressions. Differentially expressed gene (DEG) profiles as responses to stress were presented, and DEG response to saline-alkaline stress was annotated. DEG-specific responses to saline-alkaline conditions were identified through comparison with the saline stress response pathway, and specific mechanisms of saline-alkaline stress tolerance were elucidated. The transcriptome analysis of alfalfa in response to saline-alkaline stress might provide novel insights into plant saline-alkaline stress tolerance mechanism based on DEGs and reveal the related pathway of plant saline-alkaline tolerance.

Materials and Methods

Plant Materials and Treatment Conditions

Seeds of Zhaodong, supplied by the Heilongjiang Animal Science Institute, China, were germinated in pots with distilled water for 2 days in the dark, transferred to vermiculite with ½ Hoagland solution, and kept in a culture room for 4 weeks at 60% relative humidity, 22°C, and a photoperiod of 16 h light/8 h dark.

Four-week-old seedlings were treated with ½ Hoagland solution, containing 0.1 M mixed alkaline solution (Na₂CO₃: NaHCO₃, 1: 2) for 0 day (control), 1 day (short-term treatment), and 7 days (long-term treatment). The harvested whole seedlings from each treatment were frozen in liquid nitrogen and stored at –80°C for RNA-seq, qRT-PCR, and physiological analyses. Ion torrent sequencing was conducted by BGI (Shenzhen, China) in two replicates for each sample.

The same batch of whole-seedling samples was treated with 0.1 M mixed alkaline and 0.1 M NaCl for 7 days prior to qRT-PCR in order to identify the differences in response mechanisms under saline stress and saline-alkaline stress conditions.

RNA Extraction and Transcriptome Sequencing

Transcriptome sequencing was completed by BGI using Ion Proton I RNA-seq (Thermo-Fisher Scientific, Waltham, MA, USA). Total RNA was extracted from the control and samples treated with 0.1 M mixed saline-alkaline solution using One Step RNA Reagent (Bio Basic Inc., Markham, ON, Canada), according to the manufacturer’s protocol, and purified using the RNeasy Plant Mini Kit (Qiagen, Valencia, CA, USA). The whole transcriptome cDNA library was prepared using the Ion Total RNA-Seq kit V2 (Life Technologies Corp., Carlsbad, CA, USA). Double-stranded cDNA was ligated to barcoded adapters and then sequenced using an Ion PI^TM Chip (Ion torrent, Life Technologies, Carlsbad, CA, USA) by BGI. Raw data processing, adapter sequence removal, base-calling, and quality value calculations were performed using Torrent Suite^TM 4.0 (Life Technologies, Carlsbad, CA, USA). Quality reads were obtained by trimming the raw reads at a minimum PHRED score Q = 20 (Shu et al., 2015).

Reads Assembly and Functional Annotation

Total reads obtained by Ion torrent sequencing were assembled by Trinity assembler (Grabherr et al., 2011), and redundant reads were removed using iAssembler (Zheng et al., 2011). The contigs that presented zero transcript expression (RPKM) were adjusted to 1 in order to avoid negating the subsequent calculation of fold change. The differentially expressed transcripts were identified based on fold changes (FC ≥ 2 or FC ≤ 0.5) and statistical analysis (p ≤ 0.001) using edgeR package in R (R team, Vienna, Austria).

DEGs were annotated using two methods: (1) Reads generated by transcriptome analysis were mapped onto M. truncatula, a species with a completely sequenced genome (Mt4.0 data release; Tang et al., 2014) highly homologous to alfalfa, using Bowtie 0.12.7 (Langmead et al., 2009); (2) Functional annotation of unigenes assembled by Trinity was conducted using BLASTx search against the Nr and Swiss-Prot databases and M. truncatula, Arabidopsis thaliana, and Glycine max as reference sequences at a cut-off E-value of 10e-5. The unigenes were assigned to known functional groups and biological processes using Gene Ontology (GO) tools (Ye et al., 2006)¹.

Transcription factors were searched and identified from unigenes using iTAK². Venny tool was used for drawing and comparison of Venn diagrams (Oliveros, 2007–2015).

Quantitative Reverse Transcription PCR (qRT-PCR) Analysis

Total RNA was extracted from the same group of samples that were used in transcriptome analysis using an RNA extraction kit (TianGen Biotech, Beijing, China). RNA was reverse-transcribed into cDNA using the PrimeScript™ RT reagent kit (TaKaRa, Shiga, Japan). Twenty primers (20 random genes and 20 stress-induced genes) were used in real-time RT-PCR, and data were normalized using the actin gene.

The reaction mixture (20 μl) contained 10 μl of SYBR Green RealTime PCR Master Mix (Toyobo, Osaka, Japan), 2 μl of cDNA template, and 0.5 μM of each forward and reverse primer. The amplification was run in an ABI 7300 sequencer with the following cycling parameters: initial denaturation at 94°C for 30 s, followed by 45 cycles at 94°C for 15 s, 55°C for 30 s, and 72°C for 30 s. Real-time PCR was carried out in triplicate to confirm the accuracy of the results. All the relative expression levels were log₂-transformed.

Measurement of Physiological Parameters under Saline-Alkaline Stress

The levels of H₂O₂ and ${\text{O}}_2^{-}$ were measured as described by Mukherjee and Choudhuri (1983) and Liu and Pang (2010), respectively. To measure the ${\text{O}}_2^{-}$ content, plant material was extracted in potassium phosphate buffer (pH 7.8) and then, centrifuged at 10,000 rpm for 5 min. After incubating the supernatant at 25°C for 20 min, 17 mM sulfanilamide and 7 mM a-naphthylamine were added, and the mixture was incubated at 25°C for additional 20 min. The absorbance was read at 530 nm, and the ${\text{O}}_2^{-}$ content was calculated against the standard curve prepared with NaNO₂. To determine the H₂O₂ content, plant tissue was extracted with cold acetone (4°C) and mixed with 0.1% TiCl₄ and ammonia. The mixture was centrifuged at 10,000 rpm for 10 min. The precipitate was dissolved in 2 mM H₂SO₄ and the absorbance was read at 415 nm; the H₂O₂ content was calculated using the standard curve.

Peroxidase (POD) activity was measured according to the protocol described by Gong (2001). The concentration of malondialdehyde (MDA) was measured as described by Hodges et al. (1999). Plant tissue was extracted in precooled acetocaustin and centrifuged at 10,000 rpm for 10 min. The supernatant was mixed with thiobarbituric acid and transferred to a boiling water bath for 15 min. After centrifuging at 10,000 rpm for 10 min, the concentration of MDA was calculated from the absorbance read at 532 nm and 450 nm.

The chlorophyll content was measured as previously described (Gong et al., 2013). Aldo/keto reductase (AKR) activity was determined as described by Wang Q. et al. (2014). Root tissues exposed to saline stress or saline-alkaline stress for 7 days were extracted in cold potassium phosphate buffer (pH 7.4) and centrifuged at 12,000 rpm for 10 min. The supernatant was added to a mixture of 5 mM D-glucuronic acid and dimethylsulfoxide, which was previously incubated at 30°C for 5 min, and 10 mg mol⁻¹ NADPH. The absorbance was read at 532 nm every minute for 3 min. The activity of AKR was calculated as follows: u/g = ΔA × V/(6.22 × M).

Results

Ion Torrent Sequencing and Reads Assembly

After removing low-quality adaptor and barcode sequences, 16,209,292; 14,444,593; and 19,940,365 draw reads were obtained from the control, the short-term treatment, and the long-term treatment, respectively. Of these, 15,774,180; 14,109,229; and 19,355,510, respectively, were clean reads (Table 1). The percentage of clean reads was 97.67% and 97.24% in the short-term treatment and the long-term treatment, respectively, indicating a high sequencing quality. Contigs from all samples were assembled into 53,853 unigenes that were used to generate a non-redundant unigene library for further analysis of the M. sativa transcriptome under saline-alkaline stress.

TABLE 1

TABLE 1. De novo transcriptome assembly of contig sequences of alfalfa exposed to saline-alkaline stress for different length of time.

Identification of DEG and Functional Classification

Functional annotation of non-redundant unigenes against the Nr database and the Swiss-Prot database using BLASTx resulted in 2,286 (1 day treatment) and 2,233 (7 days treatment) DEGs (log₂ > 1.5, p < 1e-5; Supplementary Table S1). In total, 1,561 upregulated DEGs and 725 downregulated DEGs were identified in the short-term treatment, whereas 1,599 upregulated DEGs and 634 downregulated DEGs were identified in the long-term treatment (Figure 1).

FIGURE 1

FIGURE 1. Venn diagram of the number of differentially expressed genes (DEGs) (log₂ > 1.5) in different time points in de novo transcriptome (false discovery rate < 0.025). (A) Distribution of all DEG (B) up-regulated DEG; and (C) down-regulated DEG annotated in samples treated with alkaline-solution for 1 day (red circles) and 7 days (blue circles).

qRT-PCR Conformation of Ion Torrent Sequencing

To validate transcriptome results, 40 genes (20 random and 20 stress response genes) differentially expressed in the short-term treatment and the long-term treatment were selected. The qRT-PCR results were consistent with the transcriptome sequencing results (Table 2).

TABLE 2

TABLE 2. Genes used for validation of next generation sequencing (NGS) data by quantitative reverse transcriptase PCR (qRT-PCR).

The most highly upregulated genes were delta-1-pyrroline-5-carboxylate synthetase (Medtr4g020110.1), encoding one of the key enzymes of proline biosynthesis, glutathione S-transferase (Medtr1g090060.1) that participates in the non-enzymatic elimination pathway of ROS, and flavonol synthase (Medtr5g055690.1) that is a key enzyme of flavonoid biosynthesis. The most highly downregulated genes were the light-harvesting protein (Medtr4g099340.1) and the light-harvesting complex I chlorophyll A/B-binding protein (Medtr6g011880.1) that play an important role in the capture of photons in the light pathway, the auxin-binding protein ABP19a (Medtr8g020590.1), and F-box (Medtr7g115270.1).

The plotting of RT-PCR and Ion torrent transcriptome data (Figure 2) showed a high correlation coefficient (R² = 0.732, P < 0.05 in samples treated with saline-alkaline solution for 1 day and R² = 0.756, P < 0.05 in samples treated with saline-alkaline solution for 7 days) between the two data sets, whole data see Supplementary Table S2.

FIGURE 2

FIGURE 2. Validation of next generation sequencing (NGS data by quantitative reverse transcriptase PCR (qRT-PCR) in samples treated with alkaline-saline solution for 1 day and 7 days.

GO Analysis

The GO cluster analysis of DEGs obtained from the two treatments was assembled using WEGO. The expression of genes involved in binding, antioxidant activity, membrane, and transport was similar between the two treatments. Subcategories of the biological process category, such as biopolymer modification, phosphate metabolism, lipid biosynthesis, and cellular macromolecule catabolism, were changed significantly between the two treatments. In the category of molecular functions, transferase activity, kinase activity, phosphotransferase activity, and structural constituent of ribosome had significantly different expression levels. In the category of cell components, ribosome, intracellular non-membrane-bound organelle, and cytosol were significantly different between the two treatments (Figure 3; Supplementary Table S3).

FIGURE 3

FIGURE 3. Histogram of gene ontology (GO) categories in response to 1 day and 7 days of treatment with saline-alkaline stress. The bars represent log₂-transformed ratio between category portions in up-regulated and down-regulated gene sets.

Gene Differential Expression in Short-Term and Long-Term Treatments

Gene ontology cluster analysis of the 53,853 unigenes showed that under alkaline-saline stress conditions, the expression levels of DEGs as well as the differential expression of various gene types were modified (Figure 4; Supplementary Table S4).

FIGURE 4

FIGURE 4. Hierarchical cluster analysis of differentially expressed genes based on expression pattern of 53,853 genes derived from the control, 1 day and 7 days treatments.

In the molecular function category, the expression levels of genes clustered in GO categories, such as protein kinase, structural molecule activity, transferase, and POD, were altered significantly. Among these, protein kinase, phosphotransferase, and transferase had higher expression levels in the short-term treatment than in the long-term treatment. These results indicated that the metabolic processes were mediated by kinase, and that the efficiency of intracellular transporter was greater at on day 1 of treatment than at day 7 of treatment. The activity of transcription factor was more enhanced at day 7 than at day 1 of treatment, whereas the expression level of POD was high in both treatments.

Genes involved in the structural constituent of the ribosome subcategory did not alter their expression in the short-term treatment; however, they were highly induced in the long-term treatment, showing that more ribosomes were assembled to translate proteins. These results suggested that intracellular homeostasis does not form in plant cells that are adapted to saline-alkaline stress, unless new homeostasis is established to meet the high expression levels of the induced genes.

Lipid biosynthetic process, a subcategory of the biology process category, had significantly more expressed genes in the long-term treatment than in the short-term treatment. The DEGs in the lipid biosynthetic process subcategory were focused on the expression and modification of phospholipids at day 1 of treatment and the transport and metabolism of lipids at day 7 of treatment with the sustained expression of inositol.

Differential Expression of Peroxidase under Saline-Alkaline Stress

In the transcriptome category, BLASTx against M. truncatula genome using de novo assembly clustered PODs into three classes: POD family protein, cationic POD, and class III POD. In addition, a few genes clustered into respiratory burst oxidase-like protein D, POD family protein, and cationic POD, which belong to class I POD and were induced by the short-term treatment and upregulated in the long-term treatment. For instance, c67928_g1_i1 (similar to Medtr2g008730.1) and c70153_g1_i3 (similar to Medtr1g066380.1) were both strongly induced and upregulated by stress, reaching the expression levels of 4.766 and 4.252, respectively. In the present study, the expression of the vast majority of class I POD was lower in the long-term treatment than in the short-term treatment, yet still higher than in the control. In contrast, the member of class III POD c69915_g1_i1 (similar to Medtr7g072510.1) was upregulated under stress, increasing the levels 4.53-fold at day 7 of treatment compared with the levels observed at day 1 of treatment.

Although the majority of PODs were upregulated under stress conditions, a different expression tendency was observed in few PODs. For instance, c41991_g1_i1 (similar to Medtr4g046713.1) and c60608_g1_i1 (similar to Medtr6g043460.1) were downregulated at day 1 of treatment. However, the metabolic pathway and specific functions of these genes remain unknown because of the limited annotation information.

Flavonoids Metabolism-Related Genes Were Induced

The pathway of flavonoid biosynthesis and metabolism is well understood (Winkel-Shirley, 2002). Flavonoids participate in ROS elimination as electron donors (Hernández et al., 2009). In addition, flavonoids bind metal ions and thus reduce the damage caused by metal ions (Pourcel et al., 2006). Flavonoids are also signal molecules that interact with transcription factors, kinases, and histones to regulate gene expression at different levels (Ithal and Reddy, 2004).

In transcriptome data, the expression of several rate-limiting enzymes and modified genes in the flavonoid metabolic pathway was significantly altered. The key enzymes of the flavonoid biosynthesis pathway, including the flavonol synthase c68976_g2_i1 (similar to Medtr5g055690.1) and the chalcone synthase c69319_g1_i1 (similar to Glyma.08g109400.1), as well as the pathway intermediate enzyme 4-p-coumaric acyl CoA ligase c80508_g1_i1 (similar to Medtr1g073180.1) were upregulated under stress. p-Coumaric acyl CoA is the substrate in the synthesis of chalcone, which limits the rate of flavonoid biosynthesis. Furthermore, the flavonoid glycosyltransferase c71058_g2_i1 (similar to Medtr8g083290.1) and the 2-hydroxy isoflavone dehydratase c69844_g1_i1 (similar to Medtr1g105020.1) were similarly upregulated. Glycosylation increases the water solubility of flavonoids and induces other modifications such as acylation (Schwinn et al., 1997). Acylation enhances the stability of flavonoids in the cytoplasm. A previous study confirmed the presence of acylated flavonoids in the chloroplast, where they prevent peroxidation of the thylakoid membrane on the oxidizing side of the Photosystem II (PSII). Isoflavone, a substrate of various flavonoids, is formed by dehydration of 2-hydroxy isoflavone in a reaction catalyzed by 2-hydroxy isoflavone dehydratase. The upregulation of 2-hydroxy isoflavone dehydratase may suggest a high synthesis rate of isoflavone (Hernández et al., 2009).

Downregulation of the Light Pathway

Plants reduce the photosynthesis rate and delay their growth under both salinity stress (Zhang and Shi, 2013) and alkaline stress (Li et al., 2010). A previous study showed that the synthesis of the photosynthetic pigments is blocked by abiotic stressors such as salinity (Zhang and Shi, 2013) and freezing temperatures. In the present study, the light-harvesting complex I chlorophyll a/b-binding protein c73857_g2_i1 (similar to Medtr6g011880.1) was significantly downregulated in response to saline-alkaline stress, suggesting an outstanding inhibition of the plant light-harvesting system. Moreover, the light-harvesting proteins c66633_g1_i1 (similar to Medtr3g101670.1) and c73420_g1_i1 (similar to Medtr4g099340.1) were downregulated, suggesting the inhibition of the photosynthetic system. Nevertheless, a previous study on the transcriptome analysis of Eleusine coracana subjected to salinity stress reported the low efficiency of PSII, similar to the behavior of the chlorophyll a/b binding protein and the light-harvesting proteins under alkaline-saline stress conditions (Rahman et al., 2014).

The expression level of ribulose bisphosphate (RuBP) carboxylase (similar to Medtr7g007120.1), a key enzyme of carbon assimilation in the light pathway, declined and remained low under stress. RuBP carboxylase catalyzes the carboxylation of RuBP, initiating the process of carbon assimilation. The downregulation of RuBP carboxylase indicates that the light pathway is inhibited under saline-alkaline stress and may reduce the generation of ROS. Charge separation of the photon-excited PSII reaction center Chl (P680) is suppressed, leading to the accumulation of excited triplet P680, which reacts with O₂ producing singlet O₂ (¹O₂) (Gong et al., 2013).

Differential Expression of Lipid Metabolism

The lipid metabolism subcategory mainly clustered in GO categories related to synthesis and modification of the calcium signal pathway in the short-term treatment. However, this cluster also included the synthesis and transport of lipids in addition to the calcium signal pathway in the long-term treatment.

The calcium signal pathway plays an important role in the biotic and abiotic stress response (Xiong et al., 2002). Various calcium-dependent stress response proteins are regulated by this pathway as well as the transcription factors are regulated by calcium and calcium binding promoters. The myo-Inositol 1-phosphate synthase c52088_g1_i1 (similar to Medtr3g087590.2), which is a component of the calcium signal pathway, was upregulated, increasing the expression levels four-fold in the long-term treatment compared with those in the control and short-term treatment. Similarly, other components of the calcium pathway, such as Medtr2g095220.1, Medtr7g038700.2, and Glyma.08G109400.1.p maintained high expression levels in the long-term treatment, suggesting that the calcium signal pathway remains active during continuous stress. In a similar study on salt stress, phosphate inositol was upregulated significantly when treated with NaCl for 7 days (Postnikova et al., 2013).

Differentially Expressed Transcription Factors in Transcriptome

Transcription factors regulate the expression of genes involved in plant cell metabolism and play an important role in response to environmental stress, including stress response pathways, growth processes, morphogenesis, and responses to various stimuli. De novo assembly data contained 2,375 transcription factors homologous to M. truncatula. Moreover, the expression levels of 150 transcription factors (6.31% of all transcription factors) were significantly modified (log₂ > 1.5) under stress (Figure 5). Of these, 109 transcription factors responded to saline-alkaline stress at day 1 of treatment and 96 at day 7 of treatment, including the major transcription factor families that respond to abiotic stress such as MYB, WRKY, NAC, AP2/EREBP, bHLH, and bZIP (Table 3; Supplementary Table S5).

FIGURE 5

FIGURE 5. Transcription factors differentially expressed in samples after 1 day and 7 days of saline-alkaline stress.

TABLE 3

TABLE 3. Transcription factors differentially expressed in samples after 1 day and 7 days of saline-alkaline stress.

The MYB transcription factor family presented the most significant changes under continuous stress. Of 179 MYB transcription factors that responded to saline-alkaline stress, 10 MYB were altered significantly at day 1 of treatment and 11 at day 7 of treatment. Although the majority of MYB family transcription factors showed stable upregulation trends in the long-term treatment, a few genes, including c70259_g1_i1 (similar to Medtr3g110028.2) and c30131_g1_i1 (similar to Medtr5g010650.1), were downregulated. In previous studies, MYBs were strongly upregulated in response to salinity in Arabidopsis thaliana (Xu et al., 2015) and Glycine max (Liao et al., 2008). In the transcriptome analysis of alfalfa, MYB was the secondary transcription factor family with the strongest response to salinity stress after AP2 (Postnikova et al., 2013). The expression of MYBs in alfalfa induced by saline-alkaline stress appeared similar to that under salinity stress, suggesting the important role of MYB in alfalfa saline-alkaline tolerance.

Of the 87 differentially expressed WRKY family transcription factors, 16 exhibited significant differential expression in the short-term treatment: 15 were upregulated and only c73314_g1_i1 (similar to Medtr7g109600.1) was downregulated. Among the transcription factors in the long-term treatment, c62895_g1_i1 (similar to Medtr3g095080.1), c67877_g1_i1 (similar to Medtr8g092010.1), c70411_g1_i1 (similar to Medtr4g132430.1), and c81562_g1_i4 (similar to Medtr5g094430.1) were moderately downregulated, except for c81562_g1_i4 that was upregulated. Among them, both WRKY33-like and WRKY75-like transcription factors were highly induced under saline-alkaline stress. A previous study revealed that WRKY33 is highly induced by salt stress and closely related to ABA sensitivity (Jiang and Deyholos, 2009), increasing Arabidopsis NaCl tolerance. WRKY75 was reported to be involved in regulating a nutrient starvation response and root development (Devaiah et al., 2007) and, through its co-expression with a GSTU family member, to limit plant growth and development as well as the level of GSTU expression. Additionally, Rushton et al. (2010) reported that the WRKY family of transcription factors was induced by salinity stress in Arabidopsis thaliana, confirming the important role of WRKY gene family in biological processes and responses to abiotic stress.

The NAC family comprises plant-specific transcription factors is one of the largest transcription factor families in plants involved in responses to many biotic and abiotic stresses (Nakashima et al., 2007, 2012). The members of the NAC transcription factor family showed similar expression in both treatments, whereas the expression of 11 NAC members was significantly modified in both treatments. An AtNAC47-like gene was highly upregulated under saline-alkaline stress, but downregulated under salt, drought, and ABA treatments. The high expression levels of NAC members indicated the important role of NAC in plant saline-alkaline stress tolerance.

The bHLH family of transcription factors was strongly induced by saline-alkaline stress. c77005_g1_i1 (similar to Medtr5g030770.1) was continuously upregulated until day 7 of treatment. The bHLHs are involved in the construction of the light pathway, in flavonoid biosynthesis, and specific root morphology.

The primary family of transcription factors was highly upregulated in the short-term treatment and moderately downregulated in the long-term treatment. Some transcription factors showed successive upregulation under stress, suggesting that these DETFs play an important role in plant tolerance to saline-alkaline stress.

Non-annotated Genes

A total of 19,465 genes (36.14% of all unigenes) with complete open reading frame (ORF) were not annotated in the GenBank, and their expression levels were altered under stress. For instance, the putative protein Medtr1g100627 was upregulated 7.49-fold in the long-term treatment compared with the control, whereas the putative protein Medtr2g075340 was downregulated 7.53-fold in the short-term treatment. Several putative transmembrane genes, such as Medtr8g069283, Medtr2g022540, Medtr2g075940, Medtr1g014110, and Medtr5g095840, were altered significantly under stress, suggesting that plant saline-alkaline stress tolerance includes many unknown mechanisms.

Specifically Expressed Genes under Saline-Alkaline Stress Compared with Saline Stress

The qRT-PCR results showed that the expression of eight genes significantly changed under saline and saline-alkaline stress for 7 days (Figure 6). These differences might suggest that plant tolerance to saline-alkaline environment involved a different mechanism from that in the tolerance to saline stress (Supplementary Table S6).

FIGURE 6

FIGURE 6. Expression pattern of genes in response to saline-alkaline stress and saline stress after 7 days of treatment. a, inositol transporter; b, DNA binding protein; c, raffinose synthase; d, ferritin; e, ACC synthetase; f, NAC transcription factor; g, glutathione S-transferase; h, xyloglucan endotransglucosylase.

A group of genes showed higher expression in response to alkaline-saline stress than to salt stress. For instance, xyloglucan endotransglucosylase/hydrolase (similar to Medtr8g102220.1) was not enhanced under saline stress, but it was increased 2.08-fold by saline-alkaline stress, suggesting that the synthesis of the cell wall component was under stress.

Some genes were induced significantly by saline stress, but showed only a mild response to saline-alkaline stress under the same concentration of sodium. Raffinose synthase (Medtr4g115330.1) was downregulated 2.08-fold under saline stress, whereas the downregulation was –0.11 under saline-alkaline stress. The lower downregulation of raffinose synthase indicated that the raffinose synthase pathway is involved in the response to saline-alkaline stress in alfalfa. Similar regulation was also observed in the ferritin Medtr7g069980.1, suggesting that ferritin is always maintained at a certain level under saline-alkaline stress in the root because of its important role in the regulation of intracellular ion balance (Masuda et al., 2013).

Additional transcription factors responded differently under the two stresses. For instance, a member of the AP2 (Medtr4g100650.1) transcription factor family and an ERF transcription factor (Medtr7g085220.1) were downregulated under saline stress, but upregulated under saline-alkaline stress.

Physiology Assay under Saline-Alkaline Stress

The level of ROS is indicative of the physiological conditions in plants under saline-alkaline stress. The accumulation of ${\text{O}}_2^{-}$ showed a mild increase at day 1 of treatment, but it was markedly increased at day 7 of treatment. In contrast, the level of H₂O₂ gradually increased from day 0 to day 7 of treatment (Figures 7A,B). Furthermore, the activity of POD, which is strongly induced by ROS, was highly induced by saline-alkaline stress at day 1 of treatment and had a slightly higher activity at day 7 of treatment (Figure 7D). The content of MDA reveals the level of lipid peroxidation and indicates the damage to the membrane under stress. The accumulation of MDA had no marked change between day 0 and day 1 of treatment, but it increased at day 7 of treatment (Figure 7C).

FIGURE 7

FIGURE 7. Physiology assay of alfalfa under saline-alkaline stress. (A) content of ${\text{O}}_2^{-}$; (B) content of H₂O₂; (C) content of MDA; (D) activity of POD; (E) content of chlorophyll.

The level of chlorophyll decreased under saline-alkaline stress. The total chlorophyll continued degrading (Figure 7E), suggesting a dramatic inhibition of chlorophyll synthesis due to saline-alkaline stress.

Changes in AKR activity were determined after exposure to saline and saline-alkaline stress for 7 days. The activity of AKR was markedly different under both stresses, and it was higher under saline-alkaline stress than under saline stress (Figure 8).

FIGURE 8

FIGURE 8. Activity of AKR in root part under saline stress and saline-alkaline stress respectively.

Discussion

Transcriptome analysis showed that the component of the ROS scavenging pathway was dramatically upregulated under saline-alkaline stress conditions. Flavonoids, which are potential non-enzymatic protective substances against reactive oxygen (Hernández et al., 2009), have a strong response to sodic stress. In addition, the downregulation of RuBP carboxylase suggested that the photosynthesis pathway was inhibited under saline-alkaline stress (Li et al., 2010; Bu et al., 2012). These results indicated that the antioxidant capacity is one of the major mechanisms of tolerance to saline-alkaline stress in alfalfa. Furthermore, other potential stress response pathways may exist. Eight different expression genes were screened from transcriptome for their dramatically different or even opposite expression under the saline and saline-alkaline stress conditions. A member of NAC transcription factors and AKR were downregulated under salt stress, but upregulated under saline-alkaline stress.

Reactive oxygen species accumulate in plant cells during biotic and abiotic stress conditions. Saline-alkaline stress with high pH exacerbate the accumulation of ROS (Li et al., 2010; Gong et al., 2014), such as superoxide anions, H₂O₂, hydroxyl radicals, and singlet oxygen, leading to oxidative stress that can damage cellular components such as DNA, lipids, proteins, and sugars (Pandhair and Sekhon, 2006). ROS homeostasis in plants is regulated by a complex machinery of enzymatic and non-enzymatic antioxidants through the enzymatic pathway composed of PODs (Sugimoto et al., 2014). The accumulation of MDA was not changed markedly between day 0 and day 1 of treatment, but it dramatically increased after 7 days of treatment, suggesting an oxidative damage to the plasma membrane. The increasing level of ROS and their effect on membrane indicated that plant tissue was exposed to oxidative stress during the entire treatment period.

Flavonoids are a large family of plant secondary metabolites, which contain strong antioxidants (Pourcel et al., 2006) such as ascorbate (vitamin C) and α-tocopherol (vitamin E). It has long been recognized that salt stress induces flavonoid accumulation (Chalker-Scott, 1999; Winkel-Shirley, 2002). Intracellular flavonoids located in the plasma membrane, chloroplast, vacuole, and nuclei are strongly induced by oxidative stress (Hernández et al., 2009). The expression of flavonoid pathway genes is closely related to the content of intracellular ROS. The substantive activity of flavonoid metabolism suggests that flavonoids play an important role in the saline-alkaline stress tolerance mechanism. In addition, the inhibition of the light pathway was observed.

A decrease in the chlorophyll content suggests a decline of the photosynthetic rate, whereas the downregulation of the light-harvesting protein and RuBP carboxylase suggests the inhibition of the plant photosynthetic rate (Gong et al., 2013). A decline in the plant photosynthetic rate may decrease the accumulation of endogenous ROS and reduce the pressure on the POD system caused by excessive intracellular ROS. Moreover, the inhibition of the plant growth rate and developmental processes can improve plant tolerance to saline-alkaline stress.

Peroxidases are widespread in bacteria, fungi, animals, and plants. They are divided into three categories: class I PODs, intracellular PODs that remove intracellular ROS; class II POD, mainly exist in fungi; and class III PODs, extracellular PODs involved in auxin metabolism, cell elongation, lignin metabolism, and resistance to pathogens (Guida et al., 2014). The activity of POD was highly induced at day 1 of treatment and remained constant until day 7 of treatment, revealing the important role of PODs in response to long-term saline-alkaline stress. Class III PODs show higher expression than other PODs; they participate in ROS elimination when seedlings are exposed to long-term saline-alkaline stress. Previous studies characterized the function of class III PODs, and revealed that they are involved in many biological processes (Almagro et al., 2009), including cell elongation, stress defense, seed germination, and especially polymerization of lignin (Shigeto et al., 2015). The differential response of various species of PODs under saline-alkaline conditions suggests that PODs related to intercellular interaction are more important than those respond to ROS under long-term saline-alkaline stress and allow the resuming of plant growth.

In the present study, alfalfa employed different response mechanisms under saline and saline-alkaline stress at the same concentration of sodium ion (0.1 M) and showed a more dramatic change in common upregulated response genes and milder effect in downregulated genes. Genes involved in stress response pathways showed different responses under the two stresses.

The overexpression of xyloglucan endotransglucosylase/hydrolase promotes the cell wall development and induces cell proliferation and differentiation in the seed germination stage (Miedes et al., 2013). The high pH that accompanies saline-alkaline stress may directly damage the cell wall and the membrane; this damage may be reduced by the high expression of xyloglucan endotransglucosylase/hydrolase.

Aldo/keto reductase is a large family involved in plant stress defense. Most AKRs act by reducing aldehydes and ketones to their respective alcohols using NADPH, or rarely NADH, as a cofactor (Penning, 2004). A previous study indicated that AKRs are downregulated under 100-mM saline stress (Postnikova et al., 2013), but upregulated in response to saline-alkaline stress. The higher content of AKR could rapidly eliminate aldehydes and ketones that were generated by cellular metabolic disorders and high pH. These results indicate that detoxification and anti-oxidation play an important role in alfalfa response to saline-alkaline stress, similar to most extreme environmental factors such as heat, cold, and high level of oxidation stress.

Glutathione S-transferase presented the opposite trend in the two different treatments; it was downregulated under saline stress by −1.07, but upregulated under saline-alkaline stress by 4.07. Glutathione S-transferase plays an important role in the antioxidant pathway; the high expression of glutathione S-transferase is indicative of the high levels of active oxygen in the plant root. A similar trend was detected in the NAC transcription factor (similar to Medtr7g085220.1); it was downregulated under saline stress and upregulated under saline-alkaline stress. NAC is a plant specific transcription factor family with diverse roles in development, stress regulation, and response to biotic and abiotic stresses. Previous studies suggested that the NAC transcription factor participates in the regulation of POD synthesis and may be co-expressed with flavonol synthase (Nakashima et al., 2007). The specific response of NAC may play a role in the cellulose synthase and GST pathway regulation. Future studies on NAC identification may uncover a novel mechanism of plant saline-alkaline stress response. Additionally, DEGs with different expression in a variety of stress conditions may reveal the core mechanism and relative pathways of plant saline-alkaline tolerance. Moreover, further studies will help elucidate how to improve the saline-alkaline stress tolerance of sensitive alfalfa cultivars.

Author Contributions

C-HG and Y-MA designed the study; Y-MA, L-LS, and Y-RL performed the experiment; Y-JS supplied bioinformatics analysis platform; Y-MA and Y-JS performed the data analysis; C-HG provided scientific expertise; Y-MA wrote the manuscript.

Funding

MOST 863 project (2013AA102607-5) C-HG. Key Scientific and Technological Project of Heilongjiang Province of China (GA15B105-1) C-HG. Natural and Science Foundation of China (NO.31470571) C-HG. Graduate Innovation Fund of Harbin Normal University (HSDBSCX2015-10) Y-MA. Research supported by the National Major Project for Cultivation of Transgenic Crops (2016ZX08004-002-003) C-HG.

Conflict of Interest Statement

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.

Supplementary Material

The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fpls.2016.00931

Footnotes

1. ^wego.genomics.org.cn/cgi-bin/wego/index.pl
2. ^http://bioinfo.bti.cornell.edu/cgi-bin/itak/index.cgi

References

Almagro, L., Gomez-Ros, L. V., Belchi-Navarro, S., Bru, R., Barcelo, A. R., and Pedreno, M. A. (2009). Class III peroxidases in plant defence reactions. J. Exp. Bot. 60, 377–390. doi: 10.1093/jxb/ern277

PubMed Abstract | CrossRef Full Text | Google Scholar

Bu, N., Li, X. M., Li, Y. Y., Ma, C. Y., Ma, L. J., and Zhang, C. (2012). Effects of Na2CO3 stress on photosynthesis and antioxidative enzymes in endophyte infected and non-infected rice. Ecotoxicol. Environ. Saf. 78, 35–40. doi: 10.1016/j.ecoenv.2011.11.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Chalker-Scott, L. (1999). Environmental significance of anthocyanins in plant stress responses. Photochem. Photobiol. 70, 1–9. doi: 10.1111/j.1751-1097.1999.tb01944.x

CrossRef Full Text | Google Scholar

Devaiah, B. N., Karthikeyan, A. S., and Raghothama, K. G. (2007). WRKY75 transcription factor is a modulator of phosphate acquisition and root development in Arabidopsis. Plant Physiol. 143, 1789–1801. doi: 10.1104/pp.106.093971

PubMed Abstract | CrossRef Full Text | Google Scholar

Ge, Y., Li, Y., Zhu, Y. M., Bai, X., Lv, D. K., Guo, D., et al. (2010). Global transcriptome profiling of wild soybean (Glycine soja) roots under NaHCO3 treatment. BMC Plant Biol. 10:153. doi: 10.1186/1471-2229-10-153

PubMed Abstract | CrossRef Full Text | Google Scholar

Gong, P. B. (2001). In: Principles and Techniques of Plant Physiological Biochemical Experiment. Beijing: Higher Education Press, 121–158.

Gong, B., Li, X., Bloszies, S., Wen, D., Sun, S., Wei, M., et al. (2014). Sodic alkaline stress mitigation by interaction of nitric oxide and polyamines involves antioxidants and physiological strategies in Solanum lycopersicum. Free Radic. Biol. Med. 71, 36–48. doi: 10.1016/j.freeradbiomed.2014.02.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Gong, B., Wen, D., VandenLangenberg, K., Wei, M., Yang, F., Shi, Q., et al. (2013). Comparative effects of NaCl and NaHCO3 stress on photosynthetic parameters, nutrient metabolism, and the antioxidant system in tomato leaves. Sci. Hortic. 157, 1–12. doi: 10.1016/j.scienta.2013.03.032

CrossRef Full Text | Google Scholar

Grabherr, G. M., Haas, B. J., Yassour, M., Levin, J. Z., Thompson, D. A., Amit, I., et al. (2011). Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 29, 644–652.

Google Scholar

Guida, V., Cantarella, M., Chambery, A., Mezzacapo, M. C., Parente, A., Landi, N., et al. (2014). Purification and characterization of novel cationic peroxidases from Asparagus acutifolius L. with biotechnological applications. Mol. Biotechnol. 56, 738–746. doi: 10.1007/s12033-014-9752-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Hernández, I., Alegre, L., Breusegem, F. V., and Munne-Bosch, S. (2009). How relevant are flavonoids as antioxidants in plants. Trends Plant Sci. 14, 125–132. doi: 10.1016/j.tplants.2008.12.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Hodges, D. M., DeLong, J. M., Forney, C. F., and Prange, R. K. (1999). Improving the thiobarbituric acid-reactive-substances assay for estimating lipid peroxidation in plant tissues containing anthocyanin and other interfering compounds. Planta 207, 604–611. doi: 10.1007/s004250050524

CrossRef Full Text | Google Scholar

Irsigler, A. S. T., Costa, M. D. L., and Fontes, E. P. B. (2007). Expression profiling on soybean leaves reveals integration of ER and osmotic-stress pathways. BMC Genomics 8:431. doi: 10.1186/1471-2164-8-431

PubMed Abstract | CrossRef Full Text | Google Scholar

Ithal, N., and Reddy, A. R. (2004). Rice flavonoid pathway genes, OsDfr and OsAns, are induced by dehydration, high salt and, A. B. A., and contain stress responsive promoter elements that interact with the transcription activator, OsC1-MYB. Plant Sci. 166, 1505–1513.

Google Scholar

Jiang, Y. Q., and Deyholos, M. K. (2006). Comprehensive transcriptional profiling of NaCl-stressed Arabidopsis roots reveals novel classes of responsive genes. BMC Plant Biol. 6:25. doi: 10.1186/1471-2229-6-25

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, Y. Q., and Deyholos, M. K. (2009). Functional characterization of Arabidopsis NaCl-inducible WRKY25 and WRKY33 transcription factors in abiotic stresses. Plant Mol. Biol. 69, 91–105. doi: 10.1007/s11103-008-9408-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Kawasaki, S., Borchert, C., Deyholos, M., Wang, H., Brazille, S., Kawai, K., et al. (2001). Gene expression profiles during the initial phase of salt stress in rice. Plant Cell 13, 889–905. doi: 10.1105/tpc.13.4.889

PubMed Abstract | CrossRef Full Text | Google Scholar

Kawaura, K., Mochida, K., Yamazaki, Y., and Ogihara, Y. (2006). Transcriptome analysis of salinity stress responses in common wheat using a 22k oligo-DNA microarray. Funct. Integr. Genomics 6, 132–142. doi: 10.1007/s10142-005-0010-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Langmead, B., Trapnell, C., Pop, M., and Salzberg, S. L. (2009). Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 10:R25. doi: 10.1186/gb-2009-10-3-r25

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, H., Luo, X. Y., and Wang, D. K. (1996). Report on breeding of novel varieties of alfalfa. Heilongjiang J. Anim. Sci. Technol. 1, 3–7.

Li, R. L., Shi, F. C., Kukuda, K., and Yang, Y. L. (2010). Effects of salt and alkali stresses on germination, growth, photosynthesis and ion accumulation in alfalfa (Medicago sativa L.). Soil Sci. Plant Nutr. 56, 725–733. doi: 10.1111/j.1747-0765.2010.00506.x

CrossRef Full Text | Google Scholar

Liao, Y., Zou, H. F., and Chen, S. Y. (2008). Soybean GmMYB76, GmMYB92, and GmMYB177 genes confer stress tolerance in transgenic Arabidopsis plants. Cell Res. 18, 1047–1060. doi: 10.1038/cr.2008.280

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, F., and Pang, S. J. (2010). Stress tolerance and antioxidant enzymatic activities in the metabolism of the reactive oxygen species in two intertidal red algae grateloupia turuturu and palmaria palmate. J. Exp. Mar. Biol. Ecol. 328, 82–87. doi: 10.1016/j.jembe.2009.11.005

CrossRef Full Text | Google Scholar

Masuda, H., Kobayashi, T., Ishimaru, Y., Takahashi, M., Aung, M. S., Nakanishi, H., et al. (2013). Iron-biofortification in rice by the introduction of three barley genes participated in mugineic acid biosynthesis with soybean ferritin gene. Front. Plant Sci. 4:132. doi: 10.3389/fpls.2013.00132

PubMed Abstract | CrossRef Full Text | Google Scholar

Miedes, E., Suslov, D., Vandenbussche, F., Kenobi, K., Ivakov, A., Van Der Straeten, D., et al. (2013). Xyloglucan endotransglucosylase/hydrolase (XTH) overexpression affects growth and cell wall mechanics in etiolated Arabidopsis hypocotyls. J. Exp. Bot. 64, 2481–2497. doi: 10.1093/jxb/ert107

PubMed Abstract | CrossRef Full Text | Google Scholar

Mukherjee, S. P., and Choudhuri, M. A. (1983). Implications of water stress-induced changes in the levels of endogenous ascorbic acid and hydrogen peroxide in vigna seedlings. Physiol. Plant. 58, 166–170. doi: 10.1111/j.1399-3054.1983.tb04162.x

CrossRef Full Text | Google Scholar

Nakashima, K., Takasaki, H., Mizoi, J., Shinozaki, K., and Yamaguchi-Shinozaki, K. (2012). NAC transcription factors in plant abiotic stress responses. Biochim. Biophys. Acta 1819, 97–103. doi: 10.1016/j.bbagrm.2011.10.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Nakashima, K., Tran, L. S., Van Nguyen, D., Fujita, M., Maruyama, K., Todaka, D. et al. (2007). Functional analysis of a NAC-type transcription factor OsNAC6 involved in abiotic and biotic stress-responsive gene expression in rice. Plant J. 51, 617–630. doi: 10.1111/j.1365-313X.2007.03168.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Oliveros, J. C. (2007–2015). Venny. An Interactive Tool for Comparing Lists with Venn’s Diagrams. Available at: http://bioinfogp.cnb.csic.es/tools/venny/index.html

Google Scholar

Pandhair, V., and Sekhon, B. S. (2006). Reactive oxygen species and antioxidants in plants. J. Plant Biochem. Biotechnol. 15, 71–78. doi: 10.1007/BF03321907

CrossRef Full Text | Google Scholar

Peng, Y. L., Gao, Z. W., Gao, Y., Liu, G. F., Sheng, L. X., Wang, D. L., et al. (2008). Eco-physiological characteristics of alfalfa seedlings in response to various mixed salt-alkaline stresses. J. Integr. Plant Biol. 50, 29–39. doi: 10.1111/j.1744-7909.2007.00607.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Penning, T. M. (2004). Aldo-keto reductases and formation of polycyclic aromatic hydrocar-bon o-quinones. Methods Enzymol. 378, 31–67. doi: 10.1016/S0076-6879(04)78003-9

CrossRef Full Text | Google Scholar

Postnikova, O. A., Shao, J. T., and Nemchinov, L. G. (2013). Analysis of the alfalfa root transcriptome in response to salinity stress. Plant Cell Physiol. 54, 1041–1055. doi: 10.1093/pcp/pct056

PubMed Abstract | CrossRef Full Text | Google Scholar

Pourcel, L., Routaboul, J. M., Cheynier, V., Lepiniec, L., and Debeaujon, I. (2006). Flavonoid oxidation in plants: from biochemical properties to physiological functions. Trends Plant Sci. 12, 29–36. doi: 10.1016/j.tplants.2006.11.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Rahman, H., Jagadeeshselvam, N., Valarmathi, B., Sachin, B., Sasikala, R., Senthil, N., et al. (2014). Transcriptome analysis of salinity responsiveness in contrasting genotypes of finger millet (Eleusine coracana L.) through RNA-sequencing. Plant Mol. Biol. 85, 485–530. doi: 10.1007/s11103-014-0199-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Rushton, P. J., Somssich, I. E., Ringler, P., and Shen, Q. J. (2010). WRKY transcription factors. Plant Sci. 15, 247–258. doi: 10.1016/j.tplants.2010.02.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Schwinn, K. E., Davies, K. M., Deroles, S. C., Markham, K. R., Miller, R. M., Marie Bradley, J., et al. (1997). Expression of an Antirrhinum majus UDP-glucose: flavonoid-3-O-glucosyltransferase transgene alters flavonoid glycosylation and acylation in lisianthus (Eustoma grandiflorum Grise.). Plant Sci. 125, 53–61. doi: 10.1016/S0168-9452(97)04624-4

CrossRef Full Text | Google Scholar

Shigeto, J., Itoh, Y., Hirao, S., Ohira, K., Fujita, K., and Tsutsumi, Y. (2015). Simultaneously disrupting AtPrx2, AtPrx25 and AtPrx71 alters lignin content and structure in Arabidopsis stem. J. Integr. Plant Biol. 57, 349–356. doi: 10.1111/jipb.12334

PubMed Abstract | CrossRef Full Text | Google Scholar

Shu, Y. J., Zhang, J., Ao, Y., Song, L. L., and Guo, C. H. (2015). Analysis of the Thinopyrum elongatum transcriptome under water deficit stress. Int. J. Genomics 2015:8.

PubMed Abstract | Google Scholar

Sugimoto, M., Oono, Y., Gusev, O., Matsumoto, T., Yazawa, T., Levinskikh, M. A., et al. (2014). Genome-wide expression analysis of reactive oxygen species gene network in Mizuna plants grown in long-term spaceflight. BMC Plant Biol. 14:4. doi: 10.1186/1471-2229-14-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, Y. P., Wang, F. W., Wang, N., Dong, Y., Liu, Q., Zhao, L., et al. (2013). Transcriptome exploration in Leymus chinensis under saline-alkaline treatment using 454 pyrosequencing. PLoS ONE 8:e53632. doi: 10.1371/journal.pone.0053632

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang, H., Krishnakumar, V., Bidwell, S., Rosen, B., Chan, A., Zhou, S., et al. (2014). An improved genome release (version Mt4.0) for the model legume Medicago truncatula. BMC Genomics 15:312. doi: 10.1186/1471-2164-15-312

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Q., Ye, T., Ma, Z., Chen, R., Xie, T., and Yin, X. (2014). Characterization and site-directed mutation of a novel aldo-keto reductase from Lodderomyces elongisporus NRRL YB-4239 with high production rate of ethyl (R)-4-chloro-3-hydroxybutanoate. J. Ind. Microbiol. Biotechnol. 41, 1609–1616. doi: 10.1007/s10295-014-1502-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, C., Gao, C. Q., Wang, L. Q., Zheng, L., Yang, C. P., and Wang, Y. C. (2014). Comprehensive transcriptional profiling of NaHCO3-stressed Tamarix hispida roots reveals networks of responsive genes. Plant Mol. Biol. 84, 145–157. doi: 10.1007/s11103-013-0124-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Winkel-Shirley, B. (2002). Biosynthesis of flavonoids and effects of stress. Curr. Opin. Plant Biol. 5, 218–223. doi: 10.1016/S1369-5266(02)00256-X

CrossRef Full Text | Google Scholar

Wong, C. E., Li, Y., and Moffatt, B. A. (2006). Transcriptional profiling implicates novel interactions between abiotic stress and hormonal responses in Thellungiella, a close relative of Arabidopsis. Plant Physiol. 140, 1437–1450. doi: 10.1104/pp.105.070508

PubMed Abstract | CrossRef Full Text | Google Scholar

Xiong, L. M., Schumaker, K. S., and Zhu, J. K. (2002). Cell signaling during cold, drought, and salt stress. Plant Cell 14, S165–S183.

Google Scholar

Xu, D. H., and Tuyen do, D. T. (2012). Genetic studies on saline and sodic tolerances in soybean. Breed. Sci. 61, 559–565. doi: 10.1270/jsbbs.61.559

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, R., Wang, Y., Zheng, H., Lu, W., Wu, C., Huang, J., et al. (2015). Salt-induced transcription factor MYB74 is regulated by the RNAdirected DNA methylation pathway in Arabidopsis. J. Exp. Bot. 66, 5997–6008.

PubMed Abstract | Google Scholar

Yang, Q., Chen, Z. Z., Zhou, X. F., Yin, H. B., Li, X., Xin, X. F., et al. (2009). Overexpression of SOS (salt overly sensitive) genes increases salt tolerance in transgenic Arabidopsis. Mol. Plant 2, 22–31. doi: 10.1093/mp/ssn058

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, Z., and Wang, B. S. (2015). Present status of saline soil resources and countermeasures for improvement and utilization in China. Shandong Agric. Sci. 47, 125–130.

Ye, J., Fang, L., Zheng, H., Zhang, Y., Chen, J., Zhang, Z., et al. (2006). WEGO: a web tool for plotting GO annotations. Nucleic Acids Res. 34, W293–W297. doi: 10.1093/nar/gkl031

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, J. L., and Shi, H. Z. (2013). Physiological and molecular mechanisms of plant salt tolerance. Photosynth. Res. 115, 1–22. doi: 10.1007/s11120-013-9813-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, X., Takano, T., and Liu, S. (2006). Identification of a mitochondrial ATP synthase small subunit gene (RMtATP6) expressed in response to salts and osmotic stresses in rice (Oryza sativa L.). J. Exp. Bot. 57, 193–200. doi: 10.1093/jxb/erj025

PubMed Abstract | CrossRef Full Text | Google Scholar

Zheng, Y., Zhao, L., Gao, J., and Fei, Z. (2011). iAssembler: a package for de novo assembly of Roche-454/Sanger transcriptome sequences. BMC Bioinformatics 12:453. doi: 10.1186/1471-2105-12-453

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, J. K. (2001). Plant salt tolerance. Trends Plant Sci. 6, 66–71. doi: 10.1016/S1360-1385(00)01838-0

CrossRef Full Text | Google Scholar