Skip to main content
  • Research article
  • Open access
  • Published:

Identification of C3H2C3-type RING E3 ubiquitin ligase in grapevine and characterization of drought resistance function of VyRCHC114

Abstract

Background

RING is one of the largest E3 ubiquitin ligase families and C3H2C3 type is the largest subfamily of RING, which plays an important role in plant growth and development, and growth and responses to biotic and abiotic stresses.

Results

A total of 143 RING C3H2C3-type genes (RCHCs) were discovered from the grapevine genome and separated into groups (I-XI) according to their phylogenetic analysis, and these genes named according to their positions on chromosomes. Gene replication analysis showed that tandem duplications play a predominant role in the expansion of VvRCHCs family together. Structural analysis showed that most VvRCHCs (67.13 %) had no more than 2 introns, while genes clustered together based on phylogenetic trees had similar motifs and evolutionarily conserved structures. Cis-acting element analysis showed the diversity of VvRCHCs regulation. The expression profiles of eight DEGs in RNA-Seq after drought stress were like the results of qRT-PCR analysis. In vitro ubiquitin experiment showed that VyRCHC114 had E3 ubiquitin ligase activity, overexpression of VyRCHC114 in Arabidopsis improved drought tolerance. Moreover, the transgenic plant survival rate increased by 30 %, accompanied by electrolyte leakage, chlorophyll content and the activities of SOD, POD, APX and CAT were changed. The quantitative expression of AtCOR15a, AtRD29A, AtERD15 and AtP5CS1 showed that they participated in the response to drought stress may be regulated by the expression of VyRCHC114.

Conclusions

This study provides valuable new information for the evolution of grapevine RCHCs and its relevance for studying the functional characteristics of grapevine VyRCHC114 genes under drought stress.

Background

To survive in a changing environment, post-translational modification of proteins often occurs when plants perceive and transmit internal or external signals. The acetylation, methylation, phosphorylation, and ubiquitination of proteins are the main types of post-translational modification, which play a key role in different plant development stages and plant-environment interactions. The process of classifying intracellular proteins under the action of a variety of special enzymes, and specifically modifying the screened target proteins, is called ubiquitination [1]. In eukaryotic cells, the ubiquitin sfigystem is complex and mainly involving ubiquitin (a small molecule protein), intact 26 S proteasome, ubiquitin-activating enzyme (E1), ubiquitin-binding enzyme (E2), and ubiquitin-ligase (E3) [2]. The inactivated ubiquitin-dependent ATP is first activated by E1 through the thioester bond formed between the C-terminal of ubiquitin and the cysteine residue of E1; then the ubiquitin signal connected to E1 is transferred to the acetylcysteine of E2. In the next step, the ubiquitin linked to E2 is transferred directly or indirectly to the lysine residue of the target protein via E3. It is noteworthy that E3 ubiquitin ligase is the main factor determining the binding of specific protein during ubiquitination process [3], it can repeatedly add ubiquitin to the substrate protein, so that the target protein is degraded by the 26 S protease [4].

E3 ubiquitin ligases can be divided into 9 categories based on specific conserved domains: RING, HECT, U-box, F-box, cullin, BTB, DDB, RBX and SKP. RING E3 ligases protein has a conserved RING domain, which can provide residence sites for E2 and specific substrates and enable E2-bound ubiquitin molecules to transfer to the host protein, thus completing the ubiquitination process. In the RING domain, there are eight conserved amino acids (Cys or His) located in the center of the three-dimensional protein structure, which can combine with two zinc ions to help stabilize the entire structure. According to the types of conserved amino acid sites, they are divided into different subfamily. Among them, RING C3H2C3 is the largest subfamily. The RING conserved domain sequence of this family member is Cys-X2-Cys-X(9–39)-Cys-X(1–3)-His-X(2–3)-His-X2-Cys-X(4–48)-Cys-X2-Cys, X is any amino acid. In recent years, a growing number of studies have shown that the RING E3 ligase gene also figure prominently in abiotic stress responses of plants [5]. SpRing is a RING-type E3 ubiquitin ligase located in endoplasmic reticulum and participates in salt stress signal transmission in wild tomato variety Solanum pimpinellifolium ‘PI365967’. In addition, SpRing is silenced by virus-induced gene silencing, resulting in increased sensitivity of wild tomato to salt stress. Overexpression of SpRing in Arabidopsis can improve its salt tolerance [6]. SDIR1 (SALT AND DROUGHT-INDUCED REALLY INTERESTING NEW GENE FINGER1) is a RING-type E3 ubiquitin ligase that regulates the salt stress response and ABA signaling in Arabidopsis by degrading the target protein SDIRIP1 (SDIR1-INTERACTING PROTEIN1). The downstream transcription factor ABI5 (ABA-INSENSITIVE5) is regulated by SDIRIP1, and overexpression of ABI5 increases salt tolerance [7]. The E3 ubiquitin ligase OsHTAS (Oryza sativa HEAT TOLERANCE AT SEEDLING STAGE) regulates the stomatal opening state of leaves by regulating ROS homeostasis, thus improving the basal heat resistance of leaves. It involves two pathways, ABA-mediated [8]. In Arabidopsis, CHYR1 (CHYZINC-FINGERANDRINGPROTEIN1) encodes the RING-type E3 ubiquitin ligase which interacts with the related protein kinase KINASE2 (SnRK2) and can be phosphorylated by SnRK2.6 on its Thr-178 residues. When mediated by ABA, CHYR1 promotes the production of reactive oxygen species (ROS), stomatal closure, and drought tolerance in plants [9]. The capsicum annular E3 ubiquitin ligase, CaAIRF1 (Capsicum annuum ADIP1 INTERACTING RING FINGER PROTEIN1), can interact with protein phosphatase CaADIP1 and positively regulate ABA signaling pathway to improve drought tolerance [10]. In Zea mays, ZmXerico1 encodes a RING-type E3 ligase, which can regulate the stability of ABA8’-hydroxylase protein and thereby enable control of the dynamic balance of ABA, hence, expression of ZmXerico1 endows maize plants with ABA sensitivity and improves their water use efficiency under drought stress [11]. Furthermore, Arabidopsis AtAIRP1, AtAIRP2, AtAIRP3 and CaAIR1 jointly encode an E3 ubiquitin ligase, which regulates drought responses by regulating ABA signaling transduction, the expression of these genes increases ABA-mediated stomatal closure [12,13,14,15]. Collectively, the above studies show that E3 ligase plays a crucial role in response to abiotic stress.

Grapevine (Vitis vinifera L.) is a major cash crop, whose cultivated varieties have a total worldwide output of nearly 70 million tons of the fruit berries from more 7 million hectares of harvested land [16]. This plant is mainly grown to produce table grapes, fruit juices, and wine [17]. Most grapevine producing areas in the world incur seasonal droughts. According to global climate modeling, droughts will intensify in the near future. Drought can adversely affect the growth and development of grapevines, because under drought stress the concentration of cytokinin in grape stems decreases, vegetative and reproductive growth are inhibited [18]. When grapevines are in full bloom, drought stress will also affect the pollination process, which decreases the fruit setting rate and affects the size of the individual fruit berries produced [19]. With worsening water shortages, drought stress is likely to become a key factor impacting grapevine and wine production worldwide [20]. Therefore, it is of great significance to grapevine production and breeding to study the drought resistance of wild grapevine plants as this could uncover the molecular mechanisms enabling them to withstand drought effects. Vitis yeshanensis is a wild grapevine plant native to arid areas of China, whose morphological characteristics indicate adaptability to arid environments in many aspects [21]. Several studies have shown that wild Vitis yeshanensis has stronger drought resistance than other cultivars [22, 23].

The RING-type gene family has been found in more and more plant species, and its importance in plant stress responses and growth and development has been recognized, but RING-type genes have not been fully identified in grapevine. It is reported that the RING type E3 ubiquitin ligase is involved in grapevine stress and growth, but few studies have investigated the involvement of E3 ubiquitin ligase in regulating grapevine response to drought stress. We assume that RCHC protein may mediate the ubiquitination of key factors during grape drought stress to regulate plant drought resistance. This study aimed to characterize the RING-type E3 ubiquitin ligase in grapevine’s genome and its relevance for drought stress. Genome-wide identification of C3H2C3 genes, the largest subfamily of grapevine RING-type, was carried out, coupled to their phylogenetic analysis, gene structure analysis, chromosome mapping, gene replication analysis, and cis-acting element analysis in gene promoter regions. We also quantified the expression levels of these genes under simulated drought treatment, 136 RCHC genes were found to be expressed and 52 DEGs, which 8 DEGs at least 3 stages. The VyRCHC114 gene was confirmed by RT-qPCR, and then the ubiquitin ligase activity of the gene was verified. The function of the gene under drought conditions was elucidated using Arabidopsis transgenic plants. Our study provides an important basis for the involvement of RCHC protein in the regulation of grape ubiquitination under drought stress.

Results

Genome-wide identification of RING C3H2C3 type finger proteins in grapevine

The results of the Hidden Markov Model (HMM) were analyzed, and the gene sequences were extracted and given to SMART, CDD, and Pfam for domain authentication. From this, 143 VvRCHC genes were obtained by comparing and screening genes with eight conservative metal ligands, and the alignment members were not abandoned. The physicochemical properties of 143 VvRCHCs were identified (Table 1). The number of amino acids encoded by the 143 VvRCHCs ranged from 70 (VvRCHC50) to 763 (VvRCHC98). For these genes, the molecular weights of their products varied from 7.83 kDa to 83.58 kDa, while their isoelectric points varied from 3.88 to 9.95.

Table 1 Detailed information of all 143 VvRCHCs identified in grapevine genome

Analysis of VvRCHCs in the C3H2C3 domain

The typical RING domain is considered to be an octahedral group of metal-bound cysteine and its residues, which can chelate two zinc ions in a spherical cross-supported structure, in which the metal ligands 1 and 3, 2 and 4, each bind to one zinc ion. This structure requires a certain distance between adjacent metal ligands, it being variable between ml2 ~ ml3 and ml6 ~ ml7. We calculated statistics for this distance between adjacent metal ligands (Table S2). It was found that, except those between ml2 ~ ml3 and ml6 ~ ml7, the distances between other metal ligands were constant, while those from ml2 to ml3 spanned 11 to 24 amino acids, and for ml6 ~ ml7 the distance varied from 8 to 14 amino acids. The 143 VvRCHCs C3H2C3 domains have two amino acids between ml1 ~ ml2 and ml5 ~ ml6, while ml3 ~ ml4 contains one amino acid, ml7 ~ ml8 contains two amino acids as does ml4 ~ ml5. To understand whether these RING C3H2C3 structural domains are conserved apart from their eight special metal ligands, their comparative analysis was conducted (Fig. S1). This revealed that some amino acids in the structural domain of RING C3H2C3 have a typical position bias (Fig. 1a). In the C3H2C3 type RING region, the ml2 located in front of amino acid residues is the most common Ile (I) or Val (V); likewise, the phenylalanine (Phe, F) residue is typically before ml5, the leucine residue (Leu, L) is always next to ml2, and the aspartic acid (Asp, D) residue is usually positioned after ml6, while the tryptophan residue (Trp, W) is usually the fourth following ml6. Notably, a very conservative proline (P) was found situated after ml7. According to the RING-type C3H2H3 domain schematic diagram, two pairs of metal ligands bind to a zinc ion (Fig. 1b). The total amino acid length of the C3H2C3 domain per VvRCHC gene and the corresponding number of different lengths were calculated: the vast majority of these were 41 and 42, accounting for 88.8 % of all genes (Fig. 1c).

Fig. 1
figure 1

VvRCHCs C3H2H3 domain analysis. a Sequence of C3H2H3 domain. b Schematic diagram of zinc ion binding to the C3H2C3 structure of the VvRCHCs. Cysteine and histidine are metal ligands marked numerically in ellipses. The hexagon represents a zinc ion. The number on the line between metal ligands is the number of amino acids between them. c The number of different amino acids in the VvRCHCs in the metal ligand region

Phylogenetic analysis of VvRCHCs

To infer the evolutionary relationships of grapevine VvRCHCs, phylogenetic analysis of RCHC protein sequences of Arabidopsis, tomato, and grapevine were constructed (using the Maximum Likelihood method). According to the phylogenetic analysis, these 180 genes can be divided into 6 subgroups: I ~ VI (Fig. 2). Group I has the least number of members, only 12, and the group of the largest number of members is group III, while the RCHC gene of Arabidopsis thaliana or tomato is found in each group. It is worth noting that more RCHC genes of Arabidopsis thaliana and tomato are gathered in group VI. Most of the RING-type C3H2C3 genes of grapevine display some homology to RCHC genes of Arabidopsis or tomato. In addition, some gene pairs showed high similarity in different groups, which were confirmed in the distance of evolutionary relationship, the location of RING conserved domain and the length of protein sequence. For instance, SlATL33 and VvRCHC62, SlATL46 and VvRCHC108, SlATL51 and VvRCHC110, AtBRH1 and VvRCHC116, AtRHA1A and VvRCHC13, AtSDIR1 and VvRCHC97, AtRHC1A and VvRCHC59 etc. Next, a phylogenetic tree containing only 143 VvRCHC protein sequences was constructed. To facilitate their study and analysis, the 143 members were divided into 6 groups (I ~ VI) according to the classification and phylogenetic analysis of Fig. 3a, from which 27 pairs of genes with high homology were found. Based on their color-coded names, the VvRCHCs were then divided into six groups according to the number of conserved amino acids in their protein sequence.

Fig. 2
figure 2

Phylogenetic analysis of RCHC protein in grapevine and other plants. Via the ML method, MEGA7.0 was used to construct phylogenetic trees of grapevine, Arabidopsis, and tomato (Solanum lycopersicum) with RING-type C3H2H3 proteins. The number of bootstrap repeats was n = 1000. Displayed are the percentages of bootstrap scores greater than 50 %. Conservative domains, grouping, and information on different species are shown in the form of a, b, c

Fig. 3
figure 3

Phylogenetic tree, gene domain, and structure analysis of VvRCHCs in grapevine. a The phylogenetic tree of VvRCHCs was constructed using the ML method. Different background colors represent different grouping branches. b Domain analysis of VvRCHCs proteins. At the bottom of the line, different colored squares represent different types of conserved amino acid sequences and based on MEME analysis. The modules of different colors above the line represent the functional domains that have been identified. c Genetic structure of VvRCHCs, the CDS sequence is represented by a blue square/rectangle, the introns by black lines

Characterization of the motifs and gene structure of VvRCHCs

To further understand the diversity in motif composition between VvRCHCs, the MEME analysis of VvRCHC proteins from groups I to VI was carried out. From this, 12 conserved motifs were identified in the VvRCHC protein, named motif 1 to motif 12 (Fig. 3b), in which motif 1 and motif 2 is found in almost every VvRCHCs, this motif combines to form the eight most important metal ligand (Cys-Cys-Cys-His-His-Cys-Cys-Cys) structures of every VvRCHC gene. Importantly, there are 13 such structures in some genes, such as PA, CUE, DUF1117, zinc_ribbon_9, and zf-CHY, among others. These structures domain could be relevant for the function of VvRCHCs. The sequence information of motif 1 ~ 12 is presented in Table 2; Fig. 3d (motif data). We next analyzed the exons, introns, and several key structures of VvRCHCs (Fig. 3c). Most VvRCHCs (67.13 %) had no more than 2 introns, with a maximum of 19 introns in VvRCHC29 and none intron in 57 VvRCHCs (Fig. S3). The longest intron length was found in VvRCHC141.

Table 2 Motif data information in MEME analysis of VvRCHCs

According to the phylogenetic analysis of VvRCHCs (Fig. 3a), 45 pairs of genes can be found in the evolutionary tree. The results of the MEME and gene structure analyses of these gene pairs were also similar (Fig. 3b and c). For example, the conserved motifs in the protein sequences of VvRCHC44/64 are highly similar, and the structure type and length are also similar, such as for VvRCHC94/96, VvRCHC38/97, VvRCHC18/78, VvRCHC28/67 and VvRCHC11/107, to name a few. Unexpectedly, the MEME analysis of VvRCHC55/127, VvRCHC105/133, and VvRCHC13/116 gene pairs gave near identical results to those from the gene structure analysis, revealing a remarkably similar protein sequence length, gene structure length and the intron number among them. We thus speculate these four gene pairs may perform similar functions in grapevine plants.

Chromosomal localization and gene replication analysis of VvRCHCs

According to the location of VvRCHCs in the grapevine genome, 143 VvRCHCs were placed on 20 chromosomes (Fig. 4a), albeit unevenly distributed among them. Imprinting of the VvRCHCs was found in each chromosome of grapevine, but the number of genes on different chromosomes varied. The most found were 12 VvRCHCs on chromosome 11, the 11 VvRCHCs were identified on chromosome 1,7,13 and 18. Further, we also observed that these most of these VvRCHCs are likely distributed at both ends of the chromosome, leaving only a small portion of them in its middle part. Gene replication events include tandem replication and segmental replication, which are very vital for expanding the number of members of the gene family. To clarify the amplification mechanism of VvRCHCs during their evolution, we studied their potential repetitive events of VvRCHCs. According to the intraspecific alignment of 143 VvRCHCs, 9 pairs of genes, 7 and 2, were respectively identified as associated with tandem or segmental replication events. Among the 9 pairs of gene events, the tandem repeat frequency between chromosomes 1 was the highest, there were six tandem replication events, moreover, one pair of genes on chromosomes 3 identified as tandem replication genes. These results suggested that the main replication event mode of grapevine VvRCHCs family is via tandem replication; hence, it could have played a crucial role in the amplification of VvRCHCs during their evolutionary history.

Fig. 4
figure 4

Chromosome location, gene replication, and collinear correlation. a Location information of VvRCHCs on grapevine chromosomes. Colored boxes and line connections represent the segmental repetitive gene pairs. b Collinear correlations of RCHCs in grapevine, tomato, and Arabidopsis. Each colored square is a chromosome (serial number). The black-and-white square shows the density of genes in each chromosome. Each line connects homologous genes; colored lines correspond to RCHCs, while gray lines denote other genes

To explore the selection of grapevine VvRCHCs in terms of their repetition and differentiation, the non-synonymous (Ka), synonymous (Ks), and Ka/Ks of each duplicated VvRCHCs were calculated. Among the 9 pairs of repetitive genes in grapevine, the Ka/Ks values of one pair were less than 0.5, while the average Ka/Ks value was 0.325. It is worth noting that 8 pairs had Ka/Ks values less than 0.5, indicating that most of the repeated grapevine VvRCHCs were under negative selection during evolution (Table 3). Figure 4b shows that grapevine, Arabidopsis, and tomato all retained similar RCHC genes in their evolutionary history. It is worth noting the absence of homologous genes with VvRCHC29 in tomato, but their presence in Arabidopsis, which may have arisen from gene deletions in the process of evolution, given that the same genes are VvRCHC11, VvRCHC38, VvRCHC107, VvRCHC119, and VvRCHC137. Nonetheless, two or more RCHC genes in Arabidopsis and tomato were found homologous to one VvRCHC gene; for example, VvRCHC89 and Solyc07g053850.3/Solyc12g005470.2 and AT4G28370/AT2G20650, as well as those of VvRCHC1, VvRCHC32, VvRCHC97, VvRCHC104, VvRCHC118, and VvRCHC142. Hence, these genes may be parallel gene pairs and the putative source of amplifications of RCHC genes during evolution.

Table 3 Ka/Ks analysis and divergence time estimated for grapvine duplicated VvRCHCs paralogs

Cis-acting element analysis in VvRCHCs promoter

To further investigate the transcriptional regulation of VvRCHCs, cis-acting elements in the 2000 bp region upstream of the VvRCHCs’ codon was predicted. The predicted cis-acting elements can be divided into seven categories according to their functions: namely, light response (32), hormone response (11), growth and development response (9), stress response (6), enhanced promoter cis-acting (6), binding site cis-acting (6), and other functional cis-acting (2) elements. Most promoters of grapevine VvRCHCs contained the CAAT-box or TATA-box, which are involved in the enhanced promoter cis-acting elements. In addition, 127 VvRCHCs promoters harbored the stress response element ARE, more than half of the promoters of the VvRCHCs having the hormone response elements ABRE, TGACG-motif, CGTCA-motif, and over half of the VvRCHCs also featured the G-box, GT1-motif, and Box 4 in their promoters (Table S3). In the 2000 kb region upstream of VvRCHCs, discovered many different functions of cis element, in addition to the common cis element with light response and enhanced the promoter, also found that the more growth and adversity stress related cis element, this suggests that VvRCHCs may be widely participating in various life activities of plant.

It is known that the RING gene play a key role in plant growth and response to abiotic stresses. Accordingly, the cis-acting elements related to abiotic stress, growth and hormone regulation were focused upon here. The respective locations of the five major acting elements associated with hormone response, binding sites, growth and development, and stress of our concern, on the promoter of the VvRCHCs (Fig. S4a) were determined. To accurately identify the stress-related elements, we focused on four kinds (anaerobic induction, injury response, low temperature, drought response (Fig. S4b), low temperature response, defense and stress response), whose locations are also depicted. In addition, we counted the number of major elements related to stress, growth and development, and hormone responses in the VvRCHC gene promoter (Fig. S4b). Evidently, concerning growth and development, the number of O2-sites is the largest, there are 5 promoters of VvRCHC6 and 4 promoters of VvRCHC40. In terms of stress, the number of ARE is very large, found in 89 % of the VvRCHCs promoters, moreover, 5 of the most promoters of VvRCHC14 and VvRCHC81 occurred. In terms of hormone response, the number of ABRE is dominant, found in 64 % of the VvRCHCs promoters, moreover, 9 of the most promoters of VvRCHC3 and VvRCHC16 occurred. Surprisingly, 22 of the VvRCHC74 gene promoters were found and 11 of the VvRCHC128 gene promoters were found. These results suggest that VvRCHCs may be associated with cis-acting elements of different functions; in other words, these genes may be regulated by these elements and thereby influence related plant life activities.

Expression analysis of VvRCHCs in roots of two grapevine rootstocks with different drought sensitivity

To investigate differential VvRCHCs’ expression between plants having contrasting drought-resistant genes (101.14 vs. M4) under drought stress and their potential functioning, the grapevine RNA-Seq transcriptome database of the published dataset was used [24]. We checked the expression of 143 VvRCHCs, of them, a total of 136 VvRCHCs expression.To understand the expression of these VvRCHCs under the drought treatment, we used the ratio of WS (Water Stress) to WW (Well-Watered) gene expression of the two genotypes to draw an expression heat map, expression values are reported as log2(WS/WW) fold change (Fig. 5a), the differential multiple matrixes of these VvRCHCs is recorded in Table S4. However, more than 60 % of the VvRCHCs in the two genotypes were highly expressed under the imposed drought. To screen out the key genes, in each time period of the treatment, the gene that conforms to |log2 (WS/WW)| > 1 is considered a differential gene, and the Venn diagram was made using the differentially screened genes of the drought-tolerant genotype M4 at different times (Fig. 5b). By looking at the different genes in each period, there are finally 8 genes that are different in three periods. To robustly verify the gene expression levels, the expression patterns of these 8 genes were verified by qRT-PCR (Fig. 6), whose pattern basically conformed to the trend shown in Fig. 5c. The VyRCHC114 gene was significantly down-regulated at 2 days, with a strong downward trend of the drought treatment. The VyRCHC66, VyRCHC68, VyRCHC69 and VyRCHC95 genes had a similar expression trend, being slightly up-regulated at 2 days of drought, but strongly down-regulated thereafter. These results suggested eight key genes are probably involved in regulating the plant response to drought.

Fig. 5
figure 5

Differential heat map of VvRCHCs expression in plants under drought stress conditions. a Heat maps of two different genotypes, based on their log2(WS/WW) values from the RNA-Seq data set, under drought stress and normal conditions at times T1–T4. b The Venn diagram of DEGs obtained from the analysis of the expression of M4 genotype at different periods. c Selected eight candidate genes in genotype M4 at different stages of WS/WW heat map

Fig. 6
figure 6

Expression of 8 candidate genes were screened in qRT-PCR for plants under drought stress and control conditions. The x-axis represents the different days during the treatment and the y-axis the relative levels of a gene’s expression. Each treatment group had three biological repeats whose averages are plotted with the standard deviation. The asterisks indicate the significant level (* P < 0.05, ** P < 0.01)

Identification of E3 ubiquitin ligase activity of VyRCHC114

To clarify whether VyRCHC114 has E3 ubiquitin ligase activity, we conducted an in vitro ubiquitin activity assay, achieved by using purified MBP-VyRCHC114 fusion protein mixed with ubiquitin, E1, and E2 and by western blotting with the MBP antibody. Ubiquitin molecules were detected on the fusion protein linked by MBP antibody (Fig. 7a). This same method was used to detect ubiquitin antibody tags. The VyRCHC114 protein was detected in the fusion protein linked by the ubiquitin antibody, which indicated it had E3 ligase activity.

Fig. 7
figure 7

E3 ubiquitin ligase activity of VyRCHC114. a Determination of E3 ubiquitin ligase activity of VyRCHC114; an immunoblot analysis was performed with the ubiquitin antibody (right) and MBP antibody (left). b Determination of E3 ubiquitin ligase activity of VyRCHC114 mutants; an immunoblot analysis was performed with ubiquitin antibody (right) and MBP antibody (left). c Schematic diagram of VyRCHC114 C3H2C3 domain and putative mutation sites. C328S and N355A affect a non-conserved site of the VyRCHC114 C3H2C3 domain. Mutations in C320S and H341A affect the ubiquitin activity of VyRCHC114

We know that the RING-C3H2C3 type protein can form a RING structure for ubiquitin regulation, but this process depends on the interaction between the eight conserved metal ligands. To further illustrate whether and how E3 ligase activity of VyRCHC114 depends on these conserved metal ligands, as shown in Fig. 7c, we selected four different amino acid sites for mutation (two key conservative and two non-conservative metal ligand sites). Four corresponding proteins (C320S, C328S, H341A, N355A) were obtained, and their ubiquitin activity in vitro was tested by the same method. After the immuno-blotting analysis of MBP antibody and ubiquitin antibody, evidently the two mutant proteins C320S and H341A lost their E3 ubiquitin ligase activity due to mutations at key sites, but the two mutant proteins C328S and N355A maintained theirs (Fig. 7b). The unprocessed original image is in Fig. S6. These results indicated these conserved metal ligand sites are crucial factors for demonstrating the VyRCHC114 ligase activity.

Overexpression of VyRCHC114 enhances Arabidopsis drought tolerance

To clarify the effects of VyRCHC114’s role in plant responses to drought, we selected transgenic Arabidopsis (OE #2, #5, #13) with high expression levels of the VyRCHC114 gene for subsequent experiments (Fig. 8b). After 15 days of drought imposed upon wild plants and transgenic plants, followed by normal watering for 6 days, phenotype observations revealed that plants overexpressing VyRCHC114 had significantly improved the drought tolerance (Fig. 8a). Further, on average, more than 70 % of the plants overexpressing VyRCHC114 survived the drought stress, which was significantly higher than the 30 % survival rate of the EV-transformed group (Fig. 8c).

Fig. 8
figure 8

VyRCHC114 overexpression (OE) enhances drought resistance in Arabidopsis. a Phenotypes of three transgenic and an EV-transformed Arabidopsis lines after 15 days of drought stress and a 6-day recovery period. b, c Relative expression levels of VyRCHC114 and survival of transgenic and EV-transformed Arabidopsis plants. Data are the mean ± SD (standard deviation). The asterisk, (*) and (**), indicate that OEs and EV-transformed groups were significantly different at P < 0.05 and P < 0.01 (Student’s t-test)

To understand the relationship between plant growth and drought resistance, electrolyte leakage rates (Fig. 9a) and chlorophyll content (Fig. 9b) were both measured. These were similar between VyRCHC114-overexpressed and EV-transformed plants in the non-stress treatment, but after 8 days of drought stress, the electrolyte permeability of the former was significantly lower than the latter, while the chlorophyll content was significantly higher in overexpressing than EV-transformed plants. Additionally, the changes in photosynthesis under drought stress were further analyzed by measuring potential photosynthetic efficiency (Fig. 9c) and capacity storage capacity (Fig. 9d). Each was not significantly different from EV-transformed and VyRCHC114 overexpression plants under non-stress; however, Fv/Fm was significantly higher in the latter than the former at 4 days, and especially at 7 days, of drought stress. At 4 days, energy storage capacity of VyRCHC114-overexpressed plants was not significantly different from that of EV-transformed plants, but at 7 days of drought stress, that of the former exceeded the latter. Hence, these results suggest that VyRCHC114 can enhance the drought resistance of plants by participating in the regulation of photosynthesis.

Fig. 9
figure 9

Physiological indices of the EV-transformed and overexpressing (OE) Arabidopsis plants after drought stress. a Electrolyte leakage, b chlorophyll content, c PSII maximal photochemical efficiency (Fv/Fm), d net photosynthetic rate of leaves (Pn) were evaluated. Data are the mean ± SD (standard deviation). The asterisk, (*) and (**), indicates that OEs and EV-transformed groups were significant different at P < 0.05 and P < 0.01 (Student’s t-test)

Many studies have shown that antioxidant enzymes can influence plants’ drought tolerance. Common antioxidant enzymes are ascorbate peroxidase (APX), superoxide dismutase (SOD), peroxidase (POD), and catalase (CAT), so we examined their activity. As Fig. 10 shows, under non-stress conditions, the activity of these antioxidant enzymes was similar between the plants, whereas when drought stressed for 4 and 7 days, the activities of APX (Fig. 10a), SOD (Fig. 10b), POD (Fig. 10c) and CAT (Fig. 10d) were significantly higher in plants overexpressing VyRCHC114 than those EV-transformed. Taken together, these data indicate VyRCHC114 may also improve drought tolerance by elevating antioxidant enzyme activity.

Fig. 10
figure 10

VyRCHC114-overexpressed and EV-transformed plants’ activity of various antioxidant enzymes in Arabidopsis. a Under drought stress for 0, 4, and 7 days, are ascorbate peroxidase (APX), b superoxide dismutase (SOD), c peroxidase (POD), d catalase (CAT) activities of overexpressed (OE) and EV-transformed plants were determined. Data are mean ± SD (standard deviation). The asterisk, (*) and (**), indicates that OEs and EV-transformed groups were significantly different at P < 0.05 and P < 0.01 (Student’s t-test)

AtCOR15a, AtERD15, AtP5CS1, and AtRD29A are known to be key genes for regulating plant responses to drought stress. So we quantified their expression of imposed drought. As expected, when non-stressed, there was no significant difference between plants overexpressing VyRCHC114 overexpression and those EV-transformed. By contrast, under drought stress, all four genes were significantly higher in VyRCHC114-overexpressed plants than in those EV-transformed (Fig. 11).

Fig. 11
figure 11

Relative expression levels of drought resistance genes in transgenic and EV-transformed Arabidopsis after drought stress. a AtCOR15a, b AtERD15, c AtP5CS1, d AtRD29A. Data are the mean ± SD (standard deviation). The asterisk, (*) and (**), indicates that overexpressed (OEs) and EV-transformed groups of plants were significantly different at P < 0.05 and P < 0.01 (Student’s t-test)

Discussion

The RING C3H2C3 gene family has since been identified with many plant species [25,26,27,28]. Related studies have shown that RING genes are involved in a variety of biological processes, growth and development and hormonal responses, as well as plant responses to abiotic stresses [29]. However, for grapevine, the RING C3H2C3 gene had not yet been identified in its whole genome, with few reports available on its relevance for grapevine growth and developmental regulation or response to abiotic stress. In our study, we analyzed the whole genome of grapevine for the RING C3H2C3 gene family members. Using the criteria of whether the eight conserved metal ligands are present, a total of 143 non-redundant RING C3H2C3 genes were thus identified. Studies have shown that grapevine’s genome size is about 0.5 times that of tomato, containing 0.75 times as many genes as tomato [30, 31]. According to the known RING C3H2C3genes in tomato, the genes account for 0.58 in grapevine, which lies between the multiples of genome length and the number of genes [26].

Many RING C3H2C3 of E3 ubiquitin ligases belong to the ATL gene family [32]. According to Arabidopsis and tomato RING C3H2C3 genes, we divided grapevine’s RING C3H2C3 genes into six categories (I ~ VI) (Fig. 2). Each group has Arabidopsis or tomato in the same branch. This shows that grapevine genes have sequence similarity with Arabidopsis thaliana and tomato genes. Gene replication can arise from fragment replication, tandem replication, transposable events, and even whole genome replication, which not only provide the evolutionary potential for species to produce new functional traits but also are a main driving force for species differentiation [33, 34]. In the identification of gene families from many species, gene replication events have proven instrumental in their expansion [35]. Studies have shown that tandem replication often occurs in widely and fast evolving gene families, a good example being Nucleotide Binding Sites Leucine Rich Repeat (NBS-LRR) resistance families [36]. Segmental replication is more common in slow evolving gene families, like the MYB gene family [36]. There are 7 pairs of genes in VvRCHCs that are duplicated in tandem, and 2 pairs of genes are duplicated in segments. Tandem replication may be the cause of VvRCHCs expansion, like the WRKY family genes in the autopolyploid Saccharum spontaneum [37]. The collinearity analysis of VvRCHCs with Arabidopsis and tomato explains the homology relationship between grape RCHC gene and tomato is closer. In addition, 54 VvRCHCs were found to be homologous to genes in Arabidopsis and tomato, these genes may be preserved by the ancestors of dicots.

Cis-acting elements in gene promoter regions may be critical for gene regulation. Plant hormone regulation, growth and development, and stress were common upstream of different VyRCHCs. This situation is in fact rather common in RING genes of all species [25,26,27,28]. The ubiquitin-proteasome system has been implicated in the control of the ABA response at different points of the ABA pathway [38]. ABRE is ubiquitous in the promoter of VvRCHCs (Fig. S4), Ubiquitin ligase SDIR1 regulates stress-responsive abscisic acid signal by interacting with ABRE abscisic acid response element [39]. There are ABRE elements in the promoters of VvRCHC3,VvRCHC16 and VvRCHC74, which may be induced by the regulation of abscisic acid. ABA, GA, ethylene, trauma, drought, heat stress, and pathogen response elements are present in the promoter region of OsRING genes of rice plants, for which pathogen infection, SA, ABA, JA, and ethephon (ET) treatments could induce target genes expression to different degrees [40]. A similar analysis of RING gene ZmRHCP1 was recently done in maize [41]. Similarly, there are at least seven types of hormone regulatory binding elements in the promoters of VvRCHC39, VvRCHC65, VvRCHC126 and VvRCHC129, they may be responding to a variety of hormone regulation.

According to the analysis of RNG-Seq data set (Fig. 5a), more than 60 % of VyRCHCs were significantly up-regulated or down-regulated under drought stress, indicating those genes may play a key role in how grapevine responds to drought. Studies have revealed the molecular mechanism of many circular genes involved in the drought stress. For example, in Arabidopsis, XERICO, SDIR1, AtAIRP1, AtAIRP2, AtAIRP3 and AtAIRP4 has been found to play a key role in the drought response of plants. In addition, E3 ubiquitin ligase atrzf1 mutation increased the proline content of Arabidopsis and improved drought tolerance [42]. GpDSR7 encodes an E3 ubiquitin ligase, which overexpressed in Arabidopsis increased its tolerance to drought stress [43]. In our study, focused on genes which were significantly up-regulated or down-regulated at four time periods during the drought treatment according to previous screening methods [44], as they more likely to play a key role in grapevine’s drought stress response (Fig. 5b and c). Comparing the RNA-Seq dataset with the RT-qPCR data (Fig. 6), it was found that the VyRCHC114continues to be downregulated. Hence, we postulated the VyRCHC114 may possess E3 ubiquitin ligase activity and as a negative regulatory factor to respond to drought stress. Verifying this, we detected that VyRCHC114 has E3 ubiquitin ligase activity (Fig. 7) and the VyRCHC114 gene is overexpressed in Arabidopsis (Fig. 8). The results show the overexpression of VyRCHC114 gave Arabidopsis drought resistance. This is not in line with expectations. Use the two reference models in Fig. S5 to explain the experimental results. Pattern 1: After VyRCHC114 degrades target protein A, protein B, which is functionally redundant with A, is strongly activated, which gives plants stronger drought resistance. Genes with redundant functions often come from the same gene family, they generally have similar conserved domains and participate in various life activities together [45]. the five MIR172 have redundancy in the regulation of Arabidopsis meristem size, stem elongation and flowering [46]. AtUBP12 and AtUBP13 have functional redundancy in plant immunity, circadian clock, and photoperiod flowering regulation [47, 48]. Pattern 2: After VyRCHC114 degrades target protein A, protein B, which competes with A, is not inhibited, and continues to regulate downstream gene expression. This kind of competitive relationship is often accompanied by a complicated regulatory network. There are studies describing similar that ERF4 and MYB52 regulate downstream gene expression in an opposite manner by antagonizing each other’s DNA-binding ability through a physical interaction [49]. Fig. S5 show interesting and complex networks. The results that are not in line with expectations have prompted more attention to the target protein of VyRCHC114. The preliminary pattern diagram (Fig. S5) gives us confidence. In addition, overexpression of VyRCHC114 caused changes in many related indexes, including antioxidant enzyme activity, photosynthesis rate, active oxygen metabolism and drought resistance gene expression. The results also indirectly prove that the targeted degradation substrate of the gene may be a key regulatory factor in the process of drought stress, and its degradation strongly activates other elements to resist drought stress, thus giving plants stronger drought resistance. It’s going to be an interesting story. In the next work, we will identify the substrate protein of VyRCHC114, further study its regulatory pathway, sequence the transcriptome of Vitis yeshanensis under drought stress, the co-expression networks under drought stress, and compare the differences of drought resistance of other cultivated species.

Drought stress greatly impacts the photosynthesis of plants, by affecting their photosynthetic rates and carbon metabolic pathways [50]. A lowered rate of photosynthesis can lead to excessive accumulation of reactive oxygen species (ROS), leading to cytotoxicity, membrane lipid peroxidation, and even cell death which can be countered by antioxidant enzymes as a form of plant defense [51, 52]. The overexpression of maize E3 ubiquitin ligase gene in transgenic tobacco can reportedly improve the drought resistance of tobacco [53]. Not only that, other abiotic stresses may also be regulated by photosynthesis, thus enabling plants to adapt to stress conditions [54]. According to our results, VyRCHC114 overexpressing plants maintained a strong photosynthetic rate and energy storage capacity while under drought stress. The reason for this may be an increase in their chloroplast content, pointing to VyRCHC114’s possible involvement in the regulation of chlorophyll biosynthesis pathway as an E3 ubiquitin ligase. Nonetheless, we also examined the expression of genes known to play a major role in drought stress responses to plants [55,56,57]. Our results revealed that the expression levels of these genes were significantly higher in VyRCHC114-overexpressed than EV-transformed Arabidopsis plants. Moreover, antioxidant system may be involved in plant abiotic stress tolerance mediated by the E3 ubiquitin ligase. Here, we provide physiological evidence that VyRCHC114 heterologous expression enhances drought resistance by increasing the activity of antioxidant enzymes, which can scavenge for and eliminate ROS to indirectly reduce membrane damage.

Conclusions

VyRCHC may act as an E3 ligase to mediate substrate degradation through the ubiquitin proteasome mechanism. This interaction may cause the target protein to be labeled by ubiquitin signaling, which leads to proteasome degradation. Since VyRCHC114 likely represents a new class of positive /negative regulatory factors of the drought signal pathway, however, positive/negative depends on the regulatory characteristics of the target protein, but we think the degraded protein is a positive regulator of drought signaling, so that more of this substance may activate drought signaling. Therefore, VyRCHC114 may improve the water retention ability and antioxidant defense of plants by regulating their chlorophyll content and antioxidant system, thus participating in drought stress response. So far, however, the target protein of the plant VyRCHC114 gene has not been determined, nor is the mechanism of augmented SOD, POD, APX and CAT activities clearly understood. In future work, we will focus on the identification of VyRCHC114 target proteins under drought stress and activation mechanisms of the antioxidant system in VyRCHC114-transgenic plants.

Methods

Plant materials

The grapevine variety Vitis yeshanensis was sampled from the field, located in the grape germplasm resource garden of Northwest A&F University. Annual plants are selected for treatment, and the treatment method is the same [58]. Root samples were taken at 0d, 2d, 4d, 6d, 8d and 10d respectively for subsequent experiments. Transgenic and wild type (WT) plants of Arabidopsis thaliana ecotype Columbia (Col-0) plants were grown in vermiculite: perlite (1:1, v/v) mix in plastic pots in a growth chamber. Arabidopsis plants were grown in a soil mix of peat moss, perlite, and vermiculite (3:1:1, v/v/v) under a 12-h/12-h day/night cycle at 25℃ with 60 % relative humidity. For the drought stress treatment, plants were transformed with an empty vector (EV) or to overexpress VyRCHC114 (OE#2, OE#5, or OE#13 lines) [59], all of which were grown on individual MS medium plates for 7 days before transplantation into soil, there were 5 strains in the control and three transgenic lines. This experiment mainly followed previous research methods, albeit slightly modified [60]. After 3 weeks, these plants received a 12-day drought stress treatment (no water provided), after which they were re-watered and their survival recorded 6 days later. All experiments were repeated three times.

Identification of RING-type C3H2C3 genes in the grapevine genome

To identify the C3H2H3 type of RING, the most recent grapevine genome file in the Ensembl Plants Database (http://plants.ensembl.org/index.html) was downloaded and used. The grapevine RING C3H2C3 candidates were identified based on the HMM profiles (PF13639 and PF12678) with an e-value cutoff of 0.01. The screened proteins screened were given to Pfam (http://pfam.xfam.org/search/) and SMART ( ttp://smart.embl-heidelberg.de/), e value less than 0.01. According to the results of PFAM, SMART database protein domain identification, extract RCHC conservative domain sequences in the VyRCHC protein sequence, and use CLUSTALX 2.0 to perform multi-sequence alignment, see if conserved Cys-Cys-Cys-His-His-Cys-Cys-Cys. Finally, 143 proteins that ultimately meet the conditions. The physicochemical properties of each RING-type C3H2C3 protein were predicted using the Protparam online tool (https://web.expasy.org/protparam/). The 143 VvRCHCs were named according to their positional information on the chromosomes.

Bioinformatics analysis of VyRCHCs family

CLUSTALX 2.0 software was used to perform a multiple-sequence alignment of the 143 grapevine genes and the 18 tomato and 19 Arabidopsis RCHC protein sequences; it was also used to manually remove any untrusted gaps at both sequence ends. A phylogenetic tree was generated in MEGA 7.0 using the ML (maximum likelihood) method and bootstrapping with n = 1000 replicates, with all other settings set to their default values; the online EVOLVIEW (https://www.evolgenius.info/evolview/#login/) tool carried out the tree’s visualization. The online program Gene Structure Display Server 2.0 (http://gsds.cbi.pku.edu.cn/) was used to identify the genetic structure the VvRCHCs. Using the MEME online program (http://MEME.nbcr.net/meme/introduction.html), the VvRCHC protein sequences could be analyzed under these parameters: an optimal motif width of 6 ~ 35 and a maximum number of motifs of 12. According to the annotated positions in grapevine genome data, the 143 grapevine VvRCHCs were located on 20 chromosomes. By referring to previous studies, BLASTN was used to compare the CDS sequences of VvRCHCs in grapevine and tomato (e-value = 1 × 10− 10, homology > 75 %). The tandem repeat gene pairs and segment repeat gene pairs of VyRCHCs were also identified [61, 62]. Further, the Ka/Ks ratio between repetitive genes pairs can be used to infer the selection pressure in the process of genome evolution. Next, the MCScanX program (e-value: 1 × 10− 10, num alignments: 5) was used to detect the collinear region between VvRCHCs in grapevine and tomato/Arabidopsis; any collinear gene pair of VvRCHCs was marked with red and green lines. The cis-elements were identified from the upstream 2-kb promoter sequences of the VvRCHCs after submitting them to PlantCARE (http://bioinformatics.psb.ugent.be/webtools/plantcare/html) [63], to obtain their image display, the resulting XML file was uploaded to TBtools [64].

Expression analysis of VyRCHCs in grapevine under drought stres

To analyze the grapevine RCHC genes’ expression levels under drought stress, we obtained from the NCBI database (registration number: SRA110531) two different drought resistance genes (101.14 and M4) which were compared under two different treatments WS(Water Stress) and WW (Well-Watered) in roots and in different periods (T1–T4: 2d、4d、7d、10d) RNA-Seq data set [24]. Based on the expression values of RING C3H2C3 in the roots of the two genotypes, we calculated the log2(WS/WW) values (fold-change) in each time period (Table S4). The R package ‘pheatmap’ was used to produce a heatmap for this data.

RNA extraction and quantitative real-time PCR (qRT-PCR)

The qRT-PCR primers were designed using Primer Premier software (version 5.0). The RNA from Arabidopsis and grapevine (Vitis yeshanensis) leaves was extracted using the Spectrum Plant Total RNA Kit (Sigma-Aldrich, Beijing, China), after which reverse transcription of RNA into cDNA was done using the Prime Script RT Reagent Kit (Takara, Dalian, China). The qRT-PCR was performed in an IQ5 real-time PCR detection system (Bio-Rad Laboratories, Hercules, CA, USA) with SYBR Premium EX Taq II (Takara, Dalian, China). The reaction volume was 25 µl. The relative expression level corresponding to β-TUB4 and ubiquitin1 was calculated by using the 2−ΔΔCt method [65]; each reaction was prepared in triplicate and repeated three times. Primer sequence information in Table S1.

E3 ubiquitin ligase activity assay

The open reading frame (ORF) of VyRCHC114 and the different site mutants C320S, H341A, C328S, and N355A were separately cloned into the SalI/KpnI site of the pMAL-c5X vector (New England Biolabs UK Ltd, Hitchin, UK). According to the manufacturer’s instructions, the pMAL protein fusion and purification system (New England Biolabs) was used to purify the fusion protein. Ubiquitination activity was then measured that according to the method described above [66], albeit with the following modifications made: 250 ng of purified E3 (MBP-VyRCHC114, C320S, H341A, C328S, and N355A) in the ubiquitination buffer (50 mM Tris-HCl (pH 7.5), while the other reagents and steps used were the same. Primer sequence information in Table S1.

Physiological analysis of drought stress response of transgenic Arabidopsis

To determine the water loss rate, 10 leaves were detached from 3-week-old transgenic and WT plants and immediately weighed. The samples were then placed on dry filter paper at a relative humidity of 40–45 % at room temperature and weighed over a time course. Leaves were sampled after dehydration to detect cell death, electrolyte leakage, malondialdehyde, antioxidant enzyme activity. The leaves collected before dehydration were used as a negative control.

For chlorophyll content measurements, approximately 0.05 g of fresh leaf material was placed in 5 ml of 96 % ethanol and incubated at 4◦C in the dark overnight. The absorbance of the extracted pigments was measured at 665 and 649 nm using a spectrophotometer (Hitachi Limited, Tokyo, Japan) and the chlorophyll content was calculated as previously described [58].

Relative electrolyte leakage was measured as previously described [67], as was MDA content [66]. In addition, superoxide dismutase (SOD), peroxidase (POD), catalase (CAT), ascorbate peroxidase (APX) enzyme activities were extracted from 0.5 g leaves from abiotic stress treated plants as well as control plants, and measured as described by [68].

Statistical analysis

All the above experiments by SPSS software (version 21.0) were employed to analyze the statistically significant differences of the gene expression levels by ANOVA with Duncan’s multiple range test. All experiments were repeated three times as independent analyses.

Availability of data and materials

All data generated and analyzed during this study are included in this published article. To identify the C3H2H3 type of RING, the most recent grapevine genome file in the Ensembl Plants Database ( http://plants.ensembl.org/index.html ) was downloaded and used. Expression data of C3H2H3 type of RING genes in grapevine used in this study can be accessed via the NCBI SRA database with accession numbers of SRA110531. Figure S1 Schematic diagram of C3H2C3 conserved protein sequence alignment of VvRCHCs. Figure S2 The original tree of Fig. 2. Figure S3 Number of introns in VvRCHCs. Figure S4 Analysis of cis-acting elements in the VvRCHCs. Figure S5. Two model diagrams (VyRCHC114 is involved in drought resistance). Figure S6 VyRCHC114 in vitro ubiquitin gel imprinting (uncut). Table S1 The sequences of the primers used in these experiments. Table S2 The distance between conserved metal ligands in the C3H2H3 domain of 143 VvRCHCs. Table S3 Functions of the cis-acting elements that found in the promoter region of each of VvRCHCs. Table S4 Differential multiple expression matrix of VvRCHCs.

Abbreviations

RCHCs:

RING C3H2C3-type genes

DEGs:

Differentially expressed genes

qRT-PCR:

Quantitative real-time reverse transcription PCR

RNA-seq:

RNA-sequencing

WS:

Water stress

WW:

Well-watered

ML:

Maximum likelihood

WT:

Wild type

HMM:

Hidden Markov Model

References

  1. Sadanandom A, Bailey M, Ewan R, Lee J, Nelis S. The ubiquitin–proteasome system: central modifier of plant signalling. New Phytol. 2012;196(1):13–28.

    Article  CAS  PubMed  Google Scholar 

  2. Sharma B, Joshi D, Yadav PK, Gupta AK, Bhatt TK. Role of ubiquitin-mediated degradation system in plant biology. Front Plant Sci. 2016;7:806.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Kelley DR. E3 ubiquitin ligases: key regulators of hormone signaling in plants. Mol Cell Proteomics. 2018;17(6):1047–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Stone SL. The role of ubiquitin and the 26S proteasome in plant abiotic stress signaling. Front Plant Sci. 2014;5:135.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Shu K, Yang W. E3 ubiquitin ligases: ubiquitous actors in plant development and abiotic stress responses. Plant Cell Physiol. 2017;58(9):1461–76.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Qi S, Lin Q, Zhu H, Gao F, Zhang W, Hua X. The RING finger E3 ligase SpRing is a positive regulator of salt stress signaling in salt-tolerant wild tomato species. Plant Cell Physiol. 2016;57(3):528–39.

    Article  CAS  PubMed  Google Scholar 

  7. Zhang H, Cui F, Wu Y, Lou L, Liu L, Tian M, Ning Y, Shu K, Tang S, Xie Q. The RING finger ubiquitin E3 ligase SDIR1 targets SDIR1-INTERACTING PROTEIN1 for degradation to modulate the salt stress response and ABA signaling in Arabidopsis. Plant Cell. 2015;27(1):214–27.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Liu J, Zhang C, Wei C, Liu X, Wang M, Yu F, Xie Q, Tu J. The RING finger ubiquitin E3 ligase OsHTAS enhances heat tolerance by promoting H2O2-induced stomatal closure in rice. Plant Physiol. 2016;170(1):429–43.

    Article  CAS  PubMed  Google Scholar 

  9. Ding S, Zhang B, Qin F. Arabidopsis RZFP34/CHYR1, a ubiquitin E3 ligase, regulates stomatal movement and drought tolerance via SnRK2. 6-mediated phosphorylation. Plant Cell. 2015;27(11):3228–44.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Lim CW, Baek W, Lee SC. The pepper RING-type E3 ligase CaAIRF1 regulates ABA and drought signaling via CaADIP1 protein phosphatase degradation. Plant Physiol. 2017;173(4):2323–39.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Brugière N, Zhang W, Xu Q, Scolaro EJ, Lu C, Kahsay RY, Kise R, Trecker L, Williams RW, Hakimi S. Overexpression of RING domain E3 ligase ZmXerico1 confers drought tolerance through regulation of ABA homeostasis. Plant Physiol. 2017;175(3):1350–69.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  12. Ryu MY, Cho SK, Kim WT. The Arabidopsis C3H2C3-type RING E3 ubiquitin ligase AtAIRP1 is a positive regulator of an abscisic acid-dependent response to drought stress. Plant Physiol. 2010;154(4):1983–97.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Cho SK, Ryu MY, Seo DH, Kang BG, Kim WT. The Arabidopsis RING E3 ubiquitin ligase AtAIRP2 plays combinatory roles with AtAIRP1 in abscisic acid-mediated drought stress responses. Plant Physiol. 2011;157(4):2240–57.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Kim JH, Kim WT. The Arabidopsis RING E3 ubiquitin ligase AtAIRP3/LOG2 participates in positive regulation of high-salt and drought stress responses. Plant Physiol. 2013;162(3):1733–49.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Park C, Lim CW, Baek W, Lee SC. RING type E3 ligase CaAIR1 in pepper acts in the regulation of ABA signaling and drought stress response. Plant Cell Physiol. 2015;56(9):1808–19.

    Article  CAS  PubMed  Google Scholar 

  16. Yu Y, Bian L, Yu K, Yang S, Zhang G, Guo D. Grape (Vitis davidii) VdGATA2 functions as a transcription activator and enhances powdery mildew resistance via the active oxygen species pathway. Sci Hortic. 2020;267:109327.

    CAS  Google Scholar 

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

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  18. Hardie WJ, Considine JA. Response of grapes to water-deficit stress in particular stages of development. Am J Enol Viticult. 1976;27(2):55–61.

    Article  Google Scholar 

  19. dos Santos TP, Lopes CM, Rodrigues ML, de Souza CR, Maroco JP, Pereira JS, Silva JR, Chaves MM. Partial rootzone drying: effects on growth and fruit quality of field-grown grapevines (Vitis vinifera). Funct Plant Biol. 2003;30(6):663–71.

    Article  PubMed  Google Scholar 

  20. Chaves MM, Flexas J, Pinheiro C. Photosynthesis under drought and salt stress: regulation mechanisms from whole plant to cell. Ann bot. 2009;103(4):551–60.

    Article  CAS  PubMed  Google Scholar 

  21. Liu J, Zhao F, Guo Y, Fan X, Wang Y, Wen Y. The ABA receptor-like gene VyPYL9 from drought-resistance wild grapevine confers drought tolerance and ABA hypersensitivity in Arabidopsis. Plant Cell Tiss Org. 2019;138(3):543–58.

    Article  CAS  Google Scholar 

  22. Yang Y, He M, Zhu Z, Li S, Xu Y, Zhang C, Singer SD, Wang Y. Identification of the dehydrin gene family from grapevine species and analysis of their responsiveness to various forms of abiotic and biotic stress. BMC Plant Biol. 2012;12(1):140.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Jing Z, Wang X. Genetic relationship between Chinese wild Vitis species and American and European cultivars based on ISSR markers. Biochem Syst Ecol. 2013;46:120–6.

    Article  CAS  Google Scholar 

  24. Corso M, Vannozzi A, Maza E, Vitulo N, Meggio F, Pitacco A, Telatin A, D’Angelo M, Feltrin E, Negri AS, et al. Comprehensive transcript profiling of two grapevine rootstock genotypes contrasting in drought susceptibility links the phenylpropanoid pathway to enhanced tolerance. J Exp Bot. 2015;66(19):5739–52.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Li Y, Wu B, Yu Y, Yang G, Wu C, Zheng C. Genome-wide analysis of the RING finger gene family in apple. Mol Genet Genomics. 2011;286(1):81.

    Article  CAS  PubMed  Google Scholar 

  26. Yang L, Miao M, Lyu H, Cao X, Li J, Li Y, Li Z, Chang W. Genome-wide identification, evolution, and expression analysis of RING finger gene family in solanum lycopersicum. Int J Mol Sci. 2019;20(19):4864.

  27. Qanmber G, Yu D, Li J, Wang L, Ma S, Lu L, Yang Z, Li F. Genome-wide identification and expression analysis of Gossypium RING-H2 finger E3 ligase genes revealed their roles in fiber development, and phytohormone and abiotic stress responses. J Cot Res. 2018;1(1):1–17.

    Article  Google Scholar 

  28. Alam I, Yang YQ, Wang Y, Zhu ML, Wang HB, Chalhoub B, Lu YH. Genome-wide identification, evolution and expression analysis of RING finger protein genes in Brassica rapa. Sci Rep. 2017;7:40690.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Hua Z, Vierstra RD. The cullin-RING ubiquitin-protein ligases. Annu Rev Plant Biol. 2011;62:299–334.

    Article  CAS  PubMed  Google Scholar 

  30. Jaillon O, Aury J, Noel B, Policriti A, Clepet C, Casagrande A, Choisne N, Aubourg S, Vitulo N, Jubin C. The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla. Nature. 2007;449(7161):463.

    Article  CAS  PubMed  Google Scholar 

  31. Consortium TG. The tomato genome sequence provides insights into fleshy fruit evolution. Nature. 2012;485(7400):635.

    Article  CAS  Google Scholar 

  32. Serrano M, Parra S, Alcaraz LD, Guzmán P. The ATL gene family from Arabidopsis thaliana and Oryza sativa comprises a large number of putative ubiquitin ligases of the RING-H2 type. J Mol Evol. 2006;62(4):434–45.

    Article  CAS  PubMed  Google Scholar 

  33. Ren R, Wang H, Guo C, Zhang N, Zeng L, Chen Y, Ma H, Qi J. Widespread whole genome duplications contribute to genome complexity and species diversity in angiosperms. Mol Plant. 2018;11(3):414–28.

    Article  CAS  PubMed  Google Scholar 

  34. Moore RC, Purugganan MD. The early stages of duplicate gene evolution. Proc Nat Acad Sci USA. 2003;100(26):15682–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Wang N, Liu Y, Cong Y, Wang T, Zhong X, Yang S, Li Y, Gai J. Genome-wide identification of soybean U-Box E3 ubiquitin ligases and roles of GmPUB8 in negative regulation of drought stress response in Arabidopsis. Plant Cell Physiol. 2016;57(6):1189–209.

    Article  CAS  PubMed  Google Scholar 

  36. Cannon SB, Mitra A, Baumgarten A, Young ND, May G. The roles of segmental and tandem gene duplication in the evolution of large gene families in Arabidopsis thaliana. BMC Plant Biol. 2004;4(1):10.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Li Z, Hua X, Zhong W, Yuan Y, Wang Y, Wang Z, Ming R, Zhang J. Genome-wide identification and expression profile analysis of WRKY family genes in the autopolyploid saccharum spontaneum. Plant Cell Physiol. 2019;61(3):616–30.

    Article  CAS  Google Scholar 

  38. Park J, Nguyen KT, Park E, Jeon J-S, Choi G. DELLA proteins and their interacting RING finger proteins repress gibberellin responses by binding to the promoters of a subset of gibberellin-responsive genes in Arabidopsis. Plant Cell. 2013;25(3):927–43.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Zhang Y, Yang C, Li Y, Zheng N, Chen H, Zhao Q, Gao T, Guo H, Xie Q. SDIR1 is a RING finger E3 ligase that positively regulates stress-responsive abscisic acid signaling in Arabidopsis. Plant Cell. 2007;19(6):1912–29.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Meng XB, Zhao WS, Lin RM, Wang M, Peng YL. Molecular cloning and characterization of a rice blast-inducible RING-H2 type zinc finger gene: full length research paper. DNA Seq. 2006;17(1):41–8.

    Article  CAS  PubMed  Google Scholar 

  41. Li WL, Sun Q, Li WC, Yu YL, Zhao M. Meng ZD: Characterization and expression analysis of a novel RING-HC gene, ZmRHCP1, involved in brace root development and abiotic stress responses in maize. J Integr Agr. 2017;16(9):1892–9.

  42. Ju HW, Min JH, Chung MS, Kim CS. The atrzf1 mutation of the novel RING-type E3 ubiquitin ligase increases proline contents and enhances drought tolerance in Arabidopsis. Plant Sci. 2013;203:1–7.

    Article  PubMed  CAS  Google Scholar 

  43. Li M, Li Y, Zhao J, Liu H, Jia S, Li J, Zhao H, Han S, Wang Y. GpDSR7, a novel E3 Ubiquitin ligase gene in Grimmia pilifera is involved in tolerance to drought stress in Arabidopsis. PloS One. 2016;11(5):e0155455.

  44. Ji XR, Yu YH, Ni PY, Zhang GH, Guo DL. Genome-wide identification of small heat-shock protein (HSP20) gene family in grape and expression profile during berry development. BMC Plant Biol. 2019;19(1):433.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  45. Shanks CM, Hecker A, Cheng CY, Brand L, Collani S, Schmid M, Schaller GE, Wanke D, Harter K, Kieber JJ. Role of BASIC PENTACYSTEINE transcription factors in a subset of cytokinin signaling responses. Plant J. 2018;95(3):458–73.

    Article  CAS  PubMed  Google Scholar 

  46. Lian H, Wang L, Ma N, Zhou CM, Han L, Zhang TQ, Wang JW. Redundant and specific roles of individual MIR172 genes in plant development. PLoS Biol. 2021;19(2):e3001044.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Ewan R, Pangestuti R, Thornber S, Craig A, Carr C, O’Donnell L, Zhang C, Sadanandom A. Deubiquitinating enzymes AtUBP12 and AtUBP13 and their tobacco homologue NtUBP12 are negative regulators of plant immunity. New Phytol. 2011;191(1):92–106.

    Article  CAS  PubMed  Google Scholar 

  48. Cui X, Lu F, Li Y, Xue Y, Kang Y, Zhang S, Qiu Q, Cui X, Zheng S, Liu B, Xu X, Cao X. Ubiquitin-specific proteases UBP12 and UBP13 act in circadian clock and photoperiodic flowering regulation in Arabidopsis. Plant Physiol. 2013;162(2):897–906.

  49. Ding A, Tang X, Yang D, Wang M, Ren A, Xu Z, Hu R, Zhou G, O’Neill M, Kong Y. ERF4 and MYB52 transcription factors play antagonistic roles in regulating homogalacturonan de-methylesterification in Arabidopsis seed coat mucilage. Plant Cell. 2021;33(2):381–403.

    Article  PubMed  Google Scholar 

  50. Teng K, Li J, Liu L, Han Y, Du Y, Zhang J, Sun H, Zhao Q. Exogenous ABA induces drought tolerance in upland rice: the role of chloroplast and ABA biosynthesis-related gene expression on photosystem II during PEG stress. Acta physiol plant. 2014;36(8):2219–27.

    Article  CAS  Google Scholar 

  51. Apel K, Hirt H. Reactive oxygen species: metabolism, oxidative stress, and signal transduction. Annu Rev Plant Biol. 2004;55:373–99.

    Article  CAS  PubMed  Google Scholar 

  52. Park JJ, Yi J, Yoon J, Cho LH, Ping J, Jeong HJ, Cho SK, Kim WT, An G. OsPUB15, an E3 ubiquitin ligase, functions to reduce cellular oxidative stress during seedling establishment. Plant J. 2011;65(2):194–205.

    Article  CAS  PubMed  Google Scholar 

  53. Liu J, Xia Z, Wang M, Zhang X, Yang T, Wu J. Overexpression of a maize E3 ubiquitin ligase gene enhances drought tolerance through regulating stomatal aperture and antioxidant system in transgenic tobacco. Plant Physiol Biochem. 2013;73:114–20.

    Article  CAS  PubMed  Google Scholar 

  54. Le Hir R, Castelain M, Chakraborti D, Moritz T, Dinant S, Bellini C. At bHLH68 transcription factor contributes to the regulation of ABA homeostasis and drought stress tolerance in Arabidopsis thaliana. Physiol Plant. 2017;160(3):312–27.

    Article  PubMed  CAS  Google Scholar 

  55. Hsieh T, Li C, Su R, Cheng C, Tsai Y, Chan M. A tomato bZIP transcription factor, SlAREB, is involved in water deficit and salt stress response. Planta. 2010;231(6):1459–73.

    Article  CAS  PubMed  Google Scholar 

  56. Yu D, Zhang L, Zhao K, Niu R, Zhai H, Zhang J. VaERD15, a transcription factor gene associated with cold-tolerance in Chinese Wild Vitis amurensis. Front Plant Sci. 2017;8:297.

    Article  PubMed  PubMed Central  Google Scholar 

  57. Ma L, Li Y, Chen Y, Li X. Improved drought and salt tolerance of Arabidopsis thaliana by ectopic expression of a cotton (Gossypium hirsutum) CBF gene. Plant Cell Tiss Org. 2015;124(3):583–98.

    Article  CAS  Google Scholar 

  58. Hou H, Jia H, Yan Q, Wang X. Overexpression of a SBP-Box Gene (VpSBP16) from Chinese wild vitis species in Arabidopsis improves salinity and drought stress tolerance. Int J Mol Sci. 2018;19(4):940.

    Article  PubMed Central  CAS  Google Scholar 

  59. Yu Y, Bian L, Yu K, Yang S, Zhang G, Guo D. Grape (Vitis davidii) VdGATA2 functions as a transcription activator and enhances powdery mildew resistance via the active oxygen species pathway. Sci Hortic. 2020;267:109327.

    Article  CAS  Google Scholar 

  60. Guo R, Qiao H, Zhao J, Wang X, Tu M, Guo C, Wan R, Li Z, Wang X. The grape VlWRKY3 gene promotes abiotic and biotic stress tolerance in transgenic Arabidopsis thaliana. Front Plant Sci. 2018;9:545.

    Article  PubMed  PubMed Central  Google Scholar 

  61. Moore RC, Purugganan MD. The early stages of duplicate gene evolution. Proc Nat Acad Sci U S A. 2003;100(26):15682–7.

  62. Ni P, Ji X, Guo D. Genome–wide identification, characterization, and expression analysis of GDSL-type esterases/lipases gene family in relation to grape berry ripening. Sci Hortic. 2020;264:109162.

    Article  CAS  Google Scholar 

  63. Lescot M, Déhais P, Thijs G, Marchal K, Moreau Y, Van de Peer Y, Rouzé P, Rombauts S. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Chen C, Chen H, Zhang Y, Thomas H, Frank M, He Y, Xia R. TBtools - an integrative toolkit developed for interactive analyses of big biological data. Mol Plant. 2020;13(8):1194–202.

    Article  CAS  PubMed  Google Scholar 

  65. Livak KJ, Schmittgen T. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001;25(4):402–8.

    Article  CAS  PubMed  Google Scholar 

  66. Yu Y, Xu W, Wang J, Wang L, Yao W, Yang Y, Xu Y, Ma F, Du Y, Wang Y. The Chinese wild grapevine (Vitis pseudoreticulata) E3 ubiquitin ligase Erysiphe necator-induced RING finger protein 1 (EIRP1) activates plant defense responses by inducing proteolysis of the VpWRKY11 transcription factor. New Phytol. 2013;200(3):834–46.

    Article  CAS  PubMed  Google Scholar 

  67. Zhang L, Zhao G, Xia C, Jia J, Liu X, Kong X. A wheat R2R3-MYB gene, TaMYB30-B, improves drought stress tolerance in transgenic Arabidopsis. J Exp Bot. 2012;63(16):5873–85.

    Article  CAS  PubMed  Google Scholar 

  68. Chen R, Wu P, Cao D, Tian H, Chen C, Zhu B. Edible coatings inhibit the postharvest berry abscission of table grapes caused by sulfur dioxide during storage. Postharvest Biol Tec. 2019;152:1–8.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

We would like to thank all the colleagues in our laboratory for providing useful discussions and technical assistance.

Funding

This work was financially supported by grants from the Natural Science Foundation of China (no. NSFC:31701893 and U1504321), Zhongyuan Science and Technology Innovation Leaders (no. 194200510007), and Innovative Research Team in Henan University of Science & Technology, China (no. 2015TTD003), and Program for Innovative Research Team (in Science and Technology) in University of Henan Province (21IRTSTHN021). The Funding bodies were not involved in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Author information

Authors and Affiliations

Authors

Contributions

Y.H.Y. and S.D.Y. conceived the original screening and research plans. S.D.Y. performed most of the experiments and finished the writing of the paper. D.L.G. participated in the modification of the paper L.B., K.K.Y., X.X.M. and G.H.Z. assisted in sorting out the data. W.R.X. and W.K.Y. provided guidance for experiments. Y.H.Y. and S.D.Y. wrote the article with contributions of all the authors. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Dalong Guo.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors have no conflict of interest to report.

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: Figure S1.

Schematic diagram of C3H2C3 conserved sequence alignment of VvRCHCs. Cys(C) and His(H) amino acids were added on a blue and pink background. The C3H2C3 conserved amino acid sequence length of these genes is shown later in the sequence.

Additional file 2: Figure S2.

The original tree of Fig. 2.

Additional file 3: Figure S3.

Number of introns in VvRCHCs.

Additional file 4: Figure S4.

Analysis of Cis-acting elements in the VvRCHCs. a A list of VvRCHCs to facilitate correspondence. b Color-coded numbers of cis-acting elements of the three major types of VvRCHCspromoters. c Different types of cis-actingelements are represented by different colored squares and their position on each VvRCHC gene promoter.

Additional file 5: Figure S5.

Two model diagrams (VyRCHC114 is involved in drought resistance).

Additional file 6: Figure S6.

VyRCHC114 in vitro ubiquitin gel imprinting (uncut). Corresponding to the original image of gel imprinting in Fig.7A and B.

Additional file 7: Table S1.

The sequences of the primers used in these experiments.

Additional file 8: Table S2.

The distance between conserved metal ligands in the C3H2H3 domain of 143 VvRCHCs.

Additional file 9: Table S3.

Functions of the cis-acting elements that found in the promoter region of each of VyRCHCs.

Additional file 10: Table S4.

The number of VvRCHCs expressed in the roots of two different genotypes was expressed as log2(WS/WW), and T1-T4 represented different drought stress time.

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

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Yu, Y., Yang, S., Bian, L. et al. Identification of C3H2C3-type RING E3 ubiquitin ligase in grapevine and characterization of drought resistance function of VyRCHC114. BMC Plant Biol 21, 422 (2021). https://0-doi-org.brum.beds.ac.uk/10.1186/s12870-021-03162-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12870-021-03162-8

Keywords