Next Article in Journal
Production of Terpene Trilactones from Cell and Organ Cultures of Ginkgo biloba
Next Article in Special Issue
The Variability of Berry Parameters Could Be an Indicator of the Potential Quality of the Vineyard
Previous Article in Journal
Introduced Herbivores Threaten the Conservation Genetics of Two Critically Endangered Single-Island Endemics, Crambe sventenii and Pleudia herbanica
Previous Article in Special Issue
Rootstock Effects on Fruit Yield and Quality of ‘BRS Tainá’ Seedless Table Grape in Semi-Arid Tropical Conditions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Transcriptional Landscape of Berry Skin in Red and White PIWI (“Pilzwiderstandsfähig”) Grapevines Possessing QTLs for Partial Resistance to Downy and Powdery Mildews

by
Francesco Scariolo
1,†,
Giovanni Gabelli
1,†,
Gabriele Magon
1,2,
Fabio Palumbo
1,2,
Carlotta Pirrello
3,
Silvia Farinati
1,2,
Andrea Curioni
1,2,
Aurélien Devillars
1,2,
Margherita Lucchin
1,2,
Gianni Barcaccia
1,2 and
Alessandro Vannozzi
1,2,*
1
Department of Agronomy, Food, Natural Resources, Animals and Environment (DAFNAE), University of Padova, Agripolis, 35020 Legnaro, Italy
2
Interdepartmental Centre for Research in Viticulture and Enology, University of Padua, Via XXVIII Aprile, 31015 Conegliano, Italy
3
Research and Innovation Centre, Fondazione Edmund Mach, Via E. Mach 1, 38098 San Michele all’Adige, Italy
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Plants 2024, 13(18), 2574; https://doi.org/10.3390/plants13182574
Submission received: 8 August 2024 / Revised: 6 September 2024 / Accepted: 8 September 2024 / Published: 13 September 2024

Abstract

:
PIWI, from the German word Pilzwiderstandsfähig, meaning “fungus-resistant”, refers to grapevine cultivars bred for resistance to fungal pathogens such as Erysiphe necator (the causal agent of powdery mildew) and Plasmopara viticola (the causal agent of downy mildew), two major diseases in viticulture. These varieties are typically developed through traditional breeding, often crossbreeding European Vitis vinifera with American or Asian species that carry natural disease resistance. This study investigates the transcriptional profiles of exocarp tissues in mature berries from four PIWI grapevine varieties compared to their elite parental counterparts using RNA-seq analysis. We performed RNA-seq on four PIWI varieties (two red and two white) and their noble parents to identify differential gene expression patterns. Comprehensive analyses, including Differential Gene Expression (DEGs), Gene Set Enrichment Analysis (GSEA), Weighted Gene Co-expression Network Analysis (WGCNA), and tau analysis, revealed distinct gene clusters and individual genes characterizing the transcriptional landscape of PIWI varieties. Differentially expressed genes indicated significant changes in pathways related to organic acid metabolism and membrane transport, potentially contributing to enhanced resilience. WGCNA and k-means clustering highlighted co-expression modules linked to PIWI genotypes and their unique tolerance profiles. Tau analysis identified genes uniquely expressed in specific genotypes, with several already known for their defense roles. These findings offer insights into the molecular mechanisms underlying grapevine resistance and suggest promising avenues for breeding strategies to enhance disease resistance and overall grape quality in viticulture.

1. Introduction

Vitis vinifera L. is considered one of the ancient crops in human history, with evidence of its consumption tracing back to the era when humans were still collectors [1]. This species belongs to the Vitis L. genus and includes the wild subspecies V. vinifera subsp. sylvestris and the cultivated subspecies V. vinifera subsp. sativa or vinifera. The belief that cultivated subspecies derive from the domestication of wild individuals is based on the genetic similarity between sylvestris and vinifera, as well as the presence of intermediate genotypes [2].
Until the first half of the 19th century, European viticulture relied on the cultivation of own-rooted grapevines. This changed dramatically when European vineyards encountered a crisis brought on by the introduction of native American pathogens and pests. European grapevines, which had never faced these threats, were particularly vulnerable. The root systems of European grapevines were highly susceptible to phylloxera infestation, a pest to which wild American grapevines were resistant [3]. In addition to pests such as phylloxera, other alien pathogens arrived from the New World, such as powdery mildew (PM) and downy mildew (DM), to impact viticulture. These challenges significantly affected grapevine production, necessitating new strategies and innovations in vineyard management for maintaining yield and health. The causal agents of DM and PM, respectively, the oomycete Plasmopara viticola and ascomycete Erysiphe necator, are obligate, biotrophic, and polycyclic pathogens that infect all green parts of their host plants [4,5]. P. viticola thrives in conditions of high relative humidity and mild temperatures [6,7], while E. necator optimally grows at 85% relative humidity and 26 °C [8]. Under favorable conditions, these pathogens undergo multiple cycles of clonal reproduction, causing severe damage to fruit quality and yield. To combat the impact of DM and PM, fungicides are frequently applied in vineyards. However, the use of these agrochemicals poses potential risks to human health and the environment [9]. Although V. vinifera is susceptible to these pathogens, American and Asian Vitis species show complete or partial resistance, impacting various stages of the pathogen’s life cycle. This involves reducing infection frequency, slowing tissue colonization, extending the latent period, and decreasing spore production, all without completely halting the infection [10]. The genetic basis of resistance traits in grapevine is typically identified using biparental mapping, where genetic maps are created by crossing two heterozygous parents and analyzing their segregating markers in the progeny (pseudo-testcross mapping strategy). An alternative method involves selfing populations. For quantitative trait loci (QTL) analyses, progenies are grouped based on inherited genotypes, and their phenotypes are compared to identify significant associations between traits and allelic variants [10].
To date, phenotyping and genotyping approaches to identify QTLs associated with resistance to pathogens such as P. viticola and E. necator have been applied to over 50 populations, with the development of several mapping populations produced by crossings between V. vinifera and different wild Vitis spp. [11,12,13,14], backcross individuals [15,16,17,18], and hybrid accessions [19,20,21,22,23]. Only one genome-wide association study (GWAS) incorporating both cross-generated individuals and natural grape cultivars has been conducted [24], and only two studies have utilized pedigree information to characterize resistance QTLs [25,26].
In general, progress in identifying resistance (R) loci has been steadily advancing. To date, 33 R loci against P. viticola (Rpv) and 15 against E. necator (Ren or Run) have been identified [10]. Of these, only a small number have been genomically characterized, including Run1, Rdv1, and Rpv1 [23,24,27,28].
Since the advent of first-generation hybrids, extensive efforts have been dedicated to integrating resistance traits from American or Asian wild species into “elite” genetic material. This has led to hybrids combining disease resistance from wild grapevine species with desirable traits of traditional European V. vinifera varieties. Known as PIWI varieties (from the German “pilzwiderstandsfähig”, meaning fungus resistant), their number continues to grow through breeding programs. Currently, over 100 PIWIs are recognized, encompassing both established and newly developed cultivars. Ongoing research and breeding endeavors aim to enhance disease resistance, improve grape quality, and adapt to diverse climates.
Resistant varieties are now widely cultivated globally. Among the top resistant wine grape varieties, Concord is the most prevalent, covering 0.24% of the global vineyard area, primarily in Brazil, the USA, and Japan. Bianca follows with 0.22%, mainly in Russia, Hungary, and Moldova. Other notable varieties include Seyval Blanc, Regent, Chambourcin, Villard Noir, and Baco Noir, each accounting for 0.06–0.02% of the global vineyard area and cultivated across all continents [29,30].
Molecular biology and genomics have assisted in selecting resistant varieties by identifying R loci and establishing genotyping mapping populations. Moreover, the advent of omics technologies has further aided in understanding the molecular mechanisms behind resistance responses from a molecular perspective. On the other hand, high-throughput transcriptomic technologies have allowed a better and more comprehensive understanding of the transcriptional changes occurring during the grapevine response to DM or PM, thus pinpointing key defense regulators to be used as candidate markers of grapevine resistance. Resources aimed at this purpose could include genome-wide transcriptomic studies investigating the response of susceptible or resistant genotypes [31,32,33,34], the comparison between resistant and susceptible varieties [35,36,37,38,39,40,41], or the response of susceptible genotypes treated or not with resistance inducers [42,43].
Most of these studies were performed on leaf tissues pre- and post-pathogen inoculation. This is since, for example, concerning PM, the period of fruit susceptibility is in general small while green tissues like leaves are much more exposed to pathogen attack. After the onset of véraison, grape berries develop ontogenetic defense mechanisms that make them generally more resistant to PM and DM compared to green tissues. This process is part of the plant’s natural evolution to protect mature fruits, which are crucial for reproduction and seed dispersal.
Like other fleshy fruits, grapevine berries are complex organs composed of various tissues that develop according to both common and unique patterns. Grapes are divided into three distinct parts: exocarp (skin), mesocarp (flesh), and seeds [44]. Many of the relevant traits for table or wine grape quality seem to be specifically related to only one of the pericarp tissues. The grapevine skin consists of an outer epidermis (a single cell layer) and an inner hypodermis (a layer of one to seventeen cells). It is the specific site where terpenoid biosynthesis and carotenoid catabolism occur, and where cells lose their chlorophyll and transform into modified plastids during the véraison stage [45]. Additionally, it is the site of anthocyanin pigment synthesis during ripening, primarily responsible for berry and wine color [46,47]. Furthermore, the skin serves as an interface between the internal and external environments, shielding the latter from pathogen invasion and abiotic stresses. This requires a specialized composition of cell walls and layers of protective waxes [48,49]. In the literature, comparative transcriptome studies have been conducted on grape skin. These studies aim to better understand the genetic response related to specific traits influenced by both skin and flesh properties. For example, crispness and firmness are linked to mesocarp cell turgor pressure and skin resistance during various stages of development [50]. However, knowledge about the transcriptional landscape of the skin in resistant PIWI grape varieties is still limited compared to their elite parents.
The aim of this study was to investigate the transcriptional profile of the exocarp in mature berries from four PIWI varieties (two red and two white), and their noble parental varieties, to uncover molecular mechanisms that may contribute to the greater basal resistance observed in interspecific varieties compared to elite varieties. By employing tools such as WGCNA (Weighted Gene Co-expression Network Analysis) and tau analysis, we sought to identify genotype-specific gene clusters. Our findings highlight an increase in vesicular activity in the resistant varieties, which may play a role in their enhanced resistance.

2. Results and Discussion

2.1. Global Transcriptome Analysis of PIWI and Elite Varieties and Identification of Differentially Expressed Genes in Tolerant vs. Susceptible Red and White Varieties

RNA-seq on the NovaSeq 6000 platform produced a total of 760 million 2 × 150 bp reads, with an average of approximately 42 million reads per sample (Table S1). After filtering, 750 million reads were retained for further analysis and mapped to the 12X.v2 reference genome assembly using the VCost.v3 gene annotations [51]. On average, 88% of the reads mapped uniquely to the reference genome, equating to 36.5 million reads per sample. Summarized read count data are presented in Table S2. After removing annotations with fewer than 10 counts per sample, a principal component analysis (PCA) and a heatmap was generated to evaluate the concordance of each replicate within its group and to explore the similarities and differences between varieties.
Replicates were organized in the heatmap by variety and clustered by color (Figure 1A), indicating greater similarity among the white varieties compared to the red ones. Similarly, in the PCA plot (Figure 1B), the first principal component (PC1) accounts for 43.0% of the variance, effectively separating the red and white varieties. The second principal component (PC2), explaining 21.0% of the variance, distinguished between the elite and PIWI varieties.
Using the DESeq2 dds function, differentially expressed genes (DEGs) were identified based on an adjusted p-value (padj) < 0.05 and a log2 fold change (log2FC) > |2|. DEGs from the white and red subsets were analyzed separately, as well as single comparisons between one PIWI and its elite parent (Supplementary Table S3). Analyses revealed distinct gene expression patterns between groups. The number of up- and downregulated DEGs was plotted for each dataset (Figure 1C). Notably, the highest number of DEGs (1430) was observed in the DOWN dataset SN vs. SB suggesting a significant downregulation of genes in the resistant variety compared to its parent. Conversely, the lowest number of DEGs (804) was found in the comparison of both red PIWI varieties (CC and CV) to the noble parental CS, indicating fewer changes in gene expression between these red varieties and their parent. In general, considering the single comparisons between the PIWI and relative elite parent, the white varieties exhibited a greater number of DEGs than the red varieties both in terms of up- and downregulation. This suggests that the white PIWI varieties might undergo more extensive transcriptional reprogramming compared to their red counterparts, possibly contributing to their differential basal resistance profiles. This aligns with previous findings that white PIWIs exhibit better adaptability to variable and stressful conditions, as well as earlier maturation [52].
Further investigation of the intersections between DEGs identified in single comparisons was performed using an upset plot (Figure 1D) showing both the shared and non-shared DEGs. In particular, the number of non-shared DEGs in the red varieties is very low compared to the numbers observed in the white varieties. This observation gains further significance when considering that on a genome-wide scale, the white varieties are more transcriptionally correlated with each other, as indicated in Figure 1A. Shifting to the analysis of DEGs shared among multiple samples, the most represented groups are those shared among the comparisons involving PIWI varieties that have the same noble parent. Specifically, these are the DEGs that are up- and downregulated in the comparisons between CV, CC, and CS, as well as those in the comparisons between SN, SR, and SB. This finding indicates that PIWI varieties derived from the same noble parent share similar transcriptional responses, reflected in the up- and downregulated DEGs. The shared DEGs among CV, CC, and CS suggest a conserved set of gene expression changes inherited from CS. Similarly, the DEGs shared among SN, SR, and SB highlight a common transcriptional profile influenced by Sauvignon blanc. This pattern suggests that the noble parent plays a significant role in shaping the transcriptional landscape of PIWI offspring, which could be leveraged in breeding programs to predict and enhance desirable traits [27,53].
Gene Set Enrichment Analyses of DEGs identified in pairwise comparisons indicate several ontological categories which appear to be enriched in the skin of resistant varieties in comparison to their noble parent (Table S4). Considering categories enriched with a significance level of p > 0.01, a notable observation is that, despite the higher number of DEGs, white PIWIs exhibited fewer enriched categories (58 in SN and 54 in SR) compared to the red ones (176 in CC and 177 in CV). Moreover, these varieties shared very few ontological categories, mainly related to “organic metabolic process” (GO:0005737) and “oxoacid metabolic process” (GO:0044281). In contrast, varieties with the same noble parent showed a much higher number of shared categories, especially the red berry varieties CC and CV, which shared almost all enriched categories (28.3%, 175 out of 178). The white berry varieties showed a lower number of shared categories (31.4%, 27 categories out of 86).
Among the significant categories, we focused on those with the highest fold change in enrichment to avoid categories that, although highly significant (low p-value), were too general to provide specific biological insights. Using this approach, we observed that both the CC and CV genotypes exhibit increased activity related to membrane and vesicle transport. This includes the “late endosome to vacuole transport via multivesicular body sorting pathway” (GO:0032511, FC~3), “endosome transport via multivesicular body sorting pathway” (GO:0032509, FC~3), “multivesicular body sorting pathway” (GO:0071985, FC~3), “AP-type membrane coat adaptor complex” (GO:0030119, FC~3), and “protein localization to chloroplast and vacuole” (GO:0072598 and GO:0072665). Conversely, the white varieties appeared to be enriched in categories related to the biosynthesis and metabolism of organic acids’ metabolism including the “organic hydroxy compound metabolic process” (GO:1901615), “monocarboxylic acid metabolic process” (GO:0032787), “oxoacid metabolic process organic acid metabolic process (GO:0043436), carboxylic acid biosynthetic process” (GO:0046394), and “carboxylic acid metabolic process” (GO:0019752). Also interesting is the enrichment of the “tripeptide transport” category. Looking at the enriched “unique” ontological categories in the DEGs of SR or SN, none stood out for their uniqueness, as many were associated with general processes related to primary or secondary metabolism (Table S4).

2.2. K-Means Associated Weighted Gene Co-Expression Network Analysis Indicates Gene Clusters Associated with PIWI Genotypes and Tolerance

Within plant biology, a plethora of seminal studies have harnessed WGCNA as a valuable tool, evident in diverse plant species including pineapple [54], strawberry [55], pear [56], and grapevine [57,58,59,60]. These studies highlighted the efficacy of WGCNA in unraveling the intricate regulatory networks and uncovering functionally related gene modules within plant systems. In our study, we sought to apply an enhanced version of WGCNA by integrating k-means clustering as an additional step. This powered approach, developed by Botía et al. [61], aimed to improve the precision of gene clustering within the co-expression network analysis compared to traditional WGCNA. Thus, our dataset comprising 23,847 genes was examined using an adapted version of the WGCNA R-package to detect gene co-expression modules. A soft-thresholding power β of 30 was applied to the matrix to establish a scale-free network. Modules were defined as groups of genes tightly connected to each other, exhibiting high correlation within the same module. Initially, a minimum module size of 30 was set, and modules with strongly correlated eigengenes (threshold: 0.25) were merged. A more stringent examination was conducted with a stricter cutoff (0.10). However, given the extensive number of modules obtained (103), these data are provided as supplementary data unless explicitly cited in the main text (Table S5). The eigengene, serving as the principal component, encapsulates the expression profiles of genes within a specific module. Twenty-six unique modules were delineated (threshold: 0.25), each denoted by a specific color and visualized in a hierarchical clustering dendrogram, illustrating their interrelations via individual module eigengenes (Figure 2A). The gene content of these modules varied from a minimum of 255 genes for the light green module to a maximum of 1355 for the blue one (Figure 2B). Subsequently, a correlation analysis was conducted between these 26 modules, the 6 varieties under investigation (CS, CC, CV, SB, SN, and SR), and 2 additional traits that distinguish the genotypes under study: the tolerance/susceptibility to pathogens (T/S) and the grape berry color (GC) (Figure 2C). The objective was to pinpoint gene modules distinctly associated with the diverse genotypes considered, to the tolerance/susceptibility, and to the color. For each variety/trait, at least one highly specific module was identified (correlation p-value < 0.01; Figure 2C, Table S6). However, in certain instances, multiple modules exhibited significant correlations, both positive and negative, with the same variety/trait or multiple varieties/traits displayed significant correlations with the same module. For example, CC was associated with the brown, dark green, light yellow, and royal blue modules, whereas the yellow module was simultaneously associated with SB, CC, and T/S. The most notable correlations between the module eigengene (ME) and variety were observed between the brown module and CC (r = 1, p = 6 × 10−19), the dark red module and CV (r = 0.99, p = 3 × 10−14), and finally, the black module and SR (r = 0.98, p = 6 × 10−12). The most notable correlations between T/S and GC traits were with the yellow (r = −0.99, p = 1 × 10−15) and tan (r = 0.95, p = 2 × 10−9) modules.
To delve deeper into the genetic composition of modules exhibiting notable correlations with the studied varieties/traits, two distinct network metrics, the gene significance (GS) and the module membership (MM), were employed. Module membership (MM) quantifies the degree of correlation between a gene’s expression profile and the respective module eigengene. Conversely, gene significance (GS) serves as an additional network parameter, often represented as the negative logarithm of a p-value, offering an estimation of a gene’s biological relevance. The greater the absolute value of GSi, the greater the biological significance of the i-th gene. In essence, genes demonstrating higher GS and MM values hold increased significance with respect to the phenotypical trait [58]. Consequently, a specific module exhibiting significantly interconnected MM or GS values associated with the tolerance (susceptibility) trait may play a pivotal biological role in resistance mechanisms. Among the considered modules, the brown (CC; cor = 0.99, p < 1 × 10−200), purple (SN; cor = 0.89, p < 1 × 10−200), black (SR; cor = 0.97, p < 1 × 10−200), dark red (CV; cor = 0.97, p < 1 × 10−200), blue (SB; cor = 0.97, p < 1 × 10−200), orange (CS; cor = 0.96, p < 1 × 10−200), yellow (tolerance; cor = −0.98, p < 1 × 10−200), and cyan (GC; cor = 0.86, p < 1 × 10−200) modules displayed the strongest correlations between MM and GS (Figure S1). A Gene Set Enrichment Analysis (GSEA) was performed to ascertain the occurrence of enriched categories related to biological process (BP), molecular function (MF), or cell component (CC) ontologies in all modules obtained by the analysis (Tables S7 and S8). One particularly intriguing aspect pertains to those modules displaying significant correlations with both the T/S trait and one or more examined genotypes, totaling 5 under a 0.25 cutoff and increasing to 18 using a 0.10 cutoff (Table S6). These modules consistently demonstrate a distinct behavior: when associated to the T/S and a susceptible genotype (CS or SB), the correlations exhibit inverse tendencies, whereas their associations with the tolerance trait and a resistant variety consistently display aligned correlations (Figure 3A).
No significant ontologies indicating enrichment were identified within the royal blue module, which exhibited associations with both CC and T/S. Conversely, when examining the tan module, positively correlated with both CV and T/S, a substantial number of enriched categories were evident (Figure 3B). As a general observation, most enriched genes within the top 10 enriched categories ordered by fold change (FC) collectively represent the machinery involved in various stages of vesicle trafficking, including cargo selection, vesicle formation, and membrane fusion. This observation is consistent with the enrichment analyses performed on the DEGs in the pairwise correlations. The fact that two independent analyses highlighted the same category of genes supports the notion that vesicle transport might be a process that specifically distinguishes the red PIWI varieties, particularly CC, and can be associated, at least based on the WGCNA, with tolerance.
Of particular interest was the enrichment observed in a gene class associated with specific adaptor complexes’ proteins (APs) and with the formation of clathrin-coated vesicles. Adaptor protein (AP) complexes are evolutionarily conserved vesicle transport regulators that recruit coat proteins, membrane cargoes and coated vesicle accessory proteins and finely tune their sorting within the cell [62]. In mammalian cells, the AP-1, AP-2, and AP-3 complexes operate in concert with the scaffolding molecule clathrin to generate clathrin-coated vesicles (CCVs), which are involved in various cellular functions such as receptor-mediated endocytosis, protein sorting, and the transport of molecules within the cell. Clathrin-mediated endocytosis (CME) is the primary mechanism by which eukaryotic cells internalize extracellular or membrane-bound cargoes and play crucial roles in plant–microbe interactions [63,64,65]. CME is essential for immune responses mediated by pattern recognition receptor (PRR) kinases, facilitating the internalization and degradation of activated extracellular PRRs in the vacuole [66]. Examples include PRRs such as PEP RECEPTOR 1/2 (PEPR1/2) recognizing endogenous plant peptides, EF-TU receptor (EFR) recognizing bacterial Elongation Factor TU, FLAGELLIN-SENSING 2 (FLS2) recognizing flagellin (flg-22) as a pathogen-associated molecular pattern (PAMP), and the Cf-4 receptor-like protein recognizing the Cladosporium fulvum Avr4 avirulence effector [66,67,68]. Another group of highly represented genes including Vitvi07g00530, Vitvi05g00632, Vitvi19g01853, Vitvi19g00421, Vitvi13g01627, and Vitvi15g00311 encode for subunits of the cytosolic proteins Sec23, Sec24, and Sec31, which are constituents of the COPII complex involved in intracellular protein transport from the endoplasmic reticulum (ER) to the Golgi apparatus [69]. Additionally, Vitvi07g00179, a RANGAP1 ortholog, regulates RAN GTPase activity, which indirectly influences vesicle trafficking (Figure 4).
Although a significant correlation was only found between the tan module and the CV genotype, resistant varieties show a higher expression of vesicle cargo and clathrin-related genes compared to susceptible ones in all the PIWIs analyzed, except for SN (Figure 3C). The resistant varieties exhibit an enrichment of genes involved in vesicular trafficking activity not only between the plasma membrane and internal organelles such as the nucleus, Golgi, and endoplasmic reticulum but also among these organelles. This heightened translocation activity may be attributable to various biological processes already associated with resistance, such as (i) the transport of defense compounds including defense-related molecules such as phytoalexins, antimicrobial peptides, and secondary metabolites; (ii) the regulation of antioxidant systems, important for detoxifying reactive oxygen species (ROS) and protecting cells from oxidative damage; (iii) the maintenance of a higher level of antioxidant capacity in the berry skin, which can contribute to its resilience against biotic and abiotic stresses; (iv) the modulation and strengthening of cell wall integrity; (v) the regulation of hormone signaling networks in plants, including those mediated by salicylic acid (SA), jasmonic acid (JA), and ethylene (ET), which play key roles in plant defense responses. The resistant variety may exhibit constitutive expression of vesicle cargo and clathrin-related genes to modulate hormone signaling pathways in the berry skin, priming it for rapid and effective defense activation upon pathogen attack [70,71,72,73,74,75,76].
Another interesting observation, which could seem counterintuitive at first glance, is the enrichment of defense-related genes in the blue module, which is associated to the susceptible variety SB and is anticorrelated with the T/S trait (Figure 3B,C). Over the 78 defense-related genes in the blue module are RPS genes, R protein PRF disease resistance genes, pathogen-related proteins, MLA10, etc. There could be several reasons why we detected an enrichment of defense genes in the susceptible grapevine variety compared to the resistant one: (i) the susceptible variety might have evolved to maintain a higher baseline expression of defense genes as a proactive strategy against potential threats. This could be due to historical interactions with pathogens or environmental stresses in its natural habitat; (ii) the susceptible variety might have deficiencies in other aspects of its defense system, leading to a compensatory upregulation of certain defense genes to mitigate its vulnerability. In contrast, the resistant variety may have more robust defense mechanisms overall, requiring less reliance on constitutive expression of defense genes; (iii) the constitutive expression of defense genes could also be influenced by environmental factors such as soil composition, climate, or microbial communities present in the vineyard. The susceptible variety might perceive these factors as potential threats and maintain higher levels of defense gene expression as a preemptive response; (iv) there could be trade-offs between growth/development and defense in the susceptible variety, where resources are allocated to defense at the expense of other physiological processes. This could result in the constitutive expression of defense genes even in the absence of infection.

2.3. PIWI Varieties Show Absolutely and Highly Specific Genes Which Are Involved in Disease Resistance

In addition to the WGCNA analysis, we employed an algorithm commonly utilized in transcriptomic investigations involving animals or humans. However, it had previously been applied in the analysis of tissue-specific genes within the P. noir flower, as described by [58]. This algorithm, denoted as the tau (τ) algorithm, was employed to assess the tissue-specificity level of each predicted gene within a given genome [77]. Following the quantile normalization of 23,847 genes, which were selected based on their expression levels, and the subsequent creation of BIN profiles, the τ algorithm was utilized to assign a value ranging from zero (indicating constitutive expression across all or most tissues) to one (indicating absolute specificity for a particular tissue) to each gene. The distribution of τ values across the entire gene set is depicted in Figure 5A.
In summary, 906 genes exhibited a high degree of specificity (referred to as highly specific genes, HSGs, with τ > 0.85), and among them, 570 were identified as absolutely specific genes (absolutely specific genes, ASGs, with τ = 1). It is important to note that the τ value only characterizes the “specificity” of a gene. To ascertain the specific tissue to which a gene is exclusive, we computed the τ expression fractions (τef). Among the various grapevine varieties studied, SR displayed the highest count of HSG (246) and ASG (157). Conversely, SB exhibited the lowest number of HSG and ASG values (69 and 22, respectively), as illustrated in Figure 5B. Figure 5C reports a heatmap showing the Z-score normalized expression of ASGs in all biological replicates for each variety. As a general observation, the elite grapevine varieties Cabernet sauvignon and Sauvignon blanc displayed relatively lower counts of ASG and HSG genes. Detailed lists of HSG and ASG for each grapevine variety can be found in Table S9.
As illustrated in the scatter plot depicted in Figure 5D, the expression levels of genes tend to decrease as specificity (tau) increases. This implies that many highly specific genes associated with a particular tissue might exhibit low expression levels within that tissue and nearly absent expression in others. Prioritizing genes based on both their expression and specificity proves advantageous, pinpointing a cohort of genes that display significant specificity to the tissue, boast ample expression levels (beneficial for lab experimentation), and demonstrate minimal expression in other tissues, thereby mitigating off-target effects. To accomplish this, we compiled quantile-normalized expression data and specificity metrics for all genes identified through tau analysis. Utilizing this dataset, we devised a scoring mechanism wherein each gene was assigned a score ranging from zero to two. This score was computed by summing its tau expression fraction value (τef) and its expression value normalized within the zero to one range. Subsequently, for each tissue, we identified the top 10 genes based on their score values, representing genes that exhibit both heightened expression and specificity (Figure 5E). While the complete list of these genes is available in Table S10, Table 1 explicitly showcases the foremost gene ranked highest for optimal specificity within each distinct tissue. The tau expression fractions of the 10 optimum genes are visually depicted in Supplementary Figure S2. The absolute specificity of certain genes for a particular genotype may be attributed to a mapping issue due to variations in the coding sequence affecting read mapping and consequently the RNA-seq results. Unfortunately, genomic sequences of the resistant varieties are not available; however, sequences of Cabernet sauvignon and Sauvignon blanc are accessible (www.grapegenomics.com). To dispel any doubts regarding mapping, the coding sequences (CDSs) of the top 10 significantly expressed genes of each resistant variety were conserved and their CDS in the elite varieties were aligned with those of the PN40024 reference genome. Although such a comparison does not provide information about the lack of expression of a particular gene in other resistant varieties, it can offer insights into the absence of expression in elite ones. For instance, the GDSL lipase gene (Vitvi07g02026), exclusively expressed in Cabernet volos, shows no polymorphism in the noble parent C. sauvignon or S. blanc. The lack of expression of this gene in these varieties is thus not attributable to polymorphism in the coding region but to a different level of basal expression dictated by distinct regulatory mechanisms.
Several genes among the top 10 ASGs in resistant varieties warrant particular attention. One such gene is Vitvi07g02026, which is specifically expressed in CC and encodes a GDSL lipase (VvGELP21). GDSL esterases/lipases are a subclass of lipolytic enzymes that play crucial roles in plant growth and development, stress response, and pathogen defense. Numerous studies have demonstrated the involvement of this gene class in defense mechanisms against biotic stresses. For instance, Ji et al. (2023) demonstrated the role of a GDSL esterase/lipase, GELP1, in defending apple leaves against Colletotrichum gloeosporioides infection [78]. Additionally, the overexpression of AtGDSL1 enhances resistance to Sclerotinia sclerotiorum in oilseed rape by modulating salicylic acid (SA)-dependent and jasmonic acid (JA)-dependent pathways. This overexpression leads to increased phosphatidic acid accumulation and the activation of downstream stress response pathways following Sclerotinia infection [79]. In Cabernet volos, Vitvi11g01637 encodes for a MLA R protein (MLA10), a class of protein involved in PM resistance in barley [80].
In Sauvignon nepis, several ASGs can be associated with pathogen resistance. Vitvi09g01181 encodes for a HcrVf1 protein which is known for its role in providing resistance against apple scab, a fungal disease [81]. The specific pression of an HcrVf2-like gene in SN berries implies that this gene may play a crucial role in the plant’s defense mechanisms, being involved in recognizing and responding to pathogen attacks and thereby enhancing the plant’s ability to resist diseases.
Vitvi18g02399 is an Avr9 elicitor response factor. Genes belonging to the same family were found to be upregulated in PM infection in infected fruits of the susceptible V. vinifera cv. Carignan [34] and were reported to regulate the hypersensitive response in grapevine leaves [82,83]. Another interesting gene found to be exclusively expressed in the berry skin of SN is Vitvi13g02352, encoding for the disease-resistance protein RGA2 (Resistance Gene Analog 2), which was found to play roles in conferring resistance to downy mildew [84]. As an NBS-LRR protein, RGA2 detects specific pathogen-associated molecular patterns or effector proteins, triggering a defense response. Upon recognition, RGA2 activates signaling cascades that produce reactive oxygen species (ROS) and calcium ions, leading to localized cell death (hypersensitive response) to limit pathogen spread. Additionally, it induces the production of antimicrobial compounds and strengthens cell walls, enhancing overall plant immunity.
In conclusion, beyond specific metabolic pathways such as those related to vesicular transport identified in resistant red varieties, we have identified a series of genes (for a complete list, see Table S10) associated with resistance responses and expressed exclusively in resistant varieties. These genes could become targets for genetic improvement techniques, such as cisgenics or genome editing, of elite susceptible varieties, pending functional characterization.

3. Materials and Methods

3.1. Plant Material and Sample Collection

To examine the transcriptional landscape of the exocarp in ripe berries from PIWI varieties and their noble parents, we conducted an RNA-seq analysis on 18 berry samples collected from field-grown plants, located in a graveyard in Fossalon di Grado (Udine, Italy). The plants studied included two red-berried PIWI varieties (Cabernet cortis, CC, and Cabernet volos, CV), two white-berried PIWI varieties (Sauvignon rytos, SR, and Sauvignon nepis, SN), and their respective noble parents: Cabernet sauvignon (CS) and Sauvignon blanc (SB). The berries used for the transcriptomic analysis were collected at full technological maturity, corresponding to E-L Stage 38 of the Eichhorn–Lorenz scale. Maturity was determined based on sugar content (°Brix) and the visual assessment of full ripeness, ensuring that the fruit had reached complete development before sampling. Moreover, no disease symptoms were observed on plants or berries during sampling. We collected samples in three biological replicates, with each replicate comprising 30 berries from three different field-grown plants per variety. The PIWI varieties considered in this study possess QTLs mediating partial resistance to grapevine DM and PM. The cultivar Cabernet cortis (CC) harbors the combination of Rpv3.3 and Rpv10 loci which concerns the resistance to DM, and Ren3 and Ren9 loci for the resistance to PM. It was obtained in 1984 by Norbert Becker by crossing Cabernet sauvignon and the Solaris variety. The DM and PM resistance loci were transmitted from a Solaris cv initially introgressed from V. amurensis, a wild species of the Asian Vitis gene pool. The cultivar Cabernet volos (CV) was obtained by the University of Udine and Institute of Applied Genetics (IGA) in 2002, crossing Cabernet sauvignon (CS) and Kozma 23-2—a Vitis interspecific crossing between Bianca and Sremski karlovci 77 4-5. The cultivars Sauvignon nepis (SN) and Sauvignon rytos (SR) were obtained by the University of Udine and Institute of Applied Genetics (IGA) in 2002, crossing Sauvignon blanc (SB) and Bianca—an interspecific hybrid whose lineage is made up of V. vinifera, V. rupestris, V. lincecumii, V. labrusca, and V. berlandieri.

3.2. RNA Purification, Library Preparation, and Sequencing

Berries were sectioned, and their skins were dissected using scalpels while being kept frozen with liquid nitrogen. The isolated exocarps were processed for RNA purification using the Spectrum™ Plant Total RNA Kit (Sigma-Aldrich, St. Louis, MO, USA) following the manufacturer’s instructions. For RNA extracts of low quality due to contaminants, further purification was performed using lithium chloride (LiCl) precipitation. This involved the addition of a 2.5 M LiCl solution to the extracted RNA, followed by overnight incubation at 4 °C. The next day, samples were centrifuged at 14,000 RPM for 15 min at 4 °C. The supernatant was carefully removed, and the samples were washed with 500 µL of 70% ethanol prepared in DEPC (diethyl pyrocarbonate) water. Nanodrop (Thermo Fisher Scientific, Wilmington, DE, USA) (quantity and purity) and Labchip GX (PerkinElmer, Waltham, MA, USA) (quality) were used to assess the quality and quantity of the isolated RNA. Then, cDNA libraries were constructed using ≥600 ng RNA per sample with the Hieff NGS Ultima Dual-mode mRNA Library Prep Kit for Illumina (Yeasen Biotechnology (Shanghai) Co., Ltd., Shanghai, China), following the manufacturer’s recommendations. The effective concentration of the library was determined to be at least 2 nM. The resultant sequencing libraries were then submitted for sequencing on the Illumina NovaSeq 6000 (Illumina, San Diego, CA, USA) platform, in PE150 mode. All raw reads were deposited in the NCBI SRA database with accession numbers SRRXXXX–SRRXXXXX.

3.3. RNA-seq Analysis

Read trimming and quality filtering were performed using the fastp v.0.20.1 software with default parameters [85]. Filtered reads were mapped on the 12x.v2 assembly of the PN40024 reference genome using the VCost.v3 annotation [51]. The reads alignment and counting were conducted using STAR software version 2.7.11a [86] using the following parameters: outMultimapperOrder = Random, outSAMmultNmax = 10, outWigStrand = Stranded, outWigNorm = None, quantMode = GeneCounts, twopassMode = Basic. The resulting counts were normalized with R/DESeq2 v.1.34.0 [87] with the method of the median of ratios.

3.4. Differential Gene Expression Analysis

Statistical analysis of the mapped reads was performed using the reads count as computed by the aligner STAR, normalized through the median of ratios method, and the R package DESeq2 [87]. Once the normalized reads counts and the DESeq dataset were created, a first filtering was performed to remove transcripts with less than 10 mapped reads before DEGs identification. Differentially expressed genes were identified considering two variables: (i) the grapevine variety typology (“noble” or “PIWI”), and (ii) the berry color (“White” and “Red”). The DEGs analysis was performed using the dds function from the DESeq2 package, and those with an adjusted p-value < 0.05. This first screening was performed for the overall dataset, and for the sub-datasets, focused on red and white varieties, which compared PIWI versus elite varieties separately (e.g., Cabernet volos vs. Cabernet sauvignon; Sauvignon nepis vs. Sauvignon blanc). Principal component analyses (PCAs) and heatmaps were computed through vst and pheatmap functions, respectively, to verify the concordance between replicates of each variety and to investigate the relationships between them. Upset plots to represent the private and shared DEGs between different groups were created using the upsetR package in RStudio [88].
The Gene Ontology enrichment analysis was performed using the R package gprofiler2 (v.0.2.3) using the default parameters and no custom genes background, since the gene excluded from the analyses had no sensible expression. The functional annotation of Vitis vinifera is present in the g:Profiler database under the name “vvinifera” [89].

3.5. K-Means Corrected Weighted Gene Co-Expression Network Analysis

The normalized expression indexes of all the samples were used to build a co-expression network and a first hierarchical clustering of the genes using the R package WGCNA v.1.72.5 [90]. The following parameters were passed to the function net: TOMType = “signed”, minModuleSize = 25, power = 30. The optimal power value was defined with the function pickSoftThreshold of the same package. To obtain two different levels of sensibility (i.e., smaller and more uniform modules or bigger and more generally descriptive ones), two different values of the parameter mergeCutHeight were adopted, 0.1 and 0.25, respectively. All the other parameters were kept to the default setting. The genes modules were then used as inputs for a k-means clustering algorithm from the package CoExpNets, wrapped by the function applyKMeans [61]. The parameters adopted were the defaults except for n.iteration = 50 and excludeGrey = TRUE. A linear correlation between the eigengene of each module and a series of logical vectors describing the studied character of the samples (i.e., the cultivar, the tolerance/resistance, and the color of the berries) was used to identify the modules associated and therefore of interest. The dedicated functions in the WGCNA package were used to further investigate these modules in terms of gene significance (GS, the linear correlation between the expression vector of the genes and the logical vector of the studied characters) and module membership (MM, a measure of centrality of the gene in the module co-expression network). A Gene Ontology enrichment analysis was performed on the genes composing the modules with the same procedure described in the previous paragraph.

3.6. Identification of Highly Specific Genes (HSG) and Absolutely Specific Genes (ASG) by Means of Tau Analysis

The index τ (tau) was used as an indicator of the tissue-specificity level of the expression of a gene [77]. Its range is from zero (for widely expressed genes) to one (for tissue/stage-specific genes). To calculate it, the tispec R package v0.99.0 was used. The reads counts were normalized using the TPM method with the R package bioinfokit v2.0.3. All the genes with less than 1 TPM in any tissue were then filtered out. To allow comparisons between tissues, the values were logarithm transformed, setting all the negative values to 0, and a quantile normalization, assigning to every gene a BIN value from 0 (lowest expression) to 10 (highest expression). The τ index of every gene was computed as follows:
τ = i = 1 N ( 1 x i ) N 1
where N is the number of tissues and xi is the expression value normalized by the highest expression. Genes expressed in a single tissue, with a τ of 1, are here referred as absolutely specific genes (ASGs), while we define as highly specific genes (HSGs) genes with a τ higher or equal to 0.85. A second index, a function of the normalized expression of the gene, was called τef and calculated as follows:
τ e f = τ q n m a x
where qn is the quantile normalized expression and max is the highest quantile normalized expression. For every gene in every tissue, the sum of τ and τef was called the τ-score and used as a proxy of the relevance of the gene in that tissue.

4. Conclusions

In this study, we used a transcriptomic approach to identify key metabolic pathways and specific genes that differentiate the basal defense mechanisms of resistant and susceptible grapevine varieties. Instead of focusing on infected versus uninfected plants, we analyzed the molecular baseline to uncover critical differences that prime certain varieties for enhanced tolerance. The berry skin, as the first barrier against pathogens, was chosen for this study due to its central role in defense processes, including pathogen recognition, signaling, and the activation of protective responses.
Our analyses revealed gene clusters associated with resistance and susceptibility, with a significant upregulation of vesicular transport mechanisms in red PIWI varieties. These processes, involving vesicle formation and cargo transport, have been linked to pathogen defense across various plant species. Their role in enhancing resistance likely lies in the more efficient trafficking of antimicrobial compounds or signaling molecules, providing a mechanistic insight into the improved resistance profile of these varieties.
We also identified a set of highly specific genes (HSGs) and absolutely specific genes (ASGs) uniquely expressed in PIWI varieties, likely playing key roles in basal defense. The exclusive presence of these genes in resistant varieties underscores their importance in adaptation to environmental stress and pathogen pressure. As PIWI varieties stem from breeding programs incorporating wild genomic contributions, these wild alleles likely harbor regulatory sequences or genes that modulate basal defense. This could explain the superior resistance observed in PIWI varieties.
While our study provides valuable insights into resistance mechanisms, further validation is needed. Some of the identified genes, especially those linked to vesicular transport and defense responses, present strong candidates for functional characterization through knock-out or knock-down experiments. Such studies would confirm their precise role in resistance and could lay the groundwork for future breeding strategies.
Additionally, our findings emphasize the complex nature of quantitative trait loci (QTLs) associated with resistance traits. Identifying regulatory networks and variety-specific gene expressions is crucial for understanding the genetic basis of resistance. The genes identified in this study could serve as promising candidates for Marker-Assisted Selection (MAS) in breeding programs, allowing for the introgression of resistance-associated traits into elite grapevine cultivars.
Looking forward, these results open avenues for future molecular studies and genetic manipulation aimed at producing highly resistant grapevine varieties. Beyond traditional breeding, these genes could be targeted by advanced biotechnological approaches such as assisted evolution or cisgenesis. These methods would enable the transfer of beneficial traits from PIWI varieties into elite cultivars without introducing foreign DNA. By adopting these techniques, the viticulture industry could reduce its dependence on chemical treatments and enhance crop resilience in the face of evolving pathogen threats.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/plants13182574/s1, Figure S1. Gene significance vs. module membership in selected modules; Figure S2. Optimum gene set top 10 genes; Table S1. Sequencing results read quality; Table S2. Mean normalized counts in all samples; Table S3. Differentially expressed genes in all pairwise comparisons; Table S4. Gene Set Enrichment Analysis on DEGs. Table S5. WGCNA module gene lists; Table S6. Module–genotype trait associations; Table S7. GSEA on WGCNA modules, cutoff 0.25; Table S8. GSEA on WGCNA modules cutoff 0.10; Table S9. Tau values of genes in all tissues considered in this study; Table S10. Score values of genes in all tissues considered in this study.

Author Contributions

Conceptualization, A.V. and M.L.; methodology, A.V.; software, F.S. and G.G.; data curation, F.S. and G.G.; writing—original draft preparation, A.V., F.S., G.G., G.M., A.D. and C.P.; writing—review and editing, A.V., F.P. and S.F.; visualization, A.V.; supervision, A.V., M.L., G.B. and A.C.; project administration, A.C.; funding acquisition, G.B. and A.C. All authors have read and agreed to the published version of the manuscript.

Funding

The study was carried out within the Agritech National Research Center and received funding from the European Union Next-Generation EU (PIANO NAZIONALE DI RIPRESA E RESILI-ENZA (PNRR) MISSIONE 4 COMPONENTE 2, INVESTIMENTO 1.4-D.D. 1032 17/06/2022, CN00000022, Spoke 4, PI: Gianni Barcaccia). This study represents an original paper within Spoke 4 (Task 4.1.1.) “Next-generation genotyping and -omics technologies for the molecular prediction of multiple resilient traits in crop plants”. Moreover, the research was possible thanks to the PRID 2022 project funded by the University of Padova titled “Le uve da ibridi resistenti di vite (varietà PIWI): caratterizzazione trascrittomica e proteomica ed effetti sulla stabilità e la potenziale allergenicità dei vini”.

Data Availability Statement

RNA-seq data have been submitted in the SRA database under the accession PRJNA1132586.

Acknowledgments

The authors would like to thank Sina Abdollahi who was of great help in the collection and dissection of the berry samples and in the pre-processing steps for RNA-seq analyses.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Guzmán-Ardiles, R.E.; Pegoraro, C.; da Maia, L.C.; Costa de Oliveira, A. Genetic Changes in the Genus Vitis and the Domestication of Vine. Front. Plant Sci. 2022, 13, 1019311. [Google Scholar] [CrossRef] [PubMed]
  2. Grassi, F.; Arroyo-Garcia, R. Editorial: Origins and Domestication of the Grape. Front. Plant Sci. 2020, 11, 1176. [Google Scholar] [CrossRef] [PubMed]
  3. Forneck, A.; Walker, M.A.; Blaich, R.; Yvon, M.; Leclant, F. Interaction of Phylloxera (Daktulosphaira Vitifoliae Fitch) with Grape (Vitis Spp.) in Simple Isolation Chambers. Am. J. Enol. Vitic. 2001, 52, 28–34. [Google Scholar] [CrossRef]
  4. Gessler, C.; Pertot, I.; Perazzolli, M. Plasmopara viticola: A Review of Knowledge on Downy Mildew of Grapevine and Effective Disease Management. Phytopathol. Mediterr. 2011, 50, 3–44. [Google Scholar]
  5. Gadoury, D.M.; Cadle-Davidson, L.; Wilcox, W.F.; Dry, I.B.; Seem, R.C.; Milgroom, M.G. Grapevine Powdery Mildew (Erysiphe Necator): A Fascinating System for the Study of the Biology, Ecology and Epidemiology of an Obligate Biotroph. Mol. Plant Pathol. 2012, 13, 1–16. [Google Scholar] [CrossRef] [PubMed]
  6. Caffi, T.; Gilardi, G.; Monchiero, M.; Rossi, V. Production and Release of Asexual Sporangia in Plasmopara viticola. Phytopathology 2013, 103, 64–73. [Google Scholar] [CrossRef]
  7. Mouafo-Tchinda, R.A.; Beaulieu, C.; Fall, M.L.; Carisse, O. Effect of Temperature on Aggressiveness of Plasmopara Viticola f. Sp. Aestivalis and P. Viticola f. Sp. Riparia from Eastern Canada. Can. J. Plant Pathol. 2021, 43, 73–87. [Google Scholar] [CrossRef]
  8. Carroll, J.E.; Wilcox, W.F. Effects of Humidity on the Development of Grapevine Powdery Mildew. Phytopathology 2003, 93, 1137–1144. [Google Scholar] [CrossRef]
  9. Nicolopoulou-Stamati, P.; Maipas, S.; Kotampasi, C.; Stamatis, P.; Hens, L. Chemical Pesticides and Human Health: The Urgent Need for a New Concept in Agriculture. Front. Public Health 2016, 4, 148. [Google Scholar] [CrossRef]
  10. Possamai, T.; Wiedemann-Merdinoglu, S. Phenotyping for QTL Identification: A Case Study of Resistance to Plasmopara Viticola and Erysiphe Necator in Grapevine. Front. Plant Sci. 2022, 13, 930954. [Google Scholar] [CrossRef]
  11. Marguerit, E.; Boury, C.; Manicki, A.; Donnart, M.; Butterlin, G.; Némorin, A.; Wiedemann-Merdinoglu, S.; Merdinoglu, D.; Ollat, N.; Decroocq, S. Genetic Dissection of Sex Determinism, Inflorescence Morphology and Downy Mildew Resistance in Grapevine. Theor. Appl. Genet. 2009, 118, 1261–1278. [Google Scholar] [CrossRef] [PubMed]
  12. Pap, D.; Riaz, S.; Dry, I.B.; Jermakow, A.; Tenscher, A.C.; Cantu, D.; Oláh, R.; Walker, M.A. Identification of Two Novel Powdery Mildew Resistance Loci, Ren6 and Ren7, from the Wild Chinese Grape Species Vitis Piasezkii. BMC Plant Biol. 2016, 16, 170. [Google Scholar] [CrossRef] [PubMed]
  13. Sapkota, S.; Chen, L.-L.; Yang, S.; Hyma, K.E.; Cadle-Davidson, L.; Hwang, C.-F. Construction of a High-Density Linkage Map and QTL Detection of Downy Mildew Resistance in Vitis Aestivalis-Derived ‘Norton’. Theor. Appl. Genet. 2019, 132, 137–147. [Google Scholar] [CrossRef]
  14. Lin, H.; Leng, H.; Guo, Y.; Kondo, S.; Zhao, Y.; Shi, G.; Guo, X. QTLs and Candidate Genes for Downy Mildew Resistance Conferred by Interspecific Grape (V. Vinifera L. × V. Amurensis Rupr.) Crossing. Sci. Hortic. 2019, 244, 200–207. [Google Scholar] [CrossRef]
  15. Pauquet, J.; Bouquet, A.; This, P.; Adam-Blondon, A.-F. Establishment of a Local Map of AFLP Markers around the Powdery Mildew Resistance Gene Run1 in Grapevine and Assessment of Their Usefulness for Marker Assisted Selection. Theor. Appl. Genet. 2001, 103, 1201–1210. [Google Scholar] [CrossRef]
  16. Riaz, S.; Tenscher, A.C.; Ramming, D.W.; Walker, M.A. Using a Limited Mapping Strategy to Identify Major QTLs for Resistance to Grapevine Powdery Mildew (Erysiphe Necator) and Their Use in Marker-Assisted Breeding. Theor. Appl. Genet. 2011, 122, 1059–1073. [Google Scholar] [CrossRef]
  17. Mahanil, S.; Ramming, D.; Cadle-Davidson, M.; Owens, C.; Garris, A.; Myles, S.; Cadle-Davidson, L. Development of Marker Sets Useful in the Early Selection of Ren4 Powdery Mildew Resistance and Seedlessness for Table and Raisin Grape Breeding. Theor. Appl. Genet. 2012, 124, 23–33. [Google Scholar] [CrossRef]
  18. Karn, A.; Zou, C.; Brooks, S.; Fresnedo-Ramírez, J.; Gabler, F.; Sun, Q.; Ramming, D.; Naegele, R.; Ledbetter, C.; Cadle-Davidson, L. Discovery of the REN11 Locus From Vitis Aestivalis for Stable Resistance to Grapevine Powdery Mildew in a Family Segregating for Several Unstable and Tissue-Specific Quantitative Resistance Loci. Front. Plant Sci. 2021, 12, 733899. [Google Scholar] [CrossRef]
  19. Welter, L.J.; Göktürk-Baydar, N.; Akkurt, M.; Maul, E.; Eibach, R.; Töpfer, R.; Zyprian, E.M. Genetic Mapping and Localization of Quantitative Trait Loci Affecting Fungal Disease Resistance and Leaf Morphology in Grapevine (Vitis vinifera L.). Mol. Breed. 2007, 20, 359–374. [Google Scholar] [CrossRef]
  20. Bellin, D.; Peressotti, E.; Merdinoglu, D.; Wiedemann-Merdinoglu, S.; Adam-Blondon, A.-F.; Cipriani, G.; Morgante, M.; Testolin, R.; Di Gaspero, G. Resistance to Plasmopara Viticola in Grapevine ‘Bianca’ Is Controlled by a Major Dominant Gene Causing Localised Necrosis at the Infection Site. Theor. Appl. Genet. 2009, 120, 163–176. [Google Scholar] [CrossRef]
  21. Venuti, S.; Copetti, D.; Foria, S.; Falginella, L.; Hoffmann, S.; Bellin, D.; Cindrić, P.; Kozma, P.; Scalabrin, S.; Morgante, M.; et al. Historical Introgression of the Downy Mildew Resistance Gene Rpv12 from the Asian Species Vitis Amurensis into Grapevine Varieties. PLoS ONE 2013, 8, e61228. [Google Scholar] [CrossRef] [PubMed]
  22. Vezzulli, S.; Malacarne, G.; Masuero, D.; Vecchione, A.; Dolzani, C.; Goremykin, V.; Mehari, Z.H.; Banchi, E.; Velasco, R.; Stefanini, M.; et al. The Rpv3-3 Haplotype and Stilbenoid Induction Mediate Downy Mildew Resistance in a Grapevine Interspecific Population. Front. Plant Sci. 2019, 10, 234. [Google Scholar] [CrossRef] [PubMed]
  23. Ricciardi, V.; Crespan, M.; Maddalena, G.; Migliaro, D.; Brancadoro, L.; Maghradze, D.; Failla, O.; Toffolatti, S.L.; De Lorenzis, G. Novel Loci Associated with Resistance to Downy and Powdery Mildew in Grapevine. Front. Plant Sci. 2024, 15, 1386225. [Google Scholar] [CrossRef] [PubMed]
  24. Sargolzaei, M.; Maddalena, G.; Bitsadze, N.; Maghradze, D.; Bianco, P.A.; Failla, O.; Toffolatti, S.L.; De Lorenzis, G. Rpv29, Rpv30 and Rpv31: Three Novel Genomic Loci Associated With Resistance to Plasmopara Viticola in Vitis Vinifera. Front. Plant Sci. 2020, 11, 562432. [Google Scholar] [CrossRef] [PubMed]
  25. Di Gaspero, G.; Copetti, D.; Coleman, C.; Castellarin, S.D.; Eibach, R.; Kozma, P.; Lacombe, T.; Gambetta, G.; Zvyagin, A.; Cindrić, P.; et al. Selective Sweep at the Rpv3 Locus during Grapevine Breeding for Downy Mildew Resistance. Theor. Appl. Genet. 2012, 124, 277–286. [Google Scholar] [CrossRef]
  26. Peressotti, E.; Dolzani, C.; Poles, L.; Banchi, E.; Stefanini, M.; Salamini, F.; Velasco, R.; Vezzulli, S.; Riaz, S.; Walker, M.A.; et al. A first pedigree-based analysis (PBA) approach for the dissection of disease resistance traits in grapevine hybrids. Acta Hortic. 2015, 113–121. [Google Scholar] [CrossRef]
  27. Vezzulli, S.; Gramaje, D.; Tello, J.; Gambino, G.; Bettinelli, P.; Pirrello, C.; Schwandner, A.; Barba, P.; Angelini, E.; Anfora, G.; et al. Genomic Designing for Biotic Stress Resistant Grapevine. In Genomic Designing for Biotic Stress Resistant Fruit Crops; Kole, C., Ed.; Springer International Publishing: Cham, Switzerland, 2022; pp. 87–255. ISBN 978-3-030-91801-9. [Google Scholar]
  28. Tello, J.; Ibáñez, J. Review: Status and Prospects of Association Mapping in Grapevine. Plant Sci. 2023, 327, 111539. [Google Scholar] [CrossRef]
  29. Pirrello, C.; Magon, G.; Palumbo, F.; Farinati, S.; Lucchin, M.; Barcaccia, G.; Vannozzi, A. Past, Present, and Future of Genetic Strategies to Control Tolerance to the Main Fungal and Oomycete Pathogens of Grapevine. J. Exp. Bot. 2023, 74, 1309–1330. [Google Scholar] [CrossRef]
  30. Zanghelini, J.A.; Bogo, A.; Dal Vesco, L.L.; Gomes, B.R.; Mecabô, C.V.; Herpich, C.H.; Welter, L.J. Response of PIWI Grapevine Cultivars to Downy Mildew in Highland Region of Southern Brazil. Eur. J. Plant Pathol. 2019, 154, 1051–1058. [Google Scholar] [CrossRef]
  31. Wu, J.; Zhang, Y.; Zhang, H.; Huang, H.; Folta, K.M.; Lu, J. Whole Genome Wide Expression Profiles of Vitis Amurensisgrape Responding to Downy Mildew by Using Solexa Sequencing Technology. BMC Plant Biol. 2010, 10, 234. [Google Scholar] [CrossRef]
  32. Legay, G.; Marouf, E.; Berger, D.; Neuhaus, J.-M.; Mauch-Mani, B.; Slaughter, A. Identification of Genes Expressed during the Compatible Interaction of Grapevine with Plasmopara Viticola through Suppression Subtractive Hybridization (SSH). Eur. J. Plant Pathol. 2011, 129, 281–301. [Google Scholar] [CrossRef]
  33. Vannozzi, A.; Dry, I.B.; Fasoli, M.; Zenoni, S.; Lucchin, M. Genome-Wide Analysis of the Grapevine Stilbene Synthase Multigenic Family: Genomic Organization and Expression Profiles upon Biotic and Abiotic Stresses. BMC Plant Biol. 2012, 12, 130. [Google Scholar] [CrossRef] [PubMed]
  34. Pimentel, D.; Amaro, R.; Erban, A.; Mauri, N.; Soares, F.; Rego, C.; Martínez-Zapater, J.M.; Mithöfer, A.; Kopka, J.; Fortes, A.M. Transcriptional, Hormonal, and Metabolic Changes in Susceptible Grape Berries under Powdery Mildew Infection. J. Exp. Bot. 2021, 72, 6544–6569. [Google Scholar] [CrossRef] [PubMed]
  35. Figueiredo, A.; Fortes, A.M.; Ferreira, S.; Sebastiana, M.; Choi, Y.H.; Sousa, L.; Acioli-Santos, B.; Pessoa, F.; Verpoorte, R.; Pais, M.S. Transcriptional and Metabolic Profiling of Grape (Vitis vinifera L.) Leaves Unravel Possible Innate Resistance against Pathogenic Fungi. J. Exp. Bot. 2008, 59, 3371–3381. [Google Scholar] [CrossRef] [PubMed]
  36. Figueiredo, A.; Monteiro, F.; Fortes, A.M.; Bonow-Rex, M.; Zyprian, E.; Sousa, L.; Pais, M.S. Cultivar-Specific Kinetics of Gene Induction during Downy Mildew Early Infection in Grapevine. Funct. Integr. Genom. 2012, 12, 379–386. [Google Scholar] [CrossRef] [PubMed]
  37. Polesani, M.; Bortesi, L.; Ferrarini, A.; Zamboni, A.; Fasoli, M.; Zadra, C.; Lovato, A.; Pezzotti, M.; Delledonne, M.; Polverari, A. General and Species-Specific Transcriptional Responses to Downy Mildew Infection in a Susceptible (Vitis vinifera) and a Resistant (V. riparia) Grapevine Species. BMC Genom. 2010, 11, 117. [Google Scholar] [CrossRef] [PubMed]
  38. Malacarne, G.; Vrhovsek, U.; Zulini, L.; Cestaro, A.; Stefanini, M.; Mattivi, F.; Delledonne, M.; Velasco, R.; Moser, C. Resistance to Plasmopara Viticolain a Grapevine Segregating Population Is Associated with Stilbenoid Accumulation and with Specific Host Transcriptional Responses. BMC Plant Biol. 2011, 11, 114. [Google Scholar] [CrossRef]
  39. Amrine, K.C.H.; Blanco-Ulate, B.; Riaz, S.; Pap, D.; Jones, L.; Figueroa-Balderas, R.; Walker, M.A.; Cantu, D. Comparative Transcriptomics of Central Asian Vitis Vinifera Accessions Reveals Distinct Defense Strategies against Powdery Mildew. Hortic. Res. 2015, 2, 15037. [Google Scholar] [CrossRef]
  40. Li, X.; Wu, J.; Yin, L.; Zhang, Y.; Qu, J.; Lu, J. Comparative Transcriptome Analysis Reveals Defense-Related Genes and Pathways against Downy Mildew in Vitis amurensis Grapevine. Plant Physiol. Biochem. 2015, 95, 1–14. [Google Scholar] [CrossRef]
  41. Eisenmann, B.; Czemmel, S.; Ziegler, T.; Buchholz, G.; Kortekamp, A.; Trapp, O.; Rausch, T.; Dry, I.; Bogs, J. Rpv3–1 Mediated Resistance to Grapevine Downy Mildew Is Associated with Specific Host Transcriptional Responses and the Accumulation of Stilbenes. BMC Plant Biol. 2019, 19, 343. [Google Scholar] [CrossRef]
  42. Perazzolli, M.; Moretto, M.; Fontana, P.; Ferrarini, A.; Velasco, R.; Moser, C.; Delledonne, M.; Pertot, I. Downy Mildew Resistance Induced by Trichoderma Harzianum T39 in Susceptible Grapevines Partially Mimics Transcriptional Changes of Resistant Genotypes. BMC Genom. 2012, 13, 660. [Google Scholar] [CrossRef] [PubMed]
  43. Pagliarani, C.; Moine, A.; Chitarra, W.; Meloni, G.R.; Abbà, S.; Nerva, L.; Pugliese, M.; Gullino, M.L.; Gambino, G. The Molecular Priming of Defense Responses Is Differently Regulated in Grapevine Genotypes Following Elicitor Application against Powdery Mildew. Int. J. Mol. Sci. 2020, 21, 6776. [Google Scholar] [CrossRef] [PubMed]
  44. Lund, S.T.; Bohlmann, J. The Molecular Basis for Wine Grape Quality—A Volatile Subject. Science 2006, 311, 804–805. [Google Scholar] [CrossRef] [PubMed]
  45. Hardie, W.J.; Aggenbach, S.J.; Jaudzems, V.G. The Plastids of the Grape Pericarp and Their Significance in Isoprenoid Synthesis. Aust. J. Grape Wine Res. 1996, 2, 144–154. [Google Scholar] [CrossRef]
  46. Hichri, I.; Barrieu, F.; Bogs, J.; Kappel, C.; Delrot, S.; Lauvergeat, V. Recent Advances in the Transcriptional Regulation of the Flavonoid Biosynthetic Pathway. J. Exp. Bot. 2011, 62, 2465–2483. [Google Scholar] [CrossRef]
  47. Ali, M.B.; Howard, S.; Chen, S.; Wang, Y.; Yu, O.; Kovacs, L.G.; Qiu, W. Berry Skin Development in Norton Grape: Distinct Patterns of Transcriptional Regulation and Flavonoid Biosynthesis. BMC Plant Biol. 2011, 11, 7. [Google Scholar] [CrossRef]
  48. Muganu, M.; Bellincontro, A.; Barnaba, F.E.; Paolocci, M.; Bignami, C.; Gambellini, G.; Mencarelli, F. Influence of Bunch Position in the Canopy on Berry Epicuticular Wax during Ripening and on Weight Loss during Postharvest Dehydration. Am. J. Enol. Vitic. 2011, 62, 91–98. [Google Scholar] [CrossRef]
  49. Chang, B.-M.; Keller, M. Cuticle and Skin Cell Walls Have Common and Unique Roles in Grape Berry Splitting. Hortic. Res. 2021, 8, 168. [Google Scholar] [CrossRef]
  50. Wada, H.; Matthews, M.A.; Shackel, K.A. Seasonal Pattern of Apoplastic Solute Accumulation and Loss of Cell Turgor during Ripening of Vitis Vinifera Fruit under Field Conditions. J. Exp. Bot. 2009, 60, 1773–1781. [Google Scholar] [CrossRef]
  51. Canaguier, A.; Grimplet, J.; Di Gaspero, G.; Scalabrin, S.; Duchêne, E.; Choisne, N.; Mohellibi, N.; Guichard, C.; Rombauts, S.; Le Clainche, I.; et al. A New Version of the Grapevine Reference Genome Assembly (12X.v2) and of Its Annotation (VCost.V3). Genom. Data 2017, 14, 56–62. [Google Scholar] [CrossRef]
  52. Brighenti, E.; Souza, A.L.K.; Brighenti, A.F.; Stefanini, M.; Trapp, O.; Gardin, J.P.P.; Caliari, V.; Dalbó, M.A.; Welter, L.J. Field Performance of Five White Pilzwiderstandsfahige (PIWI) Cultivars in the South of Brazil. Acta Hortic. 2019, 115–122. [Google Scholar] [CrossRef]
  53. Bavaresco, L.; Vercesi, A.; Belvini, P.; Dalla Costa, L.; Fogal, J.; Marcon, L.; Masaro, L.; Pascarella, G.; Pastro, M.; Pizzolato, A.; et al. Agronomic Performance of 21 New Disease Resistant Winegrape Varieties Grown in Northeast Italy. VITIS J. Grapevine Res. 2023, 81–87. [Google Scholar] [CrossRef]
  54. Wang, L.; Li, Y.; Jin, X.; Liu, L.; Dai, X.; Liu, Y.; Zhao, L.; Zheng, P.; Wang, X.; Liu, Y.; et al. Floral Transcriptomes Reveal Gene Networks in Pineapple Floral Growth and Fruit Development. Commun. Biol. 2020, 3, 500. [Google Scholar] [CrossRef] [PubMed]
  55. Shahan, R.; Zawora, C.; Wight, H.; Sittmann, J.; Wang, W.; Mount, S.M.; Liu, Z. Consensus Coexpression Network Analysis Identifies Key Regulators of Flower and Fruit Development in Wild Strawberry. Plant Physiol. 2018, 178, 202–216. [Google Scholar] [CrossRef] [PubMed]
  56. Li, X.; Liu, L.; Ming, M.; Hu, H.; Zhang, M.; Fan, J.; Song, B.; Zhang, S.; Wu, J. Comparative Transcriptomic Analysis Provides Insight into the Domestication and Improvement of Pear (P. pyrifolia) Fruit. Plant Physiol. 2019, 180, 435–452. [Google Scholar] [CrossRef] [PubMed]
  57. Guo, D.-L.; Wang, Z.-G.; Pei, M.-S.; Guo, L.-L.; Yu, Y.-H. Transcriptome Analysis Reveals Mechanism of Early Ripening in Kyoho Grape with Hydrogen Peroxide Treatment. BMC Genom. 2020, 21, 784. [Google Scholar] [CrossRef]
  58. Vannozzi, A.; Palumbo, F.; Magon, G.; Lucchin, M.; Barcaccia, G. The Grapevine (Vitis vinifera L.) Floral Transcriptome in Pinot Noir Variety: Identification of Tissue-Related Gene Networks and Whorl-Specific Markers in Pre- and Post-Anthesis Phases. Hortic. Res. 2021, 8, 200. [Google Scholar] [CrossRef]
  59. Zheng, H.; Dong, Y.; Nong, H.; Huang, L.; Liu, J.; Yu, X.; Zhang, Y.; Yang, L.; Hong, B.; Wang, W.; et al. VvSUN May Act in the Auxin Pathway to Regulate Fruit Shape in Grape. Hortic. Res. 2022, 9, uhac200. [Google Scholar] [CrossRef]
  60. Wu, J.; Zhang, F.; Liu, G.; Abudureheman, R.; Bai, S.; Wu, X.; Zhang, C.; Ma, Y.; Wang, X.; Zha, Q.; et al. Transcriptome and Coexpression Network Analysis Reveals Properties and Candidate Genes Associated with Grape (Vitis vinifera L.) Heat Tolerance. Front. Plant Sci. 2023, 14, 1270933. [Google Scholar] [CrossRef]
  61. Botía, J.A.; Vandrovcova, J.; Forabosco, P.; Guelfi, S.; D’Sa, K.; Hardy, J.; Lewis, C.M.; Ryten, M.; Weale, M.E. The United Kingdom Brain Expression Consortium An Additional K-Means Clustering Step Improves the Biological Features of WGCNA Gene Co-Expression Networks. BMC Syst. Biol. 2017, 11, 47. [Google Scholar] [CrossRef]
  62. Wang, P.; Siao, W.; Zhao, X.; Arora, D.; Wang, R.; Eeckhout, D.; Van Leene, J.; Kumar, R.; Houbaert, A.; De Winne, N.; et al. Adaptor Protein Complex Interaction Map in Arabidopsis Identifies P34 as a Common Stability Regulator. Nat. Plants 2023, 9, 355–371. [Google Scholar] [CrossRef] [PubMed]
  63. Beck, M.; Heard, W.; Mbengue, M.; Robatzek, S. The INs and OUTs of Pattern Recognition Receptors at the Cell Surface. Curr. Opin. Plant Biol. 2012, 15, 367–374. [Google Scholar] [CrossRef] [PubMed]
  64. Fan, L.; Li, R.; Pan, J.; Ding, Z.; Lin, J. Endocytosis and Its Regulation in Plants. Trends Plant Sci. 2015, 20, 388–397. [Google Scholar] [CrossRef] [PubMed]
  65. Oliveira-Garcia, E.; Tamang, T.M.; Park, J.; Dalby, M.; Martin-Urdiroz, M.; Rodriguez Herrero, C.; Vu, A.H.; Park, S.; Talbot, N.J.; Valent, B. Clathrin-Mediated Endocytosis Facilitates the Internalization of Magnaporthe Oryzae Effectors into Rice Cells. Plant Cell 2023, 35, 2527–2551. [Google Scholar] [CrossRef] [PubMed]
  66. Mbengue, M.; Bourdais, G.; Gervasi, F.; Beck, M.; Zhou, J.; Spallek, T.; Bartels, S.; Boller, T.; Ueda, T.; Kuhn, H.; et al. Clathrin-Dependent Endocytosis Is Required for Immunity Mediated by Pattern Recognition Receptor Kinases. Proc. Natl. Acad. Sci. USA 2016, 113, 11034–11039. [Google Scholar] [CrossRef]
  67. Spallek, T.; Beck, M.; Khaled, S.B.; Salomon, S.; Bourdais, G.; Schellmann, S.; Robatzek, S. ESCRT-I Mediates FLS2 Endosomal Sorting and Plant Immunity. PLOS Genet. 2013, 9, e1004035. [Google Scholar] [CrossRef]
  68. Postma, J.; Liebrand, T.W.H.; Bi, G.; Evrard, A.; Bye, R.R.; Mbengue, M.; Kuhn, H.; Joosten, M.H.A.J.; Robatzek, S. Avr4 Promotes Cf-4 Receptor-like Protein Association with the BAK1/SERK3 Receptor-like Kinase to Initiate Receptor Endocytosis and Plant Immunity. New Phytol. 2016, 210, 627–642. [Google Scholar] [CrossRef]
  69. Chung, K.P.; Zeng, Y.; Jiang, L. COPII Paralogs in Plants: Functional Redundancy or Diversity? Trends Plant Sci. 2016, 21, 758–769. [Google Scholar] [CrossRef]
  70. Kwon, C.; Bednarek, P.; Schulze-Lefert, P. Secretory Pathways in Plant Immune Responses. Plant Physiol. 2008, 147, 1575–1583. [Google Scholar] [CrossRef]
  71. Kulich, I.; Žárský, V. Autophagy-Related Direct Membrane Import from ER/Cytoplasm into the Vacuole or Apoplast: A Hidden Gateway Also for Secondary Metabolites and Phytohormones? IJMS 2014, 15, 7462–7474. [Google Scholar] [CrossRef]
  72. Ruano, G.; Scheuring, D. Plant Cells under Attack: Unconventional Endomembrane Trafficking during Plant Defense. Plants 2020, 9, 389. [Google Scholar] [CrossRef] [PubMed]
  73. Kurt-Celebi, A.; Colak, N.; Hayirlioglu-Ayaz, S.; Kostadinović Veličkovska, S.; Ilieva, F.; Esatbeyoglu, T.; Ayaz, F.A. Accumulation of Phenolic Compounds and Antioxidant Capacity during Berry Development in Black ‘Isabel’ Grape (Vitis vinifera L. x Vitis labrusca L.). Molecules 2020, 25, 3845. [Google Scholar] [CrossRef] [PubMed]
  74. Phua, S.Y.; De Smet, B.; Remacle, C.; Chan, K.X.; Van Breusegem, F. Reactive Oxygen Species and Organellar Signaling. J. Exp. Bot. 2021, 72, 5807–5824. [Google Scholar] [CrossRef] [PubMed]
  75. Cao, J.; Wang, C.; Hao, N.; Fujiwara, T.; Wu, T. Endoplasmic Reticulum Stress and Reactive Oxygen Species in Plants. Antioxidants 2022, 11, 1240. [Google Scholar] [CrossRef] [PubMed]
  76. Kısaca, G.; Gazioglu Sensoy, R.I. Phenolic Contents, Organicacids and Antioxidant Capacities of Twenty Grape (Vitis vinifera L.) Cultivars Having Different Berry Colors. Food Meas. 2023, 17, 1354–1370. [Google Scholar] [CrossRef]
  77. Kryuchkova-Mostacci, N.; Robinson-Rechavi, M. A Benchmark of Gene Expression Tissue-Specificity Metrics. Brief. Bioinform. 2017, 18, 205–214. [Google Scholar] [CrossRef]
  78. Ji, Z.; Wang, M.; Zhang, S.; Du, Y.; Cong, J.; Yan, H.; Guo, H.; Xu, B.; Zhou, Z. GDSL Esterase/Lipase GELP1 Involved in the Defense of Apple Leaves against Colletotrichum Gloeosporioides Infection. Int. J. Mol. Sci. 2023, 24, 10343. [Google Scholar] [CrossRef]
  79. Ding, L.-N.; Li, M.; Guo, X.-J.; Tang, M.-Q.; Cao, J.; Wang, Z.; Liu, R.; Zhu, K.-M.; Guo, L.; Liu, S.-Y.; et al. Arabidopsis GDSL1 Overexpression Enhances Rapeseed Sclerotinia Sclerotiorum Resistance and the Functional Identification of Its Homolog in Brassica Napus. Plant Biotechnol. J. 2020, 18, 1255–1270. [Google Scholar] [CrossRef]
  80. Wei, F.; Wing, R.A.; Wise, R.P. Genome Dynamics and Evolution of the Mla (Powdery Mildew) Resistance Locus in Barley. Plant Cell 2002, 14, 1903–1917. [Google Scholar] [CrossRef]
  81. Joshi, S.G.; Schaart, J.G.; Groenwold, R.; Jacobsen, E.; Schouten, H.J.; Krens, F.A. Functional Analysis and Expression Profiling of HcrVf1 and HcrVf2 for Development of Scab Resistant Cisgenic and Intragenic Apples. Plant Mol. Biol. 2011, 75, 579–591. [Google Scholar] [CrossRef]
  82. Rowland, O.; Ludwig, A.A.; Merrick, C.J.; Baillieul, F.; Tracy, F.E.; Durrant, W.E.; Fritz-Laylin, L.; Nekrasov, V.; Sjölander, K.; Yoshioka, H.; et al. Functional Analysis of Avr9/Cf-9 Rapidly Elicited Genes Identifies a Protein Kinase, ACIK1, That Is Essential for Full Cf-9–Dependent Disease Resistance in Tomato. Plant Cell 2005, 17, 295–310. [Google Scholar] [CrossRef]
  83. Toth, Z.; Winterhagen, P.; Kalapos, B.; Su, Y.; Kovacs, L.; Kiss, E. Expression of a Grapevine NAC Transcription Factor Gene Is Induced in Response to Powdery Mildew Colonization in Salicylic Acid-Independent Manner. Sci. Rep. 2016, 6, 30825. [Google Scholar] [CrossRef]
  84. Wang, P.; Liu, C.; Wang, D.; Liang, C.; Zhao, K.; Fan, J. Isolation of Resistance Gene Analogs from Grapevine Resistant to Downy Mildew. Sci. Hortic. 2013, 150, 326–333. [Google Scholar] [CrossRef]
  85. Chen, S. Ultrafast One-Pass FASTQ Data Preprocessing, Quality Control, and Deduplication Using Fastp. iMeta 2023, 2, e107. [Google Scholar] [CrossRef] [PubMed]
  86. Dobin, A.; Davis, C.A.; Schlesinger, F.; Drenkow, J.; Zaleski, C.; Jha, S.; Batut, P.; Chaisson, M.; Gingeras, T.R. STAR: Ultrafast Universal RNA-Seq Aligner. Bioinformatics 2013, 29, 15–21. [Google Scholar] [CrossRef] [PubMed]
  87. Love, M.I.; Huber, W.; Anders, S. Moderated Estimation of Fold Change and Dispersion for RNA-Seq Data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef] [PubMed]
  88. Conway, J.R.; Lex, A.; Gehlenborg, N. UpSetR: An R Package for the Visualization of Intersecting Sets and Their Properties. Bioinformatics 2017, 33, 2938–2940. [Google Scholar] [CrossRef]
  89. Kolberg, L.; Raudvere, U.; Kuzmin, I.; Vilo, J.; Peterson, H. Gprofiler2—An R Package for Gene List Functional Enrichment Analysis and Namespace Conversion Toolset g:Profiler. F1000Research 2020, 9, ELIXIR-709. [Google Scholar] [CrossRef]
  90. Langfelder, P.; Horvath, S. WGCNA: An R Package for Weighted Correlation Network Analysis. BMC Bioinform. 2008, 9, 559. [Google Scholar] [CrossRef]
Figure 1. (A) Correlation matrix heatmap showing the Euclidean distance between samples based on normalized data obtained from 18 RNA-seq samples constituted of berry skin tissues of the CC, CV, CS, SR, SN, and SB varieties in the ripening (R) phase. A darker color indicates a stronger correlation. (B) PCA on normalized data obtained from 18 RNA-seq samples. Colors indicate different varieties considered. (C) The histogram shows the number of upregulated and downregulated DEGs in white and red PIWI varieties compared to their respective noble parents (SB for white and CS for red). It includes both cumulative comparisons of all PIWI varieties of the same color against their parental variety, as well as individual comparisons (e.g., SR vs. SB). (D) Upset plots visualizing the intersections amongst different groups of DEGs identified in pairwise comparisons. Single points indicate a private DEG identified in each group, whereas 2 to n dot plots indicate DEGs shared by 2 to n groups.
Figure 1. (A) Correlation matrix heatmap showing the Euclidean distance between samples based on normalized data obtained from 18 RNA-seq samples constituted of berry skin tissues of the CC, CV, CS, SR, SN, and SB varieties in the ripening (R) phase. A darker color indicates a stronger correlation. (B) PCA on normalized data obtained from 18 RNA-seq samples. Colors indicate different varieties considered. (C) The histogram shows the number of upregulated and downregulated DEGs in white and red PIWI varieties compared to their respective noble parents (SB for white and CS for red). It includes both cumulative comparisons of all PIWI varieties of the same color against their parental variety, as well as individual comparisons (e.g., SR vs. SB). (D) Upset plots visualizing the intersections amongst different groups of DEGs identified in pairwise comparisons. Single points indicate a private DEG identified in each group, whereas 2 to n dot plots indicate DEGs shared by 2 to n groups.
Plants 13 02574 g001
Figure 2. K-means-corrected WGCNA. (A) Cluster dendrogram of module eigengenes. Branches of the dendrogram group together eigengenes that are positively correlated. The merge threshold was set to 0.25: modules under this value were merged due to their similarity in expression profiles. (B) Bar graph showing the distribution of genes over the twenty-six modules identified. (C) Module-variety/trait association analysis. The heatmap shows the correlation between modules and varieties/traits. Each row corresponds to a module, whereas each column corresponds to a specific trait. The correlation coefficient between a given module and tissue type is indicated by the color of the cell at the row–column intersection and by the text inside the cells (squared boxes indicate significant p-values). Red and blue indicate positive and negative correlations, respectively. CC, Cabernet cortis; SN, Sauvignon nepis; SR, Sauvignon rytos; CV, Cabernet volos; SB, Sauvignon blanc; CS, Cabernet sauvignon; T/S, tolerance/susceptibility; GC, grape color. (D) Scatterplots of gene significance (GS) vs. module membership (MM) in the brown module associated with Cabernet cortis (CC). Genes highly significantly associated with a trait are often also the most important (central) elements of modules associated with the trait. (E) Heatmap visualizing gene expression within the brown module across all biological replicates of the six considered varieties, normalized using Z-scores.
Figure 2. K-means-corrected WGCNA. (A) Cluster dendrogram of module eigengenes. Branches of the dendrogram group together eigengenes that are positively correlated. The merge threshold was set to 0.25: modules under this value were merged due to their similarity in expression profiles. (B) Bar graph showing the distribution of genes over the twenty-six modules identified. (C) Module-variety/trait association analysis. The heatmap shows the correlation between modules and varieties/traits. Each row corresponds to a module, whereas each column corresponds to a specific trait. The correlation coefficient between a given module and tissue type is indicated by the color of the cell at the row–column intersection and by the text inside the cells (squared boxes indicate significant p-values). Red and blue indicate positive and negative correlations, respectively. CC, Cabernet cortis; SN, Sauvignon nepis; SR, Sauvignon rytos; CV, Cabernet volos; SB, Sauvignon blanc; CS, Cabernet sauvignon; T/S, tolerance/susceptibility; GC, grape color. (D) Scatterplots of gene significance (GS) vs. module membership (MM) in the brown module associated with Cabernet cortis (CC). Genes highly significantly associated with a trait are often also the most important (central) elements of modules associated with the trait. (E) Heatmap visualizing gene expression within the brown module across all biological replicates of the six considered varieties, normalized using Z-scores.
Plants 13 02574 g002
Figure 3. Modules contemporaneously associated with both tolerance/susceptibility and one or more grapevine varieties analyzed. (A) Table showing the orientation of correlations in all varieties/traits considered (CC, Cabernet cortis; SN, Sauvignon nepis; SR, Sauvignon rytos; CV, Cabernet volos; SB, Sauvignon blanc; CS, Cabernet sauvignon; T/S, tolerance/susceptibility; GC, grape color). Green arrows indicate a positive correlation between the specific module and the trait/genotype. Red arrows indicate a negative association between the specific module and the trait/genotype considered. (B) Gene Set Enrichment Analyses of the tan and blue modules showing the top 10 enriched categories based on fold change. The threshold p-value was set to 0.01 (C) Heatmap visualizing gene expression within the blue and tan modules across all biological replicates of the six considered varieties, normalized using Z-scores.
Figure 3. Modules contemporaneously associated with both tolerance/susceptibility and one or more grapevine varieties analyzed. (A) Table showing the orientation of correlations in all varieties/traits considered (CC, Cabernet cortis; SN, Sauvignon nepis; SR, Sauvignon rytos; CV, Cabernet volos; SB, Sauvignon blanc; CS, Cabernet sauvignon; T/S, tolerance/susceptibility; GC, grape color). Green arrows indicate a positive correlation between the specific module and the trait/genotype. Red arrows indicate a negative association between the specific module and the trait/genotype considered. (B) Gene Set Enrichment Analyses of the tan and blue modules showing the top 10 enriched categories based on fold change. The threshold p-value was set to 0.01 (C) Heatmap visualizing gene expression within the blue and tan modules across all biological replicates of the six considered varieties, normalized using Z-scores.
Plants 13 02574 g003
Figure 4. Vesicle transport pathways in plants. COP-II vesicles mediate cargo transport from the ER to the cis-Golgi, while COP-I traffics the cargo from the Golgi to the ER and intra-Golgi as well. Clathrin-mediated endocytosis (CME) is the primary mechanism by which eukaryotic cells internalize extracellular or membrane-bound cargoes and it plays crucial roles in plant–microbe interactions Clathrin-coated vesicles (CCVs) are involved in the flow of cargo from the plasma membrane and trans-Golgi network to endosomes and retromers. Grapevine genes found to be enriched in the tan module are indicated in proximity to the related transport pathway.
Figure 4. Vesicle transport pathways in plants. COP-II vesicles mediate cargo transport from the ER to the cis-Golgi, while COP-I traffics the cargo from the Golgi to the ER and intra-Golgi as well. Clathrin-mediated endocytosis (CME) is the primary mechanism by which eukaryotic cells internalize extracellular or membrane-bound cargoes and it plays crucial roles in plant–microbe interactions Clathrin-coated vesicles (CCVs) are involved in the flow of cargo from the plasma membrane and trans-Golgi network to endosomes and retromers. Grapevine genes found to be enriched in the tan module are indicated in proximity to the related transport pathway.
Plants 13 02574 g004
Figure 5. Identification of absolutely specific genes in different grapevine varieties. (A) Distribution of the variety-specificity tau parameter over the 23,847 genes considered. (B) Bar graph showing the distribution of absolutely specific genes (ASG; tau = 1) and highly specific genes (HSG; tau > 0.85) over the six varieties considered. (C) Heatmap illustrating the expression of ASG in all biological replicates of the six varieties considered (Z-score normalized). (D) Scatterplot illustrating the relation/negative correlation r = −0.78) between specificity (tau) and expression in Sauvignon nepis. Blue dots represent all genes considered in the analysis, orange dots represent ASG in S. nepis, and red dots indicate the top optimal genes for S. nepis based on the score value. (E) Heatmap showing the expression of the top 10 optimal genes identified over the six varieties considered.
Figure 5. Identification of absolutely specific genes in different grapevine varieties. (A) Distribution of the variety-specificity tau parameter over the 23,847 genes considered. (B) Bar graph showing the distribution of absolutely specific genes (ASG; tau = 1) and highly specific genes (HSG; tau > 0.85) over the six varieties considered. (C) Heatmap illustrating the expression of ASG in all biological replicates of the six varieties considered (Z-score normalized). (D) Scatterplot illustrating the relation/negative correlation r = −0.78) between specificity (tau) and expression in Sauvignon nepis. Blue dots represent all genes considered in the analysis, orange dots represent ASG in S. nepis, and red dots indicate the top optimal genes for S. nepis based on the score value. (E) Heatmap showing the expression of the top 10 optimal genes identified over the six varieties considered.
Plants 13 02574 g005
Table 1. List of the ten best ranking genes based on a score value (0–2) corresponding to the sum of the quantile normalized expression of a given gene and its tau expression factor.
Table 1. List of the ten best ranking genes based on a score value (0–2) corresponding to the sum of the quantile normalized expression of a given gene and its tau expression factor.
Gene ID V3Gene ID V2TauScoreMean ExpFunctional Annotation
Cabernet sauvignon
Vitvi01g02070VIT_01s0127g009100.951.26556.55AERO1
Vitvi10g00903-1.001.24923.71Leucin-rich repeat protein kinase
Vitvi01g02281VIT_01s0010g040101.001.1829.77Unknown protein
Vitvi10g02153-1.001.1335.16-
Vitvi07g00496VIT_07s0005g023100.931.13312.47EXPA17
Vitvi05g02072-0.931.13011.93-
Vitvi10g00183-1.001.1254.63-
Vitvi10g02415VIT_00s2472g000100.851.11227.96Enhancer of mR-decapping protein 4
Vitvi01g02068-1.001.1113.84-
Vitvi10g02416-1.001.0983.25-
Cabernet cortis
Vitvi07g02026VIT_07s0130g002001.001.423238.36VvGELP21-Lipase GDSL
Vitvi19g00082VIT_19s0014g010601.001.417220.86Sesquiterpene synthase
Vitvi09g01530VIT_09s0002g019800.961.325113.34Myosin-like protein XIK
Vitvi11g01266VIT_11s0052g012300.931.322153.64Xyloglucan endotransglucosylase/hydrolase 23
Vitvi09g01648-0.961.30991.87-
Vitvi08g02288VIT_08s0007g045800.861.303319.66UGT73C2 (UDP-glucosyl transferase 73C2)
Vitvi19g00324VIT_19s0014g040001.001.28540.05Curculin (mannose-binding) lectin
Vitvi12g02451VIT_12s0134g006501.001.27434.63Anthocyanin 5-aromatic acyltransferase
Vitvi19g01982VIT_19s0014g051400.951.27164.65-
Vitvi15g00285VIT_15s0045g002701.001.26530.67Serine/threonine-protein phosphatase BSL3
Cabernet volos
Vitvi14g00668VIT_14s0036g009900.861.4321814.13Polyubiquitin (UBQ4)
Vitvi11g01637VIT_11s0052g002700.851.346652.97R protein MLA10
Vitvi11g00879VIT_11s0065g000400.891.331321.15CYP706A12
Vitvi04g00345VIT_04s0008g040000.871.295256.39Unknown
Vitvi03g01478VIT_03s0038g042300.881.240103.63Dihydroflavonol 4-reductase
Vitvi08g02374VIT_08s0007g077600.931.17521.69Polygalacturonase PG1
Vitvi16g01677-1.001.1698.42-
Vitvi08g00789VIT_08s0058g006501.001.1617.72Aldose reductase
Vitvi11g01568VIT_11s0065g007401.001.1426.04A -phase-promoting complex subunit 8
Vitvi01g01642VIT_01s0010g035501.001.1405.97Nuclear transcription factor Y sub-B related
Sauvignon blanc
Vitvi04g00029VIT_04s0008g003700.851.15657.39Clavata1 receptor kinase (CLV1)
Vitvi06g01648VIT_06s0004g025500.851.15052.71Kiwellin Ripening-related protein grip22
Vitvi04g00021VIT_04s0008g003000.851.14750.92Clavata1 receptor kinase (CLV1)
Vitvi09g01948-0.851.14045.92HcrVf2 protein
Vitvi03g00460VIT_03s0063g010000.851.12739.01Blue (type 1) copper domain
Vitvi07g01769VIT_07s0031g008500.871.09420.80Patatin
Vitvi01g01852VIT_01s0011g009900.871.08017.09RPM1
Vitvi10g00005VIT_10s0116g001500.871.07816.85Receptor kinase RK20-1
Vitvi00g02077VIT_00s0895g000100.871.06113.52Glucan 1,3-beta-glucosidase
Vitvi16g02124VIT_00s0294g001001.001.0602.29BR insensitive 1 receptor kinase 1
Sauvignon nepis
Vitvi09g01181VIT_09s0018g007801.001.506710.79HcrVf1 protein
Vitvi18g02399VIT_18s0089g010400.971.474647.26Avr9 elicitor response
Vitvi10g01863VIT_10s0003g035300.971.427382.88Lupeol synthase
Vitvi10g01875VIT_10s0003g036500.971.427381.69Beta-amyrin synthase
Vitvi13g02352VIT_13s0139g001900.921.408529.73Disease resistance protein RGA2
Vitvi12g02393VIT_12s0059g017900.911.356321.51Caffeic acid O-methyltransferase
Vitvi02g00721VIT_02s0012g016100.871.355506.49Beta-1,3-gluca -se precursor
Vitvi03g01757-1.001.31463.03-
Vitvi03g00910VIT_03s0167g000500.921.298140.99Conca-valin A lectin
Vitvi16g00665VIT_16s0022g004200.961.29481.21SRG1 oxidoreductase
Sauvignon rytos
Vitvi01g01410-1.001.32667.07-
Vitvi08g00957VIT_08s0040g009200.871.280207.61Glutathione S-transferase 25 GSTU7
Vitvi14g00080VIT_14s0060g009901.001.26931.94Unknown
Vitvi13g02566VIT_13s0156g003901.001.26731.10Myb family
Vitvi18g03265VIT_18s0089g010000.951.25651.80F-box family protein
Vitvi15g01230-1.001.22919.32-
Vitvi15g01425VIT_15s0021g014501.001.22418.06No hit
Vitvi13g01636VIT_13s0158g000501.001.22317.85Serine carboxypeptidase
Vitvi10g01830VIT_10s0003g024201.001.21716.27SRG1oxidoreductase
Vitvi17g00462-1.001.21616.08-
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Scariolo, F.; Gabelli, G.; Magon, G.; Palumbo, F.; Pirrello, C.; Farinati, S.; Curioni, A.; Devillars, A.; Lucchin, M.; Barcaccia, G.; et al. The Transcriptional Landscape of Berry Skin in Red and White PIWI (“Pilzwiderstandsfähig”) Grapevines Possessing QTLs for Partial Resistance to Downy and Powdery Mildews. Plants 2024, 13, 2574. https://doi.org/10.3390/plants13182574

AMA Style

Scariolo F, Gabelli G, Magon G, Palumbo F, Pirrello C, Farinati S, Curioni A, Devillars A, Lucchin M, Barcaccia G, et al. The Transcriptional Landscape of Berry Skin in Red and White PIWI (“Pilzwiderstandsfähig”) Grapevines Possessing QTLs for Partial Resistance to Downy and Powdery Mildews. Plants. 2024; 13(18):2574. https://doi.org/10.3390/plants13182574

Chicago/Turabian Style

Scariolo, Francesco, Giovanni Gabelli, Gabriele Magon, Fabio Palumbo, Carlotta Pirrello, Silvia Farinati, Andrea Curioni, Aurélien Devillars, Margherita Lucchin, Gianni Barcaccia, and et al. 2024. "The Transcriptional Landscape of Berry Skin in Red and White PIWI (“Pilzwiderstandsfähig”) Grapevines Possessing QTLs for Partial Resistance to Downy and Powdery Mildews" Plants 13, no. 18: 2574. https://doi.org/10.3390/plants13182574

APA Style

Scariolo, F., Gabelli, G., Magon, G., Palumbo, F., Pirrello, C., Farinati, S., Curioni, A., Devillars, A., Lucchin, M., Barcaccia, G., & Vannozzi, A. (2024). The Transcriptional Landscape of Berry Skin in Red and White PIWI (“Pilzwiderstandsfähig”) Grapevines Possessing QTLs for Partial Resistance to Downy and Powdery Mildews. Plants, 13(18), 2574. https://doi.org/10.3390/plants13182574

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop