Skip to main content

Transcriptomic profiling of the salt-stress response in the wild recretohalophyte Reaumuria trigyna



Reaumuria trigyna is an endangered small shrub endemic to desert regions in Inner Mongolia. This dicotyledonous recretohalophyte has unique morphological characteristics that allow it to tolerate the stress imposed by semi-desert saline soil. However, it is impossible to explore the mechanisms underlying this tolerance without detailed genomic information. Fortunately, newly developed high-throughput sequencing technologies are powerful tools for de novo sequencing to gain such information for this species.


Two sequencing libraries prepared from control (C21) and NaCl-treated samples (T43) were sequenced using short reads sequencing technology (Illumina) to investigate changes in the R. trigyna transcriptome in response to salt stress. Among 65340 unigenes, 35495 (52.27%) were annotated with gene descriptions, conserved domains, gene ontology terms, and metabolic pathways with a cut-off E-value of 10-5. These included 44 Gene Ontology (GO) terms, 119 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, and 25 Clusters of Orthologous Groups families. By comparing the transcriptomes from control and NaCl-treated plants, 5032 genes showed significantly differences in transcript abundance under salt stress (false discovery rate ≤ 0.001 and |log2Ratio| ≥ 1). These genes were significantly enriched in 29 KEGG pathways and 26 GO terms. The transcription profiles indicated that genes related to ion transport and the reactive oxygen species scavenging system were relevant to the morphological and physiological characteristics of this species. The expression patterns of 30 randomly selected genes resulted from quantitative real-time PCR were basically consistent with their transcript abundance changes identified by RNA-seq.


The present study identified potential genes involved in salt tolerance of R. trigyna. The globally sequenced genes covered a considerable proportion of the R. trigyna transcriptome. These data represent a genetic resource for the discovery of genes related to salt tolerance in this species, and may be a useful source of reference sequences for closely related taxa. These results can also further our understanding of salt tolerance in other halophytes surviving under sodic stress.


Soil salinity is one of the most significant abiotic stresses limiting plant growth, productivity, and geographical distribution [1, 2]. Approximately 20% of the world’s cultivated lands, and nearly 50% of irrigated lands, are affected by salinity [3, 4]. High salt levels impose two stresses on plants: osmotic stress, resulting from the lowered water availability because of the high osmotic pressure in the soil; and ionic stress arising from solute imbalance, in which there are increased levels of Na+ and Cl- in the cytosol and changes to the intracellular K+/Na+ ratio [5, 6]. Halophytes are the natural inhabitants of highly saline soils, and have evolved resistance strategies including efficient control of the uptake and compartmentalization of salt ions, synthesis of organic ‘compatible’ solutes, and unique morphological structures such as succulent leaves, salt glands, and bladders [7]. Such remarkable adaptations to salty soil make halophytes excellent candidates to investigate salt-tolerance mechanisms and to identify effective salt-response genes in this type of flora [8].

Studies on the salt-tolerance mechanisms of halophytes have been carried out recently in several model plants, such as Thellungiella halophila, Mesembryanthemum crystallinum, Suaeda, and Populus. Excellent progress has been made in understanding the physiological and molecular mechanisms underlying stress tolerance. However, knowledge about halophytes and their physiology is rather limited with respect to their wide taxonomic occurrence. It remains unclear whether halophytes have unique tolerance mechanisms, and whether such mechanisms are taxonomically linked, and/or have evolved to respond to interactions between salinity and other environmental variables. Therefore, studies on several species from different salt-stressed environments should be performed to answer these questions [9].

Reaumuria trigyna (Reaumuria Linn genus, family Tamaricaceae) is an endangered dicotyledonous shrub with the features of a recretohalophyte. It is regarded as a living fossil owing to its Tethys Ocean origin [1013]. This species is endemic to the Eastern Alxa–Western Ordos area (106°27’E–111°28’E, 39°13’N–40°52’N, elevation 1500–2100 m), a salinized desert in Inner Mongolia, China. The area is characterized by its high soil salinity (up to 0.7% salts), hyper-drought conditions (annual average precipitation of 140.9–302.2 mm), and low temperature (annual average temperature of 6.0–9.2°C) [14]. To adapt to the semi-desert and salty soil environment, R. trigyna has developed a distinct morphology that is characterized by succulent leaves and sunken stomata. These features are typical of a recretohalophyte—a halophyte with salt excretion glands. According to our previous studies, the concentrations of Na+, Cl-, and SO42- in the soil where R. trigyna grows are higher than those of K+, Mg2+, and Ca2+. The Na+ concentration in the soil is as high as 0.6 mg/g, with a Na+/K+ ratio of approximately 16. R. trigyna excretes Na+ and Cl- via its multicellular salt glands, but little K+ is excreted, with a Na+/K+ ratio of only 0.44. When seedlings were treated with 400 mM NaCl, R. trigyna showed a rapid increase in glutathione (GSH) from nearly 4.77 to 9.62 nmol·g-1 FW, and showed significant increases in the activities of superoxide dismutase (SOD) and peroxidase (POD) (P < 0.05) [15]. As a typical recretohalophyte growing in this area, R. trigyna could be considered as a reference species to represent those that survive in the severely adverse environment of the Ordos desert. Knowledge about the unique adaption strategies employed by R. trigyna will provide valuable references to the salt-tolerance mechanisms in these species and possibly in other halophytes growing in different environments. However, these objectives are difficult to achieve without detailed genetic and sequence information for this species. Fortunately, powerful high throughput sequencing technologies are now available to sequence the genome of plants de novo. Therefore, it is now possible to obtain the required sequences of R. trigyna to address important questions about its physiology, such as its mechanisms of salt-tolerance.

Next-generation sequencing (NGS) technology has developed rapidly in recent years. It not only provides rapid, cost-effective, and comprehensive analyses of complex nucleic acid populations for model plants or species closely related to model plants, but also provides opportunities to analyze non-model plants whose genomes have never been sequenced [1619]. This technology has been used widely in comparative transcriptomics to identify differences in transcript abundance among different cultivars, organs, and different treatment conditions [2022]. In the present study, we generated a transcriptome dataset to provide genetic information to explore the salt-tolerance mechanisms of R. trigyna using the Illumina HiSeq™ 2000 platform. We compared transcriptomes of salt-stressed and control plants to identify genes showing transcriptional changes, and identified the functions of the transcripts and the KEGG pathways that showed changes. We speculate that the assembled, annotated transcriptome sequences and transcript abundance patterns will provide a valuable genetic resource for further investigations of the molecular mechanisms of salt tolerance in this species, and possibly in other recretohalophytes.


Sequencing output and assembly

Two sequencing libraries were prepared from control (C21) and NaCl-treated samples (T43) to investigate the transcriptomic responses to salt-stress in R. trigyna. In total, 26.51 million raw reads were generated from C21 and 28.17 million raw reads were generated from T43 (Sequence Read Archive accession number: SRP008320; Table 1). The raw reads had an average length of 90 bp. Among all the raw reads, more than 92% had Phred-like quality scores at the Q20 level (an error probability of 1%). These were used for de novo assembly. As a result, 459601 (C21) and 476229 (T43) contigs, and 97752 (C21) and 101931 (T43) scaffolds, were generated. After removal of redundancy, 68076 (C21) and 71194 (T43) unigenes were clustered using the Gene Indices Clustering Tools (GICT) program (Table 1, deposited in the Transcriptome Shotgun Assembly Sequence database. Accession numbers are summarized in Additional files 1 and 2). The random distribution of reads in the above-assembled unigenes indicated that the 3’ ends of all assembled unigenes contained relatively fewer reads compared with other positions within these unigenes, which showed greater and more even distributions (Figure 1).

Table 1 Summary of sequencing and assembly results
Figure 1

Random distribution of sequencing reads in assembled unigenes. X-axis represents relative position of sequencing reads in the assembled sequences. Orientation of assembled unigenes is from the 5’ end to the 3’end. Y-axis indicates number of reads. Left: random distribution of C21 reads mapped to C21-unigenes. Right: random distribution of T43 reads mapped to T43-unigenes.

Using the CAP3 assembler, we obtained 65340 all-unigenes with an average length of 562 bp and an N50 of 770 bp by combining C21 (25.6 million) and T43 (27.3 million) clean reads (Table 1). Of the 65340 unigenes, 34745 were clustered from the assembled C21 (43506) and T43 (44558) unigenes. The remaining 30595 sequences were singletons with 14764 from C21 and 15831 from T43 unigenes. The length distributions of the unigenes from both C21 and T43 were similar with sequences in the 200–300 bp range making up 21.55% and 23.2% of C21 and T43 unigenes, respectively. The length of all-unigenes was dramatically increased following further assembly of the C21 and T43 unigenes, resulting in assembled sequences with lengths of > 300 bp (Figure 2). Ultimately, 54331 gap-free sequences were obtained by assembling all of the C21- and T43-unigenes.

Figure 2

Length distribution of C21-, T43- and all-unigenes, and gap distribution of all unigenes. A, B, C: Length distribution of C21-unigenes, T43-unigenes, and all-unigenes. Values in colored square brackets indicate range of unigene length. Unigene numbers and percentages (total number of unigenes in a certain length range) are shown next to pie chart. Total number of C21-, T43-, and all-unigenes was 68076, 71194, and 65340, respectively. D: Gap distribution of all unigenes. Numbers in square brackets represent ratio of N number for one unigene.

Functional annotation of assembled unigenes

All assembled high-quality unigenes were first blasted against the NCBI non-redundant (nr) database using BLASTX with a cut-off E-value of 10-5. Of the 65340 all-unigenes, 35271 (53.9%) returned at least one match at the E-value > 10-5. Because of the lack of genome and EST information for R. trigyna or closely related taxa, 46.2% of the unigenes did not match to known genes in the database. Based on sequence homology, 13552 unigenes were categorized into 44 Gene Ontology (GO) terms (Figure 3, Additional file 3). In each of the three main GO classifications, i.e., biological process, cellular component, and molecular function, “cellular process”, “metabolic process”, “cell”, “cell part” and “binding” terms were dominant among the returned terms. Many of the identified unigenes were classified in the categories “response to stimulus”, “organelle” and “catalytic activity”, whereas only a few genes belonged to “biological adhesion”, “cell killing”, “locomotion”, “viral reproduction”, “virion” and “virion part”. Interestingly, 894 unigenes were classified into the category “transporter activity” and 71 unigenes into the category “antioxidant activity”.

Figure 3

GO functional annotation. Unigenes with best BLAST hits were aligned to GO database. All 13552 unigenes were assigned to at least one GO term and were grouped into three main GO categories and 44 sub-categories. Right Y-axis represents number of genes in a category. Left Y-axis indicates percentage of a specific category of genes in each main category.

To further evaluate the integrity of our transcriptome library and the effectiveness of our annotation process, unigene sequences were subjected to Clusters of Orthologous Groups (COG) classification. Out of 35271 nr hits, 10968 sequences showed a COG classification (Figure 4, Additional file 3). Among the 25 COG categories, the cluster for “general function prediction only” (3179, 28.98%) was the largest group, followed by “replication, recombination and repair” (1681, 15.32%), “transcription” (1621, 14.78%), “signal transduction mechanisms” (1257, 11.46%), and “carbohydrate transport and metabolism” (1089, 9.93%). The categories “nuclear structure” (8, 0.07%) and “extracellular structures” (7, 0.06%) had the fewest corresponding genes.

Figure 4

COG classification of putative proteins. All 10968 putative proteins showing significant homology to those in COG database were functionally classified into 25 molecular families. Y-axis indicates number of unigenes in a specific functional cluster.

The 35271 annotated sequences were mapped to the reference canonical pathways in the Kyoto Encyclopedia of Genes and Genomes (KEGG). Among those, 15995 unigenes were assigned to 119 KEGG pathways (Additional file 3). The pathways most strongly represented by mapped unigenes were “metabolic pathways” (3459 unigenes), “biosynthesis of secondary metabolites” (1920 unigenes), and “plant-pathogen interaction” (1110 unigenes), indicating that these were the active pathways in R. trigyna. Finally, 38197 unigenes were annotated in the transcriptome of R. trigyna, of which 35459 were generated by BLAST analysis and 2738 were annotated using the ESTscan program.

Identification and annotation of potential differentially expressed genes (DEGs)

The analysis showed that 32697 unigenes were up-regulated and 31997 were down-regulated upon treatment of R. trigyna seedlings with salt solution. In total, 5032 DEGs were identified in both libraries, including 2370 up- and 2662 down-regulated unigenes (Figure 5). Among those DEGs responding to salt-stress treatment, 3719 were annotated and 1313 had no BLAST hits.

Figure 5

Identification of DEGs between C21 and T43. DEGs were determined using a threshold of FDR ≤ 0.001 and |log2Ratio| ≥ 1. Red spots represent up-regulated DEGs and green spots indicate down-regulated DEGs. Those shown in blue are unigenes that did not show obvious changes in salt-stressed R. trigyna.

Enrichment analysis was conducted to clarify the biological functions of the identified DEGs. The results indicated that 1947 DEGs were enriched in 27 GO terms, using a corrected P-value ≤ 0.05 as the threshold (Additional file 4). Among these GO categories, “oxidoreductase activity” (171 DEGs), “catalytic activity” (742 DEGs) and “response to stress” (155 DEGs) were significantly enriched among DEGs compared with the whole transcriptome background. In total, 2086 DEGs were enriched in 33 metabolic pathways (q-value ≤ 0.05) (Additional file 4). The most enriched metabolic pathways were “metabolic pathways” (483 DEGs), “biosynthesis of secondary metabolites” (351 DEGs), and “plant-pathogen interaction” (198 DEGs). The pathways of phenylpropanoid biosynthesis, flavonoid biosynthesis, stilbenoid, diarylheptanoid, and gingerol biosynthesis were very active, while those associated with ascorbate, aldarate, and nitrogen metabolism were strongly suppressed in salt-treated plants (Figure 6, Additional file 5).

Figure 6

KEGG enrichment analysis showed significantly enriched biosynthesis and metabolic pathways. Y-axis indicates numbers of DEGs enriched in KEGG pathways. X-axis represents DEGs enriched in secondary metabolites biosynthesis and metabolism pathways. Black bar represents up-regulated DEGs; red bar indicates down-regulated DEGs.

Genes related to ion transport

Among the unigenes, we identified 135, 72, and 25 genes related to regulation of K+ uptake, H+ pumping, and Na+ efflux, respectively (Figure 7, Table 2, Additional file 6). Clearly, genes associated with K+-transport formed the largest proportions of screened genes, suggesting that these genes played important role in reestablishing Na+/K+ homeostasis in salt-stressed R. trigyna. Among 135 K+ transport genes, 42 were DEGs, including 24 up- and 18 down-regulated genes in response to salt-stress. KUP (unigene14170, unigene26051) and CNGC (unigene3965, unigene7183) transcripts were present at high levels in both libraries, but were down-regulated by nearly 2-fold under salt-stress. Seven HKT1 genes, which are responsible for Na+ influx, were salt-responsive, of which three (unigene4890, unigene27586, and unigene20631) were suppressed by more than 2-fold. We identified 16, 44, and 12 genes encoding plasma membrane H+-ATPases (PM-H+-ATPases), vacuolar H+-ATPases (V-H+-ATPases) and H+-pyrophosphatases (V-H+-PPases), respectively. Out of 44 V-ATPase homologs, four showed down-regulated transcription under salt treatment while the other 40 did not. Transcripts of six PM-H+-ATPase homologs were highly abundant in both libraries. Two PM-H+-ATPases (unigene15650 and unigene248) were up-regulated by salt, showing increases in their Reads Per Kilobase of exon model per Million mapped reads (RPKM) values from 143.1 to 376.9 and from 79.8 to 240.3, respectively. Five PPases showed high transcript levels in both libraries, but their transcriptions did not change under salt-stress. Among 35 Na+ efflux genes, one PM-Na+/H+ antiporter gene (unigene798 or SOS1B) showed moderate transcript levels and two showed low transcript abundance. Transcription of unigene798 was slightly up-regulated by salt-stress (0.3-fold increase in abundance). Among 12 genes encoding V-Na+/H+ antiporter family proteins, four (unigene20634, unigene5272, unigene8445, and unigene752) were highly abundant. The transcript levels of unigene20634 and unigene5272 were up-regulated by nearly 2-fold, whereas those of unigene8445 and unigene752 were unchanged under salt-stress. There were abundant Nha (multisubunit Na+/H+ antiporter, unigene16859) transcripts in both libraries (RPKM values of 223.0 in C21 and 311.2 in T43).

Figure 7

Transcripts related to ion transport and ROS scavenging system. Columns at left of Y-axis show genes with down-regulated transcription under salt stress; columns at right of Y-axis show genes with up-regulated transcription under salt stress. Gene numbers and categories are shown next to columns. A: Transcriptional characteristics of genes related to ion transport. HKT, high-affinity K+ transporter; CNGC, cyclic nucleotide-gated channel; KUP, K+ ion transmembrane transporter; AKT, K+ channel; HAK, high affinity K+ transporter; CHX, cation H+ exchanger; KEA, K+ efflux antiporter; KOC, outward rectifying K+ channel; SOS1, salt overly sensitive 1; NHX, Na+/H+ exchanger; NhaP, NhaP-type Na+/H+ antiporter; V-ATPase, V-H+-ATPase; V-PPase, V-H+-PPase; P-ATPase, PM-H+-ATPase. B: Transcriptional characteristics of genes related to ROS scavenging system. GLR, glutaredoxin; APX, ascorbate peroxidase; MDAR, monodehydroascorbate reductase; DHAR, dehydroascorbate reductase; GR, glutathione reductase; GST, glutathione S-transferase; GPX glutathione peroxidase; POD, peroxidases; GLP, germin-like protein; CAT, catalase; PEX11, peroxisomal biogenesis factor11; Trx, thioredoxin; PrxR, peroxiredoxin; SOD, superoxide dismutase. GST is used in both the glutathione ascorbate cycle and GPX pathway.

Table 2 DEGs and highly transcribed genes related to ion transport *

Genes related to the ROS scavenging system

In total, 298 unigenes were predicted to encode enzymes related to reactive oxygen species (ROS) scavenging (Table 3, Additional file 7). There were 130, 84, 63, and 7 unigenes categorized into the GSH-ascorbate cycle, the peroxiredoxin/thioredoxin (PrxR/Trx) pathway, the glutathione peroxidase (GPX) pathway, and the catalase (CAT) pathway, respectively. The largest group was genes encoding thioredoxins (Trxs), glutathione S-transferases (GSTs), glutaredoxins (GLRs), and PODs (Figure 7). We detected 14 genes encoding SODs; three were up-regulated and 11 were down-regulated by salt treatment. In the GSH-ascorbate ROS removal pathway, seven ascorbate peroxidase (APX) genes (unigene44406, unigene47495, unigene17786, unigene13095, unigene57465, unigene25879, unigene10954), two GLR genes (unigene33573, unigene7889), one monodehydroascorbate reductase gene (unigene33304), one dehydroascorbate reductase gene (unigene42162), and two glutathione reductase genes (unigene8218, unigene3408) showed relatively abundant transcript levels (RPKM values of > 200). We identified five PrxRs, which play a role in the PrxR/Trx pathway in the antioxidant defense system. Two PrxR genes (unigene63275 and unigene61053) showed high transcript levels in both control and salt-treated plants (RPKM values of > 200). Out of 79 Trxs, five transcripts (unigene12493, unigene17142, unigene7940, unigene34328, and unigene2658) were expressed at high levels in both libraries. In the GPX pathway, there were 7 unigenes encoding GPXs but 53 unigenes encoding GSTs. Ten GST showed increased transcription under salt-stress, with RPKM values ranging from 102.6 to 477.2 in the T43 library but from 62.6 to 313.4 in the C21 library. In addition, 11 genes belonging to the Tau family of GSTs showed increased transcription under salt stress.

Table 3 DEGs and highly transcribed genes related to ROS scavenging system

Experimental validation

PCR amplification showed that all qPCR primers produced only single fragments of the expected lengths (155–350 bp, a 100% success rate). The qPCR results for 30 selected unigenes showed general agreement with their transcript abundance changes determined by RNA-seq, suggesting the reliability of the transcriptomic profiling data. However, moderate discrepancies between the expression levels and RPKM values were observed in ten unigenes, i.e., unigene7220, unigene20645, unigene5040, unigene65175, unigene55260, unigene65219, unigene42056, unigene19588, unigene41893 and unigene39 (Figure 8). On the other hand, amplifications of longer sequences (> 500 bp) had a lower success rate. Among the 30 unigenes, 26 unigenes yielded fragments > 500 bp and the other four failed to produce amplification products, yielding a success rate of 86.7%. All positive clones from the validation studies were sequenced by Sanger sequencing, and the results completely confirmed the Illumina results.

Figure 8

Expression pattern validation of selected genes by qPCR. Changes in transcript levels of 30 selected genes as detected by qPCR. A, B, and C represent 10 up-regulated, 10 down-regulated, and 10 unigenes with changes in transcript abundance of < 2-fold, respectively. X-axis shows -fold change in transcript abundance of unigenes. Black bar indicates transcript abundance changes calculated by the RPKM method. Red bar with associated standard error bar represents relative expression level determined by qPCR using 2-∆∆CT method. Results represent mean standard deviations (± SD) of three experimental replicates.


Construction of an informative transcriptome dataset for R. trigyna

For many non-model species, there is no background genomic information available for researchers to conduct comprehensive investigations into the genetic mechanisms underlying their unique features. Therefore, the newly developed NGS technique has been used widely to explore genomic solutions to important physiological questions. Gene annotation is an important element of NGS in which biological information is attached to the predicted genes or unigenes. A high proportion of unigenes with high-confidence BLASTX similarity to protein sequences from annotated gene catalogs of other plant species is considered to indicate the integrity of transcript sequences assembled from Illumina short-read data [23]. This assumption was also verified by the results generated in the present study. In R. trigyna, 54.27% of 65340 unigenes were annotated by BLAST analysis and functional bioinformatics analyses (e.g., GO, Swiss-Prot, and KEGG). Overall, the top five species with BLAST hits to annotated unigenes were Arabidopsis thaliana, Oryza sativa, Arabidopsis lyrata, Populus trichocarpa, and Vitis vinifera, species for which the annotations of their genomes are comprehensive and largely accepted. This suggested that the sequences of the R. trigyna unigenes generated in the present study were assembled and annotated correctly.

The lengths of the assembled sequences are crucial in determining the level of significance of a BLAST match [20]. Out of all the assembled unigenes, 66.79% had lengths between 300 and 500 bp, among which 37.1% had significant BLAST hits in the public databases. There were 1224 unigenes showing strong homology with the sequences hit in the database (E-value < 1.0E-50). Among 21697 unigenes with sequence lengths > 500 bp, 87.0% had significant BLAST scores and 59.2% showed E-values less than 1.0E-50. Out of all the unannotated unigenes (28991), 70.65% were longer than 500 bp, demonstrating that the lack of annotation for these unigenes was not because of a shorter sequence length but because of a genuine lack of hits to sequences in the database. Therefore, we can speculate that these unannotated unigenes represent a specific genetic resource of R. trigyna, which warrants further investigation.

In this study, we also used other strategies to enhance the effectiveness of the short reads assembly, apart from the use of bioinformatics methods. These included preparation of high-quality RNA samples, removal of dirty raw reads, BLAST assembly of unigenes using multiple databases, and the large sample population for sequencing [20, 24]. First, RNA was isolated from sterile R. trigyna seedlings to minimize the risk of contamination by foreign RNase. Second, three seedlings were used to extract RNA samples, not only to reduce sample bias, but also to ensure comprehensive coverage of the R. trigyna transcriptome. Last, two sequencing libraries (C21 and T43) were merged to generate longer sequences and to increase the sequencing depth. These strategies were an effective way to improve the quality of assembly and annotation of assembled unigenes [23].

In summary, an extensive and diverse expressed gene catalog, representing a large proportion of the transcribed genes in R. trigyna, was successfully sampled in the present study. The gene catalog provided a comprehensive understanding of the gene transcription profiles of R. trigyna, and laid a solid foundation for further study of salt-tolerance mechanisms and identification of novel genes in this species.

Transcriptome comparison identified genes related to salt-stress in R. trigyna

Gene transcription and/or expression is often compared among different developmental stages, among different plant organs, or among plants under different growth conditions [2022]. In this study, many genes showing transcriptional changes under salt stress were identified by comparing NaCl-treated seedlings with untreated controls. There were 64694 unigenes showing differences in transcript abundance, and 5032 were defined as DEGs using the thresholds of false discovery rate (FDR) ≤ 0.001 and |log2Ratio| ≥ 1. GO clustering analysis suggested the potential biological functions of these DEGs. For example, many DEGs were enriched in GO terms such as “oxidoreductase activity”, “catalytic activity” and “response to stress”. This information will be useful to elucidate salt-tolerance mechanisms and to find new salt-stress-related genes specific to R. trigyna. KEGG enrichment analysis identified significantly enriched metabolic pathways or signal transduction pathways involving DEGs. Surprisingly, genes in the phenylpropanoid biosynthesis and flavonoid biosynthesis pathways were transcribed at much high levels under salt stress than in the control in R. trigyna, which was not expected from the results of previous studies. This result suggested that genes in these two pathways may play vital roles in resisting oxidation and maintaining membrane integrity.

An outstanding advantage of NGS is the detection of transcripts present in low copy numbers. Among the 5032 identified DEGs, almost half were detected as low-abundance transcripts that were up-regulated after exposure to salt-stress. For example, only 14 reads from C21 could be mapped to unigene24841, encoding flavanone 3-hydroxylase (F3H), compared with 470 reads from T43. Nearly 30-fold up-regulation suggested that this gene may act as a key element in the salt-stress response. Overall, transcriptome comparison analysis provided sufficient information to study the salt-responsive mechanisms and genes related to salt-tolerance related in R. trigyna.

Ion transport genes are important for salt tolerance in R. trigyna

Ion transport is a crucial element in response to salt stress in plants. This is particularly true in halophytes [25]. It has been shown that halophytes are able to grow under extreme salinity conditions because of their anatomical and morphological adaptations and/or their avoidance mechanisms [25, 26]. R. trigyna is a representative recretohalophyte with unique morphological characteristics, such as salt glands and succulent leaves, allowing it to adapt to the salinized conditions of the Ordos desert. Our previous study showed that the salt glands of this species functionally excrete salt ions under normal and salt-stressed conditions. However, the amount of Na+ and Cl- excreted significantly increased under salt treatment. In addition, R. trigyna shows a strong ability to uptake K+ from barren soil. Therefore, it is possible that the ion transport mechanisms and their gene expression patterns in this plant differ from those of other plant species to some extent.

Once Na+ is taken up into the cell, ATPases (PM-H+-ATPases, V-H+-ATPases, and V-H+-PPases) are induced to create the driving force for Na+ transport. This results in not only extrusion of Na+ into the external environment by Na+/H+ antiporters in the plasma membrane, but also in compartmentalization into vacuoles by tonoplast Na+/H+ antiporters, which are essential for reestablishing cellular ion homeostasis in salt-stressed plants [25, 2730]. In this study, five V-H+-PPases were detected in R. trigyna, and were proposed to generate a proton electrochemical gradient. Taking the morphological characteristics of this plant into account, the largest share of the driving force generated by such enzymes is likely to be consumed by the succulent leaves and the salt glands, which include two mesophyll-like collecting cells with obvious vacuoles, and six secretory cells [31]. All mesophyll cells in leaves need V-H+-ATPases and V-H+-PPases to generate an electrochemical gradient to compartmentalize excess Na+ ions into the tonoplast, as a mechanism to cope with physiological drought conditions and ion toxicity [3, 25]. When Na+ is excreted via secretory cells, both V-H+-ATPases and V-H+-PPases are required to produce an H+ electrochemical gradient to pump Na+ into collecting cells. PM-H+-ATPases generate sufficient driving force to exchange Na+ ions out of the cells. In addition, higher H+-pumping activity may be indispensable to maintain the high concentration of K+ in salt-stressed R. trigyna. The activation of PM-H+-ATPases in salinized Populus euphratica led to repolarization or hyperpolarization of the plasma membrane and thus decreased NaCl-induced K+ loss through depolarization-activated outward rectifying K+ channels (DA-KOCs) [32, 33]. In R. trigyna, 10 out of 17 RtKOCs showed high homology with those of Populus, demonstrating similar functions of these genes. The analysis of transcripts related to K+ uptake showed that the outstanding potassium uptake capacity of R. trigyna was probably conferred by enhanced KUP/HAK (potassium transporter) and AKT K+ uptake systems, among which the number of up-regulated unigenes (39) far exceeded that of down-regulated unigenes (16) under salt stress. In summary, proton electrochemical gradient-dependent K+ maintenance, ions secretion, compartmentalization, and enhanced K+ uptake systems may represent important components in reestablishing ion homeostasis in R. trigyna.

On the other hand, we detected 12 NHXs that were probably responsible for Na+ sequestration. The significantly up-regulated unigene20634 and 5272 showed high homology to the Citrus reticulata NHX1 and a Salsola komarovii Na+/H+ antiporter, respectively. The slightly down-regulated unigene8445 and 752 had a highly significant BLAST hit to the Na+/H+ antiporter of M. crystallinum and the NHX1 of Tetragonia tetragonioides, respectively. These four unigenes also showed strong homology to the NHX2 of A. thaliana, suggesting that these AtNHX2-like proteins play an important role in avoiding or mitigating the deleterious effects of high Na+ levels in the cytosol or in regulating intravacuolar K+ and pH, which has been demonstrated in Arabidopsis[34]. Our prediction is also supported to some degree by the studies of succulent leaves of the halophytes Suaeda salsa and Halostachys caspica[35].

In R. trigyna, we identified only one moderately expressed and up-regulated PM Na+/H+ antiporter (SOS1B, unigene798), named RtSOS1. The transcript profile of this gene was consistent with those of ThSOS1 and PeSOS1, which remained relatively constant or were slightly up-regulated under salt-tress [36, 37]. Escherichia coli complementation experiments with PeNhaD1 rescued salt-exposed bacteria by lowering salt accumulation [38], suggesting that the P. euphratica gene PeNhaD1 also functions as an Na+⁄H+ antiporter. Interestingly, a highly abundant transcript (unigene16859, RtNha1) encoding an Nha protein was identified in our dataset. This may play a role in Na+ exclusion in R. trigyna. At present, it is still unknown how the moderately up-regulated RtSOS1 achieves effective efflux of excess salt ions under extreme salt-stress, and whether any other gene products are also involved in this process. This requires further investigations in the near future.

ROS scavenging plays a key role in salt-stress response in R. trigyna

Salt stress causes a rapid increase in ROS, including superoxide radicals (·O2-), hydrogen peroxide (H2O2), and hydroxyl radicals (·OH), which perturb cellular redox homeostasis and result in oxidative damage to many cellular components and structures [39, 40]. According to our previous studies, ·O2- increased rapidly in R. trigyna seedlings treated with 400 mM NaCl, as did the concentrations of other antioxidants, for example, GSH. At the same time, increased activities of antioxidant enzymes (SOD, POD) were also detected. This suggested that R. trigyna probably has a similar ROS scavenging system to that of other plant species, and that this system can be enhanced to increase its antioxidative ability. As the first line of defense against oxidative damage, SODs are usually induced by salinity to rapidly dismutate ·O2- into oxygen and H2O2, which is subsequently removed through various pathways [41]. In this study, 14 SODs were identified. Two of them (unigene23472, unigene8340) were transcribed at high levels, and showed up-regulated transcription under salt stress. Therefore, they were probably associated with enhanced SOD activity and responsible for converting the increased ·O2- into H2O2. On the other hand, increased GSH could be consumed in the ascorbate-GSH cycle and the GPX cycle [39]. In the ascorbate-GSH cycle, the oxidized ascorbate (i.e., monodehydroascorbate) can be converted into dehydroascorbate (DHA); DHA is then reduced to ascorbate at the expense of GSH, yielding oxidized GSH (i.e., glutathione disulfide). In the GPX pathway, GPX can reduce H2O2 to the corresponding hydroxyl compounds using GSH. However, to fulfill its function as an antioxidant, GSH must be catalyzed by GSTs, which have GPX activity and can use GSH to reduce organic hydroperoxides of fatty acids and nucleic acids to the corresponding monohydroxy alcohols [42]. Interestingly, in our dataset, there were 35 salt-induced GSTs genes that likely participated in such catalytic processes, especially the 11 Tau family GSTs [43]. This finding was similar to the results of other studies on salt-stressed plants [44, 45]. Therefore, the interaction between GSH and GST may be an important factor in the ROS scavenging system of R. trigyna. In addition, we identified 84 genes encoding components of the PrxR/Trx defense system, which employs a thiol-based catalytic mechanism to reduce H2O2 and is regenerated using Trxs as electron donors [46].


This study describes a platform in which publicly available information was matched to the R. trigyna transcriptome using NGS technology. The substantially assembled sequences represented a considerable portion of the transcriptome of this plant. Transcriptome comparison identified the transcripts and metabolic pathways that play significant roles in the response to salt stress. These results will prompt studies on the molecular mechanisms of salt resistance, facilitate the analysis of expression profiles of genes related to salt tolerance, and accelerate the discovery of specific stress-resistance genes in R. trigyna. The information provided here can also further extend the knowledge on the salt tolerance of halophytes that survive in high sodic stress in areas such as the semi-desert saline areas in Ordos, Inner Mongolia, China.


Plant materials and experimental treatment

Seeds of R. trigyna were collected in September 2009 from the Eastern Alxa-Western Ordos area in Inner Mongolia, China. Intact plump seeds were selected, immersed in 10% sodium hypochlorite for 15 min, and then rinsed three times with sterilized double-distilled water. The seeds were germinated in a 150-ml conical flask containing 40 ml MS medium in the dark for 72 h. Germinated seeds were then grown on the same medium at 25°C under 70% relative humidity and a 16-h light/8-h dark photoperiod for 15 d. When seedlings were approximately 10 cm high, they were transferred into a test tube containing 50 ml ½-strength Hoagland’s medium, and further cultured for another 4 weeks with a change of medium every 2 days. Three healthy seedlings of a similar size were selected for further analysis. Prior to NaCl treatment, leaf samples were collected for RNA extraction from each seedling as controls. The roots were then immersed in Hoagland’s medium containing NaCl. The NaCl concentration was increased stepwise from 100 mM to 400 mM at a rate of 100 mM per step every 8 h. Leaves were collected before each increase in NaCl concentration. Leaves were also collected after the seedlings were treated in 100 mM NaCl for 0.5 h, 1 h, 2 h, 4 h, and 8 h. The leaf samples were immediately snap-frozen in liquid nitrogen and stored at −80°C until analysis.

RNA preparation

Total RNA was extracted with Plant Plus RNA reagent (DP437, Tiangen, Beijing) according to the manufacturer’s instructions. The extracted RNA was treated with RNase-free DNase I (TaKaRa Bio Inc., Otsu, Shiga, Japan) for 45 min at 37°C to remove residual DNA. RNA quality was evaluated using the Agilent 2100 Bioanalyzer with a minimum RNA integrity number value of 8. In this study, two sequencing libraries were prepared, the C21 library from control samples and the T43 library from salt-stressed samples.

cDNA library construction and sequencing

Beads with Oligo (dT) (Illumina) were used to isolate poly (A) + RNA from 20 μg of each RNA pool. Fragmentation buffer (Ambion, Austin, TX, USA) was added to break mRNA into short fragments. First-strand cDNA was synthesized using N6 random hexamers (Illumina Inc., San Diego, CA, USA) and SuperScript® III reverse transcriptase (Invitrogen, Grand Island, NY, USA). The second-strand cDNA was synthesized with DNA polymerase I (Invitrogen, Grand Island, NY, USA). Short fragments were purified with a QIAquick PCR extraction kit (Qiagen, Hilden, Germany) and resolved with EB buffer for end-repair and poly (A) addition, and then linked to sequencing adapters. Agarose gel electrophoresis and PCR amplification were used to select suitable fragments. Two cDNA libraries were constructed and sequenced on the Illumina HiSeq™ 2000 platform. The sequences were base-called and quality-checked by the Illumina data processing pipeline.

De novoassembly of sequences and sequence clustering

The raw reads were filtered to obtain high-quality clean reads prior to assembly. This was performed by removing adaptor sequences, duplicated sequences, reads containing more than 5% “N” (i.e., ambiguous bases in reads), and reads in which more than 50% of the bases showed a Q-value ≤ 5. Transcriptome de novo assembly was carried out with the short reads assembly program, SOAPdenovo, with the default settings, except that the K-mer value was 25-mer [47]. The different K-mer sizes were assessed; the 29-mer yielded the best assembly for the desired application and so it was used to construct the de Bruijn graph. Contigs without ambiguous base reads were obtained by conjoining K-mers in an unambiguous path. The reads were then mapped to contigs to construct scaffolds with paired-end information. Paired-end reads were used again for gap filling of scaffolds to obtain unigenes. To evaluate the depth of sequence coverage, all usable reads were realigned to the unigenes with SOAPaligner [48]. To minimize sequence redundancy, the scaffolds were clustered using the GICT [49]. The clustering output was passed to the sequence clustering software CAP3 assembler [50] for multiple alignment and generation of consensus sequences. The sequences that did not reach the threshold set and did not fall into any assembly were assigned as singletons.

Bioinformatics analysis of sequencing data

To assign predicted gene descriptions for the assembled unigenes, they were aligned against the plant protein dataset of nr, Swiss-Prot/Uniprot protein database, and COG databases, respectively, using BLASTX with a significance threshold of E-value ≤ 10-5. To identify the best BLASTX hits from the alignments, putative gene names, ‘CDS’, and predicted proteins of the corresponding assembled sequences were produced. The orientation of sequences was also derived from BLAST annotations. For sequences aside from those obtained from BLAST searches, the EST Scan program was used to predict the ‘CDS’ and its orientation. Functional categorization by GO terms [51] was performed based on the best BLASTX hits from the nr database using BLAST2GO software [52] according to molecular function, biological process, and cellular component ontologies with an E-value threshold of 10-5. The data were statistically analyzed using WEGO software [53]. The pathway assignments were carried out by sequence searches against the KEGG database [54], also using the BLASTX algorithm with an E-value threshold of 10-5.

Identification and functional annotation of DEGs

To identify DEGs between the two samples, the transcript abundance of all assembled unigenes were analyzed using the RPKM method [55]. With reference to “the significance of digital gene expression profiles” [56], a rigorous algorithm has been developed to identify DEGs between two samples. The probability of gene A being transcribed equally between two samples was calculated using the following formula:

p y | x = N 2 N 1 y x + y ! x ! y ! 1 + N 2 N 1 x + y + 1

where the P-value corresponds to the differential gene expression test; N1 and N2 represent the total clean tag number of C21 and T43 samples, respectively; and x and y denote the tag numbers of the gene of interest present in C21 and T43, respectively.

The FDR is a method to determine the threshold P-value in multiple tests and analysis through manipulating the FDR value. FDR ≤ 0.001 and the absolute value of log2Ratio ≥ 1 were used as the threshold to identify DEG. Functional enrichment analysis including GO and KEGG were performed using the following ultra-geometric test to find which DEGs were significantly enriched in GO terms (P-value ≤ 0.05) and metabolic pathways (q-value ≤ 0.05) compared with the whole transcriptome background.

Ρ = 1 i = 0 m 1 M i N M n i N n

Validation and expression pattern analysis

To experimentally validate the transcriptional abundance results from sequencing and computational analysis, 30 unigenes (10 up-regulated, 10 down-regulated, and 10 unigenes with no significant changes in transcript abundance) were selected for RT-PCR and qPCR analysis. Reverse transcription reactions were performed using SuperScript® III Reverse Transcriptase (Invitrogen, Grand Island, NY, USA) with approximately 5 μg total RNA following the manufacturer’s instructions. Primers for RT-PCR and qPCR were designed using Primer Premier 5 and Beacon Designer 7.0 software, respectively (shown in Additional file 8). β-Actin was used as the internal control gene. qPCR was performed on a Qiagen Rotor-gene Q realtime PCR platform (Qiagen, Hilden, Germany) using SYBR-Green real-time PCR mix (TaKaRa Bio Inc., Otsu, Shiga, Japan) to detect transcript abundance. The amplification was achieved by the following PCR protocol: first denaturation 95°C for 30 s, then 40 cycles of denaturation at 95°C for 5 s, annealing and extension at 55°C for 30 s. The relative expression levels of the selected unigenes normalized to β-Actin was calculated using 2-∆∆Ct method. All reactions were performed with three replicates, and data were analyzed using Rotor-gene Q series software.



Next-Generation Sequencing


Differentially Expressed Genes


Reads Per Kilobase of exon model per Million mapped reads


False Discovery Rate


Gene ontology


Kyoto Encyclopedia of Genes and Genomes


Clusters of Orthologous Groups


Potassium Channel


High-Affinity K+ Transporter


Cyclic Nucleotide-Gated Channel


K+ Ion Transmembrane Transporter


High Affinity K+ Transporter


Cation H+ Exchanger


K+ Efflux Antiporter


Outward Rectifying K+ Channel


Salt Overly Sensitive 1


Na+/H+ Exchanger


Nhap-Type Na+/H+ Antiporter


V-H+-ATPase, vacuolar H+-ATPase


V-H+-PPase, H+-Pyrophosphatase


PM-H+-ATPase, Plasma Membrane H+-ATPase


Reactive Oxygen Species




Superoxide Radicals, H2O2, Hydrogen Peroxide


Hydroxyl Radicals




Ascorbate Peroxidase


Monodehydroascorbate Reductase


Dehydroascorbate Reductase


Glutathione Reductase


Glutathione S-Transferase


Glutathione peroxidase




Germin-Like Protein




Peroxisomal Biogenesis Factor11






Superoxide Dismutase

A. lyrata:

Arabidopsis lyrata subsp. lyrata

A. thaliana:

Arabidopsis thaliana

A. hypogaea:

Arachis hypogaea

B. gymnorhiza:

Bruguiera gymnorhiza

C. sinensis:

Camellia sinensis

C. chinense:

Capsicum chinense

C. quinoa:

Chenopodium quinoa

C. limon:

Citrus limon

C. unshiu:

Citrus unshiu

C. reticulata:

Citrus reticulata

C. sativus:

Cucumis sativus

D. carota:

Daucus carota

E. guineensis:

Elaeis guineensis

E. camaldulensis:

Eucalyptus camaldulensis

G. max:

Glycine max

G. hirsutum:

Gossypium hirsutum

H. caspica:

Halostachys caspica

H. brasiliensis:

Hevea brasiliensis

H. vulgare:

Hordeum vulgare

I. batatas:

Ipomoea batatas

I. lacteal:

Iris lacteal

J. curcas:

Jatropha curcas

K. foliatum:

Kalidium foliatum

L. bicolor:

Limonium bicolor

L. chinensis:

Litchi chinensis

M. pumila:

Malus pumila

M. pusilla:

Malva pusilla

M. truncatula:

Medicago truncatula

M. crystallinum:

Mesembryanthemum crystallinum

N. nucifera:

Nelumbo nucifera

N. alata:

Nicotiana alata

N. benthamiana:

Nicotiana benthamiana

N. tabacum:

Nicotiana tabacum

O. sativa:

Oryza sativa

P. acutifolius :

Phaseolus acutifolius

P. sativum:

Pisum sativum

P. euphratica:

Populus euphratica

P. trichocarpa:

Populus trichocarpa

P. juliflora:

Prosopis juliflora

P. persica:

Prunus persica

P. tenuiflora:

Puccinellia tenuiflora

P. communis:

Pyrus communis

R. australe:

Rheum australe

R. communis:

Ricinus communis

S. komarovii:

Salsola komarovii

S. lycopersicum:

Solanum lycopersicum

S. nigrum:

Solanum nigrum

S. tuberosum:

Solanum tuberosum

S. oleracea:

Spinacia oleracea

S. salsa:

Suaeda salsa

T. androssowii:

Tamarix androssowii

T. hispida:

Tamarix hispida

T. tetragonioides:

Tetragonia tetragonioides

T. octandrum:

Trichostigma octandrum

V. radiate:

Vigna radiate

V. unguiculata:

Vigna unguiculata

V. vinifera:

Vitis vinifera

Z. aethiopica:

Zantedeschia aethiopica

Z. mays:

Zea mays

Z. jujuba:

Ziziphus jujuba

Z. marina Zostera marina:



  1. 1.

    Allakhverdiev SI, Sakamoto A, Nishiyama Y, Inaba M, Murata N: Ionic and osmotic effects of NaCl-induced inactivation of photosystems I and II in Synechococcus sp. Plant Physiol. 2000, 123: 1047-1056. 10.1104/pp.123.3.1047.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  2. 2.

    Zhu JK: Plant salt tolerance. Trends Plant Sci. 2001, 6: 66-71.

    Article  CAS  PubMed  Google Scholar 

  3. 3.

    Tuteja N: Mechanisms of high salinity tolerance in plants. Method Enzymol. 2007, 428: 419-438.

    Article  CAS  Google Scholar 

  4. 4.

    Yeo A: Predicting the interaction between the effects of salinity and climate change on crop plants. Sci Hortic-amsterdam. 1998, 78: 159-174. 10.1016/S0304-4238(98)00193-9.

    Article  Google Scholar 

  5. 5.

    Blumwald E, Aharon GS, Apse MP: Sodium transport in plant cells. BBA-Biomembranes. 2000, 1465: 140-151. 10.1016/S0005-2736(00)00135-8.

    Article  CAS  PubMed  Google Scholar 

  6. 6.

    Conde A, Chaves MM, Geros H: Membrane transport, sensing and signaling in plant adaptation to environmental stress. Plant Cell Physiol. 2011, 52: 1583-1602. 10.1093/pcp/pcr107.

    Article  CAS  PubMed  Google Scholar 

  7. 7.

    Breckle SW: Salinity, halophytes and salt affected natural ecosystems. Salinity: environment-plants-molecules. Edited by: Läuchli A, Lüttge U. 2004, Netherlands: Springer, 53-77.

    Google Scholar 

  8. 8.

    Türkan I, Demiral T: Recent developments in understanding salinity tolerance. Environ Exp Bot. 2009, 67: 2-9. 10.1016/j.envexpbot.2009.05.008.

    Article  Google Scholar 

  9. 9.

    Flowers TJ, Colmer TD: Salinity tolerance in halophytes. New Phytol. 2008, 179: 945-963. 10.1111/j.1469-8137.2008.02531.x.

    Article  CAS  Google Scholar 

  10. 10.

    Zhao YZ: Vascular plants in Ordos Plateau. 2006, Hohhot: Inner Mongolia University Press

    Google Scholar 

  11. 11.

    Ma YQ: Flora Intramongolica (Tomus 3). 1989, Hohhot: Typis Intramongolicae Popularis

    Google Scholar 

  12. 12.

    Delectis Florae Reipublicae Popularis Sinicae Agendae Academiae Sinicae Edita: Flora Reipublicae Popularis Sinicae (Tomus 50). 1988, Beijing: Science Press

    Google Scholar 

  13. 13.

    Yang C, Wang YC, Liu Q, Zhang YF, Zhang YJ: Conservation Biology of Tetraena mongolia. 2002, Beijing: Science Press

    Google Scholar 

  14. 14.

    Li B: Natural resources and environment study of Inner Mongolia Ordos Plateau. 1990, Beijing: Science Press

    Google Scholar 

  15. 15.

    Xue Y: PhD thesis. Studies on salt tolerance mechanisms of Reaumuria trigyna Maxim. native to Alxa. 2008, Inner Mongolia University: College of Life Sciences

    Google Scholar 

  16. 16.

    Schuster SC: Next-generation sequencing transforms today’s biology. Nat Methods. 2008, 5: 16-18. 10.1038/nmeth1156.

    Article  CAS  PubMed  Google Scholar 

  17. 17.

    Mardis ER: The impact of next-generation sequencing technology on genetics. Trends Genet. 2008, 24: 133-141. 10.1016/j.tig.2007.12.007.

    Article  CAS  PubMed  Google Scholar 

  18. 18.

    Metzker ML: Sequencing technologies-the next generation. Nat Rev Genet. 2009, 11: 31-46.

    Article  Google Scholar 

  19. 19.

    Bräutigam A, Gowik U: What can next generation sequencing do for you? Next generation sequencing as a valuable tool in plant research. Plant Biology. 2010, 12: 831-841. 10.1111/j.1438-8677.2010.00373.x.

    Article  PubMed  Google Scholar 

  20. 20.

    Shi CY, Yang H, Wei CL, Yu O, Zhang ZZ, Jiang CJ, Sun J, Li YY, Chen Q, Xia T: Deep sequencing of the Camellia sinensis transcriptome revealed candidate genes for major metabolic pathways of tea-specific compounds. BMC Genomics. 2011, 12: 131-10.1186/1471-2164-12-131.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  21. 21.

    Zenoni S, Ferrarini A, Giacomelli E, Xumerle L, Fasoli M, Malerba G, Bellin D, Pezzotti M, Delledonne M: Characterization of transcriptional complexity during berry development in Vitis vinifera using RNA-Seq. Plant Physiol. 2010, 152: 1787-1795. 10.1104/pp.109.149716.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  22. 22.

    Li PH, Ponnala L, Gandotra N, Wang L, Si YQ, Tausta SL, Kebrom TH, Provart N, Patel R, Myers CR, Reidel EJ, Turgeon R, Liu P, Sun Q, Nelson T, Brutnell TP: The developmental dynamics of the maize leaf transcriptome. Nat Genet. 2010, 42: 1060-1067. 10.1038/ng.703.

    Article  CAS  PubMed  Google Scholar 

  23. 23.

    Hao DC, Ge GB, Xiao PG, Zhang YY, Yang L: The first insight into the tissue specific Taxus transcriptome via Illumina second generation sequencing. PLoS One. 2011, 6: e21220-10.1371/journal.pone.0021220.

    PubMed Central  Article  CAS  Google Scholar 

  24. 24.

    Logacheva MD, Kasianov AS, Vinogradov DV, Samigullin TH, Gelfand MS, Makeev VJ, Penin AA: De novo sequencing and characterization of floral transcriptome in two species of buckwheat (Fagopyrum). BMC Genomics. 2011, 12: 30-10.1186/1471-2164-12-30.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  25. 25.

    Zhu JK: Regulation of ion homeostasis under salt stress. Curr Opin Pant Biol. 2003, 6: 441-445. 10.1016/S1369-5266(03)00085-2.

    Article  CAS  Google Scholar 

  26. 26.

    Flowers TJ, Hajibagheri MA, Clipson NJW: Halophytes. Q Rev Biol. 1986, 61: 313-337. 10.1111/j.1469-185X.1986.tb00657.x.

    Article  Google Scholar 

  27. 27.

    Blumwald E: Sodium transport and salt tolerance in plants. Curr Opin Cell Biol. 2000, 12: 431-434. 10.1016/S0955-0674(00)00112-5.

    Article  CAS  PubMed  Google Scholar 

  28. 28.

    Parida AK, Das AB: Salt tolerance and salinity effects on plants: a review. Ecotoxicol Environ Saf. 2005, 60: 324-349. 10.1016/j.ecoenv.2004.06.010.

    Article  CAS  PubMed  Google Scholar 

  29. 29.

    Chinnusamy V, Zhu JK: Salt stress signaling and mechanisms of plant salt tolerance. Genet Eng. 2006, 27: 141-177. 10.1007/0-387-25856-6_9.

    Article  CAS  Google Scholar 

  30. 30.

    Niu XM, Bressan RA, Hasegawa PM, Pardo JM: Ion homeostasis in NaCl stress environments. Plant Physiol. 1995, 109: 735-742.

    PubMed Central  CAS  PubMed  Google Scholar 

  31. 31.

    Xue Y, Wang YC: Study on characters of ions secretion from Reaumuria trigyna. J Desert Res. 2008, 28: 437-442.

    Google Scholar 

  32. 32.

    Chen S, Polle A: Salinity tolerance of Populus. Plant Biology. 2010, 12: 317-333.

    Article  CAS  PubMed  Google Scholar 

  33. 33.

    Sun J, Dai SX, Wang RG, Chen SL, Li NY, Zhou XY, Lu CF, Shen X, Zheng XJ, Hu ZM, Zhang ZK, Song J, Xu Y: Calcium mediates root K+/Na+ homeostasis in poplar species differing in salt tolerance. Tree Physiol. 2009, 29: 1175-1186. 10.1093/treephys/tpp048.

    Article  CAS  PubMed  Google Scholar 

  34. 34.

    Bassil E, Tajima H, Liang YC, Ohto M, Ushijima K, Nakano R, Esumi T, Coku A, Belmonte M, Blumwald E: The Arabidopsis Na+/H+ antiporters NHX1 and NHX2 control vacuolar pH and K+ homeostasis to regulate growth, flower development, and reproduction. Plant Cell. 2011, 23: 3482-3497. 10.1105/tpc.111.089581.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  35. 35.

    Guan B, Hu YZ, Zeng YL, Wang Y, Zhang FC: Molecular characterization and functional analysis of a vacuolar Na+/H+ antiporter gene (HcNHX1) from Halostachys caspica. Mol Biol Rep. 2010, 38: 1889-1899.

    Article  PubMed  Google Scholar 

  36. 36.

    Wu YX, Ding N, Zhao X, Zhao MG, Chang ZQ, Liu JQ, Zhang LX: Molecular characterization ofPeSOS1: the putative Na+/H+ antiporter of Populus euphratica. Plant Mol Biol. 2007, 65: 1-11. 10.1007/s11103-007-9170-y.

    Article  CAS  PubMed  Google Scholar 

  37. 37.

    Taji T, Seki M, Satou M, Sakurai T, Kobayashi M, Ishiyama K, Narusaka Y, Narusaka M, Zhu JK, Shinozaki K: Comparative genomics in salt tolerance between Arabidopsis and Arabidopsis-related halophyte cress using Arabidopsis microarray. Plant Physiol. 2004, 135: 1697-1709. 10.1104/pp.104.039909.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  38. 38.

    Ottow EA, Polle A, Brosche M, Kangasjärvi J, Dibrov P, Zörb C, Teichmann T: Molecular characterization of PeNhaD1: the first member of the NhaD Na+/H+ antiporter family of plant origin. Plant Mol Biol. 2005, 58: 75-88. 10.1007/s11103-005-4525-8.

    Article  CAS  PubMed  Google Scholar 

  39. 39.

    Mittler R: Oxidative stress, antioxidants and stress tolerance. Trends Plant Sci. 2002, 7: 405-410. 10.1016/S1360-1385(02)02312-9.

    Article  CAS  PubMed  Google Scholar 

  40. 40.

    Miller G, Shulaev V, Mittler R: Reactive oxygen signaling and abiotic stress. Physiol Plant. 2008, 133: 481-489. 10.1111/j.1399-3054.2008.01090.x.

    Article  CAS  PubMed  Google Scholar 

  41. 41.

    Zhang H, Han B, Wang T, Chen SX, Li HY, Zhang YH, Dai SJ: Mechanisms of plant salt response: insights from proteomics. J Proteome Res. 2012, 11: 49-67. 10.1021/pr200861w.

    Article  PubMed  Google Scholar 

  42. 42.

    Dixon DP, Lapthorn A, Edwards R: Plant glutathione transferases. Genome Biol. 2002, 3: 3004.1-3004.10.

    Article  Google Scholar 

  43. 43.

    Kilili KG, Atanassova N, Vardanyan A, Clatot N, Al-Sabarna K, Kanellopoulos PN, Makris AM, Kampranis SC: Differential roles of Tau class glutathione S-transferases in oxidative stress. J Biol Chem. 2004, 279: 24540-24551. 10.1074/jbc.M309882200.

    Article  CAS  PubMed  Google Scholar 

  44. 44.

    Pang QY, Chen SX, Dai SJ, Chen YZ, Wang Y, Yan XF: Comparative proteomics of salt tolerance in Arabidopsis thaliana and Thellungiella halophila. J Proteome Res. 2010, 9: 2584-2599. 10.1021/pr100034f.

    Article  CAS  PubMed  Google Scholar 

  45. 45.

    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-3958. 10.1021/pr200535r.

    Article  CAS  PubMed  Google Scholar 

  46. 46.

    Dietz KJ: Peroxiredoxins in plants and cyanobacteria. Antioxid Redox Sign. 2011, 15: 1129-1159. 10.1089/ars.2010.3657.

    Article  CAS  Google Scholar 

  47. 47.

    Li RQ, Zhu HM, Ruan J, Qian WB, Fang XD, Shi ZB, Li YR, Li ST, Shan G, Kristiansen K, Li SG, Yang HM, Wang J, Wang J: De novo assembly of human genomes with massively parallel short read sequencing. Genome Res. 2010, 20: 265-272. 10.1101/gr.097261.109.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  48. 48.

    Li RQ, Li YR, Kristiansen K, Wang J: SOAP: short oligonucleotide alignment program. Bioinformatics. 2008, 24: 713-714. 10.1093/bioinformatics/btn025.

    Article  CAS  PubMed  Google Scholar 

  49. 49.

    Pertea G, Huang XQ, Liang F, Antonescu V, Sultana R, Karamycheva S, Lee Y, White J, Cheung F, Parvizi B, Tsai J, Quackenbush J: TIGR Gene Indices clustering tools (TGICL): a software system for fast clustering of large EST datasets. Bioinformatics. 2003, 19: 651-652. 10.1093/bioinformatics/btg034.

    Article  CAS  PubMed  Google Scholar 

  50. 50.

    Huang XQ, Madan A: CAP3: A DNA sequence assembly program. Genome Res. 1999, 9: 868-877. 10.1101/gr.9.9.868.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  51. 51.

    Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G: Gene Ontology: tool for the unification of biology. Nat Genet. 2000, 25: 25-29. 10.1038/75556.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  52. 52.

    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-3676. 10.1093/bioinformatics/bti610.

    Article  CAS  PubMed  Google Scholar 

  53. 53.

    Ye J, Fang L, Zheng HK, Zhang Y, Chen J, Zhang ZJ, Wang J, Li ST, Li RQ, Bolund L, Wang J: WEGO: a web tool for plotting GO annotations. Nucleic Acids Res. 2006, 34 (suppl 2): W293-W297.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  54. 54.

    Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, Katayama T, Kawashima S, Okuda S, Tokimatsu T, Yamanishi Y: KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008, 36 (suppl 1): D480-D484.

    PubMed Central  CAS  PubMed  Google Scholar 

  55. 55.

    Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B: Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008, 5: 621-628. 10.1038/nmeth.1226.

    Article  CAS  PubMed  Google Scholar 

  56. 56.

    Audic S, Claverie JM: The significance of digital gene expression profiles. Genome Res. 1997, 7: 986-995.

    CAS  PubMed  Google Scholar 

Download references


We thank the Beijing Genomics Institute at Shenzhen for Illumina sequencing and raw data analysis. This work was supported by the National Natural Science Foundation of China (30760052, and 31140020), the Natural Science Foundation key project of Inner Mongolia Autonomous Region (200711020501, and 2012ZD05), and the Excellent Doctoral Dissertation training program of Inner Mongolia University.

Author information



Corresponding author

Correspondence to Ying-chun Wang.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

ZHD cultivated plants; performed RNA extraction, RT-PCR, and qPCR; analyzed, interpreted, and submitted data; and drafted the manuscript. LLZ and ZQ participated in the design of the study and data analysis. JW and ZG participated in plant management and data submission to the database. SBW participated in analysis and interpretation of data, and critically revised the manuscript. YCW conceived of the study, participated in its design and data interpretation, and revised the manuscript critically. All authors read and approved the final manuscript.

Electronic supplementary material

Additional file 1: Deposited C21-unigenes and TSA accession numbers. Sequences with no gap and length > 200 bp were selected from sample C21 and submitted to TSA database. (TXT 1 MB)

Additional file 2: Deposited T43-unigenes and the TSA accession numbers. Sequences with no gap and length > 200 bp were selected from sample T43 and submitted to TSA database. (TXT 1 MB)

Functional annotation of all assembled unigenes, including GO, COG, and KEGG analyses.

Additional file 3: All-unigene sequences were searched against protein databases (nr, KEGG, COG) using BLASTX (E-value ≤ 10-5). (XLS 2 MB)

Summary and functional annotation of identified DEGs. Unigenes with absolute value of |log

Additional file 4: 2 Ratio| ≥ 1 and FDR ≤ 0.001 were identified as DEGs. GO and KEGG analysis of DEGs were based on cut-off E-value of ≤ 10-5. (XLS 106 KB)

Additional file 5: Summary of DEGs enriched in KEGG pathways. Pathways and backbone gene numbers are shown in table. The q-value of all these pathways was ≤ 0.05. (DOCX 25 KB)

Additional file 6: Transcripts related to ion transport. (XLS 99 KB)

Additional file 7: ROS scavenging system-associated transcripts. (XLS 122 KB)

Additional file 8: Primers used for experimental validation.(XLSX 12 KB)

Authors’ original submitted files for images

Rights and permissions

Reprints and Permissions

About this article

Cite this article

Dang, Z., Zheng, L., Wang, J. et al. Transcriptomic profiling of the salt-stress response in the wild recretohalophyte Reaumuria trigyna. BMC Genomics 14, 29 (2013).

Download citation


  • Reaumuria trigyna
  • Recretohalophyte
  • Transcriptome
  • Illumina sequencing
  • Salt-stress response