Skip to main content

Genome-wide identification of WRKY family genes and their response to cold stress in Vitis vinifera



WRKY transcription factors are one of the largest families of transcriptional regulators in plants. WRKY genes are not only found to play significant roles in biotic and abiotic stress response, but also regulate growth and development. Grapevine (Vitis vinifera) production is largely limited by stressful climate conditions such as cold stress and the role of WRKY genes in the survival of grapevine under these conditions remains unknown.


We identified a total of 59 VvWRKYs from the V. vinifera genome, belonging to four subgroups according to conserved WRKY domains and zinc-finger structure. The majority of VvWRKYs were expressed in more than one tissue among the 7 tissues examined which included young leaves, mature leaves, tendril, stem apex, root, young fruits and ripe fruits. Publicly available microarray data suggested that a subset of VvWRKYs was activated in response to diverse stresses. Quantitative real-time PCR (qRT-PCR) results demonstrated that the expression levels of 36 VvWRKYs are changed following cold exposure. Comparative analysis was performed on data from publicly available microarray experiments, previous global transcriptome analysis studies, and qRT-PCR. We identified 15 VvWRKYs in at least two of these databases which may relate to cold stress. Among them, the transcription of three genes can be induced by exogenous ABA application, suggesting that they can be involved in an ABA-dependent signaling pathway in response to cold stress.


We identified 59 VvWRKYs from the V. vinifera genome and 15 of them showed cold stress-induced expression patterns. These genes represented candidate genes for future functional analysis of VvWRKYs involved in the low temperature-related signal pathways in grape.


Plants have a variety of defense mechanisms to protect themselves from adverse environmental effects. Families of transcription factors are involved in these processes by functioning to reorganize gene expression patterns. The WRKY family is among them and plays key roles in modulating genes expression during plant defense in response to pathogens [1, 2]. The WRKY transcription factors were first identified in sweet potato (SPF1) as DNA binding proteins [3]. Two similar genes (ABF1 and ABF2) were found in wheat during germination [4]. Subsequently, Rushton et al. [5] reported the identification and characterization of WRKY1, WRKY2 and WRKY3 from parsley (Petroselinum crispum) and proposed these genes belong to a gene family. This gene family was named WRKY due to a conserved region (WRKYGQK) that was identified in the N-terminal amino acid sequence of all the members [4, 5]. Further studies showed that the conserved WRKY domain had other forms such as WRKYGKK and WRKYGEK [6], or the WRKY domain could be replaced by WKKY, WKRY, WSKY, WIKY, WRIC, WRMC, WRRY or WVKY [7, 8].

According to variation in WRKY domain and a zinc finger motif in the C-terminus, WRKY proteins were divided into four groups [9, 10]. WRKY proteins with two WRKY domains composed group I. Groups II and III were characterized by a single WRKY domain. Group II WRKY proteins were further subdivided into five or more subgroups based on short conserved structural motifs while group III proteins contained a variant zinc-finger which ends with HXC. Finally, group IV WRKY proteins contained the WRKY domain, but lack a complete zinc-finger structure in the C-terminus. WRKY proteins usually functioned as transcriptional regulators via binding to W-boxes (TTGACC/T) in the promoter regions of down-stream genes and clusters of W-boxes had an amplified effect [35, 1115]. However, some other studies have found that some WRKY proteins bind to the PRE4 element (TGCGCTT), SURE element (TAAAGATTACTAATAGGAA) or SURE-like element and the WK box (TTTTCCAC) [2].

WRKY proteins have been found to play essential roles in pathogen defense in response to bacteria [16, 17], fungi [18, 19], and viruses [20, 21]. Evidence also supported that WRKY transcription factors were involved in modulating gene expression in plants during abiotic stresses such as cold [22, 23], salt [24, 25] and drought [2628]. Besides roles in response to biotic and abiotic stress, WRKY proteins were also implicated in processes that modulate plant developmental processes such as morphogenesis of trichomes and embryos, senescence, dormancy, and metabolic pathways [2].

Grape is one of the most important fruit crops worldwide. The productivity of grapevines is largely limited by disease pressure and stressful fluctuations in environmental conditions. Due to their essential role in the early response to pathogens and abiotic stresses, several WRKY genes were intensively studied in grape. VvWRKY1 and VvWRKY2, isolated from grape (V. vinifera cv. Cabernet Sauvignon) berries, were found to potentially participate in defending against fungal pathogens [18, 29]. VvWRKY1 was found involving in enhanced protection against Botrytis cinerea by transactivating the VvLTP1 promoter [30], and VvWRKY2 may regulate lignification and response to biotic or abiotic stresses in grapevine [31]. VpWRKY1 and VpWRKY2, isolated from Chinese wild V. pseudoreticulata, may contribute to resistance to powdery mildew (Erysiphe necator) and tolerance to salt and cold stresses in grape [32]. VpWRKY3 was found to be involved in pathogen defense and also interact with the salicylic acid, ethylene, and abscisic acid signal pathways [33]. Transgenetic Arabidopsis plants expressing VvWRKY11, isolated from ‘Beifeng’, an interspecific cultivar of V. thunbergii × V. vinifera, showed increased dehydration tolerance [34]. Its homologous gene, VpWRKY11, was found to serve as a negative regulator of disease resistance [35]. Although several individual WRKY genes have been identified in grapevine, the WRKY gene family in grapevine remains wholly uncharacterized.

Based on our previous transriptome analysis, we found that some WRKY genes respond to cold stress in different patterns in V. amurensis (a cold hardy grapevine species) and V. vinifera cv. Muscat Hamburg [36]. VvWRKY14 (GSVIVT01015952001) and VvWRKY12 (GSVIVT01012682001) were found up-regulated over 30 fold in V. amurensis after being subjected to cold stress but up-regulated to a lesser extent in V. vinifera. In contrast, the expression of VvWRKY43 (GSVIVT01030258001) was up-regulated in V. vinifera (26 fold) while expression remained low in V. amurensis. These different gene expression patterns in response to cold stress may be contribute to the distinctive cold hardiness between the two species. To further characterize how WRKY genes respond to freezing stress of grapevine, we initiated this study to identify the entire WRKY gene family in grapevine based on the published 12× V. vinifera cv. Pinot noir (PN40024) genome sequences [37]. A phylogenetic tree was constructed for identified WRKY proteins and the gene expression patterns in different tissues of V. vinifera were detected by RT-PCR. WRKY genes responding to biotic and abiotic stresses were cross-evaluated by using public gene-chip databases. Additionally, real time RT-PCR was used to detect the expression level of VvWRKYs under cold treatment and exogenous ABA. A comparative analysis was conducted to identify VvWRKYs that may participate in cold signal transduction pathways in V. vinifera using microarray data in public databases, our previously reported transcriptome data and qRT-PCR analysis conducted in this study.


Identifying of WRKY transcription factors in V. viniferagenome

A total of 64 transcripts in the V. vinifera genome sequence were identified as possible members of the WRKY family. Five transcripts were excluded due to a lack of the conserved WRKY domain in the predicted amino acid sequences. The remaining 59 transcripts were named from VvWRKY1 to VvWRKY59 according to their order in the V. vinifera genomic sequence (Table 1). As for the previously published six WRKY proteins in grapes [18, 2935], each amino acid sequence was downloaded and BLASTp was used to find its corresponding WRKY loci in the V. vinifera genome.

Table 1 Identified WRKY genes in 12× V. vinifera ‘Pinot Noir’ genome

The putative genome location of each VvWRKY in the grape genome was shown in Additional file 1: Figure S1. Fifty-eight of the VvWRKYs could be mapped to 18 of the 19 grape chromosomes, with no VvWRKYs found on chromosome 3. VvWRKY4 was putatively located on the ‘Chromosome Unknown’. WRKY transcription factors were not evenly distributed across the chromosomes of the grape genome. There were most abundant on Chromosome 4 (8 VvWRKYs) and chromosome 7 (7 VvWRKYs) and least abundance on Chromosome 5 and 18 (1 VvWRKY).

Categorization of VvWRKYs basis on conserved WRKY domains

The disposition of structural domains in amino acid sequences is an important clue to analyze the evolution and relationship between highly divergent sequences [38]. The relationships among the 59 WRKY proteins were investigated through constructing phylogenetic trees based on multiple alignments of the predicted amino acid sequences of the WRKY domains. As shown in Figure 1, we classified the 59 VvWRKY proteins into four large groups according to the results of the phylogenetic analyses. The models of conserved amino acid sequences of WRKY domain and zinc-finger structure in four groups were shown in Additional file 2: Figure S2. Twelve of the WRKY proteins contained two complete WRKY domains and a C2H2-type zinc finger motif. These proteins constituted group I. The N-terminal WRKY domain (NTWD) and C-terminal WRKY domain (CTWD) of VvWRKY27, VvWRKY41 and VvWRKY56 were clustered into a same clade in group I. According to Eulgem et al. [9] and by using WRKY proteins in Arabidopsis as references, 39 VvWRKY in group II were categorized into five subgroups. Three members were found in subgroup IIa, 8 in IIb, 16 in IIc, 6 in IId and 6 in IIe. Group II was divided into two parts. Subgroup IIa, IIb and IIc showed a close relationship with Group III WRKY proteins. And subgroups IId and IIe belonged to a separate clade which was closely related to group IV. Subgroup IIc showed higher divergence than the other subgroups. There were also 6 WRKY proteins in group III, and 2 in group IV which lacked a complete zinc-finger structure.

Figure 1

Phylogenetic tree of VvWRKY s. The unrooted phylogenetic tree of WRKY domains was constructed with MEGA5.1 program with the neighbor-joining method. The numbers beside the branches represent bootstrap values based on 1000 replications. The name of groups (I, II, III and IV) and subgroup (a–e) were shown at the outside of the circle. The WRKY named with suffix -N or -C indicated the N-terminal WRKY domain (NTWD) or the C-terminal WRKY domain (CTWD) in one VvWRKY with two WRKY domains. AtWRKYs were used as reference to categorize VvWRKYs.

RT-PCR based transcription levels detection of VvWRKYsin different tissues

To investigate if the putative VvWRKYs were expressed and assess their transcription levels in grape, we examined the expression of these genes in different grape tissues. Among all VvWRKYs, we successfully designed and verified 58 primer pairs representing all candidate VvWRKYs except for VvWRKY38 (Figure 2). All transcripts can be detected at least in one tissue. Nineteen VvWRKYs (including VvWRKY02, 11, 12, 13, 14, 17, 20, 24, 28, 30, 33, 34, 35, 36, 39, 41, 42, 48 and 52) were found expressed in all tissues used. Six VvWRKYs (VvWRKY05, 09, 22, 40, 44 and 58) were found only expressed in young tissues. VvWRKY05 was expressed in the stem apex and young fruit. VvWRKY40 was found in stem apex, young fruit and root. VvWRKY09, 22, 44 and 58 were detected in young leaf, stem apex, young fruit and root.

Figure 2

RT-PCR analyses of presence of VvWRKY transcripts in seven grape tissues. YL: young leaf; ML: mature leaf; T: tendril; S: stem apex; YF: young fruit; RF: ripe fruit; R: root. VvMDH and VvACT were used as control.

Gene-chip based expression analysis of 26 VvWRKYsunder various stresses

Although we identified WRKY transcription factors from the V. vinifera genome, functions for these genes in response to abiotic and biotic stress remain unknown. Using microarray results from publically available data, it was possible to find gene expression data from multiple experimental conditions for several of the grapevine WRKY genes. We carefully checked the genes on the ‘GeneChip Vitis vinifera (Grape) Genome Array’ (Affymetrix) and a total of 26 VvWRKYs were found on this chip. Microarray data related to salinity, water-deficit, PEG, cold, ABA and pathogen stresses were downloaded and their corresponding probes and the CV (coefficent of variation of the corresponding treatment means) of these genes in each of the microarray experiments were listed in Additional file 3: Table S1. If the expression of a probe set (gene) is affected by some of the treatments in an experiment, it shows a higher CV (more fluctuation); and vice versa. According to the data, the CV of 20 of the 26 VvWRKYs were over 5% in at least one experiment. The highest CV appeared in VvWRKY57 (up to 36%) associated with compatible viral diseases in berry experiment in V. vinifera cv. Cabernet Sauvignon. VvWRKY03, 06, 08, 28, and 55 responded to both abiotic and pathogens stresses while VvWRKY21, 39, 48 seemed to respond primarily to pathogens stresses.

To test the correlation between the expression patterns of 26 VvWRKYs and their phylogenetic relationship, a hierarchical cluster analysis was performed using the 11 stress related experimental datasets (Figure 3). Red, black and green elements in the matrix indicate up-, no change- and down-regulated expression of WRKY transcription factors, respectively. From the heat map, twenty-six genes were clustered into four clades. Carefully analyzing the cluster of expression data in response to abiotic stresses experiments and comparing this with the VvWRKYs phylogenetic tree, we found that genes with close phylogenetic relationship were classified into the same clade during hierarchical cluster analysis. The most obvious evidence can be found in clades 3 with 5 WRKY subgroup IId genes (including VvWRKY07, 24, 39, 48 and 57), which show similar expression patterns in response to salt, PEG and cold stresses. Clade 1 contained three WRKY group I genes and two group IIC genes. Clade 2 was mainly composed by WRKY group I and IIC and contains a majority of cold stress-related VvWRKYs (Also shown in Additional file 4: Table S2). Clade 4 only had one gene and that gene was from WRKY group III.

Figure 3

Cluster analyses of VvWRKY s from 16 k Affymetrix V. vinifera gene-chip data in PLEXdb database. The relative expression values of 26 VvWRKYs responding to different abiotic stresses (salinity, water deficit, cold) were used in analysis. Red, black and green elements in the matrix indicate up-regulated, no change and down-regulated WRKY genes, respectively. Those genes can be classified into four groups according to expression patterns, which were shown in different color with its group IDs that coincide with Figure 1. The red color was used to emphasize the VvWRKYs that changed in expression over 2 fold under cold treatment.

Real-time RT-PCR based expression analysis of VvWRKYs under cold treatment in V. vinifera

To examine the response of VvWRKYs under cold stress in grape, we examined the transcription levels of VvWRKYs in shoot apices of ‘Muscat Hamburg’ under cold-treatment (4°C). VvWRKY05, 21, 32 and 40 were excluded from cold-treated experiment since their Ct value of amplification curve were over 35 cycles in the templates of normal and cold-treated shoot apex. Detected VvWRKYs can be classified into four groups according to expression patterns as shown in Figure 4 and Additional file 5: Figure S3: A) sustained up-regulated during cold treatment (22 genes, Figure 4A), B) changed above 2 fold with irregular pattern (9 genes, Figure 4B), C) sustained down-regulated (5 genes, Figure 4C) and D) no significant difference (18 genes, as shown in Additional file 5: Figure S3). The relative expressions of 36 genes (Figure 4A, B and C) were significantly different as cold treatment. The greatest increase in expression (nearly 30 fold) was found in VvWRKY55 at 48 h cold treatment. VvWRKY18 and VvWRKY46 had the largest up-regulation of greater than 6 fold at 8 hours after cold treatment. While VvWRKY18 was degraded after 24 hours, the expression of VvWRKY46 demonstrated both up and down regulated with a spike of expression at 48 hours after intensive degradation at 24 hours.

Figure 4

qRT-PCR assays of the expression patterns of VvWRKY s under cold and exogenous ABA treatments. The default expression value for each gene was 1 at 0 hours before treatment. A, B and C represent the subgroups with different expression patterns in cold treatment and D represents the genes that up-regulated over 2 fold in ABA treatment. A: sustained up-regulated genes in cold treatment; B: genes that changes over 2 fold but without significant tendency in cold treatment; C: sustained down-regulated genes in cold treatment; D: up-regulated genes in exogenous ABA treatment. VvWRKYs that accumulated in both cold and exogenous ABA treatments were underlined. One-Way ANOVA analysis was used to test the impact of timing of cold treatment. When the effects were significantly different, we examined the difference between treatments using post hoc multiple comparisons (LSD, p < 0.05). All data analyses were conducted using IBM SPSS Statistics 20, and the results were displayed through a, b, c and d.

Exogenous ABA induced accumulation of VvWRKYs in V. Vinifera

To illustrate how the VvWRKYs respond to ABA and whether the cold stress related VvWRKYs may participate in the ABA-dependent cold signal pathway, ABA treated grapevine apices were examined using qRT-PCR. VvWRKY12, 29, and 46 were excluded from this experiment due to their higher Ct value (Figure 4D and Additional file 6: Figure S4). Among the 55 VvWRKYs we detected, twelve VvWRKYs were expressed over 2-fold greater within 2 h of exogenous ABA treatment (Figure 4D). After statistical analyses of qRT-PCR results, 7 of them were evaluated to significantly change during exogenous ABA treatments. Transcripts of VvWRKY35 showed the greatest increase in expression at 0.5 h after ABA treatments. Six other genes showed increases in expression 1 h after exogenous ABA treatment (Figure 4D).

When the data from the cold and ABA experiments were compared, 6 of 7 genes (VvWRKY, 19, 28, 35, 42, 50 and 55) that were up-regulated during exogenous ABA treatment were also up-regulated under cold treatment (Figure 4A and B, marked by underline). Two genes (VvWRKY55, 28) were greatly up-regulated, over 10 fold. The expression levels of the rest of the 44 VvWRKYs were lower than 2-fold and not significantly changed during exogenous ABA treatments (Additional file 6: Figure S4).

Identification of candidate cold-stress related VvWRKYs

Previously we reported the changes of the transcriptome during cold-treatments in ‘Muscat Hamburg’ and identified 14 cold-stress related VvWRKYs (we reported 16 VvWRKYs but subsequent annotation of these genes allowed us to exclude two genes that do not belong to the WRKY gene family)[36]. Gene-chip based methods also allowed to identify 10 cold-stress related VvWRKYs[39]. In order to overcome the deficiencies of determining gene expression from a single technological approach and obtain more reliable results, we compared the data from three different methods. Fourteen VvWRKYs from our previous transcriptome analysis, ten from publically available gene-chip based data and 36 genes from qRT-PCR results (this study) were used. The results were summarized in Figure 5 and Additional file 4: Table S2. Three VvWRKYs (VvWRKY12, 28, 55) showed identical expression patterns and were found up-regulated over 10 fold in at least one time-point under cold-treatment by qRT-PCR (Figure 4A). A total of 12 VvWRKYs were confirmed by two experimental methods (Figure 5A and B). VvWRKY56 was identified as up-regulated gene under cold treatment only in the gene-chip studies. Twenty-two genes that were characterized by qRT-PCR were not supported by the other studies. It is worth mentioning that down-regulated VvWRKYs under cold-treatment were only identified by qRT-PCR based method.

Figure 5

An overview of cold stress-related VvWRKYs in three sets of data. A: The Venn diagram of the cold stress-related VvWRKYs obtained from qRT-PCR, transcriptome and gene-chip data. B: The VvWRKYs that were found in more than one type of experimental data. Green color in forms indicated VvWRKYs induced by exogenous ABA.


WRKY family in grape

Considering the important roles that WRKY transcription factors play during plant development and in response to various stresses, it is not surprising that we identified so many family members in grapevine. Previously, 74 WRKY genes were found in Arabidopsis[2], 55 in cucumber [40], 102 in rice [2], 47 in castor bean [41], 86 in Brachypodium distachyon[42] and 136 in maize [43]. Here we identified 59 candidate WRKY proteins in V. vinifera and categorized them into four groups.

Group I WRKY proteins

When compared with WRKY family groups, WRKYs in primitive plant ancestors Giatdia lamblia, Dictyostelium discoideum and Chlamydomonas reinhardtii closely resembled Vitis group I [7, 38]. In our study, two domains of VvWRKYs in group 1 were closely related. A BLASTp search of EuGene.1100010359 from an ancient alga species (Ostreococcus sp. RCC809) which has a single WRKY domain allowed us to identify 9 corresponding WRKY homologs in grape and 8 of these belonged to group I by MAP VIEW (Plant Genome Duplication Database) [44]. These data support the hypothesis that the dual WRKY domains present in members of group I may be derived from a single WRKY domain duplication [6, 7].

Group II WRKY proteins

Group II was divided into three parts: subgroup IIa + IIb, subgroup IIc and subgroup IId + IIe (Figure 1). Subgroup IIa + IIb belong to the same clade and is sister to the WRKYs in group I. Interestingly, the presumed function of CTWDs in group I for sequence-specific DNA binding [9] were more similar to the single WRKY domain members in group II and III than to the NTWDs of group I, This result may indicate that subgroup IIa + IIb evolved from group I by domain structure loss of the group I NTWD.

Group III WRKY proteins

Group III in the phylogenetic tree was most closely related to the very large subgroup IIc, which was separately into four clades and seemed to indicate an expansion of the gene family. A thorough search of the Plant Transcription Factor Database ( indicated that the earliest evolutionary occurrences of group III genes were those found in ferns (Selaginella moellendorffii). There was no evidence of any sequenced plant species that only contain members of group I and III but we found in some species with only members of group I and II, for example in mosses (Physcomitrella patens) [1], and some gymnosperms (Pinus taeda) ( We speculated that group III may have evolved from group II, particularly IIc. As group III WRKYs in Arabidopsis responded to diverse biotic stresses [45], group III members may indicate adaptation of early plants to the stressful conditons associated with the colonization of land and subsequent increase in biotic pathogen pressures.

Group IV WRKY proteins

We found that group IV WRKY proteins, which were characterized by the loss of the zinc-finger domain, were in the same clade as subgroups IId + IIe. VvWRKY02 and VvWRKY57 were duplicated gene pairs according to a whole genome analysis of grapevine gene duplications [46]. This might suggest an origin of group IV from subgroups IId + IIe. Group IV proteins were considered non-functional due to the loss of the zinc-finger domain [10]. However, these genes of group IV can be found in all higher plant species as well as in algaes (Bathycoccus prasinos: Bathy17g02050). Furthermore, some genes were expressed in rice (OsWRKY56) [10] as were two genes identified in this study (VvWRKY02 and VvWRKY29). Therefore, it remains questionable whether group IV WRKYs have biological function in plants.

VvWRKYsparticipate in development and stress-related signal pathways

WRKY genes were found to be expressed in many tissues and seem to be involved in regulating plant developmental and physiological processes. Transcriptomic analysis of senescence in the flag leaf of wheat demonstrated that WRKY transcription factors are greatly up-regulated during the senescence process [47]. OsWRKY78 was found to be up-regulated in elongating stems and knockdown mutations in this gene cause plants to produce a semi-dwarf and small seed phenotype caused by reducing cell length [48]. Moreover, the transcription of GhWRKY15 was observed abundant in the roots and stems of tobacco and transgenic overexpression lines of these plants displayed faster elongation at the earlier shooting stages [49]. Here the expression of 15 VvWRKYs (Figure 2) can be detected in all grape tissues we used, which may indicate its fundamental roles in different cell-types in grape. Similar to expression patterns observed in other plant species, VvWRKYs were found to be expressed in young tissues such as young leaf, shoot apex, tendril and young fruit.

Several numbers of VvWRKYs were found activated in more than one type of stress condition (Figure 3 and also Additional file 3: Table S1). VpWRKY3, homologous to VvWRKY55 was observed to be up-regulated in response to many different sources of stress, including pathogen exposure, salicylic acid, ethylene, cold and drought stress [32]. VvWRKYs that were up-regulated in response to more than two types of stresses (e.g. pathogen and drought) supported the occurrence of cross-talk between signal transduction pathways in response to different stress conditions in plants [50].

Phylogenetic relationships between VvWRKY genes suggested that there may be conserved responses of these genes to salt exposure, PEG and cold-stress (Figure 3). All members of group IId clustered into one clade with similar expression pattern during these three stress conditions, suggesting the function of these VvWRKY proteins may relate to the structures of WRKY domains. Subgroup IId was identified as a novel CaM-binding transcription factor family in plants and their conserved structural motif was a Ca2+-dependent CaM-binding domain [51]. Thus the placement of the WRKYs in the phylogenetic tree may also help to predict function of new members that belong to certain gene family.

VvWRKYsthat participate in the cold related signal transduction in grape

Three different experimental methods were combined to robustly analyze the response of VvWRKY genes to cold stress (Figure 5 and Additional file 4: Table S2). Results from qRT-PCR demonstrated the greatest number of cold stress-related VvWRKYs (36) while gene-chip based methods identified the least, 10 VvWRKYs. This difference may be attributed to the method used but is also likely due to differences in the treatment conditions between experiments. During Digital Gene Expression profile (DGE) analysis [36], plant material was obtained from 4 h cold treatment at 4°C, whereas in our pRT-PCR experiment, we used samples collected at several different time periods (at 8 h, 24 h and 48 h after cold treatment at 4°C). Additionally, multiple matched tags were excluded from the final analysis performed by Xin et al. [36], which may have reduced the number of identified cold related VvWRKYs. Finally, gene-chip based methods may bias results due to a lower number of genes with corresponding probes related to the WRKY proteins (only 26 WRKY). By integrating the data from different methods, we obtained more reliable results and a total of 15 candidate cold tolerance VvWRKYs (Figure 5) were identified during our investigation.

According to previous studies, the transcriptional control of plant responses to cold stress can be divided into ABA-dependent and ABA-independent signal pathways [52]. The results of our study also indicated that 15 putative cold stress-related VvWRKYs can be divided into two groups according to their responses to exogenous ABA. Three VvWRKYs (VvWRKY28, 42 and 55) may participate in an ABA-dependent signal pathway and other 12 in ABA-independent pathway. WRKY transcription factors have been identified as key components in the ABA signaling pathways [8, 53]. In rice, OsWRKY24, 51, 71 and 72 are induced by (ABA) in aleurone cells. OsWRKY24 and 45 were functional as negative regulators in ABA induction of the HVA22 promoter-beta-glucuronidase construct, while OsWRKY72 and 77 synergistically interacted with ABA to activate this reporter construct [10]. It is still unknown how WRKYs participate in the cold stress-related signal pathway and what relationship these genes have with C-repeat Binding Factor genes (CBFs), which are critical transcription factors responsible for cold tolerance in plant [54].

The reliability of the identified 15 cold–related VvWRKYs was also supported by homologous genes in other species. STHP-64, which showed high similarity with VvWRKY43, was not present in leaves until November and December in Solanum dulcamara[55]. WRKY38, a homolog gene of VvWRKY14, was transiently accumulated when leaves and roots were exposure to low temperature in barley [56]. BcWRKY46 showed higher similarity with VvWRKY33 and responded to low temperatures in Pak-choi. Constitutive expression of BcWRKY46 reduced the freezing susceptibility in transgenic tobacco [57]. The transcription level of VvWRKY55 was up-regulated robust under cold treatment. Its homolog gene, WRKY71 was found in banana with a similar expression pattern [58]. All these VvWRKYs mentioned above were confirmed by at least two set of experiment methods, which provided appropriate candidates to illustrate the roles of WRKY protein under low temperature-related signal pathways in grape.

Although low-temperature related WRKYs were isolated in several species, the mechanism of how WRKYs respond to cold signals and regulate the expression of downstream genes is still largely unknown. Further work is needed to elucidate the function of these important genes in low-temperature related signal pathways. Previously we reported the different expression patterns of WRKYs in V. amurensis, a cold-hardness species. The WRKY genes identified here from V. vinifera may accelerate the functional analysis of this gene family in V. amurensis. The comprehensive analysis of cold stress-related WRKYs in two different Vitis species with contrasting cold hardiness phenotypes would certainly help to illustrate the function of WRKY genes in conveying cold hardiness in grapevine.


In summary, a total of 59 VvWRKYs in the V. vinifera genome were identified. The VvWRKYs were unevenly distributed in 18 of the 19 chromosomes. WRKY domain based phylogenetic analysis allowed categorizing 59 VvWRKYs into four large groups. A majority of VvWRKYs were found expressed in more than one tissue in V. vinifera. Gene-chip based data analysis suggested that a subset of VvWRKYs was activated in respond to diverse biotic and abiotic stresses. The transcription level of 36 VvWRKY genes changed over 2 fold after cold induction. A comparative analysis of qRT-PCR results, gene-chip based data and transcriptome analysis allowed us to identify 15 VvWRKYs that show identical expression patterns during cold treatment at least in two kinds of analyses. These studies not only increase our knowledge of WRKY family, but also provide candidate genes for future functional analysis of VvWRKYs involved in the low temperature-related signal pathways in grape.


Identification of WRKYgenes in the grape genome

Candidate WRKY proteins were identified from the 12X V. vinifera cv. Pinot noir genome (quasi-homozygous line PN40024, Full-length amino acids sequences of all WRKY proteins in Arabidopsis thaliana ( were used as query sequences. A BLASTp search was performed and E-value of e−6 was used as the threshold [59]. Candidate WRKY proteins were manually confirmed [60] by searching for WRKY domains in the candidate amino acids sequences using SWISS-MODEL ( and the results were shown in Table 1.

Phylogenetic analysis of WRKY family

Multiple alignments of the amino acid sequences of 73 WRKY domains from V. vinifera were performed using CLUSTALW by MEGA5.1 [61]. Twelve Arabidopsis WRKY domains from different WRKY groups were used as references to categorize the WRKY proteins from grape. The GenBank accession numbers of those AtWRKYs are AtWRKY01: ABJ17102, AtWRKY11: AEE85928.1, AtWRKY14: AAP21276.1, AtWRKY18: AAM78067, AtWRKY21: AAB63078.1, AtWRKY27: ABH04558, AtWRKY28: AEE84006, AtWRKY31: AEE84546.1, AtWRKY38: AED93044.1, AtWRKY41: AEE82969, AtWRKY43: AEC10646.1, AtWRKY45: ABD57509.1, AtWRKY49: AAQ62425.1. The parameters used during alignment were: protein weight matrix: Gonnet series; negative matrix: on; gap open penalty: 10; gap extension penalty: 0.20; delay divergent sequences: 30; residue-specific gap penalties: on; hydrophilic penalties: on; gap separation distance: 0; end gap separation penalty: on. An unrooted phylogenetic tree was constructed using Neighbor-Joining (NJ) methods and bootstrapped with 1,000 iterations to help identify WRKY protein groups.

Plant materials

‘Muscat Hamburg’ (V. vinifera) was obtained from the Institute of Botany, the Chinese Academy of Sciences. Tissues of young leaves, mature leaves, tendril, stem apex, root, young fruits and ripe fruits were collected from the vineyard in July, 2012. Cold and exogenous ABA treatment experiments were performed on tissue culture seedlings of ‘Muscat Hamburg’ according to Li et al. [62]. Briefly, seedlings were cultured on 1/2 B5 medium with 30 g/L of sucrose in a growth chamber under 16-h light/8-h dark photoperiod at 26°C. Cold treatments were performed in another growth chamber with the same parameters except for temperature (4°C). Seedlings with five well developed leaves were used and the shoot apex with one well developed leaf was collected at 0 hour (h, used as control), 8 h, 24 h and 48 h. Seedlings with five well developed leaves were transplanted in 1/2 B5 nutrient solution. Exogenous ABA treatments were performed after one week under normal culture conditions with 100 μM ABA and the shoot apex with one well developed leaves were collected at 0 h (used as control), 0.5 h, 1 h and 2 h after treatments. Three independent replicates were collected for each time point and frozen in liquid nitrogen. Samples were then stored at - 80°C for the following RNA isolation.

Expression patterns analysis of VvWRKYsby RT-PCR

Total RNA was isolated from collected samples using Plant Total RNA Isolation kit (Tiandz Inc; Beijing, China). RNase-free DNase (RQ1, Promega) was used to degrade DNA from total RNA. cDNA was synthesized by the SuperScript III Reverse Transcriptase (Invitrogen) with Oligo(dT)18 (Promega) according to the manufacturer’s instructions. Primer pairs (Additional file 7: Table S3) for VvWRKYs were designed by Primer 3 ( and tested by NCBI Primer BLAST. Two genes, β-actin (GenBank accession: EC969944; sense primer: 5′-CTTGCATCCCTCAGCACCTT-3′; antisense primer: 5′-TCCTGTGGACAATGGATGGA-3′) and malate dehydrogenase gene (MDH; GenBank accession: EC921711; sense primer: 5′-CCATGCATCACCCACAA-3′; antisense primer: 5′-GTCAACCATGCTACTGTCAAAACC-3′) were used as positive control for RT-CR [63]. Three biological replicate and 35 cycles for each reaction were performed. PCR products were detected by agarose gel electrophoresis with 2.5% gel concentration.

Gene-chip based expression pattern analysis of VvWRKYs

We explored the expression profiles of VvWRKYs using publically available data from the 16 k Affymetrix V. vinifera gene-chip stored at PLEXdb (Plant Expression Database) [64] to explore the response of VvWRKYs during biotic and abiotic stresses in grape. The different studies and datasets that were included in these analyses were: A) a short term abiotic stress experiment in ‘Cabernet Sauvignon’ [39], B) a long-term salt and water stress study [65]; C) a study examining gene expression associated with compatible viral diseases in grapevine cultivars [66]; D) an experiment designed to examine the powdery mildew-induced transcriptome in a susceptible grapevine ‘Cabernet Sauvignon’ [67]; E) the complimentary dataset of the powdery mildew-induced transcriptome of a resistant grapevine ‘Norton’ [67]; F) a study of gene expression in grapevine in response to Bois noir infection [68]; G) a study of the grape skin transcriptome of berries grown on an exogenous abscisic acid treated vine [69]; H) the complimentary dataset of the grape skin transcriptome in the berries cultured in vitro and treated with exogenous ABA [69]; and lastly, I) a gene expression study associated with compatible viral diseases in the berry [70]. In our comparative analysis, we divided these experiments into either abiotic or biotic stresses related datasets. For each microarray experiment, the Affymetrix MAS5.0 normalized data were used for calculations of the fold change of differentially expressed genes. Probe sets corresponding to the putative VvWRKYs were identified at PLANEX ( and completed via PLEXdb blast tool. Comparisons of WRKY expression level from gene-chip data for the short term abiotic stress treatment in ‘Cabernet Sauvignon’ was performed using Cluster 3.0 and JavaTreeview.

Quantitative RT-PCR

Total RNA was isolated from cold and exogenous ABA treated shoot apices following the cDNA synthesis methods mentioned above. Synthesized cDNA was diluted 1:10 with ddH2O, and the quantitative RT-PCR reaction mixture contained 5 μl of 2 × SYBR Green I Master Mix (Roche, USA), 2.6 μL ddH2O, 0.2 μL of 10 μM solution of each primer and 2 μL diluted template cDNA. Reaction specificities for each primer pair was tested using qRT-PCR melting curve analysis. The experiment was carried out using a StepOnePlus real-time PCR Instrument (Applied Biosystems). Transcription levels of each VvWRKY was normalized against the average of β-actin, MDH (as mentioned above) and glyceraldehyde-3-phosphatedehydrogenase (GAPDH: CB973647; sense primer: 5′-TTCTCGTTGAGGGCTATTCCA-3′; antisense primer: 5′-CCACAGACTTCATCGGTGACA-3′) [63]. Each sample had three biological and two technical replicates to ensure the accuracy of results, and RNA samples with the same reverse-transcription (without Reverse Transcriptase) and dilution procedure were used as negative controls. The Ct values and the real-time PCR efficiencies were obtained using Lin-RegPCR [71] and the normalized relative quantities and standard errors for each sample were calculated by qbaseplus [72]. The relative expression level of each VvWRKY in different templates was calculated based on normalized relative quantities. We used One-Way ANOVA analysis to test the impact of timing of cold treatment. When the effects were significantly different, we examined the difference between treatments using post hoc multiple comparisons (LSD, p < 0.05). All data analyses were conducted using IBM SPSS Statistics 20.



Abscisic acid.


  1. 1.

    Ulker B, Somssich IE: WRKY transcription factors: from DNA binding towards biological function. Curr Opin Plant Biol. 2004, 7 (5): 491-498. 10.1016/j.pbi.2004.07.012.

    Article  PubMed  Google Scholar 

  2. 2.

    Rushton PJ, Somssich IE, Ringler P, Shen QXJ: WRKY transcription factors. Trends Plant Sci. 2010, 15 (5): 247-258. 10.1016/j.tplants.2010.02.006.

    Article  CAS  PubMed  Google Scholar 

  3. 3.

    Ishiguro S, Nakamura K: Characterization of a cDNA encoding a novel DNA-binding protein, SPF1, that recognizes SP8 sequences in the 5′ upstream regions of genes coding for sporamin and beta-amylase from sweet potato. Mol Gen Genet. 1994, 244 (6): 563-571.

    Article  CAS  PubMed  Google Scholar 

  4. 4.

    Rushton PJ, Macdonald H, Huttly AK, Lazarus CM, Hooley R: Members of a new family of DNA-binding proteins bind to a conserved cis-element in the promoters of alpha-Amy2 genes. Plant Mol Biol. 1995, 29 (4): 691-702. 10.1007/BF00041160.

    Article  CAS  PubMed  Google Scholar 

  5. 5.

    Rushton PJ, Torres JT, Parniske M, Wernert P, Hahlbrock K, Somssich IE: Interaction of elicitor-induced DNA-binding proteins with elicitor response elements in the promoters of parsley PR1 genes. Embo J. 1996, 15 (20): 5690-5700.

    PubMed Central  CAS  PubMed  Google Scholar 

  6. 6.

    Ross CA, Liu Y, Shen QXJ: The WRKY gene family in rice (Oryza sativa). J Integr Plant Biol. 2007, 49 (6): 827-842. 10.1111/j.1744-7909.2007.00504.x.

    Article  CAS  Google Scholar 

  7. 7.

    Zhang YJ, Wang LJ: The WRKY transcription factor superfamily: its origin in eukaryotes and expansion in plants. BMC Evol Biol. 2005, 5: 1-

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  8. 8.

    Rushton DL, Tripathi P, Rabara RC, Lin J, Ringler P, Boken AK, Langum TJ, Smidt L, Boomsma DD, Emme NJ, Chen X, Finer JJ, Shen QJ, Rushton PJ: WRKY transcription factors: key components in abscisic acid signalling. Plant Biotechnol J. 2012, 10 (1): 2-11. 10.1111/j.1467-7652.2011.00634.x.

    Article  CAS  PubMed  Google Scholar 

  9. 9.

    Eulgem T, Rushton PJ, Robatzek S, Somssich IE: The WRKY superfamily of plant transcription factors. Trends Plant Sci. 2000, 5 (5): 199-206. 10.1016/S1360-1385(00)01600-9.

    Article  CAS  PubMed  Google Scholar 

  10. 10.

    Xie Z, Zhang ZL, Zou XL, Huang J, Ruas P, Thompson D, Shen QJ: Annotations and functional analyses of the rice WRKY gene superfamily reveal positive and negative regulators of abscisic acid signaling in aleurone cells. Plant Physiol. 2005, 137 (1): 176-189. 10.1104/pp.104.054312.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  11. 11.

    Eulgem T, Rushton PJ, Schmelzer E, Hahlbrock K, Somssich IE: Early nuclear events in plant defence signalling: rapid gene activation by WRKY transcription factors. Embo J. 1999, 18 (17): 4689-4699. 10.1093/emboj/18.17.4689.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  12. 12.

    Laloi C, Mestres-Ortega D, Marco Y, Meyer Y, Reichheld JP: The Arabidopsis cytosolic thioredoxin h5 gene induction by oxidative stress and its W-box-mediated response to pathogen elicitor. Plant Physiol. 2004, 134 (3): 1006-1016. 10.1104/pp.103.035782.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  13. 13.

    Yang PZ, Chen ZX: A family of dispersed repetitive DNA sequences in tobacco contain clusters of W-box elements recognized by pathogen-induced WRKY DNA-binding proteins. Plant Sci. 2001, 161 (4): 655-664. 10.1016/S0168-9452(01)00454-X.

    Article  CAS  Google Scholar 

  14. 14.

    Yamamoto S, Nakano T, Suzuki K, Shinshi H: Elicitor-induced activation of transcription via W box-related cis-acting elements from a basic chitinase gene by WRKY transcription factors in tobacco. Bba-Gene Struct Expr. 2004, 1679 (3): 279-287. 10.1016/j.bbaexp.2004.07.005.

    Article  CAS  Google Scholar 

  15. 15.

    Wang Z, Zhu Y, Wang LL, Liu X, Liu YX, Phillips J, Deng X: A WRKY transcription factor participates in dehydration tolerance in Boea hygrometrica by binding to the W-box elements of the galactinol synthase (BhGolS1) promoter. Planta. 2009, 230 (6): 1155-1166. 10.1007/s00425-009-1014-3.

    Article  CAS  PubMed  Google Scholar 

  16. 16.

    Mukhtar MS, Deslandes L, Auriac MC, Marco Y, Somssich IE: The Arabidopsis transcription factor WRKY27 influences wilt disease symptom development caused by Ralstonia solanacearum. Plant J. 2008, 56 (6): 935-947. 10.1111/j.1365-313X.2008.03651.x.

    Article  CAS  PubMed  Google Scholar 

  17. 17.

    Tao Z, Liu HB, Qiu DY, Zhou Y, Li XH, Xu CG, Wang SP: A Pair of Allelic WRKY Genes Play Opposite Roles in Rice-Bacteria Interactions. Plant Physiol. 2009, 151 (2): 936-948. 10.1104/pp.109.145623.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  18. 18.

    Marchive C, Mzid R, Deluc L, Barrieu F, Pirrello J, Gauthier A, Corio-Costet MF, Regad F, Cailleteau B, Hamdi S, Lauvergeat V: Isolation and characterization of a Vitis vinifera transcription factor, VvWRKY1, and its effect on responses to fungal pathogens in transgenic tobacco plants. J Exp Bot. 2007, 58 (8): 1999-2010. 10.1093/jxb/erm062.

    Article  CAS  PubMed  Google Scholar 

  19. 19.

    Gallou A, Declerck S, Cranenbrouck S: Transcriptional regulation of defence genes and involvement of the WRKY transcription factor in arbuscular mycorrhizal potato root colonization. Funct Integr Genomics. 2012, 12 (1): 183-198. 10.1007/s10142-011-0241-4.

    Article  CAS  PubMed  Google Scholar 

  20. 20.

    Yang PZ, Chen CH, Wang ZP, Fan BF, Chen ZX: A pathogen- and salicylic acid-induced WRKY DNA-binding activity recognizes the elicitor response element of the tobacco class I chitinase gene promoter. Plant J. 1999, 18 (2): 141-149. 10.1046/j.1365-313X.1999.00437.x.

    Article  CAS  Google Scholar 

  21. 21.

    Wang YN, Dang FF, Liu ZQ, Wang X, Eulgem T, Lai Y, Yu L, She JJ, Shi YL, Lin JH, Chen CC, Guan DY, Qiu A, He SL: CaWRKY58, encoding a group I WRKY transcription factor of Capsicum annuum, negatively regulates resistance to Ralstonia solanacearum infection. Mol Plant Pathol. 2013, 14 (2): 131-144. 10.1111/j.1364-3703.2012.00836.x.

    Article  PubMed  Google Scholar 

  22. 22.

    Zhou XF, Wang GD, Sutoh K, Zhu JK, Zhang WX: Identification of cold-inducible microRNAs in plants by transcriptome analysis. Bba-Gene Regul Mech. 2008, 1779 (11): 780-788.

    CAS  Google Scholar 

  23. 23.

    Zou CS, Jiang WB, Yu DQ: Male gametophyte-specific WRKY34 transcription factor mediates cold sensitivity of mature pollen in Arabidopsis. J Exp Bot. 2010, 61 (14): 3901-3914. 10.1093/jxb/erq204.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  24. 24.

    Jiang YQ, Deyholos MK: Functional characterization of Arabidopsis NaCl-inducible WRKY25 and WRKY33 transcription factors in abiotic stresses. Plant Mol Biol. 2009, 69 (1–2): 91-105.

    Article  CAS  PubMed  Google Scholar 

  25. 25.

    Wu XL, Shiroto Y, Kishitani S, Ito Y, Toriyama K: Enhanced heat and drought tolerance in transgenic rice seedlings overexpressing OsWRKY11 under the control of HSP101 promoter. Plant Cell Rep. 2009, 28 (1): 21-30. 10.1007/s00299-008-0614-x.

    Article  CAS  PubMed  Google Scholar 

  26. 26.

    Ren XZ, Chen ZZ, Liu Y, Zhang HR, Zhang M, Liu QA, Hong XH, Zhu JK, Gong ZZ: ABO3, a WRKY transcription factor, mediates plant responses to abscisic acid and drought tolerance in Arabidopsis. Plant J. 2010, 63 (3): 417-429. 10.1111/j.1365-313X.2010.04248.x.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  27. 27.

    Jiang YJ, Liang G, Yu DQ: Activated expression of WRKY57 confers drought tolerance in Arabidopsis. Mol Plant. 2012, 5 (6): 1375-1388. 10.1093/mp/sss080.

    Article  CAS  PubMed  Google Scholar 

  28. 28.

    Shen HS, Liu CT, Zhang Y, Meng XP, Zhou X, Chu CC, Wang XP: OsWRKY30 is activated by MAP kinases to confer drought tolerance in rice. Plant Mol Biol. 2012, 80 (3): 241-253. 10.1007/s11103-012-9941-y.

    Article  CAS  PubMed  Google Scholar 

  29. 29.

    Mzid R, Marchive C, Blancard D, Deluc L, Barrieu F, Corio-Costet MF, Drira N, Hamdi S, Lauvergeat V: Overexpression of VvWRKY2 in tobacco enhances broad resistance to necrotrophic fungal pathogens. Physiol Plant. 2007, 131 (3): 434-447. 10.1111/j.1399-3054.2007.00975.x.

    Article  CAS  PubMed  Google Scholar 

  30. 30.

    Laquitaine L, Gomes E, Francois J, Marchive C, Pascal S, Hamdi S, Atanassova R, Delrot S, Coutos-Thevenot P: Molecular basis of ergosterol-induced protection of grape against Botrytis cinerea: Induction of type I LTP promoter activity, WRKY, and stilbene synthase gene expression. Mol Plant Microbe In. 2006, 19 (10): 1103-1112. 10.1094/MPMI-19-1103.

    Article  CAS  Google Scholar 

  31. 31.

    Guillaumie S, Mzid R, Mechin V, Leon C, Hichri I, Destrac-Irvine A, Trossat-Magnin C, Delrot S, Lauvergeat V: The grapevine transcription factor WRKY2 influences the lignin pathway and xylem development in tobacco. Plant Mol Biol. 2010, 72 (1–2): 215-234.

    Article  CAS  PubMed  Google Scholar 

  32. 32.

    Li HE, Xu Y, Xiao Y, Zhu ZG, Xie XQ, Zhao HQ, Wang YJ: Expression and functional analysis of two genes encoding transcription factors, VpWRKY1 and VpWRKY2, isolated from Chinese wild Vitis pseudoreticulata. Planta. 2010, 232 (6): 1325-1337. 10.1007/s00425-010-1258-y.

    Article  CAS  PubMed  Google Scholar 

  33. 33.

    Zhu ZG, Shi JL, Cao JL, He MY, Wang YJ: VpWRKY3, a biotic and abiotic stress-related transcription factor from the Chinese wild Vitis pseudoreticulata. Plant Cell Rep. 2012, 31 (11): 2109-2120. 10.1007/s00299-012-1321-1.

    Article  CAS  PubMed  Google Scholar 

  34. 34.

    Liu HY, Yang WL, Liu DC, Han YP, Zhang AM, Li SH: Ectopic expression of a grapevine transcription factor VvWRKY11 contributes to osmotic stress tolerance in Arabidopsis. Mol Biol Rep. 2011, 38 (1): 417-427. 10.1007/s11033-010-0124-0.

    Article  CAS  PubMed  Google Scholar 

  35. 35.

    Yu Y, Xu W, Wang J, Wang L, Yao W, Yang Y, Xu Y, Ma F, Du Y, Wang Y: The Chinese wild grapevine (Vitis pseudoreticulata) E3 ubiquitin ligase Erysiphe necator-induced RING finger protein 1 (EIRP1) activates plant defense responses by inducing proteolysis of the VpWRKY11 transcription factor. New Phytol. 2013, 200 (3): 834-846. 10.1111/nph.12418.

    Article  CAS  PubMed  Google Scholar 

  36. 36.

    Xin H, Zhu W, Wang L, Xiang Y, Fang L, Li J, Sun X, Wang N, Londo JP, Li S: Genome wide transcriptional profile analysis of Vitis amurensis and Vitis vinifera in response to cold stress. PLos One. 2013, 8 (3): e58740-10.1371/journal.pone.0058740.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  37. 37.

    Jaillon O, Aury JM, Noel B, Policriti A, Clepet C, Casagrande A, Choisne N, Aubourg S, Vitulo N, Jubin C, Vezzi A, Legeai F, Hugueney P, Dasilva C, Horner D, Mica E, Jublot D, Poulain J, Bruyère C, Billault A, Segurens B, Gouyvenoux M, Ugarte E, Cattonaro F, Anthouard V, Vico V, Fabbro CD, Alaux M, Gaspero GD, Dumas V, et al: The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla. Nature. 2007, 449 (7161): 463-467. 10.1038/nature06148.

    Article  CAS  PubMed  Google Scholar 

  38. 38.

    Wu KL, Guo ZJ, Wang HH, Li J: The WRKY family of transcription factors in rice and Arabidopsis and their origins. DNA Res. 2005, 12 (1): 9-26. 10.1093/dnares/12.1.9.

    Article  CAS  PubMed  Google Scholar 

  39. 39.

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

    Article  CAS  PubMed  Google Scholar 

  40. 40.

    Ling J, Jiang WJ, Zhang Y, Yu HJ, Mao ZC, Gu XF, Huang SW, Xie BY: Genome-wide analysis of WRKY gene family in Cucumis sativus. BMC Genomics. 2011, 12: 471-10.1186/1471-2164-12-471.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  41. 41.

    Li HL, Zhang LB, Guo D, Li CZ, Peng SQ: Identification and expression profiles of the WRKY transcription factor family in Ricinus communis. Gene. 2012, 503 (2): 248-253. 10.1016/j.gene.2012.04.069.

    Article  CAS  PubMed  Google Scholar 

  42. 42.

    Tripathi P, Rabara RC, Langum TJ, Boken AK, Rushton DL, Boomsma DD, Rinerson CI, Rabara J, Reese RN, Chen XF, Rohila JS, Rushton PJ: The WRKY transcription factor family in Brachypodium distachyon. BMC Genomics. 2012, 13: 270-10.1186/1471-2164-13-270.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  43. 43.

    Wei KF, Chen J, Chen YF, Wu LJ, Xie DX: Molecular phylogenetic and expression analysis of the complete WRKY transcription factor family in maize. DNA Res. 2012, 19 (2): 153-164. 10.1093/dnares/dsr048.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  44. 44.

    Lee TH, Tang H, Wang X, Paterson AH: PGDD: a database of gene and genome duplication in plants. Nucleic Acids Res. 2013, 41 (Database issue): D1152-D1158.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  45. 45.

    Kalde M, Barth M, Somssich IE, Lippok B: Members of the Arabidopsis WRKY group III transcription factors are part of different plant defense signaling pathways. Mol Plant Microbe Interact. 2003, 16 (4): 295-305. 10.1094/MPMI.2003.16.4.295.

    Article  CAS  PubMed  Google Scholar 

  46. 46.

    Wang NA, Xiang Y, Fang LC, Wang YJ, Xin HP, Li SH: Patterns of gene duplication and their contribution to expansion of gene families in grapevine. Plant Mol Biol Report. 2013, 31 (4): 852-861. 10.1007/s11105-013-0556-5.

    Article  CAS  Google Scholar 

  47. 47.

    Gregersen PL, Holm PB: Transcriptome analysis of senescence in the flag leaf of wheat (Triticum aestivum L.). Plant Biotechnol J. 2007, 5 (1): 192-206. 10.1111/j.1467-7652.2006.00232.x.

    Article  CAS  PubMed  Google Scholar 

  48. 48.

    Zhang CQ, Xu Y, Lu Y, Yu HX, Gu MH, Liu QQ: The WRKY transcription factor OsWRKY78 regulates stem elongation and seed development in rice. Planta. 2011, 234 (3): 541-554. 10.1007/s00425-011-1423-y.

    Article  CAS  PubMed  Google Scholar 

  49. 49.

    Yu FF, Huaxia YF, Lu WJ, Wu CG, Cao XC, Guo XQ: GhWRKY15, a member of the WRKY transcription factor family identified from cotton (Gossypium hirsutum L.), is involved in disease resistance and plant development. BMC Plant Biol. 2012, 12: 144-10.1186/1471-2229-12-144.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  50. 50.

    Knight H, Knight MR: Abiotic stress signalling pathways: specificity and cross-talk. Trends Plant Sci. 2001, 6 (6): 262-267. 10.1016/S1360-1385(01)01946-X.

    Article  CAS  PubMed  Google Scholar 

  51. 51.

    Park CY, Lee JH, Yoo JH, Moon BC, Choi MS, Kang YH, Lee SM, Kim HS, Kang KY, Chung WS, Lim CO, Cho MJ: WRKY group IId transcription factors interact with calmodulin. Febs Lett. 2005, 579 (6): 1545-1550. 10.1016/j.febslet.2005.01.057.

    Article  CAS  PubMed  Google Scholar 

  52. 52.

    Zhang JZ, Creelman RA, Zhu JK: From laboratory to field. Using information from Arabidopsis to engineer salt, cold, and drought tolerance in crops. Plant Physiol. 2004, 135 (2): 615-621. 10.1104/pp.104.040295.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  53. 53.

    Yan L, Liu ZQ, Xu YH, Lu K, Wang XF, Zhang DP: Auto- and cross-repression of three Arabidopsis WRKY transcription factors WRKY18, WRKY40, and WRKY60 negatively involved in ABA signaling. J Plant Growth Regul. 2013, 32 (2): 399-416. 10.1007/s00344-012-9310-8.

    Article  Google Scholar 

  54. 54.

    Medina J, Catala R, Salinas J: The CBFs: three Arabidopsis transcription factors to cold acclimate. Plant Sci. 2011, 180 (1): 3-11. 10.1016/j.plantsci.2010.06.019.

    Article  CAS  PubMed  Google Scholar 

  55. 55.

    Huang T, Duman JG: Cloning and characterization of a thermal hysteresis (antifreeze) protein with DNA-binding activity from winter bittersweet nightshade, Solanum dulcamara. Plant Mol Biol. 2002, 48 (4): 339-350. 10.1023/A:1014062714786.

    Article  CAS  PubMed  Google Scholar 

  56. 56.

    Mare C, Mazzucotelli E, Crosatti C, Francia E, Stanca AM, Cattivelli L: Hv-WRKY38: a new transcription factor involved in cold- and drought-response in barley. Plant Mol Biol. 2004, 55 (3): 399-416. 10.1007/s11103-004-0906-7.

    Article  CAS  PubMed  Google Scholar 

  57. 57.

    Wang F, Hou XL, Tang J, Wang Z, Wang SM, Jiang FL, Li Y: A novel cold-inducible gene from Pak-choi (Brassica campestris ssp chinensis), BcWRKY46, enhances the cold, salt and dehydration stress tolerance in transgenic tobacco. Mol Biol Rep. 2012, 39 (4): 4553-4564. 10.1007/s11033-011-1245-9.

    Article  CAS  PubMed  Google Scholar 

  58. 58.

    Shekhawat UKS, Ganapathi TR, Srinivas L: Cloning and characterization of a novel stress-responsive WRKY transcription factor gene (MusaWRKY71) from Musa spp. cv. Karibale Monthan (ABB group) using transformed banana cells. Mol Biol Rep. 2011, 38 (6): 4023-4035. 10.1007/s11033-010-0521-4.

    Article  CAS  PubMed  Google Scholar 

  59. 59.

    Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25 (17): 3389-3402. 10.1093/nar/25.17.3389.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  60. 60.

    Arnold K, Bordoli L, Kopp J, Schwede T: The SWISS-MODEL workspace: a web-based environment for protein structure homology modelling. Bioinformatics. 2006, 22 (2): 195-201. 10.1093/bioinformatics/bti770.

    Article  CAS  PubMed  Google Scholar 

  61. 61.

    Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011, 28 (10): 2731-2739. 10.1093/molbev/msr121.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  62. 62.

    Li J, Wang N, Xin H, Li S: Overexpression of VaCBF4, a Transcription Factor from Vitis amurensis, Improves Cold Tolerance Accompanying Increased Resistance to Drought and Salinity in Arabidopsis. Plant Mol Biol Report. 2013, 31 (6): 1518-1528. 10.1007/s11105-013-0627-7.

    Article  CAS  Google Scholar 

  63. 63.

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

  64. 64.

    Dash S, Van Hemert J, Hong L, Wise RP, Dickerson JA: PLEXdb: gene expression resources for plants and plant pathogens. Nucleic Acids Res. 2012, 40 (Database issue): D1194-D1201.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  65. 65.

    Cramer GR, Ergul A, Grimplet J, Tillett RL, Tattersall EAR, Bohlman MC, Vincent D, Sonderegger J, Evans J, Osborne C, Quilici D, Schlauch KA, Schooley DA, Cushman JC: Water and salinity stress in grapevines: early and late changes in transcript and metabolite profiles. Funct Integr Genomics. 2007, 7 (2): 111-134. 10.1007/s10142-006-0039-y.

    Article  CAS  PubMed  Google Scholar 

  66. 66.

    Espinoza C, Vega A, Medina C, Schlauch K, Cramer G, Arce-Johnson P: Gene expression associated with compatible viral diseases in grapevine cultivars. Funct Integr Genomics. 2007, 7 (2): 95-110. 10.1007/s10142-006-0031-6.

    Article  CAS  PubMed  Google Scholar 

  67. 67.

    Fung RWM, Gonzalo M, Fekete C, Kovacs LG, He Y, Marsh E, McIntyre LM, Schachtman DP, Qiu WP: Powdery mildew induces defense-oriented reprogramming of the transcriptome in a susceptible but not in a resistant grapevine. Plant Physiol. 2008, 146 (1): 236-249.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  68. 68.

    Albertazzi G, Milc J, Caffagni A, Francia E, Roncaglia E, Ferrari F, Tagliafico E, Stefani E, Pecchioni N: Gene expression in grapevine cultivars in response to Bois Noir phytoplasma infection. Plant Sci. 2009, 176 (6): 792-804. 10.1016/j.plantsci.2009.03.001.

    Article  CAS  Google Scholar 

  69. 69.

    Koyama K, Sadamatsu K, Goto-Yamamoto N: Abscisic acid stimulated ripening and gene expression in berry skins of the Cabernet Sauvignon grape. Funct Integr Genomics. 2010, 10 (3): 367-381. 10.1007/s10142-009-0145-8.

    Article  CAS  PubMed  Google Scholar 

  70. 70.

    Vega A, Gutierrez RA, Pena-Neira A, Cramer GR, Arce-Johnson P: Compatible GLRaV-3 viral infections affect berry ripening decreasing sugar accumulation and anthocyanin biosynthesis in Vitis vinifera. Plant Mol Biol. 2011, 77 (3): 261-274. 10.1007/s11103-011-9807-8.

    Article  CAS  PubMed  Google Scholar 

  71. 71.

    Ruijter JM, Ramakers C, Hoogaars WM, Karlen Y, Bakker O, van den Hoff MJ, Moorman AF: Amplification efficiency: linking baseline and bias in the analysis of quantitative PCR data. Nucleic Acids Res. 2009, 37 (6): e45-10.1093/nar/gkp045.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  72. 72.

    Hellemans J, Mortier G, De Paepe A, Speleman F, Vandesompele J: qBase relative quantification framework and software for management and automated analysis of real-time quantitative PCR data. Genome Biol. 2007, 8 (2): R19-10.1186/gb-2007-8-2-r19.

    PubMed Central  Article  PubMed  Google Scholar 

Download references


The authors thank Professor Yuepeng Han (Key Laboratory of Plant Germplasm Enhancement and Specialty Agriculture, Wuhan Botanical Garden, The Chinese Academy of Sciences) for critical review of this manuscript. This work was supported by the National Natural Science Foundation of China (NSFC Accession No.: 31130047, 31000902).

Author information



Corresponding authors

Correspondence to Shaohua Li or Haiping Xin.

Additional information

Competing interest

The authors declare that they have no competing interests.

Authors’ contributions

HPX, LNW, SHL and JPL designed and oversaw the research. LNW, LCF, XMS, LYS, ZCL and NW performed the research. LNW and WZ performed bioinformatics analysis, including gene identification and microarray data analysis. LNW, HPX, JPL and SHL wrote the article. All authors read and approved the final manuscript.

Electronic supplementary material


Additional file 1: Figure S1: Chromosomal location of 57 VvWRKYs. VvWRKY03 was located on ‘chromosome 1 random’ and VvWRKY04 was located on ‘chromosome unknown’. Neither was shown here. (JPEG 535 KB)


Additional file 2: Figure S2: The models of conserved amino acid sequences of WRKY domain and zinc-finger structure in four groups. The numbers behind the charts indicated gene numbers in each group. (JPEG 41 KB)


Additional file 3: Table S1: The coefficient of variation of the corresponding treatment means (CV) and probe set IDs of VvWRKYs in 9 experiments. A higher CV means the expression of the probe set (gene) is affected by treatments in an experiment. Five VvWRKYs that didn’t show any changes in any treatments are marked by green color. (XLS 12 KB)


Additional file 4: Table S2: Cold stress-related VvWRKYs obtained in one of three experimental methods. Yellow, red and blue forms represent genes obtained via qRT-PCR, gene-chip data and transcriptome data respectively. Exogenous ABA induced VvWRKYs were shown in green color in form. (DOCX 16 KB)


Additional file 5: Figure S3: Quantitative RT-PCR assays of the expression level of 18 VvWRKYs under cold treatment. The transcription level of these genes didn’t show significant changes during cold treatment in V. vinifera.(JPEG 862 KB)


Additional file 6: Figure S4: Quantitative RT-PCR assays of the expression patterns of 44 VvWRKYs under exogenous ABA treatment. The transcription level of these genes didn’t show significant changes during exogenous ABA treatment in V. vinifera.(JPEG 2 MB)

Additional file 7: Table S3: The primers used for expression pattern analysis for VvWRKYs. (DOCX 17 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is 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 credited. The Creative Commons Public Domain Dedication waiver ( ) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Wang, L., Zhu, W., Fang, L. et al. Genome-wide identification of WRKY family genes and their response to cold stress in Vitis vinifera. BMC Plant Biol 14, 103 (2014).

Download citation


  • WRKY transcription factor family
  • Grapevine
  • Biotic and abiotic stress
  • Cold stress