Skip to main content

Transcriptomic analysis of grape (Vitis viniferaL.) leaves during and after recovery from heat stress



Grapes are a major fruit crop around the world. Heat stress can significantly reduce grape yield and quality. Changes at the molecular level in response to heat stress and subsequent recovery are poorly understood. To elucidate the effect of heat stress and subsequent recovery on expression of genes by grape leaves representing the classic heat stress response and thermotolerance mechanisms, transcript abundance of grape (Vitis vinifera L.) leaves was quantified using the Affymetrix Grape Genome oligonucleotide microarray (15,700 transcripts), followed by quantitative Real-Time PCR validation for some transcript profiles.


We found that about 8% of the total probe sets were responsive to heat stress and/or to subsequent recovery in grape leaves. The heat stress and recovery responses were characterized by different transcriptional changes. The number of heat stress-regulated genes was almost twice the number of recovery-regulated genes. The responsive genes identified in this study belong to a large number of important traits and biological pathways, including cell rescue (i.e., antioxidant enzymes), protein fate (i.e., HSPs), primary and secondary metabolism, transcription factors, signal transduction, and development. We have identified some common genes and heat shock factors (HSFs) that were modulated differentially by heat stress and recovery. Most HSP genes were upregulated by heat stress but were downregulated by the recovery. On the other hand, some specific HSP genes or HSFs were uniquely responsive to heat stress or recovery.


The effect of heat stress and recovery on grape appears to be associated with multiple processes and mechanisms including stress-related genes, transcription factors, and metabolism. Heat stress and recovery elicited common up- or downregulated genes as well as unique sets of responsive genes. Moreover, some genes were regulated in opposite directions by heat stress and recovery. The results indicated HSPs, especially small HSPs, antioxidant enzymes (i.e., ascorbate peroxidase), and galactinol synthase may be important to thermotolerance of grape. HSF30 may be a key regulator for heat stress and recovery, while HSF7 and HSF1 may only be specific to recovery. The identification of heat stress or recovery responsive genes in this study provides novel insights into the molecular basis for heat tolerance in grape leaves.


Most crop plants are exposed to heat stress during certain stages of their life cycle. Heat stress, defined as the temperature above a normal optimum, is expected to become a major issue in reducing crop production in coming years due to global warming [1]. Grape is a popular cultivated fruit throughout the world and represents one of the most important crops with highly valued products such as juices, liquors and wines [2]. The grape species Vitis vinifera makes up most of the grape production in the world. However, grape production and quality often fluctuate due to various environmental factors. Temperature has been broadly considered as a major determining factor. Studies show that crop production is severely limited by temperature stresses around the world [3]. In many regions, the maximum midday air temperature can reach 40°C and above, which can destroy grape berry ripening [4]. In addition, crop cultivation in sheltered conditions (e.g., greenhouses and hoop houses) is common in many regions. These conditions can further increase temperature due to inadequate air circulation. Studies indicated that temperatures above 35°C generally reduce photosynthesis in grape leaves [5]. Extreme temperatures may endanger berry quality and economic returns [4, 6]. This is expected to get worse with more frequent high temperature stress due to climate change [3, 7].

In the past, studies of response and adaptation of grape to high temperatures have focused mostly on grape morphological and physiological changes including photosynthesis, respiration, cell membrane stability, hormone changes and antioxidant systems [6, 813]. With the availability of the grape genome sequence [14, 15], study of the functional genomics of grapes has become possible [3, 16]. Transcriptomic analysis represents one of these major research opportunities. Gene expression is tissue- and development-specific. To date, transcriptomic studies in grapes has been primarily focused on berry development and water stress responses [1729]. There are few transcriptomic reports of the effect of heat stress on grape.

Transcriptomic studies of heat stress effects on Arabidopsis, rice, tobacco, potato, tomato and sunflower have been reported [3034]. Lim et al. [30] found that acclimating Arabidopsis thaliana suspension cells at a moderate heat enhanced heat resistance. Expression of 165 genes changed, especially those of heat shock proteins (HSPs). With cDNA microarrays and RT-qPCR techniques, Frank et al. [33] found that HSP70, HSP90, and heat shock transcription factors(HSF) HSFA2 and HSFA3 were important to tomato microspore resistance to heat stress. However, data on the molecular mechanisms involved in heat stress responses and thermotolerance in grape leaves are very limited.

Some studies have shown that the recovery process from heat stress in plants is very important to survival [35]. The degree of recovery from stress is a direct index of plant heat tolerance. Thus, there should be some differences between the recovery mechanisms and the direct heat response mechanisms in plants [36]. Gu et al. [37] found that sucrose synthase, calmodulin, certain peptides, and aquaporin genes in Populus were transcriptionally-activated only during recovery from salt stress.

In this study, we used oligonucleotide microarrays and quantitative Real-Time PCR (qRT-PCR) to identify genes in grape leaves with altered transcript accumulation during heat stress and after recovery to provide clues to the function of these genes during heat stress and subsequent recovery.


Expression and validation of probes sets responsive to heat stress or recovery in grape leaves

Transcriptomic profiles of gene expression variation in Vitis vinifera cv. Cabernet Sauvignon leaves in response to heat treatment and subsequent recovery was quantitatively assessed using the Affymetrix Grape Genome Array with 15,700 probe sets. Array data were averaged for three biological replicates and filtered as described in Materials and Methods. Using the filtering criteria, 1282 (about 8% of total probe sets) were significantly affected by heat or recovery and were further analyzed. Heat stress and recovery affected transcript levels in different ways. Among those probe sets showing differential expression during heat stress (dChip, P < 0.001), 247 were upregulated, 697 were downregulated, while 150 were induced and 353 were repressed after recovery compared to their corresponding control levels (Figure 1). A total of 12 probe sets were commonly upregulated by heat stress and recovery, whereas 179 and 114 were specifically induced by both treatments, respectively. A total of 73 probe sets were commonly downregulated by both treatments, whereas 600 and 224 were specifically repressed by heat stress and recovery, respectively. Moreover, 24 probe sets were repressed by heat stress but induced by recovery, and 56 probe sets were upregulated by heat stress and downregulated by recovery. Cluster analysis also indicated that some genes were up- or downregulated under heat stress and after recovery (Figure 2). In order to validate the results obtained with the microarray analyses, we carried out qRT-PCR assays on 12 cDNA sequences using gene-specific primers (Additional file 1). The qRT-PCR profiles were analyzed on three biological replicates. Linear regression analyses displayed highly significant correlations (r = 0.982) between qRT-PCR and microarray results for the 12 evaluated genes (see Additional file 2), confirming the validity of the microarray results.

Figure 1

Venn diagram of transcripts (both identified and unknown) that were up- or downregulated by heat stress or recovery. The " + " and "-" indicate up- and downregulated transcripts, respectively. A total of 1282 transcripts were significantly (dChip, P < 0.001) affected by heat or recovery. 179: unique upregulated transcripts by heat stress; 114: unique upregulated transcripts by the recovery; 600: unique downregulated transcripts by heat stress; 224: unique downregulated transcripts by the recovery; 12: commonly upregulated transcripts by heat stress and recovery; 73: commonly downregulated transcripts by heat stress and recovery; 24: downregulated by heat stress but upregulated by recovery; 56: upregulated by heat stress but downregulated by recovery.

Figure 2

Average linkage hierarchical clustering analysis of the Log 2 transformed values of fold changes of the 821 annotated genes during heat stress and after the following recovery. HS-1, HS-2 and HS-3 represent three replications of heat stress treatments; RC-1, RC-2 and RC-3 represent three replications of the recovery.

Functional analysis of probes sets responsive to heat stress or recovery in grape leaves

There were a large number of genes which were responsive to heat stress or recovery that do not match any genes with known functions. Eight hundred and twenty-one probe sets were assigned to annotated genes and ESTs in PLEXdb ( Functional categories of these probe sets based on MIPS (Munich Information Center for Protein sequences) are shown in Figure 3. Furthermore, we compared the common probe sets between heat stress and recovery, and analyzed probe sets unique to both treatments, based on functional classification. The changes of some genes such as HSPs were similar to those in previous studies on microarray analysis of heat stress in rice, Arabidopsis, tomato, barley and potato [30, 31, 33, 3841].

Figure 3

Functional classification of heat stress and recovery–responsive transcripts. Transcripts with a mean absolute expression ratio of at least 2.0 (linear scale) and a P<0.001 in a t-test for significance were classified into the categories shown, based on shared putative function and/or common structural motifs. HS heat stress; RC recovery.

Comparative analysis of common responsive probe sets between heat stress and subsequent recovery

There were 40 annotated genes that were upregulated by heat treatment, but were downregulated following recovery. They are listed in groups according to their putative involvement in different cellular events (see Additional file 3). Sixty percent of them were HSPs/molecular chaperones and may be critical to heat tolerance. The genes of high and middle molecular HSPs included HSP62.6, HSC71.5, HSP70, HSP 101, HSP80.1, low molecular HSPs genes included CPN10, HSP11.1, HSP15.7, HSP16.1, HSP17.5, HSP17.6, HSP18.3, HSP18.6, HSP22, HSP25.7, HSP37.1 and HSP40. Interestingly, gene HSF30 was upregulated 11-fold by heat stress, but downregulated more than 10-fold following recovery compared with the corresponding controls. In addition, genes related to cell rescue including Hin1, Stil and galactinol synthase changed. Three calmodulin-related proteins that suppress posttranscriptional gene silencing in plants were also involved in the process. Also affected by heat stress were genes of ethylene-responsive transcriptional co-activator-like protein, metabolism-associated enzymes, 3-hydroxy-3-methylglutaryl coenzyme A (acetyl-CoA pathway) and ER1(pyridoxin biosynthesis protein), development -related enzymes, and ripening regulated protein.

There were a total of 10 annotated genes upregulated by heat stress and subsequent recovery (see Additional file 4). HSP20 showed a larger degree of up-regulation by heat stress than by recovery. The other genes were related to metabolism, transcription, cell rescue and transposable elements.

Additional file 5 shows 19 genes which were downregulated by heat stress but upregulated during subsequent recovery. Genes related to cell rescue included peroxidase 42 and alcohol dehydrogenase 7. A number of genes were related to metabolism; for example, genes involved in starch biosynthesis (ADP-glucose pyrophosphorylase large subunit 1), polysaccharide degradation (beta-galactosidase BG1), and polyphenol oxidase (tyrosine metabolism). Aquaporin and triose phosphate were affected and are related to transport regulation. In addition, the affected gene F14J22.4 is involved in signal transduction. No genes in the regulatory factor categories were identified.

Additional file 6 shows 53 genes that were downregulated by heat stress and subsequent recovery. Twelve of these genes were related to cell rescue including stress responses. Two Class IV chitinase transcripts were downregulated about 30-fold by heat treatment, but only downregulated 3-fold after the recovery. Other downregulated genes were involved in metabolism, signal transduction, transcript regulation and cell transport. Secondary metabolism genes that were downregulated included stilbene synthase, flavonol 3-O-glucosyltransferase 2, caffeic acid O-methyltransferase. Regulatory genes including WRKY transcription factor-b, Myb-related transcription factor MybB1-2, SPF1 protein (transcription activation), and NAC domain protein NAC1 were also downregulated.

Analysis of unique probe sets responsive to heat stress in grape leaves

In Additional file 7, a total of 78 upregulated probe sets identified as unique to the heat stress response are listed. Stress-responsive genes affected included those for one salt-inducible protein, one wound-induced protein, a dehydroascorbate reductase, and a cytosolic ascorbate peroxidase (APX). There were nine HSPs/molecular chaperones influenced, including HSP21 which was upregulated 43-fold. Metabolism genes affected included tetrahydropicolinate succinylase, pectinesterase, homocysteine S-methyltransferase 2, putative ripening-related protein sterol 4-alpha-methyl-oxidase and adenine phosphoribosyl transferase. Four genes involved in signal transduction were affected including those for GTP-binding, a ADP-ribosylation factor-like protein, a receptor associated protein (transmembrane signal transduction) and strubbelig receptor family 3. Regulatory genes responding included genes for zinc finger protein-like, Myb-related protein, and triptychon. Metabolism genes affected included those for fatty acid biosynthesis.

There were 600 downregulated probe sets uniquely responsive to heat stress (Figure 1), of which 412 were identified (see Additional file 7). Stress response genes included those for chitinase, pathogenesis-related protein, osmotic, and cold stress response proteins. A total of 129 genes related to primary and secondary metabolism were downregulated by heat treatment. The transcript level of 13 genes involved in respiration and photosynthesis declined. Transcription regulatory genes that were down-regulated included those for zinc finger protein, pathogenesis-related gene transcriptional activator PTI5, Myb-related transcription factor VlMYBB1-2, ethylene response factor ERF3a, bHLH106, and MYC transcription factor. Down-regulated genes involved in transport regulation included cation, carbohydrate, amino acid, protein and nucleotide transporters.

Analysis of unique responsive probe sets to recovery from heat stress

There were 72 upregulated probe sets unique to the recovery response (see Additional file 8). Stress genes that were upregulated included adenylyl-sulfate reductase precursor, Fe-SOD, late embryogenesis abundant protein D-29, and putative pathogenesis-related protein. Upregulated genes involved in protein fate only included U-Box protein and aspartyl protease family protein, but no HSPs were observed. There were 27 genes related to metabolism that were upregulated, involving carboxylate metabolism, biosynthesis of phenylalanine, phosphate metabolism, polysaccharide degradation, and fatty acid biosynthesis. There was only one signal transduction gene, a receptor protein kinase and a protein ralf-like 34 that were upregulated. There were 17 transporter genes upregulated after the recovery, including a nitrate transporter, a tonoplast intrinsic protein and several ABC transporters. Transcript regulatory genes upregulated including a squamosa promoter binding-like protein, a Myc-like anthocyanin regulatory protein, a putative Zinc finger protein, and a translation initiation factor IF2.

In addition, 137 downregulated probe sets were unique to the recovery response (see Additional file 8). Stress- responsive proteins that were downregulated were mainly defence proteins. Downregulated genes involved in protein fate included HSPs, HSFs (HSF1, HSF7), and protein kinase domain-containing protein. Downregulated genes related to metabolism numbered 34, and were involved in biosynthesis of glutamine, stilbenes, flavonoids and metabolism of proline and carbohydrates. Some downregulated genes were related to Ca2+ mediated signal transduction, and some were regulatory genes including NAC domain proteins, WRKY transcription factors, and Zinc finger DNA-binding protein.


High temperature disturbs cellular homeostasis in plants and can lead to severe retardation of growth and development and even death. Plants exposed to heat stress exhibit a characteristic set of cellular and metabolic responses. Based on the results detailed above, the discussion focuses on the response of grape leaves to heat stress and subsequent recovery from the following aspects.

HSPs and antioxidant enzymes

One typical response to heat stress is an accelerated transcription of a set of stress or protein fate-related genes, such as those encoding HSPs, which are the main constituents of the heat shock response [1]. The major HSPs belong to five structurally distinct classes: Hsp100, Hsp90, Hsp70, Hsp60 and small HSPs (sHSPs) [42]. In this study, heat stress increased the expression of various HSPs. The transcript level of these HSPs decreased to the control level or even lower after the subsequent recovery (see Additional files 3 and 7). Under heat stress, many proteins in the cell are subject to denaturation, and HSPs function as molecular chaperones to provide protection. After the recovery, these HSPs are degraded. Moreover, a large amount of HSPs in cells resulted in the decline of transcript level of HSPs genes [43]. Wang et al. [44] reported that the expression of protein HSP21 was induced by heat stress in grapevine leaves, and was lower than that of the control after the recovery. In plants, sHSPs form a more diverse subfamily than other HSPs/chaperones with respect to sequence similarity, cellular location and function. The sHSPs are not themselves able to refold non-native proteins, but have a high capacity to bind non-native proteins, possibly through hydrophobic interactions, and to stabilize and prevent non-native aggregation, thereby facilitating their subsequent refolding by ATP-dependent chaperones such as the DnaK system or ClpB/DnaK. Increasing evidence suggests strong correlation between sHSP accumulation and plant tolerance to stress [45]. Chloroplast and mitochondrial sHSPs are considered to play an important role in heat tolerance [46, 47]. In this study, HSP22 located in mitochondria was highly expressed under heat stress and was downregulated after recovery. HSP21 located in the chloroplast was highly expressed (43-fold) by heat stress, but its transcription level declined to the control level after the following recovery (see Additional file 7). This result shows that some sHSPs may have important effects on heat tolerance of grapevines.

Among the other HSPs, over expression of HSP101 in Arabidopsis had a positive effect on growth after recovery [36]. It was recently found that an HSP101 homologue in Arabidopsis was involved in conferring thermotolerance to chloroplasts during heat stress [48]. Genetic analysis in Arabidopsis indicated that HSP101 interacts with the sHSP chaperone system to re-solubilize protein aggregates after heat stress [49]. Recently, it was demonstrated that the transcript level of HSP101 increased in maturing tomato pollen grains in response to heat stress [33, 50]. In the present study, HSP101 in grape leaves was upregulated by heat stress, exhibiting a 9-fold elevated expression, but it was downregulated 6-fold after the subsequent recovery. HSP70 has essential functions in preventing aggregation and in assisting refolding of non-native proteins under both normal and stress conditions [51]. Some family members of HSP70 are consistently expressed and are often referred to as HSC70. These members are often involved in assisting the folding of de novo synthesized polypeptides and the import/translocation of precursor proteins. In this study, the expression of HSC70 was induced by heat stress and declined after the subsequent recovery, in agreement with previous results.

In heat stress studies, increasing attention has being paid to the generation of reactive oxygen species (ROS) and the cellular antioxidant defense systems. ROS levels are controlled by a network of enzymes and metabolites, including superoxide dismutases (SOD), ascorbate peroxidase (APX), guaiacol peroxidase (GPX) and dehydroascorbate reductase (DHAR). APX plays a pivotal role in ROS metabolism. It catalyzes the reduction of hydrogen peroxide to water by using ascorbate as a specific electron donor [52]. APX appears to be regulated by HSFA2 [53]. Previous studies also indicated that it is involved in survival from high light stress [54]. The results of the present study showed that expression of APX, peroxidase 42 and DHAR were upregulated by 5-, 3.27- and 3.18-fold, respectively, indicating that these genes may have an important role in grape leaves in response to heat stress.


Temperature is one of the most active environmental factors affecting all plant metabolic activities, including amino acid and carbohydrate metabolism [55]. Secondary metabolites are involved in resistance against heat shock [56]. In this study, galactinol synthase transcript level increased significantly (about 50-fold) in response to heat stress and declined during recovery. Heat shock induced the production of sugars including raffinose and galactinol [57]. The raffinose family of oligosaccharides has been implicated in the scavenging of hydroxyl radicals [58], and in protecting liposomes from desiccation through direct sugar–membrane interactions in soybean [59]. Galactinol synthase catalyzes the first step in the synthesis of raffinose polysaccharide, and is also regulated by HSFA2 in Arabidopsis, linking heat shock proteins and raffinose metabolism [60]. Weston et al. [61] reported that expression of AtGolS1 was observed even at optimal temperature, and was upregulated during heat stress. Several other probe sets representing genes related to carbohydrate metabolism were regulated by heat stress in our experiment, such as invertase, UDP-glucose dehydrogenase, 6-phosphate dehydrogenase, sucrose synthase, amylase and trehalose-6-phosphate synthase/phosphatase (see Additional files 3 and 7). These enzymes are important in sugar metabolism, which were downregulated by heat stress in this study. In addition, lipid metabolism was also inhibited by heat stress, with effects on lipase, L-asparaginase, lipoxygenase and fatty acid hydroperoxide lyase.

Transcription factors (TFs)

HSFs play important roles in both basal and acquired thermotolerance through binding to cis-acting regulatory elements called heat shock element (HSEs) in the promoter region of HSP genes [62]. In general, plant HSFs are divided into three classes, HSFAs, HSFBs and HSFCs [63]. HSF30 belongs to HSFA2 which was the dominant HSF in thermotolerant cells [39, 64] and was highly heat shock upregulated in mature tomato microspores [33]. The transcript levels of the HSFA2 target genes (e.g. HSP101, HSP70, HSP22, HSP17.6, and HSP15.7) were highly correlated with those of HSFA2 in Pro35S:HSFA2 Arabidopsis plants, and the induction of HSFA2 target genes was strongly reduced under heat stress in HSFA2 knockout Arabidopsis plants [60]. In this study, HSF30 exhibited increased expression levels (11-fold) under heat stress and was downregulated (14.3-fold) after the subsequent recovery corresponding to 16 HSPs (see Additional file 3). The present results indicated that HSF30 may also play an important regulatory role in the thermotolerance of grape. HSF7 is an HSFb2b and was strongly induced by heat stress in Arabidopsis, maize and tomato [64]. However, HSF7 was weakly upregulated (1.7-fold) by heat stress, but showed significant down-regulation (3.7-fold) after recovery in the present study. This indicated that HSF7 may play an important role in reducing expression of HSPs after recovery in grape leaves. In addition, HSF1 governs the expression of HSPs and regulates thermotolerance [33, 65]. HSF1 belongs to the HSFA1a group, which was identified as a master regulator of thermotolerance in tomato [64]. The synthesis of members of HSP100, HSP90, HSP70, HSP60 and sHSPs under heat stress in the HSF1-RNAi strains of Chlamydomonas was dramatically reduced or completely abolished [66]. In the present study, HSF1 showed similar sign of change as HSF7. Therefore, HSF7 and HSF1 may play important roles in the recovery process from heat stress in grape leaves.

Qin et al. [67] reported that ethylene-responsive transcriptional co-activator (ERTCA) gene is upregulated more than 8-fold in all heat treated wheat leaves. Over-expression of ERTCA enhanced heat stress tolerance of Arabidopsis [68]. Heat stress induced expression of ERTCA rapidly in the sensitive genotype of tomato [38]. In the present study, the strong induction of ERTCA expression by heat stress in grape leaves provided another piece of evidence for its role in heat tolerance. Many other transcription factor genes were also affected by heat stress and recovery in grape leaves, although their roles in heat tolerance are not clear. Interestingly, the majority of heat response transcription factors genes were responsive to the heat treatment. Most of them were downregulated. Basic leucine zipper (bZIP) transcription factors play a role in plant pathogen responses, light signaling, and ABA and abiotic stress signaling [69]. Our microarray data revealed that bZIP transcription factors were heat-regulated. Some were upregulated, such as Zinc finger protein-like, some were downregulated, such as B-box type Zinc finger-containing protein. GATA-type Zinc finger protein was upregulated only by the recovery treatment (see Additional file 8). Plant transcription factors WRKYs have been reported in both biotic and abiotic stress responses [70]. In this study, three WRKY transcription factors were downregulated by heat stress, and two WRKY transcription factors were downregulated by the recovery (see Additional files 7 and 8). These results suggested that several WRKY factors could be involved in the heat response of grape leaves. Our data indicated that expression of 3 transcription factors with NAC domain genes were also heat–regulated. Plant-specific NAC family transcription factor has a conserved NAC domain at the N-terminal of the protein and has been implicated in plant development [71]. It was recently reported that several NAC transcription factors were also involved in biotic and abiotic stress response [72].

Signal transduction components

In our study, 4 genes for receptor-like kinases (RLKs) were regulated by heat stress (see Additional file 7). RLK1 is induced by wounding, pathogen attack, and salicylic acid [73]. Recent work indicated that RLK1 plays an important role in abscisic acid (ABA) signal transduction [74]. Our results suggested that heat signal transduction in grape leaves shared, at least in part, some common pathways with other biotic, abiotic and ABA stress signaling through these RLKs. Protein phosphorylation and dephosphorylation have been reported in heat signal transduction. Indeed, protein kinases and phosphatases with altered expression formed the largest group of genes. In addition, many protein phosphatases showed differential expression, suggesting that protein post-translation modification occurred during the heat response of grapevine leaves.

Calcium is a universal signaling molecule in both animals and plants, and the transient increase of Ca2+ level during heat stress is well-documented in plants [11, 75, 76]. Heat shock triggered cytosolic Ca2+ bursts, which is transferred by Ca2+ binding proteins (CBP) such as calmodulin (CaM), CaM-related proteins, Ca2+-dependent protein kinases (CDPK), and calcineurin B-like protein (CBL), and then unregulated the expression of HSPs, due to the dependence of the final step in HSF- mediated HSP expression on a Ca2+ signal [77, 78]. In the present analysis, candidate genes encoding the components of calcium- or calmodulin-mediated signal pathways, including calnexin, calmodulin, and CDPKs, were also heat- or recovery-regulated, suggesting a role of Ca2+-mediated signals in the heat stress response. Plant survival after severe environmental stress largely depends on the efficiency of recovery mechanisms. Among the genes activated after recovery, we found a calmodulin gene. qRT-PCR indicated a 4.28-fold induction of calmodulin transcript accumulation after the early response to salt stress [79].


The data presented here provides genome-wide expression profiles of grape leaves under heat stress and subsequent recovery. Affymetrix Grape Genome Array and qRT-PCR techniques were used to identify heat stress- and recovery-regulated genes which represented classic heat stress responsive and thermotolerance mechanisms. The present study highlighted the significant contribution and fundamental roles of transcriptional control in stress responses in grape leaves. We found that about 8% of total probe sets were responsive to heat stress and subsequent recovery in grape leaves. Heat stress and recovery responses displayed differential gene expression changes. The number of heat stress-regulated genes was almost twice that of recovery-regulated genes. The responsive genes identified in this study belong to a large number of important factors and biological pathways, including those for cell rescue (i.e., antioxidant enzymes), protein fate (i.e., HSPs), primary and secondary metabolism, transcription factors, signal transduction and development. Of particular interest, HSPs especially sHSPs, APX and galactinol synthase, may be very important to thermotolerance of grape leaves; HSF30 may be a key regulator for heat stress and recovery;, and, HSF7 and HSF1 may mainly function after recovery. These results provide novel insight into the grape leaf response to heat stress and have great implications for further studies on gene function annotation and molecular breeding.


Plant materials and treatments

Stem cuttings of ‘Cabernet Sauvignon’ (Vitis vinifera L.) were rooted in pots containing a mixture of 4 peat moss: 6 perlite (V/V) and grown in a greenhouse under mist conditions. When the cuttings were rooted, they were repotted into larger pots, and grown for about 10 weeks in a greenhouse at 70–80% relative humidity. During the daytime, temperature was at most 25°C and the maximum photosynthetic active radiation (PAR) was about 1,000 μmol m-2 s-1. During night time, the temperature was maintained above 18°C. Young grapevines with identical growth (10 leaves) were acclimated for two days in a controlled environment room (70 - 80% relative humidity, 25/18°C day/night cycle and PAR at 800 μmol m-2 s-1) and divided into two groups. On the following day (the first day of the experiment, Day 1), one group of grapevines was kept at 25/18°C day/night in the controlled environment room as the control. The other group was treated at 45°C in another controlled environment room (except for temperature, the other conditions were the same as the control) from 9:00 am to 14:30 pm. The stressed grapevines were then allowed to recover at 25°C rapidly (from 45°C to 25°C for at about 15 min). Then all conditions were the same as the control until 9:30 am on day 2. Leaf samples of the treatment and control were taken at 14:30 pm on Day 1 and 9:00 am on day 2 for transcriptomic analyses. Three independently replicated experiments were executed.

RNA extraction, amplification, labeling and hybridization

Total RNA were extracted from grape leaves using Trizol reagent (Invitrogen, Carlsbad, CA) according to the manufacturer’s instructions, and digested with DNase I at 37°C for 15 min to remove any contaminating DNA. The RNA was cleaned up with RNeasy Kit (Qiagen, Hilden, Germany) and the quantities and qualities were determined by spectrophotometry and 1% formaldehyde denaturing gel electrophoresis. The samples with bright bands of ribosomal 28S to 18S RNA in a ratio >1.5:1 were used for microarray analysis. The Affymetrix GeneChip Vitis vinifera (Grape) Genome Array, which contains 15,700 probe sets to cover 14,000 V. vinifera transcripts and 1,700 transcripts from other Vitis species, was used for microarray analysis. Hybridization, data capture, and analysis were performed by CapitalBio Corporation (Beijing, China), a service provider authorized by Affymetrix Inc. (Santa Clara, CA). Briefly, 200 ng of total RNA was used for cDNA synthesis, and produce biotin-tagged cRNA with MessageAmpTM Premier RNA Amplification Kit (Ambion). A 10 μg fragmented cRNA, with control oligo B2 and eukaryotic hybridization controls (bioB, bioC, bioD, cre), was hybridized to each GeneChip array at 45°C for 16 h (Affymetrix Gene Chip Hybridization Oven 640) according to the manufacturer’s instructions. After hybridization, the GeneChip arrays were washed, and then stained with streptavidin phycoerythrinonan (SAPE) with an Affymetrix Fluidics Station 450 followed by scanning with an Affymetrix GeneChip Scanner 3000 7 G.

Microarray data processing

Microarray image files (CEL) were loaded into a DNA-Chip Analyzer (dChip) software package [80]. The model-based expression values were calculated in dChip, followed by normalization via the program default method (Invariant Set Normalization). An unpaired two-group comparison for all probe sets was performed. Genes were determined to have altered expression levels in the treated samples versus control samples based on the following criteria: (1) P-value < 0.001; (2) Fold change Treatment/Control >2 or <0.5. The lower confidence bound (LCB) of the 90% confidence interval of the fold changes was used [81]; (3) An absolute difference between the means of the expression levels of the two groups was greater than 50; (4) Present calls in both samples are larger than 20%. The reliability of the comparison criteria was assessed by checking the FDR (False Discovery Rate) when permuting samples 200 times. Genes that satisfied all of the above criteria were chosen for further analysis. The gene annotation and functional categories were performed using data at PLEXdb ( and based upon Deluc et al. [82]. All data have been submitted to the PLEXdb, and the accession number is VV40.


Total RNA extraction was the same as that used for microarray analysis described above. The total RNA was treated with DNase I (Promega) to avoid DNA contamination. One microgram of RNA was reverse transcribed using the Superscript II reverse transcriptase (Invitrogen) with an oligo(dT)15 primer according to the manufacturer’s instructions (Tiangen Biotech, Beijing, China). qRT-PCR experiments were conducted using Real Master Mix (SYBR Green) (Tiangen Biotech, Beijing, China). Reactions were carried out on a MX 3000 multicolor real-time detection system. The following standard thermal profile was used for all PCR experiments: 94°C for 2 min; 40 cycles of 95°C for 15 s, 56°C for 18 s and 68°C for 20 s. Fluorescence signals were captured at the end of each cycle, and the melting curve analysis was performed from 68°C to 95°C. Gene-specific primers were designed using the Primer5 software (see Additional file 1). The amplification of 18S rRNA gene sequence (GQ849399) was used as the internal control to normalize all the data [83, 84]. Analyses of qRT-PCR data used the classic (1 + E)-ΔΔCT method (CT is the threshold cycles of one gene, E is the amplification efficiency). ΔCT is equal to the difference in threshold cycles for target (X) and reference (R) (CT,X-CT,R), while the ΔΔCT is equal to the difference of ΔCT for control (C) and treatment (T) (ΔCT,T-ΔCT,C). The amplification system (e.g., primer and template concentrations) was properly optimized, and the efficiency was close to 1. So the amount of target, normalized to an endogenous reference and relative to a calibrator, is given by:

Amount of target =2-ΔΔCT.


  1. 1.

    Wahid A, Gelani S, Ashraf M, Foolad MR: Heat tolerance in plants: An overview. Environmental and Experimental Botany. 2007, 61: 199-223. 10.1016/j.envexpbot.2007.05.011.

    Article  Google Scholar 

  2. 2.

    Vivier MA, Pretorius IS: Genetically tailored grapevines for the wine industry. Trends Biotechnol. 2002, 20: 472-478. 10.1016/S0167-7799(02)02058-9.

    PubMed  CAS  Article  Google Scholar 

  3. 3.

    Cramer GR, Urano K, Delrot S, Pezzotti M, Shinozaki K: Effects of abiotic stress on plants: a systems biology perspective. BMC Plant Biol. 2011, 11: 163. 10.1186/1471-2229-11-163.

    PubMed  PubMed Central  Article  Google Scholar 

  4. 4.

    Schultz HR: Abiotic stress ecophysiology and grape functional genomics. Climate change and world viticulture. Cost Action 858 Workshop: Vineyard under environmental constraints: adaptations to climate change. 2007, Poland: University of Lodz.

    Google Scholar 

  5. 5.

    Kriedemann PE: Photosynthesis in vine leaves as a function of light intensity, temperature, and leaf age. Vitis. 1986, 7: 213-220.

    Google Scholar 

  6. 6.

    Howell GS: Sustainable grape productivity and the growth-yield relationship: A review. Am J Enology and Viticulture. 2001, 52: 165-174.

    Google Scholar 

  7. 7.

    van Leeuwen C, Friant P, Chone X, Tregoat O, Koundouras S, Dubourdieu D: Influence of climate, soil, and cultivar on terroir. Am J Enology and Viticulture. 2004, 55: 207-217.

    Google Scholar 

  8. 8.

    Sepu´lveda G, Kliewer WM: Stomatal response of three grapevine cultivars (Vitis vinifera L.) to high temperature. Am J Enol Vitic. 1986, 37: 44-52..

    Google Scholar 

  9. 9.

    Caprio JM, Quamme HA: Weather conditions associated with grape production in the Okanagan Valley of British Columbia and potential impact of climate change. Can J Plant Sci. 2002, 82: 755-763. 10.4141/P01-160.

    Article  Google Scholar 

  10. 10.

    Mori K, Goto-Yamamoto N, Kitayama M, Hashizume K: Loss of anthocyanins in red-wine grape under high temperature. J Exp Bot. 2007, 58: 1935-1945. 10.1093/jxb/erm055.

    PubMed  CAS  Article  Google Scholar 

  11. 11.

    Wang LJ, Li SH: Salicylic acid-induced heat or cold tolerance in relation to Ca2+ homeostasis and antioxidant systems in young grape plants. Plant Sci. 2006, 170: 685-694. 10.1016/j.plantsci.2005.09.005.

    CAS  Article  Google Scholar 

  12. 12.

    Wang LJ, SH L: Heat acclimation induced acquired heat tolerance and cross adaptation in different grape cultivars: Relationships to photosynthetic energy partitioning. Functional Plant Biology. 2009, 36: 516-526. 10.1071/FP09008.

    CAS  Article  Google Scholar 

  13. 13.

    Luo HB, Ma L, Xi HF, Duan W, Li SH, Loescher W, Wang JF, Wang LJ: Photosynthetic responses to heat treatments at different temperatures and following recovery in grapevine (Vitis amurensis L.) leaves. PLoS One. 2011, 6: e23033. 10.1371/journal.pone.0023033.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  14. 14.

    Jaillon O, Aury JM, Noel B, Policriti A, Clepet C, Casagrande A, Choisne N, Aubourg S, Vitulo N, Jubin C, et al: The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla. Nature. 2007, 449: 463-467. 10.1038/nature06148.

    PubMed  CAS  Article  Google Scholar 

  15. 15.

    Velasco R, Zharkikh A, Troggio M, Cartwright DA, Cestaro A, Pruss D, Pindo M, FitzGerald LM, Vezzulli S, Reid J, et al: A high quality draft consensus sequence of the genome of a heterozygous grapevine variety. PLoS One. 2007, 2: e1326. 10.1371/journal.pone.0001326.

    PubMed  PubMed Central  Article  Google Scholar 

  16. 16.

    Cramer GR: Abiotic stress and plant responses from the whole vine to the genes. Australian J Grape and Wine Res. 2010, 16: 86-93.

    CAS  Article  Google Scholar 

  17. 17.

    Deluc LG, Grimplet J, Wheatley MD, Tillett RL, Quilici DR, Osborne C, Schooley DA, Schlauch KA, Cushman JC, Cramer GR: Transcriptomic and metabolite analyses of Cabernet Sauvignon grape berry development. BMC Genomics. 2007, 8: 429. 10.1186/1471-2164-8-429.

    PubMed  PubMed Central  Article  Google Scholar 

  18. 18.

    Deluc LG, Decendit A, Papastamoulis Y, Merillon JM, Cushman JC, Cramer GR: Water deficit increases stilbene metabolism in Cabernet Sauvignon berries. J Agric Food Chem. 2011, 59: 289-297. 10.1021/jf1024888.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  19. 19.

    Cramer GR, Ergul A, Grimplet J, Tillett RL, Tattersall EAR, Bohlman MC, Vincent D, Sonderegger J, Evans J, Osborne C, et al: Water and salinity stress in grapevines: early and late changes in transcript and metabolite profiles. Funct Integr Genomics. 2007, 7: 111-134. 10.1007/s10142-006-0039-y.

    PubMed  CAS  Article  Google Scholar 

  20. 20.

    Grimplet J, Deluc LG, Tillett RL, Wheatley MD, Schlauch KA, Cramer GR, Cushman JC: Tissue-specific mRNA expression profiling in grape berry tissues. BMC Genomics. 2007, 8: 187. 10.1186/1471-2164-8-187.

    PubMed  PubMed Central  Article  Google Scholar 

  21. 21.

    Mathiason K, He D, Grimplet J, Venkateswari J, Galbraith DW, Or E, Fennell A: Transcript profiling in Vitis riparia during chilling requirement fulfillment reveals coordination of gene expression patterns with optimized bud break. Funct Integr Genomics. 2009, 9: 81-96. 10.1007/s10142-008-0090-y.

    PubMed  CAS  Article  Google Scholar 

  22. 22.

    Sweetman C, Deluc LG, Cramer GR, Ford CM, Soole KL: Regulation of malate metabolism in grape berry and other developing fruits. Phytochemistry. 2009, 70: 1329-1344. 10.1016/j.phytochem.2009.08.006.

    PubMed  CAS  Article  Google Scholar 

  23. 23.

    Tattersall EAR, Grimplet J, DeLuc L, Wheatley MD, Vincent D, Osborne C, Ergul A, Lomen E, Blank RR, Schlauch KA, et al: Transcript abundance profiles reveal larger and more complex responses of grapevine to chilling compared to osmotic and salinity stress. Funct Integr Genomics. 2007, 7: 317-333. 10.1007/s10142-007-0051-x.

    PubMed  CAS  Article  Google Scholar 

  24. 24.

    Tillett RL, Wheatley MD, Tattersall EAR, Schlauch KA, Cramer GR, Cushman JC: The Vitis vinifera C-repeat binding protein 4 (VvCBF4) transcriptional factor enhances freezing tolerance in wine grape. Plant Biotechnol J. 2012, 10: 105-124. 10.1111/j.1467-7652.2011.00648.x.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  25. 25.

    Waters DLE, Holton TA, Ableff EM, Lee LS, Henry RJ: cDNA microarray analysis of developing grape (Vitis vinifera cv. Shiraz) berry skin. Funct Integr Genomics. 2005, 5: 40-58. 10.1007/s10142-004-0124-z.

    PubMed  CAS  Article  Google Scholar 

  26. 26.

    Waters DLE, Holton TA, Ablett EM, Lee LS, Henry RJ: The ripening wine grape berry skin transcriptome. Plant Sci. 2006, 171: 132-138. 10.1016/j.plantsci.2006.03.002.

    CAS  Article  Google Scholar 

  27. 27.

    Zamboni A, Di Carli M, Guzzo F, Stocchero M, Zenoni S, Ferrarini A, Tononi P, Toffali K, Desiderio A, Lilley KS, et al: Identification of putative stage-specific grapevine berry biomarkers and omics data integration into networks. Plant Physiol. 2010, 154: 1439-1459. 10.1104/pp.110.160275.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  28. 28.

    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  CAS  PubMed Central  Article  Google Scholar 

  29. 29.

    Perrone I, Pagliarani C, Lovisolo C, Chitarra W, Roman F, Schubert A: Recovery from water stress affects grape leaf petiole transcriptome. Planta. 2012, 235 (6): 1383-1396. 10.1007/s00425-011-1581-y.

    PubMed  CAS  Article  Google Scholar 

  30. 30.

    Lim CJ, Yang KA, Hong JK, Choi AS, Yun DJ, Hong JC, Chung WS, Lee SY, Cho MJ, Lim CO: Gene expression profiles during heat acclimation in Arabidopsis thaliana suspension-culture cells. J Plant Res. 2006, 119: 373-383. 10.1007/s10265-006-0285-z.

    PubMed  CAS  Article  Google Scholar 

  31. 31.

    Yamakawa H, Hirose T, Kuroda M, Yamaguchi T: Comprehensive expression profiling of rice grain filling-related genes under high temperature using DNA microarray. Plant Physiol. 2007, 144: 258-277. 10.1104/pp.107.098665.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  32. 32.

    Hewezi T, Leger M, Gentzbittel L: A comprehensive analysis of the combined effects of high light and high temperature stresses on gene expression in sunflower. Ann Bot. 2008, 102: 127-140. 10.1093/aob/mcn071.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  33. 33.

    Frank G, Pressman E, Ophir R, Althan L, Shaked R, Freedman M, Shen S, Firon N: Transcriptional profiling of maturing tomato (Solanum lycopersicum L.) microspores reveals the involvement of heat shock proteins, ROS scavengers, hormones, and sugars in the heat stress response. J Exp Bot. 2009, 60: 3891-3908. 10.1093/jxb/erp234.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  34. 34.

    Wei K, Cheng F, Dong H, Zhang Q, Liu K, Cao Z: Microarray analysis of gene expression profile related to grain storage metabolism in rice endosperms as affected by high temperature at filling stage. Scientia Agricultura Sinica. 2010, 43: 1-11.

    CAS  Google Scholar 

  35. 35.

    Galmes J, Medrano H, Flexas J: Photosynthetic limitations in response to water stress and recovery in Mediterranean plants with different growth forms. New Phytol. 2007, 175: 81-93. 10.1111/j.1469-8137.2007.02087.x.

    PubMed  CAS  Article  Google Scholar 

  36. 36.

    Vinocur B, Altman A: Recent advances in engineering plant tolerance to abiotic stress: achievements and limitations. Curr Opin Biotechnol. 2005, 16: 123-132. 10.1016/j.copbio.2005.02.001.

    PubMed  CAS  Article  Google Scholar 

  37. 37.

    Gu RS, Fonseca S, Puskas LG, Hackler L, Zvara A, Dudits D, Pais MS: Transcript identification and profiling during salt stress and recovery of Populus euphratica. Tree Physiol. 2004, 24: 265-276. 10.1093/treephys/24.3.265.

    PubMed  CAS  Article  Google Scholar 

  38. 38.

    Bita CE, Zenoni S, Vriezen WH, Mariani C, Pezzotti M, Gerats T: Temperature stress differentially modulates transcription in meiotic anthers of heat-tolerant and heat-sensitive tomato plants. BMC Genomics. 2011, 12: 384. 10.1186/1471-2164-12-384.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  39. 39.

    Chauhan H, Khurana N, Agarwal P, Khurana P: Heat shock factors in rice (Oryza sativa L.): genome-wide expression analysis during reproductive development and abiotic stress. Mol Genet Genomics. 2011, 286: 171-187. 10.1007/s00438-011-0638-8.

    PubMed  CAS  Article  Google Scholar 

  40. 40.

    Ginzberg I, Barel G, Ophir R, Tzin E, Tanami Z, Muddarangappa T, de Jong W, Fogelman E: Transcriptomic profiling of heat-stress response in potato periderm. J Exp Bot. 2009, 60: 4411-4421. 10.1093/jxb/erp281.

    PubMed  CAS  Article  Google Scholar 

  41. 41.

    Mangelsen E, Kilian J, Harter K, Jansson C, Wanke D, Sundberg E: Transcriptome analysis of high-temperature stress in developing barley caryopses: early stress responses and effects on storage compound biosynthesis. Mol Plant. 2011, 4: 97-115. 10.1093/mp/ssq058.

    PubMed  CAS  Article  Google Scholar 

  42. 42.

    Huang B, Xu C: Identification and characterization of proteins associated with plant tolerance to heat stress. J Integr Plant Biol. 2008, 50: 1230-1237. 10.1111/j.1744-7909.2008.00735.x.

    PubMed  CAS  Article  Google Scholar 

  43. 43.

    Didomenico BJ, Bugaisky GE, Lindquist S: Heat -shock and recovery are mediated by different translational mechanisms. Proceedings of the National Academy of Sciences of the United States of America-Biological Sciences. 1982, 79: 6181-6185. 10.1073/pnas.79.20.6181.

    CAS  Article  Google Scholar 

  44. 44.

    Wang LJ, Fan L, Loescher W, Duan W, Liu GJ, Cheng J-S, Luo HB, Li SH: Salicylic acid alleviates decreases in photosynthesis under heat stress and accelerates recovery in grapevine leaves. BMC Plant Biol. 2010, 10: 34. 10.1186/1471-2229-10-34.

    PubMed  PubMed Central  Article  Google Scholar 

  45. 45.

    Wang W, Vinocur B, Shoseyov O, Altman A: Role of plant heat-shock proteins and molecular chaperones in the abiotic stress response. Trends Plant Sci. 2004, 9: 244-252. 10.1016/j.tplants.2004.03.006.

    PubMed  CAS  Article  Google Scholar 

  46. 46.

    Heckathorn SA, Downs CA, Sharkey TD, Coleman JS: The small, methionine-rich chloroplast heat-shock protein protects photosystem II electron transport during heat stress. Plant Physiol. 1998, 116: 439-444. 10.1104/pp.116.1.439.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  47. 47.

    Sanmiya K, Suzuki K, Egawa Y, Shono M: Mitochondrial small heat-shock protein enhances thermotolerance in tobacco plants. FEBS Lett. 2004, 557: 265-268. 10.1016/S0014-5793(03)01494-7.

    PubMed  CAS  Article  Google Scholar 

  48. 48.

    Myouga F, Motohashi R, Kuromori T, Nagata N, Shinozaki K: An Arabidopsis chloroplast-targeted Hsp101 homologue, APG6, has an essential role in chloroplast development as well as heat-stress response. Plant J. 2006, 48: 249-260. 10.1111/j.1365-313X.2006.02873.x.

    PubMed  CAS  Article  Google Scholar 

  49. 49.

    Lee U, Wie C, Escobar M, Williams B, Hong SW, Vierling E: Genetic analysis reveals domain interactions of Arabidopsis Hsp100/ClpB and cooperation with the small heat shock protein chaperone system. Plant Cell. 2005, 17: 559-571. 10.1105/tpc.104.027540.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  50. 50.

    Pressman E, Shaked R, Firon N: Tomato response to heat stress: focus on pollen grains. Plant Stress. 2007, 1: 216-227.

    Google Scholar 

  51. 51.

    Hartl FU: Molecular chaperones in cellular protein folding. Nature. 1996, 381: 571-580. 10.1038/381571a0.

    PubMed  CAS  Article  Google Scholar 

  52. 52.

    Shi WM, Muramoto Y, Ueda A, Takabe T: Cloning of peroxisomal ascorbate peroxidase gene from barley and enhanced thermotolerance by overexpressing in Arabidopsis thaliana. Gene. 2001, 273: 23-27. 10.1016/S0378-1119(01)00566-2.

    PubMed  CAS  Article  Google Scholar 

  53. 53.

    Schramm F, Ganguli A, Kiehlmann E, Englich G, Walch D, von Koskull-Doring P: The heat stress transcription factor HsfA2 serves as a regulatory amplifier of a subset of genes in the heat stress response in Arabidopsis. Plant Mol Biol. 2006, 60: 759-772. 10.1007/s11103-005-5750-x.

    PubMed  CAS  Article  Google Scholar 

  54. 54.

    Rossel JB, Walter PB, Hendrickson L, Chow WS, Poole A, Mullineaux PM, Pogson BJ: A mutation affecting ascorbate peroxidase 2 gene expression reveals a link between responses to high light and drought tolerance. Plant Cell Environ. 2006, 29: 269-281. 10.1111/j.1365-3040.2005.01419.x.

    PubMed  CAS  Article  Google Scholar 

  55. 55.

    Kaplan F, Kopka J, Haskell DW, Zhao W, Schiller KC, Gatzke N, Sung DY, Guy CL: Exploring the temperature-stress metabolome of Arabidopsis. Plant Physiol. 2004, 136: 4159-4168. 10.1104/pp.104.052142.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  56. 56.

    Zu YG, Tang ZH, Yu JH, Liu SG, Wang W, Guo XR: Different responses of camptothecin and 10-hydroxycamptothecin to heat shock in Camptotheca acuminata seedlings. Acta Botanica Sinica. 2003, 45: 809-814.

    CAS  Google Scholar 

  57. 57.

    Busch W, Wunderlich M, Schoffl F: Identification of novel heat shock factor-dependent genes and biochemical pathways in Arabidopsis thaliana. Plant J. 2005, 41: 1-14.

    PubMed  CAS  Article  Google Scholar 

  58. 58.

    Nishizawa A, Yabuta Y, Shigeoka S: Galactinol and raffinose constitute a novel function to protect plants from oxidative damage. Plant Physiol. 2008, 147: 1251-1263. 10.1104/pp.108.122465.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  59. 59.

    Hincha DK, Zuther E, Heyer AG: The preservation of liposomes by raffinose family oligosaccharides during drying is mediated by effects on fusion and lipid phase transitions. Biochimica Et Biophysica Acta-Biomembranes. 2003, 1612: 172-177. 10.1016/S0005-2736(03)00116-0.

    CAS  Article  Google Scholar 

  60. 60.

    Nishizawa A, Yabuta Y, Yoshida E, Maruta T, Yoshimura K, Shigeoka S: Arabidopsis heat shock transcription factor A2 as a key regulator in response to several types of environmental stress. Plant J. 2006, 48: 535-547. 10.1111/j.1365-313X.2006.02889.x.

    PubMed  CAS  Article  Google Scholar 

  61. 61.

    Weston DJ, Karve AA, Gunter LE, Jawdy SS, Yang X, Allen SM, Wullschleger SD: Comparative physiology and transcriptional networks underlying the heat shock response in Populus trichocarpa, Arabidopsis thaliana and Glycine max. Plant Cell Environ. 2011, 34: 1488-1506. 10.1111/j.1365-3040.2011.02347.x.

    PubMed  CAS  Article  Google Scholar 

  62. 62.

    Kotak S, Larkindale J, Lee U, von Koskull-Doering 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.

    PubMed  CAS  Article  Google Scholar 

  63. 63.

    Baniwal SK, Bharti K, Chan KY, Fauth M, Ganguli A, Kotak S, Mishra SK, Nover L, Port M, Scharf KD, et al: Heat stress response in plants: a complex game with chaperones and more than twenty heat stress transcription factors. J Biosci. 2004, 29: 471-487. 10.1007/BF02712120.

    PubMed  CAS  Article  Google Scholar 

  64. 64.

    Klimecka M, Muszynska G: Structure and functions of plant calcium-dependent protein kinases. Acta Biochim Pol. 2007, 54: 219-233.

    PubMed  CAS  Google Scholar 

  65. 65.

    Lohmann C, Eggers-Schumacher G, Wunderlich M, Schoffl F: Two different heat shock transcription factors regulate immediate early expression of stress genes in Arabidopsis. Mol Genet Genomics. 2004, 271: 11-21. 10.1007/s00438-003-0954-8.

    PubMed  CAS  Article  Google Scholar 

  66. 66.

    Schulz-Raffelt M, Lodha M, Schroda M: Heat shock factor 1 is a key regulator of the stress response in Chlamydomonas. Plant J. 2007, 52: 286-295. 10.1111/j.1365-313X.2007.03228.x.

    PubMed  CAS  Article  Google Scholar 

  67. 67.

    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  PubMed Central  Article  Google Scholar 

  68. 68.

    Suzuki N, Bajad S, Shuman J, Shulaev V, Mittler R: The transcriptional co-activator MBF1c is a key regulator of thermotolerance in Arabidopsis thaliana. J Biol Chem. 2008, 283: 9269-9275. 10.1074/jbc.M709187200.

    PubMed  CAS  Article  Google Scholar 

  69. 69.

    Jakoby M, Weisshaar B, Droge-Laser W: bZIP transcription factors in Arabidopsis. Trends Plant Sci. 2002, 7: 106-111. 10.1016/S1360-1385(01)02223-3.

    PubMed  CAS  Article  Google Scholar 

  70. 70.

    Qiu Y, Yu D: Over-expression of the stress-induced OsWRKY45 enhances disease resistance and drought tolerance in Arabidopsis. Environmental and Experimental Botany. 2009, 65: 35-47. 10.1016/j.envexpbot.2008.07.002.

    CAS  Article  Google Scholar 

  71. 71.

    Olsen AN, Ernst HA, Lo Leggio L, Skriver K: NAC transcription factors: structurally distinct, functionally diverse. Trends Plant Sci. 2005, 10: 79-87. 10.1016/j.tplants.2004.12.010.

    PubMed  CAS  Article  Google Scholar 

  72. 72.

    Tran LSP, Nakashima K, Sakuma Y, Simpson SD, Fujita Y, Maruyama K, Fujita M, Seki M, Shinozaki K, Yamaguchi-Shinozaki K: Isolation and functional analysis of Arabidopsis stress-inducible NAC transcription factors that bind to a drought-responsive cis-element in the early responsive to dehydration stress 1 promoter. Plant Cell. 2004, 16: 2481-2498. 10.1105/tpc.104.022699.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  73. 73.

    Tarutani Y, Morimoto T, Sasaki A, Yasuda M, Nakashita H, Yoshida S, Yamaguchi I, Suzuki Y: Molecular characterization of two highly homologous receptor-like kinase genes, RLK902 and RKL1, in Arabidopsis thaliana. Biosci Biotechnol Biochem. 2004, 68: 1935-1941. 10.1271/bbb.68.1935.

    PubMed  CAS  Article  Google Scholar 

  74. 74.

    Osakabe Y, Maruyama K, Seki M, Satou M, Shinozaki K, Yamaguchi-Shinozaki K: Leucine-rich repeat receptor-like kinase1 is a key membrane-bound regulator of abscisic acid early signaling in Arabidopsis. Plant Cell. 2005, 17: 1105-1119. 10.1105/tpc.104.027474.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  75. 75.

    Zhou R, Li B, Liu H, Sun D: Progress in the participation of Ca(2+)-calmodulin in heat shock signal transduction. Progress in Natural Sci. 2009, 19: 1201-1208. 10.1016/j.pnsc.2008.12.011.

    CAS  Article  Google Scholar 

  76. 76.

    Mittler R, Finka A, Goloubinoff P: How do plants feel the heat?. Trends Biochem Sci. 2012, 37: 118-125. 10.1016/j.tibs.2011.11.007.

    PubMed  CAS  Article  Google Scholar 

  77. 77.

    Liu HT, Li B, Shang ZL, Li XZ, Mu RL, Sun DY, Zhou RG: Calmodulin is involved in heat shock signal transduction in wheat. Plant Physiol. 2003, 132: 1186-1195. 10.1104/pp.102.018564.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  78. 78.

    Saidi Y, Finka A, Goloubinoff P: Heat perception and signalling in plants: a tortuous path to thermotolerance. New Phytol. 2011, 190: 556-565. 10.1111/j.1469-8137.2010.03571.x.

    PubMed  CAS  Article  Google Scholar 

  79. 79.

    Kawasaki S, Borchert C, Deyholos M, Wang H, Brazille S, Kawai K, Galbraith D, Bohnert HJ: Gene expression profiles during the initial phase of salt stress in rice. Plant Cell. 2001, 13: 889-905.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  80. 80.

    Li C, Wong WH: Model-based analysis of oligonucleotide arrays: Expression index computation and outlier detection. Proc Natl Acad Sci USA. 2001, 98: 31-36. 10.1073/pnas.98.1.31.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  81. 81.

    Li C, Hung Wong W: Model-based analysis of oligonucleotide arrays: model validation, design issues and standard error application. Genome Biol. 2001, 2: RESEARCH0032.

    PubMed  CAS  PubMed Central  Google Scholar 

  82. 82.

    Deluc LG, Quilici DR, Decendit A, Grimplet J, Wheatley MD, Schlauch KA, Merillon JM, Cushman JC, Cramer GR: Water deficit alters differentially metabolic pathways affecting important flavor and quality traits in grape berries of Cabernet Sauvignon and Chardonnay. BMC Genomics. 2009, 10: 212. 10.1186/1471-2164-10-212.

    PubMed  PubMed Central  Article  Google Scholar 

  83. 83.

    Reid KE, Olsson N, Schlosser J, Peng F, Lund ST: An optimized grapevine RNA isolation procedure and statistical determination of reference genes for real-time RT-PCR during berry development. BMC Plant Biol. 2006, 6: 27. 10.1186/1471-2229-6-27.

    PubMed  PubMed Central  Article  Google Scholar 

  84. 84.

    Jarosova J, Kundu JK: Validation of reference genes as internal control for studying viral infections in cereals by quantitative real-time RT-PCR. BMC Plant Biol. 2010, 10: 146. 10.1186/1471-2229-10-146.

    PubMed  PubMed Central  Article  Google Scholar 

Download references


This work was supported by the National Natural Science Foundation of China (No. 30771758).

Author information



Corresponding authors

Correspondence to Li-Jun Wang or Shao-Hua Li.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

GTL, LJW, SHL and JFW carried out the genomic experiments, or conducted the study design and data analysis. WD, BHW, HGX, PGF participated in the study design and in data analysis. LJW, GTL, SHL, GC and ZWD drafted the manuscript. All authors read, revised and approved the final manuscript.

Electronic supplementary material

Additional file 1: Gene-specific primers for qRT-PCR.(DOC 48 KB)

Linear correlation analysis (r = 0.982) between qRT-PCR and microarray results for 12 genes.

Additional file 2: X; log2 fold change value from microarray data; Y: log2 fold change value from qRT-PCR data. (DOCX 41 KB)


Additional file 3: Genes upregulated during heat stress (HS) and downregulated after the subsequent recovery (RC) in grape leaves.(DOCX 23 KB)


Additional file 4: Genes upregulated during heat stress (HS) and after the subsequent recovery (RC) in grape leaves.(DOCX 19 KB)


Additional file 5: Genes downregulated during heat stress (HS) and upregulated after the subsequent recovery (RC) in grape leaves.(DOCX 20 KB)


Additional file 6: Genes downregulated during heat stress and after the subsequent recovery in grape leaves.(DOCX 24 KB)

Additional file 7: Genes upregulated or downregulated unique to heat stress in grape leaves.(DOCX 69 KB)

Additional file 8: Genes upregulated or downregulated unique to the recovery in grapevine leaves.(DOCX 40 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

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 (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and Permissions

About this article

Cite this article

Liu, GT., Wang, JF., Cramer, G. et al. Transcriptomic analysis of grape (Vitis viniferaL.) leaves during and after recovery from heat stress. BMC Plant Biol 12, 174 (2012).

Download citation


  • Heat Stress
  • Heat Tolerance
  • Heat Shock Transcription Factor
  • Subsequent Recovery
  • WRKY Transcription Factor