- Research article
- Open Access
Genome-wide identification and characterization of the 14–3-3 family in Vitis vinifera L. during berry development and cold- and heat-stress response
BMC Genomics volume 19, Article number: 579 (2018)
The 14–3-3 family of ubiquitous proteins in eukaryotes plays important roles in the regulation of various plant biological processes. However, less information is known about this family in grape fruit.
To investigate the characteristics and functions of 14–3-3 in grape, a total of 11 14–3-3 proteins were identified. Phylogenetic analysis of 14–3-3 proteins in grape (VviGRFs) with homologous proteins in Arabidopsis showed that these proteins were classified into two groups, namely, epsilon and non-epsilon groups. Epsilon group members commonly contained more introns and motifs than non-epsilon group, and some intron positions were found to be conserved between Vitis and Arabidopsis 14–3-3 genes. RNA-seq and qRT-PCR results indicated that VviGRF genes may be involved in the regulation of grape development and berry ripening. Moreover, six VviGRFs exhibited significantly up- or down-regulated expression in response to cold and heat stresses, thereby revealing their potential roles in the regulation of abiotic stress responses.
This work provides fundamental knowledge for further studies about the biological roles of VviGRFs in grape development and abiotic stress response. The present result will also be beneficial for understanding their molecular mechanisms and improving grape agricultural traits in the future.
The 14–3-3 proteins were first found in bovine brain and named according to their elution and migration pattern on 2D DEAE-cellulose chromatography and starch gel electrophoresis . This proteins, as highly conserved regulatory proteins, are ubiquitous in eukaryotes, including plants, insects, mammals, amphibians and yeasts [2,3,4,5]. Generally, these proteins exist as homo- or hetero-isoform dimers, which form dimer groove structures . The structures provide binding sites for 14–3-3 proteins to interact with their targets, such as enzymes involved in primary biosynthetic and energy metabolism (e.g. nitrate reductase, brassinosteroid receptor kinase and H+-ATPase) and signal proteins (e.g. lipoxygenase and protein kinases) [7,8,9]. With the interaction, 14–3-3 proteins function as scaffolding proteins that are involved in the regulation of diverse biological processes, including carbon and nitrogen metabolisms, tricarboxylic acid cycle and the shikimate pathway [8, 10].
To date, 8, 12, 25, 25, 13, 18 and 10 14–3-3 proteins were identified in rice, tomato, cotton, banana, Arabidopsis, soybean and Hevea brasiliensis, respectively [11,12,13,14,15,16,17]. Previous studies revealed that 14–3-3 proteins play important roles in plant growth and development. For example, most banana MaGRFs show an accumulated transcription during fruit development and post-harvest ripening . Various GhGRF interaction partners identified in cotton fibre are involved in plant development, metabolism, signalling transduction and other cellular processes. The overexpression of GhGRFl promotes cotton fibre elongation and maturation . In Arabidopsis, antisense technology revealed that AtGRF10 and AtGRF9 regulate starch accumulation . The 14–3-3 proteins are also related to various abiotic stress responses. In rice, most OsGRF expression changes in response to heat, low temperature and salt stresses . Overexpressed AtGRF6 transgenic cotton shows a stay-green phenotype and improved tolerance to drought stress .
Grape (Vitis vinifera L.) is one of the most important fruit crops worldwide; this crop has been widely cultivated because of its nutritional value and wide application . The grape production is often threatened with abiotic stresses, such as heat or low temperature, drought and salt stresses. Currently, the14–3-3 family in grape has not been studied systematically, except for two 14–3-3 genes (Vv1CS and Vv2CS) that are cloned from V. vinifera cv. ‘Cabernet Sauvignon’; nevertheless, the functions of these genes are still unclear . Several unknown grape 14–3-3-like proteins are also involved in biotic and abiotic stress responses . Thus, 14–3-3 gene family should be systematically analysed to characterise the expression during development and response to abiotic stress in grape.
In the present study, the identification, phylogeny, structure and gene location and duplication were performed on the 14–3-3 family in grape. The expression levels of 14–3-3 s in various tissues and organs of grape during different developmental stages and responses to heat and cold stresses were systematically studied in detail. The present results will provide important information for further study of the regulation mechanism of 14–3-3 s in grape berry development and under abiotic stress.
Identification of 14–3-3 proteins in grape
The 14–3-3 proteins in grape were identified according to 12X V1 version (Grape Genome Database, http://genomes.cribi.unipd.it/DATA/, Table 1). The 14–3-3 protein sequences of Oryza. sativa and A. thaliana used in this study were obtained from NCBI (http://www.ncbi.nlm.nih.gov/).
Hidden Markov Model (HMM), BLASTP program and NCBI-Conserved Domain Data (CDD) search were used to identify the 14–3-3 members in grape. HMM was constructed to scan 14–3-3 proteins in grape by using the obtained OsGRF and AtGRF protein sequences. OsGRFs, AtGRFs, together with 14–3-3 proteins published in NCBI were used as query sequences to blast against the grape protein database obtained from the 12X V1 version. Finally, 11 proteins containing 14–3-3 domain were identified as members of 14–3-3 in grape using NCBI-CDD (http://0-www.ncbi.nlm.nih.gov.brum.beds.ac.uk/Structure/cdd/wrpsb.cgi) .
Multiple sequence alignment and phylogenetic analysis
Identified 14–3-3 proteins in grape were aligned with AtGRFs using ClustalX2 program . The secondary structures of VviGRF proteins were analysed by searching in Self-Optimized Prediction Methods with Alignment (SOPMA, https://npsa-prabi.ibcp.fr/cgi-bin/npsa_automat.pl?page=npsa_sopma.html) . With the combined aligned results, typical α helices in 14–3-3 proteins were marked on their corresponding sequences (Fig. 1). Phylogenetic analysis was conducted by MEGA7 based on the aligned results . Neighbour-joining method was used with bootstrap replications of 1000 (Figs. 2 and 3a). 11 VviGRF proteins were named based on their domains and their phylogenetic relationships with AtGRFs . When a one-to-one orthology was present in the Arabidopsis 14–3-3 family, the grape 14–3-3 members were given the corresponding Arabidopsis-like name (e.g., AtGRF11 and VviGRF11). Two grape genes had the same phylogenetic distance from a single homologue in Arabidopsis were differentiated by a number (e.g., AtGRF9, VviGRF9a and VviGRF9b). As for remaining 14–3-3 members, there are no Arabidopsis genes homologous to them, thus we named them on the basis of domain analysis: both VIT_00s0199g00140 and VIT_02s0033g00780 have incomplete 14–3-3 domains, thus these two protein were named as VviGRF-like1 and –like2. Since there are only 13 AtGRFs in Arabidopsis 14–3-3 family, other five members in several adjacent subclasses share the same 14–3-3 domain were named from VviGRF14 to 15, 16, 17 and 18. The isoelectric point (PI) and molecular weight (Mol wt) of VviGRF proteins were examined on Swiss-Prot website (http://web.expasy.org/compute_pi/, Table 1).
Gene structure and motif analysis
The gene structure data used in this study were obtained from the 12X V1 genome annotation file (http://genomes.cribi.unipd.it/DATA/) and The Arabidopsis Information Resource website (TAIR, http://www.arabidopsis.org/browse/genefamily/14-3-3.jsp). Gene structure analysis was conducted using Gene Structure Display Server version 2.0 (http://gsds.cbi.pku.edu.cn/, Fig. 3b) . According to the amino acid sequences, the motifs of VviGRF proteins were analysed by using MEME program (http://meme-suite.org/tools/meme) . The motif distribution type was zero or one occurrence per sequence, and only motifs with E-value > 0.05 were present (Fig. 3b and Additional file 1).
Gene location and duplication
The location data of VviGRFs were also obtained from the genome annotation files (http://genomes.cribi.utnipd.it/DATA/). The gene location map was constructed using MapChart (Fig. 4) . Moreover, the gene duplication landscape was obtained using the MCScanX . Each duplicate segment with VviGRF genes was selected, and the syntonic map was generated using CIRCOS . The putative duplicated genes were linked by the connection lines (Fig. 5).
In this work, published data GSE36128, GSE62744, SRP018199 and SRP091989 were used for the analysis of expression profiles of VviGRF genes in 54 tissues and organs, during berry ripening, and in cold and heat stress responses, respectively [34,35,36,37]. To survey the temporal and tissue-specific expression of VviGRFs, 54 tissues and organs in V. vinifera cv Corvina were collected in GSE36128 (Additional file 2, Fig. 6). To search the expression pattern of VviGRFs in different varieties during berry development, berries at four growth periods (Pea, Touch, Soft, Harvest) in five grape varieties, namely, Sangiovese, Barbera, Negro amaro, Refosco and Primitivo were used in GSE62744 (Fig. 7a). In SRP018199, the V. vinifera cv. Muscat Hamburg seedlings were used under the gradual cooling treatment (decreased at 5 °C per hour from 24 °C to 4 °C and then 4 °C for additional 4 h) to analyse the expression profile of VviGRFs in response to cold treatment. Besides, to survey their response to heat treatment, detached leaves in V. vinifera cv ‘Jingxiangyu’ plantlets were treated under 25 °C, 35 °C, 40 °C and 45 °C.
To verify exon—intron structure, specific cDNA and DNA sequences were obtained from three biological replicates of green stem in V. vinifera cv. ‘Cabernet Sauvignon’. To analyse the 14–3-3 expression pattern in grape during berry ripening, three biological replicates represented by 15—20 grape berries from V. vinifera cv. ‘Cabernet Sauvignon’ were sampled at three developmental stages (EL33, EL35 and EL37 ). The grapevines were planted in south-to-north oriented rows by similar management conditions at the Germplasm Repository for Grapevines, Institute of Botany, the Chinese Academy of Sciences, Beijing, China (39° 54’ N, 116° 23′ E).
To analyse the 14–3-3 expression in grape under heat stress, 30-day-old leaves taken from V. davidii were heated in a water bath (25 °C for 2 h and 47 °C for 40 min) with three biological replicates. V. davidii were planted in the same condition as V. vinifera cv. ‘Cabernet Sauvignon’ described above. Heat treatment was processed according to a previous study : leaf discs (5.5 cm in diameter) cut from the detached leaves were wrapped in a wet paper and placed in a small vessel made of aluminium foil. Leaf discs were sampled when vessels floated in a temperature-controlled water bath at 25 °C for 2 h and subsequently at 47 °C for 40 min. The control plants were treated in the same condition as that in heat treatment, except for temperature controlled at 25 °C.
To study the expression profile of VviGRFs in response to cold stress, six-week-old Vitis amurensis plantlets grown on 1/2 Murashige and Skoog (pH 5.8, 0.7% agar and 1% sucrose in 120 mL conical flasks) solid medium were placed in 4 °C (cold treatment) and 24 °C (control) chambers with a 16 h light/8 h dark cycle and 100 μmol m− 2 s− 1 light intensity. Shoots with the first fully expanded leaf were sampled after 0 h, 8 h and 12 h cold treatment with three biological replicates.
RNA extraction and qRT-PCR analysis
All samples were collected then preserved in liquid nitrogen immediately, and stored at − 80 °C for RNA and DNA extraction. Total RNA and DNA were extracted using the RNAprep Pure Plant Kit (DP432, Tiangen Biotech, Beijing, China) and Plant Genomic DNA Kit (DP305–02, Tiangen Biotech, Beijing, China) respectively, according to the manufacturer’s instructions. The cDNA used for verifying exon—intron structure was synthesised using HiScript® II 1st Strand cDNA Synthesis Kit (R211–01, Vazyme Biotech Co., Nanjing, China). To amplify the complete cDNA and DNA fragments corresponding to the VviGRF genes, specific primers for each gene were designed through NCBI Primer-BLAST (http://0-www.ncbi.nlm.nih.gov.brum.beds.ac.uk/tools/primer-blast/, Additional file 3), and the PCR reactions were performed using PrimeSTAR Max DNA Polymerase (R045A, Takara Biotechnology Co., Ltd., Dalian, China) with the following thermal cycling profile: 98 °C for 3 min, 34 cycles of 98 °C for 10 s and 60 °C for 5 s and 72 °C for 50 s, 72 °C for 5 min, 12 °C for ∞. PCR products were submit to Beijing Majorbio Company, and the sequencing results can be seen in Additional file 4. The cDNA used for qRT-PCR analysis was synthesised through reverse transcription of obtained total RNA using HiScript® II Reverse Transcriptase (R223–01, Vazyme Biotech Co., Nanjing, China). The synthesised cDNA was subjected to qRT-PCR with an Opticon thermocycler (CFX Connect Real-Time System; Bio-Rad, Hercules, CA) using SYBR Green PCR master mix (Vazyme, Nanjing, China) according to the manufacturer’s instructions. In qRT-PCR, VviActin (Accession: EC969944) was used as reference gene, and specific primers for VviGRF genes were designed by using Primer-BLAST (Additional file 5). The specificity of designed primers was verified through gel electrophoresis and sequencing. PCR reactions were as follows: 95 °C for 10 min, 40 cycles of 95 °C for 10 s and 60 °C for 30 s. Each sample was prepared in three biological and technical replicates. The relative expression levels of the VviGRFs were calculated using the 2−ΔΔCt method. Statistical difference was determined by t-test (**P < 0.01, *P < 0.05, n = 3) using R program .
Identification, multiple sequences alignment, and phylogenetic analysis of grape 14–3-3 family
Nineteen 14–3-3 proteins were obtained by using HMM and BLASTP. A total of 11 proteins containing 14–3-3 domain were identified in grape after searching in NCBI-CDD. The results of sequence alignment showed that VviGRFs contained several conserved domains in their aligned sequences, and nine α helices were identified in their secondary structures (Fig. 1). Based on phylogenetic results, these proteins were classified into two groups, namely, epsilon (VviGRF9a/9b/11/12) and non-epsilon groups (VviGRF14/15/16/17/18/−like1/−like2) (Figs. 2 and 3a). VviGRF proteins contained 64 (VviGRF-like2) to 582 (VviGRF-like1) amino acid residues, their Mol wts varied from 6.91 kDa (VviGRF-like2) to 66.03 kDa (VviGRF-like1), and their PIs ranged from 4.74 (VviGRF16) to 7.57 (VviGRF-like2) (Table 1).
Gene structure and motif analysis
Figure 3b illustrates the predicted motif and gene structure of VviGRF genes. And the exon—intron structure of VviGRF genes have been verified by sequencing, which was accordant with our predicted analysis (Additional file 4). The number of exon in VviGRFs varied from one (VviGRF-like2) to seven (VviGRF9b/ 11/ 12), and nine out of 11 VviGRFs showed not less than four exons (except for VviGRF-like1 and VviGRF-like2). A search on MEME program identified seven motifs in VviGRF proteins. The number of motifs in VviGRFs ranged from 1 to 7, and the length of motifs varied from 8 (motif 6) to 50 (motifs 1–3) amino acids. Most VviGRFs contained all seven motifs, except for VviGRF-like1 and VviGRF-like2. VviGRF-like1 contained five motifs without motifs 2 and 5. VviGRF-like2 only had motif 7. Moreover, NCBI-CDD results showed that the 14–3-3 domain of VviGRF-like2 was incomplete in N-terminal, VviGRF-like1 was incomplete in C-terminal. Furthermore, VviGRF-like1 contained a 14–3-3 domain, a domain of unknown function (DUF4283) and a zinc-binding motif (Additional file 6). Associated with gene and motif structure, the position of introns was conserved in grape 14–3-3 family (Fig. 3). In both groups, motif 1 and 6 were translated by a single exon, motif 2 and motif 5 were translated from sequences modelled by “exon-intron-exon”. In epsilon group, motifs 3, 4 and 7 in all four members were translated by “exon-intron-exon” sequences. In non-epsilon group, the sequence of 6 members (except for VviGRF-like2) encoding motif 7 also shown in this model. The detail information of introns position in “exon-intron-exon” sequences was shown in Additional file 7. Between epsilon group and non-epsilon group, the position of introns was absolutely different. But in the same group, the introns were present at the same position within their encoding sequences.
Chromosome distribution and synteny analysis of VviGRFs
Eleven VviGRF genes were distributed unevenly on seven chromosomes and chrUn in grape (Table 1 and Fig. 4). The chr18 contained three VviGRF genes (VviGRF11, VviGRF14 and VviGRF17), chr1, 2, 7, 10, 14 and 19 only contained one VviGRF gene. VviGRF15 and VviGRF-like1 were located on chrUn. We also conducted a synteny analysis on VviGRF genes to investigate the duplication event occurring in the grape 14–3-3 family (Fig. 5). Three duplication events were observed between chr10 and chr19 (VviGRF9a and VviGRF9b), chr18 and chrUn (VviGRF14 and VviGRF-like1), chr7 and chr18 (VviGRF17 and VviGRF18), which evolved from segment duplication.
Tissue-specific expression of VviGRF genes
To analyse the temporal and tissue-specific expression of VviGRFs in grapevine, we assessed their expression data in 54 tissues and organs obtained from the published GSE36128 (Fig. 6). Most VviGRFs exhibited a quite moderate expression in different tissues and organs, but some genes showed very high or very low expression in certain tissues. For instance, VviGRF12 was normally lowly expressed in most tissues and organs, but not in some floral organs such as stamen, bud, flower and pollen. And in normal conditions, VviGRF15 showed a high expression level in most tissues and organs, except for leaf, petal and pollen. VviGRF-like1 was generally lowly expressed in most tissues and organs, except for senescencing leaf, bud, berry pericarp and berry skin. VviGRF-like2 exhibited a high expression level in most tissues, but lowly expressed in well-developed inflorescence, young leaf, and woody stem.
Expression pattern of VviGRFs during berry development
To determine the VviGRF expression in different grape genotypes, we analysed their expression levels in five grape varieties during berry ripening (GSE62744, Fig. 7a). Except that VviGRFs in Sangiovese showed similar expression at soft and harvest stages, the expression patterns of VviGRFs in five varieties were similar at the early (pea and touch) and late (soft and harvest) stages of ripening. VviGRF9b, VviGRF11, VviGRF14 and VviGRF16 exhibited low expression level at the early stages and high expression level at the late stages. By contrast, VviGRF9a and VviGRF15 showed high and low expression levels at the early and late stages, respectively. The expression level of VviGRF12, VviGRF17 and VviGRF18 was constantly changed during berry ripening. At early stages, these three genes expressed highly at pea stage and lowly at touch stage, whereas at late stages, they expressed highly at soft stage and lowly at harvest stage. To study their expression patterns during grape berry development (Fig. 7b), eight VviGRFs in V. vinifera cv. ‘Cabernet Sauvignon’ were selected to qRT-PCR. Five out of 8 VviGRFs showed no significant expression level at three key stages during berry ripening (EL33, EL35 and EL37). Only VviGRF18 expression was up-regulated at both EL35 and EL37 stages. VviGRF14 and VviGRF15 were only up-regulated at EL37.
Response to cold and heat stresses
RNA-seq GSE89113 was used to investigate the response of VviGRFs to heat stress (Fig. 8a). Four VviGRFs (VviGRF9b/12/14/16) indicated up-regulated expression in detached leaves under 35 °C treatment compared with that at 25 °C, whereas the expression level of three VviGRFs (VviGRF11/15/17) showed down-regulation. Among them, VviGRF15 exhibited a significantly down-regulated expression (log2 fold change = − 0.969124). Under 40 °C, the expression levels of five members (VviGRF11/12/14/17/18) and three members (VviGRF9b/15/16) were down-regulated and up-regulated at 35 °C, respectively. Among them, VviGRF12 showed significantly decreased expression, but VviGRF15 significantly increased (|log2 fold change| > 1). Under 45 °C, all of nine VviGRFs showed an increase in expression compared with that at 40 °C. According to the data published in SRP018199, GSVIVT01012207001 and GSVIVT01009141001 (corresponding VviGRF12 and VviGRF14 in 12X V1), the V. vinifera cv. Muscat Hamburg seedlings showed up-regulated expression in response to the gradual cooling treatment (decreased at 5 °C per hour from 24 °C to 4 °C and subsequently 4 °C for additional 4 h, Additional file 8).
To obtain further understanding about the putative function of VviGRFs under heat and cold stresses, the expression patterns of nine VviGRFs in V. davidii and V. amurensis were analysed by qRT-PCR (Figs. 8b and 9). Three VviGRFs (VdGRF9b/14/18) showed no significant change in expression under heat treatment (25 °C 2 h and 47 °C 40 min), whereas the expression level of six VviGRFs (VdGRF9a/11/15/16/17/−like2) were up-regulated (Fig. 8b). After cold treatment (4 °C, 0–12 h, Fig. 9), except for three VviGRFs (VaGRF9a/11/16), six VviGRFs demonstrated significant expression change. VaGRF17 and VaGRF-like2 continuously decreased and increased during cold treatment, respectively. The expression of VaGRF9b exhibited no significant change at 8 h but increased at 12 h. VaGRF15 and VaGRF18 increased at 8 h and subsequently decreased at 12 h. VaGRF14 decreased at 8 h and subsequently increased at 12 h.
In this work, we identified and characterized the 14–3-3 gene family in V. vinifera. It is composed by 11 members, a number higher than in rice , but lower respect to Arabidopsis , tomato , cotton , banana , switchgrass , and soybean . Based on amino acid sequences, gene structure and phylogenetic relationship with A. thaliana, the 11 VviGRF proteins were classified into two groups (Table 1 and Fig. 2, Fig. 3a). This result is also in accordance with previous studies on other species [11, 12, 42]. For example, five and eight AtGRF proteins were contained in epsilon and non-epsilon groups in Arabidopsis, respectively.
A total of seven motifs were found in the grape 14–3-3 family. Epsilon group members contained all seven motifs, whereas non-epsilon group members commonly possessed less motifs than those of epsilon group members (e.g., VviGRF-like1 and VviGRF-like2). Similar to the motif analysis, epsilon group members commonly possessed more exons or introns than those of non-epsilon group, which also verified by RT-PCR (Table 1, Fig. 3b, and Additional file 4). Generally, only epsilon group members contained introns within sequences encoding motifs 3 and 4.This observation revealed that different groups of phylogenetic tree presented different motifs and exon—intron structures. This finding can also be observed in other species, such as Arabidopsis and rice; epsilon group members often present two extra introns in the N- or C-terminal compared with those of non-epsilon group [11, 42]. Take Arabidopsis for example (Additional file 1), there were a total of six motifs, and both their sequences and distribution were similar with those in grapevine. And as found in grapevine, sequences encoding motifs in Arabidopsis 14–3-3 family also possessed “exon-intron-exon” model. Moreover, Additional file 9 showed that the position of introns in this model was conserved in the same phylogenetic group, which was in accordance with grapevine. And between Vitis and Arabidopsis, introns within motifs 2 and 4 located in the same position. These similarities suggested different distribution of “exon-intron-exon” sequences encoding motifs between epsilon and non-epsilon groups, and also revealed the conservation of the exon—intron and motif structure in 14–3-3 family.
In grape 14–3-3 family, VviGRF-like1 was the longest protein, which displayed a low homology with other 14–3-3 members. And different from other VviGRFs, VviGRF-like1 contained an additional unknown function (DUF4283) and a zinc-binding motif. DUF4283 contains two highly conserved tryptophan residues and possibly provides binding/guiding region for diverse range of other domains (view pfam14111 in NCBI-CDD) . Moreover, NCBI-CDD results (Additional file 6) showed that the 14–3-3 domain of VviGRF-like2 lacked motif 1/3/4/6, and it was incomplete in N-terminal; the 14–3-3 domain of VviGRF-like1 lacked motif 2/5, and it was incomplete in C-terminal. These results indicated that motifs 1/3/4/6 and 2/5 may be the necessary components for the N- and C-terminals of 14–3-3 domain, respectively.
Gene duplication events are crucial in genomic rearrangement and often result in the expansion of gene family, which includes tandem, segment and transposition duplication . In grape 14–3-3 family, three pairs of genes evolved from segment duplication, but no tandem duplication event occurred (Fig. 5), indicating that the segment duplication may be the dominant gene duplication on the expansion of this family. And as for the remaining five VviGRFs, they may evolve in an early divergence time or be obtained from gene transposition.
To further understand the putative function of VviGRFs, the expression pattern of VviGRFs was searched for reference in the Vespucci database , the RNA-seq and expression-array data were also analysed. Fig. 6 reveals that VviGRFs were expressed in all investigated tissues and organs, and several VviGRFs showed tissue-specific expression in different tissues (e.g., VviGRF12/15/−like1/−like2). This phenomenon is also supported by some 14–3-3 genes in switchgrass, tomato, Arabidopsis and barley microspore embryogenesis [41, 45,46,47].
During berry ripening, VviGRFs in Sangiovese expressed similar pattern at soft and harvest stages, which may be caused by the varietal character or close sampling time at these two stages. Except this, VviGRFs in different varieties exhibited similar expression patterns between the early and late stages of berry ripening. In accordance with the results, VviGRF14, VviGRF15 and VviGRF18 displayed significant change in expression levels during EL33 (early stage), EL35, and EL37 stages (late stage). However, the up-regulated or down-regulated tendency of VviGRF15 differed in Fig. 7a and b, which may be impacted of variety difference. Taken together, some VviGRF genes presented changed expression during berry ripening, and they may be involved in the regulation of grape development. Accordingly, previous studies also demonstrated that plant 14–3-3 proteins are involved in fruit development and ripening processes. For example, most MaGRF transcriptions are significantly accumulated during fruit development and postharvest ripening in banana . Transcription profiles of switchgrass suggested that four PvGRFs may be involved in regulating lignin metabolism, and PvGRFr may participate in flower development . In Arabidopsis, RNA interference revealed the fundamental role of epsilon group members in regulating PIN polarity and plant development .
In the process of heating-up temperature (RNA-seq SRP091989, Fig. 8a), the expression of VdGRFs changed significantly, and nine genes presented up-regulated expression at 45 °C. In accordance with this result, six VdGRFs exhibited up-regulated expression in response to heat stress (Fig. 8b). Similar to this result, most 14–3-3 genes also display up-regulated expression under heat treatment in rice and Brachypodium distachyon L. [11, 49]. In V. amurensis, six VaGRF expression levels significantly changed in response to cold treatment, and five of six VaGRFs up-regulated or down-regulated after 8 h treatment, thereby revealing their potential roles in cold stress response (Fig. 9). In accordance with previous studies [11, 50, 51], two AtGRFs, five OsGRFs and two TaGRFs also indicated that cold stress induced significant expression changes in Arabidopsis, rice and wheat, respectively. VviGRF15, VviGRF-like2 and VviGRF17 showed significantly change in expression level under both cold and heat treatments, which indicated their potential roles in multiple-abiotic-stress responses. In summary, VviGRFs were important genes involved in cold and heat stress responses.
In grape, a total of 11 VviGRFs were identified and classified into epsilon and non-epsilon groups according to their phylogenetic relationship with A. thaliana. Protein and gene structure, together with gene duplication event analyses were performed to determine the characteristics and potential functions of grape 14–3-3 family. Published transcription data and qRT-PCR were analysed to gain further understanding regarding the putative function of VviGRF genes. The expression level of three VviGRFs showed significantly changes during berry ripening. Six VviGRFs expression altered under abiotic stress. The present results revealed the potential roles of VviGRF genes in berry development and abiotic stress response. This work provides new insight into grape 14–3-3 family and a basis for further study on its functional mechanism.
Conserved Domain Data
Gene Structure Display Server
Plasma membrane proton pumping
Hidden Markov Model
- Mol wt:
Fu H, Subramanian RR, Masters SC. 14-3-3 proteins: structure, function, and regulation. Annu Rev Pharmacol. 2000;40:617–47.
Chung H-J, Sehnke PC, Ferl RJ. The 14-3-3 proteins: cellular regulators of plant metabolism. Trends Plant Sci. 1999;4:367–71.
Aitken A, Collinge D, Van Heusden B, Isobe T, Roseboom P, Rosenfeld G, Soll J. 14-3-3 proteins: a highly conserved, widespread family of eukaryotic proteins. Trends Biochem Sci. 1992;17:498–501.
Bunney TD, De Boer AH, Levin M. Fusicoccin signaling reveals 14-3-3 protein function as a novel step in left-right patterning during amphibian embryogenesis. Development. 2003;130:4847–58.
van Heusden GPH, Wenzel TJ, Lagendijk EL, De Steensma H, van den Berg JA. Characterization of the yeast BMH1 gene encoding a putative protein homologous to mammalian protein kinase II activators and protein kinase C inhibitors. FEBS Lett. 1992;302:145–50.
Gökirmak T, Denison FC, Laughner BJ, Paul A-L, Ferl RJ. Phosphomimetic mutation of a conserved serine residue in Arabidopsis thaliana 14-3-3ω suggests a regulatory role of phosphorylation in dimerization and target interactions. Plant Physiol Bioch. 2015;97:296–303.
Ferl RJ. 14-3-3 proteins and signal transduction. Annu Rev Plant Biol. 1996;47:49–73.
Huber SC, MacKintosh C, Kaiser WM. Metabolic enzymes as targets for 14-3-3 proteins. Plant Mol Biol. 2002;50:1053–63.
Bai M-Y, Zhang L-Y, Gampala SS, Zhu S-W, Song W-Y, Chong K, Wang Z-Y. Functions of OsBZR1 and 14-3-3 proteins in brassinosteroid signaling in rice. P Natl Acad Sci USA. 2007;104:13839–44.
Diaz C, Kusano M, Sulpice R, Araki M, Redestig H, Saito K, Stitt M, Shin R. Determining novel functions of Arabidopsis 14-3-3 proteins in central metabolic processes. BMC Syst Biol. 2011;5:192.
Yao Y, Du Y, Jiang L, Liu J-Y. Molecular analysis and expression patterns of the 14-3-3 gene family from Oryza sativa. BMB Rep. 2007;40:349–57.
Xu WF, Shi WM. Expression profiling of the 14-3-3 gene family in response to salt stress and potassium and iron deficiencies in young tomato (Solanum lycopersicum) roots: analysis by real-time RT–PCR. Ann Bot-London. 2006;98:965–74.
Sun G, Xie F, Zhang B. Transcriptome-wide identification and stress properties of the 14-3-3 gene family in cotton (Gossypium hirsutum L.). Funct Integative Genomic. 2011;11:627–36.
Li M, Ren L, Xu B, Yang X, Xia Q, He P, Xiao S, Guo A, Hu W, Jin Z. Genome-wide identification, phylogeny, and expression analyses of the 14-3-3 family reveal their involvement in the development, ripening, and abiotic stress response in banana. Front Plant Sci. 2016;7:1442.
Li X, Dhaubhadel S. Soybean 14-3-3 gene family: identification and molecular characterization. Planta. 2011;233:569–82.
TAIR: The Arabidopsis Information Resource. http://www.arabidopsis.org/browse/genefamily/14-3-3.jsp. Accessed 20 Sept 2017.
Yang Z-P, Li H-L, Guo D, Tang X, Peng S-Q. Identification and characterization of the 14-3-3 gene family in Hevea brasiliensis. Plant Physiol Bioch. 2014;80:121–7.
Zhang Z-T, Zhou Y, Li Y, Shao S-Q, Li B-Y, Shi H-Y, Li X-B. Interactome analysis of the six cotton 14-3-3s that are preferentially expressed in fibres and involved in cell elongation. J Exp Bot. 2010;61:3331–44.
Sehnke PC, Chung H-J, Wu K, Ferl RJ. Regulation of starch accumulation by granule-associated plant 14-3-3 proteins. P Natl Acad Sci USA. 2001;98:765–70.
Yan J, He C, Wang J, Mao Z, Holaday SA, Allen RD, Zhang H. Overexpression of the Arabidopsis 14-3-3 protein GF14λ in cotton leads to a “stay-green” phenotype and improves stress tolerance under moderate drought conditions. Plant Cell Physiol. 2004;45:1007–14.
Myles S, Boyko AR, Owens CL, Brown PJ, Grassi F, Aradhya MK, Prins B, Reynolds A, Chia J-M, Ware D. Genetic structure and domestication history of the grape. P Natl Acad Sci USA. 2011;108:3530–5.
Farace G, Visconti S, Marchive C, Aducci P, Marra M, Trossat-Magnin C. Characterization, functional validation and gene expression patterns of two 14-3-3 isoforms from Vitis vinifera. Vitis. 2015;54:33–40.
Liu G-T, Ma L, Duan W, Wang B-C, Li J-H, Xu H-G, Yan X-Q, Yan B-F, Li S-H, Wang L-J. Differential proteomic analysis of grapevine leaves by iTRAQ reveals responses to heat stress and subsequent recovery. BMC Plant Biol. 2014;14:110.
Marchler-Bauer A, Bo Y, Han L, He J, Lanczycki CJ, Lu S, Chitsaz F, Derbyshire MK, Geer RC, Gonzales NR. CDD/SPARCLE: functional classification of proteins via subfamily domain architectures. Nucleic Acids Res. 2016;45:D200–3.
Larkin MA, Blackshields G, Brown N, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R. Clustal W and Clustal X version 2.0. Bioinformatics. 2007;23:2947–8.
Geourjon C, Deleage G. SOPMA: significant improvements in protein secondary structure prediction by consensus prediction from multiple alignments. Bioinformatics. 1995;11:681–4.
Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Molecular Biol Evol. 2016;33:1870–4.
Grimplet J, Adam-Blondon A-F, Bert P-F, Bitz O, Cantu D, Davies C, Delrot S, Pezzotti M, Rombauts S, Cramer GR. The grapevine gene nomenclature system. BMC Genomics. 2014;15:1077.
Hu B, Jin J, Guo A-Y, Zhang H, Luo J, Gao G. GSDS 2.0: an upgraded gene feature visualization server. Bioinformatics. 2014;31:1296–7.
Bailey TL, Boden M, Buske FA, Frith M, Grant CE, Clementi L, Ren J, Li WW, Noble WS. MEME SUITE: tools for motif discovery and searching. Nucleic Acids Res. 2009;37:W202–8.
Voorrips R. MapChart: software for the graphical presentation of linkage maps and QTLs. J Hered. 2002;93:77–8.
Wang Y, Tang H, DeBarry JD, Tan X, Li J, Wang X, Lee TH, Jin H, Marler B, Guo H. MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 2012;40:e49.
Krzywinski M, Schein J, Birol I, Connors J, Gascoyne R, Horsman D, Jones SJ, Marra MA. Circos: an information aesthetic for comparative genomics. Genome Res. 2009;19:1639–45.
Xin H, Zhu W, Wang L, Xiang Y, Fang L, Li J, Sun X, Wang N, Londo JP, Li S. Genome wide transcriptional profile analysis of Vitis amurensis and Vitis vinifera in response to cold stress. PLoS One. 2013;8:e58740.
Fasoli M, Dal Santo S, Zenoni S, Tornielli GB, Farina L, Zamboni A, Porceddu A, Venturini L, Bicego M, Murino V. The grapevine expression atlas reveals a deep transcriptome shift driving the entire plant into a maturation program. Plant Cell. 2012;24:3489–505.
Palumbo MC, Zenoni S, Fasoli M, Massonnet M, Farina L, Castiglione F, Pezzotti M, Paci P. Integrated network analysis identifies fight-club nodes as a class of hubs encompassing key putative switch genes that induce major transcriptome reprogramming during grapevine development. Plant Cell. 2014;26:4617–35.
Jiang J, Liu X, Liu G, Liu C, Li S, Wang L. Integrating omics and alternative splicing reveals insights into grape response to high temperature. Plant Physiol. 2017;173:1502–18.
Coombe B. Growth stages of the grapevine: adoption of a system for identifying grapevine growth stages. Aust J Grape Wine R. 1995;1:104–10.
Xu H, Liu G, Liu G, Yan B, Duan W, Wang L, Li S. Comparison of investigation methods of heat injury in grapevine (Vitis) and assessment to heat tolerance in different cultivars and species. BMC Plant Biol. 2014;14:156.
Team RC. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2014. p. 2014.
Wu S, Yan H, Zhang A, Huang L, Yin G, Lee S. Identification and characterization of the 14-3-3 gene family in switchgrass. Genet Mol Res. 2016;15:15048688.
DeLille JM, Sehnke PC, Ferl RJ. The Arabidopsis 14-3-3 family of signaling regulators. Plant Physiol. 2001;12:35–8.
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:10.
Moretto M, Sonego P, Pilati S, Malacarne G, Costantini L, Grzeskowiak L, Bagagli G, Grando MS, Moser C, Engelen K. VESPUCCI: exploring patterns of gene expression in grapevine. Front Plant Sci. 2016;7:633.
Maraschin S, Louwerse JD, Lamers GE, Spaink HP, Wang M. Tissue-specific expression of 14-3-3 isoforms during barley microspore and zygotic embryogenesis. Acta Biol Cracov Ser Bot. 2003;45:103–6.
Wei-Feng X, Wei-Ming S, Feng Y. Temporal and tissue-specific expression of tomato 14-3-3 gene family in response to phosphorus deficiency. Pedosphere. 2012;22:735–45.
Daugherty CJ, Rooney MF, Miller PW, Ferl RJ. Molecular organization and tissue-specific expression of an Arabidopsis 14-3-3 gene. Plant Cell. 1996;8:1239–48.
Keicher J, Jaspert N, Weckermann K, Möller C, Throm C, Kintzi A, Oecking C. Arabidopsis 14-3-3 epsilon members contribute to polarity of PIN auxin carrier and auxin transport-related development. eLife. 2017;6:e24336.
Cao H, Xu Y, Yuan L, Bian Y, Wang L, Zhen S, Hu Y, Yan Y. Molecular characterization of the 14-3-3 gene family in Brachypodium distachyon L. reveals high evolutionary conservation and diverse responses to abiotic stresses. Front Plant Sci. 2016;7:1099.
Leyva A, Jarillo JA, Salinas J, Martinez-Zapater JM. Low temperature induces the accumulation of phenylalanine ammonia-lyase and chalcone synthase mRNAs of Arabidopsis thaliana in a light-dependent manner. Plant Physiol. 1995;108:39–46.
Meng X, Chen X, Wang Y, Xiao R, Liu H, Wang X, Ren J, Li Y, Niu H, Wang X, et al. Characterization and subcellular localization of two 14-3-3 genes and their response to abiotic stress in wheat. Sheng Wu Gong Cheng Xue Bao. 2014;30:232–46.
Average C. The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla. Nature. 2007;449:463–7.
We thank Erpeng Zhang, Yingying Gao and Fengrui Ren for revising the manuscript, and Renkun Tang for help with the R software.
This work was supported by the grants from the National Science Foundation of China (31572090 and 31772266), Agricultural Breeding Project of Ningxia Hui Autonomous Region (NXNYYZ20150203) and Hundred Talents of Chinese Academy of Sciences. The design of this study and the manuscript writing were supported by NXNYYZ20150203 and Hundred Talents of Chinese Academy of Sciences. And the collection, analysis, and interpretation of data used in this study were supported by 31572090 and 31772266.
Availability of data and materials
The datasets supporting the conclusions of this article are included within the article and its Additional files.
Ethics approval and consent to participate
Vitis vinifera cv. ‘Cabernet Sauvignon’ is the famous wine grape widely cultivated in the world. V. davidii and V. amurensis are the important wild plant species, and they are famous for their cold or heat resistance, respectively. These three cultivars are not listed in the appendices I, II and III of the Convention on the Trade in Endangered Species of Wild Fauna and Flora, that has been valid from 4 October 2017 (https://cites.org/eng/app/appendices.php). The plant samples used in this study were collected from Beijing Key Laboratory of Grape Sciences and Enology, Laboratory of Plant Resources, Institute of Botany, the Chinese Academy of Sciences, Beijing, China. Collection of plant materials complied with the institutional, national and international guidelines. No specific permits were required.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Phylogenetic relationship, motif and gene structure analysis of AtGRF proteins. (A) Phylogenetic tree of 14–3-3 family in Arabidopsis. (B) Gene structure and motif analysis of 14–3-3 members in Arabidopsis. The black blocks represent the untranslated region (UTR), the gray blocks represent exons, and the black lines represent introns. (C) The amino acid sequences of motifs in AtGRF proteins. Arrows showed the intron position appeared in “exon-intron-exon” sequences. Red arrow means the epsilon group, green arrow means the non-epsilon group. The detail information for the intron position was shown in Additional file 9. (TIF 2369 kb)
The primer sequence for verifying exon—intron structure of seven VviGRF genes. (DOC 30 kb)
Comparison of sequencing results and exon—intron structure of seven VviGRF genes. (A) cDNA Sequence alignment of VviGRFs. ‘_sequencing’ means cDNA of each VviGRF cloned from green stem in Cabernet Sauvignon. ‘_cDNA’ means reference sequences of each VviGRF cDNA in Point Noir. (B) DNA Sequence alignment of VviGRFs. ‘_sequencing’ means DNA of each VviGRF cloned from green stem in Cabernet Sauvignon. ‘_DNA’ means reference sequences of each VviGRF DNA in Point Noir. Sequence alignment analysis was performed by using DNAMAN version7. Sequences of seven VviGRFs in Cabernet Sauvignon were homologous with their reference sequence in Point Noir (Identity > 94%). (C) The exon—intron structure of seven VviGRFs. CS means sequencing gene structure in Cabernet Sauvignon, PN means reference gene structure in Point Noir. (DOC 18109 kb)
The primer sequence for qRT-PCR. (DOC 31 kb)
The conserved domain analysis of VviGRF proteins. 14–3-3 domain was present in blue strips, and the jagged edges means incomplete N- or C-terminal. (TIF 679 kb)
Motif sequences of VviGRF proteins. Black arrow shows the introns position appeared in “exon-intron-exon” sequences. (TIF 5120 kb)
Expression profile of 14–3-3 genes in V. vinifera cv. Muscat of Hamburg under cold stress from published data SRP018199. (DOC 28 kb)
Motif sequences of AtGRF proteins. Black arrow shows the introns position appeared in “exon-intron-exon” sequences. (TIF 6446 kb)
Detail information of grape 14–3-3 family. (DOC 39 kb)
About this article
Cite this article
Cheng, C., Wang, Y., Chai, F. et al. Genome-wide identification and characterization of the 14–3-3 family in Vitis vinifera L. during berry development and cold- and heat-stress response. BMC Genomics 19, 579 (2018) doi:10.1186/s12864-018-4955-8
- Abiotic stress