Skip to main content

Evolutionary history of the C-repeat binding factor/dehydration-responsive element-binding 1 (CBF/DREB1) protein family in 43 plant species and characterization of CBF/DREB1 proteins in Solanum tuberosum



Plants are easily affected by temperature variations, and high temperature (heat stress) and low temperature (cold stress) will lead to poor plant development and reduce crop yields. Therefore, it is very important to identify resistance genes for improving the ability of plants to resist heat stress or cold stress by using modern biotechnology. Members of the C-repeat binding factor/Dehydration responsive element-binding 1 (CBF/DREB1) protein family are related to the stress resistance of many plant species. These proteins affect the growth and development of plants and play vital roles during environmental stress (cold, heat, drought, salt, etc.). In this study, we identified CBF/DREB1 genes from 43 plant species (including algae, moss, ferns, gymnosperms, angiosperms) by using bioinformatic methods to clarify the characteristics of the CBF/DREB1 protein family members and their functions in potato under heat and cold stresses.


In this study, we identified 292 CBF/DREB1 proteins from 43 plant species. However, no CBF/DREB1 protein was found in algae, moss, ferns, or gymnosperms; members of this protein family exist only in angiosperms. Phylogenetic analysis of all the CBF/DREB1 proteins revealed five independent groups. Among them, the genes of group I do not exist in eudicots and are found only in monocots, indicating that these genes have a special effect on monocots. The analysis of motifs, gene duplication events, and the expression data from the PGSC website revealed the gene structures, evolutionary relationships, and expression patterns of the CBF/DREB1 proteins. In addition, analysis of the transcript levels of the 8 CBF/DREB1 genes in potato (Solanum tuberosum) under low-temperature and high-temperature stresses showed that these genes were related to temperature stresses. In particular, the expression levels of StCBF3 and StCBF4 in the leaves, stems, and roots significantly increased under high-temperature conditions, which suggested that StCBF3 and StCBF4 may be closely related to heat tolerance in potato.


Overall, members of the CBF/DREB1 protein family exist only in angiosperms and plays an important role in the growth and development of plants. In addition, the CBF/DREB1 protein family is related to the heat and cold resistance of potato. Our research revealed the evolution of the CBF/DREB1 family, and is useful for studying the precise functions of the CBF/DREB1 proteins when the plants are developing and are under temperature stress.


Plants are continuously affected by various environmental factors such as temperature variation, high salinity, and drought during their growth process. These stresses can lead to many undesirable consequences, such as hindered plant growth and reduced crop yields. Plant responses to abiotic stress are complex and controlled by a variety of mechanisms, including physiological, biochemical, and molecular regulator mechanisms. Among them, molecular regulation may be related to stress perception, signal transduction, gene expression and final metabolic modification to prevent cell damage under abiotic stress [1].

Many genes related to plant stress resistance, such as those related to membrane stability [2], antifreeze protein genes [3], osmo-regulatory genes [4], and antioxidant enzyme activity genes [5], have been isolated, cloned and identified in plants. However, stress resistance is a quantitative trait controlled by multiple genes in plants, and it is difficult to cause significant effects by transferring only a single stress gene to plants [6]. Therefore, the study of transcription factors (TFs) that can improve the stress resistance of plants by regulating the expression of multiple stress-responsive genes has become a popular research.

The AP2/ERF superfamily is one of the largest TF families, and its members are involved in the regulation of growth and development and biotic and abiotic stress responses in plants [7, 8]. This family includes all the genes encoding at least one APETALA2 (AP2) domain and can be further divided into ethylene response factor (ERF), AP2, RAV (Related to ABI3/VP1), and soloists. Among these, the ERF family can be further divided into two subfamilies based on the amino acid sequence of their DNA binding domain: the CBF/DREB (C-repeat binding factor/Dehydration-responsive element-binding) subfamily (group A) and the ERF subfamily (group B) [9]. The DREB subfamily consists of six subgroups, of which the A1 subgroup includes CBF/DREB1 transcription factors [10]. The CBF/DREB1 gene is a type of plant-specific transcription factor that was first discovered in Arabidopsis thaliana (Arabidopsis) [11]. Members of the CBF/DREB1 family regulate the expression of most stress-inducible genes that are independent of the ABA pathway, and binds specifically to the CRT/DRE sequence (TACCGACAT) within the promoters of stress-inducible genes to activate the expression of the target genes, improving the resistance of plants to stress [12].

Members of the CBF/DREB1 subfamily exhibit different response patterns under various abiotic stresses, including low-temperature, drought, and high salinity [13,14,15,16], and transgenic methods using these genes have shown that these factors play important roles in abiotic stress signalling and tolerance [12]. In Arabidopsis, six CBF/DREB1 genes have been characterized [17]. Among them, AtCBF1, AtCBF2, and AtCBF3 are induced by cold stress, whereas AtCBF4, AtDREB1E, and AtDREB1F are induced by osmotic stresses such as drought and salt [10, 18,19,20]. The cold-resistance of tobacco (Nicotiana tabacum) and Arabidopsis can be enhanced because of the increased content of proline and soluble sugars by overexpressing the AtCBF3 gene [21, 22]. In rice (Oryza sativa), OsDREB1A is induced by salt and cold stresses [23, 24], whereas OsDREB1B and OsDREBL are induced by cold stress [25]. In Arachis hypogaea, PNDREB1 is strongly upregulated under low-temperature conditions [26]. The expression of PpCBF1 is also induced by NaCl, cold, drought, and ABA treatments in Physcomitrella patens [27]. Under cold stress, BjDREB1B and BoCBF/DREB1 are highly expressed in Brassica juncea and Brassica oleracea (kale), respectively [28, 29]. Overexpressing AtCBF3 in maize not only enhances cold tolerance, but also enhances drought and salt tolerance [30]. In transgenic Arabidopsis, salt tolerance and freezing tolerance are enhanced by overexpressing JcDREB from the woody biodiesel plant Jatropha curcas [31]. Overexpression of CkDREBa from Caragana korshinski confers drought and salinity tolerance to transgenic chrysanthemum [32]. MbDREB1 from dwarf apple (Malus baccata) increases plant tolerance to low temperature, drought, and salt stress via both ABA-dependent and ABA-independent pathways [33]. Transgenic Arabidopsis and rice plants overexpressing OsDREB1A are tolerant to low temperature, high salinity, and drought [13, 23]. Overexpression of CBF1/DREB1B in tomato (Solanum lycopersicum) enhanced the chilling and drought tolerance of transgenic tomato plants [34, 35].

Since release of the whole-genome sequences of increased numbers of plant species, many CBF/DREB1 transcription factors have been identified in species such as Arabidopsis, rice, wheat (Triticum aestivum), barley (Hordeum vulgare), and soybean (Glycine max) [10, 36,37,38,39,40]. However, no researchers have described a detailed and complete evolutionary history of these genes. In this study, we identified CBF/DREB1 genes from the complete genome sequences of 43 plant species (including lower plants and higher plants, such as algae, moss, ferns, gymnosperms, angiosperms) by using bioinformatic methods, such as analysis of phylogenetic relationships, motifs, gene duplications, gene characteristics and gene structure. In addition, potato is one of the most important crops species in the world, ranking fourth after rice, wheat and maize, and plays a significant economic role. As the main research object of our research team, we analysed Gene Ontology (GO) annotations and transcriptional profiles of the CBF/DREB1 proteins in potato. Based on the results, the evolution of the CBF/DREB1 genes can be revealed in high detail, and some useful information can be gained for further research on potato.


Identification of CBF/DREB1 genes in 43 plant species

To study the evolution of the CBF/DREB1 family, we used the CBF/DREB1 full-length protein sequence of Arabidopsis [17] as a query to search the available protein sequences from 42 other species, including algae, mosses, ferns, gymnosperms, and angiosperms. There are two conserved amino acid sequences, PKK/RPAGRxKFxETRHP and DSAWR upstream and downstream, respectively, of the AP2 domains of the CBF/DREB1 genes [41] (Fig. 1). Based on the sequence structural features and the domains of the CBF/DREB1 family, we removed incorrect and redundant sequences prior to further analysis of the remaining CBF/DREB1 sequences that could be considered putative genes. We ultimately identified a total of 292 proteins in 43 plant species and named them accordingly (Additional file 1a).

Fig. 1

Amino acid sequence alignment of CBF/DREB1 proteins from Arabidopsis thaliana. The characteristic domains such as the PKK/RPAGRxKFxETRHP, DSAWR, and AP2 domains are indicated in the figure

The number of CBF/DREB1 proteins ranged from 0 to 18 (Table 1). Medicago truncatula (Mt) and Carica papaya (Cp) contained 18 and 2 CBF/DREB1 proteins, respectively, which were greater and fewer than those in other plant species. Moreover, we obtained the genome length of 43 species from the NCBI database (Additional file 1a) and found that the number of CBF/DREB1 proteins was not related to the genome length.

Table 1 CBF/DREB1 proteins identified in 43 sequenced plant genomes

Analysis of conserved motifs and gene duplication events

Prediction of protein motifs is an essential component of protein analysis. In this study, we used the MEME online tool to identify the motifs of 292 CBF/DREB1 proteins (the picture, which is not shown here, is too large to show clearly in the manuscript). (Additional files 2, 3). Twenty putative motifs were identified in these proteins. The motifs of the CBF/DREB1 family members were found to be conserved. The vast majority of proteins (more than 180 proteins) contained 8 or 9 motifs, and they were arranged in the same order. More than 80 proteins had 7 or 10 motifs, and a only few proteins had 5, 6 or 11 motifs. Additionally, OsCBF7 was the only protein with 4 motifs in rice (Additional file 1a).

To further understand the evolutionary relationship of CBF/DREB1 members, predicting of the molecular evolution rate can be used to clarify the gene evolution process [42]. In this study, we identified the tandem and segmental duplications in the 43 species (Table 1, Additional file 1 a, b), and the gene duplication events were identified in 28 species, based on the available information (in fact, due to incomplete information on gene location and on the PGDD website, we could not identify all the gene duplication events). In general, there were fewer tandem duplications than segmental duplication events. In addition, there were more CBF/DREB1 gene duplications in eudicots than in monocots. Most species only had tandem or segmental duplications, but some species had both. Tandem duplications did not occur in soybean, but we identified 20 segmental duplications, which was the highest number among all the proteins. Eucalyptus grandis contained the most tandem duplication events-9 pairs. We then calculated the Ka/Ks value for each duplication event (Additional file 1b), and the results showed that nearly every gene pair evolved with a Ka/Ks value less than one. However, in Lactuca sativa, the only gene pair had a Ka/Ks value of 1.1288.

Analysis of the phylogenetic tree

A phylogenetic tree of the CBF/DREB1 members in the 43 plant species was constructed on the basis of the similarities of their protein sequences (Fig. 2). As shown in Fig. 2, the 292 identified CBF/DREB1 proteins from angiosperms could be roughly divided into five groups based on the phylogenetic tree. We analysed the distribution of the CBF/DREB1 proteins of the 43 species in different subgroups and found that the more closely related the species, the more similar the protein classification was (Table 1). Six CBF/DREB1 proteins in Arabidopsis were distributed in group IV and group V, which is consistent with the findings of previous studies [43, 44]. Group I and group IV had more members than did the other groups, each with more than 70 proteins. Group II contained the fewest CBF/DREB1 protein members. The CBF/DREB1 proteins in monocots are found only in groups I and II, and no group I CBF/DREB1 proteins were identified in eudicots. This suggested that the proteins of group I may have a unique role in monocots.

Fig. 2

Neighbour-joining tree comprising CBF/DREB1 genes of 43 plant species. Bootstrap values were calculated in 1000 replications by using MEGA 7. The analysed CBF/DREB1 genes were distributed in five main groups (including Group I, Group II, Group III, Group IV and Group V), which were marked with different colours. Some CBF/DREB1 families are indicated with different marks. The black star indicates rice (Oryza sativa), the white stat indicates Arabidopsis, the black triangle indicates tomato (Solanum lycopersicum), and the black tick indicates potato (Solanum tuberosum). At the nodes, the gray dots indicate that the bootstrap values were less than or equal to 40; the yellow dots show that the bootstrap values were between 40 and 80, and the red dots represent bootstrap values that were greater than 80 but less than or equal to 100

Gene Ontology annotations and analysis of the RNA-sequencing data of StCBF proteins

To better understand the biological processes affected by CBF/DREB1 proteins, we selected the potato CBF/DREB1 proteins for GO analysis using the NCBI database [42] (Fig. 3, Additional file 4). In terms of molecular functions and biological processes, all eight genes had identical functions. These genes not only were involved in 7 molecular functions (such as DNA-binding transcription factor activity and transcription regulator activity), but also were involved in nearly 50 biological processes, including regulation of metabolic processes, RNA biosynthetic processes, regulation of cellular biosynthetic processes, etc. (Fig. 3a, b). This indicated that the StCBF genes play an important role in potato growth and development. In addition, predictions of the cellular components of the genes showed that, in addition to StCBF4, seven other genes were involved in component development of intracellular, organelle, membrane-bounded organelle, and so on. StCBF5 and StCBF6 appeared to be involved in the component development of the membrane, and StCBF5 was also involved in the integral and intrinsic components’ development of membrane (Fig. 3c, Additional file 4).

Fig. 3

Information from Gene Ontology (GO) annotations. a Molecular function. b Biological process. c Cellular component. The numbers on the abscissa indicate the number of predicted proteins

We processed the RNA-seq database information and generated a heatmap (Figs. 4, 5). Figure 4 shows that the expression levels of the StCBF genes in all stages of potato development were not very high. Moreover, StCBF1, StCBF5, StCBF6 and StCBF8 were expressed in the petiole, and StCBF1 was clearly expressed in the root and tuber sprout. StCBF6 and StCBF7 were distributed in the tuber peel, and the expression of StCBF7 was detected in the mature tubers. The plants were treated with 10 different conditions, including phytohormone treatments (6-benzylaminopurine, BAP; indole-3-acetic acid, IAA; abscisic acid, ABA; gibberellic acid, GA3); abiotic stress (heat, drought, and salinity); and biotic stress (pathogen; DL-β-amino-n-butyric acid, BABA; acibenzolar-S-methyl, BTH) [45]. Under the 10 different treatments, most genes showed downregulated expression or little change in expression. However, in response to the BAP, IAA, GA3 NaCl (salinity), mannitol (drought) and BABA treatments, the expression of StCBF6 increased significantly. StCBF1 also increased slightly under drought, and salinity stresses. Under heat stress, no genes were upregulated, and few genes were upregulated in response to the other abiotic stress treatments (salt and drought).

Fig. 4

Heatmap based on the RNA-seq information of 14 different tissues of potato. The RNA-seq information was not processed. The heatmap shows CBF/DREB1 gene expression across 14 different tissues throughout the entire potato life cycle: including flower, leaf, stolon, petiole, root, young tuber, mature tuber, shoot apex, tuber peel, tuber cortex, tuber sprout, tuber pith, and stamen tissues. In the heat map, high expression is shown in red and low expression is shown in black

Fig. 5

Expression profiles of 8 CBF/DREB1 genes under 10 different conditions in potato. Abiotic stress, biotic stress, and phytohormone treatment data and control data were obtained from the PGSC database. Transcripts were measured via RNA-Seq technology. The stress and hormone data were compared with the control data. In the heat map, upregulated expression is shown in red and downregulated expression is shown in green

CBF/DREB1 gene expression profiles in response to temperature variations

To demonstrate the change in CBF/DREB1 protein accumulation in different tissues (leaf, stem, and root tissues) in potato under low-temperature (4 °C) and high-temperature (35 °C) treatments, qRT-PCR analysis was used to conduct expression analysis (Fig. 6, Additional file 5). We mainly focused on the CBF/DREB1 genes of potato.

Fig. 6

CBF/DREB1 gene expression in potato leaves, stems, and roots at 35 °C and 4 °C. a Leaf, b Stem, c Root, d Comparison of gene expression in the three tissues based on that in the leaves. The X-axes show 8 representative CBF/DREB1 genes, and the y-axes show scales of relative expression levels. The Ef1a gene was used as a reference transcript. Leaf, stem, and root tissues were sampled from the same parts of control and experimental plants. The quantitative data were measured by taking three biological replicates and two technical replicates, and the relative expression level of each gene was calculated using the 2−ΔΔCt method

In the leaves (Fig. 6a), under low-temperature and high-temperature stresses, the expression levels of StCBF1, StCBF2, StCBF5, StCBF6, and StCBF8 did not change much, but the expression levels of StCBF3 and StCBF4 increased significantly. The expression of StCBF7 increased under high-temperature conditions. In the stems (Fig. 6b), the transcript levels of all genes decreased under low-temperature stress. However, notable upregulation of 4 genes (StCBF1-4) occurred under high temperature conditions, and the expression of StCBF7 decreased slightly, the expression of the other three genes did not significantly increase. In the roots (Fig. 6c), the expression levels of 8 genes increased significantly under high temperature. Under low-temperature conditions, the expression levels of StCBF4 and StCBF7 obviously increased, while the expression levels of StCBF3 barely increased, and the expression levels of other genes showed a decreased.

In addition, on the basis of the expression of genes in the leaves under normal conditions, we compared the expression levels of the 8 aforementioned genes in the 3 tissues (Fig. 6d). Except for StCBF7, the expression levels of 7 genes in the stems were higher than those in the leaves, and the transcript levels of StCBF4 and StCBF7 in the roots were lower than those in the leaves. Taken together, the results showed that StCBF genes in potato were expressed more in the roots and stems and less in the leaves, which is consistent with the data obtained from the PGSC website (Additional file 6).


CBF/DREB1 proteins have been identified in many species, such as wheat, rice, Chinese cabbage (Brassica rapa), and kale [43, 44, 46,47,48,49,50]. We combined these results, and we searched the CBF/DREB1 family again using the methods described in this article. For example, we found that Bra010463 (BrDREB1C1) [43] did not contain the PKK/RPAGRxKFxETRHP signature sequence, so we omitted it. In addition, we found that no CBF/DREB1 proteins were identified in algae, mosses, ferns, and the gymnosperms Picea abies. CBF/DREB1 proteins were identified only in angiosperms. According to the evolutionary relationships, the CBF/DREB1 protein did not evolve until the appearance of angiosperms.

The motifs of 292 identified proteins were analysed by the MEME website (Additional file 2). Previous studies have shown that certain conserved motifs to play functional and/or structural roles in active proteins [51]. Most proteins had 7–10 motifs, but OsCBF7 (LOC_Os06g06970.1, OsDREB1D) had only 4 motifs, which was significantly different from the numbers in the other proteins. Research has shown that OsCBF7 is a unique DREB homologous gene, and its expression cannot be detected with or without any stress in rice [23]. Yeast one-hybrid experiments showed that OsDREB1D not only can form complexes with proteins that have dehydration response elements/C-repeat motifs, but also can be combined with low-temperature response element (LTRE) sequence. When OsDREB1D was overexpressed in Arabidopsis thaliana, the transgenic plants obtained the cold and high-salt tolerance, and were also insensitive to ABA. These results indicated that the functions of OsDREB1D may be similar functions to those of other CBF/DREB1 transcription factors, and that the expression of OsDREB1D in rice may be controlled by a special mechanism for the redundancy of function [52].

Genomic duplications are the essential force driving the origin and evolution of species. However, across evolutionary history, most of these duplications have disappeared or have been silenced, with the remaining few playing a strong role in positive or purifying selection [53]. The substitution rate, known as the Ka/Ks (non-synonymous/synonymous) ratio, is an indicator for estimating selective pressure [54], and is commonly used to analyse the molecular evolutionary rate, and to investigate the evolutionary direction and selective strength of a coding sequence [55]. In this study, nearly every gene pair was subject to purification selection (Ka/Ks < 1), and only one gene pair had a Ka/Ks value greater than 1 in L. sativa, experienced positive selection pressure in evolution (Additional file 1b) [53]. The number of CBF/DREB1 proteins and duplication events suggested that the increase in CBF/DREB1 proteins may be due to numerous duplication events in angiosperms. Plants evolve with changes in the environment, and positive selection promotes the exchange of gene functions to survive. Therefore, selective patterns can partially explain the evolutionary patterns of genes [42].

The use of phylogenetic tree is an established method for the examination of the structure and function of a protein family, and the functional relationship can be understood by the phylogeny [56]. We specifically labelled the CBF/DREB1 proteins of potato, tomato, Arabidopsis and rice. Arabidopsis and rice are common model crop species used in research (Fig. 2). Tomato and potato are related species, and potato is our main research species. The role of a gene in potato can be inferred based on the function of that gene in A. thaliana, in tomato, or in rice. For example, overexpression of the AtCBF1 gene in Arabidopsis can induce the expression of the COR gene and increase the freezing tolerance of non-acclimated Arabidopsis plants [16]. Overexpression of the AtCBF3 gene of Arabidopsis in potato can improve the heat resistance of potato by regulating the expression of many genes involved in resistance, increasing the accumulation of soluble substances, improving the photosynthetic performance and antioxidant capacity, and reducing and eliminating excess ROS [57]. In rice, when plants are exposed to low temperatures, the transcripts of OsCBF2, OsCBF8 and OsCBF9 in the roots increase rapidly, and the expression level of OsCBF2 in the leaves increases under salt stress, suggesting that these genes may be involved in the low-temperature response or salt-stress response [50]. As a result, these CBF/DREB1 proteins in potato and tomato might have a similar function or a similar effect on development and responses to various stresses and stimuli.

All BoCBF genes in kale respond to abiotic stresses [49]. In addition, all CBF/DREB1 genes in Chinese cabbage are highly upregulated during cold (4 °C) treatment, and some of them are also responsive to salt (250 mM NaCl), drought (air drying), and ABA (100 mM) treatments [43]. GmDREB1B;1 can activate the expression of numerous soybean-specific stress-responsive genes and enhance ABRE (ABA response element)-mediated gene expression in an ABA-independent manner under a variety of abiotic stress conditions [58]. Therefore, we can also speculate on the function of the corresponding CBF/DREB1 gene in various species based on the existing results, which can provide a clear idea and direction for research.

On the basis of the data obtained from the PGSC database, we measured and analysed the expression levels of the 8 CBF/DREB1 genes at different stages of potato development and under various stresses [59]. The data related to the CBF/DREB1 genes from the PGSC database were analysed in double haploid potato material (DM) or heterozygous diploid material (RH) [45]. The change in the expression levels of the CBF/DREB1 proteins in potato under abiotic stress was different from the results of other studies on CBF/DREB1 genes [43, 49, 50, 57, 58]. The different results may be related to the different species or genotypes uesd. The 10 treatments also did not include cold stress. Therefore, to better understand the response of CBF/DREB1 proteins in potato to temperature stress, we analysed the expression patterns of the eight aforementioned genes of potato under cold (4 °C) and heat (35 °C) stresses using the tetraploid cultivar “Desiree”.

The expression profile analysis showed that under temperature stress, the expression levels of all 8 CBF/DREB1 proteins in potato changed, and the change in transcript levels was also different in different tissues. However, Fig. 6 shows that all genes responded to temperature variations in the leaves, stems, and roots. Interestingly, the expression levels of StCBF3 and StCBF4 significantly increased in the leaves, stems, and roots at 35 °C, which indicated that StCBF3 and StCBF4 were sensitive to high temperature and are worthy of further study.

The DREB subfamily consists of six subgroups [10]. In previous studies, it was found that the StCBF1 gene of the A1 subgroup in potato could improve the cold and freezing tolerance of Arabidopsis [60]. The StDREB2 gene of the A2 subgroup was confirmed to be involved in the plant response to several treatments (salt, ABA, etc.) [61]. The StDREB gene of the A4 subgroup plays an important role in the fungal invasion response, and can also improve salt tolerance in potato [62, 63]. Overexpression of the StDREB1 and StDREB2 genes in potato, neither of which is not a member of the A1 subgroup, improves oxidative stress damage by improving plant growth, and proline and antioxidant production, thereby increasing tolerance to cadmium (Cd) [64]. In addition, in other species, the CBF/DREB1 gene plays a very important role in improving plant stress resistance [28, 29, 31, 43, 57, 58]. However, to date, the CBF/DREB1 gene has not been studied in potato. This paper analysed the gene expression patterns of the members of the A1 subgroup of the DREB family in potato, providing an improved theoretical basis and data to further explore the heat tolerance of plants.


We performed an evolutionary analysis of the CBF/DREB1 protein family in the plant kingdom to reveal the gene structure, phylogenetic relationships, and evolution of the CBF/DREB1 proteins in each group. We identified 292 CBF/DREB1 proteins from 43 species, and CBF/DREB1 proteins were found only in angiosperms. Moreover, the genes of group I and group II were found in only monocots. All the genes of groups II-V were found in eudicots. Therefore, we can speculate that the function of these genes is essential for angiosperms. There were no genes of group I identified in eudicots, therefore, we conclude that group I genes play a crucial role in monocots. We also analysed the gene mapping, motif, gene duplication event, phylogenetic tree, GO annotation and expression data from the PGSC website. These results are important for understanding the properties and functions of CBF/DREB1 proteins, and provide a basis for studying the evolutionary relationships of CBF/DREB1 proteins in plants. Our analysis also showed that all eight CBF/DREB1 proteins in the leaves, stems, and roots of potato are related to high-temperature and low-temperature stresses. In particular, the expression levels of StCBF3 and StCBF4 in the leaves, stems, and roots were upregulated significantly under high-temperature conditions. These two genes may play an important role in the heat tolerance of potato. The above research results increase our understanding of the evolutionary relationships of the CBF/DREB1 family and lay a foundation for further functional identification of the CBF/DREB1 proteins in potato. In subsequent studies, we will try to determine the function of CBF/DREB1 proteins in terms of the biological processes that affect crop resistance and improve the response to various stresses through molecular biological techniques.


Identification of CBF/DREB1 Genes in Various Species

The genes, proteins, and coding sequences of 42 species were downloaded from the Phytozome (, PGSC (, and the NCBI ( database. Six Arabidopsis CBF/DREB1 protein sequences that were obtained from TAIR database ( were used as queries to perform a protein search against those of 42 species proteins in the database with a strict E value (< 1e−10) requirement by using BLAST 2.6.0 [65]. All candidate CBF/DREB1 protein sequences were screened by using the Conserved Domains Database (CDD) ( and were aligned using Molecular Evolutionary Genetics Analysis 7 (MEGA 7) software.

Sequence analysis of CBF/DREB1 genes/proteins

Information on the length of the protein sequences and chromosome location was obtained from the Phytozome, PGSC, NCBI, and other databases (Additional file 1). Multiple Em for Motif Elicitation website (MEME; was used to calculate the motifs of the protein sequences, with the number of motifs = 20.

Phylogenetic tree construction and gene duplication analyses

A phylogenetic tree comprising the members of the CBF/DREB1 family from 43 species was constructed using MEGA 7 software. The phylogenetic trees were produced using the neighbour-joining (NJ) method with the parameters of the Jones-Taylor-Thornton (JTT) model and 1000 replicates for bootstrap analysis [66]. The EvolView online tool ( evolview/#login) was used to draw and manage the phylogenetic trees.

Gene tandem duplication events of the CBF/DREB1 genes were analysed following the methods of Gu et al. [67]. The major criterion was that the length of an aligned sequence must cover > 75% of the longer gene, and the similarity of the aligned regions must be > 75%. Two genes were considered a tandem pair if they were located on the same chromosome and were separated by no more than 10 unrelated genes [68]. The segmental duplication pairs were analysed by using the Plant Genome Duplication Database (PGDD; Ka and Ks values were calculated by DNAsp software.

GO annotation and expression pattern analysis

Blast2GO software was used to analyse the Gene Ontology (GO) terms. The full-length amino acid sequences were uploaded into the program, and the NCBI database was chosen as the reference to for analysing the molecular functions, cellular components, and biological processes.

The expression data of the responses of CBF/DREB1 genes to different stresses were obtained from the PGSC database, and the transcript abundance was represented by using fragments per kilobase million (FPKM) values [69].

Plant materials and growth conditions

The tetraploid potato (Solanum tuberosum L.) cultivar “Diseree” was used in this study. Potato plantlets were cultured in MS media supplemented with 10 g/L agar and 20 g/L sucrose (pH 5.8) and were maintained at 25 ± 1 °C under 10,000 lx for 16 h and 20 ± 1 °C under 0 lx for 8 h. After 2 months, we selected seedlings of the same size for experimentation and divided them into three groups. The first group was subjected to a temperature of 35 °C for 24 h, the second group was subjected to a temperature of 4 °C for 24 h, and the third group was maintained under normal conditions as the control. We subsequently sampled the root, stem and leaf tissues (100 mg) from the same parts of the control plants and experimental plants. All of the samples were flash frozen in liquid nitrogen and then stored at − 80 °C prior to utilization.

RNA extraction and real-time qRT-PCR analysis of CBF/DREB1 gene expression

An RNA simple Total RNA Kit (Code No. RP3301, BioTeke, Beijing, China) was used to extract RNA. The elongation factor 1-a (Ef1a) gene was chosen as the reference gene [70]. Specific primers were designed using Primer Premier 5 software. First-strand cDNAs were synthesized from 1 μg of RNA with a PrimeScript™ RT reagent Kit (Code No. RR047A, TaKaRa, Dalian, China) in conjunction with gDNA Eraser in a 20 μL reaction volume. Real-time PCR was set up on the basis of a 2 × Plus SYBR real-time PCR mixture (Code No. PR7702, BioTeke, Beijing, China) and was performed on a QuantStudio™ 7 Flex Real-time PCR System (Applied Biosystems Inc., America) in a 10 μL reaction volume. The relative expression level of each gene was calculated using the 2−ΔΔCt method [71].

Availability of data and materials

All data generated or analyzed during this study are included in this article.



C-repeat binding factor/Dehydration responsive element-binding 1


Transcription factors




Ethylene response factor


Related to ABI3/VP1


Gene Ontology




Indole-3-acetic acid


Abscisic acid


Gibberellic acid


DL-β-amino-n-butyric acid




  1. 1.

    Agarwal PK, Agarwal P, Reddy MK, Sopory SK. Role of DREB transcription factors in abiotic and biotic stress tolerance in plants. Plant Cell Rep. 2006a;12(25):1263–74.

    Article  CAS  Google Scholar 

  2. 2.

    Kodama H, Nishiuchi T, Seo S, Ohashi Y, Iba K. Possible involvement of protein phosphorylation in the wound-responsive expression of Arabidopsis plastid ω-3 fatty acid desaturase gene. Plant Sci. 2000;2(155):153–60.

    Article  Google Scholar 

  3. 3.

    Graham LA, Liou YC, Walker VK, Davies PL. Hyperactive antifreeze protein from beetles. Nature. 1997;6644(388):727–8.

    Article  CAS  Google Scholar 

  4. 4.

    Ashraf M, Foolad MR. Roles of glycine betaine and proline in improving plant abiotic stress resistance. Environ Exp Bot. 2007;2(59):206–16.

    Article  CAS  Google Scholar 

  5. 5.

    Kornyeyev D, Logan BA, Payton P, Allen RD, Holaday AS. Enhanced photochemical light utilization and decreased chilling-induced photoinhibition of photosystem II in cotton overexpressing genes encoding chloroplast-targeted antioxidant enzymes. Physiol Plant. 2001;3(113):323–31.

    Article  Google Scholar 

  6. 6.

    Warren GJ. Cold stress: Manipulating freezing tolerance in plants. Curr Biol. 1998;15(8):514–6.

    Article  Google Scholar 

  7. 7.

    Mizoi J, Shinozaki K, Yamaguchi-Shinozaki K. AP2/ERF family transcription factors in plant abiotic stress responses. Biochim Biophys Acta. 1819;2011(2):86–96.

    Google Scholar 

  8. 8.

    Wessler SR. Homing into the origin of the AP2 DNA binding domain. Trends Plant Sci. 2005;2(10):54–6.

    Article  CAS  Google Scholar 

  9. 9.

    Sakuma Y, Liu Q, Dubouzet JG, Abe H, Shinozaki K, Yamaguchi-Shinozaki K. DNA-binding specificity of the ERF/AP2 domain of Arabidopsis DREBs, transcription factors involved in dehydration- and cold-inducible gene expression. Biochem Biophys Res Commun. 2002;3(290):998–1009.

    Article  CAS  Google Scholar 

  10. 10.

    Dietz KJ, Vogel MO, Viehhauser A. AP2/EREBP transcription factors are part of gene regulatory networks and integrate metabolic, hormonal and environmental signals in stress acclimation and retrograde signalling. Protoplasma. 2010;1–4(245):3–14.

    Article  CAS  Google Scholar 

  11. 11.

    Stockinger EJ, Gilmour SJ, Thomashow MF. Arabidopsis thaliana CBF1 encodes an AP2 domain-containing transcriptional activator that binds to the C-repeat/DRE, a cis-acting DNA regulatory element that stimulates transcription in response to low temperature and water deficit. Proc Natl Acad Sci USA. 1997;3(94):1035–40.

    Article  Google Scholar 

  12. 12.

    Lata C, Prasad M. Role of DREBs in regulation of abiotic stress responses in plants. J Exp Bot. 2011;14(62):4731–48.

    Article  CAS  Google Scholar 

  13. 13.

    Ito Y, Katsura K, Maruyama K, Taji T, Yamaguchi-Shinozaki K. Functional analysis of rice DREB1/CBF-type transcription factors involved in cold-responsive gene expression in transgenic rice. Plant Cell Physiol. 2006;1(47):141–53.

    Article  CAS  Google Scholar 

  14. 14.

    Liang CL, Li YN, Zhang XP, Song Y, Jia XD. Immunotoxicologic assessment of genetically modified drought-resistant wheat T349 with GmDREB1. Zhonghua Yu Fang Yi Xue Za Zhi. 2012;6(46):556–60.

    Google Scholar 

  15. 15.

    Hong JP, Kim WT. Isolation and functional characterization of the Ca-DREBLP1 gene encoding a dehydration-responsive element binding-factor-like protein 1 in hot pepper (Capsicum annuum L. cv. Pukang). Planta. 2005;6(220):875–88.

    Article  CAS  Google Scholar 

  16. 16.

    Jaglo-Ottosen KR, Gilmour SJ, Zarka DG, Schabenberger O, Thomashow MF. Arabidopsis CBF1 overexpression induces COR genes and enhances freezing tolerance. Science. 1998;5360(280):104–6.

    Article  Google Scholar 

  17. 17.

    Thomashow MF, Gilmour SJ, Stockinger EJ, Jaglo-Ottosen KR, Zarka DG. Role of the Arabidopsis CBF transcriptional activators in cold acclimation. Physiol Plant. 2001;112:171–5.

    CAS  Article  Google Scholar 

  18. 18.

    Haake V, Cook D, Riechmann JL, Pineda O, Thomashow MF, Zhang JZ. Transcription factor CBF4 is a regulator of drought adaptation in Arabidopsis. Plant Physiol. 2002;2(130):639–48.

    Article  CAS  Google Scholar 

  19. 19.

    Nakashima K, Ito Y, Yamaguchi-Shinozaki K. Transcriptional regulatory networks in response to abiotic stresses in Arabidopsis and grasses. Plant Physiol. 2009;1(149):88–95.

    Article  CAS  Google Scholar 

  20. 20.

    Shinwari ZK, Nakashima K, Miura S, Kasuga M, Seki M, Yamaguchi-Shinozaki K, Shinozaki K. An Arabidopsis gene family encoding DRE/CRT binding proteins involved in low-temperature-responsive gene expression. Biochem Biophys Res Commun. 1998;1(250):161–70.

    Article  Google Scholar 

  21. 21.

    Gilmour SJ, Sebolt AM, Salazar MP, Everard JD, Thomashow MF. Overexpression of the Arabidopsis CBF3 transcriptional activator mimics multiple biochemical changes associated with cold acclimation. Plant Physiol. 2000;124:1854–65.

    CAS  Article  Google Scholar 

  22. 22.

    Liu RM, Li FL, Hu GF, Hu BZ. Cold resistance of transgenic tobacco with CBF3 gene. J Northeast Agric Univ. 2011;1(42):119–23.

    Google Scholar 

  23. 23.

    Dubouzet J, Sakuma Y, Ito Y, Kasuga M, Dubouzet E, Miura S, Seki M, Shinozaki K, Yamaguchi-Shinozaki K. OsDREB genes in rice, Oryza sativa L., encode transcription activators that function in drought-, high-salt-, and cold-responsive gene expression. Plant J. 2003;33:751–63.

    CAS  Article  Google Scholar 

  24. 24.

    Agarwal P, Agarwal P, Reddy M, Sopory S. Role of DREB transcription factors in abiotic and biotic stress tolerance in plants. Plant Cell Rep. 2006b;25:1263–74.

    CAS  Article  Google Scholar 

  25. 25.

    Chen J, Dong Y, Wang Y, Liu Q, Zhang J, Chen S. An AP2/EREBP-type transcription-factor gene from rice is cold-inducible and encodes a nuclear-localized protein. Theor Appl Genet. 2003;107:972–9.

    CAS  Article  Google Scholar 

  26. 26.

    Mei Z, Wei L, Yu-Ping B, Zi-Zhang W. Isolation and identification of PNDREB1, a new DREB transcription factor from peanut (Arachis hypogaea L). Acta Agronom Sin. 2009;35:1973–80.

    Article  CAS  Google Scholar 

  27. 27.

    Liu N, Zhong N, Wang G, Li L, Liu X, He Y, Xia G. Cloning and functional characterization of PpDBF1 gene encoding a DRE-binding transcription factor from Physcomitrella patens. Planta. 2007;226:827–38.

    CAS  Article  Google Scholar 

  28. 28.

    Hadi F, Gilpin M, Fuller MP. Identification and expression analysis of CBF/DREB1 and COR15 genes in mutants of Brassica oleracea var. botrytis with enhanced proline production and frost resistance. Plant Physiol Biochem. 2011;49:1323–32.

    CAS  Article  Google Scholar 

  29. 29.

    Cong L, Chai TY, Zhang YX. Characterization of the novel gene BjDREB1B encoding a DRE-binding transcription factor from Brassica juncea L. Biochem Biophys Res Commun. 2008;4(371):702–6.

    Article  CAS  Google Scholar 

  30. 30.

    Al-Abed D, Madasamy P, Talla R, Goldman S, Rudrabhatla S. Genetic engineering of maize with the DREB1A/CBF3 gene using split-seed explants. Crop Sci. 2007;6(47):2390–402.

    Article  CAS  Google Scholar 

  31. 31.

    Tang M, Liu X, Deng H, Shen S. Over-expression of JcDREB, a putative AP2/EREBP domain-containing transcription factor gene in woody biodiesel plant Jatropha curcas, enhances salt and freezing tolerance in transgenic Arabidopsis thaliana. Plant Sci. 2011;181:623–31.

    CAS  Article  Google Scholar 

  32. 32.

    Chen S, Cui X, Chen Y, Gu C, Miao H, Gao H, Chen F, Liu Z, Guan Z, Fang W. CgDREBa transgenic chrysanthemum confers drought and salinity tolerance. Environ Exp Bot. 2011;74:255–60.

    CAS  Article  Google Scholar 

  33. 33.

    Yang W, Liu XD, Chi XJ, Wu CA, Li YZ, Song LL, Liu XM, Wang YF, Wang FW, Zhang C. Dwarf apple MbDREB1 enhances plant tolerance to low temperature, drought, and salt stress via both ABA-dependent and ABA-independent pathways. Planta. 2010;233:219–29.

    Article  CAS  Google Scholar 

  34. 34.

    Hsieh T, Lee J, Yang P, Chiu L, Charng Y, Wang Y, Chan M. Heterology expression of the Arabidopsis C-repeat/dehydration response element binding factor 1 gene confers elevated tolerance to chilling and oxidative stresses in transgenic tomato. Plant Physiol. 2002;129:1086–94.

    CAS  Article  Google Scholar 

  35. 35.

    Hsieh T, Lee J, Charng Y, Chan M. Tomato plants ectopically expressing Arabidopsis CBF1 show enhanced resistance to water deficit stress. Plant Physiol Biochem. 2002;130:618–26.

    CAS  Google Scholar 

  36. 36.

    Gil-Humanes J, Pistón F, Martín A, Barro F. Comparative genomic analysis and expression of theAPETALA2-like genes from barley, wheat, and barley-wheat amphiploids. BMC Plant Biol. 2009;1(9):66–78.

    Article  CAS  Google Scholar 

  37. 37.

    Nakano T, Suzuki K, Fujimura T, Shinshi H. Genome-wide analysis of the ERF gene family in Arabidopsis and rice. Plant Physiol. 2006;2(140):411–32.

    Article  CAS  Google Scholar 

  38. 38.

    Zhuang J, Chen JM, Yao QH, Xiong F, Sun CC, Zhou XR, Zhang J, Xiong AS. Discovery and expression profile analysis of AP2/ERF family genes from Triticum aestivum. Mol Biol Rep. 2011;2(38):745–53.

    Article  CAS  Google Scholar 

  39. 39.

    Zhang G, Ming C, Chen X, Xu Z, Shan G, Li LC, Li A, Guo J, Long M, Ma Y. Phylogeny, gene structures, and expression patterns of the ERF gene family in soybean (Glycine max L.). J Exp Bot. 2008;15(59):4095–107.

    Article  CAS  Google Scholar 

  40. 40.

    Zhuang J, Cai B, Peng RH, Zhu B, Jin XF, Xue Y, Gao F, Fu XY, Tian YS, Zhao W. Genome-wide analysis of the AP2/ERF gene family in Populus trichocarpa. Biochem Biophys Res Co. 2008;371:468–74.

    CAS  Article  Google Scholar 

  41. 41.

    Song X, Li Y, Hou X. Genome-wide analysis of the AP2/ERF transcription factor superfamily in Chinese cabbage (Brassica rapa ssp. pekinensis). BMC Genomics. 2013.

  42. 42.

    Zhang C, Zhang L, Wang D, Ma H, Liu B, Shi Z, Ma X, Chen Y, Chen Q. Evolutionary history of the glycoside hydrolase 3 (GH3) family based on the sequenced genomes of 48 plants and identification of jasmonic acid-related GH3 proteins in Solanum tuberosum. Int J Mol Sci. 2018;7(19):1850–65.

    Article  CAS  Google Scholar 

  43. 43.

    Lee SC, Lim MH, Yu JG, Park BS, Yang TJ. Genome-wide characterization of the CBF/DREB1 gene family in Brassica rapa. Plant Physiol Biochem. 2012;61:142–52.

    CAS  Article  Google Scholar 

  44. 44.

    Wang Z, Wang F, Tang J, Huang Z, Xiong A, Hou X. C-repeat binding factor gene family identified in non-heading Chinese cabbage is functional in abiotic and biotic stress response but different from that in Arabidopsis. Acta Physiol Plant. 2014;12(36):3217–29.

    Article  CAS  Google Scholar 

  45. 45.

    Li W, Dong J, Cao M, Gao X, Wang D, Liu B, Chen Q. Genome-wide identification and characterization of HD-ZIP genes in potato. Gene. 2019;697:103–17.

    CAS  Article  Google Scholar 

  46. 46.

    Mohseni S, Che H, Djillali Z, Dumont E, Nankeu J, Danyluk J, Gulick P. Wheat CBF gene family: identification of polymorphisms in the CBF coding sequence. Genome. 2012;12(55):865–81.

    Article  CAS  Google Scholar 

  47. 47.

    Ryu JY, Hong SY, Jo SH, Woo JC, Park CM. Molecular and functional characterization of cold-responsive C-repeat binding factors from Brachypodium distachyon. BMC Plant Biol. 2014;1(14):15.

    Article  Google Scholar 

  48. 48.

    Skinner JS, Zitzewitz JV, Szucs P, Marquez-Cedillo L, Filichkin T, Amundsen K, Stockinger EJ, Thomashow MF, Chen THH, Hayes PM. Structural, functional, and phylogenetic characterization of a LargeCBFGene family in Barley. Plant Mol Biol. 2005;4(59):533–51.

    Article  CAS  Google Scholar 

  49. 49.

    Thamilarasan SK, Park JI, Jung HJ, Nou IS. Genome-wide analysis of the distribution of AP2/ERF transcription factors reveals duplication and CBFs genes elucidate their potential function in Brassica oleracea. BMC Genomics. 2014;1(15):422.

    Article  CAS  Google Scholar 

  50. 50.

    Cao Y, Wang J, Guo L, Xiao K. Identification, characterization and expression analysis of transcription factor (CBF) genes in rice (Oryza sativa L.). Front Agric China. 2008;2:253–61.

    Article  Google Scholar 

  51. 51.

    Yamasaki K, Kigawa T, Seki M, Shinozaki K, Yokoyama S. DNA-binding domains of plant-specific transcription factors: structure, function, and evolution. Trends Plant Sci. 2013;5(18):267–76.

    Article  CAS  Google Scholar 

  52. 52.

    Zhang Y, Chen C, Jin X-F, Xiong A-S, Peng R-H, Hong Y-H, Yao Q-H, Chen J-M. Expression of a rice DREB1 gene, OsDREB1D, enhances cold and high-salt tolerance in transgenic Arabidopsis. BMB Rep. 2009;8(42):486–92.

    Article  Google Scholar 

  53. 53.

    Lynch M, Conery JS. The evolutionary fate and consequences of duplicate genes. Science. 2000;5494(290):1151–5.

    Article  Google Scholar 

  54. 54.

    Zhang Z, Li J, Zhao XQ, Wang J, Wong KS, Yu J. KaKs_Calculator: calculating Ka and Ks through model selection and model averaging. Genom Proteom Bioinform. 2006;4(4):259–63.

    CAS  Article  Google Scholar 

  55. 55.

    Li J, Zhang Z, Vang S, Yu J, Wong KS, Wang J. Correlation between Ka/Ks and Ks is related to substitution model and evolutionary lineage. J Mol Evol. 2009;4(68):414–23.

    Article  CAS  Google Scholar 

  56. 56.

    Sze H, Geisler M, Murphy AS. Linking the evolution of plant transporters to their functions. Front Plant Sci. 2014;4:547–8.

    Article  Google Scholar 

  57. 57.

    Dou H, Xv K, Meng Q, Li G, Yang X. Potato plants ectopically expressing Arabidopsis thaliana CBF3 exhibit enhanced tolerance to high-temperature stress. Plant Cell Environ. 2014;1(38):61–72.

    Google Scholar 

  58. 58.

    Kidokoro S, Watanabe K, Ohori T, Moriwaki T, Maruyama K, Mizoi J, Htwe NMPS, Fujita Y, Sekita S, Shinozaki K, Yamaguchi-Shinozaki K. Soybean DREB 1/ CBF-type transcription factors function in heat and drought as well as cold stress-responsive gene expression. Plant J. 2015;3(81):505–18.

    Article  CAS  Google Scholar 

  59. 59.

    Massa AN, Childs KL, Buell CR. Abiotic and biotic stress responses in group Phureja DM1–3 516 R44 as measured through whole transcriptome sequencing. Plant Genome. 2013;3(6).

  60. 60.

    Li J, Wang Y, Yu B, Song Q, Liu Y, Chen THH, Li G, Yang X. Ectopic expression of StCBF1 and ScCBF1 have different functions in response to freezing and drought stresses in Arabidopsis. Plant Sci. 2018;270:221–33.

    CAS  Article  Google Scholar 

  61. 61.

    Bouaziz D, Pirrello J, Ben Amor H, Hammami A, Charfeddine M, Dhieb A, Bouzayen M, Gargouri-Bouzid R. Ectopic expression of dehydration responsive element binding proteins (StDREB2) confers higher tolerance to salt stress in potato. Plant Physiol Bioch. 2012;60:98–108.

    CAS  Article  Google Scholar 

  62. 62.

    Charfeddine M, Bouaziz D, Charfeddine S, Hammami A, Ellouz ON, Bouzid RG. Overexpression of dehydration-responsive element-binding 1 protein (DREB1) in transgenic Solanum tuberosum enhances tolerance to biotic stress. Plant Biotechnol Rep. 2015;2(9):79–88.

    Article  Google Scholar 

  63. 63.

    Bouaziz D, Pirrello J, Charfeddine M, Hammami A, Jbir R, Dhieb A, Bouzayen M, Gargouri-Bouzid R. Overexpression ofStDREB1 transcription factor increases tolerance to salt in transgenic potato plants. Mol Biotechnol. 2013;3(54):803–17.

    Article  CAS  Google Scholar 

  64. 64.

    Charfeddine M, Charfeddine S, Bouaziz D, Messaoud RB, Bouzid RG. The effect of cadmium on transgenic potato (Solanum tuberosum) plants overexpressing the StDREB transcription factors. Plant Cell Tiss Org. 2016;3(128):521–41.

    Google Scholar 

  65. 65.

    Wang Y, Deng D, Bian Y, Lv Y, Xie Q. Genome-wide analysis of primary auxin-responsive Aux/IAA gene family in maize (Zea mays L.). Mol Biol Rep. 2010;8(37):3991–4001.

    Article  CAS  Google Scholar 

  66. 66.

    Liu B, Zhao S, Wu X, Wang X, Nan Y, Wang D, Chen Q. Characterization of phosphate transporter genes in potato. J Biotechnol. 2017;264:17–28.

    CAS  Article  Google Scholar 

  67. 67.

    Gu Z, Cavalcanti A, Chen FC, Bouman P, Li WH. Extent of gene duplication in the genomes of Drosophila, nematode, and yeast. Mol Biol Evol. 2002;3(19):256–62.

    Article  Google Scholar 

  68. 68.

    Jiang SY, González JM, Ramachandran S. Comparative Genomic and transcriptomic analysis of tandemly and segmentally duplicated genes in rice. PLoS ONE. 2013;5(8):e63551.

    Article  CAS  Google Scholar 

  69. 69.

    Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, Salzberg SL, Wold BJ, Pachter L. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010;5(28):511–5.

    Article  CAS  Google Scholar 

  70. 70.

    Nicot N, Hausman JF, Hoffmann L, Evers D. Housekeeping gene selection for real-time RT-PCR normalization in potato during biotic and abiotic stress. J Exp Bot. 2005;421(56):2907–14.

    Article  Google Scholar 

  71. 71.

    Liu F, Xu Y, Jiang H, Jiang C, Du Y, Gong C, Wang W, Zhu S, Han G, Cheng B. Systematic identification, evolution and expression analysis of the Zea mays PHT1 gene family reveals several new members involved in root colonization by Arbuscular Mycorrhizal Fungi. Int J Mol Sci. 2016;6(17):930–47.

    Article  CAS  Google Scholar 

Download references


We thank to Prof. Chen for suggesting the funding. We thank to Northwest A&F University for suggesting the laboratory.


This work was funded mainly by the National Key Research and Development Program of China (2018YFD0200805); the Project of Science and Technology from Shaanxi Province (No. 2017ZDXM-NY-004); and major collaborative innovation projects for production, education and research in Yangling demonstration zone (No. 2016CXY-05). In addition, this work was partially supported by State Key Laboratory of Crop Stress Biology in Arid Areas, China. The funding body played no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Author information




WL, DW and QC designed the experiments. WL, QC, DW, HL, MY and YC carried out the experiments. WL and MY analysed the experimental results. WL, MY, and DW analysed the data and developed the analysis tools. WL wrote the manuscript, and HL and QC revised the manuscript. All authors read and approved the manuscript.

Corresponding authors

Correspondence to Haibin Lu or Dongdong Wang or Qin Chen.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

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.

Select information about all CBF/DREB1 proteins. a: The protein length, chromosomal location, gene name, and group of each gene. b: Duplication events.

Additional file 2.

Evolutionary relationships (left), and motif predictions (right) of CBF/DREB1 genes in 43 plant species. The evolutionary history was inferred using the N-J method via MEGA 7. Bootstrap values of 1000 replications were used. The motifs, numbered 1–20, are displayed in different coloured boxes.

Additional file 3.

Logos of 20 motifs for all CBF/DREB1 proteins.

Additional file 4.

GO annotation information for each CBF/DREB1 member in potato.

Additional file 5.

Specific primers used for qRT–PCR analysis.

Additional file 6.

RNA-seq data of CBF/DREB1 proteins in potato from PGSC database.

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

Li, W., Chen, Y., Ye, M. et al. Evolutionary history of the C-repeat binding factor/dehydration-responsive element-binding 1 (CBF/DREB1) protein family in 43 plant species and characterization of CBF/DREB1 proteins in Solanum tuberosum. BMC Evol Biol 20, 142 (2020).

Download citation


  • CBF/DREB1 family
  • Heat resistance
  • Cold resistance
  • Potato (Solanum tuberosum)
  • Abiotic stress