Chinese wild V. quinquangularis cv. “83-4-96,” “Danfeng-2” and “Shang-24,” V. amurensis cv. “Tonghua-3,” V. davidii cv. “Tangwei,” V. vinifera cv. “Muscat Rose,” and V. labrusca × V. vinifera cv. “Kyoho” genotypes were obtained from the Grape Repository (34°20′N, 108°24′E) of the Northwest A&F University, Yangling, Shaanxi, China. All grapevines were trained to a T-trellis with a planting density of 1.5 m between plants in 2.5 m rows. Berries were judged to be ripe based on data from the previous season's ripening dates and soluble solids levels. Approximately 300 berries total were randomly selected from each cultivar. Approximately 80–100 berries were harvested from three separate locations, respectively, with each group of 80–100 berries considered as one replicate, resulting in three total replicates in each case.

Stilbenes were extracted according to methods described by Sun (Sun et al., 2006) with slight modifications. In brief, 2 g of berry tissue was extracted in 10 ml 80% (v/v) methanol. The samples were then treated with ultrasound for 30 min and kept at 4°C in darkness for 12 h. The extracts were centrifuged at 6000 rpm for 15 min at 4°C, and the pellets re-extracted twice more. The three resulting supernatants were then combined and vacuum-dried using a rotary evaporator, RE-52AA (Shanghai Jinpeng Analytical Instruments Co. Ltd, Shanghai, China). The dried samples were dissolved in 1 mL of pure methanol, filtered through a 0.22 μm filter and stored at −20°C pending further analysis.

Quantitative analysis of stilbenes was performed using a Waters 600E-2487 HPLC system (Waters, USA) equipped with an Agilent ZORBAX SB-C18 column (5 μm, 4.6 × 250 mm). The column temperature was set at 25°C and detection of stilbenes was achieved at 306 nm with a flow rate of 0.8 ml/min. The solvent system consisted of 1 min in isocratic 20% aqueous acetonitrile, 30 min in a linear gradient from 20 to 75% aqueous acetonitrile, 2 min in a linear gradient from 75 to 100% acetonitrile, 3 min in isocratic 100% acetonitrile, 1 min in a linear gradient from 100 to 20% aqueous acetonitrile, and 10 min in isocratic 20% aqueous acetonitrile. Trans-resveratrol was purchased from Sigma-Aldrich while cis-resveratrol was obtained from trans-resveratrol by conversion using UV-C irradiation (Wang et al., 2013b).

Mature grape berries were irradiated with UV-C (254 nm, Spectroline, Model ZQJ-254, output 300 mW/cm²) after harvest for 10 min at a distance of 15 cm. Negative controls consisted of non-irradiated material. Following irradiation, both UV-treated and control plants were transferred to a Robert manual incubator (Model PRX-350D), and incubated in the dark at 25°C with a relative humidity of 80% for 60 h. After each incubation period (0, 4, 12, 24, 36, 48, and 60 h) a subset of berries was collected, frozen in liquid nitrogen and stored at −80°C. Each treatment was performed with at least three independent replicates, and each replicate consisted of randomized plant material.

Two biological replicates of each treatment were used for all RNA-Seq experiments. Total RNA was extracted from whole berries of Chinese wild V. amurensis cv. Tonghua-3 using an SDS/phenol method (Zhang et al., 2003). RNA quality and quantity were monitored on a 1.2% denatured agarose gel and with a NanoDrop 1000 Spectrophotometer (Thermo Scientific, Wilmington, DE, USA), respectively. Treated berries were defined as those that had been irradiated with UV-C and frozen in liquid nitrogen 4 or 24 h post-treatment, while non-irradiated berries at the same time points were used as controls.

Approximately 10 μg cDNA from each of the eight samples (two biological replicates for each treatment) were sequenced in the Genomics Facility at Cornell University's Biotechnology Resource Center (http://www.biotech.cornell.edu/brc/genomics-facility). Strand-specific RNA-Seq libraries were constructed as previously described (Zhong et al., 2011) and sequenced on an Illumina HiSeq 2000 system in the single-end mode. The length of the reads was 100 bp. RNA-Seq reads were first aligned to ribosomal RNA sequences using Bowtie (Langmead et al., 2009) and the aligned reads were removed. The resulting filtered reads were then aligned to the grape genome sequence using TopHat (Trapnell et al., 2009). Following alignment, the count of mapped reads from each sample was derived and normalized to RPKM (reads per kilobase of exon model per million mapped reads). Differentially expressed genes (DEGs) between UV-C treated and control samples at each time point were identified using the DESeq 1.8.3 package (Anders and Huber, 2010) with the raw count data. Raw P-values were adjusted for multiple testing using a false discovery rate (FDR; Benjamini and Hochberg, 1995). Genes with an FDR value of < 0.05 and a fold-changes > 2.0 were considered to be DEGs. Correlations between the 4 h and 24 h time points, along with the untreated control samples, were investigated using cluster analysis with the Cluster 3.0 software.

Multi-array log2 transformation, normalization, and t-tests were performed to identify DEGs between any two selected groups using the geWorkbench genomic analysis software suite (http://www.geworkbench.org) (Floratos et al., 2010). Hierarchical clustering was performed using Pearson's metrics and total linkage algorithms. Gene Ontology and Pathway analysis of co-expressed genes was performed using Plant MetGenMAP (http://bioinfo.bti.cornell.edu/cgi-bin/MetGenMAP/home.cgi) (Joung et al., 2009). To adjust the P-value, a hypergeometric test with the Benjamini and Hochberg false discovery rate (FDR) was performed using the default parameters.

First-strand cDNA was synthesized from 1 mg DNase-treated (TaKaRa Biotechnology) total RNA with the PrimeScript™ II 1st Strand cDNA Synthesis Kit (TaKaRa Biotechnology, Dalian, China). Quantitative real-time RT-PCR was conducted using SYBR green (TaKaRa Biotechnology) on an IQ5 real-time PCR machine (Bio-Rad, CA, USA). Each reaction was carried out in triplicate with a reaction volume of 25 μl. Cycling parameters were 95°C for 30 s, 40 cycles of 95°C for 5 s, and 60 °C for 30 s. For dissociation curve analysis, a program including 95°C for 15 s, followed by a constant increase from 60°C to 95°C, was included after the PCR cycles. The expression of grape ACTIN1 (GenBank Accession number AY680701), amplified with primers F (5′ - GAT TCT GGT GAT GGT GTG AGT - 3′) and R (5′ - GAC AAT TTC CCG TTC AGC AGT - 3′), as well as grape GAPDH (GenBank Accession number CB973647), amplified with primers F (5′ - TTC TCG TTG AGG GCT ATT CCA - 3′) and R (5′ - CCA CAG ACT TCA TCG GTG ACA - 3′), were used as internal controls. Relative expression levels were analyzed using the IQ5 software and the normalized-expression method.

A stilbene synthase (VaSTS) promoter region was isolated by PCR using Chinese wild V. amurensis genomic DNA as a template. Genomic DNA was extracted with phenol/chloroform, precipitated with ethanol and then stored at −40°C for subsequent cloning. Primary PCR with the gene-specific F (5′ - GAA TTC AAT AAA TCT AAT AAA TAT T - 3′) and R (5′- GAC GTC AAT CAG ACT GGT AGA TAC AG - 3′) primers was performed to amplify the randomly selected VaSTS gene (GSVIVT01010565001) promoter, which was then ligated into the cloning vector pMD 19-T (TaKaRa Biotechnology, Dalian, China) and transformed into Escherichia coli strain DH5a (TaKaRa Biotechnology, Dalian, China). A positive candidate clone was sequenced at Invitrogen Biotechnology (Invitrogen Biotechnology, Shanghai, China), and the promoter sequence was analyzed using the PlantCARE (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/) database (Lescot et al., 2002).

Two expression vectors, pC0380GUS and pC35SGUS (Xu et al., 2010), were used for transient expression assays. The 1,570bp upstream region of VaSTS was amplified by PCR as described above, with PstI and EcoRI restriction enzyme sites introduced at the 5′ ends of each primer, respectively. The VaSTS promoter fragment was then double-digested with PstI/EcoRI, and ligated into the corresponding site of the pC0380GUS vector, which had been double-digested with the same restriction enzymes. The fusion construct was verified by sequencing and introduced into Agrobacterium tumefaciens strain GV3101 using the freeze-thaw method.

Agrobacterium-mediated transient expression assays were performed as previously described (Sparkes et al., 2006) utilizing the promoter-GUS fusion constructs. Fully expanded, infiltrated tobacco leaves were collected from intact tobacco (Nicotiana tabacum) plants harboring each construct and subsequently used to evaluate GUS activity via both histochemical and qRT-PCR analyses. Forty-eight hours following A. tumefaciens infiltration, tobacco leaves were treated with UV-C by irradiating (254 nm, Spectroline, Model ZQJ-254, output 300 mW/cm²) for 15 min at a distance of 15 cm, while controls consisted of non-irradiated material. Leaves were sampled at the indicated times for qRT-PCR analysis, or 24 h after A. tumefaciens inoculation for GUS assays.

GUS transcript levels in UV-C treated or untreated tobacco leaves were quantified by two-step qRT-PCR. Total RNA was isolated from tobacco leaf samples using the E.Z.N.A. ® Plant RNA Kit (Omega Bio-tek, USA, R6827-01), and the RNA concentration was adjusted to 500 ng for reverse transcription. The reactions were carried out in triplicate in a final volume of 25 μl using an iCycler iQ5 thermal cycler (Bio-Rad, Hercules, CA, USA), and a SYBR Premix Ex^TM TaqII kit (Takara). Cycling parameters were 95°C for 30 s, followed by 40 cycles of 95°C for 5 s and 62°C for 30 s. The tobacco ubiquitin gene (UBI, GenBank Accession number U66264) was used as the internal control. GUS transcript levels were calculated using the normalized-expression method. The following primers were used for real-time PCR amplifications: 5′ - ATTATGCGGGCAA CGTCTGGTATCAG - 3′ and 5′ - CAT CGGCTTCA AATGGCGTAGC - 3′ for GUS, 5′ - ATGAACGC TGGCGGCATG CTTA - 3′ and 5′ - AGATCTGCATTC CTCCCCTCGCTA - 3′ for UBI.

Histochemical and quantitative GUS assays were carried out according to the procedure published by Jefferson (1987). Histochemical staining of the leaves was performed as previously described (Xu et al., 2010). Quantitative GUS assays, using 4-methyl umbelliferyl glucuronide (MUG; Sigma-Aldrich) as a substrate, were quantified using a Hitachi 850 fluorescence spectrophotometer (Hitachi, Tokyo, Japan), and protein concentrations were measured using the protein–dye binding assay (Bradford, 1976) and a Nicolet Evolution 300 UV-VIS spectrophotometer (Thermo Electron Corp., Madison, WI, USA) with bovine serum albumin as the standard. GUS activity was expressed as nM of 4-methylumbelliferone (4-MU; Sigma-Aldrich) generated per minute per milligram of soluble protein. Data were subjected to one-way ANOVA for comparison of the means, and one-sided paired t-test were performed to assess significant differences between mock and treatment conditions using the SigmaPlot 10.0 statistical software package (Ashburn, VA, USA). Differences at P < 0.05 were considered significant.

Berries from a total of 7 grape accessions, corresponding to five wild species, were analyzed for resveratrol content (Table 1). The highest levels of resveratrol were found in V. amurensis Tonghua-3 (2.19 μg/g FW at maturity); therefore, this genotype was chosen for further study.

Following irradiation treatment of mature berries, trans-resveratrol levels were significantly higher by 12 h post-treatment than non-treated samples and continued to increase steadily until 60 h post-treatment, when the concentration reached 56.76 μg/g FW (Figure 1). Conversely, the content of cis-resveratrol at each time point measured did not change significantly in response to the UV-C treatment (Figure 1).

The transcriptome of mature grape berries was analyzed in response to UV-C treatment 4 and 24 h post-irradiation. At these two time points, trans-resveratrol levels were approximately equal to, and 4-fold greater than, those of the untreated control, respectively. Low quality, adaptor and barcode sequences, as well as possible rRNA reads were removed to obtain clean reads (Table S1). The percentage of reads mapping to the grape genome was approximately 70%, with the exception of repeat 1 of the untreated 4 h sample (4h-CK1) (60.12%) and repeat 2 of the UV-C treated 4 h sample (4h-T2) (50.34%). As confirmation of the robustness of the RNA-Seq dataset, we found that the two biological replicates were in good agreement with respect to gene expression levels, with 0.95 < R² < 0.99 (Table S2).

Correlations between the 4 and 24 h post-treatment time points, along with the untreated control samples, were investigated using cluster analysis (Figure 2; Eisen et al., 1998). As expected, the 4h-CK and 24 h-CK samples (untreated controls) exhibited fewer differences in gene expression between them than did 4h-T and 24h-T samples (UV-C irradiated) (Figure 2). Interestingly, there were even more substantial differences between the two UV treated samples and their untreated controls, respectively, which correlated well with our resveratrol content results.

Gene expression levels were measured as RPKM values and genes exhibiting an FDR < 0.05 and a fold-change in transcript expression levels > 2.0 between UV-C treated and untreated controls were regarded as DEGs. We identified 785 genes with increased transcript abundance and 97 genes with decreased transcript abundance in grape berries 4 h after UV-C treatment compared to untreated controls, with an even larger number of genes (819 with increased transcript abundance and 239 with decreased transcript abundance) 24 h after UV-C treatment. Of the 785 genes with increased transcript abundance at 4 h, 401 exhibited differences in expression at 24 h post-treatment, whereas only 54 down-regulated genes displayed differences in expression between 4 and 24 h post-treatment (Figure S2, Tables S3, S4).

A total of 662 of the 881 DEGs identified at 4 h post-treatment, and 824 of the 1058 DEGs identified at 24 h post-treatment, were assigned to at least one term within the GO categories of “molecular function,” “cellular component,” or “biological process.” These DEGs were further classified into 100 functional subcategories (Figure 3; Table S5).

Of the genes showing increased expression 4 h post-irradiation in the “biological process” category, “response to stress” and “metabolic processes” were the most represented terms, indicating that extensive metabolic activity was taking place in the grape berries following UV-C treatment. Within the “molecular function category,” the “catalytic activity,” “protein binding” and “binding” terms were prominently represented, while the “membrane” and “plasma membrane” terms dominated the “cellular component” category (Figures S3, S4; Tables S6, S7).

We further identified GO terms within the “biological process” category that were over-represented (p < 0.05) both 4 and 24 h after UV-C treatment (Tables 2, 3). These GO terms served as indicators of significantly different biological processes occurring within berries following UV-C treatment, compared to their untreated control counterparts. GO terms such as “metabolic processes,” “cellular processes,” “cellular metabolic processes” and “biosynthetic processes” were enriched in the annotations of both sets of transcripts from the two time points, suggesting that a set of longer term biological activities are required to repair UV-C damage. However, we also noted a few striking differences between enriched GO terms from the two time points. For example, GO terms related to “protein binding,” “nucleotide binding,” “purine nucleotide binding” and “adenyl nucleotide binding” were highly enriched in DEGs 4 h post-treatment, but not at 24 h post-treatment (Tables 2, 3).

To better understand the functional roles of UV-C responsive genes involved in different metabolic pathways, we identified biochemical pathways affected by UV-C irradiation based on our expression profiling analysis. A total of 51 and 76 biochemical pathways were significantly influenced by UV-C treatment (p-value < 0.05) at 4 and 24 h post-treatment, respectively, (Figure 6 and Table S8). These pathways encompassed the biosynthesis or degradation of diverse metabolites including hormones, sugars and polysaccharides, amino acids, fatty acids and lipids, and secondary metabolites.

Secondary metabolic pathways, such as the biosynthesis of flavonoids, phenylpropanoids (Figure 4), chlorogenic acid (Figure S5), free phenylpropanoid acid, acetyl-CoA (from pyruvate), as well as acetoin biosynthesis III and the pathway leading to the biosynthesis of phenylalanine, tyrosine, and tryptophan, all had associated genes whose expression was significantly altered 4 h and 24 h following UV-C treatment. In the case of the stilbenoid and flavonoid biosynthetic pathways, 22 genes encoding stilbene synthases (STS) exhibited 8- to 210-fold increases in their expression levels 4 h after UV-C treatment, and 29- to 64-fold increases in expression 24 h post-irradiation (Figure 4). Several genes encoding enzymes involved in phenylpropanoid biosynthesis were also highly up-regulated by the UV-C treatment (Figure 4; Tables S3, S4). These enzymes included phenylalanine ammonia-lyase (PAL) proteins, which catalyze the first step in the biosynthesis of phenylpropanoids, cinnamate 4-hydroxylase (C4H), and 4-coumarate-CoA ligases (4CL) (Figure 4). Differentially expressed UV-C responsive MYB transcription factors (transcription factors; TFs) were also identified in this study (Figure 5), with a total of 9 MYB TFs identified 4 h and 24 h following UV-C treatment, including 5 with higher and 4 with lower expression.

To further validate our RNA-Seq expression profile data, we performed qRT-PCR assays of 20 selected DEGs (Figure 7 and Figure S6, Table S9). Correlations between RNA-Seq and qRT-PCR results were evaluated by comparing fold-change measurements, whereby scatterplots were generated using the log₂ fold-changes obtained using the RNA-Seq and qRT-PCR methods (Table S10). A close correlation (R² = 0.72) was observed between the log₂-fold changes measured by RNA-Seq and qRT-PCR (Figure S7).

To investigate potential regulatory cis-acting elements within a 1.5 kb region upstream (GenBank Accession number KU376402) of one of VaSTS genes, the sequence was analyzed using the PlantCARE database. This investigation revealed that the promoter contained typical TATA and CAAT boxes, as well as several light-regulated elements, such as a 3-AF1 binding site, G-box, AT1-motif, GAG-motif, GT1-motif, and 4cl-CMA2b element, which are homologous to cis-acting elements previously identified in plants. We also identified other cis-regulatory elements, such as a HSE motif, MBS motif, AAGAA-motif, GCN4_motif, AC-I element, Skn-1 motif and circadian element in the promoter region (Table 4). Most of the predicted motifs are known to be involved in responses to environmental stress, which suggests that VaSTS plays a role in defense. Based on the results of the analysis shown in Table 4, it can also be inferred that VaSTS may be a light-inducible gene.

To test the light-inducible activity of the VaSTS promoter, a 1570 bp fragment upstream from the transcriptional start site was fused to the GUS reporter gene. The construct was introduced into tobacco leaves and these were tested for GUS activity after UV-C irradiation treatment (Figure 8). A CaMV35S::GUS (pC35SGUS) construct was used as a positive control, and a promoter-less construct (pC0380GUS) served as the negative control. The response of the VaSTS promoter was initially analyzed by histochemical GUS staining 24 h after UV-C treatment. As shown in Figure 8A, compared with levels in untreated leaves, GUS activity was significantly greater in UV-C treated leaves. The CaMV35S promoter showed no significant induction after UV-C treatment, while the negative control and wild-type only exhibited low background levels of GUS activity regardless of the treatment (Figure S8).

To further confirm that the VaSTS promoter contributes to the differential transcriptional regulation of expression following UV-C treatment, the promoter-GUS fusion constructs were used in an Agrobacterium-mediated transient tobacco expression system and subsequently analyzed using qRT-PCR. This data showed kinetic and dynamic changes of promoter-driven GUS expression in response to UV-C irradiation at the transcriptional level. Compared with levels in untreated leaves, GUS expression was greatly increased in UV-C treated leaves (Figure 8B). None of the wild-type tobacco leaves treated with UV-C or untreated leaves showed any detectable GUS expression, whereas high levels of GUS expression were found in light-induced and untreated leaves infiltrated with the CaMV35S::GUS construct. In summary, there was a low and comparable level of basal GUS expression in the transformed tobacco leaves at 0 h, which transiently increased in the UV-C-treated leaves at 0–6 h. A peak in GUS mRNA levels in transformed tobacco leaves harboring the VaSTS promoter constructs and treated with UV-C irradiation was observed 6 h after treatment, and thereafter GUS transcript levels decreased to relatively low levels. No significant change in expression was seen in the untreated plants.

Phenolic compounds in grapes can be divided into two groups: flavonoids and non-flavonoids (Teissedre et al., 1996). Resveratrol (trans-3, 5, 4′-trihydroxystilbene) is a non-flavonoid polyphenol with a wide range of beneficial biological properties. To date, resveratrol has been found in at least 72 plant species, distributed among 31 genera and 12 families (Counet et al., 2006), and a number of these are components of the human diet, such as cranberry, peanut, cocoa, and especially grape and its products, such as wine. To examine the biochemical basis of stilbene accumulation in grape berries, we compared their profiles in mature fruit from various accessions following UV-C treatment, with a focus on resveratrol accumulation, and found that V. amurensis cv. Tonghua-3 berries showed the highest resveratrol contents. Furthermore, resveratrol levels in this accession were found to increase between 4 and 60 h post-treatment (Figure 2), which is consistent with previous findings (Douillet-Breuil et al., 1999; Cantos et al., 2001; Liu et al., 2010; Crupi et al., 2013; Wang et al., 2013a).

For each RNA-Seq sample, approximately 60% of the reads could be assigned to grape genes and were used for gene expression profiling (Table S1). There was a high correlation among biological replicates, indicating the robustness of the RNA-Seq dataset (Table S2). Moreover, a very close correlation was observed between relative expression levels based on RNA-Seq and qRT-PCR data (Figure 7; Table S10), suggesting a high level of consistency between the two analytical techniques.

We identified 881 and 1058 genes that exhibited changes in transcript abundance at 4 and 24 h post-UV-C treatment, respectively (Figure S2). Analysis of these DEGs indicated that UV-C irradiation influenced many defense and general stress response pathways, and promoted the expression of genes involved in UV-C protection or the alleviation of UV-C damage; results that are congruent with previous studies of grapevine following UV-B treatment (Pontin et al., 2010).

Increased expression of genes with a putative function in the response to biotic stimuli was a common feature of UV-C-irradiated berries at both the 4 and 24 h post-treatment time points (Tables 2, 3), which is consistent with the fact that UV-C elicited responses are known to share gene activation and signal transduction pathways commonly associated with biotic stress (Yao et al., 2012; Maharaj et al., 2014). Accordingly, several genes putatively involved in pathogen signal transduction and defense responses were also found to have higher expression in treated berries, such as those encoding NBS-LRR type disease resistance proteins (GSVIVT01018889001 and GSVIVT01030098001), an Avr9/Cf-9 rapidly elicited protein (GSVIVT01032581001) and syringolide-induced proteins (GSVIVT01022165001) (Tables S3, S4). These overlapping responses following UV-C exposure, pathogenesis and wounding stimuli may result from common signaling molecules mediating wound/ defense responses, suggesting the use of shared components, and possibly a mechanism of cross-tolerance. Thus, the effects of UV-C light on the induction of biotic stress associated-genes may improve plant resistance to pathogens and pests.

Plant hormones play important roles in diverse growth and developmental processes, as well as in various plant responses to biotic and abiotic stress. In particular, ethylene, auxins, abscisic acid (ABA) and brassinosteroids have been associated with UV-acclimated phenotypes (Jansen, 2002; Ulm et al., 2004; Hectors et al., 2007; Berli et al., 2010), although other hormones, such as cytokinins (CTKs), gibberellins (GAs), jasmonates (JAs) and salicylic acid (SA) may also be involved (Piotrowska and Bajguz, 2011). In the present study, the expression of genes involved in the biosynthesis of JA, SA, and ethylene from methionine, were significantly altered following UV-C irradiation (Figures S9–S11; Table S8).

Genes involved in JA biosynthesis were up-regulated both 4 and 24 h after UV-C treatment (Figure S9), which correlates well with the observation that JA biosynthetic genes are responsive to various types of stress (Doornbos et al., 2011; Hind et al., 2011; Carvalhais et al., 2013; Van der Does et al., 2013). Similarly, we found that the expression of genes related to both SA and ethylene biosynthesis were up-regulated in response to UV-C treatment (Figures S10, S11; Table S8), which may result in an increase in SA and ethylene production. Both ethylene and SA are known to play key roles in biotic stress signaling in response to pathogen infection (Glazebrook, 2005), which again indicates an overlap between the various stress response pathways. Ethylene response factor (ERF) proteins, which belong to a large family of TFs that bind to the promoter of many ethylene inducible defense-related genes, are an integral component of this mechanism (Hao et al., 1998). We observed that the expression of several annotated ERF genes (GSVIVT 01034563001, GSVIVT01004798001, GSVIVT01013913001, GSVIVT01028050001 and GSVIVT01021146001) was up-regulated following UV-C treatment, both at 4 h and 24 h post-irradiation (Figure 5). In addition to its important role in pathogen defense, ethylene also plays a vital role in grape berry development and ripening, including the regulation of genes involved in anthocyanin biosynthesis and accumulation (Chervin et al., 2004). Our results therefore suggest that ethylene may also participate in the modulation of gene expression responses to UV-C radiation (Table S8 and Figure S10).

Many plant TFs are known to be rapidly induced by UV-C and likely play roles in UV-C-induced responses. In this study, we observed that the expression of TFs belonging to the AP2/EREBP family showed the most substantial changes in response to UV-C treatment, followed by MYB, WRKY, and NAC TFs (Figure 5; Tables S3, S4).

MYB genes are unique to eukaryotes (Yang et al., 2003; Jiang et al., 2004) and although their function in animals is restricted to the control of cell division and differentiation (Weston, 1998), they play important roles in plants during various physiological processes (Martin and Paz-Ares, 1997). Their expression is modulated by various hormones (Yanhui et al., 2006) and they have been shown to participate in processes related to epidermal cell fate (Ramsay and Glover, 2005), seed development (Penfield et al., 2001), responses to drought (Abe et al., 1997) and cold (Zhu et al., 2005; Agarwal et al., 2006), pathogen-disease resistance (Lee et al., 2001; Vailleau et al., 2002; Mengiste et al., 2003) and sucrose related responses (Teng et al., 2005), among other functions. Our data indicated that many MYB TFs were transcriptionally induced (between 3- and 30-fold) both 4 and 24 h post-UV-C irradiation. In particular, the gene encoding the MYB protein R2R3 Myb14 (GSVIVT01028328001), which may act as a direct activator of the STS genes, was found to be highly up-regulated following UV-C treatment; a result that is in agreement with previous studies (Höll et al., 2013; Fang et al., 2014; Wang et al., 2015). It should be noted that the expression of MYB14 was not induced by UV irradiation in other studies (Pontin et al., 2010; Carbonell-Bejerano et al., 2014; Suzuki et al., 2015); a difference that may reflect the use of different cultivars, tissues, and experimental conditions, highlighting the complexity of the regulatory networks controlled by TFs.

Both NAC and WRKY TFs also play important roles in plants during development, growth and responses to environmental stress stimuli (Olsen et al., 2005b; Guo et al., 2014; Song and Nan, 2014). NAC domain (Petunia NAM and Arabidopsis ATAF1/2 and CUC2) proteins are plant-specific TFs that are expressed in various developmental stages and tissues, although details of their interactions with DNA and with other proteins are still limited (Olsen et al., 2005a,b), while WRKY proteins are known to contain a DNA binding domain bearing an invariant WRKYGQK sequence (after which the domain was named; Eulgem et al., 2000). As was the case for members of the MYB TF family, genes encoding members of the WRKY and NAC families were found to be transcriptionally induced following UV-C treatment (Figure 5), suggesting a role for these TFs in protection against UV-C radiation.

It is known that the biosynthesis of secondary metabolites, such as flavonoids and other UV-C absorbing phenolic compounds, which accumulate in the vacuoles of epidermal cells, is an important factor in UV-C acclimation. Grapevine plantlets exposed to UV-C showed an up-regulation of genes annotated as being involved in secondary metabolism, especially those associated with the phenylpropanoid pathway (Figure 6; Table S8). Indeed, most PAL genes, as well as genes encoding C4H and 4CL, were up-regulated in response to UV-C treatment (Figure 4).

In addition, the resveratrol producing pathway in grape is known to be dependent upon the function of STS genes. In this study, we found 22 STS genes which were actively responsive to UV-C irradiation in grape, with up-regulation occurring in every case. In grapevine, chalcone synthase (CHS) substrates can also be used by STS enzymes for the production of stilbenes, which are classified as phytoalexins because of their role in plant defense mechanisms against fungal pathogens (Langcake and Pryce, 1977; Jeandet et al., 2002). Interestingly, we found that CHS genes (GSVIVT01000521001 and GSVIVT01032968001), which are also involved in flavonoid biosynthesis, were down-regulated after UV-C irradiation (Figure 4). In agreement with our results for the expression of genes encoding STS and CHS (Figure 4; Tables S3, S4), it has previously been found that CHS genes were down-regulated and genes encoding STS were up-regulated after UV-C treatment (Xi et al., 2014; Suzuki et al., 2015). These results suggest that since STS and CHS share a substrate, there may be a competitive and/or inhibitory relationship between them in response to UV-C exposure, which may in turn play a role in resveratrol accumulation in grape berries.

Previous studies have shown that the promoters of STS genes contain several important elements, including Box-W1 and MBS motifs (Xu et al., 2011). Box-W1 (TTGACC), which is similar to the conserved core sequence of the fungal elicitor-responsive element found in Petroselinum crispum, is specifically recognized by a salicylic acid-induced WRKY DNA-binding protein (Rushton et al., 1996), while the MBS element is a binding consensus site for a MYB homolog (Urao et al., 1993). As such, it is possible that specific WRKY and/or MYB TFs that are up-regulated following UV-C treatment may lead to positive regulation of STS genes in grape. However, another possibility, which is not mutually exclusive, is that UV-C radiation is able to stimulate certain STS promoter elements directly through light-responsive elements (Table 4).

In this study, we utilized a VaSTS promoter-GUS fusion construct and transient expression assays in tobacco to determine whether the UV-responsiveness of this gene was due to elements within its promoter. Quantitative real-time RT-PCR data of GUS expression patterns (Figure 8B) and GUS activity data (Figure 8A) suggested that this promoter exhibited significantly higher levels of activity when leaves were exposed to UV-C irradiation. This is in agreement with previous findings in which a grapevine STS promoter isolated from V. vinifera cv. Optima was shown to be activated in response to biotic (Liswidowati et al., 1991; Bais et al., 2000; Keller et al., 2003) and abiotic stresses, such as wounding (Jeandet et al., 1997). In addition, our in silico analysis indicated that the VaSTS promoter contained several light-responsive elements, including a GAG-motif, which is part of the rbcA conserved DNA module array (rbcA-CMA1) involved in light responsiveness, and an L-box, which is a cis-acting regulatory element associated with the A- and P-box (Logemann et al., 1995). These results, along with our UV-C treatment data, indicate that the expression of the VaSTS gene is regulated via motifs in its promoter that are responsive to UV-C light irradiation.

In conclusion, our transcriptome analysis revealed that UV-C irradiation had regulatory effects on defense and general stress response pathways, as well as on genes involved in UV-C protection or the alleviation of UV-C damage, which likely contributes to the acclimation of plants to UV-C exposure. We have also identified several potential structural and regulatory genes, as well as other factors, which may play a role in the accumulation of resveratrol in grape berries following UV-C exposure. These results provide an important framework for the further elucidation of precise mechanism by which resveratrol biosynthesis is induced by UV-C and other types of stress.

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.