- Research article
- Open Access
Transcriptomic network analyses of leaf dehydration responses identify highly connected ABA and ethylene signaling hubs in three grapevine species differing in drought tolerance
BMC Plant Biologyvolume 16, Article number: 118 (2016)
Grapevine is a major food crop that is affected by global climate change. Consistent with field studies, dehydration assays of grapevine leaves can reveal valuable information of the plant’s response at physiological, transcript, and protein levels. There are well-known differences in grapevine rootstocks responses to dehydration. We used time-series transcriptomic approaches combined with network analyses to elucidate and identify important physiological processes and network hubs that responded to dehydration in three different grapevine species differing in their drought tolerance.
Transcriptomic analyses of the leaves of Cabernet Sauvignon, Riparia Gloire, and Ramsey were evaluated at different times during a 24-h controlled dehydration. Analysis of variance (ANOVA) revealed that approximately 11,000 transcripts changed significantly with respect to the genotype x treatment interaction term and approximately 6000 transcripts changed significantly according to the genotype x treatment x time interaction term indicating massive differential changes in gene expression over time. Standard analyses determined substantial effects on the transcript abundance of genes involved in the metabolism and signaling of two known plant stress hormones, abscisic acid (ABA) and ethylene. ABA and ethylene signaling maps were constructed and revealed specific changes in transcript abundance that were associated with the known drought tolerance of the genotypes including genes such as VviABI5, VviABF2, VviACS2, and VviWRKY22. Weighted-gene coexpression network analysis (WGCNA) confirmed these results. In particular, WGCNA identified 30 different modules, some of which had highly enriched gene ontology (GO) categories for photosynthesis, phenylpropanoid metabolism, ABA and ethylene signaling. The ABA signaling transcription factors, VviABI5 and VviABF2, were highly connected hubs in two modules, one being enriched in gaseous transport and the other in ethylene signaling. VviABI5 was distinctly correlated with an early response and high expression for the drought tolerant Ramsey and with little response from the drought sensitive Riparia Gloire. These ABA signaling transcription factors were highly connected to VviSnRK1 and other gene hubs associated with sugar, ethylene and ABA signaling.
A leaf dehydration assay provided transcriptomic evidence for differential leaf responses to dehydration between genotypes differing in their drought tolerance. WGCNA proved to be a powerful network analysis approach; it identified 30 distinct modules (networks) with highly enriched GO categories and enabled the identification of gene hubs in these modules. Some of these genes were highly connected hubs in both the ABA and ethylene signaling pathways, supporting the hypothesis that there is substantial crosstalk between the two hormone pathways. This study identifies solid gene candidates for future investigations of drought tolerance in grapevine.
Boyer  reviewed the impact of the environment on crop production and highlighted the need for crops better suited to these environments. Much has been learned since then, yet our understanding of plant responses to abiotic and biotic stresses is very incomplete. Drought signaling within plants is a complex process involving many different signaling cascades . A rapid assay was developed to assess the physiological response of different grapevine genotypes to dehydration . This assay is a simple approach that can determine differences in dehydration sensitivity at the physiological and Omic levels.
Climate change is expected to affect water and land availability [4, 5]. Rootstocks are used in viticulture because they can confer pest or drought resistance, alter vigor to the scion or the fruit-bearing portion of the plant, thus, impacting fruit quality; rootstocks are vital in most viticultural regions [6, 7]. Much research has focused on the scion-rootstock relationship [8–11], but there is little research on the rootstock response to abiotic conditions.
Three different Vitis genotypes were shown previously to have differences in their dehydration sensitivity [3, 6, 12]. Cabernet Sauvignon (Vitis vinifera L.) along with two North American Vitis species commonly used as rootstocks, Ramsey (Vitis champinii Planch., a naturally occurring hybrid between Vitis candicans Engelm. and Vitis rupestris Scheele) and Riparia Gloire (Vitis riparia Michx.). Ramsey, a drought tolerant genotype, originates from hot, dry regions of Texas. Riparia Gloire originates from wet, riparian areas and is drought sensitive [3, 6, 12].
Transcriptomic analyses allow one to have a “holistic” snapshot of the plant’s transcriptional response to a changing environment . A time-series transcriptomic analysis allows one to begin to elucidate the sensitivities of the response and the primary or secondary responses. Co-expression analyses allow one to identify networks, genes that have high connectivity or correlation with each other. A particularly powerful approach is the weighted coexpression network analysis (WGCNA) method [13–16]. This analysis can identify clusters (modules) of genes with high biological meaning. It can also identify those genes with high connectivity or module membership, which are essentially hub genes. Similar to airport hubs, if one hub is not functioning, the whole system can slow down or become chaotic. Therefore we consider hub genes important or essential genes for the proper functioning of the system or organism.
Abscisic acid (ABA) and ethylene are two important hormones that regulate abiotic stress responses in plants [2, 17]. In a preliminary survey of more than 30 genotypes, we found variation in the increase in transcript abundance of VviNCED3, the rate-limiting step in ABA biosynthesis, in response to rapid dehydration. We hypothesize that the transcriptomic responses to rapid dehydration between these grapevine genotypes are different and may involve ABA signaling. In this time-series transcriptomic study, we identify significant genes by standard and network (WGCNA) analysis methods; a number of these genes are involved in ABA and ethylene signaling and correlate with the relative drought tolerance between the genotypes. In particular, ABI5, a transcription factor gene normally associated with ABA regulation of germination, was highly sensitive and increased in Ramsey leaves, the most drought tolerant of the three genotypes, early in the response to dehydration, but there was little effect of dehydration on abscisic acid insensitive 5 (ABI5) transcription factor expression in Riparia Gloire leaves, the most drought sensitive of the genotypes. This gene and other hub genes are identified as “solid” candidates for future drought tolerance research.
Plant material and experimental conditions
Three grapevine genotypes (Vitis vinifera cv. Cabernet Sauvignon clone 8, Vitis riparia cv. Riparia Gloire, and Vitis champinii cv. Ramsey (a naturally occurring hybrid between Vitis candicans and Vitis rupestris)) were pruned to two shoots and grown in 13.3 liter pots containing 10 L SuperSoil® potting mix supplemented with slow release fertilizer (5-10-10). The original vine cuttings of Riparia Gloire and Ramsey were obtained from Dr. Andrew Walker at the University of California, Davis, CA, USA. The original vine cuttings of Cabernet Sauvignon vines were obtained from Inland Desert Nursery in Benton City, WA, USA. Vines were grown in a greenhouse with supplemental sodium vapor lamp lighting (day/night cycles of 16 h/8 h light (minimum 400 μE m−2 s−1) and 28 °C/18 °C). Fully-developed leaves were subjected to dehydration as previously described . Leaves were removed from dehydration boxes at 1, 2, 4, 8, and 24 h and frozen immediately in liquid nitrogen. Control leaves were taken from the second shoot of the intact plant at the corresponding daytime of the dehydration assay to account for circadian effects on transcript abundance.
RNA extraction and microarray hybridization
Frozen leaves were ground using a Retsch MM 301 ball mill  for 1 min at 30 revolutions s−1. Total RNA was extracted from approximately 100 mg of tissue using a cetyl trimethylammonium bromide (CTAB)-based method [19, 20]. Extracts were treated with DNase (Qiagen RNeasy Plant Kit, ) according to manufacturer’s instructions. RNA quality and quantity were assessed with a Nanodrop ND-1000 spectrophotometer (ThermoFisher Scientific, Waltham, MA) and an Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA) according to the manufacturer’s instructions. Microarrays were hybridized by MOgene (St. Louis, MO, USA) using the NimbleGen microarray 090818 Vitis exp HX12 (Roche, NimbleGen Inc., Madison, WI, USA) according to the manufacturer’s instructions.
All microarrays were analyzed as one set as previously described [22–24]. One of ninety arrays (Rip4S2) exhibited considerable spatial variation, and thus was excluded. Two other arrays (Cab4S3 and Cab8S1) had notable statistically significant outliers across replicates when compared to other arrays, and were thus also excluded.
A simple 3-way analysis of variance (ANOVA) was performed on normalized (log-transformed) and quality-controlled processed data to determine which probesets on the array were differentially expressed with statistical significance across Genotype, Treatment, and Time, and the 2-way and 3-way interaction of these effects. The processed and normalized expression values were not normally distributed, thus an extension of the Kruskal Wallis rank sum test was used for the ANOVA . A multiple testing correction was applied to the p-values of the ANOVA , and any probeset with a significant Genotype x Treatment or Genotype x Treatment x Time interaction term with adjusted p-value ≤ 0.05 was examined further.
Principal component analysis (PCA) was applied to quality-controlled expression data using the covariance matrix to visualize any trends in the expression data [27–30]. The PCA (Additional file 1) showed a very clear separation between genotypes, supporting the large number of probesets with a statistically significant tissue effect (95.6 %).
Gene expression was also evaluated with WGCNA  using the following settings for the adjacency function (datExpr, power = 16, type = "signed hybrid", corFnc = "bicor", corOptions = "use = 'p', maxPOutliers = 0.1") and for the cuttreeDynamic function (dendro = geneTree, distM = dissTOM, method = "hybrid", deepSplit = 2, pamRespectsDendro = F, minClusterSize = 30); these functions have been shown to be the best approach for biologically meaningful results . WGCNA also confirmed clear separation by genotype and treatment (Additional file 2).
Functional categorization of significant transcripts was performed with the BiNGO plugin  in Cytoscape  using a gene ontology (GO) file created with the EnsemblPlants BioMart  for Vitis vinifera. Overrepresented (enriched) categories were determined using a hypergeometric test with a significance threshold at 0.05 after a Benjamini and Hochberg false discovery rate correction.
Dehydration causes massive changes in gene expression
Fully mature Cabernet Sauvignon, Riparia Gloire, and Ramsey leaves were dehydrated  to asses rapid transcriptomic changes. Briefly, the leaf was excised from the plant and placed into a dehydration box for various time points over a 24-h period. Leaf dehydration occurred in the air above a solution of NaCl in a sealed container in a growth chamber. Leaves were removed from the box and immediately frozen in liquid nitrogen at specific time points of dehydration. Control samples were taken from the same plant at the corresponding time to account for any circadian effect on transcript abundance. RNA was extracted from three experimental replicates for treatment and control samples at each time point.
In order to test the hypothesis that the response to rapid dehydration between the different genotypes is different over time at the transcript level, a 2 × 2 × 5 factorial (Genotype x Treatment x Time) experimental design was established. Transcriptomic analysis was carried out using the NimbleGen Grape Whole-Genome Microarray. A parametric ANOVA was originally performed but because the expression data were not quite normally distributed, the expression data were reanalyzed using a nonparametric ANOVA. The nonparametric ANOVA reduced the overall number of genes with differential expression by about 10 %. The nonparametric ANOVA determined that the abundance of 28,030 transcripts changed significantly with an adjusted p-value ≤ 0.05 (herein referred to as “significant” throughout the paper) with respect to genotype (Table 1, Additional file 3). Clearly there are large differences in gene expression between these species regardless of treatment. There were 18,237 and 23,656 transcripts that changed significantly with Treatment and Time, respectively; 11,436 transcripts changed significantly with respect to Genotype x Treatment; 24,543 transcripts changed significantly with genotype over time; 17,488 transcripts changed significantly with Treatment x Time, and finally 6,285 transcripts changed significantly for the Genotype x Treatment x Time term (Table 1, Additional file 3).
Two approaches were taken to analyze these large datasets: “standard analysis” and “network analysis”. With standard analysis we used a posteriori knowledge to sort through known biochemical and signaling pathways affected by dehydration. With network analysis, we took an a priori approach by using WGCNA and GO enrichment methods.
As we are interested in elucidating mechanisms of drought tolerance, we focused on the Genotype x Treatment and the Genotype x Treatment x Time interaction sets of genes. GO categories for these gene sets were determined with a custom Vitis GO file (see Methods) and analyzed for significant overrepresentation using BiNGO, a Cytoscape plug-in . However, GO enrichment analysis was not very informative with these large datasets. There were two biological processes significantly overrepresented for the Genotype x Treatment gene set (Additional file 4): translation and phenylpropanoid metabolism (more will be discussed about these categories later in the network analysis section). The Genotype x Treatment x Time gene set had one category that was significantly overrepresented: oxidation-reduction (Additional file 5). This latter interaction term did not provide any obvious clues about differences between the genotypes, so we used previous knowledge from our research to ascertain if there were differences in ABA metabolism and signaling, the hypothesis in which we were most interested.
Dehydration induces significant changes in ABA metabolism transcripts between the genotypes
The rate-limiting step in ABA biosynthesis is catalyzed by genes that encode 9-cis-epoxycarotenoid dioxygenase (NCED) [34, 35]. In Vitis there are three NCED genes, which can lead to the production of ABA . The gene symbols used are based upon the symbol used to the closest ortholog in Arabidopsis. Vitis and Arabidopsis loci for these symbols are listed in Additional file 6. A significant difference in the transcript abundance for these three genes was observed (Fig. 1). In Cabernet Sauvignon, VviNCED3 expression was slightly decreased at 1 h of dehydration while the expression in Riparia Gloire and Ramsey was increased in response to dehydration. Riparia Gloire had larger initial response than Ramsey, but after 4 h, VviNCED3 expression of Ramsey and Cabernet Sauvignon exceeded that of Riparia Gloire. These results confirmed our preliminary results that there were differences between the genotypes in the expression of VviNCED3 in response to dehydration.
VviNCED5 transcript abundance increased in Cabernet Sauvignon leaves from 2 to 8 h; there was an increase in Riparia Gloire at the 8 h time point and a slight increase in Ramsey (Fig. 1). The transcript abundance of VviNCED6 in both North American genotypes increased within 1 h of dehydration, but there was little response in Cabernet Sauvignon (Fig. 1). These data indicate a different regulation of these genes in response to dehydration as well as differences between the genotypes.
ABA action within the plant is also dependent on degradation, conjugation and transport. Significant differences in the transcript abundance of genes involved in these processes were different between the genotypes. Degradation of ABA is catalyzed by a group of cytochrome P450 enzymes known as ABA-hydroxylases, and then continues by a few non-enzymatic steps leading to the formation of phaseic acid. Two genes annotated to be ABA-hydroxylases, VviCYP707A3 and VviCYP707A4, responded differently to the dehydration (Fig. 1). The transcript abundance of VviCYP707A4 decreased throughout the experiment for all genotypes. At 4 h of dehydration, transcript abundance in both Riparia Gloire and Ramsey were at their lowest points with a log2 fold decrease greater than 4. Interestingly, VviCYP707A3 increased in expression throughout the experiment, most notably in Cabernet Sauvignon (Fig. 1).
Active ABA can also be produced through β-glucosidase, which involves the hydrolysis of an inactive form of ABA, Glc-conjugated ABA (ABA-GE), to active ABA. These enzymes are localized in the vacuole where ABA-GE is known to be stored . Previously, Zhang et al.  found significant differences in the expression of three genes in Vitis encoding β-glucosidases in Vitis vinifera cv. Muscat Hamburg berry ripening. In our study, transcript abundance of VviBGLU40 increased in Riparia Gloire leaves, most notably at 2 h of dehydration (Fig. 1). Interestingly, Zhang et al.  reported that VviBG3 transcript abundance decreased through véraison in berry samples. A similar expression profile was observed in this study in all genotypes surveyed, most notably in Riparia Gloire at 8 h of dehydration (Fig. 1).
ABA transport can also affect ABA concentrations. The transcript abundance of two ABA transporters was shown to increase significantly in response to dehydration (Fig. 1) and was different between the genotypes over time. VviABCG25 is an ATP-binding cassette (ABC) transporter that exports ABA from vascular tissues allowing ABA to reach distant guard cells [39, 40]. After 4 h of dehydration, transcript abundance increased with a log2 fold change of more than 2. All genotypes displayed a large increase in transcript abundance by 8 h of dehydration.
In Arabidopsis, AtABCG40 imports ABA directly into guard cells . In our study, transcript abundance in all genotypes increased within 1 h of dehydration (Fig. 1). Interestingly, Riparia Gloire increased nearly 5-fold indicating a massive change in transcript abundance. This gene may contribute to the dehydration and ABA sensitivity of stomatal conductance of Riparia Gloire leaves .
Together these results indicate that transcripts involved in ABA metabolism changed significantly in response to rapid dehydration and the responses between the genotypes were different, consistent with the differences in dehydration sensitivity previously observed .
ABA core-signaling response to dehydration
Downstream of ABA biosynthesis is a complex ABA signaling network involving many different genes. Recently, Lumba et al.  took a systems biology approach to create an ABA core-signaling network consisting of over 500 interactions between 138 proteins in Arabidopsis. Many different processes are represented such as proteins involved in transport, metabolism, proteolysis, calcium sensing, as well as numerous transcription factors and kinases. Vitis orthologs were compiled based on the closest orthologs identified by Gramene ( release 44 (January 2015); see list in Additional file 6. Significant differences in gene expression within the ABA core-signaling network were detected (Fig. 2).
ABA binds to receptors in the cytoplasm known as PYR/PYL/RCAR proteins [44, 45]. Evidence indicates that there are additional receptors located at the plasma membrane . Interestingly, the transcript abundance of VviPYL4 changed significantly in response to dehydration between the genotypes; transcript abundance was decreased in Cabernet Sauvignon leaves while there was a slight increase in Riparia Gloire and little response in Ramsey (Fig. 2). ABA receptors interact with specific type 2C protein phosphatases (PP2C), which inhibit the activity of serine/threonine-protein kinase 2.6 (SnRK2.6; OST1) when ABA is not present. In the presence of ABA, PP2C proteins dissociate from the kinase allowing autophosphorylation and subsequent downstream signaling. The transcript abundance of two PP2C genes, VviHAI1 and VviAHG3, increased significantly in response to dehydration between the genotypes (Fig. 2). The transcript abundance of VviOST1 significantly increased in response to dehydration, most notably in Ramsey after 24 h of dehydration (Fig. 2). Downstream targets of VviOST1 also increased significantly in response to dehydration (Fig. 2). Targets include an ABA-responsive binding elements factor (ABF) that contains an ABA-responsive element (ABRE; PyACGTGG/TC) as a conserved cis-element in the promoter region . In our study, VviABF2 transcript abundance increased significantly in response to dehydration in all the genotypes (Fig. 2). Another ABF protein that is a target of OST1 is ABI5. VviABI5 transcript abundance increased significantly in response to dehydration, most notably in Ramsey leaves as early as 1 h of dehydration (Fig. 2). Recently, Yoshida et al.  noted that ABF2 in Arabidopsis is one of four predominant AREB/ABF transcription factors downstream of SnRK2.6 (OST1) in response to various osmotic stress conditions. Interestingly, the authors did not see a significant increase in the expression of AtABI5 under their conditions. In our study, there is clear evidence for the induction and differential expression of VviABI5 in response to dehydration.
A number of genes are induced by ABA, but lack the specific binding element mentioned above including proteins known to be involved in transport. The transcript abundance of two genes, guard cell S-type anion channel, VviSLAC1, and nitrate transporter, VviNRT1, changed significantly in response to dehydration (Fig. 2). SLAC1 is required for stomatal closure under conditions of high CO2 and ABA . Interestingly, VviSLAC1 transcript abundance decreased in response to dehydration (Fig. 2). In contrast, the transcript abundance of VviNRT1 increased significantly in response to dehydration most notably in Cabernet Sauvignon and Ramsey leaves (Fig. 2). NRT1 acts as both a low- and high-affinity nitrate transporter in Arabidopsis depending on it phosphorylation status . Nitrate concentrations appear to be important because nitric oxide (NO) is an important signaling molecule in ABA-induced stomatal closure and the production of NO is mediated through nitrate reductase activity .
Two additional transcripts shown to be within the ABA core-signaling network are the ABA-inducible HVA22-like homolog D (HVA22D) and mildew resistance locus 4 (MLO4). In our study, VviHVA22D transcript abundance increased significantly, most notably in Ramsey at 24 h of dehydration (Fig. 2). In contrast, VviMLO4 transcript abundance significantly decreased in response to dehydration, again most notably in Ramsey starting at 4 h of dehydration and continuing throughout the experiment (Fig. 2).
Metabolic transcripts within the ABA core-signaling network
Multiple transcripts within the ABA core-signaling network are known to be involved in metabolism. Changes in the expression of a number of these genes were observed (Fig. 2). Transcript abundance changes occurred more noticeably in Ramsey. For example, VviGPX3 encodes a glutathione peroxidase; its transcript abundance significantly decreased more than 1.5 log2 fold at 24 h of dehydration in Ramsey (Fig. 2). In addition, VviBCAT2 (branched-chain amino acid transaminase 2), VviGAT1 (involved in amino acid transport), and VviCAD5 (cinnamyl alcohol dehydrogenase) significantly increased in expression in response to dehydration (Additional files 3 and 6). The transcript abundance of an alcohol dehydrogenase, VviADH1, increased in response to dehydration within 4 h in Cabernet Sauvignon leaves with both Riparia Gloire and Ramsey increasing later in the experiment (Fig. 2).
Kinase transcripts within the ABA core-signaling network
In addition to VviOST1, a number of other transcripts encoding kinases significantly changed in response to dehydration. For example, a number of SnRK3 kinases significantly increased in expression in response to dehydration (Fig. 2). These kinases are involved in a number of plant stress responses including cold, salt, and drought . For example, the transcript abundance of VviSnRK3.11 increased most notably in Riparia Gloire at 2 h of dehydration (Fig. 2). Another, VviSnRk3.16 increased with a peak in expression at 8 h in both Cabernet Sauvignon and Ramsey with little response in Riparia Gloire. Finally, the transcript abundance of VviSnRK3.6 increased in all genotypes surveyed, the earliest in Riparia Gloire leaves at 4 h of dehydration (Fig. 2).
Another kinase of note changing in response to dehydration is a diacylglycerol kinase (DGK), VviDGK2 (Fig. 2). DGK synthesizes phosphatidic acid (PA), which is an important lipid-signaling molecule in plants involved in both biotic and abiotic signaling pathways . The transcript abundance of VviDGK2 increased as early as 1 h of dehydration in all genotypes (Fig. 2). The largest increase was observed in Cabernet Sauvignon leaves with a log2 fold increase of nearly 2.5, with Riparia Gloire lower at 1.8, and Ramsey at 1.2.
Other kinases increasing in expression include VviETR2 and VviMAP3Kδ4. Ethylene receptor 2 (ETR2) is a member of a group of ethylene receptors, which upon binding ethylene initiate a large signaling cascade (see below). In response to dehydration, VviETR2 increased in all genotypes with Cabernet Sauvignon and Ramsey increasing within 2 h of dehydration followed by Riparia Gloire at 4 h (Fig. 2). In contrast, VviMAP3Kδ4 encodes an activated mitogen kinase, which increased in Riparia Gloire most notably at 2 h of dehydration followed by both Cabernet Sauvignon and Ramsey.
Two examples of kinases that significantly decreased in transcript abundance were a CRINKLY4 related 2 (VviCCR2) and VviWNK2 (with no lysine (K)) (Fig. 2); VviCCR2 decreased with a log2 fold change of −1.8 and −1.6 in Cabernet Sauvignon and Ramsey, respectively, with little response in Riparia Gloire. In contrast, VviWNK2 displayed a log2 fold decrease of more than 2 in Riparia Gloire at 24 h of dehydration, the lowest observed (Additional file 3). Together, these data indicate significant changes in the expression of a number of kinases involved in multiple processes.
Transcription factor transcripts within the ABA core-signaling network
There are many transcription factors in the ABA core-signaling network. An investigation into the expression of all is outside the scope of this work. Instead, interesting differences between the genotypes are highlighted. Transcription factors from multiple families are represented indicating changes in many different signaling cascades. For example, a number of transcripts from the MYB (myeloblastosis), NAC (for NAM (no apical meristem), ATAF (Arabidopsis transcription activation factor), CUC (cup-shaped cotyledon)), and AP2/ERF (APETALA2/Ethylene-Responsive Element Binding Protein) domain transcription factor families changed significantly between the genotypes in response to dehydration.
Two MYB genes in particular within the ABA core-signaling network responded differently to dehydration. The transcript abundance of VviMYB12 decreased particularly in Riparia Gloire beginning at 1 h of dehydration (Fig. 2). In Arabidopsis, MYB12 regulates flavonoid biosynthetic genes  and subsequent reactive oxygen species (ROS) scavenging leading to greater drought tolerance . Conversely, a MYB77-like gene, VviMYB77 increased in expression particularly in Cabernet Sauvignon at 4 h of dehydration. MYB77 responds to ethylene and is involved in stress memory .
The transcript abundance of another two genes classified as homeodomain leucine zipper class I transcription factors, VviHB6 and VviHB12, changed significantly between the genotypes in response to dehydration. VviBH6 transcript abundance decreased in expression, particularly in Ramsey at 24 h of dehydration (Fig. 2). In contrast, VviBH12, followed a similar pattern in all genotypes displaying a decrease in expression early followed by an increase throughout the rest of the experiment.
One of the largest groups of transcription factors in plants is the AP2/ERF superfamily. With more than 130 members in Vitis, this particular family is known to regulate many different processes such as response to biotic and abiotic stress, development, reproduction, and response to hormones . Recently, Cramer et al.  reanalyzed the phylogeny of this family in Vitis. Out of 130 family members on the Vitis microarray, 91 changed significantly in response to dehydration between the genotypes with 99 changing significantly in response to dehydration between genotypes over time (Additional files 3 and 6). This indicates that ethylene and ethylene signaling may play important roles in the dehydration response.
Recently, Dubois et al.  classified ERF6 and ERF5 as the “master regulators” of leaf growth in Arabidopsis. A number of ERF6-like transcription factors changed in response to dehydration. VviERF6L3 and VviERF6L1 responded similarly by increasing rapidly at 1 h of dehydration in Cabernet Sauvignon leaves with little response in both Ramsey and Riparia Gloire (Fig. 2). In contrast, the ERF/AP2 transcription factors, VviRAP2.4 and VviERF11, displayed a similar pattern in Riparia Gloire increasing with a peak in expression at 2 h of dehydration. These data indicate a difference in the response between genotypes for multiple AP2/ERF transcription factors.
Other hormone signaling transcripts within the ABA core-signaling network
Transcripts involved in ABA signaling are also known to interact with other hormone signaling pathways. For example, ARR5 (Arabidopsis response regulator 5) is an essential component of cytokinin signaling . The transcript abundance of the closest Vitis ortholog, VviARR5, was decreased for all genotypes (Fig. 2). The transcript abundance in Cabernet Sauvignon decreased as early as 2 h of dehydration, however, at 24 h of dehydration the greatest changes were observed in Ramsey with a log2 fold change decrease of −2.0. The transcript abundance of VviARR4 significantly decreased in response to dehydration (Fig. 2). This particular gene is known to be involved in the ethylene signaling pathway .
In plants, ethylene is synthesized from S-adenosine-L-methionine (SAM), and 1-aminocyclopropane-1-carboxylate (ACC). The conversion of SAM to ACC is catalyzed by ACC synthase (ACS), which is followed by the oxidation of ACC to ethylene catalyzed by ACC oxidase (ACO) . ACS is the rate-limiting enzyme for ethylene biosynthesis. In our study, a number of putative ACS genes changed significantly in response to dehydration (Figs. 2 and 3). For example, transcript abundance of VviACS2 increased in all genotypes, most notably in Riparia Gloire at 1 h of dehydration (Fig. 2). In both Riparia Gloire and Cabernet Sauvignon, gene expression remained high throughout the experiment while Ramsey was increased at 4 h of dehydration. Another example that has been linked to the ABA core-signaling network is VviACS6 (Fig. 2). Interestingly, this gene was differentially regulated between the genotypes. The transcript abundance increased at 1 h in Cabernet Sauvignon, followed by Riparia Gloire, and the lowest expression was in Ramsey (Fig. 2).
Together, these data indicate that there are many changes in gene expression within the ABA core-signaling network as defined by Lumba et al. . Changes in gene expression involved in processes such as transport, transcription factor expression, kinase expression, and hormone signaling, highlight a complex regulatory network that is involved in the response to dehydration.
Dehydration induces significant changes in ethylene metabolism transcripts
Similar to ABA metabolism, massive changes in ethylene metabolism and signaling were observed. Within plants, small gene families encode multiple ACS and ACO genes that are known to be regulated differently depending on environmental, developmental, and hormonal signals [62–64]. In addition to the ACS genes previously mentioned, others changed significantly in response to dehydration (see Additional file 6 for gene annotations). VviACS4 and VviACS8-like increased in expression in both Cabernet Sauvignon and Riparia Gloire with little response in Ramsey (Fig. 3). Transcript abundance of yet another, VviACS7, increased most notably in Cabernet Sauvignon leaves at 1 h of dehydration. Interestingly a negative regulator of ethylene production, ETHYLENE OVERPRODUCER 1 (VviETO1), also increased in expression, particularly in Cabernet Sauvignon and Ramsey (Fig. 3). Previously, Yoshida et al.  used yeast-two hybrid assays to establish that the Arabidopsis ETO1 interacts with AtACS5 and not other ACS proteins surveyed. There is no clear ortholog to AtACS5 in Vitis; therefore, further investigation is needed.
Downstream of ACS, transcripts for ACC oxidase also displayed a significant change in response to dehydration. In particular, both VviACO4 and another ACO-like gene (VviACOL) increased in Cabernet Sauvignon leaves at 1 h of dehydration. VviACO1 increased notably in Ramsey and 2 h and 8 h in Riparia Gloire with little response in Cabernet Sauvignon (Fig. 3).
Dehydration induces significant changes in ethylene signaling transcripts
A large number of transcripts involved in ethylene signaling significantly changed in response to dehydration (Figs. 3 and 4). The transcript abundance of a number of ethylene receptors increased significantly in response to dehydration (Fig. 3). Ethylene receptors are broken down into two subfamilies based on conserved histidine kinase domains and are localized within the endoplasmic reticulum (for review see ). Ethylene receptor 2 (VviETR2) increased notably in Cabernet Sauvignon and Riparia Gloire at 1 h of dehydration with Ramsey responding later during the treatment. Ethylene response sensor 2 (VviERS2) and VviERS1 followed a similar pattern by increasing at 2 h of dehydration in Cabernet Sauvignon and Ramsey with a later response in Riparia Gloire (Fig. 3).
Additional transcripts involved in ethylene signaling were mapped (Fig. 4) using the ATTED-II database as a template . The closest Vitis orthologs were determined according to Gramene ( release 44 (January 2015)). A detailed investigation of all transcripts within this ethylene-signaling network is outside the scope of this study, however, a few key genes are discussed here.
The transcript abundance of a number of WRKY domain transcription factors involved in ethylene signaling increased in response to dehydration (Fig. 4). The transcript abundance of VviWRKY33 and VviWRKY40 increased rapidly within 1 h of dehydration in all genotypes. The transcript abundance of VviWRKY22 also increased rapidly within 1 h of dehydration in Riparia Gloire, but decreased at later time points in Ramsey indicating differences in the regulation of this transcription factor (Fig. 4). Recently, WRKY33 has been shown to bind directly to the promoter of ACS2 and ACS6 to induce gene expression in Arabidopsis . AP2/ERF transcription factors appear to regulate the expression of WRKY40 , further indicating cross-talk between the WRKY transcription factors and ethylene signaling.
Within this network, AP2/ERF transcription factors were also observed to change significantly. For example, the transcript abundance of VviERF6L1 increased at 1 h of dehydration with little response in Riparia Gloire and Ramsey (Fig. 4). The transcript abundance of VviERF104 also increased at 1 h of dehydration, however, this increase was observed in Riparia Gloire and Ramsey with little response in Cabernet Sauvignon. Finally, The transcript abundance of VviERF2 and VviERF1L2 increased at 1 h of dehydration and remained high throughout the experiment, notably in Cabernet Sauvignon and Riparia Gloire leaves.
Notable AP2/ERF transcript responses in drought sensitive Riparia Gloire
Riparia Gloire leaves are more sensitive to rapid dehydration by closing their stomata more quickly compared to both Cabernet Sauvignon and Ramsey . In this study, multiple AP2/ERF transcription factors increased in transcript abundance at 1 h in Riparia Gloire with little or no response in the other genotypes. For example, the transcript abundance of VviERF128 increased rapidly at 1 h with a slight response in Cabernet Sauvignon and little response in Ramsey leaves (Fig. 5). According to Cramer et al.  this gene does not have a clear ortholog to Arabidopsis indicating the possibility for a unique function in Vitis. Another AP2/ERF transcript that had a similar pattern of expression was VviERF098 (Fig. 5). The closest Arabidopsis ortholog AtERF98 (At3g23230), increases ascorbic acid (AsA) biosynthesis leading to an increase in salt tolerance . AsA has a number of roles in plants including as an antioxidant, protecting the plant from reactive oxygen species (ROS), which can result in enhanced tolerance to a variety of abiotic stresses [71, 72].
The transcript abundance of other AP2/ERF transcription factors in Riparia Gloire that increased rapidly in response to dehydration included VviERF055 and VviERF022 (Fig. 5). Interestingly, VviERF055 is closely related to TRANSLUCENT GREEN (TG), an ERF family transcription factor in Arabidopsis. TG binds directly to the promoter of multiple aquaporin genes, and overexpression results in increased drought tolerance . CYTOKININ RESPONSE FACTOR 2 (VviCRF2) is another example of an AP2/ERF transcript that increases in abundance in Riparia Gloire leaves with little to no response in the other genotypes (Fig. 5). The closest Arabidopsis ortholog, cytokinin response factor 2 (CRF2, At4g23750), interacts with other CRF proteins within the cytokinin signaling pathway .
The transcript abundance of two negative regulators of ethylene signaling, VviERF9 and VviERF022, also increased rapidly at 1 h of dehydration in Riparia Gloire leaves. VviERF9 belongs to subgroup VIII in the ERF family [23, 75], which is known to be involved in transcriptional repression [76–78]. In our study, VviERF9 expression increased rapidly in Riparia Gloire, to a lesser extent in Ramsey, with no response in Cabernet Sauvignon (Fig. 5). VviERF022 in Riparia Gloire and Ramsey leaves increased in expression at 1 h of dehydration, followed by a decrease at later time points. Little is known about the function of this gene in Vitis. Nowak et al.  demonstrated the importance of AtERF022, the closest ortholog, in somatic embryogenesis in Arabidopsis. This study provides direct evidence that ERF022 regulates the biosynthesis and signaling of ethylene. These data indicate some unique changes occurred in Riparia Gloire leaves within the AP2/ERF family. Changes in both positive and negative regulators of ethylene signaling also indicate a complex regulatory network requiring further investigation.
Together, these data provide further evidence for changes in the expression of genes involved in both ABA and ethylene metabolism and signaling in response to dehydration. These changes were also observed to be different between the genotypes surveyed indicating different types of transcriptional regulation.
Network analysis: WGCNA
Weighted gene coexpression network analysis (WGCNA) was performed to get a better understanding of which genes within these complex signaling networks were the most connected hubs. Thirty modules or gene networks were detected, assigned color names, and correlated to genotype and treatment effects over time (Fig. 6). All genes were correlated with these 30 colored modules; the grey category is not a true module, but a place to put all the leftover genes not correlated well enough with one of the significant colored modules. In addition, a measure called the kME (module eigengene-based connectivity) was calculated for each gene to every module (Additional file 7). Genes with a kME score of 1 are perfectly positively correlated with that particular module (network) and scores of −1 are perfectly negatively correlated to the module. The advantage of computing kME scores for each module is that genes can act as a hub in more than one module or network. Genes with the highest kME score have the most connectivity and are therefore designated the largest hubs. The largest hubs were identified along with the top GO categories overrepresented in each module (Additional file 8). All top hub genes had kME scores of 0.92 or higher. Some of the top hub genes in each module were of unknown function, making them interesting candidates for future research (Additional files 7 and 8).
In contrast to previous GO enrichment analyses done with Genotype x Treatment or the Genotype x Treatment x Time subsets of genes, GO enrichment analysis of the individual modules identified with WGCNA provided much more significant results with more biological meaning. Some of the modules were highly enriched (based upon their p-value) in gene ontologies such as plum2 (translation), firebrick4 (photosynthesis) and yellow2 (protein amino acid phosphorylation); a few were less distinct or not significant, such as yellow3, blue, and tan4 (Additional file 8).
The most interesting modules were modules (lightsteelblue, darkolivegreen, orangered3, and yellow3) because they were correlated positively with dehydration and genotype (Fig. 6). The most interesting module is the yellow3 module, because the heat map for genotypes was most highly correlated with the genotype order of drought tolerance: Ramsey > Cabernet Sauvignon > Riparia Gloire, and there was an increase in response to dehydration. The top hub gene in the yellow3 module was a gene for an unknown protein with a TPR domain (Additional file 8; Fig. 7). TPR domains are structural motifs that facilitate protein-protein interactions and assembly of protein complexes. There were several other TPR domain protein genes in the top of this module. These genes are clear targets for future research due to lack of knowledge of these genes in any plant species.
The yellow3 module is overrepresented with a significant but relatively high p-value (1.28 × 10−4) in gas (oxygen) transport. Basically this means that either this network of genes is poorly annotated or there is a general mix of all kinds of gene functions that is not very distinct from the entire genome. Nevertheless, there is a very interesting subset of hub genes all having kME values above 0.80, including the before mentioned ABA core-signaling transcription factors, ABF2 and ABI5 (Fig. 7; Additional file 7). Other notable hub genes in this module are a SNF1-related protein kinase (VviSnRK1; a central regulator of metabolism), a cinnamyl alcohol dehydrogenase (VviCAD; phenylpropanoid metabolism), a β-carotene hydroxylase (VviBCH; carotenoid metabolism), GSDL Lipase 1 (VviGLIP1; ethylene signaling) and a universal stress protein (VviUSP; defense). The one gene that stands out in this group is VviABI5, because of the differences displayed in transcript abundance between the genotypes in their response to dehydration with Ramsey > Cabernet Sauvignon > Riparia Gloire. There are 81 hub genes in this network with a kME higher than 0.80, indicating a very complex network.
The other very interesting module is the lightsteelblue module (Additional file 7). A number of the genes with high kMEs in the yellow3 module also have high kMEs in the lightsteelblue module, including VviABI5, VviABF2, VviSnRK1, some wound-inducible genes, and genes encoding late-embryogenesis abundant (LEA) proteins or dehydrins. This module is significantly enriched in the ethylene signaling GO category. The top hub gene is an unknown gene, whose ortholog in Arabidopsis is induced by Al (Fig. 8; Additional files 7 and 8). Other top hub genes with similar but slightly different expression patterns are Gibberellin Insensitive Dwarf 1B (VviGID1B; a gibberellin receptor), an uncharacterized VviMYB, RAS-related Nuclear Protein (VviRAN1; GTPase signaling), VviERF1 (ethylene signaling), VviNCED3 (ABA biosynthesis), and Indeterminate Domain 2 (VviIDD2; a C2H2-type zinc finger protein). Interestingly, VviERF1, which is at the start of the ERF transcription factor cascade is more highly expressed in Riparia Gloire in response to dehydration (Fig. 8). Thus there is a clear overlap in ABA and ethylene signaling hubs in these two modules, which supports the argument that there is strong crosstalk between these two hormone signaling pathways.
Top hub genes from other modules with high GO enrichment categories such as translation, photosynthesis, phenylpropanoid metabolism were less interesting because the transcripts in these sets were less correlated with drought tolerance, and therefore will not be discussed further. We have only just scratched the surface of the analyses of these data. This dataset (Additional file 7) is very deep and can be exploited for years to come. Many other fascinating networks are present within this dataset but there is not enough space and time to describe them here.
There are thousands of different grapevine genotypes leading to a large pool of natural genetic diversity [80, 81]. The phylloxera devastation of the European wine industry during the 19th century resulted in the introduction of North American genotypes as rootstocks for biotic and abiotic stress resistance . With global climate change, a better understanding of the genes involved in drought tolerance will be needed. Previous work by Padgett-Johnson et al.  clearly showed differences in water-use-efficiency between different genotypes of grapevine in the field. However, studies like this require multiple years for the establishment of an experimental vineyard.
The strengths and weaknesses of the rapid dehydration assay
To complement previous field studies, Hopper et al.  developed a rapid dehydration assay that is amenable to large-scale phenotyping studies. The assay takes only a few hours and can precisely control the level of dehydration. It is simple, not requiring expensive or sophisticated equipment. It can be done at any time of the year as long as plant material is available, such as from a greenhouse. The results from the rapid dehydration assay are consistent with field assays, indicating that it is a suitable assay for evaluating dehydration tolerance. This assay is excellent as a first screen for genetic tendencies.
The screen only reflects the inherent tendencies of a leaf. The leaf is isolated from the rest of the plant and cannot sense signals from other parts of the plant. For example, in an earlier study, massive changes in proteins were detected in Cabernet Sauvignon shoots that had been gradually dehydrated . These changes occurred before any photosynthetic or growth symptoms indicating that the shoot detected a feed-forward signal from another part of the plant, such as the roots.
This rapid dehydration assay by its nature cannot look at the more gradual changes that occur in a plant in a more natural setting. Nevertheless, it was able to detect differences between genotypes that are reflective of their behavior in the field.
Sensitivity and timing of stress responses
One of the strengths of our transcriptomic approach in this study was the collection of samples along a detailed time course. This allowed detection of genotype sensitivity to dehydration. Sensitivity could reflect sensitivity to injury and thus indicate a stress-sensitive (dehydration sensitive) plant or it could indicate a plant with heightened awareness of impending drought, giving it more time to respond and thus fostering dehydration resistance. The use of genotypes differing in their drought tolerance allowed us to distinguish between these sensitivities. ABA signaling appears to be key in this response with the drought tolerant genotype, Ramsey, showing heightened responses in ABA signaling and decreased responses in ethylene signaling relative to the drought sensitive genotype, Riparia Gloire, which showed slower responses in ABA signaling and heightened responses in ethylene signaling. The intermediate genotype, Cabernet Sauvignon, was consistent with these responses, being intermediate in both ABA and ethylene signaling.
The importance of ABA and ethylene responses to dehydration
ABA is a major plant hormone involved in plant responses to dehydration. Significant differences in transcript abundance of genes involved in ABA metabolism and signaling were observed between the genotypes. Previously, Hopper et al.  found that application of ABA to leaves prior to dehydration decreased the overall amount and rate of water lost. Differences between the genotypes in their stomatal conductance after application at various concentrations of ABA also indicate differences in ABA sensitivity. Differences observed in water loss are presumably through differences in stomatal conductance and regulation of gene transcription, protein synthesis and other signaling pathways (for review see ).
Based on both the standard and network analysis approaches, significant changes in the transcript abundance of genes involved in ABA metabolism and signaling were observed (Figs. 2 and 3). For example, the transcript abundance of the three NCED genes in Vitis that catalyze the rate-limiting step in ABA biosynthesis increased in response to dehydration. Interestingly, VviNCED6, is increased in both Riparia Gloire and Ramsey with a different response in Cabernet Sauvignon. This indicates differences in the regulation of these key genes to dehydration. Recently, Lumba et al.  created an ABA interactome using a systems biology approach. Utilizing transcriptomic and proteomic data, these authors were able to create a detailed interactome with over 500 interactions, highlighting the complexity of ABA signaling. Studies like this also highlight the need and utility for transcriptomic data such as in the study presented here. Our study provides further evidence for differences in ABA metabolism and signaling at the transcript level that may be involved in dehydration adjustment for the Vitis genotypes.
Key candidate factors appear to be the transcription factors, ABI5 and ABF2, at the start of the ABA transcriptional signaling pathway. In particular, the transcriptional response of VviABI5 to dehydration in Ramsey was more sensitive than Cabernet Sauvignon, which was more sensitive than Riparia Gloire. This makes it a good target for future research.
ABI5 was first cloned in Arabidopsis . Its expression is highest in seeds, but a low level of expression is observed in vegetative tissues. The expression of some AtLEA genes is dependent upon AtABI5 expression in Arabidopsis. Note that there were several VviLEA genes within the top kMEs (>0.80) in the yellow3 module (Additional file 7). ABI5 expression is largely associated with seed development and not vegetative organs or tissues . However recent evidence has linked ABI5 to vegetative functions such as the induction of flowering , leaf senescence , drought tolerance , and stomatal movement . With these exciting new findings, it is necessary to focus more research on the function of ABI5 in grapevine vegetative organs.
The SnRK1 γ-subunit, another hub in the yellow3 module is a part of the SnRK1 complex that is a central regulator of metabolism. SnRK1 responds to sugar and ABA signaling . PP2Cs (ABI1 and AHG3) directly dephosphorylate the SnRK1 α-subunit ; thus, the inhibition of PP2C activity by ABA directly stimulates SnRK1 expression and activity, complementing the stress response through energy regulation and coordination.
Ethylene is another plant hormone involved in dehydration stress signaling (for review see ). The ethylene response may reflect dehydration adjustment, damage control or senescence. The ERF domain of ERF transcription factors contains 60 to 70 conserved amino acids and was first identified in four DNA binding proteins NtERF1-4 from Nicotiana tabacum . Group I ERF transcription factors have previously been shown to be important in mediating drought tolerance in plants [93–95]. VviERF1 had a high kME in the lightsteelblue module, which was significantly enriched in ethylene signaling gene ontologies. The closest orthologs to VviRAP2.4 and VviERF055 belong to Group I in Arabidopsis  and were both up-regulated in response to dehydration, most notably in Riparia Gloire leaves (Figs. 2 and 5).
ERF5 and ERF6 are also important transcription factors involved in drought signaling and have been called the “master regulators” in leaf growth in response to environmental changes . In Cabernet Sauvignon leaves, a number of VviERF6-like transcription factors increased in response to dehydration (Fig. 2). A similar trend in the expression of VviERF6-like genes was seen as sugar (°Brix) levels increased within the skins of ripening Cabernet Sauvignon berries . This may indicate a similar level of regulation of these genes in response to dehydration, ripening or other stress responses.
Similar to ABA, transcripts involved in ethylene metabolism and signaling also changed significantly in response to dehydration (Figs. 3 and 4). A number of transcripts were observed to be changing significantly based on the Genotype x Treatment and Genotype x Treatment x Time interaction terms indicating a different response between the genotypes in response to dehydration and over time.
There is substantial evidence that there is crosstalk between ABA and ethylene in response to dehydration [88, 96]. Overexpression of an ethylene responsive transcription factor RAV2 along with the ABA-response transcription factor ABI5 had synergistic effects on drought tolerance in cotton . Interestingly, VviRAV2 has a high kME (0.81) in the lightsteelblue module along with VviABI5 (Additional file 7). Lumba et al.  provide evidence for hormonal crosstalk with a number of genes involved in ethylene metabolism and signaling within the ABA core network. For example, ACS catalyzes the rate-limiting step in ethylene biosynthesis and ACS6 in Arabidopsis is regulated by ABA Insensitive 1 (ABI1), a negative regulator of ABA signaling . Multiple ACS genes in Vitis changed significantly between the genotypes in response to dehydration (Figs. 2 and 3).
In summary, the leaf dehydration time-series assay allowed the detection of a very large number of transcriptional changes in a coordinated fashion. Key genes were identified by a standard a posteriori analysis that involved mapping known biochemical pathways in leaves responding to dehydration (e.g. ABA and ethylene metabolism and signaling pathways). An a priori data analysis approach using WGCNA proved more powerful; it confirmed the results from the standard approach and it identified 30 distinct modules (networks), most of which had highly enriched GO categories that were biologically and functionally relevant. In addition, WGCNA enabled the identification of gene hubs in these modules that are likely to be important and significant operators within these networks.
The results from this study confirmed our hypothesis that ABA signaling pathways were different between the grapevine genotypes. Furthermore, the results indicate that the dehydration response had substantial cross-talk between the ABA and ethylene signaling pathways. Some of the genes in the ABA and ethylene signaling pathways were highly connected hubs and were correlated with drought tolerance. A number of interesting unknown genes were also identified and associated with these pathways. VviABI5 is one of the best candidate hubs for drought tolerance because its expression was more rapid and higher in the drought tolerance genotypes and because of its known functions in ABA signaling and drought tolerance. The facts that this gene is normally associated with seed dormancy, that the transcript abundance of VviABI5 in leaves increased significantly in a species (Vitis champinii) that has evolved in a hot and dry climate, and did not respond substantially in the leaves of a species (Vitis riparia) that has evolved in a cooler and wetter climate, make this gene particularly interesting.
This study provides a very rich data set that can be further explored for new discoveries. Only the tip of this data “iceberg” has been discussed here. Future research will focus on the further elucidation of the regulation of these gene networks and on gradual responses of the root to dehydration, since multiple lines of evidence indicate it may act as the initial sensor for drought stress signaling.
Boyer JS. Plant productivity and environment. Science. 1982;218(4571):443–8.
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.
Hopper DW, Ghan R, Cramer GR. A rapid dehydration leaf assay reveals stomatal response differences in grapevine genotypes. Hortic Res. 2014;1:2.
Hannah L, Roehrdanz PR, Ikegami M, Shepard AV, Shaw MR, Tabor G, Zhi L, Marquet PA, Hijmans RJ. Climate change, wine, and conservation. Proc Natl Acad Sci U S A. 2013;110(17):6907–12.
van Leeuwen C, Schultz HR, Garcia de Cortazar-Atauri I, Duchene E, Ollat N, Pieri P, Pieri P, Bois B, Goutouly JP, Quenol H, Touzard JM et al. Why climate change will not dramatically decrease viticultural suitability in main wine-producing areas by. Proc Natl Acad Sci U S A 2013. 2050;110(33):E3051–3052.
Lowe KM, Walker MA. Genetic linkage map of the interspecific grape rootstock cross Ramsey (Vitis champinii) x Riparia Gloire (Vitis riparia). Theor Appl Genet. 2006;112(8):1582–92.
Marguerit E, Brendel O, Lebon E, Van Leeuwen C, Ollat N. Rootstock control of scion transpiration and its acclimation to water deficit are controlled by different genes. New Phytol. 2012;194(2):416–29.
Gambetta GA, Manuck CM, Drucker ST, Shaghasi T, Fort K, Matthews MA, Walker MA, McElrone AJ. The relationship between root hydraulics and scion vigour across Vitis rootstocks: what role do root aquaporins play? J Exp Bot. 2012;63(18):6445–55.
Candolfi-Vasconcelos MC, Candolfi MP, Kohlet W. Retranslocation of carbon reserves from the woody storage tissues into the fruit as a response to defoliation stress during the ripening period in Vitis vinifera L. Planta. 1994;192(4):567–73.
Padgett-Johnson M, Williams LE, Walker MA. The influence of Vitis riparia rootstock on water relations and gas Exchange of Vitis vinifera cv. Carignane scion under non-irrigated conditions. Am J Enol Vitic. 2000;51(2):137–43.
Soar CJ, Dry PR, Loveys BR. Scion photosynthesis and leaf gas exchange in Vitis vinifera L. cv. Shiraz: Mediation of rootstock effects via xylem sap ABA. Aust J Grape Wine Res. 2006;12(2):82–96.
Padgett-Johnson M, Williams LE, Walker MA. Vine water relations, gas exchange, and vegetative growth of seventeen vitis species grown under irrigated and nonirrigated conditions in California. J Amer Soc Hort Sci. 2003;128(2):269–76.
Dong J, Horvath S. Understanding network concepts in modules. BMC Syst Biol. 2007;1:24.
Langfelder P, Horvath S. Eigengene networks for studying the relationships between co-expression modules. BMC Syst Biol. 2007;1:54.
Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. 2008;9:559.
Song L, Langfelder P, Horvath S. Comparison of co-expression measures: mutual information, correlation, and model based indices. BMC Bioinformatics. 2012;13:328.
Goda H, Sasaki E, Akiyama K, Maruyama-Nakashita A, Nakabayashi K, Li W, Ogawa M, Yamauchi Y, Preston J, Aoki K et al. The AtGenExpress hormone and chemical treatment data set: experimental design, data evaluation, model data analysis and data access. Plant J. 2008;55(3):526–42.
Retsch: Solutions in Milling and Sieving. [http://www.retsch.com/]
Gambino G, Perrone I, Gribaudo I. A rapid and effective method for RNA extraction from different tissues of grapevine and other woody plants. Phytochem Anal. 2008;19(6):520–5.
Tattersall EAR, Ergul A, AlKayal F, Deluc L, Cushman JC, Cramer GR. Comparison of methods for isolating high-quality RNA from leaves of grapevine. Am J Enol Vitic. 2005;56:400–6.
Rattanakon S, Ghan R, Gambetta G, Deluc L, Schlauch KA, Cramer GR. Transcriptomic analysis reveals that abscisic acid signaling differs in different organs of grapevine (Vitis vinifera L.). BMC Plant Biol 2016. (in press).
Cramer GR, Ghan R, Schlauch KA, Tillett RL, Heymann H, Ferrarini A, Delledonne M, Zenoni S, Fasoli M, Pezzotti M. Transcriptomic analysis of the late stages of grapevine (Vitis vinifera cv. Cabernet Sauvignon) berry ripening reveals significant induction of ethylene signaling and flavor pathways in the skin. BMC Plant Biol. 2014;14(1):370.
Ghan R, Van Sluyter SC, Hochberg U, Degu A, Hopper DW, Tillet RL, Schlauch KA, Haynes PA, Fait A, Cramer GR. Five omic technologies are concordant in differentiating the biochemical characteristics of the berries of five grapevine (Vitis vinifera L.) cultivars. BMC Genomics. 2015;16(1):946.
Sokal R, Rohlf F. Biometry: the principles and practice of statistics in biological research. 3rd ed. New York: W.H.Freeman; 1995.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Royal Stat Soc: Series B (Stat Methodol). 1995;57:289–300.
Knudsen S. Guide to analysis of DNA microarray data. 2nd ed. New York: John Wiley & Sons; 2004.
Baldi P, Hatfield GW. DNA microarrays and gene expression: from experiments to data analysis and modeling. Cambridge: Cambridge University Press; 2002.
Gordon AD. Classification. Boca Raton: Chapman & Hall/CRC; 1999.
Stekel D. Microarray bioinformatics. Cambridge: Cambridge University Press; 2003.
Maere S, Heymans K, Kuiper M. BiNGO: a Cytoscape plugin to assess overrepresentation of gene ontology categories in biological networks. Bioinformatics. 2005;21(16):3448–9.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.
Kinsella RJ, Kahari A, Haider S, Zamora J, Proctor G, Spudich G, Almeida-King J, Staines D, Derwent P, Kerhornou A et al. Ensembl BioMarts: a hub for data retrieval across taxonomic space. Database (Oxford). 2011;2011:bar030.
Schwartz SH, Qin X, Zeevaart JA. Elucidation of the indirect pathway of abscisic acid biosynthesis by mutants, genes, and enzymes. Plant Physiol. 2003;131(4):1591–601.
Endo A, Sawada Y, Takahashi H, Okamoto M, Ikegami K, Koiwai H, Seo M, Toyomasu T, Mitsuhashi W, Shinozaki K et al. Drought induction of Arabidopsis 9-cis-epoxycarotenoid dioxygenase occurs in vascular parenchyma cells. Plant Physiol. 2008;147(4):1984–93.
Grimplet J, Adam-Blondon AF, Bert PF, Bitz O, Cantu D, Davies C, Delrot S, Pezzotti M, Rombauts S, Cramer GR. The grapevine gene nomenclature system. BMC Genomics. 2014;15:1077.
Burla B, Pfrunder S, Nagy R, Francisco RM, Lee Y, Martinoia E. Vacuolar transport of abscisic acid glucosyl ester is mediated by ATP-binding cassette and proton-antiport mechanisms in Arabidopsis. Plant Physiol. 2013;163(3):1446–58.
Zhang G, Duan C, Wang Y, Wang Y, Ji K, Xu H, Dai S, Li Q, Chen P, Sun Y et al. The expression pattern of β-glucosidase genes (VvBGs) during grape berry maturation and dehydration stress. Plant Growth Regul. 2012;70(2):105–14.
Kuromori T, Miyaji T, Yabuuchi H, Shimizu H, Sugimoto E, Kamiya A, Moriyama Y, Shinozaki K. ABC transporter AtABCG25 is involved in abscisic acid transport and responses. Proc Natl Acad Sci U S A. 2010;107(5):2361–6.
Kuromori T, Sugimoto E, Shinozaki K. Intertissue signal transfer of abscisic acid from vascular cells to guard cells. Plant Physiol. 2014;164(4):1587–92.
Kuromori T, Shinozaki K. ABA transport factors found in Arabidopsis ABC transporters. Plant Signal Behav. 2010;5(9):1124–6.
Lumba S, Toh S, Handfield LF, Swan M, Liu R, Youn JY, Cutler SR, Subramaniam R, Provart N, Moses A et al. A mesoscale abscisic acid hormone interactome reveals a dynamic signaling landscape in Arabidopsis. Dev Cell. 2014;29(3):360–72.
Ware DH, Jaiswal P, Ni J, Yap IV, Pan X, Clark KY, Teytelman L, Schmidt SC, Zhao W, Chang K et al. Gramene, a tool for grass genomics. Plant Physiol. 2002;130(4):1606–13.
Ma Y, Szostkiewicz I, Korte A, Moes D, Yang Y, Christmann A, Grill E. Regulators of PP2C phosphatase activity function as abscisic acid sensors. Science. 2009;324(5930):1064–8.
Park SY, Fung P, Nishimura N, Jensen DR, Fujii H, Zhao Y, Lumba S, Santiago J, Rodrigues A, Chow TF et al. Abscisic acid inhibits type 2C protein phosphatases via the PYR/PYL family of START proteins. Science. 2009;324(5930):1068–71.
Pandey S, Nelson DC, Assmann SM. Two novel GPCR-type G proteins are abscisic acid receptors in Arabidopsis. Cell. 2009;136(1):136–48.
Hauser F, Waadt R, Schroeder JI. Evolution of abscisic acid synthesis and signaling mechanisms. Curr Biol. 2011;21(9):R346–355.
Yoshida T, Fujita Y, Maruyama K, Mogami J, Todaka D, Shinozaki K, Yamaguchi-Shinozaki K. Four Arabidopsis AREB/ABF transcription factors function predominantly in gene expression downstream of SnRK2 kinases in abscisic acid signalling in response to osmotic stress. Plant Cell Environ. 2015;38(1):35–49.
Vahisalu T, Kollist H, Wang Y-F, Nishimura N, Chan W-Y, Valerio G, Lamminmäki A, Brosché M, Moldau H, Desikan R et al. SLAC1 is required for plant guard cell S-type anion channel function in stomatal signalling. Nature. 2008;452(7186):487–91.
Liu KH, Tsay YF. Switching between the two action modes of the dual-affinity nitrate transporter CHL1 by phosphorylation. EMBO J. 2003;22(5):1005–13.
Yamasaki H, Sakihama Y. Simultaneous production of nitric oxide and peroxynitrite by plant nitrate reductase: in vitro evidence for the NR-dependent formation of active nitrogen species. FEBS Lett. 2000;468(1):89–92.
Coello P, Hey SJ, Halford NG. The sucrose non-fermenting-1-related (SnRK) family of protein kinases: potential for manipulation to improve stress tolerance and increase yield. J Exp Bot. 2011;62(3):883–93.
Arisz SA, Testerink C, Munnik T. Plant PA signaling via diacylglycerol kinase. Biochim Biophys Acta. 2009;1791(9):869–75.
Stracke R, Ishihara H, Huep G, Barsch A, Mehrtens F, Niehaus K, Weisshaar B. Differential regulation of closely related R2R3-MYB transcription factors controls flavonol accumulation in different parts of the Arabidopsis thaliana seedling. Plant J. 2007;50(4):660–77.
Nakabayashi R, Yonekura-Sakakibara K, Urano K, Suzuki M, Yamada Y, Nishizawa T, Nishizawa T, Matsuda F, Kojima M, Sakakibara H, Shinozaki K et al. Enhancement of oxidative and drought tolerance in Arabidopsis by overaccumulation of antioxidant flavonoids. Plant J. 2014;77(3):367–79.
Ding Y, Liu N, Virlouvet L, Riethoven JJ, Fromm M, Avramova Z. Four distinct types of dehydration stress memory genes in Arabidopsis thaliana. BMC Plant Biol. 2013;13:229.
Licausi F, Ohme-Takagi M, Perata P. APETALA2/Ethylene Responsive Factor (AP2/ERF) transcription factors: mediators of stress responses and developmental programs. New Phytol. 2013;199(3):639–49.
Dubois M, Skirycz A, Claeys H, Maleux K, Dhondt S, De Bodt S, Vanden Bossche R, De Milde L, Yoshizumi T, Matsui M et al. Ethylene Response Factor6 acts as a central regulator of leaf growth under water-limiting conditions in Arabidopsis. Plant Physiol. 2013;162(1):319–32.
Kurepa J, Li Y, Smalle JA. Proteasome-dependent proteolysis has a critical role in fine-tuning the feedback inhibition of cytokinin signaling. Plant Signal Behav. 2013;8(3):e23474.
Mira-Rodado V, Veerabagu M, Witthoft J, Teply J, Harter K, Desikan R. Identification of two-component system elements downstream of AHK5 in the stomatal closure response of Arabidopsis thaliana. Plant Signal Behav. 2012;7(11):1467–76.
Yang SF, Hoffman NE. Ethylene biosynthesis and its regulation in higher plants. Annu Rev Plant Physiol. 1984;35:155–89.
Rottmann WH, Peter GF, Oeller PW, Keller JA, Shen NF, Nagy BP, Taylor LP, Campbell AD, Theologis A. 1-Aminocyclopropane-1-carboxylate synthase in tomato is encoded by a multigene family whose transcription is induced during fruit and floral senescence. J Mol Biol. 1991;222(4):937–61.
El-Sharkawy I, Kim WS, Jayasankar S, Svircev AM, Brown DC. Differential regulation of four members of the ACC synthase gene family in plum. J Exp Bot. 2008;59(8):2009–27.
Van de Poel B, Van Der Straeten D. 1-aminocyclopropane-1-carboxylic acid (ACC) in plants: more than just the precursor of ethylene! Front Plant Sci. 2014;5:640.
Yoshida H, Nagata M, Saito K, Wang KL, Ecker JR. Arabidopsis ETO1 specifically interacts with and negatively regulates type 2 1-aminocyclopropane-1-carboxylate synthases. BMC Plant Biol. 2005;5:14.
Gallie DR. Appearance and elaboration of the ethylene receptor family during land plant evolution. Plant Mol Biol. 2015;87(4–5):521–39.
Obayashi T, Okamura Y, Ito S, Tadaka S, Aoki Y, Shirota M, Kinoshita K. ATTED-II in 2014: evaluation of gene coexpression in agriculturally important plants. Plant Cell Physiol. 2014;55(1):e6.
Li G, Meng X, Wang R, Mao G, Han L, Liu Y, Zhang S. Dual-level regulation of ACC synthase activity by MPK3/MPK6 cascade and its downstream WRKY transcription factor during ethylene induction in Arabidopsis. PLoS Genet. 2012;8(6):e1002767.
Koyama T, Nii H, Mitsuda N, Ohta M, Kitajima S, Ohme-Takagi M, Sato F. A regulatory cascade involving class II ETHYLENE RESPONSE FACTOR transcriptional repressors operates in the progression of leaf senescence. Plant Physiol. 2013;162(2):991–1005.
Zhang Z, Wang J, Zhang R, Huang R. The ethylene response factor AtERF98 enhances tolerance to salt through the transcriptional activation of ascorbic acid synthesis in Arabidopsis. Plant J. 2012;71(2):273–87.
Upadhyaya, Hemavathi CP, Akula N, Young KE, Chun SC, Kim DH, Park SW. Enhanced ascorbic acid accumulation in transgenic potato confers tolerance to various abiotic stresses. Biotechnol Lett. 2009;32(2):321–30.
Huang C, He W, Guo J, Chang X, Su P, Zhang L. Increased sensitivity to salt stress in an ascorbate-deficient Arabidopsis mutant. J Exp Bot. 2005;56(422):3041–9.
Zhu D, Wu Z, Cao G, Li J, Wei J, Tsuge T, Gu H, Aoyama T, Qu LJ. TRANSLUCENT GREEN, an ERF family transcription factor, controls water balance in Arabidopsis by activating the expression of aquaporin genes. Mol Plant. 2014;7(4):601–15.
Cutcliffe JW, Hellmann E, Heyl A, Rashotte AM. CRFs form protein–protein interactions with each other and with members of the cytokinin signalling pathway in Arabidopsis via the CRF domain. J Exp Bot. 2011;62(14):4995-5002.
Licausi F, Giorgi FM, Zenoni S, Osti F, Pezzotti M, Perata P. Genomic and transcriptomic analysis of the AP2/ERF superfamily in Vitis vinifera. BMC Genomics. 2010;11:719.
Nakano T, Suzuki K, Fujimura T, Shinshi H. Genome-wide analysis of the ERF gene family in Arabidopsis and rice. Plant Physiol. 2006;140(2):411–32.
Hiratsu K, Matsui K, Koyama T, Ohme-Takagi M. Dominant repression of target genes by chimeric repressors that include the EAR motif, a repression domain, in Arabidopsis. Plant J. 2003;34(5):733–9.
Ohta M, Matsui K, Hiratsu K, Shinshi H, Ohme-Takagi M. Repression domains of class II ERF transcriptional repressors share an essential motif for active repression. Plant Cell. 2001;13(8):1959–68.
Nowak K, Wojcikowska B, Gaj MD. ERF022 impacts the induction of somatic embryogenesis in Arabidopsis through the ethylene-related pathway. Planta. 2015;241(4):967–85.
This P, Lacombe T, Thomas MR. Historical origins and genetic diversity of wine grapes. Trends Genet. 2006;22(9):511–9.
Myles S, Boyko AR, Owens CL, Brown PJ, Grassi F, Aradhya MK, Prins B, Reynolds A, Chia JM, Ware D et al. Genetic structure and domestication history of the grape. Proc Natl Acad Sci U S A. 2011;108(9):3530–5.
Cramer GR, Van Sluyter SC, Hopper DW, Pascovici D, Keighley T, Haynes PA. Proteomic analysis indicates massive changes in metabolism prior to the inhibition of growth and photosynthesis of grapevine (Vitis vinifera L.) in response to water deficit. BMC Plant Biol. 2013;13:49.
Seki M, Umezawa T, Urano K, Shinozaki K. Regulatory metabolic networks in drought stress responses. Curr Opin Plant Biol. 2007;10(3):296–302.
Finkelstein RR, Lynch TJ. The Arabidopsis abscisic acid response gene ABI5 encodes a basic leucine zipper transcription factor. Plant Cell. 2000;12(4):599–609.
Fujita Y, Fujita M, Shinozaki K, Yamaguchi-Shinozaki K. ABA-mediated transcriptional regulation in response to osmotic stress in plants. J Plant Res. 2011;124(4):509–25.
Wang Y, Li L, Ye T, Lu Y, Chen X, Wu Y. The inhibitory effect of ABA on floral transition is mediated by ABI5 in Arabidopsis. J Exp Bot. 2013;64(2):675–84.
Sakuraba Y, Jeong J, Kang MY, Kim J, Paek NC, Choi G. Phytochrome-interacting transcription factors PIF4 and PIF5 induce leaf senescence in Arabidopsis. Nat Commun. 2014;5:4636.
Mittal A, Gampala SS, Ritchie GL, Payton P, Burke JJ, Rock CD. Related to ABA-Insensitive3(ABI3)/Viviparous1 and AtABI5 transcription factor coexpression in cotton enhances drought stress adaptation. Plant Biotechnol J. 2014;12(5):578–89.
Qiao Z, Li CL, Zhang W. WRKY1 regulates stomatal movement in drought-stressed Arabidopsis thaliana. Plant Mol Biol. 2016;91(1-2):53–65.
Jossier M, Bouly JP, Meimoun P, Arjmand A, Lessard P, Hawley S, Grahame HD, Thomas M. SnRK1 (SNF1-related kinase 1) has a central role in sugar and ABA signalling in Arabidopsis thaliana. Plant J. 2009;59(2):316–28.
Rodrigues A, Adamo M, Crozet P, Margalha L, Confraria A, Martinho C, Elias A, Rabissi A, Lumbreras V, Gonzalez-Guzman M et al. ABI1 and PP2CA phosphatases are negative regulators of Snf1-related protein kinase1 signaling in Arabidopsis. Plant Cell. 2013;25(10):3871–84.
Ohme-Takagi M, Shinshi H. Ethylene-inducible DNA binding proteins that interact with an ethylene-responsive element. Plant Cell. 1995;7(2):173–82.
Zhang JY, Broeckling CD, Blancaflor EB, Sledge MK, Sumner LW, Wang ZY. Overexpression of WXP1, a putative Medicago truncatula AP2 domain-containing transcription factor gene, increases cuticular wax accumulation and enhances drought tolerance in transgenic alfalfa (Medicago sativa). Plant J. 2005;42(5):689–707.
Saleh A, Lumbreras V, Lopez C, Kizis ED-PD, Pagès M. Maize DBF1-interactor protein 1 containing an R3H domain is a potential regulator of DBF1 activity in stress responses. Plant J. 2006;46(5):747–57.
Lin R-C, Park H-J, Wang H-Y. Role of Arabidopsis RAP2.4 in Regulating Light- and Ethylene-Mediated Developmental Processes and Drought Stress Tolerance. Mol Plant. 2008;1(1):42–57.
Wilkinson S, Davies WJ. Drought, ozone, ABA and ethylene: new insights from cell to plant to community. Plant Cell Environ. 2010;33(4):510–25.
Ludwikow A, Ciesla A, Kasprowicz-Maluski A, Mitula F, Tajdel M, Galganski L, Ziolkowski PA, Kubiak P, Malecka A, Piechalak A et al. Arabidopsis protein phosphatase 2C ABI1 interacts with type I ACC synthases and is involved in the regulation of ozone-induced ethylene biosynthesis. Mol Plant. 2014;7(6):960–76.
Barrett T, Wilhite SE, Ledoux P, Evangelista C, Kim IF, Tomashevsky M, Marshall KA, Phillippy KH, Sherman PM, Holko M et al. NCBI GEO: archive for functional genomics data sets--update. Nucleic Acids Res. 2013;41(Database issue):D991–995.
This material is based upon work that is supported by the National Institute of Food and Agriculture, U.S. Department of Agriculture, Hatch Grant (NEV00345) and was supported by a grant from the National Institute of General Medical Sciences (P20GM103440).
Availability of supporting data
DWH designed and conducted the experiments, analyzed the data and wrote the body of the paper; RG helped in the data analysis; KAS performed statistical analysis; GRC acquired the funding and supervised all aspects of the project. All authors reviewed, edited and approved the final version of the manuscript.
The authors declare that they have no competing interests.
Consent to publish
Principal component analysis (PCA) of expression data. Symbol names refer to genotype (Ram = Ramsey, Rip = Riparia and CS = Cabernet Sauvignon); the following number is the time in hours (1, 2, 4, 8, and 24), then there is a separating “.”, followed by the treatment (C = control, S = stress or dehydration), and the final number refers to the sample replicate number. (PDF 19 kb)
WGCNA sample dendrogram and trait heatmap of expression data. (PDF 29 kb)
Annotation, transcript abundance values, and statistics of all genes on the NimbleGen Grape Whole-Genome microarray. (XLSX 20763 kb)
Overrepresented GO categories of all significantly changing transcripts in response to Genotype x Treatment interaction. (XLSX 192 kb)
Overrepresented GO Categories of all significantly changing transcripts with the Genotype x Treatment x Time interaction. (XLSX 130 kb)
kME values of all transcripts in each of the gene modules determined by WGCNA. (XLSX 15417 kb)
Summary table of WGCNA modules including top hub genes and gene ontologies. (XLSX 50 kb)