Skip to main content

Identification of defense related gene families and their response against powdery and downy mildew infections in Vitis vinifera

Abstract

Background

Grapevine (Vitis vinifera) productivity has been severely affected by various bacterial, viral and fungal diseases worldwide. When a plant is infected with the pathogen, various defense mechanisms are subsequently activated in plants at various molecular levels. Thus, for substantiating the disease control in an eco-friendly way, it is essential to understand the molecular mechanisms governing pathogen resistance in grapes.

Results

In our study, we performed genome-wide identification of various defensive genes expressed during powdery mildew (PM) and downy mildew (DM) infections in grapevine. Consequently, we identified 6, 21, 2, 5, 3 and 48 genes of Enhanced Disease Susceptibility 1 (EDS1), Non-Race-specific Disease Resistance (NDR1), Phytoalexin deficient 4 (PAD4), Nonexpressor of PR Gene (NPR), Required for Mla-specified resistance (RAR) and Pathogenesis Related (PR), respectively, in the grapevine genome. The phylogenetic study revealed that V. vinifera defensive genes are evolutionarily related to Arabidopsis thaliana. Differential expression analysis resulted in identification of 2, 4, 7, 2, 4, 1 and 7 differentially expressed Nucleotide-binding leucine rich repeat receptor (NLR), EDS1, NDR1, PAD4, NPR, RAR1 and PR respectively against PM infections and 28, 2, 5, 4, 1 and 19 differentially expressed NLR, EDS1, NDR1, NPR, RAR1 and PR respectively against DM infections in V. vinifera. The co-expression study showed the occurrence of closely correlated defensive genes that were expressed during PM and DM stress conditions.

Conclusion

The PM and DM responsive defensive genes found in this study can be characterized in future for impelling studies relaying fungal and oomycete resistance in plants, and the functionally validated genes would then be available for conducting in-planta transgenic gene expression studies for grapes.

Background

Grapevine is one of the most commercially valuable fruit crops grown around the world. Nevertheless, its growth is significantly impacted by climate in the mass-producing areas, which allows for the development of various pathogenic diseases. Powdery mildew (PM) and downy mildew (DM) are two of those destructive diseases that propagate in the areas with colder temperatures, and higher relative humidity infecting grapevine [1]. The causal agent of PM is Erysiphe necator, an ascomycete fungus, while DM is caused by an oomycete, Plasmopara viticola. The epidemic of these diseases acutely diminishes both the productivity as well as quality by infecting all the green tissues in grapevine. As majority of the cultivated varieties belong to Vitis vinifera, which is susceptible to both PM and DM, it is vital to understand the molecular processes involved in grapevine resistance to powdery and downy mildew infection in order to manage these diseases in an eco-friendly manner. Under abiotic and biotic stress environments, different defensive mechanisms are triggered in plants as part of their innate immune response to protect themselves. Essentially, three main steps are involved in host defense against pathogenic attack, i.e. pathogen identification, signal transduction and gene induction; contributing to the production of molecules with antimicrobial properties. Once the emergence of PM and DM infections establishes in grapes, under optimal environmental conditions, the conidiospores of Erysiphe necator and Plasmopara viticola germinate on the grapevine leaf surface. This resulted in initial association of pathogen-associated molecular patterns (PAMPs) and grapevine defense molecules known as pattern recognition receptors (PRRs), thereby activating the first defense layer i.e. Pathogen-triggered immunity (PTI) [2]. However, the activated defense is not strong enough to deter the pathogenic augmentation in plants. In addition, due to evolution, certain pathogenic isolates have adapted to producing effector proteins in response, that are delivered within the cell, hence the PTI gets suppressed [3, 4]. During that time, another defense mechanism known as effector-triggered immunity (ETI) is activated in plants that triggers the expression of intracellular resistance (R) proteins often leading to apoptosis of the infected cell, thus, blocking any further pathogenic proliferation [5,6,7,8,9].

During plant immune response, numerous signal transduction pathways are triggered which lead to the host plant transcriptome re-programming and activation of the defense-responsive genes. Activated nucleotide binding sites leucine rich repeats (NLRs) trigger a cascade of intracellular immune responses, including an increase in reactive oxygen species (ROS) production, mitogen-activated protein kinase (MAPK) cascades, calcium spikes, phytohormone production, transcriptional reprogramming etc. [10]. Various plant hormones function as key players in activating the plant immune signaling network downstream the PTI or ETI induction [11,12,13,14]. Most often, resistance to biotrophic and hemibiotrophic pathogens is driven by salicylic acid (SA) signaling, while resistance to necrotrophic pathogens is mediated by jasmonate (JA) and ethylene (ET) signaling [15]. Since E. necator and P. viticola are both biotrophic pathogens, therefore, the defense mechanism triggered in the plant in response to infection with powdery and downy mildew might be modulated by SA signaling pathway. Once the SA signaling is initiated at the infection site, a similar response is triggered in undamaged plant tissues to provide protection against further pathogen invasion. This broad-spectrum long lasting resistance response is known as systemic acquired resistance (SAR) [16].

In the present research, we have focused on R-gene mediated defense signaling through SA accumulation. Once PAMPs or pathogen effectors are recognized in PTI and ETI, SA production is activated [17]. The R-gene mediated resistance has been previously documented to be triggered by 2 main proteins i.e. Enhanced Disease Susceptibility 1 (EDS1) and Non-Race-specific Disease Resistance (NDR1), both interacting with different types of resistance (R) proteins. EDS1 interacts with Toll/interleukin-1 receptor (TNL) type R proteins and NDR1 interacts with Coiled coil (CNL) type R proteins [18,19,20]. Recent study on Arabidopsis have demonstrated that any perturbation in the interactions between EDS1, NDR1 and R proteins can affect the activation of SA signaling pathway, thereby indicating the importance of these two main plant defense regulators [21,22,23,24]. Further, EDS1 forms heterodimers with Phytoalexin deficient 4 (PAD4) and Senescence-associated gene 101(SAG101), two proteins with which it shares sequence homology [25, 26]. These three proteins’ N-terminal domains allowed EP domains at C-terminal to bind and form interaction surfaces on the heterodimer [27].

Several other major molecules of SA mediated defense signaling includes avrPphB susceptible 2 (PBS2), required for Mla-specified resistance (RAR1), Nonexpressor of PR Gene (comprising of NPR1, NPR2, NPR3 and NPR4). RAR1, a chaperone protein is an important member of NLR mediated ETI utilized by both TIR-NLR and CC-type NLR [28]. The NPR class is a pivotal SA-mediated signal transduction regulator that is triggered after accumulation and subsequently binding of SA [29]. NPR1, NPR2, NPR3, and NPR4 are among the members of this protein family. According to the literature survey, NPR1 and NPR2 positively regulate SAR, while NPR3 and NPR4 regulate SAR negatively [30, 31]. The SA binds to NPR’s ankyrin-repeat domain, liberating it from BTB/POZ domain’s auto-inhibition. The C-terminal domain then associates with transcription factors (TFs) of WRKY and TGA class in nucleus, and further it binds with promoter sequences of Pathogenesis Related (PR) genes [32,33,34]. Recently, structural analyses of SA detection by several NPR proteins have been published, which aids in better understanding of mechanism of SA recognition by NPR proteins [35]. The ultimate players of plant defense mediated through SA are PR proteins that were described firstly in tobacco leaves infected with the tobacco mosaic virus [36]. Since then, these proteins have been identified in many monocot and dicot plants. Although these genes are wide-ranging, many of them are encoding for antimicrobial proteins [37]. PR-1 is a well-studied PR gene that is often used as a strong indicator of SA-mediated gene expression.

The importance of the role of SA or its derivatives in triggering defense-responsive proteins is well stated in literature. For instance, the plant susceptibility to pathogen increases when the plant becomes deficient in SA synthesis or accumulation [38]. It is evident in literature through experiments performed in plants like Arabidopsis that eds5 and eds16 mutants that are lacking SA accumulation are compromised in some R-gene pathways as well as SAR and basal resistance [38]. Various PR proteins that are known to induce SAR, get expressed when the plant is treated with SA or its derivatives [39] . This suggested the importance of SA in providing immunity to plants whenever the plant gets pathogenic infection. Previously, we identified 386 NLR genes in the grape genome, of which 63 such genes were found to be responsive to PM stress [40]. In the present study, we have carried out the genome-wide identification of various defense-responsive gene families (EDS1, NDR1, PAD4, NPR, RAR1 and PR). The evolutionary relatedness of identified gene families is studied with model plant Arabidopsis thaliana. The differential expression studies are conducted based on the PM and DM-responsive RNA seq data. As a result, we have identified various PM and DM-responsive defensive genes of above mentioned gene families (EDS1, NDR1, PAD4, NPR, RAR1 and PR) along with their co-expression to identify highly correlative defensive genes responsive to both PM and DM infections. The different genes obtained through co-expression analysis are functionally characterized by the Blast2GO results. Additionally, various cis-acting regulatory elements (CARE) for the identified genes have also been predicted against the PM and DM stresses.

Results

Genome-wide identification of various defense-responsive gene families in Vitis vinifera

The genome-wide identification of various classes of defensive genes in V. vinifera, namely EDS1, NDR1, PAD4, NPR, RAR1 and PR in was conducted by running standalone BLASTP between protein coding sequences of V. vinifera and defensive gene families. For analysis, we collected 38,120 protein coding sequences of V. vinifera from the NCBI genome database. Next, we retrieved 41, 127, 17, 34, 9 and 92 coding sequences of EDS1, NDR1, PAD4, NPR, RAR1 and PR proteins respectively from A. thaliana, A. lyrata, O. sativa, P. trichocarpa and V. vinifera using NCBI. The sequence alignments were done and consequently, 27, 47, 29, 928, 3 and 148 putative gene hits from EDS1, NDR1, PAD4, NPR, RAR1 and PR gene families were detected in V. vinifera. The confirmation of the presence of respective domains was done through BLASTP analysis of putative hits against the non-redundant (NR) database. Finally, 6, 21, 2, 5, 3 and 48 such sequences were obtained that were considered as EDS1, NDR1, PAD4, NPR, RAR and PR proteins after validation through bioinformatics tools and databases such as Genoscope, CDD, Pfam, InterProScan and SMART. The information regarding protein IDs and transcript IDs of all gene families is given in Additional file 1.

The identified genes were characterized in silico on the basis of gene structure, motif presence, chromosomal locations and physicochemical properties. The gene structure analysis was performed by using SPLIGN tool which provided valuable information about gene arrangement in terms of exons and introns. The results depicted that the number of introns in the EDS1, NDR1, PAD4, NPR, RAR1 and PR gene families ranged from 1 to 2, 1–2, 1–3, 1–4, 1–4, 1–4 and 1–2 respectively (Additional file 2). Further, on the basis of the localization of an intron with respect to reading frame of gene, three intron phases were depicted. All EDS1 gene sequences were present in phase 0; NDR1 in phase 0 and phase 1; PAD4 in phase 1 and phase 2; NPR, RAR1 and PR gene sequences in all the three intron phases respectively (Additional file 2). Next, the protein sequences were assessed with ‘MEME’ server to identify conserved motifs. A total of 16 conserved motifs were identified in the studied gene families (Additional file 2). Their final verification for the presence of desired motif was done through Pfam server and InterProScan. Conserved domains identified were Lipase 3 and EDS1 EP in EDS1; LEA2 domain in NDR1; EDS1 EP in PAD4; Ank2, Ank5, NPR1_like_C; NPR1_interact in NPR; CHORD domain in RAR1; CAP, Thaumatin, Glycohydro 17, Glycohydro 18, Glycohydro 19, Betv1, Barwin and Chitin_bind_1 in different classes of PR coding sequences. The block diagrams and the logos of the conserved motifs are shown in additional file 2.The identified genes were found to be distributed widely in the grape genome. All EDS1, PAD4 and RAR1 genes were located on single chromosome i.e. 17, 7 and 16; NPR genes were located on chromosomes 7, 8, 10 and 11; NDR1 and PR genes were unevenly distributed at several chromosomes. However, no gene sequences were identified on chromosome number 1 and 12 in any gene family (Additional file 2). The various parameters computed to determine physicochemical properties included amino acid length, molecular weight, isoelectric point (pI) and instability index (Additional file 1). The range of amino acid length for the studied proteins is 112–629 and the molecular weight range is between 12 and 71.3 kDa. Based on pI prediction, it was depicted that EDS1 proteins were acidic, basic and neutral in nature; all the NDR1 proteins were basic in nature; PAD4 and NPR proteins were acidic; RAR1 proteins were neutral and basic in nature; PR were mostly acidic and some are basic in nature. The results of instability index prediction suggested that most of the proteins identified were unstable under in-vitro conditions with the exception of some NDR1 and PR proteins of PR1, PR2, PR10.1, PR10.7 families that were found to be present in stable state.

Phylogenetic analysis

The evolutionary history of EDS1, NDR1, PAD4, NPR, RAR1 and PR genes was traced by conducting their phylogenetic analysis with model plant A. thaliana. Full length amino acid sequences of all genes families were used to perform multiple sequence alignment and subsequently, phylogenetic trees were constructed which were classified into different phylogenetic clades (Fig. 1). The analysis depicted that all the studied genes families of V. vinifera and A. thaliana were closely related. The phylogenetic tree of EDS1 proteins was classified into six clades, with clades II-VI comprising V. vinifera proteins exclusively and clade I carrying A. thaliana proteins. The NDR1 proteins from both plants were divided into 13 clades, with clades I, II, X and XI including primarily A. thaliana proteins and the remaining clades containing V. vinifera proteins. The PAD4, NPR, and RAR1 proteins were classified into five, four, and four clades, respectively, and the majority of the protein members of both plants were found in separate clades. Likewise, we constructed distinct phylogenetic trees for three different PR protein classes, namely PR1, PR2, and PR5. The results showed that PR1 proteins from both plants were classified into 18 clades, PR2 proteins were classified into four clades and PR3 proteins were grouped into five clades. Overall, the analysis showed that majority of the members of studied gene families from both plants were evolutionary conserved, hence clustered together and placed in distinct clades (Fig. 1).

Fig. 1
figure 1

Phylogenetic trees of defense responsive proteins of V. vinifera and A. thaliana A.) EDS1, B.) NDR1, C.) PAD4, D.) NPR, E.) RAR1, F.) PR1, G.) PR2, H.) PR5. Visualization of phylogenetic trees was done through iTOL. The colour of the labels indicate different clades, while those of leaf nodes depict the species

Identification of PM and DM-responsive defensive genes in V. vinifera

To better understand the R-gene mediated molecular mechanism of powdery and downy mildew resistance in grapevine, the identification of PM and DM-responsive defensive gene families (NLR, EDS1, NDR1, PAD4, NPR, RAR1 and PAD4) and transcription factors (TF) was done by conducting the differential gene expression analysis. We reported the occurrence of 386 NLR gene sequences in the grape genome in our previous research, including 63 such sequences that were sensitive to PM stress [40]. In this study, we found 2 new highly expressing PM-responsive and 28 DM-responsive NLR genes. Additionally, 4 4, 7, 2, 4, 1 and 7 EDS1, NDR1, PAD4, NPR, RAR1 and PR genes expressing differentially in reaction to PM infection stress and 2, 5, 4, 1 and 19 differentially expressed DM-responsive defensive genes (EDS1, NDR1, NPR, RAR1 and PR) were identified (Fig. 2). The differentially expressed genes were represented as heat maps that depicted the expression fold change. During analysis, it was observed that for the triplicates of each condition; the TPM values obtained were nearly same. Therefore, we took the cumulative mean of these values for every triplicate during heat map representation.

Fig. 2
figure 2

Heat maps illustrating the differential expression patterns for defensive genes in Vitis vinfera. (a) PM-responsive defensive genes at time points of 0 h and 36 hpi. The gradient value of color mapping ranges from − 1 to + 1 (b) DM-responsive defensive genes at time points of 0 h, 1 dpi, 2 dpi and 3 dpi. The gradient value ranges from − 1.73 to + 1.73. Green color indicates declined expression level, red color indicates elevated expression levels, black color indicates no expression. PM and DM defensive genes are clustered in two groups indicated by gene names with blue and red color

The TPM findings of NLR genes showed that most of the NLR genes were up-regulated in later stages of pathogen infection relative to initial phase during both PM and DM-infection (Fig. 2). Analysis of EDS1 and NDR1 gene families showed that all EDS1 and 5 NDR1 genes were up-regulated during PM-infection relative to control condition (Fig. 2A). Likewise, during DM-infection, all EDS1 and 4 NDR1 genes were up-regulated at 3 dpi (Fig. 2B). The analysis of PAD4 gene family revealed that 2 PAD4 genes were detected in PM infected samples, one of which was up-regulated (XP_010652916.1) while the other one was down regulated (XP_002275637.1) during PM infection as compared to control conditions (Fig. 2A). No differentially expressed PAD4 genes were identified in DM-infected sample. We found 1 PM-responsive and 1 DM-responsive differentially expressed transcripts encoding RAR1 in both PM and DM infection and both transcripts demonstrated up-regulation during pathogenic stress (Fig. 2). The DEA of NPR proteins indicated that 4 PM-responsive and 4 DM-responsive NPR (comprised of NPR1 and NPR2) genes were expressed differentially during infection and all NPR genes were shown to be steadily up-regulated from initial to later stages of both powdery and downy mildew infection (Fig. 2).

Further, the DEA of 3 transcription factors (TFs) families i.e. TGA, WRKY and NAC was done as T.Fs are known to play a significant role in stimulating the PR genes expression during R-gene mediated defense signaling after interacting with NPR genes [32,33,34]. By using Plant Transcription Factor Database (TFDB), the total number of TGA, WRKY and NAC TFs retrieved was 47, 59 and 71. Their differential expression was conducted and consequently, 20 PM-responsive and 39 DM-responsive TFs were identified. Following that, we found 7 PM-responsive and 19 DM-responsive differentially expressed PR genes, with the majority of genes showing a significant surge in expression levels during infections.

Overall, we identified 27 PM-responsive and 59 DM-responsive differentially expressed defensive genes (NLR, EDS1, NDR1, PAD4, NPR, RAR1 and PAD4) in V. vinifera accessions used in the current study (Fig. 2). The findings indicated that most genes were up-regulated during PM and DM infection relative to control conditions.

Co-expression analysis of PM and DM responsive defensive genes in V. vinifera

As co-expression analysis effectively detects correlation between genes and provides an important information about which genes are expressing simultaneously, we conducted co-expression (CE) analysis of differentially expressed PM and DM-responsive defensive genes using the software tool CoExpress 1.5.2. Consequently, we screened some highly co-expressing genes in V. vinifera that could regulate the powdery and downy mildew infection through R-gene mediated signaling cascade (Fig. 3). During PM stress, 6 NLRs, 1 EDS1, 1 NDR1, 1 PAD4, 3 NPRs (comprised of 2 NPR1 and 1 NPR2), 5 TFs and 7 PRs co-expressed with 41, 15, 14, 15, 18, 36 and 73 other transcripts respectively, of defensive gene families as studied above. However, no RAR1 was found to be co-expressed with any transcript of defensive gene family during PM stress (Fig. 3A). Likewise, when differentially expressed DM-responsive genes were analyzed for co-expression, 7 NLRs, 2 EDS1, 2 NDR1, 1 RAR1, 1 NPR2, 5 TFs and 4 PRs co-expressed with 43, 19, 18, 13, 13, 38 and 44 other transcripts respectively (Fig. 3B). The exception in this case is PAD4 that was not found to be co-expressed with any defensive gene family transcript. Altogether, we screened 24 PM-responsive and 22 DM-responsive differentially expressed genes and TFs in V. vinifera accessions, which were found to be highly co-expressing amongst each other (Fig. 3).

Fig. 3
figure 3

Co-expression analysis of defensive genes showing pairwise correlation of different protein families (NLR, EDS1, PAD4, NPR1, WRKY, TGA, PR) using  Gephi 0.91 software. (a) Co-expression amongst PM-responsive defensive genes visualized by cytoscape 3.7.1. (b) Co-expression amongst DM-responsive defensive genes was shown

qPCR expression analysis

On the basis of TPM and co-expression analysis, we picked 14 PM-responsive and 13 DM-responsive highly expressing defensive genes for real time expression analysis under various stress conditions (Fig. 4 and Fig. 5). The primer sequences, amplicon length and annealing temperature of PM and DM-responsive defensive genes are listed in Additional file 3.

Fig. 4
figure 4

Relative expression levels of PM-responsive defensive genes during PM infection and salicylic acid (SA) treatment. The fold change was calculated relative to healthy leaf samples which were assigned a value of 1. Error bars indicate the standard deviation (SD) of means. Significant differences at P < 0.001, P < 0.01, P < 0.1 have been shown with 3 (***), 2 (*) and 1 (*) asterisks

Fig. 5
figure 5

Relative expression levels of DM-responsive defensive genes during DM infection and salicylic acid (SA) treatment. The fold change was calculated relative to healthy leaf samples which were assigned a value of 1. Error bars indicate the standard deviation (SD) of means. Significant differences at P < 0.001, P < 0.01, P < 0.1 have been shown with 3 (***), 2 (*) and 1 (*) asterisks

For analysis, 3 different types of stresses were chosen i.e. PM stress, DM stress and salicylic acid (SA). In previous studies, it was seen that the exogenous application of SA or SA analogs induces SAR and restores resistance in numerous mutants compromised in signaling steps upstream of SA production [25, 41]. In the present study, as we are focusing on genes of R-gene mediated defense signaling which is known for having the involvement of SA signaling, therefore we studied its effects on resistance response against PM and DM infections. On the leaves of Thompson seedless variety, the SA was applied exogenously at 5 different time intervals i.e. 1 hpi, 2 hpi, 1 dpi, 2 dpi and 3 dpi to test the expression pattern changes of various defense-responsive genes. In addition, we also collected PM and DM infected leaf samples at various time intervals of 1 dpi for PM and 1 dpi and 3 dpi for DM infection and performed qPCR. The results of qPCR indicated that mostly the expression levels of PM-responsive defensive genes were higher during 2 dpi PM infection, 2 hpi and 3 dpi SA treatment (Fig. 4). At 1 hpi of SA treatment, the higher relative expression of only NDR1 gene was observed compared to other gene families; after 1 hpi and upto 2 hpi of SA treatment, expression levels of all genes increased; however at 1 dpi SA treatment, all genes were down regulated compared to 2 hpi; then slight increase in expression level at 2 dpi; and at 3 dpi of SA treatment, maximum up-regulation was observed and all genes were up-regulated compared to 2 dpi. Overall, it was found that PR genes were showing the maximum expression at 2 hpi; all other gene classes (NLR, EDS1, NDR1, NPR2) were showing maximum expression at 3 dpi. Additionally, for PM-infected sample at 2 dpi, the relative expression levels for all genes were higher compared to control condition. Also, the expression levels for all gene families at 2 dpi PM-infection was comparable to that of SA treated samples at 2 hpi and 3 dpi (Fig. 4).

Likewise, the expression levels of DM-responsive defensive genes were higher mainly during late infection i.e. 3 dpi DM infection, 1 hpi and 3 dpi SA treatment (Fig. 5). At 1 hpi of SA treatment, the relative expression of all genes was higher; expression level decreased from 1 hpi to 2 hpi of SA treatment for all genes; however at 1 dpi of SA treatment, there was an increase in expression levels of some of the genes (NLR, NDR1, PR2); at 2 dpi, almost similar expression levels as that at 1 dpi was observed; however, at 3 dpi, all genes were up-regulated as compared to 2 dpi condition. At 1 dpi of DM infection, the relative expression levels of NLR, NDR1, EDS1, NPR2 and TGA TFs coding genes was higher as compared to control. However, none of the PR coding genes were showing an up-regulation. At 3 dpi of DM infection, the expression levels of most of the genes increases as compared to 1 dpi condition. Additionally, at 3 dpi of DM infection, the slight increase in expression level of PR1 gene was observed (Fig. 5). Overall, if we compare an expression levels between 1 dpi and 3 dpi of DM infection, expression levels of all genes increased from 1 dpi to 3 dpi DM infection except PR2 gene that was down regulated and PR1 gene that was up-regulated at 3 dpi DM infection. Moreover, if we compare 1 dpi and 3 dpi SA treatment condition and 1 dpi and 3 dpi DM infection; almost similar trend was followed from 1 dpi to 3 dpi for all genes.

Functional annotation of PM and DM-responsive defensive genes

In order to recognize the diverse roles of proteins at the molecular level, we performed Blast2GO of highly interacting defensive genes found in the co-expression study. From this analysis, the depiction of biological process, cellular component and molecular function, the 3 key characteristics of genes was done. Various GO terms were assigned to gene sequences for determining the possible functions of PM and DM-responsive defensive genes. For both PM and DM-responsive genes, in the biological process category, we obtained maximum gene sequences that were stimulus responsive (GO:0050896). [Additional file 4A(a) and Additional file 4B(a)]. Likewise, the GO term annotations in cellular component category suggested that most of our gene sequences performed their respective functions in cell (GO:0005623), cell part (GO:0044464) and organelle (GO:0043226) [Additional file 4A(b) and Additional file 4B(b)]. The findings of molecular process category revealed that maximum number of sequenceswere involved in binding activity (GO:0005488)) [Additional file 4A(c) and Additional file 4B(c)]. Altogether, our results depicted that the PM and DM-responsive defensive genes found in our study may play an immense role in biological, cellular and molecular processes.

Prediction of cis-acting regulatory elements in PM and DM-responsive defensive genes

Predictions of cis-regulatory elements in PM and DM-responsive defensive genes derived from a co-expression study revealed the existence of two key cis-acting elements i.e. TATA and CAAT boxes in the promoter and enhancer parts of all gene sequences. In addition, 5 classes of PM and DM-responsive regulatory elements were identified, namely defense and stress responsive element, salicylic acid responsive element, wound and pathogen responsive element, stress responsive element, As-1 and TGA element. Among all, the defense and stress responsive element was located in 12 PM-responsive and 10 DM-responsive promoter sequences, salicylic acid responsive element in 15 each PM-responsive and DM-responsive sequences, wound and pathogen responsive element in 17 PM-responsive and 16 DM-responsive sequences, As-1 in 10 each PM-responsive and DM-responsive sequences and TGA motif in 6 PM-responsive and 9 DM-responsive promoter sequences (Fig. 6). Taken together, the most frequent cis elements found in maximum number of PM and DM-responsive gene sequences were from wound and pathogen responsive element category.

Fig. 6
figure 6

Assessment of promoter sequences of stress responsive defensive genes. (a) PM-stress responsive (b) DM-stress responsive cis-elements were shown in different colors. Each color represents different cis-element

Discussion

As the cultivation of grapes is vulnerable to the emergence of diseases like PM and DM, efforts are increasingly directed towards the understanding of defense mechanisms of disease resistance at the molecular level. Various signaling pathways are triggered in plants at molecular level during pathogenic infection that provide plant defense. It has been reported in literature that during invasion by biotrophic organisms, the active defense response in plants is signaled by SA [32]. In the present study, we have identified various defense-responsive genes (EDS1, NDR1, PAD4, NPR, RAR1 and PR) involved in SA signaling pathway. The members of these gene families have also been identified previously in grapevine and other plants, but there is no report of genome-wide recognition study of such gene families. The phylogenetic analysis of V. vinifera defensive gene families is conducted with model plant A. thaliana to understand their evolutionary relatedness. The findings show that the orthologous genes from both plants share a close evolutionary relatedness with each other, implying that the EDS1, NDR1, PAD4, NPR, RAR and PR genes have been conserved in different plants. The evolutionary conservation also indicated that these gene families play an important role in plant defense. Additionally, the paralogous defensive genes are found to be closely grouped, which indicated their structural homology.

Next, the changes in the expression patterns of above mentioned defensive gene families under pathogenic and non-pathogenic conditions are studied by conducting differential expression. Our analysis depicted that during infection, most of the defensive gene classes were up-regulated relative to non-pathogenic conditions, indicating the active involvement of SA-mediated defense pathway during infection with PM and DM. To confirm the TPM findings, we performed real time PCR of PM and DM infected leaf samples to test differences in relative expression between healthy and infected samples. In addition, since SA plays a prominent role in activating the defense mechanism, we gave SA treatment to leaf samples at different intervals of time. The results of real time study revealed that during SA treatment, PM infection and DM infection, most of the defense-responsive genes were up-regulated relative to control condition. This meant that anytime the plant suffers an outbreak, the innate immune response of the plant is activated through which various signaling pathways and comes into action. The expression level of defense-responsive genes is enhanced during infection as compared to non-pathogenic condition that eventually induces SAR in the plant. We found that in case of PM-responsive genes, initially, the expression level of most genes increase from 1 hpi to 2 hpi of SA treatment; however, a sudden decrease in expression levels was witnessed after 2 hpi till 2 dpi. After that, ultimately, at 3 dpi, a tremendous increase in most of the genes was visualized. In case of DM-responsive genes, initially, the defense-responsive genes were highly expressed at 1 hpi of SA treatment; after that there was a decline in expression levels; and finally at 3 dpi, most of the defensive genes were up-regulated. The possible assumption would be that following an infection, SA levels will rise in plant, leading to defense reaction. Thus, when we gave SA treatment to our samples, initially up-regulation of genes occurs at 1 h or 2 h; and later on to provide SAR to the plant, the expression levels of both PM and DM-responsive defensive genes peaks up at 3 dpi. Additionally, our results of powdery and downy mildew infection through quantitative real time PCR are in accordance with the findings of the TPM study. Most of the genes were up-regulated during PM and DM infections as compared to control conditions in both TPM and real time studies.

Past studies have shown that over-expression of these genes confers broad spectrum resistance to infected plants. For instance, one analysis in Arabidopsis showed that over-expression of EDS1 enhances the accumulation of PR1 protein that ultimately provides broad spectrum resistance to the plant [41]. NPR proteins have been reported to engage in SA signal transduction in various plants. Once activated by SA, these genes regulate the expression of PR genes by interacting with various TFs. Through differential expression analysis, we have identified TFs of 3 different classes (bZIP, WRKY and NAC), that might activate PR proteins by interacting with NPR1 and NPR2 proteins. Previously, it has been shown in Arabidopsis that the PR-1 gene expression is regulated through NPR1 gene binding to TGA transcription factors [42]. Similarly WRKY and NAC transcription factors were also found to interact with NPR1 gene in previous studies [43, 44]. Additionally, over-expression of the NPR1 gene in plants like Arabidopsis, wheat, tomato and apple provides broad spectrum resistance against fungal and bacterial pathogens [45,46,47,48]. In one study, conducted in Fuji apples, it has been shown that over-expression of NPR1 genes enhances PM resistance through SA mediated signaling pathway [49]. All these studies supported our results of expression analysis of various defense responsive genes conducted in grapes. Next, we performed the co-expression analysis of different genes and constructed their co-expression networks. Such networks provide very important information about candidate disease gene identification, regulatory genes identification, gene annotation etc. [50]. Apart from conventional co-expression networks, differential co-expression analysis is becoming more popular now-a-days [51]. This method identify genes with different co-expression pairs under diverse conditions like disease states, developmental stages, tissue forms etc. [51,52,53,54,55,56]. In the present study, we performed differential co-expression between the differentially expressed transcripts obtained from the gene expression data. All gene classes were analyzed for co-expression by taking two gene classes at a time. First, a pairwise correlation was examined for each gene class pair. Next, highly co-expressing transcripts with a threshold value of 0.99 to 1 were screened and their co-expression networks were constructed.

Next, to allocate roles to PM and DM-responsive defensive genes identified in the present research, Blast2GO was done. The study of GO annotations has showed that the maximum defense-responsive gene sequences are engaged in biological processes, accompanied by sequences in cellular component and molecular function category. The GO terms in biological process category are directly implicated in defense response which supported the functioning of proteins identified in our study in resistance mechanism. The proteins in cellular component class are found to perform their role inside cell which clearly indicated that the defensive pathway we are targeting is functional inside the cell. In the category of molecular processes, many of the proteins are involved in binding activity which supported the interacting nature of defense-responsive proteins.

In our research, we have recognized multiple regulatory elements that react to pathogenic stress, namely defense and stress responsive element (TC-rich repeats), salicylic acid responsive element (TCA element), wound and pathogen responsive element (WUN-motif and W-box), stress responsive element (GT1 Box), As-1 and TGA; they provide an important information regarding upstream control of PM and DM-responsive defensive proteins. The function of the TFs in response to PM and DM-stress has been well elucidated in previous studies. For example; the NAC TFs associate with GT1 box, bZIP bind to TGA elements and WRKY TFs associates with W-box in reaction to fungal infection [57,58,59]. In the present analysis, we have observed the involvement of these elements in defensive proteins, which clearly confirmed that these genes were up-regulated as a defense mechanism in the situation of PM and DM stress. Altogether, numerous regulatory elements are identified in the present study that could play an important role in regulating the expression of stress-responsive genes.

In the present study, we focused on molecular aspects and gene expression analysis to identify various defense-responsive genes and study their response against PM and DM infections in grapes. In future, the identified genes can be characterized for resistance against PM and DM in Arabidopsis plants. Further, the functionally validated genes will be available for transgenic studies in grapes for imparting fungal resistance.

Conclusion

In our research, grapevine genome was examined to identify various classes of defensive genes that were directly implicated in R-gene mediated defense signaling during powdery and downy mildew infections. In total, we identified 6, 21, 2, 5, 3 and 48 EDS1, NDR1, PAD4, NPR, RAR and PR genes in the grape genome amongst which 27 PM-responsive and 59 DM-responsive differentially expressed defensive genes (NLR, EDS1, NDR1, PAD4, NPR, RAR and PR) were identified in 2 different V. vinifera accessions. The gene co-expression study depicted 24 PM-responsive and 22 DM-responsive defensive genes and Transcription factors (TF) that were found to be highly co-expressing with each other. Taken together, we can summarize that defensive genes identified in the current research could be beneficial in improving grapevine resistance to powdery and downy mildew.

Methods

Genome-wide identification of various defense-responsive gene families in Vitis vinifera

For the identification of V. vinifera EDS1, NDR1, PAD4, NPR, RAR1 and PR gene families, all possible protein coding sequences were procured from NCBI genome database (ftp://ftp.ncbi.nlm.nih.gov/genomes/Vitis_vinifera/protein/) and constructed their local sequence database. Next, the previously documented protein sequences of above listed gene families from different plants namely Arabidopsis thaliana, Arabidopis lyrata, Oryza sativa, Populus trichocarpa and Vitis vinifera were obtained from NCBI (http://www.ncbi.nlm.nih.gov/protein) and their FASTA sequences were compiled into a file which was considered as a query file. To identify the putative hits in the grape genome, standalone protein BLAST (Basic Local Alignment Search Tool) was performed between the database and query file created above and the e-value used was 1e-05, respectively. The candidate gene sequences were verified by assessing the occurrence of respective domains of various gene families at different servers such as Grape Genome Browser (12X) (http://www.genoscope.cns.fr/externe/GenomeBrowser/Vitis/), conserved domains database (CDD) (http://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi) [60], InterProScan (https://www.ebi.ac.uk/interpro/search/sequence-search) [61] and SMART (http://smart.embl-heidelberg.de/) [62]. Next, the ‘In silico characterization of defense-responsive gene families was done based on their gene structure, motif presence, chromosomal locations and physicochemical properties. For the gene structure prediction, an online tool ‘Splign’ was used to detect exons and introns (http://www.ncbi.nlm.nih.gov/sutils/splign/splign.cgi/). Visualization of exon and intron locations was done by using an online tool i.e. GSDS2.0 (http://gsds.cbi.pku.edu.cn/) [63]. The presence of conserved motifs was predicted through Multiple Expectation Maximization for Motif Elicitation (MEME) Suite (http://meme-suite.org/) and Pfam database (https://pfam.xfam.org/) [64, 65]. To determine the chromosomal locations, Grape Genome Browser at Genoscope (http://www.genoscope.cns.fr/externe/GenomeBrowser/Vitis/) and Ensembl Plants portal (https://plants.ensembl.org/Vitis_vinifera/Location/Genome/) were used. The physicochemical properties of defensive gene families were depicted by using ProtParam bioinformatics tool (https://web.expasy.org/protparam/) with default parameters [66].

Multiple sequence alignment and phylogenetic analysis

Multiple sequence alignment of protein sequences of EDS1, NDR1, PAD4, NPR, RAR1 and PR gene families was performed to identify conserved residues using Clustal W program at default parameters (http://www.ebi.ac.uk/Tools/msa/clustalw2/). To study the evolutionary relatedness of EDS1, NDR1, PAD4, NPR, RAR1 and PR proteins with A. thaliana, phylogenetic trees of all gene families were constructed following the Neighbor-joining method in MEGA 7 with 1000 bootstrap replicates. The JTT matrix-based method was used to compute the evolutionary distances [67,68,69,70]. For visualization of phylogenetic trees, an Interactive Tree of Life (iTOL) online software tool (https://itol.embl.de) was used that worked by taking Newick phylogenetic tree format as an input [71].

Identification of PM and DM responsive defensive gene families in V. vinifera

To identify PM and DM responsive defensive genes, we made a FASTA file of protein coding sequences of defensive gene families (EDS1, NDR1, NPR1, RAR1, PR) and NLR genes identified in our previous research [40]. In addition, as transcription factors (TFs) are also presumed to participate in R-gene mediated defense pathway, various transcription factors (TFs) of V. vinifera belonging to 3 distinct families i.e. TGA, WRKY and NAC were retrieved from Plant Transcription Factor Database (TFDB) [72]. The RNA-seq data of two V. vinifera varieties i.e. Thompson Seedless and Pinot Noir was procured from NCBI’s Sequence Read Archive (SRA) database and the project I.Ds for these studies were SRP116308 and PRJEB24540 (http://www.ncbi.nlm.nih.gov/sra/). The transcriptomic data of Thompson Seedless was responsive to PM and was accessible at 1 time point i.e. 36 h post inoculation (hpi) (Northwest A&F University). The RNA-seq data of Pinot Noir was DM-responsive and retrieved at 3 time points i.e. 1 dpi, 2 dpi and 3 dpi [73]. Each of the derived SRA data consisted of three biological replicates. The identification of PM and DM-responsive defensive genes is based on the digital expression analysis of different groups of proteins identified in our study. The Trinity-V2.4 package’s RSEM (RNA-Seq by Expectation-Maximization) program was used to quantify the abundance of all PM and DM-responsive defensive genes as transcript per kilobase million (TPM) [74]. All three biological replicates of every condition were analyzed individually. Thereafter, EdgeR was used to calculate the differential gene expression by assigning a cut-off value of 4-fold change and P-value of 0.001. The visualization of differential gene expressions was done by constructing heat maps with Hierarchical Clustering Explorer 3.5 (http://www.cs.umd.edu/hcil/hce/) [75].

Co-expression analysis of PM and DM responsive defensive genes in V. vinifera

Next, we conducted co-expression (CE) analysis of differentially expressed defensive genes to identify a group of genes expressing simultaneously during PM and DM-infection. The software used for detecting the co-expression was CoExpress 1.5.2. This is a stand-alone software tool that is focused on the Pearson correlation coefficient (R). The parameters were set to default while running this software. (http://Bioinformatics.lu/CoExpress/). The expression values of two sets of genes were loaded at a time and a linear correlation was measured amongst them. Finally, we screened genes that are highly co-expressing with the R value of 1 or close to 1. The visualization of co-expression networks for genes was achieved using the software Gephi 0.91 [76].

Plant material and treatment

As SA is playing a pivotal role in R-gene mediated defense signaling against PM and DM-infections, we studied its effect on Thompson Seedless variety of V. vinifera. The stem cuttings were collected from National Research Centre for Grapes (NRCG), Pune. The Centre has the National Active Germplasm Site for grapes, where the grape germplasm is being maintained. The institute code number for Thompson seedless variety of V. vinifera is A37–3. The collected stem cuttings were grown in pots containing soil:soil-rite in the ratio of 2:1 in growth chamber in Department of Biotechnology, Panjab University, Chandigarh. The experiment was performed at 5 different time intervals i.e. 1 hpi, 2 hpi, 1 dpi, 2 dpi and 3 dpi in the sets of three biological replicates at each time point. The SA treatment was given at the final concentration of 100 mg/litre. The control and treated leaves were collected, snap frozen and stored in − 80 °C till further use. Additionally, we also studied the effects of PM and DM infection on the activity of defense-responsive genes identified in the present study. We collected the healthy and infected leaf samples of Thompson Seedless variety of V. vinifera from National Research Centre for Grapes (NRCG), Pune. The time points for PM and DM infected leaf samples collected were 2 dpi for PM infected samples and 1 dpi and 3 dpi for DM infected samples. During sample collection, the symptoms were not evident on the leaf surface. The collection was done on the basis of SRA data. During visible condition, the PM-infected leaf is covered with a dusty whitish powdery growth and DM appears as yellowish circular spots encircled by a brownish-yellow halo that covers the majority of the leaf surface. This experiment was also performed in the sets of three biological replicates at the respective time points.

RNA isolation, cDNA synthesis and quantitative real time PCR analysis

Based on the genes identified in the co-expression analysis, we performed quantitative real time PCR (qPCR) of PM and DM-responsive differentially co-expressed defensive genes to measure the relative expression. Leaf samples were used to extract the total RNA by following Ghawana et al. 2011 protocol and the cDNA was prepared using the Superscript III first strand cDNA synthesis kit (Invitrogen USA) [77].

We used two housekeeping genes namely actin (ACT) and elongation factor 1 (EF1) as internal reference genes for the normalization of qPCR results. The real time primers were designed by using Primer 3 software (http://primer3.ut.ee/). The qPCR experiment was performed with three replicates of healthy, PM-infected, DM-infected and SA treated leaf samples by employing Bio-Rad CFX96 Real-Time PCR detection system. The following conditions were used while performing an experiment: 95 °C for 7 min, followed by 40 cycles of 95 °C for 20 s, Tm for 20 s and 72 °C for 20 s. In order to interpret the results of qPCR, we used REST 2009 algorithm (Qiagen) and 2(−∆∆CT) method to obtain relative gene expression ratios of target genes with respect to control genes (http://www.REST.de.com/) [78, 79]. All the triplicate experiments are represented as means ± standard deviations (SDs). Analysis of variance (ANOVA) followed by Tukey’s multiple comparisons test was conducted at p-value < 0.05 to detect the significant differences in gene expressions amongst control and treatment conditions using GraphPad Prism 9.2.0 software.

Functional annotation of PM and DM responsive defensive genes in V. vinifera

For assigning biological functions to PM and DM-responsive defensive gene sequences, Blast2GO (https://www.blast2go.com/blast2go-pro/) tool was used [80]. Essentially, Blast2GO annotation requires 3 steps: BLAST to locate homologous sequences, Mapping to allocate GO terms to each hit obtained, Annotation to assign role to query sequences on the basis of their cellular positions, molecular functions and biological processes. Furthermore, knowledge regarding various protein domains or motifs may also be obtained from Blast2GO by running InterProScan.

Promoter study of PM and DM responsive defensive genes in V. vinifera

We also predicted various cis-acting regulatory elements of PM and DM-responsive defensive genes by using PlantCare database with default parameters (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/) [81]. To conduct the analysis, an upstream sequence from transcription start site was needed. Nucleotide BLAST was then done between genomic sequence and coding sequence of respective gene. Eventually, different regulatory elements were predicted that might be playing an important role in the regulation of PM and DM-responsive genes of defense pathway.

Availability of data and materials

Not applicable.

Abbreviations

ACT:

Actin

CARE:

Cis-acting Regulatory Elements

CE:

Co-Expression

CNL:

coiled-coil NLRs

DEA:

Differential Expression Analysis

DM:

Downy Mildew

DPI:

days post inoculation

EDS1:

Enhanced Disease Susceptibility

EF1:

Elongation Factor 1

ETI:

Effector Triggered Immunity

ET:

Ethylene

TPM:

Transcript per kilobase million

GO:

Gene Ontology

GSDS:

Gene Structure Display Server

HPI:

hours post inoculation

HR:

Hypersensitive Response

JA:

Jasmonate

LEA2:

Late embryogenesis abundant

MAPK:

Mitogen activated protein kinase

MEME:

Multiple Expectation Maximization for Motif Elicitation

NLR:

Nucleotide-Binding Leucine Rich Repeat Receptor

NDR1:

Non-Race-specific Disease Resistance

NPR:

Nonexpressor of PR

nr:

Non Redundant

PAD4:

Phytoalexin-deficient 4

PAMP:

Pathogen-associated Molecular Patterns

PBS2:

avrPphB susceptible 2

pI:

Isoelectric point

PM:

Powdery Mildew

PR:

Pathogenesis Related

PRR:

Pattern Recognition Receptors

PTI:

Pathogen Triggered Immunity

qPCR:

quantitative real-time PCR

R:

Resistance

RAR1:

Required for Mla-specified Resistance

RSEM:

RNA-Seq by Expectation-Maximization

SA:

Salicylic Acid

SAG101:

Senescence associated gene 101

SAR:

Systemic Acquired Resistance

SRA:

Sequence Read Archive

TF:

Transcription Factor

TNL:

TIR-type NLRs

References

  1. 1.

    Gessler C, Pertot I, Perazzolli M. Plasmopara viticola: a review of knowledge on downy mildew of grapevine and effective disease management. 2011.

  2. 2.

    Zipfel C. Plant pattern-recognition receptors. Trends Immunol. 2014;35(7):345–51. https://doi.org/10.1016/j.it.2014.05.004.

    CAS  Article  PubMed  Google Scholar 

  3. 3.

    Dangl JL, Horvath DM, Staskawicz BJ. Pivoting the plant immune system from dissection to deployment. Science. 2013;341(80):746–51.

    CAS  Article  Google Scholar 

  4. 4.

    Karasov TL, Horton MW, Bergelson J. Genomic variability as a driver of plant–pathogen coevolution? Curr Opin Plant Biol. 2014;18:24–30. https://doi.org/10.1016/j.pbi.2013.12.003.

    Article  PubMed  PubMed Central  Google Scholar 

  5. 5.

    Cui H, Tsuda K, Parker JE. Effector-triggered immunity: from pathogen perception to robust defense. Annu Rev Plant Biol. 2015;66(1):487–511. https://doi.org/10.1146/annurev-arplant-050213-040012.

    CAS  Article  PubMed  Google Scholar 

  6. 6.

    Jacob F, Vernaldi S, Maekawa T. Evolution and conservation of plant NLR functions. Front Immunol. 2013;4:297. https://doi.org/10.3389/fimmu.2013.00297.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  7. 7.

    Jones JDG, Dangl JL. The plant immune system. Nature. 2006;444(7117):323–9. https://doi.org/10.1038/nature05286.

    CAS  Article  PubMed  Google Scholar 

  8. 8.

    Martin GB, Bogdanove AJ, Sessa G. Understanding the functions of plant disease resistance proteins. Annu Rev Plant Biol. 2003;54(1):23–61. https://doi.org/10.1146/annurev.arplant.54.031902.135035.

    CAS  Article  PubMed  Google Scholar 

  9. 9.

    Nimchuk Z, Eulgem T, Holt Iii BF, Dangl JL. Recognition and response in the plant immune system. Annu Rev Genet. 2003;37(1):579–609. https://doi.org/10.1146/annurev.genet.37.110801.142628.

    CAS  Article  PubMed  Google Scholar 

  10. 10.

    Buscaill P, Rivas S. Transcriptional control of plant defence responses. Curr Opin Plant Biol. 2014;20:35–46. https://doi.org/10.1016/j.pbi.2014.04.004.

    CAS  Article  PubMed  Google Scholar 

  11. 11.

    Bari R, Jones JDG. Role of plant hormones in plant defence responses. Plant Mol Biol. 2009;69(4):473–88. https://doi.org/10.1007/s11103-008-9435-0.

    CAS  Article  PubMed  Google Scholar 

  12. 12.

    Howe GA, Jander G. Plant immunity to insect herbivores. Annu Rev Plant Biol. 2008;59(1):41–66. https://doi.org/10.1146/annurev.arplant.59.032607.092825.

    CAS  Article  PubMed  Google Scholar 

  13. 13.

    Katagiri F, Tsuda K. Understanding the plant immune system. Mol Plant-Microbe Interact. 2010;23(12):1531–6. https://doi.org/10.1094/MPMI-04-10-0099.

    CAS  Article  PubMed  Google Scholar 

  14. 14.

    Pieterse CMJ, Leon-Reyes A, Van der Ent S, Van Wees SCM. Networking by small-molecule hormones in plant immunity. Nat Chem Biol. 2009;5(5):308–16. https://doi.org/10.1038/nchembio.164.

    CAS  Article  PubMed  Google Scholar 

  15. 15.

    Pieterse CMJ, Van der Does D, Zamioudis C, Leon-Reyes A, Van Wees SCM. Hormonal modulation of plant immunity. Annu Rev Cell Dev Biol. 2012;28(1):489–521. https://doi.org/10.1146/annurev-cellbio-092910-154055.

    CAS  Article  PubMed  Google Scholar 

  16. 16.

    Vlot AC, Dempsey DA, Klessig DF. Salicylic acid, a multifaceted hormone to combat disease. Annu Rev Phytopathol. 2009;47(1):177–206. https://doi.org/10.1146/annurev.phyto.050908.135202.

    CAS  Article  PubMed  Google Scholar 

  17. 17.

    Mishina TE, Zeier J. Pathogen-associated molecular pattern recognition rather than development of tissue necrosis contributes to bacterial induction of systemic acquired resistance in Arabidopsis. Plant J. 2007;50(3):500–13. https://doi.org/10.1111/j.1365-313X.2007.03067.x.

    CAS  Article  PubMed  Google Scholar 

  18. 18.

    Day B, Dahlbeck D, Staskawicz BJ. NDR1 interaction with RIN4 mediates the differential activation of multiple disease resistance pathways in Arabidopsis. Plant Cell. 2006;18(10):2782–91. https://doi.org/10.1105/tpc.106.044693.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  19. 19.

    Wiermer M, Feys BJ, Parker JE. Plant immunity: the EDS1 regulatory node. Curr Opin Plant Biol. 2005;8(4):383–9. https://doi.org/10.1016/j.pbi.2005.05.010.

    CAS  Article  PubMed  Google Scholar 

  20. 20.

    Knepper C, Savory EA, Day B. The role of NDR1 in pathogen perception and plant defense signaling. Plant Signal Behav. 2011;6(8):1114–6. https://doi.org/10.4161/psb.6.8.15843.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  21. 21.

    Heidrich K, Wirthmueller L, Tasset C, Pouzet C, Deslandes L, Parker JE. Arabidopsis EDS1 connects pathogen effector recognition to cell compartment–specific immune responses. Science. 2011;80(334):1401–4.

    Article  Google Scholar 

  22. 22.

    Selote D, Shine MB, Robin GP, Kachroo A. Soybean NDR 1-like proteins bind pathogen effectors and regulate resistance signaling. New Phytol. 2014;202(2):485–98. https://doi.org/10.1111/nph.12654.

    CAS  Article  PubMed  Google Scholar 

  23. 23.

    Wang J, Shine MB, Gao Q-M, Navarre D, Jiang W, Liu C, et al. Enhanced disease susceptibility1 mediates pathogen resistance and virulence function of a bacterial effector in soybean. Plant Physiol. 2014;165(3):1269–84. https://doi.org/10.1104/pp.114.242495.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  24. 24.

    Bhattacharjee S, Halane MK, Kim SH, Gassmann W. Pathogen effectors target Arabidopsis EDS1 and alter its interactions with immune regulators. Science. 2011;80(334):1405–8.

    Article  Google Scholar 

  25. 25.

    Feys BJ, Moisan LJ, Newman M, Parker JE. Direct interaction between the Arabidopsis disease resistance signaling proteins, EDS1 and PAD4. EMBO J. 2001;20(19):5400–11. https://doi.org/10.1093/emboj/20.19.5400.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  26. 26.

    Feys BJ, Wiermer M, Bhat RA, Moisan LJ, Medina-Escobar N, Neu C, et al. Arabidopsis SENESCENCE-ASSOCIATED GENE101 stabilizes and signals within an ENHANCED DISEASE SUSCEPTIBILITY1 complex in plant innate immunity. Plant Cell. 2005;17(9):2601–13. https://doi.org/10.1105/tpc.105.033910.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  27. 27.

    Wagner S, Stuttmann J, Rietz S, Guerois R, Brunstein E, Bautor J, et al. Structural basis for signaling by exclusive EDS1 heteromeric complexes with SAG101 or PAD4 in plant innate immunity. Cell Host Microbe. 2013;14(6):619–30. https://doi.org/10.1016/j.chom.2013.11.006.

    CAS  Article  PubMed  Google Scholar 

  28. 28.

    Muskett PR, Kahn K, Austin MJ, Moisan LJ, Sadanandom A, Shirasu K, et al. Arabidopsis RAR1 exerts rate-limiting control of R gene–mediated defenses against multiple pathogens. Plant Cell. 2002;14(5):979–92. https://doi.org/10.1105/tpc.001040.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  29. 29.

    Durrant WE, Dong X. Systemic acquired resistance. Annu Rev Phytopathol. 2004;42(1):185–209. https://doi.org/10.1146/annurev.phyto.42.040803.140421.

    CAS  Article  PubMed  Google Scholar 

  30. 30.

    Backer R, Naidoo S, Van den Berg N. The nonexpressor of pathogenesis-related genes 1 (NPR1) and related family: mechanistic insights in plant disease resistance. Front Plant Sci. 2019;10:102. https://doi.org/10.3389/fpls.2019.00102.

    Article  PubMed  PubMed Central  Google Scholar 

  31. 31.

    Zhang Y, Cheng YT, Qu N, Zhao Q, Bi D, Li X. Negative regulation of defense responses in Arabidopsis by two NPR1 paralogs. Plant J. 2006;48(5):647–56. https://doi.org/10.1111/j.1365-313X.2006.02903.x.

    CAS  Article  PubMed  Google Scholar 

  32. 32.

    Loake G, Grant M. Salicylic acid in plant defence—the players and protagonists. Curr Opin Plant Biol. 2007;10(5):466–72. https://doi.org/10.1016/j.pbi.2007.08.008.

    CAS  Article  PubMed  Google Scholar 

  33. 33.

    Mou Z, Fan W, Dong X. Inducers of plant systemic acquired resistance regulate NPR1 function through redox changes. Cell. 2003;113(7):935–44. https://doi.org/10.1016/S0092-8674(03)00429-X.

    CAS  Article  PubMed  Google Scholar 

  34. 34.

    Eulgem T, Somssich IE. Networks of WRKY transcription factors in defense signaling. Curr Opin Plant Biol. 2007;10(4):366–71. https://doi.org/10.1016/j.pbi.2007.04.020.

    CAS  Article  PubMed  Google Scholar 

  35. 35.

    Wang W, Withers J, Li H, Zwack PJ, Rusnac D-V, Shi H, et al. Structural basis of salicylic acid perception by Arabidopsis NPR proteins. Nature. 2020:1–6.

  36. 36.

    Van Loon LC, Pierpoint WS, Boller TH, Conejero V. Recommendations for naming plant pathogenesis-related proteins. Plant Mol Biol Report. 1994;12(3):245–64. https://doi.org/10.1007/BF02668748.

    Article  Google Scholar 

  37. 37.

    van Loon LC, Rep M, Pieterse CMJ. Significance of inducible defense-related proteins in infected plants. Annu Rev Phytopathol. 2006;44(1):135–62. https://doi.org/10.1146/annurev.phyto.44.070505.143425.

    CAS  Article  PubMed  Google Scholar 

  38. 38.

    Wildermuth MC, Dewdney J, Wu G, Ausubel FM. Isochorismate synthase is required to synthesize salicylic acid for plant defence. Nature. 2001;414(6863):562–5. https://doi.org/10.1038/35107108.

    CAS  Article  PubMed  Google Scholar 

  39. 39.

    Gaffney T, Friedrich L, Vernooij B, Negrotto D, Nye G, Uknes S, et al. Requirement of salicylic acid for the induction of systemic acquired resistance. Science. 1993;80(261):754–6.

    Article  Google Scholar 

  40. 40.

    Goyal N, Bhatia G, Sharma S, Garewal N, Upadhyay A, Upadhyay SK, et al. Genome-wide characterization revealed role of NBS-LRR genes during powdery mildew infection in Vitis vinifera. Genomics. 2019.

  41. 41.

    Falk A, Feys BJ, Frost LN, Jones JDG, Daniels MJ, Parker JE. EDS1, an essential component of R gene-mediated disease resistance in Arabidopsis has homology to eukaryotic lipases. Proc Natl Acad Sci. 1999;96(6):3292–7. https://doi.org/10.1073/pnas.96.6.3292.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  42. 42.

    Zhou J-M, Trifa Y, Silva H, Pontier D, Lam E, Shah J, et al. NPR1 differentially interacts with members of the TGA/OBF family of transcription factors that bind an element of the PR-1 gene required for induction by salicylic acid. Mol Plant-Microbe Interact. 2000;13(2):191–202. https://doi.org/10.1094/MPMI.2000.13.2.191.

    CAS  Article  PubMed  Google Scholar 

  43. 43.

    Puranik S, Sahu PP, Srivastava PS, Prasad M. NAC proteins: regulation and role in stress tolerance. Trends Plant Sci. 2012;17(6):369–81. https://doi.org/10.1016/j.tplants.2012.02.004.

    CAS  Article  PubMed  Google Scholar 

  44. 44.

    Singh KB, Foley RC, Oñate-Sánchez L. Transcription factors in plant defense and stress responses. Curr Opin Plant Biol. 2002;5(5):430–6. https://doi.org/10.1016/S1369-5266(02)00289-3.

    CAS  Article  PubMed  Google Scholar 

  45. 45.

    Fitzgerald HA, Chern M-S, Navarre R, Ronald PC. Overexpression of (at) NPR1 in rice leads to a BTH-and environment-induced lesion-mimic/cell death phenotype. Mol Plant-Microbe Interact. 2004;17(2):140–51. https://doi.org/10.1094/MPMI.2004.17.2.140.

    CAS  Article  PubMed  Google Scholar 

  46. 46.

    Makandar R, Essig JS, Schapaugh MA, Trick HN, Shah J. Genetically engineered resistance to fusarium head blight in wheat by expression of Arabidopsis NPR1. Mol Plant-Microbe Interact. 2006;19(2):123–9. https://doi.org/10.1094/MPMI-19-0123.

    CAS  Article  PubMed  Google Scholar 

  47. 47.

    Lin W-C, Lu C-F, Wu J-W, Cheng M-L, Lin Y-M, Yang N-S, et al. Transgenic tomato plants expressing the Arabidopsis NPR1 gene display enhanced resistance to a spectrum of fungal and bacterial diseases. Transgenic Res. 2004;13(6):567–81. https://doi.org/10.1007/s11248-004-2375-9.

    CAS  Article  PubMed  Google Scholar 

  48. 48.

    Malnoy M, Jin Q, Borejsza-Wysocka EE, He SY, Aldwinckle HS. Overexpression of the apple MpNPR1 gene confers increased disease resistance in Malus× domestica. Mol Plant-Microbe Interact. 2007;20(12):1568–80. https://doi.org/10.1094/MPMI-20-12-1568.

    CAS  Article  PubMed  Google Scholar 

  49. 49.

    Chen X-K, Zhang J-Y, Zhang Z, Du X-L, Du B-B, Qu S-C. Overexpressing MhNPR1 in transgenic Fuji apples enhances resistance to apple powdery mildew. Mol Biol Rep. 2012;39(8):8083–9. https://doi.org/10.1007/s11033-012-1655-3.

    CAS  Article  PubMed  Google Scholar 

  50. 50.

    Van Dam S, Vosa U, van der Graaf A, Franke L, de Magalhaes JP. Gene co-expression analysis for functional classification and gene–disease predictions. Brief Bioinform. 2018;19(4):575–92. https://doi.org/10.1093/bib/bbw139.

    CAS  Article  PubMed  Google Scholar 

  51. 51.

    Amar D, Safer H, Shamir R. Dissection of regulatory networks that are altered in disease via differential co-expression. PLoS Comput Biol. 2013;9(3):e1002955. https://doi.org/10.1371/journal.pcbi.1002955.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  52. 52.

    Hu R, Qiu X, Glazko G, Klebanov L, Yakovlev A. Detecting intergene correlation changes in microarray analysis: a new approach to gene selection. BMC Bioinformatics. 2009;10(1):1–9. https://doi.org/10.1186/1471-2105-10-20.

    CAS  Article  Google Scholar 

  53. 53.

    Kostka D, Spang R. Finding disease specific alterations in the co-expression of genes. Bioinformatics. 2004;20 suppl_1:i194–i199.

  54. 54.

    Hudson NJ, Reverter A, Dalrymple BP. A differential wiring analysis of expression data correctly identifies the gene containing the causal mutation. PLoS Comput Biol. 2009;5(5):e1000382. https://doi.org/10.1371/journal.pcbi.1000382.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  55. 55.

    Xue Z, Huang K, Cai C, Cai L, Jiang C, Feng Y, et al. Genetic programs in human and mouse early embryos revealed by single-cell RNA sequencing. Nature. 2013;500(7464):593–7. https://doi.org/10.1038/nature12364.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  56. 56.

    Pierson E, Koller D, Battle A, Mostafavi S. Consortium Gte. Sharing and specificity of co-expression networks across 35 human tissues. PLoS Comput Biol. 2015;11(5):e1004220. https://doi.org/10.1371/journal.pcbi.1004220.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  57. 57.

    Gebrie SA. Promoter analysis of a powdery mildew induced Vitis vinifera NAC transcription factor gene. Adv Crop Sci Tech. 2017;5(01):2. https://doi.org/10.4172/2329-8863.1000252.

    CAS  Article  Google Scholar 

  58. 58.

    Kim HS, Delaney TP. Over-expression of TGA5, which encodes a bZIP transcription factor that interacts with NIM1/NPR1, confers SAR-independent resistance in Arabidopsis thaliana to Peronospora parasitica. Plant J. 2002;32(2):151–63. https://doi.org/10.1046/j.1365-313X.2001.01411.x.

    CAS  Article  PubMed  Google Scholar 

  59. 59.

    Wang M, Vannozzi A, Wang G, Liang Y-H, Tornielli GB, Zenoni S, et al. Genome and transcriptome analysis of the grapevine (Vitis vinifera L.) WRKY gene family. Hortic Res. 2014;1:14016.

    Article  Google Scholar 

  60. 60.

    Marchler-Bauer A, Derbyshire MK, Gonzales NR, Lu S, Chitsaz F, Geer LY, et al. CDD: NCBI’s conserved domain database. Nucleic Acids Res. 2014;43:D222–6.

    Article  Google Scholar 

  61. 61.

    Jones P, Binns D, Chang H-Y, Fraser M, Li W, McAnulla C, et al. InterProScan 5: genome-scale protein function classification. Bioinformatics. 2014;30(9):1236–40. https://doi.org/10.1093/bioinformatics/btu031.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  62. 62.

    Letunic I, Doerks T, Bork P. SMART: recent updates, new developments and status in 2015. Nucleic Acids Res. 2015;43(D1):D257–60. https://doi.org/10.1093/nar/gku949.

    CAS  Article  PubMed  Google Scholar 

  63. 63.

    Hu B, Jin J, Guo A-Y, Zhang H, Luo J, Gao G. GSDS 2.0: an upgraded gene feature visualization server. Bioinformatics. 2014;31:1296–7.

    Article  Google Scholar 

  64. 64.

    Bailey TL, Elkan C. Fitting a mixture model by expectation maximization to discover motifs in bipolymers. 1994.

  65. 65.

    Finn RD, Bateman A, Clements J, Coggill P, Eberhardt RY, Eddy SR, et al. Pfam: the protein families database. Nucleic Acids Res. 2013;42:D222–30.

    Article  Google Scholar 

  66. 66.

    Gasteiger E, Hoogland C, Gattiker A, Wilkins MR, Appel RD, Bairoch A. Protein Identification and Analysis Tools on the ExPASy Server. In: Protein identification and analysis tools on the ExPASy server. In: The proteomics protocols handbook. Springer; 2005. p. 571–607.

    Google Scholar 

  67. 67.

    Thompson JD, Gibson TJ, Higgins DG. Multiple sequence alignment using ClustalW and ClustalX. Curr Protoc Bioinforma. 2003:2–3.

  68. 68.

    Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33(7):1870–4. https://doi.org/10.1093/molbev/msw054.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  69. 69.

    Saitou N, Nei M. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987;4(4):406–25. https://doi.org/10.1093/oxfordjournals.molbev.a040454.

    CAS  Article  Google Scholar 

  70. 70.

    Felsenstein J. Confidence limits on phylogenies: an approach using the bootstrap. Evolution (N Y). 1985;39:783–91.

    Google Scholar 

  71. 71.

    Letunic I, Bork P. Interactive tree of life (iTOL) v5: an online tool for phylogenetic tree display and annotation. Nucleic Acids Res. 2021;49(W1):W293–6. https://doi.org/10.1093/nar/gkab301.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  72. 72.

    Jin J, Tian F, Yang D-C, Meng Y-Q, Kong L, Luo J, et al. PlantTFDB 4.0: toward a central hub for transcription factors and regulatory interactions in plants. Nucleic Acids Res. 2016;:gkw982.

  73. 73.

    Toffolatti SL, De Lorenzis G, Costa A, Maddalena G, Passera A, Bonza MC, et al. Unique resistance traits against downy mildew from the center of origin of grapevine (Vitis vinifera). Sci Rep. 2018;8(1):12523. https://doi.org/10.1038/s41598-018-30413-w.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  74. 74.

    Haas BJ, Papanicolaou A, Yassour M, Grabherr M, Blood PD, Bowden J, et al. De novo transcript sequence reconstruction from RNA-seq using the trinity platform for reference generation and analysis. Nat Protoc. 2013;8(8):1494–512. https://doi.org/10.1038/nprot.2013.084.

    CAS  Article  PubMed  Google Scholar 

  75. 75.

    Seo J, Gordish-Dressman H, Hoffman EP. An interactive power analysis tool for microarray hypothesis testing and generation. Bioinformatics. 2006;22(7):808–14. https://doi.org/10.1093/bioinformatics/btk052.

    CAS  Article  PubMed  Google Scholar 

  76. 76.

    Bastian M, Heymann S, Jacomy M. Gephi: an open source software for exploring and manipulating networks. Proceedings of the International AAAI Conference on Web and Social Media: In; 2009.

    Google Scholar 

  77. 77.

    Ghawana S, Paul A, Kumar H, Kumar A, Singh H, Bhardwaj PK, et al. An RNA isolation system for plant tissues rich in secondary metabolites. BMC Res Notes. 2011;4(1):85. https://doi.org/10.1186/1756-0500-4-85.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  78. 78.

    Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2− ΔΔCT method. methods. 2001;25:402–8.

    CAS  Article  Google Scholar 

  79. 79.

    Pfaffl MW, Horgan GW, Dempfle L. Relative expression software tool (REST©) for group-wise comparison and statistical analysis of relative expression results in real-time PCR. Nucleic Acids Res 2002;30:e36–e36, 9, 336, DOI: https://doi.org/10.1093/nar/30.9.e36.

  80. 80.

    Conesa A, Götz S. Blast2GO: a comprehensive suite for functional analysis in plant genomics. Int J Plant Genomics. 2008;2008:1–12. https://doi.org/10.1155/2008/619832.

    CAS  Article  Google Scholar 

  81. 81.

    Lescot M, Déhais P, Thijs G, Marchal K, Moreau Y, Van de Peer Y, et al. PlantCARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences. Nucleic Acids Res. 2002;30(1):325–7. https://doi.org/10.1093/nar/30.1.325.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

N Goyal and G Bhatia are thankful to DBT, India and ICMR, India respectively for providing junior and senior research fellowships.

Funding

Dr. Kashmir Singh is thankful to SERB, India (grant number EEQ/2016/000103) to provide research grant to support purchase of instruments, reagents, chemicals etc. to carry out this work.

Author information

Affiliations

Authors

Contributions

KS; conceived the idea, designed the experiments, analyzed the results and finalized the manuscript. N Goyal and GB, collected and analyzed the data, compiled the results. AKU provided the plant material and analyzed the results. N Goyal and N Garewal performed wetlab experiments and wrote the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Kashmir Singh.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1

Protein IDs, chromosomal locations and physicochemical properties of EDS1 (Sheet 1), NDR1 (Sheet 1), PAD4 (Sheet 1), NPR (Sheet 1), RAR1 (Sheet 1) and PR (Sheet 1) defensive genes. The table describes the protein ID, transcript ID, chromosomal position, polypeptide length, molecular weight, isoelectric point and instability index of various defensive gene families.

Additional file 2.

Gene configuration and motif study of various classes of defensive genes. The arrangement of introns and exons is defined by GSDS 2.0 server and motif analysis was conducted with Pfam database and MEME tool (A) Gene structure and conserved domains of EDS1 class of defensive genes. (B) Gene structure and conserved domains of NDR1 class of defensive genes. (C) Gene structure and conserved domains of NPR class of defensive genes. (D) Gene structure and conserved domains of PAD4 class of defensive genes. (E) Gene structure and conserved domains of PR class of defensive genes. (F) Gene structure and conserved domains of RAR1 class of defensive genes. Exons are represented by yellow boxes partitioned by thin intron lines and blue boxes reflect UTRs. Different types of conserved motifs are shown with different colored boxes.

Additional file 3.

List of primers used for qRT-PCR analysis of selected PM and DM-responsive defensive genes as well as Internal control genes.

Additional file 4.

Functional characterization of defensive genes by assigning Gene Ontology (GO) terms using Blast2GO tool. (a) Representation of (a) PM-responsive (b) DM-responsive defensive genes classified on the basis of GO terms enrichment in biological process, cellular component and molecular function categories.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Goyal, N., Bhatia, G., Garewal, N. et al. Identification of defense related gene families and their response against powdery and downy mildew infections in Vitis vinifera. BMC Genomics 22, 776 (2021). https://doi.org/10.1186/s12864-021-08081-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-021-08081-4

Keywords

  • Grapevine
  • Mildew resistance
  • Defense related genes
  • Host-pathogen interactions
  • Biotic stress