- Research article
- Open Access
An RNA sequencing transcriptome analysis of the high-temperature stressed tall fescue reveals novel insights into plant thermotolerance
© Hu et al.; licensee BioMed Central. 2014
- Received: 24 June 2014
- Accepted: 12 December 2014
- Published: 19 December 2014
Tall fescue (Festuca arundinacea Schreb.) is major cool-season forage and turf grass species worldwide, but high-temperature is a major environmental stress that dramatically threaten forage production and turf management of tall fescue. However, very little is known about the whole-genome molecular mechanisms contributing to thermotolerance. The objectives of this study were to analyzed genome-wide gene expression profiles in the leaves of two tall fescue genotypes, heat tolerant ‘PI578718’ and heat sensitive ‘PI234881’ using high-throughput RNA sequencing.
A total of 262 million high-quality paired-end reads were generated and assembled into 31,803 unigenes with an average length of 1,840 bp. Of these, 12,974 unigenes showed different expression patterns in response to heat stress and were categorized into 49 Gene Ontology functional subcategories. In addition, the variance of enrichment degree in each functional subcategory between PI578718 and PI234881 increased with increasing treatment time. Cell division and cell cycle genes showed a massive increase in transcript abundance in heat-stressed plants and more activated genes were detected in PI 578718 by Kyoto Encyclopedia of Genes and Genomes pathways analysis. Low molecular weight heat shock protein (LMW-HSP, HSP20) showed activated in two stressed genotypes and high molecular weight HSP (HMW-HSP, HSP90) just in PI578718. Assimilation such as photosynthesis, carbon fixation, CH4, N, S metabolism decreased along with increased dissimilation such as oxidative phosphorylation.
The assembled transcriptome of tall fescue could serve as a global description of expressed genes and provide more molecular resources for future functional characterization analysis of genomics in cool-season turfgrass in response to high-temperature. Increased cell division, LMW/HMW-HSP, dissimilation and antioxidant transcript amounts in tall fescue were correlated with successful resistance to high temperature stress.
- High-temperature stress
- RNA sequencing
- Cell cycle
- Energy metabolism
- Tall fescue
High-temperature stress has devastating effects on plant growth and productivity, and is a growing concern in the context of global climate change [1, 2]. Temperatures 5°C above optimal growing conditions will induce heat shock or heat stress to plants . Every year, there is substantial crop yield loss due to the elevated temperature and drought stress . The development of genetically heat-resistant plants has been considered as a promising approach for alleviating the threats from high-temperature stress .
Heat tolerance of plants is a complex multigenic process with different regulatory mechanisms sets throughout the gene network during plant development [6, 7]. Larkindale et al. [8, 9] reported that each of 45 Arabidopsis mutants had different responses to high temperature, which were involved in signaling pathways [abscisic acid (ABA), ethylene (ET), salicylic acid (SA), calcium, oxidative burst] and reactive oxygen metabolism (ascorbic acid or glutathione production, catalase). Additionally, ABA, SA or active oxygen molecule may serve as one signal to regulate gene chain in heat tolerant Arabidopsis. Although many researches on heat-stress tolerance have advanced considerably in recent years, genome-wide comparisons in plants have not yet been made using the next-generation high-throughput sequencing technologies. Determining the multigenic net regulatory mechanisms of plant response to high temperature, particularly at the transcriptomic level, will be helpful in the development of heat-tolerant species.
Tall fescue (Festuca arundinacea Schreb.) is a forage and cool-season turfgrass species grown widely in the temperate regions of the world such as in United States, China, Japan, Australia and many countries in Europe because of its agronomic importance [10–12]. However, high temperature (30–41°C) is one of the limiting factors affecting forage production and turf management of tall fescue in the south such as the Great Plains of the U.S. and south of the Yangtze Rive in China .
The genome size of tall fescue is approximately 6 × 103 Mbp and about 14 times larger than that of rice (Oryza sativa L.) . To date, complete genome sequences of tall fescue are not achieved. Although 59, 578 expressed sequence tags (ESTs) had been deposited in the NCBI GenBank database as of June 2014, molecular resouces of tall fescue are still limited because of the lack of genomic and transcriptomic information. An accelerated effort to acquire transcriptomes of tall fescue in response to high-temperature stress will be helpful in the development of heat-tolerant tall fescue cultivars. In the current study, two tall fescue genotypes, heat tolerant PI578718 and heat sensitive PI234881 identified from 120 tall fescue genotypes in our previous study (Additional file 1) , were used to investigate responses to high temperature at a global transcriptional level using an high-throughput RNA sequencing (RNA-Seq) approach. The objectives of this study were (1) to provide an initial investigation of genome-wide transcriptomic analysis; (2) to reveal novel insights into the regulation pattern of certain candidate genes in the critical metabolic pathways response to heat stress in cool-season turfgrass species.
Analyses of RNA-Seq data
Overview of the obtained RNA-Seq data from the Trinity de novo assembly program in tall fescue leaves induced by heat stress
Assembly statistic (Contig)
Total Length (bp)
Max Length (bp)
Ave Length (bp)
>N50 Reads No.
Assembly statistic (Transcript)
Total Length (bp)
Max Length (bp)
Ave Length (bp)
>N50 Reads No.
Total Length (bp)
Max Length (bp)
Ave Length (bp)
>N50 Reads No.
Using the Trinity software, a final total of 31,803 unigenes (N50 of 2584 bp) from 967,565 transcripts were generated, whose average length of 1,840 bp was longer than those reported in other plants, such as litchi (Litchi chinensis Sonn.; 601 bp) , bracken fern (Pteridium aquilinum; 547.2 bp) , tea (Camellia sinensis; 355 bp) . 15,689 (49.33%) unigenes of these 31,803 unigenes were matched with the homologs in databases with known function for tall fescue. A transcriptome atlas of two wheat genotypes response to heat treatments using microarray analyses only detected the expression of 22,720 probe sets . Other suppression subtractive hybridization (SSH)-based studies involved in high temperature response just detected the expression of 1,090 transcripts in fescue (Festuca sp.) seeds of PI 297901 and PI283316 . The unigenes obtained from a recent whole-genome-wide transcriptional analysis in two tall fescue genotypes under stress could greatly enrich the gene expression library response to high temperature in the grass family, especially in turfgrass and forage grasses. This Transcriptome Shotgun Assembly project has been deposited at DDBJ/EMBL/GenBank under the accession GBYN00000000. The version described in this paper is the first version, GBYN01000000.
GO enrichment of the tall fescue heat transcriptome
The number of functional subcategories and variance degree in these three categories increased along with increased treatment time (Additional files 2, 3, 4). PI 578718, one of the most heat-hardy wild tall fescue genotype, may have a larger number of genes enriching in GO categories. However, RNA-Seq results did not support this hypothesis. On the contrary, a greater number of genes (1,568 vs. 1,363) were seen in the heat sensitive PI234881. In addition, the gene variance degree between the two genotypes in the same GO subcategories reached the maximum at the longest treatment time. At 36 HAT, genes involved in 18 subcategories from three main GO categories in the heat-sensitive genotype had more than two times enrichment than in the heat- tolerant genotype. Furthermore, genes, just encoding abscission, death, enzyme regulator activity, pollination, cytoskeleton, Golgi apparatus, nucleolus, ribosome, structural molecule activity and thylakoid proteins, were more enriched in the heat-tolerant genotype than in the heat-sensitive genotype at 36 HAT (Figure 1b). Similar results were reported in grapevine shoot apices in response to cold stress .
COG functional annotation and classification of the tall fescue transcriptome
The unigenes variance degree between the two genotypes in the same COG categories also exhibited the maximum at the longest treatment time. At 36 HAT, unigenes involved in 18 COG categories in the heat-sensitive genotype had more than two times enrichment than in the heat- tolerant genotype (Figure 2c). At 0 HAT, unigenes just encoding ‘coenzyme transport and metabolism’ in the heat-sensitive tall fescue had more than two times enrichment than in the heat-tolerant tall fescue (Figure 2a). At 12 HAT, unigenes encoding ‘translation, ribosomal structure and biogenesis’ and ‘cytoskeleton’ proteins in the heat-sensitive tall fescue had more than two times enrichment than in the heat-tolerant tall fescue (Figure 2b). Unigenes involved in ‘chromatin structure and dynamics’ and ‘defense mechanisms’ in the heat-tolerant tall fescue had more two times enrichment than in the heat-sensitive tall fescue at 0 HAT and 12 HAT, respectively (Figure 2a, b). Under high-temperature stress, higher enrichment of the tall fescue unigenes in the most of COG categories for heat-sensitive genotypes compared with heat-tolerant ones suggests that thermotolerance was related with the changed expression of genes involved in transcriptional regulation as well as metabolic pathways .
Differential expression of transcripts are affected by high-temperature stress
Cell cycle and division are affected by high temperature stress
In response to heat stress in tall fescue, cell division control protein 6 gene (CDC6, comp_71_104860_c0_seq2) initiating DNA replication showed increased expression at 12 and 36 HAT in PI 578718. Under exogenous environment stress, S-phase cells are susceptible to chromosomal damage because of the potential replication errors . Minichromosome maintenance (MCM) and origin recognition complex (ORC) family proteins are essential components in DNA replication licensing and play a key role in preventing centrosome reduplication, particularly for MCM complex in S-phase checkpoint activation and in subsequent lesion repair [25, 26]. At 12 HAT, ORC1 (comp_159_122824_c0_seq5) and ORC5 (comp_159_115294_c0_seq3) were activated in PI578718 and PI234881, respectively, but were unaffected in grass leaves at 36 HAT. The transcript abundance of both MCM2 (comp_159_114387_c1_seq5) and MCM6 (comp_159_125142_c0_seq1) increased in the heat-stressed PI578718, but were unaffected in heat-stressed PI234881. The result suggests that the heat-tolerant PI578718 had greater ability to self-repair to protect S-phase genome stability than the heat-sensitive PI 234881.
Plant cyclins are typically classified as A, B, C, or D-types . Cyclin A protein is present in the S-phase to M phase and plays a major role in the control of DNA replication, and cyclin D regulates the G1-S transition . Two cyclins [cyclin A (CCNA) and cyclin D3 (CYCD3)] were detected for the first time in tall fescue by genome wide RNA-Seq analysis. CCNA (comp_159_115639_c1_seq1) was identified activated in heat-stressed PI578718, unaffected in heat-stressed PI234881. However, CYCD3 showed decreased transcript abundance in 12 h-treated PI578718. The 14-3-3 proteins are a family of phosphoserine/phosphothreonine binding molecules that promote cell survival and control the function of a wide array of cellular proteins . 14-3-3 epsilon (YWHAE, comp_71_104651_c2_seq1) was detected with increased expression at 36 HAT in PI 578718, but not in heat-stressed PI 234881.
One sister-chromatid cohesion gene (SCC1) plays a vital role in keeping sister chromatid cohesion at spindle checkpoint until the activation of separase by the anaphase-promoting complex/cyclosome (APC/C) . Cdc20/Fizzy and Cdh1/Fizzy-related complex are critical for the activity of multisubunit E3 ubiquitin ligase APC/C and then triggers the degradation of multiple substrates during mitosis . In this study, SCC1 (comp_159_122824_c0_seq5) and CDC20 (comp_71_114623_c0_seq3) showed increased transcript abundance at 36 HAT in PI 578718. CDC20 showed increased expression at 12 HAT in PI 234881. APC10 (comp_159_108228_c1_seq2) and E-cadherin gene (CDH1) (comp_159_134911_c0_seq34) had more transcript abundance in PI234881 than in PI578718 at 0 and 12 HAT, respectively. However, there was no significant change between PI234881 and PI578718 at 36 HAT. The results suggest that SCC1and CDC20 play more important role in protecting M-phase stability in tall fescue in response to heat stress.
Energy metabolism is affected by high temperature stress
In the tall fescue leaf, of the 473 annotated unigenes, 164 unigenes in the KEGG category “energy metabolism” were significantly affected by high temperature, mainly involved in six kinds of assimilation pathway: “photosynthesis (KO 00195)”, “photosynthesis-antenna proteins (KO 00196)”, “carbon fixation in photosynthetic organisms (KO 00710)”, “methane metabolism (KO 00680)”, “nitrogen metabolism (KO 00910)” and “sulfur metabolism (KO 00920)”, and one dissimilation “oxidative phosphorylation (KO 00190)” (Additional file 6). Assimilation and dissimilation are the two important aspects of the energy metabolism, and the balance of energy metabolism process between assimilation and dissimilation is fundamental to the maintenance of cell life, the source of its life and development [31, 32]. In “energy metabolism”, most of the genes (75.0%, 66) involved in assimilation, such as “photosynthesis”, “photosynthesis-antenna proteins” and “carbon fixation in photosynthetic organisms” showed decreased transcript abundance in two tall fescue genotypes under heat stress. In the methane, nitrogen and sulfur metabolic pathways, there were also more genes (70.5%, 31) inhibited in the two tall fescue genotypes. However, most of the genes involved in dissimilation, such as “oxidative phosphorylation” showed increased transcript abundance in the heat-tolerant PI 578718 (78.57%, 22) and only inhibited in heat-sensitive PI234881 (66.7%, 6). The result finding indicates that PI578718 obtained more energy than PI234881 in response to high temperature stress. The expression data suggest that greater dissimilation metabolism in PI578718 relative to PI512315 may result in a more stable cell division in PI578718 than in PI234881 in response to high temperature. In the three pathways methane, nitrogen and sulfur metabolism, there were more genes inhibited in PI578718 than in PI234881. One gene encoding cytochrome c oxidase subunit 2 (COX2) showed increased abundance in PI 578718, but decreased in PI 234881 under heat stress. The transcript abundance of a ferredoxin gene (petF) belonging to this pathway was dramatically increased at all treatment times in PI 578718 but decreased at 36 HAT in PI234881.
HSP and antioxidant genes are affected by high temperature stress
Plants respond to high temperature by synthesizing a set of the heat shock proteins (HSPs) to acquire thermotolerance . According to their predicted molecular weight, HSPs could be grouped into three families: high molecular weight (HMW-HSP, 80–114 kDa), HSP70 (69–71 kDa), low molecular weight (LMW-HSP, 15–30 kDa) . Many researches have proved that the increased transcript abundance of the three kinds of HSPs and their molecular chaperone was directly involved in the activation of protection mechanisms in response to heat stress [34–37]. However, in contrast with these results, very few HSP70 genes showed differential expression in the two tall fescue genotypes under heat stress (Additional file 7). Moreover, heat transfer products group (htpG) encoding the molecular chaperone HtpG protein and HSPBP1 encoding hsp70-interacting protein showed decreased transcript abundance at 12 HAT in PI578718. It is interesting to find that nine HSP20 genes showed increased transcript abundance in all treatment lines for the both genotypes, and two HSP90B genes increased abundance in heat-stressed PI 578718, which suggests that the LMW-HSP and HMW-HSP may play a more important role in enhancing thermotolerance for tall fescue. There is new evidence showing that plants have developed antioxidant enzymatic systems for scavenging antioxidant enzymes to avoid these oxidative injuries induced by abiotic stress [38, 39]. In this study, 14 antioxidant genes, in the KEGG category “peroxisome”, regulating antioxidant enzyme activity, were significantly affected in both heat-stressed genotypes and the tolerant genotype had more antioxidant genes activated than in the sensitive genotype (Additional file 8). Particularly, three SOD genes increased transcript abundance at 36 HAT in PI578718 but were unaffected in heat-stressed PI234881. The result suggests that the high expression level of antioxidant genes may be one key functional trait to identify the heat tolerance in tall fescue species.
Quantitative real-time-PCR validation of differentially expressed transcripts from RNA-Seq
In conclusion, this study provided the first large-scale transcriptome dataset in cool-season turfgrass in response to high-temperature stress. More than 262 million high-quality 100-bp paired-end reads paired-end reads were generated and assembled 31,803 unigenes in tall fescue. Furthermore, 12,974 unigenes expressed specifically in tall fescue under high-temperature stress, as well as their classification, functions, and metabolic pathways were shown for the first time. High expression levels of genes involved in LMW (HMW)-HSP, antioxidants, cell cycle and cell division and dissimilation (oxidative phosphorylation) in tall fescue respond to heat stress could play a vital role in achieving thermotolerance.
Plant materials and growth conditions
Healthy and single clonal plants of two tall fescue (Festuca arundinacea Schreb.) genotypes, heat-tolerant PI578718 and heat-sensitive PI234881 were collected from turfgrass field plots at Wuhan Botanical Garden, Chinese Academy of Sciences, Wu Han, Hubei, China. Tall fescue plants were grown in plastic pots (13 cm diameter, 11 cm deep) with a mix of sand and peat soil (1/1, v/v) in the greenhouse under natural sunlight, relative humidity of 87%, wind speed of 0.8 m · s−1 and temperatures of 22/18°C (day/night). After 3 months of establishment, plants were transferred into controlled-environment growth chambers (HP300GS-C; Ruihua Instrument, Wuhan, China), with a 14-h photoperiod, photosynthetically active radiation at 450 μmol m−2 s−1 at the canopy level, a day/night temperature of 22/18°C and 70% relative humidity. The plants were allowed to acclimate for one week before they were exposed to heat treatments. Throughout the experiment, the plants were fertilized three times a week with half-strength Hoagland’s solution  until dripping occurred.
Treatments and experimental design
The plant-pot system was weighed at 48-h intervals to determine transpiration rate (Tr) before heat treatment initiation based on the method described by Hu et al. . All plants from each tall fescue genotype were divided into two groups with similar Tr (Group I, II). Group I of tall fescues was kept at 22/18°C day/night in the controlled-environment growth chambers. Group II was subjected to heat treatments at 40°C in the controlled-environment growth chambers (except for temperature, the other conditions remained the same as the control). At the beginning of heat treatment, all plots were hand irrigated until dripping occurred; upon reaching this point, the irrigation was withheld during the experiment. Heat treatment began at 7: 00 am and it took about 45 min for the temperature to increase from 22°C to 40°C. The plants were exposed to prolonged 40°C heat stress for 36 h at 450 μmol m−2 s−1 and 70% relative humidity. The leaf samples for Illumina deep sequencing analysis were collected at 12 and 36 h after plants were exposed to 40°C, respectively. At each sampling time point, the leaf tissues from the three pots (three replicates, one pot per replicate) of each genotype pooled together as one biological replicate and froze immediately with liquid nitrogen, and then stored at −80°C for RNA-Seq analysis. A total of 12 control samples including the leaves collected at 0 h and 22°C in Group II were pooled together as the sample at 0 h for each genotype. A total of six samples in total were used for Illumina Genome Analyzer deep sequencing and 6 G data performed for each sample.
Total RNA, Poly(A) RNA Isolation, and Library Preparation
Total RNA from tall fescue leaf tissues was extracted separately using Trizol reagent (Invitrogen, Carlsbad, CA) and purified with the RNeasy Plant Mini Kit (Qiagen, Valencia, CA). The concentration and quality of the total RNA was determined by a NanoDrop 8000 spectrophotometer (NanoDrop, Wilmington, DE) and checked by running a gel electrophoresis in 1.5% (w/v) agarose gels. According to the manufacturer’s instructions (Illumina, San Diego, CA), the poly (A) mRNA was isolated from purified total RNA using biotin-Oligo (dT) magnetic beads and fragmented into small pieces using an RNA fragmentation kit. First-strand cDNA was generated from the cleaved RNA fragments using reverse transcriptase and random primers. Second-strand cDNA was synthesized using RNase H and DNA Polymerase I. Following adaptor ligation, the unsuitable fragments were removed with AMPureXP beads, 200-bp cDNA fragments were enriched by 18 cycles of PCR and checked with Pico green staining quantified with Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA). The products were sequenced on the Illumina HiSeq 2000 instrument in Shanghai Personal Biotechnology Co., Ltd.
Processing and Mapping of Illumina Reads
The RNA-Seq raw reads were processed to obtain high-quality reads by removing the adapter sequences and low-quality bases at the 3’ end, trimming low-quality bases (Q < 20) from the 5’ and 3’ ends of the remaining reads. Reads filtering out reads containing ‘N’ and greater than 25 bp were considered for analysis. Clean reads were assembled into contigs, transcripts and unigenes with Trinity software (http://trinityrnaseq.sf.net). RPKM was used to normalize the abundances of transcripts . More than a 2-fold change was used to identify the significance of different gene expression between different treatment lines.
Identification of differentially expressed unigenes
All unigenes were blastx searched against Brachypodium hereafter (Brachypodium distachyon), wheat (Triticum spp.), rice (Oryza sativa Linn.), maize (Zea mays L.), sorghum [Sorghum bicolor (L.) Moench] in KEGG and eggNOG protein databases (E-value < 10−5) and functionally annotated by Blast2GO Gene Ontology Functional Annotation Suit (E-value < 10−5) (http://www.blast2go.org/). Metabolic pathways were predicted by KEGG mapping. Putative transcription factors were identified by searching Arabidopsis Gene Regulatory Information Server (AGRIS) Database .
Hierarchical clustering analysis was carried out for all annotated transcripts from heat tolerant PI578718 and heat sensitive PI234881. The RPKM counts for each transcript were clustered with the software Cluster 3.0, and JAVA Treeview was used to view the cluster image. The results were visualized using JAVA Treeview .
Validation of RNA-Seq by qRT-PCR
First-strand cDNA was generated from 3 μg of the same total RNA used by RNA-Seq analysis. The expression level of the target genes was determined by means of quantitative Real-Time PCR (qPCR) using ABI StepOne Plus Real-Time PCR system (Applied Biosystems, Foster City, CA) and SYBR Green Real-Time PCR Master Mix (Toyobo, Osaka, Japan) in 20 μL reactions. Each reaction consisted of a 2 ng of total RNA, 0.5 μL of each primer and 10 μL master mix. The PCRs were performed in thermocycler conditions starting with 3 min at 95°C, 40 cycles of 10 s at 94°C, 20 s at 55°C, and 20 s at 72°C, followed by 5 min at 72°C. To verify the presence of a specific product, the melting curve analysis of amplification products was performed at the end of each PCR reaction. In addition, the size of each amplified DNA fragment was verified on a 1.5% (w/v) agarose–ethidium bromide gel at 100 V for 40 min in 1 × TE buffer (10 mM Tris, 1 mM EDTA). Primers (see Additional file 8) for qPCR were designed with Primer Premier software (Primer Premier v5.0; Premier Biosoft International, Palo Alto, Calif.). The comparative CT method of quantification was used to quantify the relative expression of specific genes . YT521-B gene was used as references . The experiments were repeated twice with three replicates.
To assess the correlation between different platforms, the Pearson correlations were calculated by SPSS 16.0 to compare the mRNA expression levels measured by RNA-Seq and qPCR. To assess treatment effects between different treatment lines, the variance analysis was performed using the Statistical Analysis System (SAS 9.0 for windows, SAS Institute Inc., Cary, NC). The mean separation was performed with Fisher’s least significant difference test at P < 0.05.
We thank Dr. Wang Lu (Wuhan Botanical Garden, Chinese Academy of Science) for helpful suggestions on bioinformatics analysis of RNA-Seq data. This work was supported by the National Natural Science Foundation of China (Grant #: 31071822; 31272194) and the National High Technology Research and Development Program of China (i.e., “863” plan, Grant #: 2011AA100209-2). We are also grateful to Shanghai Personal Biotechnology Co., Ltd. for RNA-Seq technical support.
- Xu C, Huang B: Root proteomic responses to heat stress in two Agrostis grass species contrasting in heat tolerance. J Exp Bot. 2008, 59: 4183-4194. 10.1093/jxb/ern258.PubMed CentralPubMedView ArticleGoogle Scholar
- Shah F, Huang J, Cui K, Nie L, Shah T, Chen C, Wang K: Impact of high-temperature stress on rice plant and its traits related to tolerance. J Agr Sci. 2011, 149: 545-556. 10.1017/S0021859611000360.View ArticleGoogle Scholar
- Zhang Y, Mian MR, Chekhovskiy K, So S, Kupfer D, Lai H, Roe BA: Differential gene expression in Festuca under heat stress conditions. J Exp Bot. 2005, 56: 897-907. 10.1093/jxb/eri082.PubMedView ArticleGoogle Scholar
- Qin D, Wu H, Peng H, Yao Y, Ni Z, Li Z, Zhou C, Sun Q: Heat stress-responsive transcriptome analysis in heat susceptible and tolerant wheat (Triticum aestivum L.) by using Wheat Genome Array. BMC Genomics. 2008, 9: 432-10.1186/1471-2164-9-432.PubMed CentralPubMedView ArticleGoogle Scholar
- Wahid A, Gelani S, Ashraf M, Foolad M: Heat tolerance in plants: an overview. Environ Exp Bot. 2007, 61: 199-223. 10.1016/j.envexpbot.2007.05.011.View ArticleGoogle Scholar
- Suzuki N, Rizhsky L, Liang H, Shuman J, Shulaev V, Mittler R: Enhanced tolerance to environmental stress in transgenic plants expressing the transcriptional coactivator multiprotein bridging factor 1c. Plant Physiol. 2005, 139: 1313-1322. 10.1104/pp.105.070110.PubMed CentralPubMedView ArticleGoogle Scholar
- Kotak S, Larkindale J, Lee U, von Koskull-Döring P, Vierling E, Scharf K-D: Complexity of the heat stress response in plants. Curr Opin Plant Biol. 2007, 10: 310-316. 10.1016/j.pbi.2007.04.011.PubMedView ArticleGoogle Scholar
- Larkindale J, Knight M-R: Protection against heat stress-induced oxidative damage in Arabidopsis involves calcium, abscisic acid, ethylene, and salicylic acid. Plant Physiol. 2002, 128: 682-695. 10.1104/pp.010320.PubMed CentralPubMedView ArticleGoogle Scholar
- Larkindale J, Hall JD, Knight MR, Vierling E: Heat stress phenotypes of Arabidopsis mutants implicate multiple signaling pathways in the acquisition of thermotolerance. Plant Physiol. 2005, 138: 882-897. 10.1104/pp.105.062257.PubMed CentralPubMedView ArticleGoogle Scholar
- Barnes RF: Importance and problems of tall fescue. Biotechnology in tall fescue improvement. Edited by: Kasperbauer MJ. 1990, 1-12. ISBN 0-8493-4891-9Google Scholar
- Wang Z, Hopkins A, Mian R: Forage and turf grass biotechnology. Crit Rev Plant Sci. 2001, 20: 573-619. 10.1016/S0735-2689(01)80005-8.View ArticleGoogle Scholar
- Wang Z-Y, Ge Y, Scott M, Spangenberg G: Viability and longevity of pollen from transgenic and nontransgenic tall fescue (Festuca arundinacea) (Poaceae) plants. Am J Bot. 2004, 91: 523-530. 10.3732/ajb.91.4.523. doi:10.3732/ajb.91.4.523PubMedView ArticleGoogle Scholar
- Mian MR, Zhang Y, Wang ZY, Zhang JY, Cheng X, Chen L, Chekhovskiy K, Dai X, Mao C, Cheung F: Analysis of tall fescue ESTs representing different abiotic stresses, tissue types and developmental stages. BMC Plant Biol. 2008, 8: 27-10.1186/1471-2229-8-27.PubMed CentralPubMedView ArticleGoogle Scholar
- Sun X, Hu L, Xie Y, Fu J: Evaluation of genotypic variation in heat tolerance of tall fescue by functional traits. Euphytica. 2014, 199: 247-260. 10.1007/s10681-014-1122-4.View ArticleGoogle Scholar
- Zhang Y, Zwonitzer J, Chekhovskiy K, May G-D, Mian M-R: A functional genomics approach for identification of heat tolerance genes in tall fescue. Mol Breed Forage Turf. 2004, 11: 87-96. 10.1007/1-4020-2591-2_8.View ArticleGoogle Scholar
- Li C, Wang Y, Huang X, Li J, Wang H, Li J: De novo assembly and characterization of fruit transcriptome in Litchi chinensis Sonn and analysis of differentially regulated genes in fruit in response to shading. BMC Genomics. 2013, 14: 552-10.1186/1471-2164-14-552.PubMed CentralPubMedView ArticleGoogle Scholar
- Der JP, Barker MS, Wickett NJ, Wolf PG: De novo characterization of the gametophyte transcriptome in bracken fern, Pteridium aquilinum. BMC Genomics. 2011, 12: 99-10.1186/1471-2164-12-99.PubMed CentralPubMedView ArticleGoogle Scholar
- Shi C, Yang H, Wei C, Yu O, Zhang Z, Jiang C, Sun J, Li Y, 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 CentralPubMedView ArticleGoogle Scholar
- Xin H, Zhu W, Wang L, Xiang Y, Fang L, Li J, Sun X, Wang N, Londo J-P, Li S: Genome wide transcriptional profile analysis of Vitis amurensis and Vitis vinifera in response to cold stress. Plos One. 2013, 8: e58740-10.1371/journal.pone.0058740.PubMed CentralPubMedView ArticleGoogle Scholar
- Han X, Zhao CZ, Hou L, Li AQ, Zhao SZ, Bi YP, An J, Zhao YX, Wan SB, Wang XJ: Transcriptome profiling of peanut gynophores revealed global reprogramming of gene expression during early pod development in darkness. BMC Genomics. 2013, 14: 517-10.1186/1471-2164-14-517.View ArticleGoogle Scholar
- Lai ZX, Lin YL: Analysis of the global transcriptome of longan (Dimocarpus longan Lour.) embryogenic callus using Illumina paired-end sequencing. BMC Genomics. 2013, 14: 561-10.1186/1471-2164-14-561.PubMed CentralPubMedView ArticleGoogle Scholar
- Lockwood BL, Sanders JG, Somero GN: Transcriptomic responses to heat stress in invasive and native blue mussels (genus Mytilus): molecular correlates of invasive success. J Exp Biol. 2010, 213: 3548-3558. 10.1242/jeb.046094.PubMedView ArticleGoogle Scholar
- Hurford R, Cobrinik D, Lee MH, Dyson N: pRB and p107/p130 are required for the regulated expression of different sets of E2F responsive genes. Gene Dev. 1997, 11: 1447-1463. 10.1101/gad.11.11.1447.PubMedView ArticleGoogle Scholar
- Zheng N, Schulman BA, Song L, Miller JJ, Jeffrey PD, Wang P, Chu C, Koepp D-M, Elledge SJ, Pagano M, Conaway RC, Conaway JW, Harper JW, Pavletich NP: Structure of the Cul1-Rbx1-Skp1-F boxSkp2 SCF ubiquitin ligase complex. Nature. 2002, 416: 703-709. 10.1038/416703a.PubMedView ArticleGoogle Scholar
- Bailis JM, Forsburg S-L: MCM proteins: DNA damage, mutagenesis and repair. Curr Opin Genet Dev. 2004, 14: 17-21. 10.1016/j.gde.2003.11.002.PubMedView ArticleGoogle Scholar
- Ferguson RL, Pascreau G, Maller J-L: The cyclin A centrosomal localization sequence recruits MCM5 and Orc1 to regulate centrosome reduplication. J Cell Sci. 2010, 123: 2743-2749. 10.1242/jcs.073098.PubMed CentralPubMedView ArticleGoogle Scholar
- De Veylder L, Beeckman T, Inze D: The ins and outs of the plant cell cycle. Nat Rev Mol Cell Bio. 2007, 8: 655-665. 10.1038/nrm2227.View ArticleGoogle Scholar
- Jang S, Liu X, Fu H, Rees H, Yepes M, Levey A, Ye K: Interaction of Akt-phosphorylated SRPK2 with 14-3-3 mediates cell cycle and cell death in neurons. J Biol Chem. 2009, 284: 24512-24525. 10.1074/jbc.M109.026237.PubMed CentralPubMedView ArticleGoogle Scholar
- Sonoda E, Matsusaka T, Morrison C, Vagnarelli P, Hoshi O, Ushiki T, Nojima K, Fukagawa T, Waizenegger IC, Peters JM: Scc1/Rad21/Mcd1 is required for sister chromatid cohesion and kinetochore function in vertebrate cells. Dev Cell. 2001, 1: 759-770. 10.1016/S1534-5807(01)00088-0.PubMedView ArticleGoogle Scholar
- Pesin JA, Orr-Weaver TL: Regulation of APC/C activators in mitosis and meiosis. Annu Rev Cell Dev Bi. 2008, 24: 475-499. 10.1146/annurev.cellbio.041408.115949.View ArticleGoogle Scholar
- Yinger J-M: Toward a theory of assimilation and dissimilation. Ethn Rac Stud. 1981, 4: 249-264. 10.1080/01419870.1981.9993338.View ArticleGoogle Scholar
- Huner N, Öquist G, Sarhan F: Energy balance and acclimation to light and cold. Trends Plant Sci. 1998, 3: 224-230. 10.1016/S1360-1385(98)01248-5.View ArticleGoogle Scholar
- Howarth C: Molecular responses of plants to an increased incidence of heat shock. Plant Cell Environ. 1991, 14: 831-841. 10.1111/j.1365-3040.1991.tb01446.x.View ArticleGoogle Scholar
- Wang D, Luthe DS: Heat sensitivity in a bentgrass variant. failure to accumulate a chloroplast heat shock protein isoform implicated in heat tolerance. Plant Physiol. 2003, 133: 319-327. 10.1104/pp.102.018309.PubMed CentralPubMedView ArticleGoogle Scholar
- Tian J, Belanger FC, Huang B: Identification of heat stress-responsive genes in heat-adapted thermal Agrostis scabra by suppression subtractive hybridization. J Plant Physiol. 2009, 166: 588-601. 10.1016/j.jplph.2008.09.003.PubMedView ArticleGoogle Scholar
- Giorno F, Wolters-Arts M, Grillo S, Scharf KD, Vriezen WH, Mariani C: Developmental and heat stress-regulated expression of HsfA2 and small heat shock proteins in tomato anthers. J Exp Bot. 2010, 61: 453-462. 10.1093/jxb/erp316.PubMed CentralPubMedView ArticleGoogle Scholar
- Li S, Li F, Wang J, Zhang W, Meng Q, Chen TH, Murata N, Yang X: Glycinebetaine enhances the tolerance of tomato plants to high temperature during germination of seeds and growth of seedlings. Plant Cell Environ. 2011, 34: 1931-1943. 10.1111/j.1365-3040.2011.02389.x.PubMedView ArticleGoogle Scholar
- Granger D, Höllwarth M, Parks D: Ischemia-reperfusion injury: role of oxygen-derived free radicals. Acta Physiol Scand Suppl. 1985, 548: 47-63.Google Scholar
- Hu T, Li H, Zhang X, Luo H, Fu J: Toxic effect of NaCl on ion metabolism, antioxidative enzymes and gene expression of perennial ryegrass. Ecotox Environ Safe. 2011, 74: 2050-2056. 10.1016/j.ecoenv.2011.07.013.View ArticleGoogle Scholar
- Hoagland D-R, Arnon DI: Circular. Calif Agr Exp Station. 1950, The water-culture method for growing plants without soi, 2, 347Google Scholar
- 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.PubMedView ArticleGoogle Scholar
- Palaniswamy SK, James S, Sun H, Lamb RS, Davuluri RV, Grotewold E: AGRIS and AtRegNet. a platform to link cis-regulatory elements and transcription factors into regulatory networks. Plant Physiol. 2006, 140: 818-829. 10.1104/pp.105.072280.PubMed CentralPubMedView ArticleGoogle Scholar
- Saldanha AJ: Java Treeview—extensible visualization of microarray data. Bioinformatics. 2004, 20: 3246-3248. 10.1093/bioinformatics/bth349.PubMedView ArticleGoogle Scholar
- Kakumanu A, Ambavaram MM, Klumas C, Krishnan A, Batlang U, Myers E, Grene R, Pereira A: Effects of drought on gene expression in maize reproductive and leaf meristem tissue revealed by RNA-Seq. Plant Physiol. 2012, 160: 846-867. 10.1104/pp.112.200444.PubMed CentralPubMedView ArticleGoogle Scholar
- Lee JM, Roche JR, Donaghy DJ, Thrush A, Sathish P: Validation of reference genes for quantitative RT-PCR studies of gene expression in perennial ryegrass (Lolium perenne L.). BMC Mol Biol. 2010, 11: 8-10.1186/1471-2199-11-8.PubMed CentralPubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.