International Journal of Genomics
Volume 2015 (2015), Article ID 265791, 8 pages
http://dx.doi.org/10.1155/2015/265791

Research Article

Analysis of the Thinopyrum elongatum Transcriptome under Water Deficit Stress

Yongjun Shu, Jun Zhang, You Ao, Lili Song, and Changhong Guo

Key Laboratory of Molecular Cytogenetics and Genetic Breeding of Heilongjiang Province, College of Life Science and Technology, Harbin Normal University, Heilongjiang 150025, China

Received 5 November 2014; Revised 13 January 2015; Accepted 20 January 2015

Academic Editor: Ferenc Olasz

Copyright © 2015 Yongjun Shu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The transcriptome of Thinopyrum elongatum under water deficit stress was analyzed using RNA-Seq technology. The results showed that genes involved in processes of amplification of stress signaling, reductions in oxidative damage, creation of protectants, and roots development were expressed differently, which played an important role in the response to water deficit. The Th. elongatum transcriptome research highlights the activation of a large set of water deficit-related genes in this species and provides a valuable resource for future functional analysis of candidate genes in the water deficit stress response.

1. Introduction

Water deficit is responsible for the greatest crop losses worldwide and is expected to worsen, heightening international interest in drought tolerance in crops [1]. Plant adaptation to water deficit is the result of many different physiological and molecular mechanisms that interact in a complex manner. Previous studies have shown that the plant response to water deficit stress involves numerous genes, which activate a series of physiological and biochemical processes to counteract the effects of the water-limited environment, including (1) the synthesis and accumulation of various osmoprotectants, (2) maintaining intracellular ion homeostasis via the expression of transporters, and (3) scavenging of reactive oxygen species (ROS) generated as a secondary effect of water deficit by detoxification enzymes [2, 3].

In addition to these physiological and biochemical effects, regulatory systems that link the sensing and signaling of environmental stress in plants also play important roles in the response to water deficit [4, 5]. The components that control and modulate stress-adaptive pathways mainly include transcription factors and protein kinases [2, 6]. Many transcription factors belonging to different transcription factor families, such as bZIP, AP2/ERF, MYB, NAC, WRKY, and zinc finger, are important regulators of the plant response to abiotic stress, and their activity can improve stress tolerance in transgenic plants [2, 3]. Moreover, these protein kinases, including calmodulin-dependent protein kinases (CDPKs), mitogen-activated protein kinases (MAPKs), receptor protein kinases (RPKs), and ribosomal protein kinases, participate in signal transduction processes in abiotic stress signaling and function as hubs in abiotic stress signaling.

Thinopyrum elongatum (syn. Lophopyrum elongatum or Agropyron elongatum), a perennial species in the tribe Triticeae and the genus Elytrigia, shares an ancestor with common wheat [7–9]. This species is easily crossed with common wheat, which makes it a good source for genetic improvement of wheat. In the past several decades, a number of genes from Th. elongatum have been introduced into common wheat to improve yield and provide resistance to wheat streak mosaic virus, barley yellow dwarf virus [10], stripe rust [11, 12], leaf rust [13], Fusarium head blight, and so on [14]. In addition to harboring pathogen resistance genes, Th. elongatum has also been used to improve tolerance to abiotic stresses, such as drought, waterlogging, and salinity, through the introduction of its chromosomes into wheat [7, 15–17]. However, the genome of Th. elongatum has not yet been published, which seriously limits the identification, characterization, and development of valuable genes in this species.

Recently, the development of next-generation sequencing (NGS) technologies and associated bioinformatics tools has provided a new method for transcriptomic research, that is, RNA-Seq [18]. RNA-Seq provides a precise way to measure transcript levels while simultaneously providing sequence information. RNA-Seq is highly efficient, reliable, and cost-effective, which makes it widely used to characterize the transcriptomes of plants, particularly nonmodel plants without published reference genomes [19–24].

In this study, we performed large-scale transcriptome sequencing of Th. elongatum under water deficit stress using ion torrent sequencing technology. We then compared the global expression profiles of Th. elongatum shoot and root tissues under control and water deficit stress conditions and identified a number of differentially expressed genes in response to water deficit stress. Gene annotation analysis of these genes provided novel insights into the response of Th. elongatum to water deficit stress, which should greatly facilitate wheat improvement in the future.

2. Materials and Methods

2.1. Plant Material and Water Deficit Treatment

Thinopyrum elongatum (PI 531718, ) seeds were kindly provided by GRIN (http://www.ars-grin.gov/), ARS, US Department of Agriculture. The seeds were pregerminated on wet filter paper in the dark at 25°C as described by Placido et al. [7]. When the coleoptiles were approximately 1 cm long, uniform seedlings were selected and transplanted to plastic pots filled with a mixture of surface soil collected from a field and washed sand. At 16 weeks after transplanting, each pot was supplied with 50 mL of water and 50 mL of half-strength Hoagland solution twice weekly (irrigated four times per week). Then, 16 weeks later, the plants were randomly divided into two groups, including the control group (watered normally as described above) and the water deficit group (supplied with only 50 mL of half-strength Hoagland solution twice weekly). After eight weeks of treatment, plant materials were collected from both the control and water deficit groups. All seedlings were grown in a greenhouse from March 2013 to September 2013 in Harbin, China. The greenhouse temperature was between 26°C and 30°C, the humidity ranged from 40% to 80%, and the light period was 06:00 to 18:00, as supplied by metal halide lamp 1 kW bulbs (Philips Lighting). The root and shoot tissue samples were separated, cleaned quickly, frozen in liquid nitrogen, and stored at −80°C for RNA isolation.

2.2. Total RNA Extraction, RNA-Seq Library Construction, and Sequencing

Frozen plant samples were ground in liquid nitrogen and total RNA was extracted using One Step RNA Reagent (Bio Basic Inc., Canada) as per the manufacturer’s protocol and purified using an RNeasy Plant Mini Kit (Qiagen, Valencia, CA). The integrity of the RNA was assessed by formaldehyde agarose gel electrophoresis. Total RNA was quantified using a NanoDrop ND-1000 spectrophotometer (Thermo Fisher Scientific, Wilmington, DE, USA) and a Bioanalyzer 2100 (Agilent Technologies, CA). RNA integrity number (RIN) values were greater than 8.0 for all samples. Ribosomal RNA depletion was carried out using a RiboMinus RNA plant kit for RNA-Seq (Life Technologies, CA). The whole-transcriptome cDNA library was prepared using an Ion Total RNA-Seq kit v2 (Life Technologies Corporation, CA). Double-stranded cDNA was ligated to barcoded adapters and sequenced by BGI-Shenzhen Ltd. (Shenzhen, China) using an Ion PI Chip (ion torrent, Life Technologies, CA). Processing of raw data, removal of adapter sequences, base-calling, and quality value calculations were performed using Torrent Suite Software 4.0 (ion torrent, Life Technologies, CA). Quality reads were obtained by trimming the raw reads at a minimum PHRED score of .

2.3. RNA-Seq Data Processing, De Novo Assembly, and Annotation

RNA-Seq reads were first processed with FASTX toolkit to remove low-quality sequences with parameters “− 33– 20– 70”. The resulting high-quality cleaned reads were assembled de novo into contigs using Trinity with the parameters “min_kmer_cov 2” [25]. To remove the redundancy of Trinity-generated contigs, the reads were further assembled de novo using iAssembler with minimum percent identity () set to 97 [26]. Blast searches of the resulting unique transcripts were performed against combined databases harboring Arabidopsis, rice, maize, and Brachypodium distachyon protein sequences with a cutoff -value of [27]. Gene ontology (GO) terms were assigned to the assembled transcripts based on the GO terms annotated to their corresponding homologs in the combined database, and the GO annotation results were explored using WEGO [28]. Annotations from MapMan were also retrieved based on homology search results [29]. Plant transcription factors (TF) and protein kinases were identified and classified into different families (or groups) using the iTAK pipeline (http://bioinfo.bti.cornell.edu/tool/itak/) [30].

2.4. Gene Expression Quantification and Differential Expression Analysis

High-quality cleaned RNA-Seq reads were aligned to the assembled Th. elongatum transcripts using the Bowtie program, allowing one mismatch [31]. Following the alignments, raw counts for each transcript and in each sample were derived and normalized to reads per kilobase of exon model per million mapped reads (RPKM). Differentially expressed genes (fold changes ≥ 2 or fold changes ≤ 0.5 and adjusted value ≤ 0.001) between normal and water deficit stress conditions were identified with the edgeR package [32]. GO terms enriched in the set of differentially expressed genes affected by water deficit stress were identified using the topGO package [33].

3. Results

3.1. Sequencing and De Novo Assembly of the Thinopyrum elongatum Transcriptome

To obtain a global view of water deficit stress-induced changes in Th. elongatum at the transcriptome level, we performed whole genome transcriptome sequencing of shoots and roots collected from control and water deficit-stressed plants using the ion proton platform [34]. After removing low-quality, adaptor, and barcode sequences, a total of 39,273,796 reads were obtained. All raw and processed data were submitted to the NCBI database (Accession numbers: SRX729803 and SRX729805-07). De novo assembly of these high-quality cleaned reads generated 169,990 unique transcripts with an average length of 550.5 bp; the longest transcript was 10,851 bp long. The length distribution of the assembled Th. elongatum unique transcripts is shown in Figure 1.

Figure 1: Length distribution of Thinopyrum elongatum unique transcripts.

3.2. Annotation of Thinopyrum elongatum Unique Transcript Sequences

The assembled Th. elongatum unique transcripts were annotated by Blast analysis against the combined databases, including Arabidopsis, rice, maize, and Brachypodium distachyon protein sequences, revealing a total of 81,061 (47.9%) unique transcripts with significant hits. Consistent with previous reports [6], the results show that the percentage of genes that could be annotated was positively correlated with the length of the genes, as shown in Figure 1. The short transcripts were annotated to fewer targets, and the longer transcripts generated more hits. Among these unique transcripts, 59,704 (35.1%) were assigned to at least one GO term in three main categories, that is, biological process, molecular function, and cellular component. We further classified these unique transcripts into different functional categories, as shown in Figure 2. The result shows that metabolic process (GO:0008152) was the most abundant group in the biological process category, followed by biological regulation (GO:0065007), while response to stimulus (GO:0050896) and response to stress (GO:0006950) were also common, which is consistent with the transcriptome data collected from Th. elongatum plants under water deficit stress. In the molecular function category, the most abundant groups included binding (GO:0005488), catalytic activity (GO:0003824), oxidoreductase activity (GO:0016491), transferase activity (GO:0016740), transporter activity (GO:0005215), and transcription regulator activity (GO:0030528). There were also transcripts classified into specific groups, such as antioxidant activity (GO:0016209), indicating that antioxidants play an important role in trapping free radicals to protect Th. elongatum from water deficit damage.

Figure 2: GO annotation results of Thinopyrum elongatum unique transcripts.

To investigate the transcriptional regulation process in detail, we used the iTAK pipeline to mine for transcription factors and protein kinases. In total, we identified 2,988 transcription factors classified into 77 different families and 3,154 protein kinases classified into 85 different families from among the Th. elongatum transcripts, shown as Figure 3. These TFs belong to many families that play important roles in the plant response to abiotic stress, such as C2H2, C3H, WRKY, MYB, SNF2, bZIP, bHLH, NAC, AUX/IAA, AP2-EREBP, CCAAT, and MADS. The protein kinases were classified into legume lectin domain kinase, leucine rich repeat kinase, DUF26 kinase, S domain kinase, GmPK6/AtMRK1 family, CDPK, SnRK, MAPK family, and so on. These protein kinases broadly participate in the regulation of gene expression, while protein kinases involved in the plant response to abiotic stress, especially the CDPK and MAPK families [35, 36], were also highly abundant in our transcriptome dataset.

(a)

(b)

Figure 3: Number of unique transcripts classified as transcription factor and protein kinase gene transcripts in Thinopyrum elongatum. (a) Transcription factor gene transcripts. (b) Protein kinase gene transcripts.

To estimate possible differences in transcript sequences between Th. elongatum and wheat, a BLASTN search was performed against the wheat transcript sequences from IWGSC. The results show that 46.63% (79,265/169,990) of unique transcripts from Th. elongatum had significant matches with wheat transcripts, most with high identity percentages, as shown in Figure 4. The remaining transcripts from Th. elongatum (53.37%) without significant matches in wheat represent Th. elongatum-specific genes, which could be beneficial for wheat improvement.

Figure 4: Sequence identity distribution of unique transcripts in Thinopyrum elongatum compared to wheat.

3.3. Differentially Expressed Genes under Water Deficit Stress

Using the edgeR Bioconductor package, we identified 1,300 and 3,604 differentially expressed transcripts from Th. elongatum shoot and root tissue, respectively, while 122 transcripts were differentially expressed in both tissues, shown as Figure 5. Among these transcripts, 2,690 were induced by water deficit stress in root tissue, while 914 were repressed in roots. There were almost three times as many upregulated transcripts as downregulated transcripts in roots. However, in shoots, 700 transcripts were induced while 600 were repressed. GO terms were assigned to all 4,782 differentially expressed transcripts, and enrichment analysis for GO annotation was performed using the topGO package; the results are shown in Table S1 in Supplementary Material available online at http://dx.doi.org/10.1155/2015/265791. As expected, GO terms in the biological process category were highly enriched, including GO:0050896 (response to stimulus), GO:0009628 (response to abiotic stimulus), GO:0006950 (response to stress), GO:0009651 (response to salt stress), and GO:0006970 (response to osmotic stress), which is consistent with previous reports in other plants. Meanwhile, GO terms in the molecular function category, such as GO:0005507 (copper ion binding), GO:0016491 (oxidoreductase), GO:0016209 (antioxidant), GO:0004784 (superoxide dismutase), GO:0022857 (transmembrane transporter), and GO:0005215 (transporter), were also highly enriched under water deficit stress in Th. elongatum, indicating that the antioxidant and transport systems play important role in protecting plants from damage due to environmental stress.

Figure 5: Venn diagram of the distribution of differentially expressed transcripts in shoot and root tissues.

We further annotated the functions of differentially expressed transcripts using MapMan. The results show that shoots and roots have different ways of responding to water deficit stress, as shown in Figure 6. In shoots, differentially expressed transcripts were more enriched in the category photorespiration. However, compared to shoot tissue, some transcripts were more abundant in roots, including those in the categories across cell wall, lipids, antioxidant (including ascorbate, glutathione, and OPP), and sucrose metabolism. Meanwhile, transcripts in the categories abiotic stress (20.2) and ascorbate and glutathione (21.2) were highly enriched, which is consistent with the GO annotation results. We analyzed the expression of transcripts in these two categories based on transcriptome data, revealing that most of these transcripts were induced under water deficit stress, as shown in Figure 7.

(a)

(b)

Figure 6: Overview of metabolic responses to water deficit stress. (a) Distribution of water deficit-responsive transcripts in shoots and (b) distribution of water deficit-responsive transcripts in roots.

Figure 7: Heatmap showing the expression profiles of transcripts involved in the response to abiotic stress and ROS metabolism.

4. Discussion

In plants, roots are often able to continue growing under water deficit stress in order to seek deeper water resources, while shoot elongation is completely inhibited due to a decline in photosynthesis. In this study, RNA-Seq technology was utilized to compare the shoot and root transcriptomes of Th. elongatum under water deficit stress to those grown under control conditions. We identified a total of 4,782 differentially expressed transcripts. Among these water deficit-responsive genes, 3,604 were detected in roots, while 1,300 were detected in shoots and only a few (122) were expressed in both roots and shoots, which is consistent with previous reports in maize and cotton [37, 38]. Placido et al. [7] examined the role of the brassinosteroid gene regulatory network in root development, finding that the drought tolerance of common wheat was improved through introgressing an alien chromosome segment from Th. elongatum. In the current study, we identified three differentially expressed transcripts involved in brassinosteroid metabolism and the brassinosteroid-mediated signaling pathway in Th. elongatum (ThUN007188, ThUN021506, and ThUN006380), which is consistent with the transcripts detected in wheat by Placido et al. [7]. We also found that GO categories GO:0044036 (cell wall macromolecule metabolic process, nine of 50 differentially expressed transcripts) and GO:0048364 (root development, six of 45 differentially expressed transcripts) were enriched in root tissues, which is consistent with the fact that brassinosteroids promote cell wall loosening, root elongation, and root development to mitigate the effect of water deficit stress on plant growth.

Our transcriptome data suggest that TFs play important roles in the response of Th. elongatum to water deficit stress, since we identified 47 differentially expressed TF genes. In plants, the bHLH TF directly regulates the GA and JA signaling pathways to induce the initiation of trichomes [39], which act as barriers to protect plants from water loss. Indeed, in this study, we identified 12 differentially expressed bHLH TF genes, including 10 expressed in shoots and three expressed in roots (with one expressed in both tissues). Moreover, nine transcripts involved in JA and GA metabolism were also differentially expressed in shoots compared to two in roots. These results suggest that the modulation of the GA and JA signaling pathways by bHLH in shoots helps protect Th. elongatum from water deficit stress. In addition to bHLH, other TFs reported to be involved in plant responses to abiotic stress, such as AP2/ERBP (9) [40, 41], MYB (8) [42–44] (6), and WRKY [45] (6), were also differentially expressed under water deficit stress, confirming that they play crucial roles in regulating transcription processes under water deficit stress.

ROS are produced in plant tissues due to the partial reduction of oxygen, for example, in the photosynthetic and respiratory electron chains, and their levels increase dramatically under environmental stress [46]. However, as ROS can cause cellular damage, they are scavenged (and generated) by oxidoreductases. Therefore, the homeostasis of ROS is determined by interactions between the ROS-producing and ROS-scavenging pathways in plants. ROS play important roles as signaling molecules that modulate many pathways [47], for example, MAPK cascades [48], and influence the activity of TFs under water deficit stress. Based on the transcriptome data, we identified 28 differentially expressed transcripts that are related to ROS metabolism, most of which were induced by water deficit stress, as shown in Figure 7. These transcripts were mainly classified into two highly represented GO categories, “oxidoreductase activity” and “antioxidant activity,” which were identified by topGO analysis. These genes encode enzymes including oxidoreductase, ascorbate, tocopherol, and glutathione, all of which scavenge harmful ROS to protect Th. elongatum from cellular damage [20, 49]. However, the detailed mechanisms of ROS metabolism are unknown. Additional studies are needed to fully elucidate the complex interactions of ROS metabolism under water deficit stress.

5. Conclusions

In the present study, thousands of transcripts were identified from Th. elongatum shoots and roots that were differentially expressed in response to water deficit stress, with three times as many transcripts in roots as in shoots. These transcripts are mainly involved in the response to stress, signal transduction, transcriptional regulation, ROS metabolism, and so on, especially the regulation of root development under water deficit. All of these genes help Th. elongatum adapt to water stress, and their introduction into cultivated wheat may improve the water deficit tolerance of this crop in the future.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This project received funding from the grant from Specialized Research Fund for the Doctoral Program of Higher Education of China (no. 20122329120001), supported by Heilongjiang Province Postdoctoral Science Foundation (no. LBH-Z14126), the National Basic Research Program of China (no. 2011CB111500), and the International S&T Cooperation Program of China (no. 2009DFA32470).

References

L. H. Comas, S. R. Becker, V. M. V. Cruz, P. F. Byrne, and D. A. Dierig, “Root traits contributing to plant productivity under drought,” Frontiers in Plant Science, vol. 4, article 442, 2013.

D. Golldack, I. Lüking, and O. Yang, “Plant tolerance to drought and salinity: stress regulating transcription factors and their functional significance in the cellular transcriptional network,” Plant Cell Reports, vol. 30, no. 8, pp. 1383–1391, 2011.

W. Wang, B. Vinocur, and A. Altman, “Plant responses to drought, salinity and extreme temperatures: towards genetic engineering for stress tolerance,” Planta, vol. 218, no. 1, pp. 1–14, 2003.

D. M. DeLaat, C. A. Colombo, A. F. Chiorato, and S. A. M. Carbonell, “Induction of ferritin synthesis by water deficit and iron excess in common bean (Phaseolus vulgaris L.),” Molecular Biology Reports, vol. 41, no. 3, pp. 1427–1435, 2014.

D. M. Neves, M. A. C. Filho, B. S. Bellete et al., “Comparative study of putative 9-cis-epoxycarotenoid dioxygenase and abscisic acid accumulation in the responses of Sunki mandarin and Rangpur lime to water deficit,” Molecular Biology Reports, vol. 40, no. 9, pp. 5339–5349, 2013.

Y. Xu, S. Gao, Y. Yang et al., “Transcriptome sequencing and whole genome expression profiling of chrysanthemum under dehydration stress,” BMC Genomics, vol. 14, no. 1, article 662, 2013.

D. F. Placido, M. T. Campbell, J. J. Folsom et al., “Introgression of novel traits from a wild wheat relative improves drought adaptation in wheat,” Plant Physiology, vol. 161, no. 4, pp. 1806–1819, 2013.

P. P. Jauhar, “Synthesis and cytological characterization of trigeneric hybrids involving durum wheat, Thinopyrum bessarabicum, and Lophopyrum elongatum,” Theoretical and Applied Genetics, vol. 84, no. 5-6, pp. 511–519, 1992.

P. P. Jauhar, “Multidisciplinary approach to genome analysis in the diploid species, Thinopyrum bessarabicum and Th. elongatum (Lophopyrum elongatum), of the Triticeae,” Theoretical and Applied Genetics, vol. 80, no. 4, pp. 523–536, 1990.

H. Sharma, H. Ohm, L. Goulart, R. Lister, R. Appels, and O. Benlhabib, “Introgression and characterization of barley yellow dwarf virus resistance from Thinopyrum intermedium into wheat,” Genome, vol. 38, no. 2, pp. 406–413, 1995.

Q. Huang, X. Li, W. Q. Chen et al., “Genetic mapping of a putative Thinopyrum intermedium-derived stripe rust resistance gene on wheat chromosome 1B,” Theoretical and Applied Genetics, vol. 127, no. 4, pp. 843–853, 2014.

L. J. Hu, G. R. Li, Z. X. Zeng, Z. J. Chang, C. Liu, and Z. J. Yang, “Molecular characterization of a wheat-Thinopyrum ponticum partial amphiploid and its derived substitution line for resistance to stripe rust,” Journal of Applied Genetics, vol. 52, no. 3, pp. 279–285, 2011.

G. M. Schachermayr, M. M. Messmer, C. Feuillet, H. Winzeler, M. Winzeler, and B. Keller, “Identification of molecular markers linked to the Agropyron elongatum-derived leaf rust resistance gene Lr24 in wheat,” Theoretical and Applied Genetics, vol. 90, no. 7-8, pp. 982–990, 1995.

S. Chen, Z. Huang, Y. Dai et al., “The Development of 7E Chromosome-Specific Molecular Markers for Thinopyrum elongatum Based on SLAF-seq Technology,” PLoS ONE, vol. 8, no. 6, Article ID e65122, 2013.

R. P. Jacoby, A. H. Millar, and N. L. Taylor, “Investigating the role of respiration in plant salinity tolerance by analyzing mitochondrial proteomes from wheat and a salinity-tolerant amphiploid (wheat lophopyrum elongatum),” Journal of Proteome Research, vol. 12, no. 11, pp. 4807–4829, 2013.

P. Monneveux, M. P. Reynolds, J. González Aguilar, and R. P. Singh, “Effects of the 7DL.7Ag translocation from Lophopyrum elongatum on wheat yield and related morphophysiological traits under different environments,” Plant Breeding, vol. 122, no. 5, pp. 379–384, 2003.

M. Taeb, R. M. D. Koebner, and B. P. Forster, “Genetic variation for waterlogging tolerance in the Triticeae and the chromosomal location of genes conferring waterlogging tolerance in Thinopyrum elongatum,” Genome, vol. 36, no. 5, pp. 825–830, 1993.

Z. Wang, M. Gerstein, and M. Snyder, “RNA-Seq: a revolutionary tool for transcriptomics,” Nature Reviews Genetics, vol. 10, no. 1, pp. 57–63, 2009.

S. R. Strickler, A. Bombarely, and L. A. Mueller, “Designing a transcriptome next-generation sequencing project for a nonmodel plant species,” The American Journal of Botany, vol. 99, no. 2, pp. 257–266, 2012.

G. Sathiyaraj, O. R. Lee, S. Parvin, A. Khorolragchaa, Y.-J. Kim, and D. C. Yang, “Transcript profiling of antioxidant genes during biotic and abiotic stresses in Panax ginseng C. A. Meyer,” Molecular Biology Reports, vol. 38, no. 4, pp. 2761–2769, 2011.

J. Wang, Q. Wang, Y. Yang et al., “De novo assembly and characterization of stress transcriptome and regulatory networks under temperature, salt and hormone stresses in Lilium lancifolium,” Molecular Biology Reports, 2014.

T. Bandyopadhyay, B. Gohain, R. Bharalee et al., “Molecular landscape of Helopeltis theivorainduced transcriptome and defense gene expression in tea,” Plant Molecular Biology Reporter, pp. 1–16, 2014.

H. Wu, D. Chen, J. Li et al., “De novo characterization of leaf transcriptomeusing 454 sequencing and development of EST-SSR markers in tea (Camellia sinensis),” Plant Molecular Biology Reporter, vol. 31, no. 3, pp. 524–538, 2013.

S. Yu, F. Zhang, Y. Yu, D. Zhang, X. Zhao, and W. Wang, “Transcriptome profiling of dehydration stress in the Chinese Cabbage (Brassica rapa L. ssp. pekinensis) by tag sequencin,” Plant Molecular Biology Reporter, vol. 30, no. 1, pp. 17–28, 2012.

B. J. Haas, A. Papanicolaou, M. Yassour et al., “De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis,” Nature Protocols, vol. 8, no. 8, pp. 1494–1512, 2013.

Y. Zheng, L. Zhao, J. Gao, and Z. Fei, “IAssembler: a package for de novo assembly of Roche-454/Sanger transcriptome sequences,” BMC Bioinformatics, vol. 12, article 453, 2011.

S. F. Altschul, T. L. Madden, A. A. Schäffer et al., “Gapped BLAST and PSI-BLAST: a new generation of protein database search programs,” Nucleic Acids Research, vol. 25, no. 17, pp. 3389–3402, 1997.

J. Ye, L. Fang, H. Zheng et al., “WEGO: a web tool for plotting GO annotations,” Nucleic Acids Research, vol. 34, pp. W293–W297, 2006.

O. Thimm, O. Bläsing, Y. Gibon et al., “MAPMAN: a user-driven tool to display genomics data sets onto diagrams of metabolic pathways and other biological processes,” Plant Journal, vol. 37, no. 6, pp. 914–939, 2004.

P. Pérez-Rodríguez, D. M. Riaño-Pachón, L. G. G. Corrêa, S. A. Rensing, B. Kersten, and B. Mueller-Roeber, “PlnTFDB: updated content and new features of the plant transcription factor database,” Nucleic Acids Research, vol. 38, no. 1, pp. D822–D827, 2009.

B. Langmead, “Aligning short sequencing reads with Bowtie,” in Current Protocols in Bioinformatics, chapter 11:unit 11.7, 2010.

M. D. Robinson, D. J. McCarthy, and G. K. Smyth, “edgeR: a bioconductor package for differential expression analysis of digital gene expression data,” Bioinformatics, vol. 26, no. 1, pp. 139–140, 2010.

A. Alexa, J. Rahnenführer, and T. Lengauer, “Improved scoring of functional groups from gene expression data by decorrelating GO graph structure,” Bioinformatics, vol. 22, no. 13, pp. 1600–1607, 2006.

M. A. Quail, M. Smith, P. Coupland et al., “A tale of three next generation sequencing platforms: comparison of Ion Torrent, Pacific Biosciences and Illumina MiSeq sequencers,” BMC Genomics, vol. 13, no. 1, article 341, 2012.

K. Moustafa, S. AbuQamar, M. Jarrar, A. J. Al-Rajab, and J. Trémouillaux-Guiller, “MAPK cascades and major abiotic stresses,” Plant Cell Reports, vol. 33, no. 8, pp. 1217–1225, 2014.

T. Asano, N. Hayashi, S. Kikuchi, and R. Ohsugi, “CDPK-mediated abiotic stress signaling,” Plant Signaling & Behavior, vol. 7, no. 7, pp. 817–821, 2012.

N. Opitz, A. Paschold, C. Marcon et al., “Transcriptomic complexity in young maize primary roots in response to low water potentials,” BMC Genomics, vol. 15, no. 1, article 741, 2014.

W. Park, B. E. Scheffler, P. J. Bauer, and B. T. Campbell, “Genome-wide identification of differentially expressed genes under water deficit stress in upland cotton (Gossypium hirsutum L.),” BMC Plant Biology, vol. 12, article 90, 2012.

T. Qi, H. Huang, D. Wu et al., “Arabidopsis DELLA and JAZ proteins bind the WD-repeat/bHLH/MYB complex to modulate gibberellin and jasmonate signaling synergy,” Plant Cell, vol. 26, no. 3, pp. 1118–1133, 2014.

J. Mizoi, K. Shinozaki, and K. Yamaguchi-Shinozaki, “AP2/ERF family transcription factors in plant abiotic stress responses,” Biochimica et Biophysica Acta—Gene Regulatory Mechanisms, vol. 1819, no. 2, pp. 86–96, 2012.

D. Kizis, V. Lumbreras, and M. Pagès, “Role of AP2/EREBP transcription factors in gene regulation during abiotic stress,” FEBS Letters, vol. 498, no. 2-3, pp. 187–189, 2001.

L. Zhang, G. Zhao, J. Jia, X. Liu, and X. Kong, “Molecular characterization of 60 isolated wheat MYB genes and analysis of their expression during abiotic stress,” Journal of Experimental Botany, vol. 63, no. 1, pp. 203–214, 2012.

L. Zhang, C. Xia, G. Zhao et al., “A novel wheat bZIP transcription factor, TabZIP60, confers multiple abiotic stress tolerances in transgenic Arabidopsis,” Physiologia Plantarum, 2014.

M. A. Hossain, J. I. Cho, M. Han et al., “The ABRE-binding bZIP transcription factor OsABF2 is a positive regulator of abiotic stress and ABA signaling in rice,” Journal of Plant Physiology, vol. 167, no. 17, pp. 1512–1520, 2010.

L. Chen, Y. Song, S. Li, L. Zhang, C. Zou, and D. Yu, “The role of WRKY transcription factors in plant abiotic stresses,” Biochimica et Biophysica Acta—Gene Regulatory Mechanisms, vol. 1819, no. 2, pp. 120–128, 2012.

G. Miller, N. Suzuki, S. Ciftci-Yilmaz, and R. Mittler, “Reactive oxygen species homeostasis and signalling during drought and salinity stresses,” Plant, Cell and Environment, vol. 33, no. 4, pp. 453–467, 2010.

S. Gilroy, N. Suzuki, G. Miller et al., “A tidal wave of signals: calcium and ROS at the forefront of rapid systemic signaling,” Trends in Plant Science, vol. 19, no. 10, pp. 623–630, 2014.

G. Xie, K. Sasaki, R. Imai, and D. Xie, “A redox-sensitive cysteine residue regulates the kinase activities of OsMPK3 and OsMPK6 in vitro,” Plant Science, vol. 227, pp. 69–75, 2014.

S. S. Gill and N. Tuteja, “Reactive oxygen species and antioxidant machinery in abiotic stress tolerance in crop plants,” Plant Physiology and Biochemistry, vol. 48, no. 12, pp. 909–930, 2010.