- Research article
- Open Access
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 Evolutionary Biology volume 20, Article number: 142 (2020)
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 .
Many genes related to plant stress resistance, such as those related to membrane stability , antifreeze protein genes , osmo-regulatory genes , and antioxidant enzyme activity genes , 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 . 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) . The DREB subfamily consists of six subgroups, of which the A1 subgroup includes CBF/DREB1 transcription factors . The CBF/DREB1 gene is a type of plant-specific transcription factor that was first discovered in Arabidopsis thaliana (Arabidopsis) . 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 .
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 . In Arabidopsis, six CBF/DREB1 genes have been characterized . 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 . In Arachis hypogaea, PNDREB1 is strongly upregulated under low-temperature conditions . The expression of PpCBF1 is also induced by NaCl, cold, drought, and ABA treatments in Physcomitrella patens . 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 . In transgenic Arabidopsis, salt tolerance and freezing tolerance are enhanced by overexpressing JcDREB from the woody biodiesel plant Jatropha curcas . Overexpression of CkDREBa from Caragana korshinski confers drought and salinity tolerance to transgenic chrysanthemum . MbDREB1 from dwarf apple (Malus baccata) increases plant tolerance to low temperature, drought, and salt stress via both ABA-dependent and ABA-independent pathways . 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  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  (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).
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.
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 . 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.
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  (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).
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) . 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).
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.
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)  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 . 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 . 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 .
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 . The substitution rate, known as the Ka/Ks (non-synonymous/synonymous) ratio, is an indicator for estimating selective pressure , and is commonly used to analyse the molecular evolutionary rate, and to investigate the evolutionary direction and selective strength of a coding sequence . 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) . 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 .
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 . 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 . 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 . 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 . 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 . 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 . 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 . 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 . 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) . 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 . 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 . The StDREB2 gene of the A2 subgroup was confirmed to be involved in the plant response to several treatments (salt, ABA, etc.) . 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) . 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 (https://phytozome.jgi.doe.gov/pz/portal.html), PGSC (https://solanaceae.plantbiology.msu.edu/pgsc_download.shtml), and the NCBI (https://www.ncbi.nlm.nih.gov) database. Six Arabidopsis CBF/DREB1 protein sequences that were obtained from TAIR database (https://www.Arabidopsis.org/) 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 . All candidate CBF/DREB1 protein sequences were screened by using the Conserved Domains Database (CDD) (https://www.ncbi.nlm.nih.gov/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; https://meme-suite.org/tools/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 . The EvolView online tool (https://www.evolgenius.info/ 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. . 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 . The segmental duplication pairs were analysed by using the Plant Genome Duplication Database (PGDD; https://chibba.agtec.uga.edu/duplication). 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 .
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 . 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 .
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
Ethylene response factor
Related to ABI3/VP1
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.
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.
Graham LA, Liou YC, Walker VK, Davies PL. Hyperactive antifreeze protein from beetles. Nature. 1997;6644(388):727–8.
Ashraf M, Foolad MR. Roles of glycine betaine and proline in improving plant abiotic stress resistance. Environ Exp Bot. 2007;2(59):206–16.
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.
Warren GJ. Cold stress: Manipulating freezing tolerance in plants. Curr Biol. 1998;15(8):514–6.
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.
Wessler SR. Homing into the origin of the AP2 DNA binding domain. Trends Plant Sci. 2005;2(10):54–6.
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.
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.
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.
Lata C, Prasad M. Role of DREBs in regulation of abiotic stress responses in plants. J Exp Bot. 2011;14(62):4731–48.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Lynch M, Conery JS. The evolutionary fate and consequences of duplicate genes. Science. 2000;5494(290):1151–5.
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.
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.
Sze H, Geisler M, Murphy AS. Linking the evolution of plant transporters to their functions. Front Plant Sci. 2014;4:547–8.
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.
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.
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).
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Select information about all CBF/DREB1 proteins. a: The protein length, chromosomal location, gene name, and group of each gene. b: Duplication events.
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.
Logos of 20 motifs for all CBF/DREB1 proteins.
GO annotation information for each CBF/DREB1 member in potato.
Specific primers used for qRT–PCR analysis.
RNA-seq data of CBF/DREB1 proteins in potato from PGSC database.
About this article
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). https://doi.org/10.1186/s12862-020-01710-8
- CBF/DREB1 family
- Heat resistance
- Cold resistance
- Potato (Solanum tuberosum)
- Abiotic stress