Next Article in Journal
Biofumigation for the Management of Fusarium graminearum in a Wheat-Maize Rotation
Next Article in Special Issue
Editorial for the Special Issue: “Fungal Pathogenicity Factors”
Previous Article in Journal
A Broad Spectrum Antiparasitic Activity of Organotin (IV) Derivatives and Its Untargeted Proteomic Profiling Using Leishmania donovani
Previous Article in Special Issue
Secreted Glycosyltransferase RsIA_GT of Rhizoctonia solani AG-1 IA Inhibits Defense Responses in Nicotiana benthamiana
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genome-Wide Identification and Characterisation of Stress-Associated Protein Gene Family to Biotic and Abiotic Stresses of Grapevine

1
College of Horticulture and Landscape Architecture, Southwest University, Chongqing 400716, China
2
Key Laboratory of Horticulture Science for Southern Mountainous Regions, Ministry of Education, Chongqing 400716, China
*
Author to whom correspondence should be addressed.
Pathogens 2022, 11(12), 1426; https://doi.org/10.3390/pathogens11121426
Submission received: 6 October 2022 / Revised: 21 November 2022 / Accepted: 23 November 2022 / Published: 27 November 2022
(This article belongs to the Special Issue Fungal Pathogenicity Factors)

Abstract

:
Grapevine is one of the earliest domesticated fruit crops and prized for its table fruits and wine worldwide. However, the concurrence of a number of biotic/abiotic stresses affects their yield. Stress-associated proteins (SAPs) play important roles in response to both biotic and abiotic stresses in plants. Despite the growing number of studies on the genomic organisation of SAP gene family in various species, little is known about this family in grapevines (Vitis vinifera L.). In this study, a total of 15 genes encoding proteins possessing A20/AN1 zinc-finger were identified based on the analysis of several genomic and proteomic grapevine databases. According to their structural and phylogenetics features, the identified SAPs were classified into three main groups. Results from sequence alignments, phylogenetics, genomics structure and conserved domains indicated that grapevine SAPs are highly and structurally conserved. In order to shed light on their regulatory roles in growth and development, as well as the responses to biotic/abiotic stresses in grapevine, the expression profiles of SAPs were examined in publicly available microarray data. Bioinformatics analysis revealed distinct temporal and spatial expression patterns of SAPs in various tissues, organs and developmental stages, as well as in response to biotic/abiotic stresses. This study provides insight into the evolution of SAP genes in grapevine and may aid in efforts for further functional identification of A20/AN1-type proteins in the signalling cross-talking induced by biotic/abiotic stresses.

1. Introduction

Grapevine (Vitis vinifera L.) is a perennial fruit crop that is cultivated globally with 8110 varieties [United States Department of Agriculture Foreign Agricultural Service 2022, United Nations FAO, 2019] [1]. In 2020, viticulture occupies an area of 7.3 million ha, with around 85 million tonnes of table grapes, 1.4 million tonnes of dried grapes, and 260 million hL of wine being produced worldwide [International Organization of Vine and Wine, OIV, 2019]. However, in the prosperous grapevine industry, there are increasing challenges from viruses, bacteria, fungi, pests and environmental stresses; for reviews see [2,3,4,5]. Under the frame of agroecology, the strategies for viticulture stress management have been shifted from traditional chemical control to antagonistic anti-microorganism screening and key resistant loci mining for genomics-based breeding [6,7]. The boundary between biotic and abiotic stress-induced plant defence responses are sometimes obscure; for instance, the jasmonic acide (JA) dependent default pathway is triggered by salt stress in grapevine. Their modulation by a parallel pathway activated by pathogen-derived factors has been approved [8]. These biotic/abiotic signals overlap in their upstream events (calcium influx, apoplastic alkalinisation and induction of JA signalling), but differ in their downstream output (stilbene metabolism, different phytoalexin synthesis genes and oxidative burst) [8]. Another example comes from V. rupestris affected by O-methylmellein, a conserved pathogenic factor. In this case, O-methylmellein as an amplifier of flg22 shares a specific signature depending on RboH, independent of calcium influx. This means it works neither as a candidate toxin, nor an elicitor activating basal immunity, but as a weapon in interspecies competition by amplifying the basal immune response [9].
Moreover, it is helpful to refer to the classic CODIT (Compartmentalization Of Decay In Trees) and recent findings on the cellular level responses to have a better understanding of host responses to biotic/abiotic stresses [10]. Where the vascular wilt systemically spreads microorganisms or structures such as spores or hyphae, and phytotoxins are mainly restricted to the vessel lumen and cells surrounding vessels, canker is hardly able to spread systemically due to the death of a portion of the vascular cambium that disables the newly formed functional xylem and phloem. Therefore, from the plant side, the resource allocation could be subjected to a trade-off between vessel occlusion and defence reaction in conjunctive living tissues, as well as the signalling transduction cascade within single cells and amongst the symplast through plasmodesmata. The worldwide spread of Grapevine Trunk Disease has been shown to be damaging depending on fluctuating environments as the fungi can be latent in the tissues for years [11]. This may hold true for the other biotic stresses. Therefore, it is interesting to find out about some signalling to diverge the abiotic stresses from the biotic stresses, or branching points upstream or downstream of certain genes.
Among these, stress-associated proteins (SAPs) have attracted much attention. SAPs are a class of zinc-finger proteins containing A20/AN1 domains, which have emerged as important components in plant stress tolerance in the past decade [12,13]. The first identified zinc finger protein was the cytokine-induced protein A20, functioning in the human umbilical vein endothelial cells [14]. This protein is induced by tumour necrosis factor (TNF) and inhibits TNF-induced apoptosis [15]. The other zinc finger domain, AN1, has been characterised for the first time in a ubiquitin-like fusion protein localised in the hemisphere of Xenopus laevis eggs and early embryos [16]. As for the C-terminal AN1 domain, it was first found in the protein encoding the hemispheric parent RNA of Xenopus laevis. The conserved sequence is CX2CX9-12CX1-2CX4CX2HX5HXC [17]. The AN1 zinc-finger domain is usually found associated with the A20 zinc-finger and such types of proteins are present in many eukaryotes [18].
Though A20/AN1 domains were discovered in animals, their roles have been gradually illuminated in plants. Early studies on SAPs have been focused on abiotic stresses at the transcriptional level and gene function. The first plant A20/AN1 protein was cloned from rice (OsiSAP1) and showed enhanced expression in response to cold, drought, salt, submersion, heavy metals and injury [19]. In addition, OsSAP9 (ZFP177) and AtSAP5 were found to be induced by heat, cold and H2O2, and to play an important role in enhancing temperature stress tolerance [20,21]. Evidence from the OsiSAP family as to salinity, osmotic pressure, and hormone and oxidative stresses were reported for a review; see [22,23]. To date, a total of 18 SAPs have been isolated from rice, 14 in Arabidopsis, and 37 in cotton. SAPs have also been identified in other horticulture plant species, such as grapevine, tomato and apple [24,25,26].
SAPs interact with many other proteins to execute their functions [27]. There is direct evidence in the cross-talk between biotic and abiotic stress signalling pathways from the first identified SAP in a plant that OsSAP1 plays a key role in basal resistance against pathogen infection [28]. With respect to grapevine, the SAP gene family has been studied under abiotic stress. The overexpression of VaSAP15 rendered the transgenic plants lower in electrolyte leakage and malondialdehyde levels and higher in antioxidant enzyme activity under cold stress [29]. However, the knowledge of SAP gene family reactions under biotic stress remains largely unknown.
In this study, we identified all VvSAPs genes and explored their evolutionary relationships, sequence features and duplication events. The tissue-specific expression of grapevine SAPs among species under biotic/abiotic stresses were analysed and their potential cross-talking signalling pathways were further discussed. Those will provide references for optimizing the management practices and for mining candidate resistant loci for genomics-based breeding. The signalling pathway induced by SAPs encoding the A20/AN1 zinc protein domain will shed light on further study on the biotic and abiotic stresses inducing signalling cross-talking in grapevine.

2. Materials and Method

2.1. Identification and Sequence Analysis

The nucleotide, protein sequences and GFF files of V. vinifera were downloaded from the grapevine genome database (Vitis_vinifera-Ensembl Genomes 51). The HMM profile of the SAP DNA-binding domain (PF02037) was downloaded from the Pfam protein families database (http://pfam.janelia.org, accessed on 20 November 2022) and used to survey all grapevine proteins. A local BLASTP search was used to query the full-length amino acid sequences of SAPs from Arabidopsis (TAIR-Home Page, arabidopsis.org, accessed on 20 November 2022). To detect conserved domains in each candidate sequence, A20/AN1 zinc-finger motifs were characterised with motif scanning using SMART (http://smart.embl-heidelberg.de/smart/, accessed on 20 November 2022), NCBI-Conserved Domain Database (NCBI Conserved Domain Search (nih.gov, accessed on 20 November 2022)) and InterPro (http://www.ebi.ac.uk/interpro/, accessed on 20 November 2022). Protein sequences containing A20 or AN1 domains were considered members of the candidate SAPs gene family in grapevine, and the coding sequences (CDS) were obtained using GFF files by TBtools1.098763 (South China Agricultural University, Guangzhou, China).
The characteristics of VvSAPs, including protein length, isoelectric point (pI), molecular weight (MW) and grand average of hydropathicity (GRAVY), were predicted by the ExPasy website (http://web.expasy.org/protparam/, accessed on 20 November 2022). The subcellular localisation was predicted by Cell-PLoc-2.0 (http://www.csbio.sjtu.edu.cn/bioinf/Cell-PLoc-2/) (accessed on 11 December 2021).

2.2. Chromosomal Location and Gene Duplication Analysis

The chromosomal positions of candidates were also obtained from V. vinifera genome annotation files (Vitis_vinifera-Ensembl Genomes 51). BLAST results and the GFF files of grapevine genomes were inputted into TBtools1.098763 to identify syntenic relationships and used to exhibit the duplication patterns. Physical chromosomal locations were graphically represented by scaling the 8 chromosomes.

2.3. Multiple Sequences Alignment and Phylogenetic Analysis

To discover the sequence characteristic, the amino acid sequences of VvSAPs were used for multiple sequence alignments through ClustalX 2.1 with default settings. To compare the evolutionary relationships and identify the subfamilies, the SAPs from grapevine, Arabidopsis and rice were used to construct the phylogenetic tree using MEGA7.0 (https://www.megasoftware.net/) (accessed on 12 March 2022) with the neighbour-joining (NJ) method. The bootstrap process was used to detect the tree with 1000 replicates. The phylogenetic tree was then visualised by online software iTol (iTOL: Interactive Tree of Life (embl.de)) (EMBL, Heidelberg, Germany).

2.4. Gene Structure, Domain Composition, Conserved Motif and Promoter Regions

The genomic DNA and corresponding CDS sequences of the putative SAPs, including coding sequence distribution, exon–intron boundaries and intron patterns, were analysed with the Gene Structure Display Server (GSDS) (http://gsds.cbi.pku.edu.cn/, accessed on 15 October 2021). Complete protein sequences were submitted to the NCBI-Conserved Domain Database (NCBI Conserved Domain Search (nih.gov, accessed on 15 December 2021)) to find A20/AN1 domain in VvSAPs. Complete amino acid sequences were analysed with MEME (http://meme.nbcr.net/meme/cgi-bin/meme.cgi, accessed on 15 December 2021) to discover conserved motifs in VvSAPs. The DNA sequences (2 kb) upstream of the initiation codon of each candidate gene were downloaded, and cis-regulatory elements were predicted with Plant CARE (http://bioinfor-matics.psb.ugent.be/webtools/plantcare/html/, accessed on 15 December 2021). The combination of gene structure, domain composition, conserved motif, promoter regions and phylogenetic tree were then generated using TBtools1.098763.

2.5. Expression Profiles Mining of Grapevine SAPs in Various Tissues, Organs and Developmental Stages

The expression profiles of VvSAP genes predicted from an analysis of the grapevine genome releases were analysed in a V. vinifera ‘Corvina’ gene expression atlas of different organs at various developmental stages. Microarray data were downloaded from the Gene Expression Omnibus (GEO) under the series entry GSE36128 (https://www.ncbi.nlm.nih.gov/geo/, accessed on 5 March 2022). Expression analyses in response to biotic and abiotic stresses were based on microarray data from Fung et al. (2008), Cramer et al. (2007) and Tattersall et al. (2007) [30,31,32] (series matrix accession numbers GSE31594, GSE31677 and GSE6404), downloaded from GEO datasets. In addition, the expression analysis dataset of biotic stress in casual agents of fungi disease was obtained from the Sequence Reads Archive (SRA) database under the accession number PRJNA692532 [33]. Data were analysed and graphically represented using the mean of Multi Experiment Viewer software. Data obtained from the V. vinifera ‘Corvina’ expression atlas were normalised based on the mean expression value of each gene in all tissues/organs analysed, and the other data were calculated as log2 fold change in treated and untreated samples.

2.6. Plant Material

A suspension cell culture line of V. rupestris expressing GFP-AtTUB6 was used in this study [34]. Cells were cultivated in a liquid medium containing 4.3 g L−1 Murashige and Skoog salts (Duchefa, Haarlem, The Netherlands), 30 g L−1 sucrose, 200 mg L−1 KH2PO4, 100 mg L−1 inositol, 1 mg L−1 thiamine and 0.2 mg L−1 2, 4-dichlorophenoxyacetic acid (2, 4-D), which was adjusted to a final pH of 5.8 and supplemented with 30 mg L−1 hygromycin B. Cells were subcultured weekly, inoculating 10 mL of stationary cells into 30 mL of fresh medium in 100 mL Erlenmeyer flasks. The cell suspensions were incubated in the dark at 26 °C on an orbital shaker (Shanghai ZhiCheng, Shanghai, China) at 150 rpm [33].

2.7. Fungi Isolation, Inoculation and Transcript Analysis

The culture filtrates of D. seriata 98.1 and N. parvum Bt-67 (isolated from Pyrénées Orientales and Colmar, France, respectively, kindly provided by Prof. Christophe Bertsch, Université de Haute-Alsace, Colmar, France) and B. dothidea Sd-18 (isolated from vineyards in Yantai, China) were generated by cultivating 3 plugs of 2×1 cm of the strain grown on solid medium containing 37 g L−1 potato dextrose broth (PDB, Solarbio, Beijing, China) with 15 g L−1 agar powder (Solarbio, Beijing, China) for 1 week in 250 mL 20 g L−1 malt-medium (Bacto malt-extract, Solarbio, Beijing, China) in a 500 mL Erlenmeyer flask in the dark at 28 °C, 200 rpm for 21 days. Further, the D. seriata 98.1 culture was filtered with sterilized filter papers to remove the big fungi clusters and the liquid was harvested in an Erlenmeyer flask. After the liquid was vacuum-filtered through a 0.22 µm filter, it was transferred into sterilized 20 mL screw capped vials and stored at −20 °C [33].
An amount of 40 µL of culture filtrate of D. seriata 98.1 and N. parvum Bt-67 were added to 10 ml V. rupestris expressing GFP-AtTUB6 suspension cells when subcultured. The cells were incubated in the dark in a constant shaker at 26 °C, 200 rpm, for 3 days. As a contrast, a short-term treatment with 2% (v/v) of the culture filtrate of N. parvum Bt-67 and B. dothidea Sd-18 for 24 h was tested. Sterile water (40 µL) was used in the mock-treated control. V. rupestris expressing GFP-AtTUB6 suspension cells were then collected by filtering through a 30 µm pore size nylon net filter in an aseptic vacuum filter system (Millipore, Burlington, MA, USA). The cells were then frozen in liquid nitrogen for RNA extraction. Two biological replicates, each comprising three technical replicates, were conducted for each treatment. RNA extraction and transcriptome sequencing and analyses refer to Zhao et al. (2021) [33].

3. Results

3.1. Identification of SAP Gene Family in Grapevine

In this study, both Hidden Markov Model (HMM) and BLASTP program were utilized to scan the grapevine database of the 12 X grapevine genome. After the removal of redundant entries, 15 SAP genes with A20/AN1 domain were identified, and numbered as VvSAP1 to VvSAP15 based on their chromosomal locations according to Vannozzi et al. (2021) [35]. The parameters used to describe the VvSAP proteins are listed in Table 1, including deduced protein length, molecular weight, isoelectric point, aliphatic index, and subcellular and chromosome localisations. The deduced length of VvSAPs ranged from 67 (VvSAP14, VvSAP15) to 293 amino acid residues (VvSAP8), whereas the pI ranged from 7.99 (VvSAP5, VvSAP7) to 10.26 (VvSAPs14). The results of the predicted subcellular localisations of these proteins revealed that all 15 proteins may be located in the nucleus, and one also in the chloroplast.
The fifteen candidate genes were subjected to SMART, NCBI-Conserved Domain Database and InterPro analysis for confirmation. All candidates contained the conserved zinc-finger domain AN1, and most also contained the zinc-finger domain A20, allowing validation of all fifteen members of the SAP gene family (Figure 1). Among the fifteen SAPs, nine encoded proteins contained both N-terminal A20 and C-terminal AN1 domains, which is the typical configuration for plant SAPs. Four SAPs contained only a single AN1 domain, while two SAPs contained two AN1 domains. Interestingly, no SAPs with only a single A20 domain were found in V. vinifera. This is similar to the findings in Arabidopsis, cotton and tomato, but not in rice.

3.2. Chromosomal Distribution of SAPs

The chromosomal locations of each gene are shown in Figure 2. Three SAPs (VvSAP6, VvSAP7 and VvSAP8) are present on chromosome 8, and two SAPs (VvSAP2 and VvSAP3, VvSAP4 and VvSAP5, VvSAP9 and VvSAP10, VvSAP12 and VvSAP13, and VvSAP14 and VvSAP15) are present on chromosome 2, chromosome 6, chromosome 13, chromosome 16 and chromosome 18, respectively. One SAP (VvSAP1, VvSAP11) is present on chromosome 1 and chromosome 14, respectively.

3.3. Multiple Sequence Alignments and Phylogenetic Analysis of SAPs

In order to understand the genetic relationship and evolutionary characteristics of grapevine SAP members, and speculate as to their biological function, an un-rooted phylogenetic tree was constructed by using aligned SAP sequences of grapevine, Arabidopsis and rice. As shown in Figure 3, the SAPs from the three species could be divided into three groups, which indicates that the plant SAPs are relatively conservative during evolution. Group I contained two AtSAPs and three VvSAPs. Group II contained seven AtSAPs, six VvSAPs and four OsSAPs. Group III was composed of five AtSAPs, fourteen OsSAPs and six VvSAPs. Group II and Group III included more than half of the total SAPs, ranking the largest grapevine SAPs groups. The grapevine SAPs categorized as Group I members were A20-AN1-type SAPs (VvSAP1 and VvSAP13) and a single AN1 domain (VvSAP12). The grapevine SAPs in group II were almost all the A20-AN1-type (VvSAP2, VvSAP5, VvSAP6, VvSAP7 and VvSAP9), with the only exception of VvSAP3 that contained a single AN1 domain. The grapevine SAPs in group III were A20-AN1-type (VvSAP4, VvSAP8, VvSAP11, VvSAP14 and VvSAP15), with the exception of VvSAP10 which contained a single A20 domain.

3.4. Conserved Motif, Typical Domain and Gene Structure Analysis

Gene structure often diverges during the evolution of multigene families. This mechanism of gene evolution facilitates the acquisition of new functions which allow organisms to adapt to various challenging environments. To further investigate the function of SAP genes in grapevine, the structural features of SAPs were analysed (Figure 1). The N-terminal of VvSAP1-VvSAP10 contained conserved motif one, the C-terminal of VvSAP1-VvSAP13 contained conserved motifs two and three, and the C-terminal of VvSAP14 and VvSAP15 contained conserved motifs seven and eight. In addition, VvSAP3, VvSAP6 and VvSAP8 also contained conserved motifs five and six.
Through the multiple sequence alignment of 15 VvSAP protein sequences (Figure 4), it was found that all of the protein sequences in the SAP family contained the AN1 conserved domain, and most of them contained the zf-A20 conserved domain. Specifically, the VvSAP1-VvSAP15 protein sequences contained the ZnF_AN1 conserved domain, and the VvSAP14 and VvSAP15 protein sequences contained the ZnF_AN1 domain, without the zf_A20 domain. The sequences of this gene family were highly conserved, and they may perform similar functions in grapevine.

3.5. Cis-Regulatory Elements in VvSAP Promoters

The sequences 2 kb upstream of the translation initiation site of each SAP were analysed using the PLANTCARE server to discover possible elements involved in mediating gene expression. The cis-regulatory elements found in the SAP promoters are shown in Supplementary Material Table S1. Various cis-elements related to plant abiotic stresses, hormone responses and pathogen defences were found in SAP promoters. These results indicate that SAPs are possibly involved in responses to both biotic and abiotic stresses, and play important roles in plant defence and hormone signalling transduction.

3.6. Expression Patterns of VvSAPs in Various Tissues, Organs and Developmental Stages

The expression patterns of VvSAPs in different tissues, organs and developmental stages were analysed (http://bar.utoronto.ca/efp_grape/cgi-bin/efpWeb.cgi, accessed on 15 March 2022. As shown in Figure 5, nearly all VvSAPs were up-regulated in both vegetative and reproductive organs, whereas VvSAP1, VvSAP12 and VvSAP13 were down-regulated. The declined extent was much stronger in vegetative organs than in reproductive organs. Among these genes, VvSAP9 and VvSAP7 have the highest transcript levels among all tissues. VvSAP4 and VvSAP8 exhibited much higher transcript abundance in pollen than in other flower structures. The consistently down-regulated VvSAP13 showed an up-regulation in seed and berry. The transcriptome data showed that several VvSAPs genes, such as VvSAP2, VvSAP3, VvSAP4, VvSAP6, VvSAP8 and VvSAP11, displayed remarkable transcriptions at certain time points during plant development.

3.7. Involvement of VvSAPs in Biotic Stresses

To clarify the expression pattern of the grapevine SAPs gene family under biotic stress, in this study, the expression levels of the grapevine SAP family in native American grapevines (V. rupestris expressing GFP-AtTUB6), native European grapevine (V. vinifera ‘Cabernet Sauvignon’) and a hybrid species (V. aestivalis × V. vinifera ‘Norton’) challenged by fungi infection were compared (Figure 6) [30,33]. VvSAP10 was shown to be highly induced in V. rupestris, expressing GFP-AtTUB6 inoculated by N. parvum Bt-67, D. seriata 98.1 and B. dothidea Sd-18 at 12 and 72 hpi, respectively (Figure 6A,B). To be specific, the expression levels of VvSAP10 were more obvious in two less virulent fungi, D. seriata 98.1 and B. dothidea Sd-18, than in N. parvum Bt-67. Moreover, the expression of VvSAP3, VsSAP6 and VvSAP11 at 72 hpi were shown to be inverse results between D. seriata 98.1 and N. parvum Bt-67 inoculation (Figure 6B). Therefore, an earlier (from 0-48 hpi) relatively expression pattern (to control) was shown (Figure 6A). Moreover, the expression levels in tissue of V. vinifera ‘Cabernet Sauvignon’ and V. aestivalis × V. vinifera ‘Norton’ without inoculation and after inoculation with Erysiphe necator conidospores were further analysed (Figure 6C). Interestingly, the expression of VvSAP1, VvSAP2 and VvSAP3 reduced during the early time frame but were resilient nearly to the initial level in both V. aestivalis × V. vinifera ‘Norton’ and V. vinifera ‘Cabernet Sauvignon’, although the relative expression levels were higher in V. aestivalis × V. vinifera ‘Norton’. The induction of VvSAP10 presented an earlier induction in V. vinifera ‘Cabernet Sauvignon’ (at 12 hpi) than in V. aestivalis × V. vinifera ‘Norton’ (at 24 hpi). However, the relative expression of VvSAP15, VvSAP4 and VvSAP14 held a similar pattern to VvSAP1, VvSAP2 and VvSAP3, although the final expression levels were higher compared to the initial (Figure 6C).

3.8. Involvement of VvSAPs in Abiotic Stresses

With regard to the grapevine distribution evolutionarily, and its cultivation status nowadays, the abiotic stresses that affect grapevine the most are drought, salt and temperature. In this study, the grapevine SAPs have been analysed in silico on transcript level based on the GEO database [31,32]. Consistent with plant responses to biotic stress, plant responses to abiotic stress were analysed in both the short time (within 24 h) and the long term (after days), respectively. VvSAP6 was the most significantly enhanced gene in all of the tested abiotic stresses. Particularly under osmotic stress (PEG 1-24 hpi, Figure 7A, and water deficit 1-16 dpi, Figure 7B), the VvSAP6 mediated salinity-induced resistance occurred later than osmotic stress (NaCl 1-24 hpi, Figure 7A, and NaCl 4-16 dpi in Figure 7B). The one with obvious expression induced by biotic stress, VvSAP10, was attractive to mine its function in abiotic stresses. Contrary to the VvSAP6, VvSAP10 showed mild changes to the osmotic and salinity, but to the temperature stress as promptly as within 8 h (5 °C, Figure 7A). Moreover, VvSAP9 showed a similar expression pattern to VvSAP6, although the extent was not as prominent as the later one (Figure 7A,B). The rest of the VvSAPs could be categorized into three groups. VvSAP15, VvSAP4 and VvSAP14 were induced continuously from hours to days in all the abiotic stresses. The expression of VvSAP5 and VvSAP8 remained fairly steady, with only a bit of down-regulation at day 12 to osmotic and salinity stresses, respectively. VvSAP1, VvSAP2 and VvSAP3 were firstly decreased at 4 hpi, but came back to the similar level as initially within 1 d and continued to increase (Figure 7A,B).

4. Discussion

Resistance to biotic and abiotic stresses is a central topic for sustainable agriculture, especially in grapevine, one of the world-wide cultivated crops between 30° and 50° N and 30° and 40° S with the highest economic output per acreage. Abiotic stresses strongly affect biotic interactions at various levels, for instance, plant morphogenesis from cells, tissues to organs, transcriptomic regulation, signalling pathways, as well as primary and secondary metabolisms. These factors combine to make an individual plant a more or less suitable host for its counterpart [36,37].
Proteins containing A20/AN1 zinc fingers have been characterised as regulating the immune responses in animals and environmental responses in plants [38]. The grapevine genome database was referred to in this study, and a total of 15 VvSAP gene members was identified and characterised (Table 1). Genome-wide survey of SAP gene family composition was further performed in other plant species. The family members of SAPs in grapevine are more than in maize (11), Arabidopsis (14) and tomato (14); while less than in rice (18), potato (17), tobacco (16) and cotton (37). The conservative domain of the VvSAP family was analysed (Figure 4). It was found that all candidates contained the conserved AN1 zinc-finger domain and mostly also contained the A20 zinc-finger domain, which is similar to findings in Arabidopsis, cotton and tomato, but not in rice. In addition, there were five SAPs in rice and four in Arabidopsis with only the AN1 zinc-finger domain. On the other hand, in the evolutionarily primitive organisms such as Saccharomyces cerevisae and Aspergillus fumigatus, only an AN1 zinc finger domain, being devoid of A20 domain, has been found, and with similar properties in cotton. This may indicate the ancient origin of the AN1 zinc finger domain compared to A20 domain [12,24]. Chromosomal positions of VvSAP were identified using Tbtools1.098763 (Figure 2).
Analysis of the evolutionary relationships of SAPs among grapevine and two other species indicated that these SAPs can be divided into three groups. The grapevine SAPs in group II were mostly with the A20-AN1-type, excepting VvSAP3 containing a single AN1 domain (Figure 3). The conserved motif distributions of VvSAP according to the evolutionary relationship were further examined. A total of 10 conserved motifs were identified, and their distributions exhibited strong evolutionary conservation (Figure 1). Gene structure can also effectively reveal evolutionary relationships among gene families. Zhou et al. (2018) reported that the intronless SAPs could reduce post-transcriptional processing, and therefore can be utilised to test immediate responses to abiotic stresses [39]. According to these findings, SAPs are highly conserved evolutionarily in plants [39]. The chromosomal positions of candidates were also obtained from V. vinifera genome annotation: Table 1. The location of SAP genes in grapevine can help to identify these genes in other species, directly or using microsynteny among them.
Plant SAP genes are thought to play a vital role in regulating a variety of biotic and abiotic stresses according to previous studies [23,40,41]. The most well-characterised plant protein in this class is the rice A20/AN1 zinc-finger domain stress associated protein 1, OsISAP1. Mukhopadhyay et al. (2004) found that the over-expression of OsISAP1 exhibits tolerances to cold, water defection and salinity stress [19]. Moreover, the OsSAPs are responsive to multiple biotic stresses. The overexpression of OsSAP1 enhances the basal resistance against pathogen infection in tobacco [28]. Our results indicated that nearly all VvSAP genes were up-regulated in both vegetative and reproductive organs, whereas VvSAP1, VvSAP12 and VvSAP13 were down-regulated. The decreased extent was much stronger in vegetative organs than in reproductive organs. VvSAP4 and VvSAP8 exhibited a much higher transcript abundance in pollen than in other flower structures (Figure 5). Further, a native American grapevine (V. rupestris), a native European grapevine (V. vinifera ‘Cabernet Sauvignon’) and a hybrid species (V. aestivalis × V. vinifera ‘Norton’) were compared respectively for their biotic stress resistance (to Botryosphaeriaceae and Erysiphaceae) (Figure 6). Given that the expression data at the quantitative level were generated by different methods (microarray Figure 5, Figure 6C and Figure 7 vs. RNAseq Figure 6A,B), the comparisons in this study are restricted within the same method. There are significant different expression patterns (nothing to consider their quantitative level) among American and European grapevine varieties. VvSAP10 was highly induced in V. rupestris. The expression levels were more obvious in two less virulent fungi. D. seriata 98.1 and B. dothidea Sd-18. than in N. parvum Bt-67 (Figure 6A,B) [11]. The induction of VvSAP10 presented an earlier induction in V. vinifera ‘Cabernet Sauvignon’ than in V. aestivalis × V. vinifera ‘Norton’. VvSAP10 showed mild changes to osmotic and salinity stresses, but strong changes to the temperature stress (Figure 7). VvSAP6 was the most significantly enhanced gene against all of the tested abiotic stresses. Especially under osmotic stress, the VvSAP6 mediated salinity-induced resistance occurred later than osmotic stress (Figure 7A,B). The expression pattern of VvSAP1, VvSAP2, VvSAP3, VvSAP5 and VvSAP8 was down-regulated in the short term and then enhanced (Figure 7A,B), which is different to the SAPs in Arabidopsis, cucumber and cotton [20,42].
Temporary water deficit stress could enhance the immunity of Arabidopsis by stomatal closure, thereby enhancing the resistance to Pseudomonas syringae pv. Tomato DC3000. In this case, a downstream gene transcription level has confirmed that the E3 ligase SINAT4 down-regulated the cysteine protein RD21A [43]. Our previous study on the actin filament (AF) marker line in grapevine indicates that guard cells act as pacemakers of defence, dominating the responses of the remaining epidermal cells. It has been shown that to the phytopathogenic bacteria P. syringae pv. Tomato DC3000, actin responses were confined to the guard cells. In contrast, upon contact with zoospores of Plasmopara viticola, not only the guard cells, but also epidermal pavement cells where no zoospores had attached, responded with the formation of a perinuclear actin basket [44].
The cold stress normally inhibits the JA-dependent signalling pathway but activates the salicylic acid (SA)-dependent signalling pathway; the high temperature brings a contrary consequence [36]. Although the gene expression levels induced by abiotic stresses are normally much milder than the biotic stress induction (Figure 6 and Figure 7), it probably can affect the activation of certain resistance genes in plants challenged by further attacks of pathogenic microorganisms [36]. Paolinelli-Alfonso et al. (2016) found that transcriptome analysis of L. theodromae in the presence of grapevine wood revealed an upregulation of genes involved in the SA pathway [45]. This upregulation was enhanced by high temperature. The high temperature-induced contradictory SA signalling pathway can be understand by taking account of the biotic induced signalling. An activation of L-tyrosine catabolism pathway could lead to inhibition of the SA pathway, favouring fungal development especially during heat stress [45].
On the other hand, high JA concentration is known to have an antagonistic effect on the SA pathway [46]. Botryosphaeriaceae are known to produce JA, which has been identified in the culture filtrate of some GTD agents, e.g., N. parvum and N. vitifusiforme [47]. Hence it could be hypothesised that, in addition to a toxic effect, JA inhibits the SA pathway involved in the stimulation of plant defences. Zhao et al. (2021) have reported a downregulation of VvWRKY70 [30], a putative central component of the SA pathway, after treatment of V. rupestris cells with D. seriata secreted compounds.
These signalling cascades are partially sharing with abiotic induced defence responses, although it is to be noted that this is a matter of timing events. Otherwise, a consequence-centred explanation may lead to an apparent discrepancy. For instance, under cold tolerance study, a transient disassembly of microtubules followed by their recovery with progressive acclimation has been approved to be a kind of cold acclimation, which is obviously contradictory if to follow the final consequence that a stable microtubules array was established for cold challenging [48]. Moreover, the above discussed JA is necessary to convey efficient disassembly of MTs under cold stress, but that is not the exclusive signal. Another example comes from the cytoskeleton (actin filaments and microtubules) early arrangement (within a time frame of 20 min vs. long-time adaptation, e.g., hours or days) to distinguish the PTI from ETI [49]. This evidence could come back to answer the basic question, that environmental mediated biotic stress symptoms are not always coincidental with the invasion of casual agents.
The above findings indicate crosstalk and convergence of mechanisms in these pathways and the existence of a general stress response, suggesting that resource allocation could be subjected to a trade-off between vessel occlusion and defence reaction in conjunctive living tissues in grapevine, as well as the signalling transduction cascade within single cells and amongst the symplasts through plasmodesmata. The principle of the CODIT lies in the compartmentalisation of the established four obstacles that restrict and decay pathogen spread [50]. The first three obstacles located in xylem per se can be interpreted as the reaction zone after pathogen invasion. Obstacle four, consisting of living cells, e.g., paratracheal parenchyma, fibres and ray, can protect the host via arming themselves with newly formed structures. The mystery of environment-related grapevine disease, such as the disruption of water flow, can probably be explained via vessel dependence systemically spread, as well as local living cell signalling transduction against these four obstacles [10]. VaSAP15 from a Chinese wild grapevine V. amurensis has been confirmed to localise in both the nucleus and cell membrane, which prompts a pending question during the penetration of fungi, that the contraction of the actin around the penetration site probably linked with the actin nucleus basket to attract the nucleus towards the penetration site [29,51]. Plant-specific class XIV kinesin KCH1 has been confirmed for its function as a tether [9,52]. Therefore, the VvSAPs are probably triggered at the early time point via an acting-dependent endocytotic of a receptor on the member.

5. Conclusions

In the present study, the systematic identification and characterisation of SAP gene family members in grapevine were studied. Their characteristics, including genomic locations, gene duplications, evolutionary relationships, conserved domains and motifs, gene structures and cis-elements in promoter regions were analysed based on bioinformatics methods. In addition, the regulations of VvSAPs in biotic and abiotic stresses were investigated based on transcriptomic analysis and discussed under the frame of signalling cross-talking at cellular and organic level. On the one hand, since plant response to abiotic stress must be coordinated with growth and development, such a host physiological status further confers susceptibility or resistance through biotic/abiotic signalling cross-talking under biotic stress; on the other hand, endophytic and rhizospheric microorganisms can enhance plant resistance to abiotic stresses, and understanding the biotic/abiotic signalling cross-talking can be helpful for the utilisation of biocontrol. This study provides essential information for elucidating the potential roles of VvSAP in response to biotic and abiotic stresses. Further studies focus on the respective roles of the different SAPs investigated via transgenic overexpression or knock-outs and can provide direct evidences and potential implications for grapevine breeding.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/pathogens11121426/s1, Table S1: The cis-elements in the VvSAP promoters.

Author Contributions

X.G. conceived the study. X.S. performed the experiments. X.S. and X.X. wrote the manuscript. X.S., X.X. and X.G. interpreted the experimental data and revised the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by projects of the National Natural Science Foundation of China 31901990, and the Natural Science Foundation of Chongqing, China cstc2021jcyj-msxmX0994.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors have no conflicts of interest related to the work described in the manuscript.

References

  1. Tu, M.; Wang, X.; Yin, W.; Wang, Y.; Li, Y.; Zhang, G.; Li, Z.; Song, J.; Wang, X. Grapevine VlbZIP30 improves drought resistance by directly activating VvNAC17 and promoting lignin biosynthesis through the regulation of three peroxidase genes. Hortic Res. 2020, 7, 150–165. [Google Scholar] [CrossRef] [PubMed]
  2. Fuchs, M. Grapevine viruses: A multitude of diverse species with simple but overall poorly adopted management solutions in the vineyard. J. Plant Pathol. 2020, 102, 643–653. [Google Scholar] [CrossRef]
  3. Mondello, V.; Songy, A.; Battiston, E.; Pinto, C.; Coppin, C.; Trotel-Aziz, P.; Clement, C.; Mugnai, L.; Fontaine, F. Grapevine trunk diseases: A review of fifteen years of trials for their control with chemicals and biocontrol agents. Plant Dis. 2018, 102, 1189–1217. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Pertot, I.; Caffi, T.; Rossi, V.; Mugnai, L.; Hoffmann, C.; Grando, M.S.; Gary, C.; Lafond, D.; Duso, C.; Thiery, D.; et al. A critical review of plant protection tools for reducing pesticide use on grapevine and new perspectives for the implementation of IPM in viticulture. Crop Prot. 2017, 97, 70–84. [Google Scholar] [CrossRef]
  5. Petrussa, E.; Braidot, E.; Zancani, M.; Peresson, C.; Bertolini, A.; Patui, S.; Vianello, A. Plant flavonoids—Biosynthesis, transport and involvement in stress responses. Int. J. Mol. Sci. 2013, 14, 14950–14973. [Google Scholar] [CrossRef] [Green Version]
  6. Yin, X.; Li, T.; Jiang, X.; Tang, X.; Zhang, J.; Yuan, L.; Wei, Y. Suppression of grape white rot caused by coniella vitis ssing the potential biocontrol agent bacillus velezensis GSBZ09. Pathogens 2022, 11, 248. [Google Scholar] [CrossRef] [PubMed]
  7. Pirrello, C.; Zeilmaker, T.; Bianco, L.; Giacomelli, L.; Moser, C.; Vezzulli, S. Mining grapevine downy mildew susceptibility genes: A Resource for genomics-based breeding and tailored gene editing. Biomolecules 2021, 11, 181. [Google Scholar] [CrossRef]
  8. Ismail, A.; Riemann, M.; Nick, P. The Jasmonate pathway mediates salt tolerance in grapevines. J. Exp. Bot. 2021, 63, 2127–2139. [Google Scholar] [CrossRef] [Green Version]
  9. Guan, P.; Schmidt, F.; Riemann, M.; Fischer, J.; Thines, E.; Nick, P. Hunting modulators of plant defence: The grapevine trunk disease fungus Eutypa lata secretes an amplifier for plant basal immunity. J. Exp. Bot. 2020, 71, 3710–3724. [Google Scholar] [CrossRef]
  10. Pouzoulet, J.; Pivovaroff, A.L.; Santiago, L.S.; Rolshausen, P.E. Can vessel dimension explain tolerance toward fungal vascular wilt diseases in woody plants? Lessons from Dutch elm disease and esca disease in grapevine. Front. Plant Sci. 2014, 12, 253–264. [Google Scholar] [CrossRef]
  11. Belair, M.; Grau, A.L.; Chong, J.; Tian, X.; Luo, J.; Guan, X.; Pensec, F. Pathogenicity factors of botryosphaeriaceae associated with grapevine trunk diseases: New developments on their action on grapevine defense responses. Pathogens 2022, 11, 951. [Google Scholar] [CrossRef] [PubMed]
  12. Gao, W.; Long, L.; Tian, X.; Jin, J.; Liu, H.; Zhang, H.; Xu, F.; Song, C. Genome-wide identification and expression analysis of stress-associated proteins (SAPs) containing A20/AN1 zinc finger in cotton. Mol. Genet. Genom. 2016, 291, 2199–2213. [Google Scholar] [CrossRef] [PubMed]
  13. Giri, J.; Dansana, P.K.; Kothari, K.S.; Sharma, G.; Vij, S.; Tyagi, A.K. SAPs as novel regulators of abiotic stress response in plants. Bioessays 2013, 35, 639–648. [Google Scholar] [CrossRef] [PubMed]
  14. Dixit, V.M.; Green, S.; Sarma, V.; Holzman, L.B.; Wolf, F.W.; O’Rourke, K.; Ward, P.A.; Prochownik, E.V.; Marks, R.M. Tumor necrosis factor-alpha induction of novel gene products in human endothelial cells including a macrophage-specific chemotaxin. J. Biol. Chem. 1990, 265, 2973–2978. [Google Scholar] [CrossRef] [PubMed]
  15. Opipari, A.W., Jr.; Boguski, M.S.; Dixit, V.M. The A20 cDNA induced by tumor necrosis factor alpha encodes a novel type of zinc finger protein. J. Biol. Chem. 1990, 265, 14705–14708. [Google Scholar] [CrossRef] [PubMed]
  16. Rebagliati, M.R.; Weeks, D.L.; Harvey, R.P.; Melton, D.A. Identification and cloning of localized maternal RNAs from Xenopus eggs. Cell 1985, 42, 769–777. [Google Scholar] [CrossRef] [PubMed]
  17. Huang, J.; Teng, L.; Li, L.; Liu, T.; Li, L.; Chen, D.; Xu, L.G.; Zhai, Z.; Shu, H.B. ZNF216 Is an A20-like and IκB Kinase γ-Interacting Inhibitor of NFκB Activation. J. Biol. Chem. 2004, 279, 16847–16853. [Google Scholar] [CrossRef] [Green Version]
  18. Evans, P.C.; Ovaa, H.; Hamon, M.; Kilshaw, P.J.; Hamm, S.; Bauer, S.; Ploegh, H.L.; Smith, T.S. Zinc-finger protein A20, a regulator of inflammation and cell survival, has de-ubiquitinating activity. Biochem. J. 2004, 378, 727–734. [Google Scholar] [CrossRef] [Green Version]
  19. Mukhopadhyay, A.; Vij, S.; Tyagi, A.K. Overexpression of a zinc-finger protein gene from rice confers tolerance to cold, dehydration, and salt stress in transgenic tobacco. Proc. Natl. Acad. Sci. USA 2004, 101, 6309–6314. [Google Scholar] [CrossRef] [Green Version]
  20. Kang, M.; Fokar, M.; Abdelmageed, H.; Allen, R.D. Arabidopsis SAP5 functions as a positive regulator of stress responses and exhibits E3 ubiquitin ligase activity. Plant Mol. Biol. 2011, 75, 451–466. [Google Scholar] [CrossRef]
  21. Huang, J.; Wang, M.M.; Jiang, Y.; Bao, Y.M.; Huang, X.; Sun, H.; Xu, D.Q.; Lan, H.X.; Zhang, H.S. Expression analysis of rice A20/AN1-type zinc finger genes and characterization of ZFP177 that contributes to temperature stress tolerance. Gene 2008, 420, 135–144. [Google Scholar] [CrossRef] [PubMed]
  22. Giri, J.; Vij, S.; Dansana, P.K.; Tyagi, A.K. Rice A20/AN1 zinc-finger containing stress-associated proteins (SAP1/11) and a receptor-like cytoplasmic kinase (OsRLCK253) Interact via A20 zinc-finger and confer abiotic stress tolerance in transgenic Arabidopsis Plants. New Phytol. 2011, 191, 721–732. [Google Scholar] [CrossRef] [PubMed]
  23. Kanneganti, V.; Gupta, A.K. Overexpression of OsiSAP8, a member of stress associated protein (SAP) gene family of rice confers tolerance to salt, drought and cold stress in transgenic tobacco and rice. Plant Mol. Biol. 2008, 66, 445–462. [Google Scholar] [CrossRef] [PubMed]
  24. Vij, S.; Tyagi, A.K. Genome-wide analysis of the stress associated protein (SAP) gene family containing A20/AN1 zinc-finger(s) in rice and their phylogenetic relationship with Arabidopsis. Mol. Genet. Genom. 2006, 276, 565–575. [Google Scholar] [CrossRef] [PubMed]
  25. Solanke, A.U.; Sharma, M.K.; Tyagi, A.K.; Sharma, A.K. Characterization and phylogenetic analysis of environmental stress-responsive SAP gene family encoding A20/AN1 zinc finger proteins in tomato. Mol. Genet. Genom. 2009, 282, 153–164. [Google Scholar] [CrossRef] [PubMed]
  26. Dong, Q.; Duan, D.; Zhao, S.; Xu, B.; Luo, J.; Wang, Q.; Huang, D.; Liu, C.; Li, C.; Gong, X. Genome-wide analysis and cloning of the apple stress-associated protein gene family reveals MdSAP15, which confers tolerance to drought and osmotic stresses in transgenic Arabidopsis. Int. J. Mol. Sci. 2018, 19, 2478. [Google Scholar] [CrossRef] [Green Version]
  27. Shukla, V.; Choudhary, P.; Rana, S.; Muthamilarasan, M. Structural evolution and function of stress associated proteins in regulating biotic and abiotic stress responses in plants. J. Plant Biochem. Biotechnol. 2021, 30, 779–792. [Google Scholar] [CrossRef]
  28. Tyagi, H.; Jha, S.; Sharma, M.; Giri, J.; Tyagi, A.K. Rice SAPs are responsive to multiple biotic stresses and overexpression of OsSAP1, an A20/AN1 zinc-finger protein, enhances the basal resistance against pathogen infection in tobacco. Plant Sci. 2014, 225, 68–76. [Google Scholar] [CrossRef]
  29. Shu, X.; Ding, L.; Gu, B.; Zhang, H.; Guan, P.; Zhang, J. A stress associated protein from Chinese wild Vitis amurensis, VaSAP15, enhances the cold tolerance of transgenic grapes. Sci. Hortic. 2021, 285, 110147–110157. [Google Scholar] [CrossRef]
  30. Fung, R.W.; Gonzalo, M.; Fekete, C.; Kovacs, L.G.; He, Y.; Marsh, E.; Mcintyre, L.M.; Schachtman, D.P.; Qiu, W. Powdery mildew induces defense-oriented reprogramming of the transcriptome in a susceptible but not in a resistant grapevine. Plant Physiol. 2008, 146, 236–249. [Google Scholar] [CrossRef]
  31. Cramer, G.R.; Ergul, A.; Grimplet, J.; Tillett, R.L.; Tattersall, E.A.; Bohlman, M.C.; Vincent, D.; Sonderegger, J.; Evans, J.; Osborne, C.; et al. Water and salinity stress in grapevines: Early and late changes in transcript and metabolite profiles. Funct. Integr. Genom. 2007, 7, 111–134. [Google Scholar] [CrossRef] [PubMed]
  32. Tattersall, E.A.R.; Grimplet, J.; Deluc, L.; Wheatley, M.D.; Vincent, D.; Osborne, C.; Erguel, A.; Lomen, E.; Blank, R.R.; Schlauch, K.A.; et al. Transcript abundance profiles reveal larger and complex responses of grapevine to chilling compared to osmotic and salinity stress. Funct. Integ. Genom. 2007, 7, 317–333. [Google Scholar] [CrossRef] [PubMed]
  33. Zhao, L.; You, S.; Zou, H.; Guan, X. Transcriptome analysis and cell morphology of Vitis rupestris cells to Botryosphaeria dieback pathogen Diplodia seriata. Genes 2021, 12, 179. [Google Scholar] [CrossRef] [PubMed]
  34. Guan, X.; Buchholz, G.; Nick, P. Tubulin marker line of grapevine suspension cells as a tool to follow early stress responses. J. Plant Physiol. 2015, 176, 118–128. [Google Scholar] [CrossRef] [PubMed]
  35. 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–219. [Google Scholar] [CrossRef]
  36. Zarattini, M.; Farjad, M.; Launay, A.; Cannella, D.; Soulie, M.C.; Bernacchia, G.; Fagard, M. Every cloud has a silver lining: How abiotic stresses affect gene expression in plant-pathogen interactions. J. Exp. Bot. 2021, 72, 1020–1033. [Google Scholar] [CrossRef]
  37. Kissoudis, C.; Clemens, V.; Visser, R.; Gerard, V. Enhancing crop resilience to combined abiotic and biotic stress through the dissection of physiological and molecular crosstalk. Front. Plant Sci. 2014, 5, 207–226. [Google Scholar] [CrossRef] [Green Version]
  38. Vij, S.; Tyagi, A.K. A20/AN1 zinc-finger domain-containing proteins in plants and animals represent common elements in stress response. Funct Integr Genom. 2008, 8, 301–307. [Google Scholar] [CrossRef]
  39. Zhou, Y.; Zeng, L.; Chen, R.; Wang, Y.; Song, J. Genome-wide identification and characterization of stress-associated protein (SAP) gene family encoding A20/AN1 zinc-finger proteins in Medicago truncatula. Arch. Biol. Sci. 2018, 70, 87–98. [Google Scholar] [CrossRef]
  40. Hozain, M.; Abdelmageed, H.; Lee, J.; Kang, M.; Fokar, M.; Allen, R.D.; Holaday, A.S. Expression of AtSAP5 in cotton up-regulates putative stress-responsive genes and improves the tolerance to rapidly developing water deficit and moderate heat stress. J. Plant Physiol. 2012, 169, 1261–1270. [Google Scholar] [CrossRef]
  41. Liu, S.; Yuan, X.; Wang, Y.; Wang, H.; Wang, J.; Shen, Z.; Gao, Y.; Cai, J.; Li, D.; Song, F. Tomato stress-associated protein 4 contributes positively to immunity against necrotrophic fungus Botrytis cinerea. Mol. Plant-Microbe Interact. 2018, 32, 566–582. [Google Scholar] [CrossRef] [PubMed]
  42. Lai, W.; Zhou, Y.; Pan, R.; Liao, L.; He, J.; Liu, H.; Yang, Y.; Liu, S. Identification and expression analysis of stress-associated proteins (SAPs) containing A20/AN1 zinc finger in cucumber. Plants 2020, 9, 400. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Liu, Y.; Wang, K.; Cheng, Q.; Kong, D.; Zhang, X.; Wang, Z.; Wang, Q.; Xie, Q.; Yan, J.; Chu, J.; et al. Cysteine protease RD21A regulated by E3 ligase SINAT4 is required for drought-induced resistance to Pseudomonas syringae in Arabidopsis. J. Exp. Bot. 2020, 71, 5562–5576. [Google Scholar] [CrossRef] [PubMed]
  44. Guan, X.; Buchholz, G.; Nick, P. Actin marker lines in grapevine reveal a gatekeeper function of guard cells. J. Plant Physiol. 2014, 171, 1164–1173. [Google Scholar] [CrossRef] [PubMed]
  45. Paolinelli-Alfonso, M.; Villalobos-Escobedo, J.M.; Rolshausen, P.; Herrera-Estrella, A.; Galindo-Sanchez, C.; Lopez-Hernandez, J.F.; Hernandez-Martinez, R. Global transcriptional analysis suggests Lasiodiplodia theobromae pathogenicity factors involved in modulation of grapevine defensive response. BMC Genom. 2016, 17, 615–634. [Google Scholar] [CrossRef] [Green Version]
  46. Sun, Q.; Rost, T.L.; Reid, M.S.; Matthews, M.A. Ethylene and not embolism is required for wound-induced tylose development in stems of grapevines. Plant Physiol. 2007, 145, 1629–1636. [Google Scholar] [CrossRef] [Green Version]
  47. Salvatore, M.; Giambra, S.; Naviglio, D.; DellaGreca, M.; Salvatore, F.; Burruano, S.; Andolfi, A. Fatty acids produced by Neofusicoccum vitifusiforme and N. parvum, fungi associated with grapevine Botryosphaeria dieback. Agriculture 2018, 8, 189. [Google Scholar] [CrossRef] [Green Version]
  48. Wang, L.X.; Sadeghnezhad, E.; Riemann, M.; Nick, P. Microtubule dynamics modulate sensing during cold acclimation in grapevine suspension cells. Plant Sci. 2019, 280, 18–30. [Google Scholar] [CrossRef] [PubMed]
  49. Guan, X.; Buchholz, G.; Nick, P. The cytoskeleton is disrupted by the bacterial effector HrpZ, but not by the bacterial PAMP flg22, in tobacco BY-2 cells. J. Exp. Bot. 2013, 64, 1805–1816. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Shigo, A.L.; Marx, H.G. Compartmentalization of Decay in Trees; US Department of Agriculture, Forest Service: Washington, DC, USA, 1977; Volume 405, pp. 1–73. Available online: https://www.fs.usda.gov/treesearch/pubs/5292 (accessed on 22 June 2022). [CrossRef]
  51. Gus-Mayer, S.; Naton, B.; Hahlbrock, K.; Schmelzer, E. Local mechanical stimulation induces components of the pathogen defense response in parsley. Proc. Natl. Acad. Sci. USA 1998, 95, 8398–8403. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Durst, S.; Hedde, P.N.; Brochhausen, L.; Nick, P.; Nienhaus, G.U.; Maisch, J. Organization of perinuclear actin in live tobacco cells observed by PALM with optical sectioning. J. Plant Physiol. 2014, 171, 97–108. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Conserved motif arrangement in grapevine SAP proteins according to the phylogenetic relationship. (A) Conserved motif distributions of the VvSAP proteins annotated with the MEME server. Ten motifs are marked by different colours; (B) analysis of conserved domains of VvSAP proteins; (C) analysis of cis-regulatory elements of VvSAP protein; (D) analysis of VvSAP protein coding.
Figure 1. Conserved motif arrangement in grapevine SAP proteins according to the phylogenetic relationship. (A) Conserved motif distributions of the VvSAP proteins annotated with the MEME server. Ten motifs are marked by different colours; (B) analysis of conserved domains of VvSAP proteins; (C) analysis of cis-regulatory elements of VvSAP protein; (D) analysis of VvSAP protein coding.
Pathogens 11 01426 g001
Figure 2. Distribution of the VvSAP genes among eight grapevine chromosomes. The scale on the left represents the length of each chromosome displayed in megabase (Mb). Tandemly duplicated genes are coloured with red.
Figure 2. Distribution of the VvSAP genes among eight grapevine chromosomes. The scale on the left represents the length of each chromosome displayed in megabase (Mb). Tandemly duplicated genes are coloured with red.
Pathogens 11 01426 g002
Figure 3. Phylogenetic analysis of SAP proteins from grapevine (V. vinifera) and other two plants (A. thaliana and O. sativa). Sequence alignment was performed with full-length SAP amino acid sequences from different plant species using MAFFT as the default parameters, and the neighbour-joining (NJ) phylogenetic tree was created by MEGA 7.0 using a bootstrap option of 1000 replications.
Figure 3. Phylogenetic analysis of SAP proteins from grapevine (V. vinifera) and other two plants (A. thaliana and O. sativa). Sequence alignment was performed with full-length SAP amino acid sequences from different plant species using MAFFT as the default parameters, and the neighbour-joining (NJ) phylogenetic tree was created by MEGA 7.0 using a bootstrap option of 1000 replications.
Pathogens 11 01426 g003
Figure 4. Alignment of conserved domains of SAP family in V. vinifera. Alignment of conserved A20/AN1 and grapevine SAP protein sequences. Identical proteins are highlighted in red boxes.
Figure 4. Alignment of conserved domains of SAP family in V. vinifera. Alignment of conserved A20/AN1 and grapevine SAP protein sequences. Identical proteins are highlighted in red boxes.
Pathogens 11 01426 g004
Figure 5. Expression levels of VvSAP genes in different tissues of grapevine. Abbreviations stand for plant developmental stages, where FS stands for fruit set, PFS for post-fruit set, V for veraison, MR for mid-ripening, R for ripening, PHWI for post-harvest withering I (1st month), BerryPericarp-PHWII for post-harvest withering II (2nd month) and BerryPericarp-PHWIII for post-harvest withering III (3rd month), Y stands for Young, FB stands for flowering begins, F stands for flowering, ML stands for mature leaf, WD stands for well developed, W stands for woody, S represents different developmental stages in Bud and Leaf, swell in bud and senescence in leaf. Black lines indicate reproductive and vegetative organs, respectively. The colour scale represents log 2 expression values.
Figure 5. Expression levels of VvSAP genes in different tissues of grapevine. Abbreviations stand for plant developmental stages, where FS stands for fruit set, PFS for post-fruit set, V for veraison, MR for mid-ripening, R for ripening, PHWI for post-harvest withering I (1st month), BerryPericarp-PHWII for post-harvest withering II (2nd month) and BerryPericarp-PHWIII for post-harvest withering III (3rd month), Y stands for Young, FB stands for flowering begins, F stands for flowering, ML stands for mature leaf, WD stands for well developed, W stands for woody, S represents different developmental stages in Bud and Leaf, swell in bud and senescence in leaf. Black lines indicate reproductive and vegetative organs, respectively. The colour scale represents log 2 expression values.
Pathogens 11 01426 g005
Figure 6. Expression of VvSAPs under biotic stresses. VvSAPs of V. rupestris expressing GFP-AtTUB6 under the infection of culture filtration of N. parvum Bt-67 and B. dothidea Sd-18 at 12 hpi (A); N. parvum Bt-67 and D. seriata 98.1 at 3 dpi (B); and VvSAPs of V. aestivalis × V. vinifera ‘Norton’ and V. vinifera ‘Cabernet Sauvignon’ under infection of conidiospores of Erysiphe necator (C). cf, culture filtration, cs conidiospores, con, control. The expression level of genes was determined based on the value of FPKM. The colour scale represents log 2 expression values.
Figure 6. Expression of VvSAPs under biotic stresses. VvSAPs of V. rupestris expressing GFP-AtTUB6 under the infection of culture filtration of N. parvum Bt-67 and B. dothidea Sd-18 at 12 hpi (A); N. parvum Bt-67 and D. seriata 98.1 at 3 dpi (B); and VvSAPs of V. aestivalis × V. vinifera ‘Norton’ and V. vinifera ‘Cabernet Sauvignon’ under infection of conidiospores of Erysiphe necator (C). cf, culture filtration, cs conidiospores, con, control. The expression level of genes was determined based on the value of FPKM. The colour scale represents log 2 expression values.
Pathogens 11 01426 g006
Figure 7. Expression of VvSAPs under abiotic stresses. Expression analysis of VvSAPs under drought, salt (match water potential to those of the water-deficit-treated plants), and cold (5 ℃) stresses in grapevine in short-term treatment (A); and long-term treatment (B). The expression level of genes was determined based on the value of RPKM. The colour scale represents log 2 expression values.
Figure 7. Expression of VvSAPs under abiotic stresses. Expression analysis of VvSAPs under drought, salt (match water potential to those of the water-deficit-treated plants), and cold (5 ℃) stresses in grapevine in short-term treatment (A); and long-term treatment (B). The expression level of genes was determined based on the value of RPKM. The colour scale represents log 2 expression values.
Pathogens 11 01426 g007
Table 1. Identification and characterisation of SAP family genes in grapevine.
Table 1. Identification and characterisation of SAP family genes in grapevine.
Gene Name.Accession NumberProteinMWpIGRAVYAliphatic IndexLocChrom
VvSAP1VIT_01s0011
g02290.t01
15417280.028.78−0.66152.6NucleusChr1
VvSAP2VIT_02s0025
g05050.t01
17218460.838.77−0.49446.05NucleusChr2
VvSAP3VIT_02s0025
g05080.t01
10010918.548.98−0.49852.7NucleusChr2
VvSAP4VIT_06s0004
g01790.t01
14216150.559.37−0.84957.68NucleusChr6
VvSAP5VIT_06s0004
g01820.t01
17218502.957.99−0.556.8NucleusChr6
VvSAP6VIT_08s0007
g03530.t01
16117192.688.92−0.366.21NucleusChr8
VvSAP7VIT_08s0007
g07950.t01
17218153.577.99−0.32768.72NucleusChr8
VvSAP8VIT_08s0007
g08360.t01
29332457.788.56−0.61458.84NucleusChr8
VvSAP9VIT_13s0064
g01210.t01
17218504.118.22−0.37161.4NucleusChr13
VvSAP10VIT_13s0064
g01220.t01
17218359.698.95−0.48657.5NucleusChr13
VvSAP11VIT_14s0066
g01880.t01
18921133.118.74−0.847.99NucleusChr14
VvSAP12VIT_16s0022
g01680.t01
23926461.398.98−0.77251.46Chloro-plast NucleusChr16
VvSAP13VIT_16s0022
g01980.t01
15216850.548.77−0.55858.36NucleusChr16
VvSAP14VIT_18s0001
g00430.t01
677980.510.26−0.67659.7NucleusChr18
VvSAP15VIT_18s0001
g01260.t01
677966.4310.2−0.6759.7NucleusChr18
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sun, X.; Xia, X.; Guan, X. Genome-Wide Identification and Characterisation of Stress-Associated Protein Gene Family to Biotic and Abiotic Stresses of Grapevine. Pathogens 2022, 11, 1426. https://doi.org/10.3390/pathogens11121426

AMA Style

Sun X, Xia X, Guan X. Genome-Wide Identification and Characterisation of Stress-Associated Protein Gene Family to Biotic and Abiotic Stresses of Grapevine. Pathogens. 2022; 11(12):1426. https://doi.org/10.3390/pathogens11121426

Chicago/Turabian Style

Sun, Xiaoye, Xue Xia, and Xin Guan. 2022. "Genome-Wide Identification and Characterisation of Stress-Associated Protein Gene Family to Biotic and Abiotic Stresses of Grapevine" Pathogens 11, no. 12: 1426. https://doi.org/10.3390/pathogens11121426

APA Style

Sun, X., Xia, X., & Guan, X. (2022). Genome-Wide Identification and Characterisation of Stress-Associated Protein Gene Family to Biotic and Abiotic Stresses of Grapevine. Pathogens, 11(12), 1426. https://doi.org/10.3390/pathogens11121426

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