Skip to main content

Identification and expression analysis of BURP domain-containing genes in jujube and their involvement in low temperature and drought response



Plant-specific BURP domain-containing genes are involved in plant development and stress responses. However, the role of BURP family in jujube (Ziziphus jujuba Mill.) has not been investigated.


In this study, 17 BURP genes belonging to four subfamilies were identified in jujube based on homology analysis, gene structures, and conserved motif confirmation. Gene duplication analysis indicated both tandem duplication and segmental duplication had contributed to ZjBURP expansion. The ZjBURPs were extensively expressed in flowers, young fruits, and jujube leaves. Transcriptomic data and qRT-PCR analysis further revealed that ZjBURPs also significantly influence fruit development, and most genes could be induced by low temperature, salinity, and drought stresses. Notably, several BURP genes significantly altered expression in response to low temperature (ZjPG1) and drought stresses (ZjBNM7, ZjBNM8, and ZjBNM9).


These results provided insights into the possible roles of ZjBURPs in jujube development and stress response. These findings would help selecting candidate ZjBURP genes for cold- and drought-tolerant jujube breeding.

Peer Review reports


Plants suffer diverse abiotic and biotic stresses, such as low temperature, drought, soil salinity, and damage from several diseases and pests, during their lifecycle. Stress influenced plant growth, development, and qualities in varying degree. In contrast, plants have also evolved certain mechanisms to settle in disadvantageous condition. Unraveling the mechanism of plant response to stress would be helpful for crops and fruit trees to gain knowledge needed for maintaining the agriculture and food security with increasing world population [1]. Recently, progresses have been made to elucidate plant response by identifying stress-related genes, metabolites, and pathways [2]. Some studies proved that a stress-responsive gene family encoding BURP domain-containing protein is plant-specific and performs important functions in plant development and stress response [3, 4]. The BURP proteins contain a BURP domain with conserved structures of CHX10CHX23–37CHX23–26CHX8W [5] and has been identified in some plant species, including rice (Oryza sativa) [6], soybean (Glycine max) [7], maize (Zea mays) [8], poplar (Populus trichocarpa) [9], cotton (Gossypium raimondii, Gossypium arboretum, and Gossypium hirsutum) [10], common bean (Phaseolus vulgaris) [11], Chinese rose (Rosa chinensis) [12], and legumes (Phaseolus vulgaris, Cicer arietinum, Cajanus cajan, and Vigna radiata) [13].

The BURP family is classified into BNM2-like, USP-like, RD22-like, and PG1β-like subfamilies; moreover, other subfamilies specifically exist in few species. Members belonging to different subfamilies has varied expression and function to maintain plant development. For instance, BnBNM2 is expressed in Brassica napus during microspore embryogenesis [14,15,16]. VfUSP participates in the early development of zygotic embryogenesis of field bean [17, 18]. Additionally, anther-specific BURP proteins, OsRAFTIN1 and RA8, were both important for anther dehiscence in rice [4, 19]. Moreover, SCB1 extensively affects seed coat formation [20]. Moreover, BURP proteins have also displayed their critical functions in fruit development and ripening. For example, VvBURP1 is involved in early fruit morphogenesis of Vitis vinifera [21], and PG1β monitors pectin solubilization and degradation during tomato maturation [22, 23].

The BURPs may be also induced by phytohormones and abiotic stress treatments, including abscisic acid (ABA), salicylic acid (SA), NaCl, cold, and drought. For instance, AtUSPL1, belonging to RD22 subfamily, suppresses ABA-mediated drought stress response [24]. BgBDC1, 2, 3, and 4 also have similar responses to salt and drought in an ABA-mediated pathway in mangroves [25]. Besides, both ABA and SA may induce all RD22-like GhBURPs in cotton [10]. Additionally, OsBURPs and most PtBURPs are significantly expressed in Oryza sativa and poplar, respectively, in response to at least one stress treatment, including cold, salt, drought, and ABA. These results indicated the potential role of BURPs in plant stress. However, the function of these BURPs remains to be further investigated. We have previously identified two BURP genes possibly associated with the variation of several cyclopeptide alkaloids in 214 jujube (Ziziphus jujuba Mill.) accessions using mGWAS analysis [26]. Both BURP genes showed consistent expression along with the metabolite accumulation. Therefore, the BURP domain-containing genes in jujube genome should be identified and investigate their potential roles in jujube development and stress response.

Jujube is an economically and ecologically important fruit tree with immense edible and medicinal values [27, 28]. Moreover, jujube is highly adaptative and tolerant to adverse environmental conditions, and wild jujube colonizes barren mountains. Thus, as jujube is a pioneer tree, revealing and studying its stress-responsive candidate genes is of vital importance. Till date, the involvement of several transcription factors, such as WRKY [29], NAC [30], and AP2/ERFs [31], have been identified in stress response and fruit development in jujube, while few downstream functional protein gene families have been investigated [32, 33]. This study reports a functional protein gene family in jujube and attempts to understand their variation and candidate functions using a series of bioinformation analysis. ZjBURP transcript expression in different tissues during development processes and in jujube leaf under various stress conditions were further studied. Our findings provide an important foundation for future functional studies of BURP genes in jujube, and will be potentially useful for jujube molecular breeding to improve stress resistance and fruit quality, particularly in cold- and drought-tolerant jujube breeding.


Identification, chromosomal location, and gene duplication of BURPs in jujube

A total of 17 putatively encoded BURP genes were identified from jujube genome. These genes were designated according to their phylogenetic relationship and previously reported homologous groups in other species. Notably, our result corrected a false induced by structural annotation in jujube reference genome, and therefore we manually identified two gene sequences as evm.model.Contig91.143 (ZjPG2) and evm.model.Contig91.146 (ZjPG3) from the genomic region. The information of 17 candidate genes, including their amino acid length, molecular weight (Mw), isoelectric point (pI), subcellular localization, and physical locations on chromosome, is summarized in Table 1. The length of ZjBURP protein sequences ranged from 215 (ZjBURP5) to 634 (ZjPG2 and ZjPG3). The Mw of ZjBURP proteins were between 24.08658 (ZjBURP5) and 70.05929 (ZjPG3), and the pI values ranged from 5.25 (ZjBURP4) to 8.93 (ZjBURP1).

Table 1 Summary of sequence bioinformation analysis of 17 BURP genes identified in jujube

The 17 ZjBURPs were unevenly located on five chromosomes (Table 1 and Fig. 1). A total of 15 BURPs located on chr 2, 5, and 6 and 1 ZjBURP was anchored on chr 7 and 10. Interestingly, six genes (ZjBNM1-6) showed close distribution on chr 6, and four genes (ZjBURP2-5) were also located near each other on chr 5. We further performed gene duplication analysis and demonstrated that two tandem duplication events (ZjBURP4 / ZjBURP5, ZjBNM1 / ZjBNM2) and one segment duplication event (ZjPG1 / ZjPG2) existed in nearby clusters.

Fig. 1
figure 1

Chromosomal location and gene duplication events of ZjBURPs. The scale bar is in million bases (Mb). The genes connected with dotted gray and red lines represent the segmental (gray) and tandem (red) gene pairs, respectively

Phylogenetic analysis and classification of the BURP gene family in jujube

To explore the evolutionary relationships of BURP proteins, a phylogenetic tree was generated by clustering 135 BURP proteins from different species, which were then classified into eight distinct subfamilies (BNM2- like, USP-like, RD22-like, PG1β-like, BURP V, BURP VI, BURP VII, and BURP VIII; Fig. 2). All 17 ZjBURPs belonged to four subfamilies, including BNM2-like, RD22-like, PG1β-like and BURP V.

Fig. 2
figure 2

The phylogenetic tree of BURP proteins and ZjBURP classification. The protein sequences from Brachypodium distachyon (10), Setaria italica (13), Oryza sativa (12), Sorghum bicolor (10), Zea mays (9), Arabidopsis thaliana (5), Glycine max (21), Cucumis sativus (6), Citrus sinensis (6), Brassica rapa (9), Populus trichocarpa (16), Vicia faba (1), and Ziziphus jujuba (17) were aligned using MUSCLE algorithm in MEGA-X, and the phylogenetic tree was constructed by MEGA-X with the neighbor-joining (NJ) method and 1,000 bootstrap replicates

Additionally, BNM2-like subfamily contained nine genes (ZjBNM1–9) adjacent to each other in the phylogenetic tree except ZjBNM6. BURP V subfamily consisted of four genes (ZjBURP2–5), which were also closely located. Two genes of PG1β-like subfamily, ZjPG2 and ZjPG3, were close except ZjPG1. The relative position of 17 genes in phylogenetic tree were extremely similar to the distribution of them on chromosomes, suggesting that the gene locations on chromosomes might be consistent with their evolutionary relationship.

Sequence alignment of BURP proteins in jujube

The signal peptides and BURP domains of these proteins in jujube were detected. The results showed that 17 and 12 BURP proteins contained the BURP domain and signal peptides, respectively (Fig. S1). Moreover, BURP proteins with close evolutionary relationships showed similar signal peptide and BURP domain composition and arrangement. The multiple sequence alignment analysis of these proteins revealed the existence of many highly conserved amino acid residues and four CH motifs. The results indicated that they were vital for the basic function of BURP family in jujube. Moreover, the C-terminus of BURP proteins in jujube could be concluded as CHX5YX4CHX25–34CHXDTX2WX6FX2LX4GX3 CHX8W with more conserved sequences between the first two CH motifs (Fig. S2).

Gene structure and conserved motifs of BURPs in jujube

To investigate BURP conservation and diversification in jujube, gene structures and conserved motifs were further analyzed (Fig. 3). Most ZjBURPs contained at least one intron, except ZjBNM1, ZjBNM2, and ZjPG1.

Fig. 3
figure 3

The phylogenetic relationships, gene structure, and conserved motifs of ZjBURPs. a The construction of the phylogenetic tree with the neighbor-joining method. b ZjBURP gene structure. c The conserved motifs identified in ZjBURP gene family

Most members from BNM2-like and PG1β-like subfamilies contained two exons, whereas genes from RD22-like and BURP V subfamilies contained 2–4 exons. These results showed that members from the same subfamily owned similar gene structure (Fig. 3a, b).

We further detected 15 conserved motifs utilizing MEME to explore the similarity and differentiation of BURP proteins in different subfamilies (Fig. 3c). The BURP domain in C-terminus contained motifs 4, 1, 8, 3, 5. Besides, five genes (ZjBNM1-5) from BNM2-like subfamily also contained motif 7, and all genes contained motif 13 in the C-terminus except those belonging to PG1β-like subfamily. Notably, motifs 6, 8, and 9 and motifs 10, 11, 14, and 15 existed only in BNM2-like and PG1β-like subfamilies, respectively. Overall, motif compositions and arrangements were particularly similar in ZjBURPs with close relationship.

Analysis of cis-elements in ZjBURP promoter regions

To investigate the possible regulatory functions of cis-elements in the putative promoter regions of ZjBURPs, stress- and phytohormone-related cis- elements in the 3,000 bp upstream of the start codons in all 17 ZjBURPs were identified (Fig. 4 and Table S1). The identified phytohormone-related elements were AuxRR-core (auxin), TGA-element (auxin), ABRE (ABA), P-box (gibberellin), TATC-box (gibberellin), GARE-motif.

Fig. 4
figure 4

Distribution of stress- and phytohormone-related cis-elements in the putative ZjBURP promoter region. The location of these cis-elements was predicted by plantCARE database

(gibberellin), CGTCA-motif (MeJA), TGACG-motif (MeJA), and TCA-element (SA). Stress-responsive regulatory elements included MBS (drought), TC-rich repeats (defense and stress responsiveness), and LTR (cold stress). All ZjBURPs, expect ZjBNM4, contained at least one putative stress-related elements. Moreover, the promoters of some ZjBURPs, such as ZjBURP4 (1 MBS, 2 LTRs and 1 TC-rich repeats) and ZjBURP2 (1 MBS, 1 LTR and 1 TC-rich repeats), contained multiple stress-responsive elements. Notably, low temperature-related elements were more common than drought- and defense-related elements, and ZjPG1 owned up to three low temperature elements. All members from RD22-like and BURP V subfamilies contained one drought-responsive element. In addition, ABA- and MeJA-related elements were more general than other phytohormone-related elements, indicating that ZjBURPs might take part in the ABA or MeJA regulatory pathway.

Transcript expression of ZjBURPs in different tissues during development

To gain insight into the potential function of ZjBURPs during jujube development, we estimated their. transcription level in different jujube tissues, including flowers, young fruit, and leaf by RT-PCR. The results showed that all ZjBURPs had various expression levels in three tissues (Fig. 5a, Fig.S3). Four genes (ZjBURP1, 2, 3, and ZjPG1) expression was high in all tissues (Fig. 5a). However, ZjBNM2 expression was extremely low in tested tissues, suggesting that it may be a redundant gene or pseudogene, or expressed only in specific organs or at specific developmental processes or under specific conditions. In addition, three BNM2-like subfamily members (ZjBNM7, 8, and 9) expression in reproductive organs, like flowers and young fruit, was higher than that in leaves, implying their possible roles in fruit development. Besides, ZjBURP4 showed dominant expression in leaves.

Fig. 5
figure 5

ZjBURP expression in three tissues and during fruit development processes. a RT-PCR analysis of ZjBURPs in three tissues. F: flower; YF: young fruit; L: leaf. b Transcriptional abundance of ZjBURPs during fruit development of ‘Junzao’. FPKM was selected to determine ZjBURP expression levels. JYF, JE, JWM, JBR, JHR, and JFR represent the young fruit, enlargement, white maturity, beginning-red, half-red, and full- red stages of ‘Junzao,’ respectively. c Transcriptional abundance of ZjBURPs during fruit development of ‘Qingjiansuanzao’. WYF, WE, WWM, WBR, WHR, and WFR represented the young fruit, enlargement, white maturity, beginning-red, half-red, and full-red stages of ‘Qingjiansuanzao,’ respectively

We further analyzed the ZjBURP expression levels during wild (Qingjiansuanzao) and cultivated jujube (Z. jujuba cv. Junzao) fruit development with in-house transcriptomic data to investigate their potential role in fruit development. ZjBNM8 was not expressed in either wild or cultivated jujube fruit (Fig. 5b, c). The results revealed that 16 ZjBURPs showed similar expression patterns, with preferential expression at young fruit and enlargement between the two accessions during fruit development, which indicated that ZjBURPs may play important roles in young fruit development and enlargement in both wild and cultivated jujube. Furthermore, transcript abundance at young fruit and enlargement varied among these genes, and six genes (ZjBNM4, ZjBNM7, ZjBURP1, ZjBURP2, ZjBURP3, and ZjPG1) had higher transcript abundance (FPKM ≥ 30) than other genes.

ZjBURP expression in wild jujube leaf in response to different stress conditions

Based on the analysis of cis-elements in promoter regions and previous reports on BURP genes in response to diverse stresses in other plants, ZjBURPs might be stress-responsive in jujube. To examine whether ZjBURPs are involved in various stresses, the transcript levels of all ZjBURPs were investigated in the leaves of wild jujube seedlings under three treatments including low temperature, salt, and drought stresses.

All ZjBURPs were more or less upregulated or downregulated under low temperature (4 ℃), revealing the possible roles of ZjBURPs in response to cold stress (Fig. 6). Furthermore, two main expression patterns were found. One expression pattern included the downregulation of ZjBURPs, including ZjBNM2, ZjBNM6, ZjPG2, and ZjPG3, at all timepoints after cold treatment. The other expression pattern involved the upregulation of ZjBURPs at certain timepoints and downregulation at other timepoints, and it could be divided into five categories: upregulation at 6 h, including ZjBNM1, ZjBNM3, ZjBNM4, ZjBNM5, and ZjPG1; upregulation at 12 h, containing ZjBURP2 and ZjBURP3; upregulation at 24 h, comprising ZjBURP4 and ZjBURP5; upregulation at 72 h, consisting of ZjBNM7, ZjBNM8, and ZjBNM9; upregulation at 6 h and 72 h, exemplified by ZjBURP1. The results showed that the differential expression characteristics of ZjBURPs existed also in the same subfamily, and those genes with close relationship might share similar expression patterns.

Fig. 6
figure 6

ZjBURP expression in low temperature. 0 h, 6 h, 12 h, 24 h, and 72 h represent the time period after the treatment. The error bars indicate the standard deviation of three biological replicates. Different letters represent a significant variation between the experimental groups (p < 0.05)

Under salt stress condition, only two genes (ZjBURP4 and ZjBURP5; both belonging to BURP V subfamily) were downregulated, whereas the remaining 15 genes were upregulated with different degrees and various expression characteristics (Fig. 7). Among the upregulated genes, four expression patterns were observed. One was upregulation at 12 h and then downregulation, such as that for ZjBNM1, ZjBNM3, ZjBNM4, ZjBNM5, ZjBURP1, and ZjBURP3. The second expression pattern was upregulation at 12 h and 48 h, such as that for ZjBNM7, ZjBNM8, ZjBNM9, ZjPG2, and ZjPG3. The third expression pattern was upregulation at 24 h, such as that for ZjBNM6, ZjBURP2, and ZjPG1. ZjBNM2 represented the last expression pattern, with upregulation at 48 h. The results revealed that all ZjBURPs could be induced by salt stress; notably, ZjBNM8 expression level was the highest, with over 60-fold upregulation.

Fig. 7
figure 7

ZjBURP gene expression under NaCl treatment. 0 h, 12 h, 24 h, 48 h, and 72 h represent the time period after NaCl treatment. The error bars indicate the standard deviation of three biological replicates. Different letters represent a significant variation between the experimental groups (p < 0.05)

Four genes (ZjBURP3, ZjBURP4, ZjBURP5 and ZjPG1) were downregulated, while other genes were upregulated in seedlings under drought stress (Fig. 8). The other 13 ZjBURPs exhibited four expression patterns. One was upregulation of ZjBURPs at 12 h and 48 h, including ZjBNM1, ZjBNM2 (slight downregulation at 24 h), ZjBNM3, ZjBNM4, and ZjBNM5. The second expression pattern was significant upregulation at 24 h, such as that for ZjBNM7, ZjBNM8, ZjBNM9, ZjBURP2, and ZjPG3. The third expression pattern was upregulation at 48 h, such as that for ZjBNM6 and ZjPG2. The last expression pattern was specific for ZjBURP1, with upregulation at 72 h. Among these genes, three genes (ZjBNM7, ZjBNM8, and ZjBNM9). were mostly induced approximately six-folds after 24 h drought stress, indicating they might perform significant functions in response to drought stress in jujube. We also noticed that six genes (ZjBNM1, ZjBNM3, ZjBNM4, ZjBNM5, ZjBURP2, and ZjBURP5) exhibited similar expression characteristics under salt and drought stress.

Fig. 8
figure 8

ZjBURP gene expression under PEG treatment. 0 h, 12 h, 24 h, 48 h, and 72 h represent the time period after PEG treatment. The error bars indicate the standard deviation of three biological replicates. Different letters represent a significant variation between the experimental groups (p < 0.05)


BURP genes have been extensively analyzed in several plants for their significant roles in plant development and stress response [18, 22, 23, 34]. However, this gene family has not been studied in jujube. This study reports a comprehensive analysis of the BURP gene family in jujube, including analysis of chromosomal location, gene duplication events, phylogeny, sequence alignment, gene structure, conserved motifs, cis-elements in the promoter region, and expression analysis. A total of 17 BURP genes were identified in ‘Junzao’ genome, and the variation of length of these sequences suggested high complexity within the ZjBURP gene family, as also reported in Medicago [35]. In addition, the genome sequences of ‘Dongzao’ and ‘Wild jujube’ have also been reported recently [36, 37]; therefore, we can obtain more insight into the synteny and evolution of BURP gene family among wild, fresh, and dry jujube. Moreover, the chromosomal location of ZjBURPs in ‘Junzao’ genome may also provide clues for the formation and evolution of this gene family in jujube. Four BURP gene clusters were distributed on three jujube chromosomes, and genes belonging to the same cluster might originate from a common ancestor and formed through gene duplication. Further analysis of gene duplication events showed that two and one BURP gene pairs were involved in tandem and segmental duplication events, respectively, indicating that both tandem and segmental duplication events served as the vital driving force during the long-term of BURP gene evolution in jujube. Similar results have been observed in G. raimondii and G. hirsutum [10]; however, in many plants, only one gene duplication event occurred in BURP gene family, such as one tandem duplication in G. arboretum [10], six tandem duplication events in Medicago [35], and one segmental duplication event in R. chinensis [12]. These results implied that BURP genes are subjected to strong positive selection in numerous plants, and conformed to different expansion models in diverse plants [38].

According to the phylogenetic tree, the BURP members from 13 plants could be divided into eight subfamilies: BNM2-like, USP-like, RD22-like, PG1β-like, BURP V, BURP VI, BURP VII, and BURP VIII. Interestingly, the BURP genes belonging to BNM2-like, USP-like, and BURP V subfamilies were only found in investigated dicotyledonous species (dicots), whereas those belonging to the BURP VI and BURP VII subfamilies contained only the BURP members from investigated monocotyledonous species (monocots). These results were consistent with the results of many previous studies [10, 38]. The PG1β-like subfamily was observed in both dicots and monocots, indicating that PG1β-like subfamily might have originated before the divergence of monocot and dicot plants, and BNM2-like, USP-like, as well as BURP V subfamilies, might have subsequently evolved separately and exert specific functions in dicots. Some results were distinguished from those of previous studies due to the number of species analyzed in the phylogenetic analysis for more reliable results of classification. Notably, BURP proteins from cotton, poplar, and jujube showed similar distributions in four subfamilies (BNM2-like, RD22-like, PG1β-like and BURP V).

The further structural analysis showed that all 17 BURP proteins in jujube included the conserved BURP domain, particularly the four CH motifs, which were also detected in BURP sequences in previously studied plants, like cotton [10], rose [12], Medicago [35], and rice [6], hinting that they may be crucial for the basic function of BURP proteins. Moreover, the component motifs and their arrangement in ZjBURPs were conserved within a subfamily and different in various subfamilies, which correspond to the results of the structure analysis of ZjBURPs. Similar findings were obtained in previous studies and might be relative with the conserved and diverse functions of BURP genes [7, 8].

To elucidate the possible functions of ZjBURPs in jujube development, their expression in different organs during development processes were analyzed. The results of expression patterns in three organs suggested that all ZjBURPs were expressed in flowers, young fruits and leaves, and some ZjBURPs with high sequence similarity also exhibited similar expression characteristics. For instance, ZjPG2 and ZjPG3, which shared similar sequences, were both transcribed in all three organs and highly expressed in flowers. Three closely clustered genes, ZjBNM7, ZjBNM8 and ZjBNM9, were highly expressed in flowers and young fruits, and weakly expressed in leaves. However, not all genes with close relationship had the same expression pattern, such as ZjBURP4 and ZjBURP5, as gene duplication pairs, showed different expression characteristics. ZjBURP4 expression was largely restricted to the leaves, whereas ZjBURP5 was also transcribed in flowers. The conserved and diversification of BURP expression patterns have also been observed in soybean [7]. Notably, ZjBNM2 transcript levels in all three tissues were extremely low, however, the degree of its upregulation under NaCl treatment was significant, indicating that it was induced by salt stress and its expression levels might be suppressed under normal conditions. Furthermore, the results of BURP expression patterns during fruit development showed that all 17 ZjBURPs were preferentially expressed during young fruit development and fruit enlargement in both ‘Junzao’ and ‘Qingjiansuanzao,’ indicating that ZjBURPs may execute important functions in young fruit development and enlargement, which was supported by the results of a previous study on a BURP protein associated with grapevine fruit development [21, 39]. In particular, several ZjBURPs were strikingly highly expressed (over 400- and 200-fold) during fruit development and enlargement, respectively. For instance, ZjBNM7 is homologous to AtUSPL1, which associates with the seed development [40]; ZjBURP1, ZjBURP2, and ZjBURP3 shared similar sequence with GhRDL1 in cotton, which is predominantly transcribed during fiber rapid elongation and promotes fiber elongation by interacting with GhEXPA1 [41]; ZjPG1 is homologous to AtPGL3, which promotes cell enlargement [42]; these ZjBURPs were highly expressed during young fruit development and fruit enlargement, indicating their key roles in regulating jujube development and enlargement, which should be verified in future functional analysis.

Besides, BURP genes are also pervasively involved in plant stress response in many plants. To investigate whether ZjBURPs were similarly responsive to stresses, we analyzed their stress-related cis-elements in the promoter region using qRT-PCR under different stresses, including low temperature, salt, and drought. The analysis of cis-elements showed that all promoters of ZjBURPs contained at least two investigated stress-related elements, implying that these genes might respond to different stresses. Furthermore, qRT-PCR analysis reports that the transcription of all ZjBURPs were induced by cold, salinity, and drought stresses, although several genes were slightly induced. Under low temperature condition, ZjPG2 and ZjPG3 expression were significantly downregulated, which were similar to the expression pattern of OsBURP12, which belongs to the same subfamily, in rice [6]. Moreover, ZjPG1 expression was similar to that of OsBURP16, which also belongs to the PG1β-like subfamily and decreases cold tolerance by enhancing PG activity and reducing the pectin content of cell wall in OsBURP16-transgenic plants [43, 44]. Therefore, we may conclude that ZjPG1 downregulation is important for jujube cold tolerance. Consistently, ZjPG1 contained the most low temperature-responsive cis-elements among all the BURP genes in jujube, indicating that it might be involved in jujube cold response through the low temperature-responsive cis-elements. Most ZjBURPs, except ZjBURP4 and ZjBURP5, were upregulated under NaCl treatment. In particular, ZjBNM7, ZjBNM8, and ZjBNM9 expression levels were elevated 12 h and 48 h after NaCl treatment, indicating that these genes could positively respond to salt stress. Under drought treatment, several ZjBURPs, including ZjBNM1, ZjBNM3, ZjBNM4, ZjBNM5, ZjBURP4, and ZjBURP5, possessed similar expression characteristics to those under salt stress. In addition, ZjPG1 was downregulated under drought treatment, which was similar to MtBURP33 and MtBURP28 downregulation in Medicago [35]. ZjBNM7, ZjBNM8, and ZjBNM9 expression levels were strikingly upregulated by drought stress, indicating that they might be key BURP genes involved in jujube drought response. This result was supported by their close relationship to AtUSPL1, which is upregulated as part of the ABA-mediated moisture stress response and involved in Arabidopsis thaliana drought tolerance [24, 34, 45, 46]. Interestingly, no drought-responsive cis-elements were found in their promoter regions, revealing that the results of qRT-PCR are not always consistent with cis-elements analysis [35], suggesting that many BURP genes with no related cis-elements identified in this study may have other regulation patterns.


In this study, BURP genes in the jujube genome were systematically analyzed. The 17 identified ZjBURPs could be classified into four subfamilies based on phylogeny, gene structure and conserved motif analysis. Gene duplication analysis indicated that both tandem and segmental duplication events might contribute to the ZjBURP gene family expansion. Expression analysis showed that all ZjBURPs were transcribed in flowers, young fruits and leaves, and four genes (ZjBURP1, ZjBURP2, ZjBURP3 and ZjPG1) were highly expressed in all three tissues. Transcriptome data revealed that ZjBURPs were preferentially expressed at young fruit and enlargement stages in both ‘Junzao’ and ‘Qingjiansuanzao,’ indicating their possible roles in regulating young fruit development and enlargement in jujube. The qRT-PCR analysis of ZjBURP genes under various stress treatments indicated that all ZjBURPs were induced by low temperature, salt and drought, even though some genes were slightly induced. We also concluded that several key BURP genes may perform significant functions in response to low temperature (ZjPG1) and drought stresses (ZjBNM7, ZjBNM8, and ZjBNM9), and their functions needs to be confirmed in future studies. This study provides a comprehensive view of the ZjBURP gene family and may serve as foundation for future breeding of cold/drought-tolerant jujube trees.


Plant material and treatment

Wild jujube seeds were collected from the Germplasm Resource Nursery (Xingtai, China) and sowed in sterilized soil after rinsing with tap water for 24 h. All seedlings were grown at 24 ℃ in a climate-controlled glasshouse (light/dark cycle: 16 h/8 h), and four-week-old seedlings with uniform height and biomass were selected for treatments. For salt and PEG6000-simulated drought stress treatment, seedlings were treated with 100 mM NaCl, 20% PEG6000, and water as control after one-day adaptation to the hydroponic conditions. The leaves were collected at 0, 12, 24, 48, and 72 h after treatments. For low temperature stress, seedlings were transferred to 4 ℃ and those grown at 24 ℃ served as a control group. The leaves were collected at 0, 6, 12, 24, and 72 h after treatment. All samples were immediately frozen in liquid nitrogen and transferred to − 80 ℃ until analysis.

Identification, chromosomal location, and gene duplication of BURP genes in jujube

The BURP protein sequences of Arabidopsis thaliana were downloaded from Swissprot database ( and used as queries for local blast analysis. Jujube genome dataset was obtained from our previous publication [47]. The blast output results were further aligned by Protein BLAST in NCBI ( and confirmed for BURP domain using Batch Web CD-Search Tool ( Next, the Mw and pI of these validated proteins were predicted using ExPASy website ( [48]. The subcellular localization was predicted by CELLO v2.5 server ( [49]. Chromosomal distribution of ZjBURPs was visualized based on their physical location from jujube genome databases. Gene duplication events were determined by MCScanX [50, 51].

Sequence alignment of BURP genes and phylogenetic analysis

All ZjBURP protein sequences were multi-aligned using DNAMAN v. (Lynnon Biosoft, CA, USA) [52]. In addition, BURP proteins from jujube and several other plants, including Brachypodium distachyon, Setaria italica, Oryza sativa, Sorghum bicolor, Zea mays, Arabidopsis thaliana, Glycine max, Cucumis sativu, Citrus sinensis, Brassica rapa, Vicia faba, and Populus trichocarpa, were multi-aligned by MEGA-X [38, 53]. A phylogenetic tree was further constructed through the neighbor-joining (NJ) method of MEGA-X with Poisson model and 1,000 bootstrap replications, and the results were finally visualized by IToL v6.5.2 ( [54].

BURP sequence analysis in jujube

Distribution of predicted signal peptides and BURP domains in jujube genes were analyzed using SignalP 4.0 server ( [55] and the Batch Web CD-Search Tool (, respectively. ZjBURP gene structures were illustrated by Gene Structure Display Server 2.0 (GSDS: [56]. The online MEME program [57] was employed to detect the conserved motifs of ZjBURP proteins with following parameters: size distribution, zero or one occurrence per sequence; motif count, 15; and motif width, 6–50 amino acids. The cis-elements in the region of 3,000 bp upstream of the transcription initiation site of ZjBURP genes were predicted by PlantCARE database ( [58].

Transcript expression of ZjBURPs by RT-PCR and qRT-PCR

ZjBURP expressions in different organs, including flowers, young fruit, and leaves of jujube were analyzed by RT-PCR. Total RNA was isolated from samples using SteadyPure Plant RNA Extraction Kit (Accurate Biotechnology; Hunan, China). RNA quality and integrity was confirmed utilizing a NanoDrop 2000 UV–vis spectrophotometer (Thermo Fisher Scientific, USA) and 1% agarose gel electrophoresis, respectively. The RNA was reversed using Evo M-MLV RT Premix (Accurate Biotechnology). The gene-specific primers for 17 ZjBURPs and Zj26S-2 (reference gene) were designed by Primer Premier 5 (Table S2). RT-PCR was carried out using 2 × Rapid Taq Master Mix (Vazyme) on an ABI 2720 Thermal Cycler (ABI; Marsiling, Singapore), and PCR products were detected by 1% agarose electrophoresis.

In addition, the relative expression levels of ZjBURPs in wild jujube leaf in response to different stresses (low temperature, salt and drought stresses) were determined by qRT-PCR using SYBR® Green Premix Pro Taq HS qPCR Kit II AG11702 on a Roche LightCycler96. The protocol was set as follows: 95 °C for 30 s; 40 cycles of 95 °C for 5 s, 58 °C for 30 s, 72 °C for 30 s; and 95 °C for 10 s, 65 °C for 60 s, 97 °C for 1 s; 37 °C for 30 s. The 2−ΔΔCt method was employed to calculate the relative expression levels of ZjBURPs from three biological replicates after normalization to Zj26S-2 expression (GenBank Accession: NC_029685; forward primer: 5′-TGGCTGAAGAATTGGGCCTT-3′, reverse primer: 5′-AGCCAAAGGAACTGCCCTAC-3′) as the internal reference gene.

Availability of data and materials

The reference jujube genome dataset could be found in the National Center for Biotechnology Information (NCBI) with the accession GCA_001835785.2. All original analysis data have been uploaded as supplementary tables or figures that could be found online. Further detail datasets generated and analyzed could be available from the corresponding author on reasonable request.



BNM2, USP, RD22 and PG1β


Fragments per kilobase of transcript per million fragments


Quantitative real-time polymerase chain reaction


  1. Gong Z, Xiong L, Shi H, Yang S, Herrera-Estrella LR, Xu G, Chao DY, Li J, Wang PY, Qin F, et al. Plant abiotic stress response and nutrient use efficiency. Sci China Life Sci. 2020;63(5):635–74.

    Article  PubMed  Google Scholar 

  2. Zhu JK. Abiotic stress signaling and responses in plants. Cell. 2016;167(2):313–24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Phillips K, Ludidi N. Drought and exogenous abscisic acid alter hydrogen peroxide accumulation and differentially regulate the expression of two maize RD22-like genes. Sci Rep-Uk. 2017;7:8821.

    Article  Google Scholar 

  4. Jeon JS, Chung YY, Lee S, Yi GH, Oh BG, An GH. Isolation and characterization of an anther-specific gene, RA8, from rice (Oryza sativa L.). Plant Mol Biol. 1999;39(1):35–44.

    Article  CAS  PubMed  Google Scholar 

  5. Hattori J, Boutilier KA, Campagne MMV, Miki BL. A conserved BURP domain defines a novel group of plant proteins with unusual primary structures. Mol Gen Genet. 1998;259(4):424–8.

    Article  CAS  PubMed  Google Scholar 

  6. Ding X, Hou X, Xie K, Xiong L. Genome-wide identification of BURP domain-containing genes in rice reveals a gene family with diverse structures and responses to abiotic stresses. Planta. 2009;230(1):149–63.

    Article  CAS  PubMed  Google Scholar 

  7. Xu HL, Li YX, Yan YM, Wang K, Gao Y, Hu YK. Genome-scale identification of soybean BURP domain-containing genes and their expression under stress treatments. Bmc Plant Biol. 2010;10:197.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Gan D, Jiang H, Zhang J, Zhao Y, Zhu S, Cheng B. Genome-wide analysis of BURP domain-containing genes in maize and sorghum. Mol Biol Rep. 2011;38(7):4553–63.

    Article  CAS  PubMed  Google Scholar 

  9. Shao Y, Wei G, Wang L, Dong Q, Zhao Y, Chen B, Xiang Y. Genome-wide analysis of BURP domain-containing genes in Populus trichocarpa. J Integr Plant Biol. 2011;53(9):743–55.

    CAS  PubMed  Google Scholar 

  10. Sun H, Wei H, Wang H, Hao P, Gu L, Liu G, Ma L, Su Z, Yu S. Genome-wide identification and expression analysis of the BURP domain-containing genes in Gossypium hirsutum. BMC Genomics. 2019;20(1):558.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Kavas M, Yildirim K, Secgin Z, Abdulla MF, Gokdemir G. Genome-wide identification of the BURP domain-containing genes in Phaseolus vulgaris. Physiol Mol Biol Plants. 2021;27(9):1885–902.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Fu L, Zhang Z, Wang H, Zhao X, Su L, Geng L, Lu Y, Tong B, Liu Q, Jiang X. Genome-wide analysis of BURP genes and identification of a BURP-V gene RcBURP4 in Rosa chinensis. Plant Cell Rep. 2022;41(2):395–413.

    Article  CAS  PubMed  Google Scholar 

  13. Chitkara P, Poddar N, Singh A, Kumar S. BURP domain-containing genes in legumes: genome-wide identification, structure, and expression analysis under stresses and development. Plant Biotechnol Rep. 2022;16:369–388.

  14. Treacy BK, Hattori J, Prudhomme I, Barbour E, Boutilier K, Baszczynski CL, Huang B, Johnson DA, Miki BL. Bnm1, a Brassica pollen-specific gene. Plant Mol Biol. 1997;34(4):603–11.

    Article  CAS  PubMed  Google Scholar 

  15. Boutilier KA, Gines MJ, Demoor JM, Huang B, Baszczynski CL, Iyer VN, Miki BL. Expression of the Bnmnap subfamily of napin genes coincides with the induction of Brassica microspore embryogenesis. Plant Mol Biol. 1994;26(6):1711–23.

    Article  CAS  PubMed  Google Scholar 

  16. Teerawanichpan P, Xia Q, Caldwell SJ, Datla R, Selvaraj G. Protein storage vacuoles of Brassica napus zygotic embryos accumulate a BURP domain protein and perturbation of its production distorts the PSV. Plant Mol Biol. 2009;71(4–5):331–43.

    Article  CAS  PubMed  Google Scholar 

  17. Chesnokov YV, Meister A, Manteuffel R. A chimeric green fluorescent protein gene as an embryogenic marker in transgenic cell culture of Nicotiana plumbaginifolia Viv. Plant Sci. 2002;162(1):59–77.

    Article  CAS  Google Scholar 

  18. Bassuner R, Baumlein H, Huth A, Jung R, Wobus U, Rapoport TA, Saalbach G, Muntz K. Abundant embryonic mRNA in field bean (Vicia faba L.) codes for a new class of seed proteins: cDNA cloning and characterization of the primary translation product. Plant Mol Biol. 1988;11(3):321–34.

    Article  CAS  PubMed  Google Scholar 

  19. Wang AM, Xia Q, Xie WS, Datla R, Selvaraj G. The classical Ubisch bodies carry a sporophytically produced structural protein (RAFTIN) that is essential for pollen development. P Natl Acad Sci USA. 2003;100(24):14487–92.

    Article  CAS  Google Scholar 

  20. Batchelor AK, Boutilier K, Miller SS, Hattori J, Bowman LA, Hu M, Lantin S, Johnson DA, Miki BL. SCB1, a BURP-domain protein gene, from developing soybean seed coats. Planta. 2002;215(4):523–32.

    Article  CAS  PubMed  Google Scholar 

  21. Fernandez L, Torregrosa L, Terrier N, Sreekantan L, Grimplet J, Davies C, Thomas MR, Romieu C, Ageorges A. Identification of genes associated with flesh morphogenesis during grapevine fruit development. Plant Mol Biol. 2007;63(3):307–23.

    Article  CAS  PubMed  Google Scholar 

  22. Zheng L, Heupel RC, DellaPenna D. The β Subunit of tomato fruit polygalacturonase lsoenzyme 1: isolation, characterization, and identification of unique structural features. Plant Cell. 1992;4:10.

    Google Scholar 

  23. Watson CF, Zheng LS, Dellapenna D. Reduction of tomato polygalacturonase beta-subunit expression affects pectin solubilization and degradation during fruit ripening. Plant Cell. 1994;6(11):1623–34.

    CAS  PubMed  PubMed Central  Google Scholar 

  24. Harshavardhan VT, Van Son L, Seiler C, Junker A, Weigelt-Fischer K, Klukas C, Altmann T, Sreenivasulu N, Baumlein H, Kuhlmann M. AtRD22 and AtUSPL1, members of the plant-specific BURP domain family involved in Arabidopsis thaliana drought tolerance. PLoS ONE. 2014;9(10):e110065.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Banzai T, Sumiya K, Hanagata N, Dubinsky Z, Karube I. Molecular cloning and characterization of genes encoding BURP domain-containing protein in the mangrove. Bruguiera gymnorrhiza Trees. 2001;16(2–3):87–93.

    Google Scholar 

  26. Zhang Z, Shi QQ, Wang B, Ma AM, Wang YK, Xue QT, Shen BQ, Hamaila H, Tang T, Qi XQ, et al. Jujube metabolome selection determined the edible properties acquired during domestication. Plant J. 2022;109(5):1116–33.

    Article  CAS  PubMed  Google Scholar 

  27. Liu M, Wang J, Wang L, Liu P, Zhao J, Zhao Z, Yao S, Stanica F, Liu Z, Wang L, et al. The historical and current research progress on jujube-a superfruit for the future. Hortic Res. 2020;7:119.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Gao QH, Wu CS, Wang M. The jujube (Ziziphus jujuba Mill.) fruit: a review of current knowledge of fruit composition and health benefits. J Agric Food Chem. 2013;61(14):3351–63.

    Article  CAS  PubMed  Google Scholar 

  29. Chen X, Chen R, Wang Y, Wu C, Huang J. Genome-wide identification of WRKY transcription factors in Chinese jujube (Ziziphus jujuba Mill.) and their involvement in fruit developing, ripening, and abiotic stress. Genes (Basel). 2019;10(5):360.

    Article  CAS  PubMed Central  Google Scholar 

  30. Li M, Hou L, Liu S, Zhang C, Yang W, Pang X, Li Y. Genome-wide identification and expression analysis of NAC transcription factors in Ziziphus jujuba Mill. reveal their putative regulatory effects on tissue senescence and abiotic stress responses. Ind Crop Prod. 2021;173:114093.

    Article  CAS  Google Scholar 

  31. Zhang Z, Li XG. Genome-wide identification of AP2/ERF superfamily genes and their expression during fruit ripening of Chinese jujube. Sci Rep-Uk. 2018;8:15612.

    Article  Google Scholar 

  32. Wang L, Wang L, Gao M, Qi C, Yang J, Li M, Ji S, Liu Z, Zhang M, Liu M. Genome-wide identification of IQ67 domain (IQD) gene families in Chinese jujube (Ziziphus jujuba Mill.) and expression profiles in response to cold stress. Scientia Horticulturae. 2022;293:110686.

    Article  CAS  Google Scholar 

  33. Wang L, Li M, Liu Z, Dai L, Zhang M, Wang L, Zhao J, Liu M. Genome-wide identification of CNGC genes in Chinese jujube (Ziziphus jujuba Mill.) and ZjCNGC2 mediated signalling cascades in response to cold stress. BMC Genomics. 2020;21(1):191.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Yamaguchishinozaki K, Shinozaki K. The plant hormone abscisic-acid mediates the drought-induced expression but not the seed-specific expression of rd22, a gene responsive to dehydration stress in Arabidopsis thaliana. Mol Gen Genet. 1993;238(1–2):17–25.

    Article  CAS  Google Scholar 

  35. Li Y, Chen X, Chen Z, Cai R, Zhang H, Xiang Y. Identification and expression analysis of BURP domain-containing genes in Medicago truncatula. Front Plant Sci. 2016;7:485.

    PubMed  PubMed Central  Google Scholar 

  36. Shen LY, Luo H, Wang XL, Wang XM, Qiu XJ, Liu H, Zhou SS, Jia KH, Nie S, Bao YT, et al. Chromosome-scale genome assembly for Chinese sour jujube and insights into its genome evolution and domestication signature. Front Plant Sci. 2021;12:773090.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Liu MJ, Zhao J, Cai QL, Liu GC, Wang JR, Zhao ZH, Liu P, Dai L, Yan G, Wang WJ, et al. The complex jujube genome provides insights into fruit tree biology. Nat Commun. 2014;5:5315.

    Article  CAS  PubMed  Google Scholar 

  38. Wang L, Wu N, Zhu Y, Song W, Zhao X, Li Y, Hu Y. The divergence and positive selection of the plant-specific BURP-containing protein family. Ecol Evol. 2015;5(22):5394–412.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Matus JT, Aquea F, Espinoza C, Vega A, Cavallini E, Dal Santo S, Canon P, Rodriguez-Hoces de la Guardia A, Serrano J, Tornielli GB, et al. Inspection of the grapevine BURP superfamily highlights an expansion of RD22 genes with distinctive expression features in berry development and ABA-mediated stress responses. PLoS One. 2014;9(10):e110372.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Van Son L, Tiedemann J, Rutten T, Hillmer S, Hinz G, Zank T, Manteuffel R, Baumlein H. The BURP domain protein AtUSPL1 of Arabidopsis thaliana is destined to the protein storage vacuoles and overexpression of the cognate gene distorts seed development. Plant Mol Biol. 2009;71(4–5):319–29.

    Article  CAS  PubMed  Google Scholar 

  41. Xu B, Gou JY, Li FG, Shangguan XX, Zhao B, Yang CQ, Wang LJ, Yuan S, Liu CJ, Chen XY. A cotton BURP domain protein interacts with alpha-expansin and their co-expression promotes plant growth and fruit production. Mol Plant. 2013;6(3):945–58.

    Article  CAS  PubMed  Google Scholar 

  42. Park J, Cui Y, Kang BH. AtPGL3 is an Arabidopsis BURP domain protein that is localized to the cell wall and promotes cell enlargement. Front Plant Sci. 2015;6:412.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Liu HH, Ma Y, Chen N, Guo SY, Liu HL, Guo XY, Chong K, Xu YY. Overexpression of stress-inducible OsBURP16, the beta subunit of polygalacturonase 1, decreases pectin content and cell adhesion and increases abiotic stress sensitivity in rice. Plant Cell Environ. 2014;37(5):1144–58.

    Article  CAS  PubMed  Google Scholar 

  44. Jin J, Duan J, Shan C, Mei Z, Chen H, Feng H, Zhu J, Cai W. Ethylene insensitive3-like2 (OsEIL2) confers stress sensitivity by regulating OsBURP16, the beta subunit of polygalacturonase (PG1beta-like) subfamily gene in rice. Plant Sci. 2020;292:110353.

    Article  CAS  PubMed  Google Scholar 

  45. Abe H, YamaguchiShinozaki K, Urao T, Iwasaki T, Hosokawa D, Shinozaki K. Role of Arabidopsis MYC and MYB homologs in drought- and abscisic acid-regulated gene expression. Plant Cell. 1997;9(10):1859–68.

    CAS  PubMed  PubMed Central  Google Scholar 

  46. Abe H, Urao T, Ito T, Seki M, Shinozaki K, Yamaguchi-Shinozaki K. Arabidopsis AtMYC2 (bHLH) and AtMYB2 (MYB) function as transcriptional activators in abscisic acid signaling. Plant Cell. 2003;15(1):63–78.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Huang J, Zhang C, Zhao X, Fei Z, Wan K, Zhang Z, Pang X, Yin X, Bai Y, Sun X, et al. The jujube genome provides insights into genome evolution and the domestication of sweetness/acidity taste in fruit trees. PLoS Genet. 2016;12(12):e1006433.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Artimo P, Jonnalagedda M, Arnold K, Baratin D, Csardi G, de Castro E, Duvaud S, Flegel V, Fortier A, Gasteiger E, et al. ExPASy: SIB bioinformatics resource portal. Nucleic Acids Res. 2012;40(W1):W597–603.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Yu CS, Lin CJ, Hwang JK. Predicting subcellular localization of proteins for Gram-negative bacteria by support vector machines based on n-peptide compositions. Protein Sci. 2004;13(5):1402–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Wang YP, Tang HB, DeBarry JD, Tan X, Li JP, Wang XY, Lee TH, Jin HZ, Marler B, Guo H, et al. MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 2012;40(7):e49.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Chen CJ, Chen H, Zhang Y, Thomas HR, Frank MH, He YH, Xia R. TBtools: An integrative toolkit developed for interactive analyses of big biological data. Mol Plant. 2020;13(8):1194–202.

    Article  CAS  PubMed  Google Scholar 

  52. Ma YP, Han YR, Feng XR, Gao HD, Cao B, Song LH. Genome-wide identification of BAM (beta-amylase) gene family in jujube (Ziziphus jujuba Mill.) and expression in response to abiotic stress. BMC Genomics. 2022;23(1):438.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Kumar S, Stecher G, Li M, Knyaz C, Tamura K. MEGA X: Molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. 2018;35(6):1547–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Letunic I, Bork P. Interactive Tree Of Life (iTOL) v5: an online tool for phylogenetic tree display and annotation. Nucleic Acids Res. 2021;49(W1):W293–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Petersen TN, Brunak S, von Heijne G, Nielsen H. SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Methods. 2011;8(10):785–6.

    Article  CAS  PubMed  Google Scholar 

  56. Hu B, Jin JP, Guo AY, Zhang H, Luo JC, Gao G. GSDS 2.0: an upgraded gene feature visualization server. Bioinformatics. 2015;31(8):1296–7.

    Article  PubMed  Google Scholar 

  57. Bailey TL, Williams N, Misleh C, Li WW. MEME: discovering and analyzing DNA and protein sequence motifs. Nucleic Acids Res. 2006;34:W369–73.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Lescot M, Dehais P, Thijs G, Marchal K, Moreau Y, Van de Peer Y, Rouze P, Rombauts S. PlantCARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences. Nucleic Acids Res. 2002;30(1):325–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references


Not applicable.


The research was supported by the National Key Research and Development Program of China (Grant No. 2018YFD1000607) and National Natural Science Foundation of China (Grant No. 32001343).

Author information

Authors and Affiliations



XL and ZZ conceived and designed the experiments. WW performed experiments, data analysis, and wrote the manuscript. ZZ and XL contributed to manuscript preparation. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Xingang Li.

Ethics declarations

Ethics approval and consent to participate

The plant material of this study protocol complies with relevant institutional, national, and international guidelines and legislation. All experimental protocols were approved by Northwest A&F University.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Supplementary Information

Additional file 1: Fig S1.

The conserved BURP domain and signal peptide of BURP proteins in jujube. The phylogenetic tree of 17 BURP proteins is shown using MEGA-X with the neighbor-joining (NJ) method. The red filled boxed represent the BURP domain, and the green filled boxes represent the signal peptide.

Additional file 2: Fig S2.

The conserved amino acid residues in BURP domain of jujube sequences. Multialignment of 17 ZjBURP proteins using DNAMAN v. The black-, pink-, and aquamarine-shaded amino acids represented the homology level with 100%, over 75%, and over 50%, respectively

Additional file 3: Table S1.

Number of cis-elements in putative promoters of ZjBURPs.

Additional file 4:

 Fig S3. The original, uncropped gel of RT-PCR analysis of ZjBURPs in three tissues. The red box denoted the region of the original gel. The bands in the F, YF, and L row represented ZjBURPs expression in flower, young fruit, and leaf, respectively. 1-9: ZjBNM1, 2, 3, 4, 5, 6, 7, 8, and 9; 10-14: ZjBURP1, 2, 3, 4, and 5; 15-17: ZjPG1, 2, and 3; 18: Zj26S-2.

Additional file 5: Table S2.

Primer pairs used in RT-PCR and qRT-PCR analysis.

Rights and permissions

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

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Wang, W., Zhang, Z. & Li, X. Identification and expression analysis of BURP domain-containing genes in jujube and their involvement in low temperature and drought response. BMC Genomics 23, 692 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Ziziphus jujuba Mill
  • BURP domain-containing protein
  • Gene expression
  • Abiotic stress
  • Fruit development