- Research article
- Open Access
Genome-wide identification and analysis of B-BOX gene family in grapevine reveal its potential functions in berry development
BMC Plant Biology volume 20, Article number: 72 (2020)
The B-BOX (BBX) proteins are the class of zinc-finger transcription factors and can regulate plant growth, development, and endure stress response. In plants, the BBX gene family has been identified in Arabidopsis, rice, and tomato. However, no systematic analysis of BBX genes has been undertaken in grapevine.
In this study, 24 grapevine BBX (VvBBX) genes were identified by comprehensive bioinformatics analysis. Subsequently, the chromosomal localizations, gene structure, conserved domains, phylogenetic relationship, gene duplication, and cis-acting elements were analyzed. Phylogenetic analysis divided VvBBX genes into five subgroups. Numerous cis-acting elements related to plant development, hormone and/or stress responses were identified in the promoter of the VvBBX genes. The tissue-specific expressional dynamics of VvBBX genes demonstrated that VvBBXs might play important role in plant growth and development. The transcript analysis from transcriptome data and qRT-PCR inferred that 11 VvBBX genes were down-regulated in different fruit developmental stages, while three VvBBX genes were up-regulated. It is also speculated that VvBBX genes might be involved in multiple hormone signaling (ABA, ethylene, GA3, and CPPU) as transcriptional regulators to modulate berry development and ripening. VvBBX22 seems to be responsive to multiple hormone signaling, including ABA, ethylene GA3, and CPPU. Some VvBBX genes were strongly induced by Cu, salt, waterlogging, and drought stress treatment. Furthermore, the expression of VvBBX22 proposed its involvement in multiple functions, including leaf senescence, abiotic stress responses, fruit development, and hormone response.
Our results will provide the reference for functional studies of BBX gene family, and highlight its functions in grapevine berry development and ripening. The results will help us to better understand the complexity of the BBX gene family in abiotic stress tolerance and provide valuable information for future functional characterization of specific genes in grapevine.
Zinc-finger transcription factor is one of the most important family in the plant kingdom, which play essential role in plant growth, development, and response to environmental stimuli [1, 2]. Zinc finger transcription factors are divided further into several subfamilies according to the structural and functional features of their individual members. The BBX proteins belong to zinc-finger transcription factors and attracted more attention in recent years due to their multiple functions. BBX proteins in plants are comprised of one or two conserved BBX domains in the N-terminus and occasionally a CCT (CONSTANS, CO-like, and TOC1) domain in the C-terminus. B-BOX motifs play significant role in protein-protein interactions and transcriptional regulation [3, 4]. The CCT domain participates in nuclear transport and transcriptional regulation [5,6,7]. In Arabidopsis, 32 BBX proteins have been identified and these members are classified into five subgroups depending on the presence of B-BOX domains along with the CCT domain . Subsequently, increasing evidence suggests that plant BBX proteins play pivotal role in diverse physiological and biochemical processes, such as flower induction [8, 9], photomorphogenesis [10, 11], shade avoidance response , carotenoid biosynthesis , and biotic and abiotic stress response [14, 15].
CONSTANS (CO)/AtBBX1 is the first investigated BBX gene in Arabidopsis, which is a central coordinator in controlling flowering time by triggering the expression of Flowering Locus T (FT) gene [16,17,18]. co mutants delay flowering significantly under a long-day, whereas co-overexpressed transgenic plants flower early in both long and short-day conditions [19,20,21]. There are other BBX genes, such as BBX4, BBX7, and BBX32, which have been found to regulate the flowering times [22,23,24]. At least ten BBX genes have been identified as a regulators of the early photomorphogenesis. Out of ten, BBX4, BBX21, BBX22, and BBX23 positively regulate plant photomorphogenesis [22, 25,26,27], whereas BBX19, BBX20, BBX24, BBX25, BBX28 and BBX32 are negative regulator of photomorphogenesis in Arabidopsis [11, 28,29,30]. For example, BBX28 negatively regulates photomorphogenesis by repressing the expression of HY5 and undergoes COP1-mediated degradation . Several BBX genes also show their functions in shade avoidance by mediating cell elongation [12, 31, 32].
The BBX proteins are also involved in abiotic stress response and hormonal signaling networks. For example, BBX18 and BBX23 are positive thermomorphogenesis regulators and the deficient mutations of BBX18 and BBX23 result in reduced thermoresponsive hypocotyl elongation . In Arabidopsis, BBX24 is initially isolated as a salt tolerant protein (STO) and increases salt tolerance activities in yeast cells . Overexpression of BBX24 in Arabidopsis enhances the growth of roots under high salinity conditions . In Chrysanthemum, CmBBX24 also enhances cold and drought tolerance besides delaying flowering time . BBX genes also reveal their functions in phytohormone signal transduction. AtBBX18 (AtDBB1a) is a positive regulator in the gibberellin (GA) signaling pathway , whereas BBX20 (AtBZS1) negatively regulates the brassinosteroid signaling network . Interestingly, BBX proteins also show their importance in the development of fruits, especially in anthocyanin and carotenoid biosynthesis. MdBBX22 (MdCOL11) is involved in MdHY5-mediated signal transduction and modulates anthocyanin accumulation in apple peel . SlBBX20 promotes chloroplast development and carotenoid accumulation in tomato by directly activating the expression PHYTOENE SYNTHASE 1 .
Grapevine is one of the most widely cultivated and commonly consumed fruit crops throughout the world due to its economic importance and essential nutrition [37, 38]. Although the BBX family has been identified in Arabidopsis, tomato, pear, and apple [3, 39,40,41], yet no comprehensive study of BBX genes in grapevine has been reported so far. With the release of the grapevine genome , we have a better possibility to systematically investigate the putative functions of BBX genes in grapevine. In this study, 24 non-redundant members of the VvBBX gene family were characterized in grapevine. Subsequently, the detailed gene structures, phylogenetic relationships, tissue expression profiles and expression profiles under different stress conditions were investigated. Our results of VvBBX genes will provide a foundation for further functional characterization of BBX genes in grapevine.
Identification of VvBBX genes in grapevine
To identify and obtain the BBX genes in the grapevine genome, the Arabidopsis BBX proteins were used as a query to search against the local grapevine genome database by using DNAtools software. Then, the hidden Markov model (HMM) profile of the B-box domain (Pfam00643) was employed to perform a global search of the grapevine genome. After analyzing the conserved domain and removing the redundant sequences, a total of 24 putative VvBBX genes were identified in grapevine. For the sake of nomenclature and consistency, these VvBBX genes were named from VvBBX1 to VvBBX24 depending on their homology to the pear BBX members, based on their similar numbers of BBX family between grapevine and pear . The detailed information of VvBBX was listed in Table 1, including gene name, protein length, chromosome location, molecular weight, theoretical isoelectric point, aliphatic index, and GRAVY. The 24 VvBBX proteins had diverse molecular length and weight, ranging from 127 (VvBBX22) to 474 (VvBBX7) in amino acid length. VvBBX22 showed the lowest value of the molecular weight (14.26 kDa), while the highest of the molecular weight (51.43 kDa) was observed in VvBBX7. Theoretical isoelectric points of these VvBBX proteins varied from 4.29 (VvBBX21) to 9.01 (VvBBX24) and the value of the aliphatic index ranged from 47.66 (VvBBX23) to 78.38 (VvBBX24). The GRAVY of all VvBBXs was less than zero, indicating the hydrophilic nature of VvBBX proteins (Table 1). The majority of VvBBX proteins were predicted to be located on the nucleus by WoLF PSORT, but a few of them may be located in other subcellular compartments, such as chloroplast and cytoplasm (Table 1).
Protein sequence and phylogenetic analysis of the VvBBX gene family
The length of VvBBX proteins varied widely from 127 to 474 amino acids. Among them, eight VvBBXs were found to contain two B-BOX domains and a conserved CCT domain. Nine members consisted of two B-BOX domains but no CCT domain. Two VvBBX contained one B-BOX plus a CCT domain, and five with only one B-BOX domain (Additional file 1: Figure S1). The protein sequence alignment and motifs logos showed that the B-Box 1 and B-Box 2 domains of the VvBBXs had similar conserved amino acid residues, and the CCT domain was also highly conserved among the VvBBX proteins (Additional file 2: Figure S2 and Additional file 3: Figure S3). The motifs logos of these domains were shown in Additional file 2: Figure S2, and their correspondence locations were illustrated in Additional file 3: Figure S3.
To explore the evolutionary relationship and functional divergence of the VvBBX members, full-length amino acid sequences of VvBBXs and pear BBXs (PbBBXs) was used to construct the phylogenetic tree by using the neighbor-joining and maximum likelihood method in MEGA, respectively (Fig. 1a and Additional file 4: Figure S4). As shown in the phylogenetic tree (Fig. 1a and Additional file 4: Figure S4), the VvBBX family was divided into five subgroups according to the phylogenetic analysis and consistent with the previous study in pear, tomato, and Arabidopsis [3, 39, 40]. Additionally, to gain a better understanding of the classification of BBX members, the phylogenetic tree of VvBBXs, together with PbBBXs, AtBBXs, and SlBBXs was also constructed by using the maximum likelihood method (Additional file 5: Figure S5). As shown in the phylogenetic tree (Additional file 5: Figure S5), all BBX proteins were also divided into five subgroups and most BBX members from grapevine and pear clustered together. Furthermore, the sequences of B-Box 1 (Fig. 1b), B-Box 2 (Fig. 1c) and CCT (Fig. 1d) domain were also employed for the phylogenetic analysis, respectively. The members of subgroups I, II and III contained both B-BOX and CCT domains except for VvBBX9 and VvBBX10, which only owned two B-BOX domains. On the contrary, the members in subgroup IV and V had no CCT domain and only contained two or one B-BOX domain(s), respectively. Eight VvBBXs in subgroup I and II contained all the three domains, and they could be divided into two different groups (Fig. 1b, c, and d). Three of them (VvBBX1, 2, 3) were clustered together in subgroup I, while the other five members from subgroup II also aligned together.
Chromosomal localization, gene duplication and gene structure analysis of VvBBXs
To understand the genomic distribution and gene duplication of VvBBX genes, 22 VvBBX genes were distributed unevenly throughout the 11 out of the 19 chromosomes and the remaining two genes (VvBBX5 and VvBBX23) had not yet been assembled to any chromosome according to the current grapevine genome (Fig. 2). Among them, the chromosomes 12 had the highest number of VvBBX genes (four), while only one VvBBX gene was localized on chromosome 5, 7, 9 and 11. Three VvBBX genes were located on each of chromosomes 1 and 19 and two VvBBX genes were distributed on chromosome 3, 4, 14 and 18, respectively (Fig. 2).
Tandem duplication and segmental duplications occurred frequently in gene families evolution and expansion. Tandem duplication usually caused gene clusters and segmental duplication might lead to scattered family members . Only one tandem duplication cluster (VvBBX11/VvBBX22) in VvBBX gene family was identified on grapevine chromosome 12. Then, five pairs of duplicated segments (VvBBX7/VvBBX8, VvBBX10/VvBBX11, VvBBX13/VvBBX14, VvBBX14/VvBBX15, and VvBBX15/VvBBX16) in VvBBX gene family were identified within the grapevine genome (Fig. 2). The result suggested that segmental duplication events may be more important than tandem duplication in the expansion of the VvBBX gene family in grapevine. Furthermore, we also calculated the value of Ka/Ks of segmental genes’ pairs, which could be used as an indicator for the selection pressure of a gene during evolution. Our results showed that all the Ka/Ks values were less than 1, indicating that the VvBBX genes primarily evolved under the influence of purifying selection (Additional file 6: Table S1).
To better insight into the evolutionary relationships of the VvBBX genes, the intron/exon structures were investigated by aligning the cDNA sequences and corresponding genomic DNA sequences. As shown in Additional file 7: Figure S6, the number of exons varied from 1 to 5. VvBBX13 and VvBBX14 contained the highest amounts of exons (5), and VvBBX22 and VvBBX24 only had one exon among the VvBBXs. Additionally, the highly similar gene structure was identified in the same group of the VvBBX genes. For example, all three VvBBXs in the group I contained two exons, and most members within the group II had four exons except for VvBBX9 (Additional file 7: Figure S6). These results suggested that exon-gain or -loss had occurred during the evolution of the VvBBX gene family and gene structure reveal the evolutionary relationship of VvBBXs.
Identification of cis-elements in the promoters of VvBBX genes
To better understand the transcriptional regulation and the gene function of VvBBX, the cis-elements in the promoter regions of the VvBBX (2 kb of genomic DNA sequence upstream of the translation start site) were used to search the PlantCARE database (Additional file 8: Figure S7, Additional file 9: Table S2). As expected, CAAT-box and TATA-box, the conventional promoter elements, were found in all the VvBBX promoters. A series of cis-elements involved in plant growth and development, phytohormone responses and stress responses were identified (Additional file 8: Figure S7, Additional file 9: Table S2). As shown in Fig. 4, the CAT-box involved in meristem expression was identified in the promoters of 12 VvBBX genes. The zein metabolism regulation element (O2 site) was found in 10 VvBBX genes. Additionally, the seed-specific regulation element (RY element) and endosperm expression regulation element (GCN4_motif) were also found in the promoters of the VvBBX genes. Among the cis-acting elements involved in hormone responses, the abscisic acid responsive element (ABRE), ethylene responsive element (ERE), the MeJA responsive element (CGTCA-motif and TGACG-motif) and gibberellin-responsive element (GARE-motif, P-box and TATC-box) were found in the promoters of 21, 18, 17 and 12 VvBBX genes, respectively. Salicylic acid responsive element (TCA-element) and auxin-responsive element (AuxRR-core and TGA-element) were also observed in 9 and 5 VvBBX genes (Additional file 8: Figure S7, Additional file 9: Table S2). In stress-related responses elements, ARE, which was the most abundant element and involved in anaerobic induction, was found VvBBX genes. Some stresses-related (low-temperature, drought and wound) cis-acting elements were also found in the promoter regions of the VvBBX genes (Additional file 8: Figure S7, Additional file 9: Table S2). In addition, a number of light response elements were also found in the promoter regions of VvBBX, including Box 4, G-Box, G-box, GATA-motif, GT1-motif, I-box, MRE, TCCC-motif and TCT-motif.
The organ-specific expression patterns analysis of VvBBX genes in grapevine
In order to investigate the putative roles of the VvBBX genes in grapevine development, the organic-specific expression patterns of VvBBXs were analyzed in the V. vinifera cv. Corvina global gene expression atlas from the GEO DataSets (GSE36128), which contained 42 various organs/tissues at different developmental stages obtained by microarray analysis (Fig. 3, Additional file 10: Table S3). As shown in Fig. 3, some VvBBX genes exhibited similar expression profiles in different organs/tissues, while other VvBBXs showed tissue-specific transcript accumulation patterns, potentially suggesting the functional divergence of VvBBX genes during grapevine growth and development. For example, four VvBBX (VvBBX2, 3, 5 and 6) genes were ubiquitously high expressed in nearly all tissues tested, whereas VvBBX16 and VvBBX17 were expressed at a very low level in all tested tissues.
Some VvBBX genes showed a very high level in specific organs/tissues. For example, three VvBBX (VvBBX4, 21 and 24) genes displayed higher expression level in pollen than other tissues, implying that they were important proteins for signaling in grapevine pollen growth and development. VvBBX1 and VvBBX22 exhibited high levels of expression in senescing leaf, which suggested that they might be involved in leaf senescing. VvBBX20 showed relatively high expression level in woody stem (Stem-W), indicating involvement in lignification of xylem and phloem during secondary growth of woody tissues. VvBBX15 was preferentially expressed in seedling, root, and seed (seed V), which suggested VvBBX15 played important roles in three tissues. Furthermore, two pair VvBBX in the segmental duplications showed similar expression patterns. For example, VvBBX7 and VvBBX8 with similar gene structure and common motif composition had low expression abundances in berries. VvBBX10 and VvBBX11 showed relatively high expression level in the stamen. (Fig. 3, Additional file 10: Table S3). Remarkably, VvBBX6, VvBBX14, and VvBBX23 were highly expressed in berries, which indicated that these genes may play an important role in berry development and ripening. These results aroused us to investigate the expression patterns of VvBBX genes during different fruit development and ripening stages.
Expression patterns of VvBBX genes during berry developmental and ripening
To understand the potential function of VvBBX genes in berry development and ripening, the transcript expression patterns of 24 VvBBX genes were investigated during three fruit developmental stages in grapevine using the expression profiles from the NCBI Gene Expression Omnibus (GEO) DataSets (GSE77218). Different members of the VvBBX genes showed distinct transcript accumulation patterns during different fruit developmental stages. As shown in Fig. 4a, VvBBX3, VvBBX6, and VvBBX14 showed relatively high expression levels during the ripening process, whereas six VvBBX genes (VvBBX5, 10, 11, 15, 16 and 17) were almost undetectable during different berry developmental stages (Fig. 4a, Additional file 11: Table S4). 11 VvBBX genes were down-regulated expression patterns, while three VvBBX genes displayed up-regulated expressions during fruit developmental (Fig. 4a, Additional file 11: Table S4). For example, the expression level of VvBBX1 was moderate and gradually decreased from veraison till to ripe stage. Similar to VvBBX1, but to a lesser extent, VvBBX18, VvBBX19, and VvBBX20 had the highest expression at green fruit stage. VvBBX22 and VvBBX23 exhibited a significant increase in expression during berry development and reached their picks at the ripening stage (Fig. 4a, Additional file 11: Table S4).
In order to validate the transcript abundance of VvBBX genes in the various developmental stages of the berry by microarray data, qRT-PCR analysis of all 18 detectable VvBBX genes was further performed at three berry development stages. As was expected, qRT-PCR results were highly consistent with the RNA-Seq data except for VvBBX13 and VvBBX24 (Fig. 4b). For example, the expression of VvBBX1 and VvBBX23 were dramatically decreased and increased during berry developmental and ripening, respectively. VvBBX18, VvBBX19, and VvBBX20 also depicted the highest expression at the green fruit stage (Fig. 4b). All these results were consistent with the data from RNA-Seq data. However, the expression profiles of VvBBX13 and VvBBX24 did not correspond with RNA-Seq data. VvBBX13 were highly expressed in ripe stage from RNA-Seq data, whereas the qRT-PCR result showed the highest expression in the veraison stage. VvBBX24 had a relatively high expression in veraison berry from qRT-PCR analysis, whereas the RNA-Seq data showed the highest expression in the green berry stage (Fig. 4b). All these results implied that some VvBBX genes showed different expression patterns and might play essential roles during fruit development and ripening.
To provide more information on the berry developmental and ripening functions of VvBBX genes in grapevine, we investigated their transcript accumulation patterns among 10 different grapevine varieties by using microarray data from the NCBI GEO DataSets (GSE62744 and GSE62745), which consists of four different fruit developmental stages (the pea-sized berry stage at 20d after flowering, the berries beginning to touch stage just prior to veraison, the berry-softening stage at the end of veraison, and the fully ripe berry stage at harvest . As shown in Fig. 5, the transcript accumulation profiles of VvBBX3 and VvBBX6 remained relatively strong but decreased gradually throughout grapevine fruit ripening, which was corresponded with the data from RNA-Seq data. VvBBX9 and VvBBX22 were preferentially expressed in End_veraison and ripe stage, and the similar expression pattern was also observed from RNA-Seq data. On the contrary, VvBBX1 and VvBBX19 were higher expression level in pea-sized berry and Pre_veraison stage and rapidly down-regulated during ripening, which was agreed with the data from RNA-Seq and qRT-PCR analysis (Fig. 4). Additionally, 12 VvBBX genes exhibited slight or no expression in fruit among all 10 grapevine varieties. All these results indicated that some VvBBX genes maybe play multiple important roles in grapevine fruit development.
VvBBX genes in response to exogenous ABA, ethylene, GA3 and CPPU hormones
Plant hormones are originally characterized as regulators in growing and developmental processes, simultaneously, the role of BBX proteins in hormonal signaling pathways is scarce in grapevine. To reveal the potential roles of the VvBBX genes in response to (abscisic acid) ABA and ethylene treatments, qRT-PCR was used to analyze the expression of the VvBBX genes, which was relatively high expression during grapevine berry ripening. As shown in Fig. 6, the expression of VvBBX1 was significantly inhibited by pre-veraison ABA and ethylene treatment throughout the entire berry development and ripening period, indicating that VvBBX1 play a negative role in the regulation of ripening (Fig. 6). Similarity, VvBBX2, VvBBX3 and VvBBX13 seemed to be suppressed at ripening stage, implying that these two genes in grapevine might be act as a negative regulator of fruit ripening. VvBBX6 and VvBBX9 were slight increase by ethylene and ABA treatment, respectively. VvBBX12 was up-regulated at green and veraison stage and VvBBX14 was down-regulated at green and ripe stage after exogenous ethylene treatment (Fig. 6), indicating that VvBBX12 and VvBBX14 might have opposite function in response to ethylene signaling pathway. It is worth noting that the transcript level of VvBBX22 was significantly upregulated at green and veraison stage by ABA and ethylene treatment (Fig. 6), indicating VvBBX22 might play a positive role in response to ABA and ethylene treatments.
Currently, GA3 and forchlorfenuron (CPPU) commonly applied to increase fruit weights, produce seedless grapes and inhibit russet development . Therefore, we further investigated the potential roles of the VvBBX genes in response to GA3 and CPPU by previous RNA-seq data . As shown in Additional file 12: Table S5, the expression of VvBBX7, VvBBX8 and VvBBX12 was significantly increased, whereas VvBBX21, VvBBX22 and VvBBX24 was significantly decreased by GA3 and CPPU treatment (Additional file 12: Table S5).
Expression patterns of VvBBX genes under different abiotic stresses
Copper (Cu), salt, waterlogging and drought are common abiotic stresses of vineyards. To understand the potential functions of the VvBBX genes in responses to these different environmental stresses, we collected RNA-seq data sets for grapevines subjected to the above four stresses (drought stress, 20 days, SRP074162, waterlogging stress, 48 h, SRP070475, salt, 48 h, see Additional file 2 in Int. J. Mol. Sci. 2018, 19, 4019, copper stress, 24 h, see Additional file 9: Table S2 in Sci. Rep. 2015, 5, 17,749). Overall, the VvBBX genes responded to waterlogging and drought stress to a greater extent than to Cu and NaCl treatment (Fig. 7, Additional file 13: Table S6). Among them, two VvBBXs (VvBBX22 and VvBBX23) were significantly induced against the Cu-treated sample and two VvBBXs (VvBBX1 and VvBBX14) were down-regulated by Cu stress (Fig. 7, Additional file 13: Table S6). Under salt stress, four VvBBX genes were up-regulated and only VvBBX8 was down-regulated. In contrast, 10 VvBBX genes were identified as differentially expressed genes under waterlogging stress, of which, only VvBBX14 and VvBBX19 were down-regulated, and the remaining genes were up-regulated (Fig. 7, Additional file 13: Table S6). In response to drought stress, nine VvBBX genes showed an increased expression pattern and six VvBBXs were more or less reduced (Fig. 7, Additional file 13: Table S6). Notably, four VvBBXs (VvBBX1, 14, 15 and 16) responded to at least three treatments. Interestingly, VvBBX22 was strongly induced at least 7-fold by all four types of stresses, suggesting that VvBBX22 might be a candidate gene for mitigating abiotic stresses. Similar to VvBBX22, but to a lesser extent, VvBBX23 was also induced dramatically by various environmental stresses (Fig. 7, Additional file 13: Table S6).
Subcellular localization of grapevine BBX proteins
The nuclear localization of transcription factors is very important for its regulatory function. Previous studies showed that BBX proteins were predominantly located on the nucleus, such as PbBBX4, PbBBX5, and PbBBX13 in pear . The subcellular localization of VvBBX proteins was firstly predicted by WoLF PSORT (Table 1). Seventeen of them had a high probability to be located in the nucleus. Four genes (VvBBX1, 9, 13 and 22) which were specifically expressed in fruit, were selected for a transient expression assay in Nicotiana benthamiana epidermal cells. As indicated in Fig. 8, green fluorescence signals from the expressed fusion VvBBX1-GFP, VvBBX9-GFP, VvBBX13-GFP, and VvBBX22-GFP were specifically distributed within the nucleus as confirmed by a mCherry-labelled nuclear marker (NF-YA4-mCherry). Notably, VvBBX13-GFP were also localized in the cytoplasm (Fig. 8), which was inconsistent with the prediction result (Table 1). These results suggested that VvBBX1, VvBBX9, and VvBBX22 were nuclear proteins, and consistent with the previous results [39, 40].
The BBX protein is one of the important transcription factors, which play an important role in regulating plant growth and development . To date, the features and functions BBX gene family have been identified in several plant species, such as Arabidopsis , rice , tomato , pear  and apple . Although a comprehensive analysis of the BBX gene family in grapevine during growth, berry development, and stress response has not been studied so far. In this study, we performed a genome-wide analysis of the BBX genes in grapevine by investigating their linkage group organization, phylogenetic relationships, gene structure, duplication events, cis-acting elements, gene expression profiles in different tissues and developmental stages, and under various stress treatments. Genome-wide analysis of the BBX genes in grapevine will set a foundation for further functional studies of this gene family for the molecular cloning of grapevine.
Evolution of the grapevine BBX gene family
The presence of BBX genes in the genome of different species from algae to monocots and dicots clearly suggests an ancient origin [4, 47]. A total of 24 VvBBX genes were identified from the grapevine genome, and the number of VvBBX genes was much smaller than that from other plants, such as 64 for Malus domestica , 32 for Arabidopsis , 30 for rice , and 29 for tomato . The difference might be due to the variable state of paralogous genes in these genomes. For example, only 5 paralogous gene pairs were found in grapevine, while 22 paralogous and 9 paralogous gene pairs from segmental duplication events were identified in apple and rice, respectively. This conclusion was further supported by the previous study in D. officinale (19 members) and P. equestris (16 members), including only two DoBBXs paralogous and two PeBBXs paralogous gene pairs, respectively . Furthermore, the composition of the BBX gene in different subgroups was also different among species (Additional file 14: Figure S8). In grapevine, the numbers of BBX members with two tandem B-BOXes plus the CCT domain, two tandem B-BOXes, Box 1 plus CCT, and B-Box 1 only were 8, 9, 2 and 5, respectively (Additional file 14: Figure S8). The corresponding numbers were 13, 8, 4, and 7 in Arabidopsis , while 8, 10, 5 and 6 were in tomato . These results demonstrated that BBX genes of different species might share a common ancestor, and underwent an independent expansion after the divergence of the monocots and the dicots .
Sequence alignment and phylogenetic analysis showed that all VvBBXs were classified into five subgroups, which are consistent with those of Arabidopsis, tomato and pear BBXs [3, 39, 40]. All VvBBX genes from subgroup I, II and IV had two B-BOX domains in grapevine, except for VvBBX11 (Fig. 1a). In contrast to animal BBX with two different types of B-BOXs, the amino acid sequences of the two B-BOX domains in grapevine were more conservative and retained the same topology . Most green algae had only a single B-box domain. However, two B-box domains were found in the unicellular green alga Chlamydomonas, implying that the B-box duplication event has taken place in some cases much before the colonization of land plants [4, 49, 50]. The rapid expansion of BBX proteins during evolution, and the fact that they are highly conserved across the plant kingdom indicates that BBX proteins might play important roles in the adaptation of land plants [4, 50].
Tandem duplication and segmental duplication events play a crucial role in the expansion of gene family members during plant evolution . To further elucidate the expansion mechanism of the VvBBX gene family in grapevine, both the tandem and segmental duplication events were analyzed. The results demonstrated that eight VvBBX genes (VvBBX7/VvBBX8, VvBBX10/VvBBX11, VvBBX13/VvBBX14, VvBBX14/VvBBX15, and VvBBX15/VvBBX16) were identified to participate in the segmental duplication (Fig. 2). Moreover, only one gene pair VvBBX11/VvBBX22 showed tandem duplication (Fig. 2), indicating that segmental duplication was more common in VvBBX genes. Taken together, these results suggested that both segmental and tandem duplications contributed to the expansion of the grapevine VvBBX gene family during their evolution. The same potential mechanism of gene family evolution was also identified in the SBP-box and WRKY gene family in grapevine [51, 52]. Furthermore, the value of Ka/Ks of segmental gene pairs was also calculated. Generally, Ka/Ks ratio greater than 1, equal to 1, and less than 1 represents positive selection, neutral selection, and negative selection, respectively. Remarkably, the Ka/Ks ratios of all grapevine gene pairs were less than 1, indicating that these gene pairs have been experiencing a markedly purifying selection during their evolution.
Potential roles of VvBBX genes in plant growth and development
Accumulating studies have shown that BBX genes are involved in multiple aspects of plant growth and development, such as seedling photomorphogenesis, shade avoidance, chlorophyll accumulation, and flower induction . The expression pattern of 24 VvBBX genes in 42 different grapevine tissues during different developmental stages were examined using an expression atlas of V.vinifera (cv. Corvina) . The gene expression analysis indicated that some VvBBX genes can be classified together according to their expression abundance in tissue-specific response of grapevine, probably reflecting their participation in a common metabolic and/or developmental process. In addition, several VvBBX genes might have a unique function in specific developmental stages.
Previous studies have also shown that BBX genes are involved in pollen tube growth. For example, the expression of PbBBX5 increased significantly in the pear pollen tubes, suggesting that it may play a critical role in pear pollen tube growth . Moreover, VvBBX4, an ortholog of PbBBX5 showed relatively high transcriptional activity in the pollen tube, indicating its potential involvement in the pollen tube growth and/or senescence. Likewise, VvBBX24 might play a similar role in grapevine, as it was also highly expressed in the pollen tube. BBX genes participate in the photoperiod pathway of flowering. In Arabidopsis, BBX6/COL5 promotes flowering under short-day by enhancing the expression of FT , whereas BBX7/COL9 negatively regulates flowering by repressing the expression of CO and FT under long day . BBX32/EIP6 represses flowering probably in a CO-independent manner under long-day . VvBBX10 and VvBBX11 were expressed at a higher level in floral tissues, signifying their putative role in the regulation of flower development. Taken together, our results proposed that VvBBX genes were not only involved in the development of pollen tubes, but also in the regulation of flower development.
AtBBX21 (also known as SALT TOLERANCE HOMOLOG 2, STH2) interacts with ELONGATED HYPOCOTYL 5 (HY5), which is considered as a key signaling regulator in photomorphogenesis, and positively regulates seedling photomorphogenesis . AtBBX21, which shows expression in dry and germinative seeds, is also involved in ABA signaling and control seed germination . In our findings, VvBBX15, the closest homolog of AtBBX21 in grapevine, showed significantly higher relative expression level in seedling and veraison seed, indicating that the VvBBX15 might be involved in seedling photomorphogenesis and seed development. On the contrary, AtBBX24 (STO) and AtBBX25 (STH) show opposite functions to AtBBX21 and suppress seedling photomorphogenesis [29, 58]. VvBBX19, which is closely related to AtBBX24 and AtBBX25, also expressed at high levels in the seedling. The result speculated that the VvBBX19 gene might be involved in grapevine photomorphogenesis and displayed opposite functions with VvBBX15. However, VvBBX19 showed more widespread and less tissue-specific transcript accumulation patterns, such as in buds, leaves, floral parts, and fruits. These results indicated that VvBBX19 genes might participate in regulatory functions during multiple growth and developmental stages. Additionally, CmBBX22 plays an important role in delaying leaf senescence and enduring drought severity in chrysanthemum . The expression of VvBBX22 was significantly higher in senescent leaf than other tissues, suggested that VvBBX22 is likely to play similar roles as CmBBX22 in regulating leaf senescence in grapevine. Taken together, all the VvBBXs were expressed at most tissues and developmental stages tested, which revealed that VvBBX genes can play multiple important roles in various developmental and biological processes.
VvBBX is likely to play important roles during grapevine berry development and ripening
Grape berry development and ripening is a complex physical and biochemical process that is controlled by various transcriptional networks and regulatory proteins, such as MYB, bHLH and MADS-box . There are growing evidences that BBX genes are involved in anthocyanin accumulation and berry ripening [56, 61, 62]. For example, AtBBX22 directly promotes the expression of the production of anthocyanin pigment 1 (PAP1) and anthocyanin biosynthesis genes by interacting with AtHY5 [56, 61]. Mutation of AtBBX22 displays lower anthocyanin accumulation by down-regulation of PAP1 expression, suggesting that AtBBX22 modulates anthocyanin synthesis by regulating the expression of MYB genes [10, 61]. Studies in apple proposed that MdBBX22 (MdCOL11) involves in MdHY5-mediated signal transduction network and positively regulates anthocyanin biosynthesis in apple peel . On the contrary, MdBBX54 (MdCOL4) suppresses anthocyanin accumulation in apple skin by inhibiting the expression of MdANS and MdUFGT, which encode genes in the anthocyanin biosynthetic pathway. Moreover, MdBBX54 indirectly inhibits the expression of MdMYB1 by interacting with MdHY5 . In rice, OsBBX14 also interacts with the OsHY5 and regulates anthocyanin biosynthesis . Furthermore, transcript accumulation of MaCOL1 in banana pulp markedly increased during natural or ethylene-induced fruit ripening, indicating that MaCOL1 is a transcription activator and might be associated with the ripening of banana fruit . PaBBX28 is highly co-expressed with ANS, CHI, F3H, and F3′H, which are deliberated as four vital genes for anthocyanin biosynthesis, suggesting that the BBX transcription factors are involved in anthocyanin biosynthesis in sweet cherry .
In grapevine, VvBBX18 and VvBBX19 share high homology with the AtBBX22 and MdBBX54, showed high expression levels in green berries (Fig. 4), indicating their potential role at the early stage of grape berry development. VvBBX1 exhibits a similar high expression with VvBBX18 and VvBBX19 in green berries, suggesting that VvBBX1 may likely to perform a function similar to VvBBX18 and VvBBX19. Interestingly, protein interaction network by STRING V11 revealed that both VvBBX18 and VvBBX19 showed interaction with the VvHY5 (Additional file 15: Figure S9 and Additional file 16: Table S7). These results are consistent with the previous studies in Arabidopsis and apple [36, 61, 63] and signifies that VvBBX18 and VvBBX19 can regulate anthocyanin accumulation or berry development by MdHY5-mediated signal transduction network. VvBBX14, the homologs of MdBBX22, exhibited the highest expression at veraison stage, implying its role in anthocyanin accumulation and berry ripening (Figs. 4, 5). Additionally, VvBBX9, VvBBX22, and VvBBX23 were preferentially expressed in ripe berries, indicating that these three genes may be involved in grapevine berry ripening process.
Potential roles of VvBBX in response to various phytohormone in grapevine
Although grapevine berries have been classified as non-climacteric fruits, several hormones including ABA and ethylene may involve in the control of ripening in grape berry. It is well known that ABA promote grape berry ripening, and ABA content in berry shows a strong increase at the end of the color turning period and during the initial stages of ripening . Exogenous ABA treatment causes an increase in berry weight, a decrease in titrable acidity, and an increase in total anthocyanin content . Futhermore, ethylene can affect the physiological processes during maturation of grapevine, including berry expansion and anthocyanin accumulation . Exogenous ethylene treatment stimulates berry coloration and enhances the expression of genes related to the anthocyanin biosynthesis . Ethylene application at veraison also leads to a berry expansion by increasing the expression of water exchange and cell wall structure genes . At present, the role of VvBBX genes in the ABA and ethylene signaling pathways is still poorly studied. In previous studies, four tomato BBX genes were induced by ETH, and all of them had the ethylene-responsive cis-element (ERE) in their promoters . In Arabidopsis, AtBBX21 coordinates with HY5 and ABI5 on the ABI5 promoter and that these transcriptional regulators work in concert to integrate light and ABA signaling . In our study, the expression of VvBBX1 was significantly decreased during the whole berry ripening period and three other VvBBX genes (i-e., VvBBX2, VvBBX3, and VvBBX13) were also suppressed at ripening stage by pre-veraison ABA and ethylene treatment, implying that these genes might be act as a negative regulator of fruit ripening in grapevine. By contrast, VvBBX22 was upregulated expression after exogenous ABA and ethylene treatment, suggesting that VvBBX22 can positively regulate fruit ripening in grapevine. All these findings indicated that these VvBBX genes might be involved in multiple hormone signaling as transcriptional regulators to modulate berry development and ripening.
In addition, VvBBX genes might be attributed to grape berry expansion and seedlessness. Multiple VvBBX genes participated in response to GA3 and CPPU treatment, which were commonly applied to increase berry size and inhibit russet development in vineyard. For example, the expression of VvBBX7, VvBBX8 and VvBBX12 was significantly increased, whereas VvBBX21, VvBBX22 and VvBBX24 was significantly decreased by GA3 and CPPU treatment (Additional file 12: Table S5). The current results indicated that the majority of VvBBX genes detected herein were up-regulated or down-regulated following the GA3 and CPPU treatment, demonstrating that BBX transcription factors might have opposite functions in grape berry expansion and seedlessness. Interestingly, the expression of VvBBX22 seemed to be responding to multiple hormone signaling, including ABA, ethylene, GA3, and CPPU, implying that VvBBX22 regulated multiple aspects of grapevine fruit ripening by modulating multiple hormone signaling pathways.
Potential roles of VvBBX genes in response to different abiotic stress in grapevine
Transcriptional regulation of stress-responsive genes is an important part of the plant response to a series of abiotic and biotic stresses. Transcription factors are master regulators that control the expression of gene clusters by binding the cis-acting element in the promoter regions of target genes. Various abiotic stresses, including Cu, salinity, waterlogging, and drought have a negative impact on grapevine growth and development. Previous studies have shown that BBX genes are involved in various abiotic stresses response [39, 41]. In this study, a series of stress-responsive cis-acting elements, such as ARE, DRE, MBS, and TC-rich, frequently occurred in the promoter regions of VvBBX genes (Additional file 9: Table S2), and are involved in drought, salt, and waterlogging. All 24 VvBBX genes possessed at least one of the stress-responsive cis-acting element, indicating their potential functions in response to abiotic stresses (Additional file 9: Table S2). In our findings, eight VvBBX (VvBBX1, 5, 15, 16, 20, 22, 23, 24) genes were induced against waterlogging stress, and all of them contained ARE elements in their promoters. Furthermore, there were four, five, ten, and fifteen VvBBX genes showing potential involvement in Cu, salt, waterlogging, and drought, respectively. The results suggested that the majority of VvBBX genes were induced or repressed at varying degrees depending on stress treatments. In addition, both VvBBX15 and VvBBX16 were co-expressed (upregulated) in salt, waterlogging and drought treatments (Fig. 7), indicating that these two genes may integrate different stress signals. The expression of VvBBX1 was enhanced by waterlogging and drought stress but was declined under Cu stress, suggesting that VvBBX1 might have different mechanisms to maintain protection against various abiotic signals. In particular, the expression of VvBBX22 was up-regulated by 7–13 fold under four abiotic stresses (Fig. 7), indicating VvBBX22 might play a vital role in response to multiple abiotic stress networks. Our result is consistent with the high expression of VvBBX22 in senescent leaf and implied that VvBBX22 might under various abiotic stresses. Characterization of VvBBX genes in response to different abiotic stresses will greatly improve our understanding of the functions and the crosstalk that occurs among different abiotic stresses signaling pathways.
In this study, 24 VvBBX genes were identified in grapevine, and the systematic and comprehensive analysis of the VvBBX gene family was performed, including conserved domain, phylogenetic relationship, gene structure, chromosome location, gene duplication, cis-acting elements, and expression pattern analysis. Numerous cis-acting elements were found in the VvBBX promoter sequences, indicating that VvBBX genes were involved in complex regulatory networks controlling development and responses to abiotic stresses. The transcription of VvBBX in different tissues and developmental stages and under various stress conditions indicated that VvBBX may have various functions in grapevine growth and development. Transcriptome and qRT-PCR analysis revealed VvBBX genes might play important roles in fruit developmental and ripening stages by modulating multiple hormone signaling pathways. Taken together, genome-wide analysis of the VvBBX will provide a solid foundation for functional analyses of BBX genes in grapevine, and further study on several BBX genes is undergoing to understand their biological functions.
Identification of BBX genes in the grapevine genome
Two different procedures were used to identify and interpret the BBX gene in grapevine as previous reports . First, all the Arabidopsis BBX protein sequences downloaded from the Arabidopsis information source (TAIR) database (http://www.arabidopsis.org) were used as queries, and grape genome database was screened by the blast p Program (E-value <1e-5). Second, the hidden Markov model (HMM) of b-box domain (pfam00643) was established by Pfam database (http://pfam.xfam.org/) and (http://genomes.cribi.unipd.it/grape/) was used to retrieve the grape genome database. Subsequently, the potential VvBBX genes in grapevine genome were further verified for the presence of the B-BOX domain by screening against the Pfam (http://pfam.sanger.ac.uk/), InterProScan (http://www.ebi.ac.uk/Tools/pfa/iprscan/) and SMART (http://smart.embl-heidelberg.de/) database. The molecular weights, isoelectric points (pI) and grand average of hydropathicity (GRAVY) of VvBBX proteins were calculated by the ExPasy website (https://web.expasy.org/protparam/). The subcellular locations of grapevine BBX proteins were predicted by WoLF PSORT (http://www.genscript.com/psort/wolf_ psort.html).
Construction of VvBBXs phylogenetic tree and sequence alignment analysis
ClustalW program (version 2.1; http://www.clustal.org/) was used for multiple sequence alignment of 24 VvBBX proteins. A phylogenetic tree based on VvBBXs protein sequence was constructed by using MEGA5.0 program neighbor-joining method and bootstrap analysis (1000 replicates) . Domains were identified with Pfam (http://pfam.xfam.org), InterProscan (http://www.ebi.ac.uk/interpro/search/sequence-search) and SMART (http:// smart.embl-heidelberg.de) programs. The motif logos of the VvBBXs were generated by submitting the sequences to the MEME website (http://meme.nbcr.net/meme/cgi-bin/meme.cgi).
Chromosomal location, gene structure, and duplication analysis
All VvBBX genes were mapped to grapevine chromosomes based on physical positions at the Grape Genome CRIBI website (http://genomes.cribi.unipd.it/) and MapInspect software was used for drawing of the map. Accordingly, the cDNA sequences and its corresponding genomic DNA sequences of VvBBX members were obtained from grape genome, and then the coding sequence and its corresponding genomic sequence were compared by using GSDs software (http://gsds.cbi.pku.edu.cn) to determine the exon-intron organization . Tandem duplicated genes were determined by detecting its physical locations on specific chromosomes and were identified as adjacent paralogous on a grape chromosome, with no more than one intervening gene. For joint analysis, download and analyze the joint blocks in the grape genome from the plant genome replication database . For duplicate pairs, Ka (nonsynonymous substitution rate) and Ks (synonymous substitution rate) and evolutionary constraint (Ka/Ks) between paralogous pairs of VvBBX genes were calculated by ParaAT and KaKs_Calculator as previous report.
Cis-element analysis for VvBBX gene promoter
The promoter sequence of 2000 bps upstream of each VvBBX genes coding regions were retrieved from CRIBI (http://genomes.cribi.unipd.it/). PlantCARE online program (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/) was used to search for assumed cis-acting element .
Expression profiles of VvBBXs in various organs and different berry developmental stages
The expression profile of VvBBX gene in a V. vinifera cv ‘Corvina’ (clone48) in different organs at different developmental stages was determined. Microarray data were attained from the NCBI gene expression omnibus (GEO) datasets under the series entry GSE36128 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE36128) . Mean of expression value of each gene in all tissues of different organs were evaluated and graphically characterized using Multi Experiment Viewer (MeV) software . Expression patterns of VvBBX genes in developmental stages of ‘fujimino’ grape fruits were acquired from the gene expression omnibus (GEO) database of NCBI (GSE77218) (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi), which is measured by RNA-sequencing (RNA-Seq) data . During the whole growing season, the fruits of three-year-old ‘fujimino’ grape trees were sampled three times in the green fruit expanding stage (40DAF or DAF40), verison stage (65DAF or DAF65), and ripe stage (90DAF or DAF90). In addition, according to the RNA sequence data downloaded from NCBI geographic data set (accession numbers GSE62744 and GSE62745), the expression of VvBBX gene in four berry development stages of 10 different grape varieties were analyzed . The 10 varieties include five red-skinned berries (Sangiovese, Barbera, Negro amaro, Refosco, and Primitivo) and five white-skinned berries (Vermentino, Garganega, Glera, Moscato Bianco, and Passerina). Berries were sampled in three replicates at four developmental stages, the pea-sized berry stage at 20d after flowering, the berries beginning to touch stage just prior to veraison (Pre_veraison), the berry-softening stage at the end of veraison (End_veraison), and the fully ripe berry stage at harvest.
Plant growth condition and different hormone treatment
The 4-years-old ‘Fujiminori’ grapevine trees grown in the standard field conditions at the Qingdao Agricultural University fruit farm, Qingdao, China, were selected as experimental material. In order to study gene expression characteristics of the VvBBX genes during berry development and ripening, grapevine berry samples were also collected at three-time points during the growing season: the green fruit expanding stage (50 DAF), veraison (70 DAF) and ripe/harvest stages (90 DAF).
Ethylene and ABA treatments (500 mg/L ethephon (ETH, an ethylene-releasing reagent) and 100 mg/L ABA in 0.02% (v/v) Tween 20 and 1% (v/v) ethanol (both used as surfactants) were sprayed on grapevine berries at green fruit expanding stage (40 DAF) and at the same time the control berries were sprayed in 0.02% (v/v) Tween 20 and 1% (v/v) ethanol. The ethephon, ABA and control solutions were applied with a handheld sprayer until run-off. During the growing season the berries samples were collected at three-time points: the green fruit expanding stage (50 DAF), veraison (70 DAF) and ripe/harvest stages (90 DAF) from treated and control plants. All samples were collected three times form each sample and were immediately frozen in liquid nitrogen and stored at − 80 °C until used. In addition, in order to study the expression profiles of VvBBXs genes, in response to GA3 and CPPU, the grapevine RNA-seq data sets were retrieved from the published supplementary data sets .
The expression of VvBBX genes under stress condition
In order to study the expression profiles of VvBBXs in response to different stress treatment (Cu, salt, waterlogging and drought stress), the grapevine RNA-seq data in response to waterlogging and drought stress were retrieved from NCBI database (SRA accession no. SRP070475 and SRP074162, respectively) (https://www.ncbi.nlm.nih.gov/sra/SRP070475 and https://www.ncbi.nlm.nih.gov/sra/?term=SRP074162) [77, 78]. RNA-seq data for expression profiles in response to Cu and salt were retrieved from published supplemental data sets [37, 79]. Summer black (hybrids of V. vinifera and V. labrusca) 2-years-old grapevine were used to detect mRNA expression of VvBBX under abiotic stresses. Cu stress of potted grapevine plants was simulated with 100 μM CuSO4 and salt stress was treated with 0.8% NaCl [37, 79]. Control plants are also treated with distilled water and waterlogging treatment was carried out by immersing the plants in the water for 48 h  and drought treatment was carried out by withholding water 20 days . Grapevine plantlets were used as a control and have grown in an appropriate conditions. All types of samples were copied three times, and the expanded leaves from the third and fourth pieces of the shooting vertices were collected from treatment and control group during deep sequencing. The analysis of RNA-seq data was based on the previous method , and the RPKM (Reads Per Kilobase per Million mapped reads) values were used to estimate the gene expression level. The heatmap of VvBBX genes was revealed using R software (http://www.bioconductor.org/).
Total RNA samples were extracted by TIANGEN RNAprep pure (Tiangen, Beijing, China) agreeing to the manufacturer’s instructions, and to remove the residual contaminated genomic DNA the RNA was processed with DNase I (RNase free; TaKaRa Biotechnology, Dalian, China). Subsequently, the first-strand cDNA was reversed from 1.5 μg total RNA using PrimeScript RTase (TaKaRa Biotechnology, Dalian, China). The primers that were used for qRT-PCR were designed by Primer 3.0 online software and the primer sequences details were presented in Additional file 17: Table S8. The qRT-PCR was performed using SYBR® Premixm Ex Taq™ (TaKaRa, Japan) with the Applied Biosystems 7500 Real-Time PCR System. All the experiments were carried out with three biological replicates. The relative levels of gene expression were determined by the 2–ΔΔCT methods with Actin (AB073011) serving as a housekeeping gene .
Subcellular localization analysis
In order to verify the subcellular localization of the VvBBX, without the stop codon of four full-length open reading frames (ORFs) of the VvBBX genes (VvBBX1, 9, 13 and 22) were amplified from the cDNAs of four-years-old ‘Fujiminori’ grapevine. Under the control of the CaMV35S promoter, the amplification product is cloned into a pCAMBIA1300 vector with green fluorescence protein (GFP). The gene-specific primers were presented in Additional file 17: Table S8. Through electroporation, four VvBBXs recombinant plasmids were transfected into the Agrobacterium tumefaciens strain EH105. Transient detection of Agrobacterium-mediated infected into the epidermal cells of Nicotiana benthamiana, as described by Sparkes et al. (2006) . The transient expression of VvBBX-GFP was observed using a laser confocal microscopes (Zeiss LSM700, Germany) and the nucleus was visualized using the mCherry-labeled nuclear marker (NF-YA4-m Cherry).
Availability of data and materials
The Arabidopsis, pear and grapevine BBX protein sequences were downloaded from the Arabidopsis information source (TAIR) database (http://www.arabidopsis.org) and GigaDB database (http://gigadb.org/site/index) and Grape Genome CRIBI website (http://genomes.cribi.unipd.it/), respectively. Microarray datas for expression profiles of V. vinifera cv ‘Corvina’ (clone48) in different organs at different developmental stages were downloaded from NCBI gene expression omnibus (GEO) datasets (accession numbers: GSE36128) (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE36128). RNA-seq data for expression profiles of ‘Fujimino’ grape fruits in developmental stages were acquired from NCBI gene expression omnibus (GEO) database (GSE77218) (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi). RNA-seq data for expression profiles of 10 different grape varieties in four berry development stages were downloaded from NCBI gene expression omnibus (GEO) database (GSE62744 and GSE62745) (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi). RNA-seq data in response to GA3 and CPPU were retrieved from the published supplementary data sets (https://doi.org/10.1016/j.scienta.2019.02.048). The grapevine RNA-seq data in response to waterlogging and drought stress were retrieved from NCBI database (SRA accession numbers: SRP070475 and SRP074162, respectively) (https://www.ncbi.nlm.nih.gov/sra/SRP070475 and https://www.ncbi.nlm.nih.gov/sra/?term=SRP074162).
Abscisic acid responsive element
Constitutively Photomorphogenic 1
Days after flowering
Flowering Locus T
Green fluorescence protein
Hidden Markov Model
Protein long hypocotyl 5
Multiple Em for Motif Elicitation
Open reading frames
The production of anthocyanin pigment 1
Reads Per Kilobase per Million mapped reads
SALT TOLERANCE HOMOLOG 2
Salt tolerance protein
Riechmann JL, Heard JE, Reuber TL, Jiang CZ, Keddie J, Adam L, et al. Arabidopsis transcription factors: genome-wide comparative analysis among eukaryotes. Science. 2000;290:2105–10.
Kiełbowicz-Matuk A. Involvement of plant C2H2-type zinc finger transcription factors in stress responses. Plant Sci. 2012;185-186:78–85.
Khanna R, Kronmiller B, Maszle DR, Coupland G, Holm M, Mizuno T, et al. The Arabidopsis B-box zinc finger family. Plant Cell. 2009;21:3416–20.
Gangappa SN, Botto JF. The BBX family of plant transcription factors. Trends Plant Sci. 2014;19:460–70.
Griffiths S, Dunford RP, Coupland G, Laurie DA. The evolution of CONSTANS-like gene families in barley, rice, and Arabidopsis. Plant Physiol. 2003;131:1855–67.
Yan H, Marquardt K, Indorf M, Jutt D, Kircher S, Neuhaus G, et al. Nuclear localization and interaction with COP1 are required for STO/BBX24 function during photomorphogenesis. Plant Physiol. 2011;156:1772–82.
Gendron JM, Pruneda-Paz JL, Doherty CJ, Gross AM, Kang SE, Kay SA. Arabidopsis circadian clock protein, TOC1, is a DNA- binding transcription factor. Proc Natl Acad Sci U S A. 2012;109:3167–72.
Gonzálezschain ND, Díazmendoza M, Zurczak M, Suárezlópez P. Potato CONSTANS is involved in photoperiodic tuberization in a graft-transmissible manner. Plant J. 2012;70(4):678–90.
Lee YS, Jeong DH, Lee DY, Yi J, Ryu CH, Kim SL, et al. OsCOL4 is a constitutive flowering repressor upstream of Ehd1 and downstream of OsphyB. Plant J. 2010;63:18–30.
Chang CS, Maloof JN, Wu SH. COP1-mediated degradation of BBX22/LZF1 optimizes seedling development in Arabidopsis. Plant Physiol. 2011;156(1):228–39.
Lin F, Jiang Y, Li J, Yan TT, Fan LM, Liang JS, et al. B-BOX DOMAIN PROTEIN28 negatively regulates photomorphogenesis by repressing the activity of transcription factor HY5 and undergoes COP1-mediated degradation. Plant Cell. 2018;30:2006–19.
Crocco CD, Holm M, Yanovsky MJ, Botto JF. Function of B-BOX under shade. Plant Signal Behav. 2011;6(1):101–4.
Xiong C, Luo D, Lin AH, Zhang CL, Shan LB, He P, et al. A tomato B-box protein SlBBX20 modulates carotenoid biosynthesis by directly activating PHYTOENE SYNTHASE 1, and is targeted for 26S proteasome-mediated degradation. New Phytol. 2019;221:279–94.
Soitamo AJ, Piippo M, Allahverdiyeva Y, Battchikova N, Aro EM. Light has a specific role in modulating Arabidopsis gene expression at low temperature. BMC Plant Biol. 2008;8:13.
Ding L, Wang S, Song ZT, Jiang YP, Han JJ, Lu SJ, et al. Two B-Box domain proteins, BBX18 and BBX23, interact with ELF3 and regulate thermomorphogenesis in Arabidopsis. Cell Rep. 2018;25:1718–28.
Putterill J, Robson F, Lee K, Simon R, Coupland G. The CONSTANS gene of Arabidopsis promotes flowering and encodes a protein showing similarities to zinc finger transcription factors. Cell. 1995;80(6):847–57.
Robson F, Costa MMR, Hepworth SR, Vizir I, Reeves PH, Putterill J, et al. Functional importance of conserved domains in the flowering-time gene CONSTANS demonstrated by analysis of mutant alleles and transgenic plants. Plant J. 2001;28:619–31.
Yang Y, Ma C, Xu Y, Wei Q, Imtiaz M, Lan H, et al. A zinc finger protein regulates flowering time and abiotic stress tolerance in chrysanthemum by modulating gibberellin biosynthesis. Plant Cell. 2014;26:2038–54.
Suarez-Lopez P, Wheatley K, Robson F, Onouchi H, Valverde F, Coupland G. CONSTANS mediates between the circadian clock and the control of flowering in Arabidopsis. Nature. 2001;410:1116–20.
Samach A, Onouchi H, Gold SE, Ditta GS, Schwarz-sommer Z, Yanofsky MF, et al. Distinct roles of CONSTANS target genes in reproductive development of Arabidopsis. Science. 2000;288:1613–6.
Valverde F, Moueadov A, Soppe W, Ravenscroft D, Samach A, Coupland G, et al. Photoreceptor regulation of CONSTANS protein in photoperiodic flowering. Science. 2004;303:1003–6.
Datta S, Hettiarachchi GH, Deng XW, Holm M. Arabidopsis CONSTANS-LIKE3 is a positive regulator of red light signaling and root growth. Plant Cell. 2006;18:70–84.
Cheng XF, Wang ZY. Overexpression of COL9, a CONSTANS-LIKE gene, delays flowering by reducing expression of CO and FT in Arabidopsis thaliana. Plant J. 2005;43:758–68.
Tripathi P, Carvallo M, Hamilton EE, Preuss S, Kay SA. Arabidopsis B-BOX32 interacts with CONSTANS-LIKE3 to regulate flowering. Proc Natl Acad Sci U S A. 2017;114:172–7.
Xu D, Jiang Y, Li J, Lin F, Holm M, Deng XW. BBX21, an Arabidopsis B-box protein, directly activates HY5 and is targeted by COP1 for 26S proteasome-mediated degradation. Proc Natl Acad Sci U S A. 2016;113:7655–60.
Xu D, Jiang Y, Li J, Holm M, Deng XW. B-box domain protein BBX21 promotes photomorphogenesis. Plant Physiol. 2018;176:2365–75.
Zhang X, Huai J, Shang F, Xu G, Tang W, Jing Y, et al. A PIF1/PIF3-HY5-BBX23 transcription factor cascade affects photomorphogenesis. Plant Physiol. 2017;174:2487–500.
Fan XY, Sun Y, Cao DM, Bai MY, Luo XM, Yang HJ, et al. BZS1, a B-box protein, promotes photomorphogenesis downstream of both brassinosteroid and light signaling pathways. Mol Plant. 2012;5:591–600.
Gangappa SN, Holm M, Botto JF. Molecular interactions of BBX24 and BBX25 with HYH, HY5 HOMOLOG, to modulate Arabidopsis seedling development. Plant Signal Behav. 2013;8:e25208.
Wang CQ, Sarmast MK, Jiang J, Dehesh K. The transcriptional regulator BBX19 promotes hypocotyl growth by facilitating COP1-mediated EARLY FLOWERING3 degradation in Arabidbopsis. Plant Cell. 2015;27:1128–39.
Crocco CD, Holm M, Yanovsky MJ, Botto JF. AtBBX21 and COP1 genetically interact in the regulation of shade avoidance. Plant J. 2010;64:551–62.
Wang HG, Zhang ZL, Li HY, Zhao X, Liu XM, Ortiz M, et al. CONSTANS-LIKE 7 regulates branching and shade avoidance response in Arabidopsis. J Exp Bot. 2013;64:1017–24.
Lippuner V, Cyert MS, Gasser CS. Two classes of plant cDNA clones differentially complement yeast calcineurin mutants and increase salt tolerance of wild-type yeast. J Biol Chem. 1996;271:12859–66.
Nagaoka S, Takano T. Salt tolerance-related protein STO binds to a Myb transcription factor homologue and confers salt tolerance in Arabidopsis. J Exp Bot. 2003;54:2231–7.
Wang QM, Zeng JX, Deng KQ, Tu XJ, Zhao XY, Tang DY, et al. DBB1a, involved in gibberellin homeostasis, functions as a negative regulator of blue light-mediated hypocotyl elongation in Arabidopsis. Planta. 2011;233:13–23.
Bai SL, Saito T, Honda C, Hatsuyama Y, Ito A, Moriguchi T. An apple B-box protein, MdCOL11, is involved in UV-B-and temperature-induced anthocyanin biosynthesis. Planta. 2014;240(5):1051–62.
Leng XP, Jia HF, Sun X, Shangguan LF, Mu Q, Wang BJ, et al. Comparative transcriptome analysis of grapevine in response to copper stress. Sci Rep. 2015;5:17749.
Leng XP, Wang PP, Wang C, Zhu XD, Li XP, Li HY, et al. Genome-wide identification and characterization of genes involved in carotenoid metabolic in three stages of grapevine fruit development. Sci Rep. 2017;7:4216.
Chu ZN, Wang X, Li Y, Yu HY, Li JH, Lu YG, et al. Genomic organization, phylogenetic and expression analysis of the B-BOX gene family in tomato. Front Plant Sci. 2016;7:1552.
Cao YP, Han YH, Meng DD, Li DH, Jiao CY, Jin Q, et al. B-BOX genes: genome-wide identification, evolution and their contribution to pollen growth in pear (Pyrus bretschneideri Rehd.). BMC Plant Biol. 2017;17:156.
Liu X, Li R, Dai YQ, Chen XS, Wang XY. Genome-wide identification and expression analysis of the B-box gene family in the apple (Malus domestica Borkh.) genome. Mol Gen Genomics. 2018;293:303–15.
Jaillon O, Aury JM, Noel B, Policriti A, Clepet C, Casagrande A, et al. The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla. Nature. 2007;449:463–7.
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.
Massonnet M, Fasoli M, Tornielli GB, Altieri M, Sandri M, Zuccolotto P, et al. Ripening transcriptomic program in red and white grapevine varieties correlates with berry skin anthocyanin accumulation. Plant Physiol. 2017;174:2376–96.
Xu YS, Hou XD, Feng J, Khalil-Ur-Rehman M, Tao JM. Transcriptome sequencing analyses reveals mechanisms of eliminated russet by applying GA3 and CPPU on ‘Shine Muscat’ grape. Sci Hortic. 2019;250:94–103.
Huang JY, Zhao XB, Weng XY, Wang L, Xie WB. The rice B-Box zinc finger gene family: Genomic identification, characterization, expression profiling and diurnal analysis. PLoS ONE. 2012;7(10):e48242.
Peers G, Niyogi KK. Pond scum genomics: the genomes of Chlamydomonas and Ostreococcus. Plant Cell. 2008;20:502–7.
Cao YP, Meng DD, Han YH, Chen TZ, Jiao CY, Chen Y, Jin Q, Cai YP. Comparative analysis of B-BOX genes and their expression pattern analysis under various treatments in Dendrobium officinale. BMC Plant Biol. 2019;19:245.
Crocco CD, Botto JF. BBX proteins in green plants: insights into their evolution, structure, feature and functional diversification. Gene. 2013;531(1):44–52.
Kenrick P, Crane PR. The origin and early evolution of plants on land. Nature. 1997;389:33–9.
Hou HM, Li J, Gao M, Singer SD, Wang H, Mao LY, et al. Genomic organization, phylogenetic comparison and differential expression of the SBP-box family genes in grape. PLoS ONE. 2013;8(3):e59358.
Wang M, Vannozzi A, Wang G, Liang YH, Tornielli GB, Zenoni S, et al. Genome and transcriptome analysis of the grapevine (Vitis vinifera L.) WRKY gene family. Hortic Res. 2014;1:16.
Fasoli M, Dal Santo S, Zenoni S, et al. The grapevine expression atlas reveals a deep transcriptome shift driving the entire plant into a maturation program. Plant Cell. 2012;24:3489–505.
Hassidim M, Harir Y, Yakir E, Kron I, Green RM. Over-expression of CONSTANS-LIKE 5 can induce flowering in short-day grown Arabidopsis. Planta. 2009;230:481–91.
Park HY, Lee SY, Seok HY, Kim SH, Sung ZR, Moon YH. EMF1 interacts with EIP1, EIP6 or EIP9 involved in the regulation of flowering time in Arabidopsis. Plant Cell Physiol. 2011;52:1376–88.
Datta S, Hettiarachchi C, Johansson H, Holm M. SALT TOLERANCE HOMOLOG2, a B-Box protein in Arabidopsis that activates transcription and positively regulates light-mediated development. Plant Cell. 2007;19:3242–55.
Xu DQ, Li JG, Gangappa SN, Hettiarachchi C, Lin F, Andersson MX, et al. Convergence of light and ABA signaling on the ABI5 promoter. PLoS Genet. 2014;10(2):e1004197.
Gangappa SN, Crocco CD, Johansson H, Datta S, Hettiarachchi C, Holm M, et al. The Arabidopsis B-BOX protein BBX25 interacts with HY5, negatively regulating BBX22 expression to suppress seedling photomorphogenesis. Plant Cell. 2013;25:1243–57.
Liu YN, Chen H, Ping Q, Zhang ZX, Guan ZY, Fang WM, et al. The heterologous expression of CmBBX22 delays leaf senescence and improves drought tolerance in Arabidopsis. Plant Cell Rep. 2019;38:15–24.
Kuhn N, Guan L, Dai ZW, Wu BH, Lauvergeat V, Gomès E, et al. Berry ripening: recently heard through the grapevine. J Exp Bot. 2014;65(16):4543–59.
Chang C, Li Y, Chen L, Chen W, Hsieh W, Shin J, et al. LZF1, a HY5-regulated transcriptional factor, functions in Arabidopsis de-etiolation. Plant J. 2008;54:205–19.
Chen J, Chen JY, Wang JN, Kuang JF, Shan W, Lu WJ. Molecular characterization and expression profiles of MaCOL1, a CONSTANS-like gene in banana fruit. Gene. 2012;496:110–7.
Fang HC, Dong YH, Yue XX, Chen XL, He NB, Hu JF, et al. MdCOL4 interaction mediates crosstalk between UV-B and high temperature to control fruit colouration in apple. Plant Cell Physiol. https://doi.org/10.1093/pcp/pcz023.
Kim DH, Park S, Lee JY, Ha SH, Lee JG, Lim SH. A rice B-Box protein, OsBBX14, finely regulates anthocyanin biosynthesis in rice. Int J Mol Sci. 2018;19(8):2190.
Guo X, Wang YT, Zhai ZF, Huang TJ, Zhao D, Peng X, et al. Transcriptomic analysis of light-dependent anthocyanin accumulation in bicolored cherry fruits. Plant Physiol Biochem. 2018;130:663–77.
Wheeler S, Loveys B, Ford C, Davies C. The relationship between the expression of abscisic acid biosynthesis genes, accumulation of abscisic acid and the promotion of Vitis vinifera L. berry ripening by abscisic acid. Aust J Grape Wine R. 2009;15:195–204.
Chervin C, El-Kereamy A, Roustan JP, LatcheA LJ, Bouzayen M. Ethylene seems required for the berry development and ripening in grape, a non-climacteric fruit. Plant Sci. 2004;167:1301–5.
El-Kereamy A, Chervin C, Roustan JP, Cheynier V, Souquet JM, Moutounet M, et al. Exogenous ethylene stimulates the long-term expression of genes related to anthocyanin biosynthesis in grape berries. Physiol Plant. 2003;119:175–82.
Chervin C, Tira-umphon A, Terrier N, Zouine M, Severac D, Roustan JP. Stimulation of the grape berry expansion by ethylene and effects on related gene transcripts, over the ripening phase. Physiol Plant. 2008;134:534–46.
Leng XP, Wei HR, Xu XZ, Ghuge SA, Jia DJ, Liu GS, Wang YZ, Yuan YB. Genome-wide identification and transcript analysis of TCP transcription factors in grapevine. BMC Genomics. 2019;20:786.
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28:2731–9.
Hu B, Jin J, Guo YA, Zhang H, Luo J, Gao G. GSDS 2.0: an upgraded gene feature visualization server. Bioinformatics. 2014;31(8):1296.
Tang H, Wang X, Bowers JE, Ming R, Alam M, Paterson AH. Unraveling ancient hexaploidy through multiply-aligned angiosperm gene maps. Genome Res. 2008;18:1944–54.
Postel D, Vanlemmens P, Gode P, Ronco G, Villa P. Plant CARE, 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:325–7.
Saeed AI, Bhagabati NK, Braisted JC, Sharov V, Howe EA, Li JW, et al. TM4 microarray software suite. Methods Enzymol. 2006;411:134–93.
Shangguan LF, Mu Q, Fang X, Zhang KK, Jia HF, Li XY, et al. RNA-sequencing reveals biological networks during table grapevine (‘Fujiminori’) fruit development. PLoS ONE. 2017;12(1):e0170571.
Haider MS, Zhang C, Kurjogi MM, Pervaiz T, Zheng T, Zhang CB, et al. Insights into grapevine defense response against drought as revealed by biochemical, physiological and RNA-Seq analysis. Sci Rep. 2017;7:13134.
Zhu XD, Li XP, Jiu ST, Zhang KS, Wang C, Fang JG. Analysis of the regulation networks in grapevine reveals response to waterlogging stress and candidate gene-marker selection for damage severity. R Soc Open Sci. 2018;5:172253.
Guan L, Haider MS, Khan N, Nasim M, Jiu ST, Fiaz M, et al. Transcriptome sequence analysis elaborates a complex defensive mechanism of grapevine (Vitis vinifera L.) in response to salt stress. Int J Mol Sci. 2018;19(12):4019.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using realtime quantitative PCR and the 2−ΔΔCT method. Methods. 2001;25(4):402–8.
Sparkes IA, Runions J, Kearns A, Hawes C. Rapid, transient expression of fluorescent fusion proteins in tobacco plants and generation of stably transformed plants. Nat Protoc. 2006;1(4):2019–25.
We extend our thanks to the reviewers and editors for their careful reading and helpful comments on this manuscript.
This work was supported by the National Key Research and Development Program of China to XP Leng (2019YFD1001405) and YB Yuan (2016YFD0400100), the High-level Scientific Research Foundation of Qingdao Agricultural University to XP Leng (665/1118011) and JG Fang (665/1119002), Shandong Provincial Natural Science Foundation, China (ZR2017YL022), and the Qingdao People’s Livelihood Science and Technology Project (#18–8–1-428-nsh). The funders had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Ethics approval and consent to participate
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.
Additional file 1 Figure S1
. Structure of the VvBBX proteins. Numbers indicate amino acid length and position of the corresponding conserved domains. The blue, red and green boxes indicate the B-Box 1, B-Box 2 and CCT domain, respectively. The scale bar represents 100 amino acids.
Additional file 2 Figure S2
. The conserved domains in the VvBBX proteins. A, B and C represent the protein alignment of the B-Box 1, B-Box 2 and CCT domain, respectively. The x-axis indicates the conserved sequences of the domain. The height of each letter indicates the conservation of each residue across all proteins. The y-axis is a scale of the relative entropy, which reflects the conservation rate of each amino acid.
Additional file 3 Figure S3
. Multiple sequence alignments of the conserved domains of the VvBBXs. Multiple sequence alignments of the B-box 1 (A), B-box 2 (B) and CCT (C) domains are shown. The sequences were aligned using DNAMAN7.0.
Additional file 4 Figure S4
. Phylogenetic analysis of BBX genes in grapevine and pear. The full-length amino acid sequences of BBX proteins from grapevine (VvBBX) and pear (PbBBX) were aligned by ClustalX, and the phylogenetic tree was constructed using the maximum likelihood method by MEGA5.0.
Additional file 5 Figure S5.
Phylogenetic analysis of BBX genes in grapevine, pear, Arabidopsis and tomato. The full-length amino acid sequences of BBX proteins from grapevine (VvBBX), pear (PbBBX), Arabidopsis (AtBBX) and tomato (SlBBX) were aligned by ClustalX, and the phylogenetic tree was constructed using the maximum likelihood method by MEGA5.0.
Additional file 6 Table S1
. The Ka/Ks ratios and divergence between paralogous VvBBX gene pairs.
Additional file 7 Figure S6
. Gene structure of the VvBBX family generated from GSDS. The yellow block means the coding sequence (CDS), the blue block means the upstream or downstream of the genes, and the black line indicates the intron. The scale bar indicates the length of the DNA sequences.
Additional file 8 Figure S7
. Promoter Cis-regulatory elements analysis of grapevine VvBBX genes. Number of each cis-acting element in the promoter region (2.0 kb upstream of the translation start site) of VvBBX genes. Based on the functional annotation, the cis-acting elements were classified into three major classes: plant growth and development, phytohormone responsive, or abiotic and biotic stresses-related cis-acting elements (detailed results shown in Additional file 6: Table S1).
Additional file 9 Table S2.
Promoter analysis of the grapevine BBX gene family.
Additional file 10 Table S3.
Expression profiles of the grapevine VvBBX genes in different organs, tissues and developmental stages.
Additional file 11 Table S4
. Expression profiles of the grapevine VvBBX genes during three fruit developmental stages.
Additional file 12 Table S5
. Expression profiles of the VvBBX genes in response to GA3 and CPPU treatment.
Additional file 13 Table S6.
Expression profiles of the VvBBX genes in response to abiotic stress.
Additional file 14 Figure S8.
BBX family members of grapevine, Arabidopsis, rice, tomato and pear.
Additional file 15 Figure S9.
Protein interaction network of grapevine VvBBX18 and VvBBX19 by STRING.
Additional file 16 Table S7.
Protein interaction network of grapevine VvBBX18 and VvBBX19 by STRING.
Additional file 17 Table S8.
The primers sequences of VvBBX genes for qRT-PCR and gene amplification.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.
About this article
Cite this article
Wei, H., Wang, P., Chen, J. et al. Genome-wide identification and analysis of B-BOX gene family in grapevine reveal its potential functions in berry development. BMC Plant Biol 20, 72 (2020). https://doi.org/10.1186/s12870-020-2239-3
- Gene expression
- Berry development
- Stress response