The Grapevine E3 Ubiquitin Ligase VriATL156 Confers Resistance against the Downy Mildew Pathogen Plasmopara viticola

Downy mildew, caused by Plasmopara viticola, is one of the most severe diseases of grapevine (Vitis vinifera L.). Genetic resistance is an effective and sustainable control strategy, but major resistance genes (encoding receptors for specific pathogen effectors) introgressed from wild Vitis species, although effective, may be non-durable because the pathogen can evolve to avoid specific recognition. Previous transcriptomic studies in the resistant species Vitis riparia highlighted the activation of signal transduction components during infection. The transfer of such components to V. vinifera might confer less specific and therefore more durable resistance. Here, we describe the generation of transgenic V. vinifera lines constitutively expressing the V. riparia E3 ubiquitin ligase gene VriATL156. Phenotypic and molecular analysis revealed that the transgenic plants were less susceptible to P. viticola than vector-only controls, confirming the role of this E3 ubiquitin ligase in the innate immune response. Two independent transgenic lines were selected for detailed analysis of the resistance phenotype by RNA-Seq and microscopy, revealing the profound reprogramming of transcription to achieve resistance that operates from the earliest stages of pathogen infection. The introduction of VriATL156 into elite grapevine cultivars could therefore provide an effective and sustainable control measure against downy mildew.


Introduction
Disease resistance in plants is controlled by multiple genetic and epigenetic mechanisms, including general physical barriers and chemical defense compounds as well as the recognition of pathogen-associated molecular patterns (PAMPs) or race-specific pathogen effectors via pattern recognition receptors (PRRs) or receptors encoded by resistance (R) genes, respectively. Signaling via such receptors results in the transcriptional and metabolic reprogramming of plant cells to improve defense and prevent colonization by pathogens [1]. The genes and metabolites involved in resistance have been determined in model plant species, including key signal transduction components and metabolic pathways that produce antimicrobial compounds. This knowledge can now be used to generate disease-resistant varieties of important crops such as grapevine (Vitis vinifera L.), facilitated V. vinifera ATL156 genes was investigated by functional analysis of the promoter regions in silico and in vivo in transgenic A. thaliana plants, and the involvement of ATL156 in resistance to downy mildew was studied by the overexpression of VriATL156 in transgenic V. vinifera plants, which reduced their susceptibility to downy mildew and induced a profound transcriptional reprogramming of defense-related genes.

Bioinformatic Analysis of the V. vinifera and V. riparia ATL156 Genes
VviATL156 (VIT_05s0077g01970) is the V. vinifera ATL gene most closely related to A. thaliana defense-related gene AtATL2 [18,20]. VviATL156 features a 954-bp coding sequence without introns, flanked by 5 and 3 untranslated regions (UTRs) of 98 and 322 bp, respectively. Our previous microarray analysis of grapevine plants infected with P. viticola revealed the weak induction of VviATL156 24 hpi in V. vinifera, whereas the V. riparia ortholog was induced earlier (12 hpi) and more strongly [17].
The differential response of the ATL156 gene in V. vinifera and V. riparia was initially investigated in silico by comparing the coding sequences and promoters (Supplementary Figures S1 and S2). Due to the lack of the V. riparia genome at that time, the VriATL156 gene was amplified from V. riparia cv. Gloire de Montpellier genomic DNA using primers designed for the corresponding V. vinifera gene. The resulting amplicon was 960 bp in length, indicating that VriATL156 is intronless like its ortholog in V. vinifera and most other ATL genes [20,23]. Comparison of the VviATL156 and the VriATL156 coding sequences revealed very high similarity, with only five single nucleotide polymorphisms (SNPs) and one additional tandem duplication of a 6-bp sequence that is actually a tandem repeat in V. vinifera and a triplet repeat in V. riparia (Supplementary Figure S1a). At the protein level, three of the five SNPs were silent and the other two resulted in changes from Asp (D) to Glu (E) at position 246 and from Thr (T) to Pro (P) at position 309 in the V. riparia polypeptide. The 6-bp tandem duplication caused the insertion of Pro (P) and Gln (Q) at positions 313 and 314, respectively, in the V. riparia polypeptide (Supplementary Figure S1b).
The 1174-bp VviATL156 promoter sequence was cloned for re-sequencing, allowing us to design primers to clone a corresponding 1186-bp VriATL156 promoter fragment. Sequence analysis showed that the two promoter regions shared 96% identity (Supplementary Figure S2). The prediction of promoter motifs in silico using PlantCARE [24] also revealed a near identical distribution and frequency of cis-regulatory elements between the two promoters, including the presence of two TATA boxes (Supplementary Figure S3a). Neither of the promoters featured pathogenesis-related cis-regulatory elements, but both contained hormone-response motifs, particularly for salicylic acid [25] and jasmonic acid [26]. Many of the predicted motifs have been associated with transcription factors involved in stress responses. The most abundant motifs in both promoters were those recognized by basic leucine zipper (bZIP) transcription factors, which are known to regulate stress and defense responses and well as light-dependent signaling. Binding sites for WRKY and MYB transcription factors were also abundant in both promoters [27]. Other motifs were identified as binding sites for basic helix loop helix (bHLH), DNA-binding with one finger (Dof) and AT-Hook family transcription factors, which are not specific for defense or stress signaling but regulate developmental processes (Supplementary Figure S3b). The high level of similarity between the promoters did not reveal any clear distinctions that would explain the differential regulation of the V. vinifera and V. riparia ATL156 genes, so we decided to test their functions directly in A. thaliana.

The V. riparia ATL156 Promoter is More Responsive to Pathogens and Hormone Treatment
The V. vinifera and V. riparia ATL156 promoters were fused to the gusA reporter gene, allowing their activity to be compared in transgenic A. thaliana plants by analyzing βglucuronidase (GUS) activity. Two independent lines for each promoter were examined under normal growth conditions throughout the developmental cycle from seedling emergence to flowering. GUS activity was weaker in lines carrying the V. riparia promoter than in lines carrying the V. vinifera promoter, as shown by the intensity of GUS staining (Supplementary Figure S4a). This was in agreement with the higher basal levels of ATL156 mRNA detected in V. vinifera compared to V. riparia (Supplementary Figure S4b) and confirmed the suitability of the heterologous A. thaliana system for the analysis of these promoters.
To assess the response of the ATL156 promoters to pathogens, transgenic A. thaliana lines expressing each promoter-reporter system were infected with avirulent or virulent strains of Pseudomonas syringae pv. tomato DC3000 (Pst), the former carrying the avirulence gene AvrB. The analysis of GUS activity revealed that both promoters were induced by the pathogen, regardless of the virulence of the strain ( Figure 1). However, there was a significant difference in the induction ratio, with the V. vinifera promoter triggering a 1.6-2-fold increase in GUS levels and the V. riparia promoter triggering a 4-5-fold increase, in both cases relative to a mock infection. Thus, the V. riparia promoter clearly showed a stronger response to the pathogen.
The V. vinifera and V. riparia ATL156 promoters were fused to the gusA reporter gene, allowing their activity to be compared in transgenic A. thaliana plants by analyzing β-glucuronidase (GUS) activity. Two independent lines for each promoter were examined under normal growth conditions throughout the developmental cycle from seedling emergence to flowering. GUS activity was weaker in lines carrying the V. riparia promoter than in lines carrying the V. vinifera promoter, as shown by the intensity of GUS staining (Supplementary Figure S4a). This was in agreement with the higher basal levels of ATL156 mRNA detected in V. vinifera compared to V. riparia (Supplementary Figure S4b) and confirmed the suitability of the heterologous A. thaliana system for the analysis of these promoters.
To assess the response of the ATL156 promoters to pathogens, transgenic A. thaliana lines expressing each promoter-reporter system were infected with avirulent or virulent strains of Pseudomonas syringae pv. tomato DC3000 (Pst), the former carrying the avirulence gene AvrB. The analysis of GUS activity revealed that both promoters were induced by the pathogen, regardless of the virulence of the strain (Figure 1). However, there was a significant difference in the induction ratio, with the V. vinifera promoter triggering a 1.6-2-fold increase in GUS levels and the V. riparia promoter triggering a 4-5-fold increase, in both cases relative to a mock infection. Thus, the V. riparia promoter clearly showed a stronger response to the pathogen. To gain more insight into the molecular basis of this phenomenon, we treated the transgenic plants with various hormones that are known to influence the regulation of ATLs [22]. Both promoters were similarly induced by two different concentrations of jasmonic acid and by 1 mM abscisic acid (ABA), but only the V. riparia promoter was induced by 1 mM salicylic acid or 100 µM ethylene ( Figure 2). Hormonal regulation may therefore facilitate the stronger response of the V. riparia promoter to pathogens. To gain more insight into the molecular basis of this phenomenon, we treated the transgenic plants with various hormones that are known to influence the regulation of ATLs [22]. Both promoters were similarly induced by two different concentrations of jasmonic acid and by 1 mM abscisic acid (ABA), but only the V. riparia promoter was induced by 1 mM salicylic acid or 100 µM ethylene ( Figure 2). Hormonal regulation may therefore facilitate the stronger response of the V. riparia promoter to pathogens.

Transgenic Grapevine Plants Can Express a Functional VriATL156 Gene
We expressed the VriATL156 coding sequence in V. vinifera under the control of the constitutive cauliflower mosaic virus 35S (CaMV 35S) promoter to determine its ability to confer resistance against P. viticola. We transformed V. vinifera cv. Shiraz callus [28] by co-cultivation with Agrobacterium tumefaciens EHA105 carrying the pART27-VriATL156 binary vector or the empty pART27 vector as a negative control. This vector carries a kanamycin resistance gene and an expression cassette for green fluorescent protein (GFP), allowing dual selection by growth on kanamycin-supplemented medium and visual selection for fluorescence (Supplementary Figure S5a). Southern blot analysis revealed a total of 14 plantlets representing six transformation events, which differ in their number of insertions (Supplementary Figure S5b). All plantlets were regenerated and propagated in vitro (Supplementary Figure S5c) before acclimatization in the greenhouse for another three seasons, except for L7 and L23, which did not survive.

Transgenic Grapevine Plants Can Express a Functional VriATL156 Gene
We expressed the VriATL156 coding sequence in V. vinifera under the control of the constitutive cauliflower mosaic virus 35S (CaMV 35S) promoter to determine its ability to confer resistance against P. viticola. We transformed V. vinifera cv. Shiraz callus [28] by cocultivation with Agrobacterium tumefaciens EHA105 carrying the pART27-VriATL156 binary vector or the empty pART27 vector as a negative control. This vector carries a kanamycin resistance gene and an expression cassette for green fluorescent protein (GFP), allowing dual selection by growth on kanamycin-supplemented medium and visual selection for fluorescence (Supplementary Figure S5a). Southern blot analysis revealed a total of 14 plantlets representing six transformation events, which differ in their number of insertions (Supplementary Figure S5b). All plantlets were regenerated and propagated in vitro (Supplementary Figure S5c) before acclimatization in the greenhouse for another three seasons, except for L7 and L23, which did not survive.
VriATL156 expression was evaluated by qPCR in transgenic lines L1, L14, L18 and L19, representing each of the remaining independent transformation events. The mean normalized expression (MNE) ranged from 1 to 6 depending on the transformation event, compared with a value of 0.13 for the endogenous VviATL156 gene in control plants. Line L1 showed the highest expression level of VriATL156, corresponding to a ~45.5-fold increase in expression compared to VviATL156 in the controls (Figure 3a). VriATL156 expression was evaluated by qPCR in transgenic lines L1, L14, L18 and L19, representing each of the remaining independent transformation events. The mean normalized expression (MNE) ranged from 1 to 6 depending on the transformation event, compared with a value of 0.13 for the endogenous VviATL156 gene in control plants. Line L1 showed the highest expression level of VriATL156, corresponding to a~45.5-fold increase in expression compared to VviATL156 in the controls (Figure 3a).
Functional analysis in four lines, one for each available transformation event, was carried out by infecting leaf disks from the third, fourth and fifth leaves of these lines with P. viticola, observed at 6 days post-infection (dpi) for sporulation [29]. The transgenic plants showed enhanced P. viticola resistance compared to vector-only controls, correlating with the VriATL156 expression (Figure 3a-c). Moreover, resistance level correlated with the abundance of total ubiquitinated proteins in transgenic lines L1 and L14, indicating the expression of a functional enzyme (Figure 3d). Infection experiments were carried out at least six times across two consecutive years on at least seven disks per line, with similar sporulation scores. Functional analysis in four lines, one for each available transformation event, was carried out by infecting leaf disks from the third, fourth and fifth leaves of these lines with P. viticola, observed at 6 days post-infection (dpi) for sporulation [29]. The transgenic plants showed enhanced P. viticola resistance compared to vector-only controls, correlating with the VriATL156 expression (Figure 3a-c). Moreover, resistance level correlated with the abundance of total ubiquitinated proteins in transgenic lines L1 and L14, indicating the expression of a functional enzyme ( Figure 3d). Infection experiments were carried out at least six times across two consecutive years on at least seven disks per line, with similar sporulation scores.

Resistance to P. viticola in Lines Expressing VriATL156 is Associated with Pathogen Growth Arrest and Enhanced Hormone Production
Tissues infected with P. viticola were analyzed by microscopy 24, 48 and 96 hpi to visualize the extent of mycelium colonization in the leaves of the two best performing lines (L1 resistant and L14 partially resistant) as well as regenerated, empty vector-transformed plants as controls. The first infection symptoms appeared on control leaves 72 hpi, and sporulation typically occurred 6 dpi. Zoospores were able to germinate at a similar frequency on all plants including controls, but subsequent mycelial growth was severely impaired in line L1, where sporulation was seldom observed (Figure 4), comparable to the outcome of infection in V. riparia [17]. We observed no evidence of necrosis or a hypersensitive response related to the infection.

Resistance to P. viticola in Lines Expressing VriATL156 is Associated with Pathogen Growth Arrest and Enhanced Hormone Production
Tissues infected with P. viticola were analyzed by microscopy 24, 48 and 96 hpi to visualize the extent of mycelium colonization in the leaves of the two best performing lines (L1 resistant and L14 partially resistant) as well as regenerated, empty vector-transformed plants as controls. The first infection symptoms appeared on control leaves 72 hpi, and sporulation typically occurred 6 dpi. Zoospores were able to germinate at a similar frequency on all plants including controls, but subsequent mycelial growth was severely impaired in line L1, where sporulation was seldom observed (Figure 4), comparable to the outcome of infection in V. riparia [17]. We observed no evidence of necrosis or a hypersensitive response related to the infection.
The measurement of hormone levels revealed that line L1 produced higher basal levels of salicylic acid and jasmonic acid, but the levels of ABA and cis-(+)-12-oxophytodienoic acid (OPDA), a precursor of jasmonic acid [30], were unaffected in the absence of infection ( Figure 5).

Transcriptional Profiling of Transgenic Grapevine Plants Expressing VriATL156 Indicates the Strong Activation of Defense-Related Signaling before Infection
To characterize the molecular phenotype of the transgenic plants expressing Vri-ATL156, we compared the most resistant line L1 to vector-only control plants by RNA-Seq analysis before infection with P. viticola and 24 hpi. We chose the 24 hpi time point because resistance is fully established 48 hpi, and transcriptional reprogramming is therefore likely to begin at an earlier time (Figure 4). We identified 2360 transcripts that showed significant differential expression (log2FC > |1|, adj. The measurement of hormone levels revealed that line L1 produced higher basal levels of salicylic acid and jasmonic acid, but the levels of ABA and cis-(+)-12-oxophytodienoic acid (OPDA), a precursor of jasmonic acid [30], were unaffected in the absence of infection ( Figure 5).  The measurement of hormone levels revealed that line L1 produced higher basal levels of salicylic acid and jasmonic acid, but the levels of ABA and cis-(+)-12-oxophytodienoic acid (OPDA), a precursor of jasmonic acid [30], were unaffected in the absence of infection ( Figure 5). Seq analysis before infection with P. viticola and 24 hpi. We chose the 24 hpi time point because resistance is fully established 48 hpi, and transcriptional reprogramming is therefore likely to begin at an earlier time (Figure 4). We identified 2360 transcripts that showed significant differential expression (log2FC > |1|, adj. p-value < 0.05) in L1 vs. control plants in the absence of infection, 1436 of which were upregulated and 924 downregulated in L1 plants, showing that the constitutive expression of VriATL156 induces profound transcriptional reprogramming (Figure 6a; Supplementary Dataset S1). Enrichment analysis revealed the prevalence of the "protein modification" (265 genes) and "reproduction" (55 genes) functional categories among the genes upregulated in line L1, with both categories mainly comprising receptor kinase-encoding genes (Figure 7a; Supplementary Dataset S6).
In the "protein modification" category, genes encoding mitogen-activated protein kinases (such as MPK3 homologous) and protein phosphatases are also present, as well as numerous proteins involved in protein degradation (such as U-box domain containing proteins). As anticipated, this group included VriATL156 itself, which was strongly expressed in the L1 plants relative to the endogenous expression of VviATL156 in the controls (log 2 FC = 4.2). Another key functional category enriched in the L1 plants was "response to biotic stimulus", including genes encoding pathogenesis-related (PR) proteins, phenylalanine ammonia lyases, stilbene synthases and trans-resveratrol di-O-methyltransferases, which are known to be associated with biotic stress responses in grapevine (Figure 7a Enrichment analysis revealed the prevalence of the "protein modification" (265 genes) and "reproduction" (55 genes) functional categories among the genes upregulated in line L1, with both categories mainly comprising receptor kinase-encoding genes ( Figure  7a; Supplementary Dataset S6).  Figure 7. Gene Ontology (GO) functional category enrichment analysis among differentially expressed genes in mock-treated transgenic line L1 (L1m) compared to mock-treated empty-vector control plants (EVm). GO categories of (a) upregulated genes and (b) downregulated genes significantly enriched were analyzed with BINGO. The size of the nodes is proportional to the number of genes in the test set which are annotated to that node while the color of the node represents the corrected p-value. The node color scale ranges from white (not significantly over-represented), yellow (p-value < 0.05) to dark orange (p-value < 10 −5 × 0.05).
Most of these genes were commonly found in the "secondary metabolic process" category, together with genes encoding laccases, endochitinases, glutathion-S-transferases, peroxidases and gdsl esterase lipases, all families well known for their role in plant defense mechanisms (Figure 7a; Supplementary Dataset S6). Only a few functional categories were enriched among the 924 downregulated genes in L1 plants, and the statistical significance of enrichment was less striking than the upregulated genes (Figure 7b; Supplementary Dataset S7). The most notable categories were "lipid metabolic process", including 15 genes encoding GDSL motif-containing lipase/esterase proteins, "secondary metabolic process" (mainly laccase genes), and "response to endogenous stimulus", one third of which were auxinrelated genes, supporting the proposed role of auxin in the susceptibility of plants to biotrophic pathogens [32], as well as several ethylene-responsive transcription factors. The few genes modulated in L1 plants following infection (L1i:L1m) were not enriched for any particular relevant functional category (Supplementary Figure S6). VriATL156 expression therefore triggers the constitutive upregulation of genes primarily involved cell signaling (e.g., receptor-like kinases) and defense responses (pathogenesis-related proteins, enzymes involved in secondary metabolism), indicating the establishment of a pre-infection resistance status in L1 plants correlating with the observed resistance phenotype.
The strong pre-infection expression level of many defense-related genes in L1 plants is reminiscent of the previously reported transcriptional comparison of V. vinifera and V. riparia, where all transcripts modulated in both species following infection with P. viticola were consistently more strongly regulated in V. riparia [16]. Accordingly, we compared the differentially expressed genes in non-infected L1 plants (this report) and in V. riparia plants infected with P. viticola  [20,22]. Moreover, in line with the role of calcium in mediating defense signal transduction and the particular role of calcium-binding proteins in response to pathogens in grapevine [33], several grapevine calmodulin-like proteins (CMLs) are commonly or specifically regulated in V. riparia and L1 plants, thus correlating with the increased plant resistance of both genotypes compared to wild-type V. vinifera. These data suggest that VriATL156 expression controlled by the constitutive CaMV 35S promoter mimics the effect of an infection response in the downy mildew-resistant grapevine species V. riparia, providing strong evidence that ATL156 is a key component in the defense response typical of resistant grapevine genotypes. and L1 plants, thus correlating with the increased plant resistance of both genotypes compared to wild-type V. vinifera. These data suggest that VriATL156 expression controlled by the constitutive CaMV 35S promoter mimics the effect of an infection response in the downy mildew-resistant grapevine species V. riparia, providing strong evidence that ATL156 is a key component in the defense response typical of resistant grapevine genotypes.

Discussion
Genetic resistance to disease is a highly valued grapevine trait because it avoids the need for chemical treatments, particularly the multiple chemical sprays required to control P. viticola. Genetic engineering allows the rapid transfer of such resistance traits to elite varieties without affecting their complex genetic background. Candidate resistance genes must be evaluated carefully, given the long and challenging process of grapevine transformation and regeneration. Several categories of genes can be considered, representing different layers of resistance: genes encoding receptors or signaling factors that regulate downstream defense responses, and genes encoding the components of those defense mechanisms [34]. We targeted the signaling layer, bypassing the specific pathogen recognition step in an attempt to engineer more durable resistance by preventing the adaptation of pathogens to individual receptors [12,35].
The candidate gene VriATL156 from V. riparia is also present in V. vinifera, with minor sequence differences, i.e., the change from polar (Thr) to non-polar (Pro), predominantly affecting the protein C-terminus, which controls substrate recognition by E3 ubiquitin ligases and may therefore specify downstream ubiquitinated targets [36]. Following the recent release of the V. riparia genome [37], the cloned sequences were also compared with

Discussion
Genetic resistance to disease is a highly valued grapevine trait because it avoids the need for chemical treatments, particularly the multiple chemical sprays required to control P. viticola. Genetic engineering allows the rapid transfer of such resistance traits to elite varieties without affecting their complex genetic background. Candidate resistance genes must be evaluated carefully, given the long and challenging process of grapevine transformation and regeneration. Several categories of genes can be considered, representing different layers of resistance: genes encoding receptors or signaling factors that regulate downstream defense responses, and genes encoding the components of those defense mechanisms [34]. We targeted the signaling layer, bypassing the specific pathogen recognition step in an attempt to engineer more durable resistance by preventing the adaptation of pathogens to individual receptors [12,35].
The candidate gene VriATL156 from V. riparia is also present in V. vinifera, with minor sequence differences, i.e., the change from polar (Thr) to non-polar (Pro), predominantly affecting the protein C-terminus, which controls substrate recognition by E3 ubiquitin ligases and may therefore specify downstream ubiquitinated targets [36]. Following the recent release of the V. riparia genome [37], the cloned sequences were also compared with those predicted from the genome assembly, retrieving an additional SNP at the 3 terminus, leading to an Arg to Gln amino acid substitution at position 316 of the corresponding proteins (Supplementary Figure S7), which would further support the observed differences between ATL156 from V. vinifera and V. riparia in the C terminal domain involved in substrate recognition. The overexpression of VriATL156 in V. vinifera driven by the CaMV 35S promoter increased resistance to P. viticola, most effectively in line L1 and to a lesser extent in line L14, correlating with the transgene expression levels. This indicates a key role for this E3 ligase in disease resistance, in line with previous work on other ATL genes and co-expression network analysis [20,22,23]. RNA-Seq analysis showed that line L1 underwent broad transcriptional reprogramming relative to the control line, particularly involving genes in functional categories related to signal transduction and defense. These results indicate that defense mechanisms are already active at a basal level in plants overexpressing VriATL156 and are not induced further by infection. Indeed, based on the functional analysis of the VriATL156 and VviATL156 promoters in the same A. thaliana genetic background, it is tempting to speculate that the higher basal activity of the V. vinifera promoter may not be sufficient to provide an effective response in the absence of strong post-infection induction, as observed in V. riparia. In contrast, the very high basal expression of VriATL156 driven by the CaMV 35S promoter in line L1 may be sufficient to reach the threshold required to achieve resistance in V. riparia following infection. All genes induced in common between uninfected L1 plants and infected controls were expressed at a much higher level in the L1 plants. An interesting parallel can be drawn with the previous report comparing gene expression in V. vinifera and V. riparia, showing that all genes commonly modulated in these species in response to infection were induced more strongly in V. riparia [17]. This suggests that infection triggers a powerful amplification of the response in resistant species.
Line L1 also differed from the control in terms of hormonal signaling, with higher levels of salicylic and jasmonic acids in the uninfected L1 plants, in parallel with the downregulation of a number of auxin-related genes. The suppression of auxin-related genes is evidence that the hormone balance tips in favor of resistance in response to the transcriptional reprogramming, based on previous reports that auxins antagonize salicylic acid responses to favor growth and promote developmental processes, associated with susceptibility to biotrophic pathogens [33].
The molecular basis of the improved resistance in the transgenic lines expressing VriATL156 is not clear because the upstream regulators and downstream targets have not been identified, as for other ATL E3 ubiquitin ligases including AtATL2 [38,39]. We previously reported that VviATL156 is co-regulated with several WRKY transcription factors across a number of plant-pathogen interactions, suggesting that one or more of these WRKY proteins may be regulators or targets of ATL156 [22]. Likewise, the E3 ubiquitin ligase EIRP1 targets the negative regulator of resistance VpWRKY11 in Vitis pseudoreticulata [40], and the E3 ubiquitin ligase ATL31 in A. thaliana is transcriptionally regulated by WRKY33 [41]. Additional evidence of ATL protein interactors includes EDR1 kinase, which negatively regulates AtATL1 to control programmed cell death and disease resistance [42], and the 14-3-3 protein targeted for degradation by ATL31 [43]. ATL31 also interacts with the SNARE protein SYP121 (involved in papilla formation), conferring resistance against powdery mildew in A. thaliana [43]. However, ATL31 does not seem to target SYP121 for degradation, probably acting as a regulator via other mechanisms [44].
The ability of ATL proteins to regulate their targets via mechanisms other than protein degradation reflects the well-established diverse effects of different modes of ubiquitination (mono-ubiquitination vs. poly-ubiquitination), different types of ubiquitin chain linkages, and different ubiquitin chain lengths. These assorted ubiquitin chains perform a range of functions, including the regulation of iron uptake, the repair of double-stranded breaks in DNA, and the regulation of auxin responses [45]. However, these diverse roles are poorly understood in plants [46,47]. Because the transgenic plants expressing VriATL156 accumulate higher levels of ubiquitinated proteins, despite the absence of treatment with a proteasome inhibitor, the proteins ubiquitinated by VriATL156 may not be targeted for degradation. Resistance may not reflect the ubiquitin-proteasome system (UPS)-mediated degradation of particular negative defense regulators, but an unknown signaling role mediated by ubiquitin attachment, generating more canonical post-translational modifications [45].
To gain more insight into the function of VriATL156, we compared the transcriptomes of line L1 and control plants. The overexpression of VriATL156 affected thousands of transcripts, resulting in profound transcriptional reprogramming that redirected the transcriptome to activate defense-related signaling. Functional category enrichment revealed that these plants appear to be primed for defense, especially at the level of "protein modification", with the upregulation of many genes encoding receptor kinases and putative E3 ubiquitin ligases, including ATLs. Some of these effects may reflect the constitutive expression of VriATL156 driven by the CaMV 35S promoter. However, the genes encoding defense effectors such as PR proteins and stilbene synthases, which are typically associated with disease resistance in grapevine, though enriched were less represented among the modulated transcripts compared with signaling components. This correlates with the lack of adverse phenotypic effects often observed when plants overexpress constitutive defenses, reflecting the damage incurred by plant tissues, for example, due to prolonged exposure to ROS [48].
In conclusion, our functional analysis has shown that VriATL156 can regulate the resistance response in grapevine, and that the overexpression of this gene using a constitutive promoter leads to the strong pre-infection modulation of multiple components of the resistance signaling cascade. Given the vast transcriptomic reprogramming of resistance signaling components in the transgenic plants, this kind of resistance may be more durable than that achieved by expressing classical R genes because pathogens do not have the opportunity to evade specific receptors. The comparison of ATL156 native promoters in two grapevine species also provides insight into the contribution of VriALT156 in the resistant species V. riparia, which could facilitate the identification of pathogen-responsive motifs suitable for cisgenic strategies or genome editing approaches.

Cloning the VriATL156 Coding Sequence
Genomic DNA extracted from the leaves of V. riparia Gloire de Montpellier plants [51] was amplified with High Fidelity DNA Polymerase (KAPA Biosystems, Roche, Basel, Switzerland) and the following primers, designed using Primer3 [52] and prepared by Eurofins Genomics (Milan, Italy): ATL2_int_For (5 -GAT CGA CTT CCG CTT TGT GCT CC-3 ) and ATL2_int_Rev (5 -GCA CAC GGC ACA TTC AG AC-3 ). The VriATL156 gene was verified by sequencing (BMR Genomics, Padua, Italy) and transferred to a pGEM vector under the control of the CaMV 35S promoter for cloning in Escherichia coli. The cassette was then transferred to pART7 [28] and finally isolated as NotI fragment and inserted at the NotI site of the modified binary vector pART27-4a (also containing an eGFP5 expression cassette driven by the A. thaliana ubiquitin 10 promoter) for grapevine transformation. The vector (10 ng) was introduced into A. tumefaciens strain EHA105 by electroporation (25 µF, 2.5 KVolt, 2000 Ohm). Transformed cells were rapidly collected with 1 mL of Lysogeny broth (LB) medium and incubated for 1 h at 28 • and 200 rpm. Then, 75-150-225 µL of cell suspension was inoculated into a solid selective medium and incubated at 28 • C for 48 h.

Cloning the VviATL156 and VriATL156 Promoter Regions
A putative promoter region extending~1.2 kb upstream of the transcriptional start site was isolated from both genes using primers ATL156for (5 -CAC CGA CCT TTT ACT ATA CAT AGG TAG-3 ) and ATL156rev (5 -TTT CTC CGA TCT GAA AGA GAT AAA G-3 ), carrying the CACC and the TT clamps (bold italics) for directional cloning in pENTR™/D-TOPO ® (Invitrogen) and transformation into E. coli Top10 competent cells.

Analysis of Promoter Responses in Transgenic A. thaliana
The cloned promoter regions of VviATL156 and VriATL156 were transferred to vector pKGWFS7 and introduced into A. tumefaciens strain GV3101 for the transformation of A. thaliana Col-0 plants by floral dipping [53]. T3 transgenic lines carrying the gusA gene under the control of the VriATL156 promoter (Vri lines 4 and 8) or the VviATL156 promoter (Vvi lines 3 and 10) were grown in Murashige and Skoog (MS) medium from 3 days to 4 weeks with a 16-h photoperiod at 24 • C. The same lines were also cultivated in soil pots for 4 weeks for infection challenge experiments with Pseudomonas syringae pv. tomato (PstDC3000). Basal promoter activity was measured after 24, 48 and 96 h, and after 1, 2, 3 and 4 weeks. Three plants from each line were sampled at each time point for GUS staining and the measurement of GUS activity.
Chemical stimuli were applied to 5-day-old transgenic plants from the four lines described above as well as vector-only controls. The chemicals were applied to plants in 96well plates containing half-strength liquid MS medium supplemented with methylsalicylic acid, methyljasmonate, abscisic acid or ethephon (ethylene donor), at indicated concentrations, and promoter activity was tested after 24 h. For the challenge studies, transgenic plants were infected with Pseudomonas syringae pv. tomato strains AvrB (avirulent on Col-0) or DC3000 (virulent) at 2 × 10 7 CFU/mL, by the infiltration of three biological replicates using a 1 mL needless syringe. Mock inoculations were carried out with 10 mM MgCl 2 .

GUS Staining and GUS Fluorimetric Assay
Whole plants or detached leaves (A. thaliana lines Vri 4.17, Vri 8.4, Vvi 3.14 and Vvi 10.7) were submerged in GUS staining solution (2.5 mM Na 3 PO 4 , pH 7.0, 0.1% Triton x-100, 100 µM X-Gluc) in a vacuum chamber for 5 min and then incubated at 37 • C for 5 h before clearing in 70% ethanol and visualization with a LEICA MZ16F stereomicroscope. The GUS fluorimetric assay was carried out as previously described (DyNA Quant Application Note 3) using a pool of three leaves from each biological replicate. Fluorescence was measured using a 3300 Nanodrop fluorimeter (Thermo Fisher Scientific, Waltham, MA, USA). Proteins were quantified in 1 mL of extract using Bradford reagent (Bio-Rad Laboratories, Hercules, CA, USA).

Transformation of V. vinifera
Pre-embryogenic callus of V. vinifera cv Shiraz was generated from anthers collected at the pre-flowering stage and initiated on callus initiation medium (PIV) medium before transfer to charcoal containing medium (GS1CA) C1 medium [54]. Transformation was carried out as previously described [55,56]. Callus was selected on C1 medium supplemented with 150 µg/mL kanamycin, and stably transformed embryos were screened for the presence of GFP, transferred to shooting medium and then deflasked as previously described [57].
Transformation was confirmed by Southern blot analysis using a probe designed to specifically recognize a sequence encompassing the final portion of the ubiquitin terminator (GFP expression cassette) and a portion of the T-DNA region next to the left border. The probe was labeled with alkaline phosphatase using the AlkPhos Direct Kit (GE Healthcare, Amersham, UK).
Transgenic plants plus controls (regenerated Shiraz vines transformed with the empty vector pART27) were regularly micropropagated in vitro and cultivated in rooting medium in a growth chamber at 27 • C with a 16-h photoperiod. Rooted plantlets were transferred to soil and cultivated in the greenhouse under shade for~3 weeks before returning to normal conditions. Plants were tested for resistance about 3 months after transfer to the soil.

Infection with P. viticola
An isolate of P. viticola harvested from a local field in the Verona province (Italy) in 2013 was propagated axenically as previously described [17]. For experimental infections, leaf disks were cut from the sterilized 3rd, 4th and 5th leaves from the apex of plant shoots and inoculated with 50 µL drops (30,000 sporangia/mL). Inoculated disks were kept in the dark overnight and inoculum droplets were removed after 24 h. Distilled water drops were used as control and collected at the same time points. Disks were incubated in a growth chamber with a 16-h photoperiod at 21-24 • C and were monitored for 6 days, until sporulation was observed on the controls. For microscopy, leaf disks were collected 24, 48 and 96 hpi, stained with 0.05% (w/v) aniline blue in 0.1% (w/v) Na 2 CO 3 (pH 10) and observed under a Leica DM/RB epifluorescence microscope (Leica Microsystems, Wetzlar, Germany) with excitation band pass filter (BP) 340-380 nm, dichroic mirror 400 nm and suppression log pass filter (LP) > 430 nm. For RNA-Seq analysis, leaf disks were collected 24 hpi, immediately frozen in liquid nitrogen and stored at -80 • C. Three independent biological replicates were artificially infected for analysis.

Analysis of Endogenous Hormone Levels
We extracted 200 mg of powdered leaf tissue in 2 mL methanol-water (7:3 v/v) containing 1% acetic acid by sonication for 15 min at 4 • C. Extractions were performed from 3 independent samples. The samples were centrifuged at 16,000× g for 5 min at 4 • C and the pellet was re-extracted with the same solvent. The supernatants were pooled and evaporated in a speed-vac to a final volume of 100 µl before the analysis of 2 µL samples by HPLC-MS/MS on an Agilent 1100 HPLC system (Agilent Technologies, Santa Clara, CA, USA) connected to an LTQ Orbitrap mass spectrometer (Thermo Fisher Scientific) equipped with a Luna-HILIC column (100 mm × 2.0 mm, 5 µm column; Phenomenex, Torrance, CA, USA). The solvent gradient was 50% A (acetonitrile containing 1% formic acid) and 50% B (water containing 1% formic acid) to 100% B over 20 min at a flow rate of 1.5 mL/min. The injection volume was 2 µL.

Quantitative Real-Time PCR (qPCR)
Grapevine leaf disks (in biological triplicate) were ground in liquid nitrogen and total RNA was isolated from 200 mg of powder per sample, using the Spectrum Plant Total RNA Kit (MilliporeSigma). The quantity, integrity and purity of the RNA were determined using a Nanodrop 2000 (Thermo Fisher Scientific) and a Bioanalyzer Chip RNA 7500 series II (Agilent Technologies). The RNA was treated with DNase from the TURBO DNA-free kit (Thermo Fisher Scientific) and 1 µg was reverse transcribed using the SuperScript III kit (Invitrogen, Thermo Fisher Scientific). qPCR was carried out using the SYBR Green PCR Master Mix kit (Applied Biosystems, Waltham, MA, USA). The ubiquitin1-encoding gene (VIT_16s0098g01190) was used as the reference for expression normalization. The list of primers is provided in Supplementary Table S1. The "Comparative Quantitation" protocol was applied using the Mx3000P Real-time qPCR system (Stratagene, Agilent Technologies). The real amplification efficiency was calculated using LinRegPCR [58]. Mean normalized expression (MNE) values were calculated according to [59]. Standard errors (SE) were calculated from three biological replicates.

Western Blot
Leaves from transgenic lines L1 and L14 and control (regenerated empty vectortransformed) V. vinifera cv. Shiraz plants were ground in liquid nitrogen and protein was extracted from 100 mg of powder as previously described [60]. Samples (equivalent to 35 µg total protein) were fractionated by SDS-PAGE and transferred to a nitrocellulose membrane (GE Healthcare) at 100 V for 1 h. Membranes were incubated with a primary anti-ubiquitin antibody 1:2000 (Enzo Life Sciences, Inc.) and a secondary anti-mouse-HRP IgG antibody 1:50,000 in Tris buffer saline medium supplemented with 0.1% Tween 20 (TBS/T) with milk 5%, for 2 h and 1.5 h, respectively, at room temperature. Finally, bands were revealed by exposing membranes to chemiluminescent substrate (Amersham™ ECL Western Blotting Detection Kit, GE Healthcare). Images were acquired with a ChemiDoc imaging system (Bio-Rad).

RNA-Seq Analysis
For each sample, an unstranded library was prepared from 2.5 µg of total RNA using the Illumina TruSeq Library Prep Kit v2 (Illumina, San Diego, CA, USA) and was sequenced on an Illumina HiSeq500 (101-bp, paired-end reads) at the Functional Genomics Center (Department of Biotechnology, University of Verona, Italy). The data were preprocessed and analyzed using high performance computing resources made available by CINECA (Class C ISCRA project IsC33_NobleRot). Raw reads were pre-processed using trimmomatic [61] to remove low-quality bases and adapter sequences. Filtered reads were mapped against the reference grapevine genome [2] using TopHat v2.0.11 [62]. Mapped reads were summarized at the gene level into a count-matrix using the htseq-count tool from the HTSeq library [63]. Gene models provided to the htseq-count tool to score reads mapping unambiguously to a single gene were obtained from the grape genome database release V1 (http://genomes.cribi.unipd.it/DATA/) [64]. Differentially expressed genes were identified and analyzed using EdgeR software [65]. The raw read counts were normalized to the gene GC content using the EDASeq package [66] and sequencing depth among samples (TMM trimmed mean of M values normalization in EdgeR). We considered as differentially expressed all the genes with a p value < 0.05 after false discovery rate correction and a Log2 fold change >|1|.