Effects of drought and salt-stresses on gene expression in Caragana korshinskii seedlings revealed by RNA-seq
- Shaofeng Li†1,
- Chengming Fan†2,
- Yan Li†2,
- Jianhui Zhang2,
- Jingshuang Sun1,
- Yuhong Chen2,
- Changyan Tian3,
- Xiaohua Su4,
- Mengzhu Lu4,
- Chengzhi Liang2 and
- Zanmin Hu2Email author
© Li et al. 2016
Received: 11 August 2015
Accepted: 29 February 2016
Published: 8 March 2016
Drought and soil salinity are major abiotic stresses. The mechanisms of stress tolerance have been studied extensively in model plants. Caragana korshinskii is characterized by high drought and salt tolerance in northwestern China; unique patterns of gene expression allow it to tolerate the stress imposed by dehydration and semi-desert saline soil. There have, however, been no reports on the differences between C. korshinskii and model plants in the mechanisms underlying their drought and salt tolerance and regulation of gene expression.
Three sequencing libraries from drought and salt-treated whole-seedling- plants and the control were sequenced to investigate changes in the C. korshinskii transcriptome in response to drought and salt stresses. Of the 129,451 contigs, 70,662 (54.12 %) were annotated with gene descriptions, gene ontology (GO) terms, and metabolic pathways, with a cut-off E-value of 10−5. These annotations included 56 GO terms, 148 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, and 25 Clusters of Orthologous Groups (COG). On comparison of the transcriptomes of the control, drought- and salt-treated plants, 1630 and 1521 contigs showed significant differences in transcript abundance under drought and salt stresses. Compared to the differentially expressed genes (DEGs) in drought- or salt-treated Arabidopsis in the database, 542 DEGs in drought-treated C. korshinskii and 529 DEGs in salt-treated samples were presumably unique to C. korshinskii. The transcription profiles revealed that genes related to transcription factors, protein kinases, and antioxidant enzymes are relevant to the tolerance of drought and salt stress in this species. The expression patterns of 38 randomly selected DEGs were confirmed by quantitative real-time PCR and were essentially consistent with the changes in transcript abundance identified by RNA-seq.
The present study identified potential genes involved in drought and salt tolerance in C. korshinskii, as well as many DEGs uniquely expressed in drought- or salt-treated C. korshinskii samples compared to Arabidopsis. To our knowledge, this study is the first exploration of the C. korshinskii transcriptome under drought and salt conditions, and these results will facilitate the discovery of specific stress-resistance-related genes in C. korshinskii, possibly leading to the development of novel plant cultivars through genetic engineering.
Environmental stresses or abiotic constraints such as salinity, drought, heavy metals and low temperature are major obstacles to plant survival and development in many places worldwide. Plant responses and adaptations to abiotic stresses are complex . Many cellular and physiological processes such as signal perception and transduction [2, 3], regulation of transcription, membrane trafficking, energy metabolism , and protein phosphorylation and dephosphorylation are involved in plants’ toleration of excessive salt and drought stresses. Several studies on drought and salt stresses using microarrays, transcriptome sequencing, and microRNA sequencing [5–7] have identified stress-inducible genes involved in water transport (aquaporins), ion transport (plasma membrane (PM) H+-ATPase, Na+/H+ exchanger or Na+/H+ antiporter), cellular membrane integrity (proline, glycine betaine, mannitol), scavenging of free oxygen radicals (superoxide dismutase, catalase and peroxidase), and protecting macromolecules (late embryogenesis abundant proteins and chaperones). Other proteins, such as regulatory proteins (transcription factors, protein kinases, protein phosphatases, and calmodulin-binding proteins), were found to be involved in signal transduction [8, 9]. ABA-dependent and ABA-independent signaling pathways were induced by plants in response to stress . The ABA-dependent transcription factors include MYC/MYB and ABA-responsive element binding/ABA-binding factor (AREB/ABF). The ABA-independent transcription factors include dehydration-responsive element-binding proteins (DREB), C-repeat/drought-responsive element (CRT/DRE), and CRT/DRE-binding factor (CBF) [11, 12]. Other transcription factors responding to abiotic stress conditions are zinc-finger proteins , basic-domain leucine-zipper (bZIP) , WRKY , and NACs .
The development of next-generation sequencing technology (NGS) has led to major progress in understanding plant responses to drought and salt stress [17–19]. Some recent efforts include tissue-selective signaling and hormone crosstalk in response to salt and osmotic stresses and, recently, the dynamics of transcriptional networks that are functional at the seedling stage in response to water stress in Gossypium arboreum . The transcriptional characteristics of genes related to the reactive oxygen species (ROS)-scavenging system are important to the salt tolerance of Reaumuria trigyna . Further, genome-wide gene-expression profiling was used to examine the molecular mechanism and physiological response of Gossypium herbaceum to drought . Thumma et al.  used RNA sequencing (RNA-seq) to study the effect of water stress on gene expression in Eucalyptus camaldulensis seedlings derived from three natural populations. RNA-seq technology and its applications have dramatically accelerated plant-genomics research, including high-throughput sequencing of non-model-plant transcriptomes and large-scale genome-wide expression analysis.
Caragana korshinskii is a leguminous shrub with highly developed root systems and strong stress tolerance [24, 25] that is commonly found in arid and semi-arid lands in northwestern China and Mongolia [26, 27]. Therefore, C. korshinskii is widely planted in China to prevent desertification and improve vegetation coverage. The effect of drought or salt stress on plant gene expression has been intensely studied in numerous species including Arabidopsis , rice , maize , wheat (Triticum aestivum L.) , and Sorghum bicolor . In C. korshinskii, some investigations have focused on growth properties, nutritive characteristics, evapotranspiration, and genetic diversity [31–33], and several drought- and salt-stress-response genes have already been characterized [34–36], but there have been no studies on the stress-signaling pathways. This lack of molecular studies is a significant obstacle to understanding the molecular mechanisms underlying drought or salt adaptation in C. korshinskii.
In this study, we generated transcriptome datasets to explore the salt- and drought-tolerance mechanisms of C. korshinskii using the Illumina HiSeq™ 2000 platform. The main objectives of this study are to identify genes showing transcriptional differences and to identify their putative functions, as well as to identify genes putatively unique to C. korshinskii compared to the model plant Arabidopsis and describe the patterns in their transcript abundance under salt- and drought-treated conditions. The assembled, annotated transcriptome data and differential expression profiles will facilitate further genetic and genomics studies on the molecular mechanisms of salt and drought tolerance in C. korshinskii and most likely in other leguminous plants.
Illumina sequencing and de novo assembly
Summary of Caragana korshinskii transcriptome assembly in this study compared with the study by Long et al. 
Mean length (bp)
Functional annotation and analysis
Among the 129,451 transcripts, the direction could be determined for 70,662. Altogether, 70,491 (54.45 %) contigs were successfully annotated in the nr, GO, KEGG, COG, Swiss-Prot and InterProScan databases. These contigs are listed in Additional file 3. Of these, 70,062 (54.12 %) of 129,451 contigs had significant matches in the nr database, 62,605 (43.36 %) contigs aligned to the KEGG database, 47,462 (36.66 %) contigs matched GO annotations, and 8321 (6.43 %) contigs were similar to proteins in the InterProScan database (Additional file 4). Because of the lack of genome and EST information for C. korshinskii, 45.55 % (58,960) of the contigs did not match any known genes in the database (Additional file 5).
Analysis of differential expression of assembled C. korshinskii transcripts under drought and salt treatments
Venn-diagram analysis revealed the specificity and overlap among differentially expressed genes (Fig. 4). Among 1630 and 1521 differentially expressed genes, numerous genes overlapped in expression in response to the salt and drought treatments. For example, in drought- and salt-treated plants, 148 upregulated transcripts overlapped and 368 downregulated transcripts overlapped. In addition, 1005 of 1521 contigs appeared only in salt-treated samples, while 1111 of 1630 contigs were expressed only in drought-treated samples.
GO annotation and statistical analyses demonstrated that 1844 tested DEGs (some contigs have more than one GO annotation) were classified into three GO ontologies and 60 terms in drought-treated plants. Of these DEGs, 843 were associated with biological processes, 419 with cellular components, and 612 with molecular function. Among these GO categories, “oxidoreductase activity” (GO:0016491,168 DEGs), “oxidation-reduction process” (GO:0055114,146 DEGs), and “response to stress” (GO:0006950, 149 DEGs) were significantly enriched among DEGs compared to the whole transcriptome background. In the molecular-function category, 2640 tested DEGs (some contigs have more than one GO annotation) were enriched in 63 terms in salt-treated samples. The overwhelming majority of the genes with differential expression between the control and salt-treated conditions were also related to “oxidoreductase activity” (GO:0016491,186 DEGs), “oxidation-reduction process” (GO:0055114,169 DEGs), and “response to stress” (GO:0006950,131 DEGs). Approximately 41.7 % of 60 GO terms in drought-treated C. korshinskii were also found in salt-treated C. korshinskii, such as oxidoreductase activity, oxidation-reduction process, electron carrier activity, and monooxygenase activity (Additional file 8).
In our analysis, 1603 tested DEGs from drought-treated samples mapped to 12 pathways (Corrected P-Value <0.05), with most mapping to “metabolic pathways” (28 DEGs), “biosynthesis of secondary metabolites” (18 DEGs), “porphyrin and chlorophyll metabolism” (6 DEGs), “starch and sucrose metabolism” (6 DEGs), and “cysteine and methionine metabolism” (5 DEGs). In salt-treated samples, 1011 DEGs were enriched in 15 metabolic pathways. The most-enriched metabolic pathways were “ribosome” (53 DEGs), “metabolic pathways” (25 DEGs), “biosynthesis of secondary metabolites” (17 DEGs), “glycolysis/gluconeogenesis” (10 DEGs), and “phenylpropanoid biosynthesis” (10 DEGs). Among these metabolic pathways, starch and sucrose metabolism, biosynthesis of secondary metabolites, metabolic pathways, pentose and glucuronate interconversions, linoleic-acid metabolism, and the ribosome were enriched in both drought- and salt-treated C. korshinskii (Additional file 10).
Functional genes related to protein kinase
We further identified 289 contigs responsible for mitogen-activated protein kinases, of which two contigs (Cko_contig_220 and Cko_contig_10124) were upregulated dramatically by drought or salt stress (Additional file 14). In addition, eight G-type lectin S-receptor-like serine/threonine-protein kinases were identified. Among these kinases, three (Cko_contig_99573, Cko_contig_39198, and Cko_contig_36276) were upregulated more than 32-fold under drought stress (Additional file 13). We also found that 12 cysteine-rich repeat receptor-like protein kinase (CRK) genes showed differential expression (Additional file 13).
Transcription factors responding to drought and salt treatments in C. korshinskii
Genes of unknown function and genes unique to C. korshinskii
In this study, many contigs had no homologous matches among database sequences and thus have no known function. The vast set of un-annotated contigs (45.55 %, 58,960 of 129,451) could be novel genes not present in the known genomes and specific to C. korshinskii (Additional file 5). Among those unknown contigs, 124 (7.61 %) of 1630 DEGs and 114 (7.50 %) of 1521 DEGs were drought- or salt-responsive, respectively, in C. korshinskii (Additional file 7).
To identify genes that might be expressed uniquely in C. korshinskii under salt and/or drought stresses, we performed a new set of BLAST searches to compare these contigs to genes expressed in Arabidopsis under salt and/or drought stresses .
A total of 295 DEGs that were downregulated in drought-treated C. korshinskii were also determined to be downregulated in drought-treated Arabidopsis (Fig. 9, Additional file 16). Some of these genes encode nitrate reductase, glutathione S-transferase, etc.; however, some genes (29 DEGs) that were strongly downregulated in drought-treated C. korshinskii did not significantly change in expression in drought-treated Arabidopsis (Additional file 16), such as peroxidase and thiol protease. Many genes (249 DEGs) that were significantly downregulated in drought-treated C. korshinskii were upregulated in drought-treated Arabidopsis (Fig. 9, Additional file 16), including genes encoding subtilisin-like protease, ubiquitin carboxyl-terminal hydrolase, etc.
A total of 188 DEGs were upregulated in salt-treated Arabidopsis and C. korshinskii, such as spore-specific catalase and aldehyde dehydrogenase; however, many genes (17 DEGs) were strongly upregulated in salt-treated C. korshinskii but remained nearly constant in expression in salt-treated Arabidopsis (Fig. 9, Additional file 16), such as plasma membrane H+-ATPase and 60S ribosomal protein. Some genes (210 DEGs) that were significantly upregulated in salt-treated C. korshinskii were slightly downregulated in salt-treated Arabidopsis (Additional file 16), including genes encoding ADP, ATP-carrier protein, and mannitol dehydrogenase.
A total of 241 DEGs were downregulated in salt-treated C. korshinskii and Arabidopsis (Fig. 9, Additional file 16), including genes encoding LRR receptor-like serine/threonine-protein kinase and beta-glucosidase; however, some genes (33 DEGs) were strongly downregulated in salt-treated C. korshinskii but did not significantly change in expression in salt-treated Arabidopsis (Additional file 16), such as glutathione S-transferase and isoflavone conjugate-specific beta-glucosidase. Many genes (242 DEGs) that were significantly downregulated in salt-treated C. korshinskii were instead upregulated in salt-treated Arabidopsis (Additional file 16), including genes encoding peptide/nitrate transporter and E3 ubiquitin-protein ligase.
Of the 148 DEGs that were upregulated in both drought- and salt-treated C. korshinskii, 51 DEGs were also upregulated in both drought- and salt-treated Arabidopsis and 31 were downregulated in both drought- and salt-treated Arabidopsis. Of the 1630 DEGs upregulated in drought-treated C. korshinskii, 53 DEGs were upregulated and 35 DEGs were downregulated in drought-treated Arabidopsis. In salt-treated Arabidopsis, 53 of these DEGs were upregulated and 33 were downregulated. Six DEGs upregulated in drought- and salt-treated C. korshinskii remained relatively constant in both drought- and salt-treated Arabidopsis (Additional file 17).
Among 368 DEGs that were downregulated in both drought- and salt-treated C. korshinskii, 88 were upregulated in both drought- and salt-treated Arabidopsis and 117 were downregulated in both drought- and salt-treated Arabidopsis. Of the 1521 DEGs downregulated in drought-treated C. korshinskii, 95 were upregulated and 121 downregulated in drought-treated Arabidopsis, while, 93 of these DEGs were upregulated and 120 downregulated in salt-treated Arabidopsis. Twelve DEGs that were downregulated in both drought- and salt-treated C. korshinskii remained relatively constant in both drought- and salt-treated Arabidopsis, except that Cko_contig_45026 remained relatively constant in drought-treated Arabidopsis but not in salt-treated Arabidopsis (Additional file 17).
In drought-treated C. korshinskii, however, 542 DEGs (33.25 %) did not match any Arabidopsis sequences in this study and thus were presumably unique to C. korshinskii. The expression levels of 17 genes changed significantly in drought-treated C. korshinskii; for example, Cko_contig_61538 was predicted to encode E3 ubiquitin-protein ligase and was upregulated 790.9-fold (Additional file 18). In addition, genes for proteins such as phosphoinositide phospholipase C (Cko_contig_42646), trypsin-protein inhibitor (Cko_contig_23762, Cko_contig_27564), and bZIP transcription factor (Cko_contig_23833) were expressed uniquely in drought-treated C. korshinskii.
In addition, 529 DEGs (34.80 %) in salt-treated C. korshinskii had no homologous sequence in the Arabidopsis datasets and were considered unique to C. korshinskii. The expression levels of 27 genes changed dramatically in salt-treated C. korshinskii, such as Cko_contig_58120, predicted to encode polyphenol oxidase, which was upregulated 313.3-fold, and Cko_contig_92720, predicted to encode serine/threonine kinase, which was upregulated 2678.4-fold (Additional file 18). Other genes, such as thiol protease (Cko_contig_105443), UDP-glucosyltransferase (Cko_contig_3136), and polyphenol oxidase (Cko_contig_58120), were found to be uniquely expressed in salt-treated C. korshinskii (Additional file 18). Although we cannot yet be certain which genes are responsible for C. korshinskii’s stress resistance, we are convinced that the majority of predicted unigenes unique to C. korshinskii represent a valuable basis on which to explore gene diversity in C. korshinskii and support drought- or salt-stress resistance in C. korshinskii.
Quantitative real-time-PCR validation of differentially expressed transcripts from RNA-seq
Cko_contig_95695, a homolog of glutathione S-transferase, was upregulated 22.93-fold under drought stress and 8.08-fold under salt stress: 31.31-fold in the drought-treated root, 4.04-fold in the salt-treated root, 13.05-fold in the drought-treated stem, 10.03-fold in the salt-treated stem, 7.08-fold in the drought-treated leaf and 1.14-fold in the salt-treated leaf (Fig. 10b). Nonetheless, moderate discrepancies between the expression levels and RPKM values were observed for three contigs: Cko_contig_12764, Cko_contig_59435, and Cko_contig_27564 (Fig. 10a and b).
In this study, we generated de novo transcriptomes of C. korshinskii under drought and salt treatment and compared the expression profiles between control and drought and salt conditions, as well as between C. korshinskii and the model plant Arabidopsis. From our analysis, we found that many genes/contigs related to protein kinases, certain transcription factors, the ROS-scavenging system and ion transport proteins, gibberellin biosynthesis, carbon metabolism and osmotic regulation, and brassinosteroid-biosynthesis pathways are involved in drought- and salt-stress responses in C. korshinskii, with some of these genes being expressed significantly between C. korshinskii and Arabidopsis, as discussed in detail below. These results provide a molecular basis for understanding the specific drought- and salt-tolerance mechanisms of C. korshinskii.
Construction of an informative transcriptome dataset for C. korshinskii
Transcriptome sequencing of C. korshinskii is fundamental to the identification of functional genes and the elucidation of the molecular mechanisms involved in drought or salt adaptation and tolerance in C. korshinskii. In a recent study in C. korshinskii , a total of 66.35 million raw sequencing reads were generated, 64.03 million clean reads were obtained, and 86,265 contigs were assembled. The percentage of annotation was 48.09 %. In our study, a total of 199.94 million raw reads were generated from the samples, and 182.93 million clean reads were assembled into 129,451 contigs (Table 1). Notably, the percentage of annotation was 54.45 %. The numbers of total raw reads, clean reads, assembled contigs and the percentage of annotation were all higher than those obtained previously , likely because we used several samples under different growth conditions for sequencing analysis, whereas only a sample under normal growth conditions was used by Long et al. . Our de novo assembly of RNA-seq improves the genome annotation of C. korshinskii.
Based on comparisons with different plant species, the top four species with BLAST hits to annotated unigenes were Glycine max, Medicago truncatula, Vitis vinifera, and Lotus japonicus, species whose genome annotations are comprehensive and largely accepted. The low rate of annotation compared with other plants such as Safflower Flowers (58 %) , Aegilops variabilis (66.46 %) , and Ma Bamboo (78.9 %)  might be due to technological limitations, such as sequencing depth, splice junctions and read length , which are common in de novo assembly and RNA-seq analysis. We can speculate that the unannotated contigs (45.55 %, 58,960 of 129,451 contigs) represent a specific genetic resource for C. korshinskii, which is worthy of further investigation. Unannotated sequences were on average much shorter than the annotated contigs (334 bp vs 975 bp); however, a considerable percentage of contigs (20.2 %, or 14,209 of 70,491) that were longer than 500 bp and that had RPKM values  over three had no known homologs in known plant species.
Protein kinases are important for drought and salt tolerance in C. korshinskii
Protein kinases phosphorylate and activate many target proteins in major signaling events in response to osmotic stress in higher plants . Microarray analysis revealed that Glycine soja G-type lectin S-receptor-like serine/threonine-protein kinase GsSRK significantly increased in expression under alkaline (NaHCO3) treatment, identifying it as an alkaline-responsive gene . In this work, Cko_contig_103451, associated with G-type lectin S-receptor-like serine/threonine-protein kinase, was upregulated 7.6-fold in drought-treated C. korshinskii whole plants and 4.6-fold in drought-treated Arabidopsis datasets. We also found that three other homologous genes (Cko_contig_99573, Cko_contig_39198, and Cko_contig_36276) were upregulated more than 32-fold in drought-treated C. korshinskii but changed little in Arabidopsis datasets (Additional file 20). These results are also consistent with the results of studies on the expression pattern of GsSRK in Glycine soja under NaCl (upregulated 4-fold in leaves and 3.5-fold in roots) and PEG treatments (upregulated 3.5-fold in leaves and 3.3-fold in roots) .
Arabidopsis thaliana L-type lectin-domain-containing receptor kinase AtLPK1 is strongly induced by stress hormones and salt treatments, and overexpression of AtLPK1 in Arabidopsis enhances salt-stress tolerance . In this study, Cko_contig_12764 had a highly significant BLAST hit to the L-type lectin-domain-containing receptor kinase and showed increased transcript abundance in salt-treated Arabidopsis datasets (upregulated 3.5-fold), but was downregulated 8-fold in salt-treated C. korshinskii. Another homologous gene (Cko_contig_49124) was upregulated 7-fold in drought-treated C. korshinskii but upregulated slightly in Arabidopsis (Additional file 20). Kempa et al.  found that a plastid-localized glycogen synthase kinase 3 GSK-3 was induced at high concentrations of NaCl, KCl and sorbitol and could modulate salt-stress tolerance and regulate carbohydrate metabolism in Medicago sativa cv. Europa. We identified 58 contigs related to glycogen synthase kinase and only one differentially expressed transcript (Cko_contig_85930), which was uniquely expressed and upregulated 7.7-fold under drought stress in C. korshinskii, with no homologous gene detected in the examined Arabidopsis datasets.
In C. korshinskii, a large number of protein kinases may play an essential role in sensing external drought and salt signals and regulating gene expression at the cellular level in response, as seen in Arabidopsis (Additional file 20); however, the genes encoding protein kinases in response to stresses (e.g., salt or drought) may differ between C. korshinskii and other plants.
Transcription factors involved in drought- and salt-stress response in C. korshinskii
Many WRKY, MYB, NAM, ATAF, and CUC (NAC) transcription factors were found to confer resistance to drought or salt stress in plants, and overexpression of these transcription factors enhanced the drought or salt tolerance of transgenic plants [7, 48–52]. Previous research has found that the expression levels of wheat (Triticum aestivum L.) TaWRKY10 were upregulated 1.5-fold under PEG treatment and 2.6-fold under salt stress . In the C. korshinskii and Arabidopsis datasets, we detected eight WRKY transcription factors that showed increased expression in both drought-treated C. korshinskii and Arabidopsis (Additional file 20), especially Cko_contig_52997 (upregulated 30.8-fold in drought-treated C. korshinskii and 12.3-fold in drought-treated Arabidopsis). The characteristics of WRKY expression in C. korshinskii and Arabidopsis were consistent with findings from research on wheat under drought stress . The expression of the Arabidopsis zinc-finger proteins Zat 11 and Zat 12 is stimulated by drought (Zat 12 was upregulated 10-fold in leaves) and salt stresses (Zat 11 was upregulated 95-fold in roots) . In this study, we found that a zinc finger protein, Cko_contig_79632, was upregulated 255.8-fold in drought-treated C. korshinskii but downregulated slightly in drought-treated Arabidopsis, and Cko_contig_85079 showed increased transcript abundance (20.0-fold increase) in salt-treated C. korshinskii but was only upregulated 0.9-fold in salt-treated Arabidopsis. Interestingly, we found that eight zinc-finger protein contigs were downregulated in drought- or salt-treated C. korshinskii. These findings support the possibility that zinc-finger proteins could play key roles as both positive and negative regulators of plant defenses .
We also identified differences in the expression profiles of some other transcription factors in C. korshinskii and Arabidopsis under the same stresses (Additional file 20); for example, some NAC, MYB, and ERFs were downregulated in drought- or salt-treated C. korshinskii. This result indicated that some MYB and ERF transcription factors might function as negative regulators of drought or salt tolerance in C. korshinskii.
Functional genes related to the ROS-scavenging system and ion transport in C. korshinskii
ROS are produced in plant cells under normal and stressful conditions  and can be removed by the antioxidant system [55, 56]. In C. korshinskii, antioxidant enzymes such as superoxide dismutase (SOD), catalase (CAT), and peroxidase (POX) protect cells and subcellular systems from drought or salt damage caused by ROS (Additional file 20). For example, Cko_contig_105070 showed strong homology to CAT and was transcribed at high levels in salt-treated C. korshinskii (180.9-fold increase) and upregulated 1.13-fold in salt-treated Arabidopsis (Additional file 20). Therefore, it is most likely associated with enhanced CAT activity and responsible for converting the increased H2O2 into H2O and O2 . Ascorbate peroxidase (APX) represents the ascorbate-dependent H2O2-scavenging mechanism in plants, which reduces H2O2 to H2O using ascorbic acid as an electron acceptor . Eleven APX proteins were uniquely expressed in treated C. korshinskii (Additional file 5) but were not detected in Arabidopsis datasets. Two of these proteins (Cko_contig_33634, Cko_contig_46972) showed increased transcript abundance under drought stress. In the glutathione peroxidase (GPX) pathway, GPX can reduce H2O2 to the corresponding hydroxyl compounds [21, 58]. Interestingly, in our datasets, 12 GPX proteins were expressed uniquely in drought- or salt-treated C. korshinskii. In particular, Cko_contig_17138 was upregulated 2.4-fold in drought-treated C. korshinskii and upregulated 79.1 % under salt stress, no homologous contigs were found in published Arabidopsis datasets (Additional file 5).
It is possible that the survival of plants under salt stress depends on the accumulation of high levels of ions in the tissue for the maintenance of osmotic balance . Cko_contig_24748, associated with a K+-transporter, was uniquely upregulated 16.5-fold under salt stress in C. korshinskii but was not found in the considered Arabidopsis datasets.
The transcript profile was consistent with the profile of K+-transporters in R. trigyna, most of which were upregulated under salt stress . This result suggests that K+-transporters were enhanced and Na+/K+ homeostasis reestablished in salt-treated C. korshinskii, as seen in R. trigyna under salt stress . It is possible that plasma-membrane H+-ATPases (PM-H+-ATPases) are involved in the mitigation of physiological disturbances imposed by salt stress . Interestingly, in this study, Cko_contig_104930, encoding plasma membrane H+-ATPases (PM-H+-ATPases), showed significant increases in expression (3.1-fold increase) in salt-treated C. korshinskii but remained relatively constant in salt-treated Arabidopsis. This finding echoes the results of other studies on salt-treated plants, such as Oryza sativa Indica cultivars . Transcripts encoding proteins associated with functions such as diterpenoid biosynthesis mechanosensitive ion channels (Cko_contig_66582) accumulated in both drought- and salt-treated C. korshinskii but were not detected in Arabidopsis datasets. These results indicate that C. korshinskii may have the specific ability to transduce these stress stimuli into electrical signals to enhance drought or salt tolerance. In summary, a small number of contigs uniquely found in C. korshinskii but not in Arabidopsis suggest that C. korshinskii may have a stress-tolerance mechanism distinct from that of the model plant Arabidopsis.
Contigs in response to gibberellin biosynthesis in C. korshinskii
Plant growth and development are regulated by gibberellins (GAs), which are regulated by numerous genes involved in GA biosynthesis and degradation. GA biosynthesis is tightly controlled by Gibberellin 20-beta-dioxygenase (GA20ox) and Gibberellin 3-beta-dioxygenase (GA3ox). Gibberellin 2-beta-dioxygenase (GA2ox) catalyzes the deactivation of bioactive GAs . In Arabidopsis, the cellular concentration of bioactive GAs was reduced via an increase in expression of GA2ox7, resulting in an accumulation of DELLA, an inhibitor of plant growth . Based on transcriptome-sequence analysis of C. korshinskii and Arabidopsis, we identified nine drought- or salt-responsive GA2ox proteins, but these contigs were not found in Arabidopsis datasets (Additional file 5). Only five GA3ox proteins and Cko_contig_63577 were overexpressed in both-drought treated C. korshinskii and in researched Arabidopsis datasets (Additional file 20). Two GA20ox proteins were detected at higher levels in drought- or salt-treated C. korshinskii, but no homologous gene was found in the Arabidopsis datasets (Additional file 5). In addition, three contigs associated with DELLA protein were downregulated in both treated C. korshinskii and Arabidopsis (Additional file 20). Thus, we propose that GAs are upregulated under salt stress for normal growth of halophytes such as C. korshinskii, with accompanying repression of the expression of DELLA proteins. Our prediction is supported to some degree by studies of a GA-deficient mutant of Arabidopsis [62, 63].
Contigs associated with carbon metabolism and osmotic regulation in C. korshinskii and Arabidopsis
Responses of carbon metabolism or starch metabolism to drought stress have been found in many plants, such as alfalfa (Medicago sativa)  and maize (Zea mays) . In C. korshinskii, seven sugar-transport proteins, 13 carbohydrate-transport proteins, one granule-bound glycogen (starch) synthase protein and two sucrose-UDP glucosyltransferase proteins were found to be differentially expressed under drought or salt stress (Additional file 5). In particular, Cko_contig_43585, encoding a sugar-transport protein, was upregulated 568.5-fold in salt-treated samples. The expression data suggest that sugar or carbon metabolism, as well as various sugar-related signaling pathways, are influenced by drought stress in C. korshinskii.
The synthesis and transport of substances involved in osmotic regulation (such as proline and other amino acids) promote salt tolerance in most plants and crops [3, 66, 67]. Previous studies have reported that ProT accumulation is increased under salt stress in several plants, such as barley and mangrove [68, 69]. Our experimental results showed an increase in expression of the proline transporter gene (ProT, Cko_contig_90699), which was uniquely upregulated 2.6-fold and 1.7-fold in drought- and salt-treated C. korshinskii, respectively, but not in Arabidopsis. This drought- or salt-treated expression pattern in C. korshinskii was similar to the findings in barley and mangrove [68, 69]. In yeast, amino-acid permeases (AAPs) complement proline-uptake deficiency and mediate the efficient transport of proline, alanine, and valine, but the expression of most members of the AAP-gene family was repressed under high-salt conditions . Surprisingly, the transcription of AAP (Cko_contig_4505) was uniquely induced (3.3-fold increase) in drought-treated C. korshinskii, and no homologous gene was found in published Arabidopsis datasets. This finding suggests that C. korshinskii had the distinct ability to synthesize and transport various osmoprotectants, such as proline, osmotin, and amino acids, to maintain cell turgor under osmotic stress, which is not seen in the model plant Arabidopsis or in other eukaryotes.
Contigs associated with brassinosteroid-biosynthesis pathways in C. korshinskii
Brassinolides regulate plant growth and development and have been known to improve plant tolerance to abiotic stresses . The results of the RNA-seq expression analysis showed that 30 contigs involved in brassinosteroid biosynthesis pathway were regulated in treated C. korshinskii (Additional file 5). Cko_contig_7817 (85A1/2 or 85A2), a homolog of brassinosteroid-6-oxidase family protein, was differentially regulated in treated C. korshinskii and Arabidopsis transcriptome datasets (Additional file 20, Additional file 21). We further found four contigs, Cko_contig_82841 (188.8.131.52), Cko_contig_8184 (184.108.40.206), Cko_contig_32733 (220.127.116.11), and Cko_contig_74975 (18.104.22.168), that are involved in starch and sucrose metabolism in drought- or salt-treated C. korshinskii, but no homologous genes were detected in Arabidopsis datasets (Additional file 5). These genes are necessary to enhance the biosynthesis of soluble sugars such as glucose, sucrose, and fructose and to maintain cell osmotic adjustment and balance (Additional files 22 and 23).
This work presents an original transcriptome-sequencing analysis of mixed RNA from C. korshinskii whole seedlings (leaves, stems, and roots) under drought and salt stresses. A total of 129,451 cDNA contigs from whole C. korshinskii plants were obtained, of which 54.51 % were annotated with gene descriptions, while the remaining 45.49 % were unknown. A large number of candidate functional genes were identified as being potentially involved in multiple mechanisms that might contribute to drought or salt tolerance, meriting further investigation. Several key pathways, such as metabolic pathways, biosynthesis of secondary metabolites, ribosome, purine metabolism and plant hormone signal transduction, and GO terms, such as metabolic process, cell or cell part, and catalytic activity, were identified as being involved in drought and salt tolerance in C. korshinskii. In addition, hundreds of DEGs appeared to be expressed uniquely in C. korshinskii under drought and salt stress, respectively, compared to the published Arabidopsis database. These unknown genes and uniquely expressed DEGs were presumably genes unique to C. korshinskii that are involved in drought or salt tolerance. This information of C. korshinskii genes will be very useful in future studies on molecular adaptations to abiotic stress in leguminous plants and will facilitate the genetic manipulation of important plants.
Plant materials, growth conditions, and stress treatments
Seeds of C. korshinskii collected from Turpan Desert Botanical Garden, Chinese Academy of Sciences were sown in plastic pots (12 × 12 cm) filled with a nonsterile sandy soil irrigated with 1 L sterile water. After germination, C. korshinskii seedlings were grown in a greenhouse with a 16 h light/8 h dark photoperiod at a temperature of 22 °C ± 1 °C and a relative humidity of 50 ± 10 % and reirrigated with 0.7 L sterile water. Two weeks later, the salt-treated plants were watered with 400 mM NaCl, the drought-treated plants were unwatered, and control plants were irrigated with 1 L sterile water. One month later, when seedlings were approximately 10-cm high, all seedlings, leaves, stems, and roots were harvested for analysis.
cDNA-library preparation and Solexa sequencing for transcriptome analysis
Total RNA was extracted from C. korshinskii whole seedlings treated with 400 mM NaCl or drought stress using a Plant RNA Extraction Kit (Autolab, China), following the manufacturer’s protocol. The concentration and quality of each RNA sample were determined using a NanoDrop 2000TM micro-volume spectrophotometer (Thermo Scientific, Waltham, MA, USA) and gel electrophoresis. The poly (A) mRNA was isolated from 10 mg total RNA using magnetic oligo(dT) beads, then divided into short fragments by fragmentation buffer (Ambion, Austin, TX, USA). The first-strand cDNA was synthesized by random-hexamer primers, and then the second-strand cDNA was synthesized using DNA polymerase I (New England Biolabs), RNase H (Invitrogen), buffer, and dNTPs. Short DNA fragments were purified with a QIAquick PCR Purification Kit (Qiagen Inc., Valencia, CA, USA), then subjected to end repair and poly(A) addition and ligated to sequencing adaptors. Suitable fragments (350–450 bp) were purified by agarose-gel electrophoresis and gathered by PCR amplification. One cDNA library (Control sample) was sequenced on a PE flow cell and three cDNA libraries (Control, drought-treated and salt-treated samples) were sequenced on a SE flow cell using the Illumina HiSeq™ 2000 platform.
De novo transcriptome assembly
Filter reads with adapter contamination.
Filter reads in which unknown nucleotides >5 %.
Filter reads in which more than 20 % of bases showed a Q-value less than 20.
High-quality clean reads obtained after the above filtering and correction steps were used for de novo assembly using the Trinity method  with an optimized k-mer length of 25 by default and with all other parameters set to default values. The assembled sequences were further incorporated into distinct contigs by using the CD-HIT  and CAP3  programs. The parameters in CAP3 were “-p 98 -o 50”. CD-HIT was also used to cluster transcripts on the basis of sequence similarity (-c sequence identity threshold >0.98) and retain only the longest transcript from each cluster; redundant contigs were excluded.
Multiple complementary approaches were adopted for alignment and analysis. All contig sequences were entered as search terms in the NCBI non-redundant nucleotide sequence (nr) database, the Kyoto Encyclopedia of Genes and Genomes (KEGG) database, the Cluster of Orthologous Groups of proteins (COG) database, and the Swiss-Prot protein database using a cut-off E-value of 10−5. Gene Ontology (GO) annotations of the contigs were determined using Blast2go (https://www.blast2go.com/)  according to the molecular function, biological process, and cellular component ontologies (http://www.geneontology.org/) with an E-value threshold of 10−5; the GO level of the GO distribution figure was three. The best alignment results were used to determine the 5′-to-3′ orientation of the contigs. Annotated contigs with no GO hits were investigated in terms of protein domains using InterProScan databases (http://www.ebi.ac.uk/interpro/). The “GetORF” program from the EMBOSS software package was used to predict the open reading frames (ORFs) of each contig . Enzyme mapping of the annotated contigs was performed by direct GO-to-enzyme mapping and used to query the KEGG to define the KEGG orthologs (KOs). The KEGG mapping tool was used to plot these KOs into a complete metabolic atlas .
Identification of differential expression transcripts (DEGs) and functional analysis
The expression of each contig was calculated by the RPKM method . The method of Audic and Claverie  was applied to determine the P-value, which corresponds to the differential transcript expression of the contigs. The false discovery rate (FDR) is used to determine the P-value threshold in multiple tests. FDR ≤0.05 and the absolute value of log2Ratio ≥1 were used as the threshold to judge the significance of differential gene expression. The differentially expressed genes were used for GO- and KEGG-enrichment analyses. GO-enrichment analysis was performed using the Blast2go software, and the P-values were calculated using the Benjamini-Hochberg (BH) correction . We selected a corrected P-value <0.05 as a threshold to determine which DEGs were significantly enriched in GO terms. In KEGG-enrichment analysis, the cellular metabolism, biochemical pathways, and potential biological behaviors of the differentially expressed genes were examined (www.genome.jp/kegg/). A corrected P-value <0.05 was used as the threshold to determine significant enrichment of the gene sets.
We performed a new set of BLAST searches to compare C. korshinskii contigs to genes expressed in Arabidopsis under salt and/or drought stresses, in order to identify genes that might be expressed uniquely in C. korshinskii under salt and/or drought stresses. The RNA-seq data of Arabidopsis can be found in ArrayExpress and the accession number was E-GEOD-48235.
Real-time quantitative reverse-transcription polymerase chain reaction (qRT-PCR) for validation and expression-pattern analysis
Total RNA was extracted from whole seedlings and organs (leaves, stems, and roots) of C. korshinskii treated with 400 mM NaCl or drought stress using a Plant RNA Extraction Kit (Autolab, China) according to the manufacturer’s protocol. Reverse-transcription reactions were performed using Superscript II reverse transcriptase (Invitrogen, Grand Island, NY, USA) according to the manufacturer’s instructions and 1 μg of total RNA and oligo(dT) primers.
Primers for qRT-PCR were designed using the Primer Premier 5.0 software (Premier Biosoft Int., Palo Alto, CA, USA). The sequences of the specific primer sets are listed in Additional file 19. The qRT-PCR analysis was performed using an ABI Prism 7500 sequence detector (Applied Biosystems, Foster City, CA, USA) and the SYBR® Premix Ex TaqTM Kit (TaKaRa, Tokyo, Japan) to detect transcript abundance. The qRT-PCR reactions were performed in 20-μl volumes containing 2 μl first-strand cDNA, 200 nM of each primer, and 10 μl of the 1 × SYBR PCR mixture with the following cycling programs: denaturation at 95 °C for 30 s, then 40 cycles of denaturation at 95 °C for 3 s and annealing and extension at 60 °C for 30 s. The specificity of the amplifications of all products were verified by melting curve analysis, which was conducted at 95 °C for 15 s, 60 °C for 60 s, and 95 °C for 15 s. Three replicates were conducted in parallel, and the results were normalized to the expression level of the constitutive 18S rRNA gene. The relative quantitative method (ΔΔCt) was used to evaluate quantitative variation.
Availability of data and materials
The RNA-seq data of Caragana korshinskii in this study was deposited into the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA) database under number SRP061143.
The RNA-seq data of Arabidopsis that was used as reference data for the comparison between Arabidopsis and C. korshinskii can be found in ArrayExpress and the accession number was E-GEOD-48235.
We thank Dr. Chaoying He (Institute of Botany, Chinese Academy of Sciences) for the helpful advice and technical assistance. We also thank Dr. Liying Song (Institute of Genetics and Developmental Biology, Chinese Academy of Sciences) and Dr. Jin Xu and Jue Ruan (Beijing Institute of Genomics, Chinese Academy of Sciences) for their sequencing work. This work was supported by a project (2011AA10020105) at the Ministry of Science and Technology of China and a project from National Natural Science Foundation of China (31400570).
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Zhu JK. Cell signaling under salt, water and cold stresses. Curr Opin Plant Biol. 2001;4:401–6.View ArticlePubMedGoogle Scholar
- Zhu JK. Salt and drought stress signal transduction in plants. Annu Rev Plant Biol. 2002;53:247–73.View ArticlePubMedPubMed CentralGoogle Scholar
- Cockburn W, Whitelam GC, Broad A, Smith J. The participation of phytochrome in the signal transduction pathway of salt stress responses in Mesembryanthemum crystallinum L. J Exp Bot. 1996;47:647–53.View ArticleGoogle Scholar
- Liska AJ, Shevchenko A, Pick U, Katz A. Enhanced photosynthesis and redox energy production contribute to salinity tolerance in Dunaliella as revealed by homology-based proteomics. Plant Physiol. 2004;136:2806–17.View ArticlePubMedPubMed CentralGoogle Scholar
- Kawaguchi R, Girke T, Bray EA, Bailey-Serres J. Differential mRNA translation contributes to gene regulation under non-stress and dehydration stress conditions in Arabidopsis thaliana. Plant J. 2004;38(5):823–39.View ArticlePubMedGoogle Scholar
- Rabbani MA, Maruyama K, Abe H, Khan MA, Katsura K, Ito Y, et al. Monitoring expression profiles of rice genes under cold, drought, and high-salinity stresses and abscisic acid application using cDNA microarray and RNA gel-blot analyses. Plant Physiol. 2003;133(4):1755–67.View ArticlePubMedPubMed CentralGoogle Scholar
- Seki M, Narusaka M, Ishida J, Nanjo T, Fujita M, Oono Y, et al. Monitoring the expression profiles of 7000 Arabidopsis genes under drought, cold and high-salinity stresses using a full-length cDNA microarray. Plant J. 2002;31(3):279–92.View ArticlePubMedGoogle Scholar
- Shinozaki K, Yamaguchi-Shinozaki K. Gene networks involved in drought stress response and tolerance. J Exp Bot. 2007;58(2):221–7.View ArticlePubMedGoogle Scholar
- Agarwal P, Agarwal P, Reddy M, Sopory S. Role of DREB transcription factors in abiotic and biotic stress tolerance in plants. Plant Cell Rep. 2006;25(12):1263–74.View ArticlePubMedGoogle Scholar
- Shinozaki K, Yamaguchi-Shinozaki K. Gene expression and signal transduction in water-stress response. Plant Physiol. 1997;115(2):327–34.View ArticlePubMedPubMed CentralGoogle Scholar
- Jaglo-Ottosen KR, Gilmour SJ, Zarka DG, Schabenberger O, Thomashow MF. Arabidopsis CBF1 overexpression induces COR genes and enhances freezing tolerance. Science. 1998;280(5360):104–6.View ArticlePubMedGoogle Scholar
- Shinozaki K, Yamaguchi-Shinozaki K, Seki M. Regulatory network of gene expression in the drought and cold stress responses. Curr Opin Plant Biol. 2003;6(5):410–7.View ArticlePubMedGoogle Scholar
- Ciftci-Yilmaz S, Morsy MR, Song L, Coutu A, Krizek BA, Lewis MW, et al. The EAR-motif of the Cys2/His2-type zinc finger protein Zat7 plays a key role in the defense response of Arabidopsis to salinity stress. J Biol Chem. 2007;282:9260–8.View ArticlePubMedGoogle Scholar
- Jakoby M, Weisshaar B, Dröge-Laser W, Vicente-Carbajosa J, Tiedemann J, Kroj T, et al. bZIP transcription factors in Arabidopsis. Trends Plant Sci. 2002;7(3):106–11.View ArticlePubMedGoogle Scholar
- Marè C, Mazzucotelli E, Crosatti C, Francia E, Stanca A, Cattivelli L. Hv-WRKY38: a new transcription factor involved in cold- and drought-response in barley. Plant Mol Biol. 2004;55(3):399–416.View ArticlePubMedGoogle Scholar
- Olsen AN, Ernst HA, Leggio LL, Skriver K. NAC transcription factors: structurally distinct, functionally diverse. Trends Plant Sci. 2005;10(2):79–87.View ArticlePubMedGoogle Scholar
- Harb A, Krishnan A, Ambavaram MMR, Pereira A. Molecular and physiological analysis of drought stress in Arabidopsis reveals early responses leading to acclimation in plant growth. Plant Physiol. 2010;154(3):1254–71.View ArticlePubMedPubMed CentralGoogle Scholar
- Ma J, Zhang M, Xiao X, You J, Wang J, Wang T, et al. Global transcriptome profiling of Salicornia europaea L. shoots under NaCl treatment. PLoS One. 2013;8(6):e65877. doi:10.1371/journal.pone.0065877.View ArticlePubMedPubMed CentralGoogle Scholar
- Skirycz A, Claeys H, De Bodt S, Oikawa A, Shinoda S, Andriankaja M, et al. Pause-and-stop: the effects of osmotic stress on cell proliferation during early leaf development in Arabidopsis and a role for ethylene signaling in cell cycle arrest. Plant Cell. 2011;23(5):1876–88.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhang X, Yao D, Wang Q, Xu W, Wei Q, Wang C, et al. mRNA-seq analysis of the Gossypium arboreum transcriptome reveals tissue selective signaling in response to water stress during seedling stage. PLoS One. 2013;8(1):e54762. doi:10.1371/journal.pone.0054762.View ArticlePubMedPubMed CentralGoogle Scholar
- Dang Z, Zheng L, Wang J, Gao Z, Wu S, Qi Z, et al. Transcriptomic profiling of the salt-stress response in the wild recretohalophyte Reaumuria trigyna. BMC Genomics. 2013;14:29.View ArticlePubMedPubMed CentralGoogle Scholar
- Ranjan A, Nigam D, Asif MH, Singh R, Ranjan S, Mantri S, et al. Genome wide expression profiling of two accession of G. herbaceum L. in response to drought. BMC Genomics. 2012;13:94.View ArticlePubMedPubMed CentralGoogle Scholar
- Thumma BR, Sharma N, Southerton SG. Transcriptome sequencing of Eucalyptus camaldulensis seedlings subjected to water stress reveals functional single nucleotide polymorphisms and genes under selection. BMC Genomics. 2012;13:364.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhang ZS, Li XR, Liu LC, Jia RL, Zhang JG, Wang T. Distribution, biomass, and dynamics of roots in a revegetated stand of Caragana korshinskii in the tengger desert, northwestern China. J Plant Res. 2009;122:109–119.View ArticlePubMedGoogle Scholar
- Zhang X. Principles and optimal models for development of Maowusu sandy grassland. Acta Phytoecologica Sinica. 1994;18:1–16.Google Scholar
- Fu HC. Caragana. In: Fu HC, editor. Flora Intramongolica. 2nd ed. Huhhot: Inner Mongolia People Press; 1989. p. 236–8.Google Scholar
- Fang XW, Turner NC, Li FM, Li WJ, Guo XS. Caragana korshinskii seedlings maintain positive photosynthesis during short-term, severe drought stress. Photosynthetica. 2011;49(4):603–9.View ArticleGoogle Scholar
- Fang Y, You J, Xie K, Xie W, Xiong L. Systematic sequence analysis and identification of tissue-specific or stress-responsive genes of NAC transcription factor family in rice. Mol Gen Genomics. 2008;280(6):547–63.View ArticleGoogle Scholar
- Niu CF, Wei W, Zhou QY, Tian AG, Hao YJ, Zhang WK, et al. Wheat WRKY genes TaWRKY2 and TaWRKY19 regulate abiotic stress tolerance in transgenic Arabidopsis plants. Plant Cell Environ. 2012;35(6):1156–70.View ArticlePubMedGoogle Scholar
- Dugas DV, Monaco MK, Olson A, Klein RR, Kumari S, Ware D, et al. Functional annotation of the transcriptome of Sorghum bicolor in response to osmotic stress and abscisic acid. BMC Genomics. 2011;12:51.View ArticleGoogle Scholar
- Li SB, Wang B, Bai YQ, Wang L. Studies on the characteristics of the main forage shrub species in Yanchi sandy land. Sci Silvae Sin. 2000;36(4):119–25.Google Scholar
- Wang XM, Gao ZY. Study on process technology and its key technique of wood based panels from desert shrubs. Chinese Wood Industry. 2003;17:11–3.Google Scholar
- Wang Z, Gao HW, Wu YQ, Han JG. Genetic diversity and population structure of Caragana korshinskii revealed by AFLP. Crop Sci. 2007;47(4):1737–43.View ArticleGoogle Scholar
- Yang DH, Song LY, Hu J, Yin WB, Li ZG, Chen YH, et al. Enhanced tolerance to NaCl and LiCl stresses by over-expressing Caragana korshinskii sodium/proton exchanger 1 (CkNHX1) and the hydrophilic C terminus is required for the activity of CkNHX1 in Atsos3-1 mutant and yeast. Biochem Biophys Res Commun. 2012;417(2):732–7.View ArticlePubMedGoogle Scholar
- Wang X, Chen X, Liu Y, Gao H, Wang Z, Sun G. CkDREB gene in Caragana korshinskii is involved in the regulation of stress response to multiple abiotic stresses as an AP2/EREBP transcription factor. Mol Biol Rep. 2011;38(4):2801–11.View ArticlePubMedGoogle Scholar
- Wang X, Dong J, Liu Y, Gao H. A novel dehydration-rsponsive element-binding protein from Caragana korshinskii is involved in the response to multiple abiotic stresses and enhances stress tolerance in transgenic tobacco. Plant Mol Biol Rep. 2010;28(4):664–75.View ArticleGoogle Scholar
- Ding Y, Liu N, Virlouvet L, Riethoven JJ, Fromm M, Avramova Z. Four distinct types of dehydration stress memory genes in Arabidopsis thaliana. BMC Plant Biol. 2013;13:229.View ArticlePubMedPubMed CentralGoogle Scholar
- Long Y, Wang Y, Wu S, Wang J, Tian X, Pei X. De novo assembly of transcriptome sequencing in Caragana korshinskii Kom. and characterization of EST-SSR markers. PLoS One. 2015;10(1):e0115805. doi:10.1371/journal.pone.0115805.View ArticlePubMedPubMed CentralGoogle Scholar
- Huang P, Ju HW, Min JH, Zhang X, Kim SH, Yang KY, et al. Overexpression of L-type lectin-like protein kinase 1 confers pathogen resistance and regulates salinity response in Arabidopsis thaliana. Plant Sci. 2013;203–204:98–106.View ArticlePubMedGoogle Scholar
- Xu DL, Long H, Liang JJ, Zhang J, Chen X, et al. De novo assembly and characterization of the root transcriptome of Aegilops variabilis during an interaction with the cereal cyst nematode. BMC Genomics. 2012;13:133.View ArticlePubMedPubMed CentralGoogle Scholar
- Liu MY, Qiao GR, Jiang J, Yang H, Xie LH, et al. Transcriptome sequencing and de novo analysis for Ma Bamboo (Dendrocalamus latiflorus Munro) using the Illumina platform. PLoS One. 2012;7:e46766.View ArticlePubMedPubMed CentralGoogle Scholar
- Novaes E, Drost DR, Farmerie WG, Pappas GJ, Grattapaglia D, et al. High-throughput gene and SNP discovery in Eucalyptus grandis, an uncharacterized genome. BMC Genomics. 2008;9:312.View ArticlePubMedPubMed CentralGoogle Scholar
- Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008;5:621–8.View ArticlePubMedGoogle Scholar
- Umezawa T, Yoshida R, Maruyama K, Yamaguchi-Shinozaki K, Shinozaki K. SRK2C, a SNF1-related protein kinase 2, improves drought tolerance by controlling stress-responsive gene expression in Arabidopsis thaliana. PNAS. 2004;101:17306–11.View ArticlePubMedPubMed CentralGoogle Scholar
- Ge Y, Li Y, Zhu YM, Bai X, Lv DK, Guo D, et al. Global transcriptome profiling of wild soybean (Glycine soja) roots under NaHCO3 treatment. BMC Plant Biol. 2010;10:153.View ArticlePubMedPubMed CentralGoogle Scholar
- Sun XL, Yua QY, Tang LL, Ji W, Bai X, Cai H, et al. GsSRK, a G-type lectin S-receptor-like serine/threonine protein kinase, is a positive regulator of plant tolerance to salt stress. J Plant Physiol. 2013;170(5):505–15.View ArticlePubMedGoogle Scholar
- Kempa S, Rozhon W, Samaj J, Erban A, Baluska F, Becker T, et al. A plastid-localized glycogen synthase kinase 3 modulates stress tolerance and carbohydrate metabolism. Plant J. 2007;49:1076–90.View ArticlePubMedPubMed CentralGoogle Scholar
- Takasaki H, Maruyama K, Kidokoro S, Ito Y, Fujita Y, Shinozaki K, et al. The abiotic stress-responsive NAC-type transcription factor OsNAC5 regulates stress-inducible genes and stress tolerance in rice. Mol Gen Genomics. 2010;284(3):173–83.View ArticleGoogle Scholar
- Jeong JS, Kim YS, Baek KH, Jung H, Ha SH, Choi YD, et al. Root-specific expression of OsNAC10 improves drought tolerance and grain yield in rice under field drought conditions. Plant Physiol. 2010;153(1):185–97.View ArticlePubMedPubMed CentralGoogle Scholar
- Hao YJ, Wei W, Song QX, Chen HW, Zhang YQ, Wang F, et al. Soybean NAC transcription factors promote abiotic stress tolerance and lateral root formation in transgenic plants. Plant J. 2011;68(2):302–13.View ArticlePubMedGoogle Scholar
- Zhou QY, Tian AG, Zou HF, Xie ZM, Lei G, Huang J, et al. Soybean WRKY-type transcription factor genes, GmWRKY13, GmWRKY21, and GmWRKY54, confer differential tolerance to abiotic stresses in transgenic Arabidopsis plants. Plant Biotechnol J. 2008;6(5):486–503.View ArticlePubMedGoogle Scholar
- Hu H, Dai M, Yao J, Xiao B, Li X, Zhang Q, et al. Overexpressing a NAM, ATAF, and CUC (NAC) transcription factor enhances drought resistance and salt tolerance in rice. Proc Natl Acad Sci U S A. 2006;103:12987–92.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang C, Deng P, Chen L, Wang X, Ma H, Hu W, et al. A wheat WRKY transcription factor TaWRKY10 confers tolerance to multiple abiotic stresses in transgenic tobacco. PLoS One. 2013;8(6):e65120.View ArticlePubMedPubMed CentralGoogle Scholar
- Mittler R, Kim Y, Song L, Coutu J, Coutu A, Ciftci-Yilmaz S, et al. Gain- and loss-of-function mutations in Zat10 enhance the tolerance of plants to abiotic stress. FEBS Lett. 2006;580:6537–42.View ArticlePubMedPubMed CentralGoogle Scholar
- Grene R. Oxidative stress and acclimation mechanisms in plants. In: Somerville CR, Meyerowitz EM, editors. The Arabidopsis Book. Rockville: American Society of Plant Biologists; 2002. doi:10.1199/tab.0036.1. http//aspb.org/publications/arabidopsis/.Google Scholar
- Alscher RG, Erturk N, Heath LS. Role of superoxide dismutases (SODs) in controlling oxidative stress in plants. J Exp Bot. 2002;53:1331–41.View ArticlePubMedGoogle Scholar
- Daud MK, Mei L, Najeeb U, Khan M, Deeba F, Raza I, et al. In vitro cadmium-induced alterations in growth and oxidative metabolism of upland Cotton (Gossypium hirsutum L.). Sci World J. 2014;2014:309409.Google Scholar
- Chen FG, Zhang S, Jiang H, Ma WJ, Korpelainen H, Li CY. Comparative proteomics analysis of salt response reveals sex-related photosynthetic inhibition by salinity in Populus cathayana cuttings. J Proteome Res. 2011;10:3944–58.View ArticlePubMedGoogle Scholar
- Wu YX, Ding N, Zhao X, Zhao MG, Chang ZQ, Liu JQ, et al. Molecular characterization of PeSOS1: the putative Na+/H+ antiporter of Populus euphratica. Plant Mol Biol. 2007;65:1–11.View ArticlePubMedGoogle Scholar
- Sahu BB, Shaw BP. Salt-inducible isoform of plasma membrane H+ATPase gene in rice remains constitutively expressed in natural halophyte, Suaeda maritima. J Plant Physiol. 2009;166(10):1077–89.View ArticlePubMedGoogle Scholar
- Richards DE, King KE, Ait-ali T, Harberd NP. How gibberellin regulates plant growth and development: a molecular genetic analysis of gibberellinsignaling. Annu Rev Plant Physiol Plant Mol Biol. 2001;52:67–88.View ArticlePubMedGoogle Scholar
- Magome H, Yamaguchi S, Hanada A, Kamiya Y, Oda K. The DDF1 transcriptional activator upregulates expression of a gibberellin-deactivating gene, GA2ox7, under high-salinity stress in Arabidopsis. Plant J. 2008;56:613–26.View ArticlePubMedGoogle Scholar
- Magome H, Yamaguchi S, Hanada A, Kamiya Y, Oda K. Dwarf and delayed-flowering 1, a novel Arabidopsis mutant deficient in gibberellin biosynthesis because of overexpression of a putative AP2 transcription factor. Plant J. 2004;37:720–9.View ArticlePubMedGoogle Scholar
- Naya L, Ladrera R, Ramos J, González EM, Arrese-Igor C, Minchin FR, et al. The response of carbon metabolism and antioxidant defenses of alfalfa nodules to drought stress and to the subsequent recovery of plants. Plant Physiol. 2007;144(2):1104–14.View ArticlePubMedPubMed CentralGoogle Scholar
- Kakumanu A, Ambavaram MMR, Klumas C, Krishnan A, Batlang U, Myers E, et al. Effects of drought on gene expression in maize reproductive and leaf meristem tissue revealed by RNA-Seq. Plant Physiol. 2012;160:846–67.View ArticlePubMedPubMed CentralGoogle Scholar
- Yancey PH, Clark ME, Hand SC, Bowlus RD, Somero GN. Living with water-stress evolution of osmolyte systems. Science. 1982;217:1214–22.View ArticlePubMedGoogle Scholar
- Chang B, Yang L, Cong W, Zu Y, Tang Z. The improved resistance to high salinity induced by trehalose is associated with ionic regulation and osmotic adjustment in Catharanthus roseus. Plant Physiol Biochem. 2014;77(2):140–148.View ArticlePubMedGoogle Scholar
- Ueda A, Shi WM, Sanmiya K, Shono M, Takabe T. Functional analysis of salt-inducible proline transporter of barley roots. Plant Cell Physiol. 2001;42:1282–9.View ArticlePubMedGoogle Scholar
- Waditee R, Hibino T, Tanaka Y, Nakamura T, Incharoensakdi A, et al. Functional characterization of betaine/proline transporters in betaine-accumulating mangrove. J Biol Chem. 2002;277:18373–82.View ArticlePubMedGoogle Scholar
- Rentsch D, Hirner B, Schmelzer E, Frommer WB. Salt stress-induced proline transportem and salt stress-repressed broad specificity amino acid permeases identified by suppression of a yeast amino acid permease-targeting mutant. Plant Cell. 1996;8:1437–46.View ArticlePubMedPubMed CentralGoogle Scholar
- Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29:644–52.View ArticlePubMedPubMed CentralGoogle Scholar
- Li W, Godzik A. Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics. 2006;22:1658–9.View ArticlePubMedGoogle Scholar
- Huang X, Madan A. CAP3: a DNA sequence assembly program. GenomeRes. 1999;9:868–77.Google Scholar
- Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21:3674–6.View ArticlePubMedGoogle Scholar
- Rice P, Longden I, Bleasby A. EMBOSS: the European molecular biology open software suite. Trends Genet. 2000;16:276–7.View ArticlePubMedGoogle Scholar
- Okuda S, Yamada T, Hamajima M, Itoh M, Katayama T, Bork P, et al. KEGG atlas mapping for global analysis of metabolic pathways. Nucleic Acids Res. 2008;36:W423–426.View ArticlePubMedPubMed CentralGoogle Scholar
- Audic S, Claverie JM. The significance of digital gene expression profiles. Genome Res. 1997;7:986–95.PubMedGoogle Scholar
- Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B. 1995;57:289–300.Google Scholar