Genome-wide analysis of the DREB family genes and functional identification of the involvement of BrDREB2B in abiotic stress in wucai (Brassica campestris L.)
BMC Genomics volume 23, Article number: 598 (2022)
Dehydration responsive element binding protein (DREB) is a significant transcription factor class known to be implicated in abiotic stresses. In this study, we systematically conducted a genome-wide identification and expression analysis of the DREB gene family, including gene structures, evolutionary relationships, chromosome distribution, conserved domains, and expression patterns. A total of 65 DREB family gene members were identified in Chinese cabbage (Brassica rapa L.) and were classified into five subgroups based on phylogenetic analysis. Through analysis of the conserved domains of BrDREB family genes, only one exon existed in the gene structure. Through the analysis of cis-acting elements, these genes were mainly involved in hormone regulation and adversity stress. In order to identify the function of BrDREB2B, overexpressed transgenic Arabidopsis was constructed. After different stress treatments, the germination rate, root growth, survival rate, and various plant physiological indicators were measured. The results showed that transgenic Arabidopsis thaliana plants overexpressing BrDREB2B exhibited enhanced tolerance to salt, heat and drought stresses. Taken together, our results are the first to report the BrDREB2B gene response to drought and heat stresses in Chinese cabbage and provide a basis for further studies to determine the function of BrDREBs in response to abiotic stresses.
Abiotic stresses, such as high and low temperatures, drought, and salinization, have severe adverse effects on plant growth and development [1, 2]. At the transcription level, transcription factors (TFs) can specifically bind to cis-acting elements in the promoter region, thereby regulating the expression of plant stress response genes. Therefore, the key role of transcription factors in plants under environmental stress has attracted increasing research attention [3, 4].
The AP2/ERF (APETALA2/ethylene-responsive factor) superfamily, which is one of the largest groups of TFs families, is involved in the regulation of plant developmental processes and biotic and abiotic stress responses [5, 6]. The genes of the superfamily are characterized by encoding at least one APETALA2 (AP2) domain and can be further separated into ERF, AP2, DREB, RAV and Soloist families. The ERF family contains a single AP2 domain, while the AP2 gene family contains one or two AP2 domains [7, 8]. The RAV (Related to ABI3 and VP1) family encodes proteins encoding a single AP2 domain plus an additional B3 domain . There is another subfamily with a single AP2 domain called soloist .
DREB transcription factors are a subgroup of the AP2/ERF (APETALA2/ethylene-responsive factor) gene family, which contain a single conserved AP2 domain that specifically interacts with the dehydration responsive element (DRE) (core motif A/GCCGAC) [10,11,12]. DREB transcription factors bind to the DRE and the C-repeat responsive (CRT) cis-element in the promoter region of stress resistance genes to regulate the expression of a series of downstream genes and enhance the resistance of plants to different abiotic stresses . Recently, a number of DREB homologous genes have been identified from different plants, such as Arabidopsis [13, 14], rice , wheat , and soybean .
Based on the sequence identities of the DNA-binding domains, DREB proteins can be divided into six subgroups (A1-A6) . DREB members in different subgroups play different roles in plants. In subgroup A2, DREB2A and DREB2B, as the main functional genes in DREB2s, have been reported to play an important role in resistance to dehydration, high salt, and heat shock . In Arabidopsis, DREB2A and DREB2B were identified and actively responded to dehydration and high salinity stress . In mung bean (Vigna radiata L.), VrDREB2A was induced by drought and salt stresses, and heterologous expression of VrDREB2A enhanced the tolerance to drought and high salt stress without growth retardation in Arabidopsis . In mustard (Brassica juncea L), BjDREB2 playing a pivotal role in ABA-independent gene expression under drought stress .
Wucai (Brassica campestris L. ssp. chinensis var. rosularis), a variant of non-heading Chinese cabbage, is a semi-hardy vegetable widely grown in the Yangtze-Huai River basin [23, 24]. Recently, abiotic stresses such as high temperature, freezing, and drought have greatly affected the yield and quality of wucai. Therefore, improvement of resistance to abiotic stress in wucai is desirable. There are few studies on the molecular mechanism of stress resistance, the expression regulation of genes related to stress resistance, and the stress signal transduction mechanisms. In this study, the genome-wide identification and mining of DREB transcription factor family genes in Chinese cabbage were systematically performed. We also determined the gene structure, evolution relationships, cis-acting element evolution relationship, and chromosome distribution of the Chinese cabbage DREB gene family in detail. In addition, we cloned and characterized the BrDREB2B gene in wucai. Furthermore, the overexpression of BrDREB2B enhanced the drought, high-salt, and heat tolerance of transgenic Arabidopsis plants. Our study provides a theoretical foundation to understand the role of BrDREBs in resistance to abiotic stress and a basis for further study of the biological functions of DREB genes.
Materials and methods
Genome-wide identification of BrDREB genes in chinese cabbage
The Hidden Markov Model (HMM) profile of the AP2 domain (PF00847) was downloaded from the protein family database (http://pfam.xfam.org/) with an E-value < 0.001, and was used to search the Chinese cabbage genome database (http://brassicadb.org/brad/) [25,26,27,28]. There is only one AP2 domain in the DREB subfamily member protein, and the 14th and 19th positions are valine and glutamic acid, respectively. The NCBI-CDD database (https://www.ncbi.nlm.nih.gov/Structure/bwrpsb/bwrpsb.cgi) and SMART website (http://smart.embl-heidelberg.de/) were used for domain prediction and identified members of the DREB family of Chinese cabbage. The physical and chemical properties of the DREB family gene were analysed using the Expasy website (https://web.expasy.org/compute_pi/) , and the subcellular locations of DREB family members were predicted using the WoLF PSORT website (https://wolfpsort.hgc.jp) .
Conserved motifs and cluster analysis of the BrDREB gene family
We predicted the conserved motifs of the BrDREB protein using the MEME website (http://meme-suite.org/tools/meme). The identified Chinese cabbage DREB amino acid sequence and the reported Arabidopsis DREB amino acid sequence were aligned using MEGA7.0 software , and unrooted phylogenetic trees were constructed by the neighbor-joining (NJ) method. We used the bootstrap method in the phylogeny test, and the number of bootstrap replications was set to 1000.
Analysis of cis-acting elements of the BrDREB gene family
The sequence upstream (2000 bp) of the transcriptional start site was used to identify cis-elements in the promoter sequences of the DREB family genes in B. rapa. The cis-acting elements were predicted and analyzed using the PlantCARE website (http://bioinformatics.psb.ugent.be/webtools/plantcare/htmL/)  and were visualized with TBtools software .
Chromosome distribution and collinearity analysis
The chromosome distribution of the BrDREB genes were extracted from the B. rapa genome annotation gff3 file for collinearity analysis by TBtools software. Analyses of interspecies synteny between BrDREBs and AtDREBs were investigated using MCScan X .
Plant materials and stress treatments
Wucai (Brassica campestris L.) cultivar WS-1, a high-generation inbred line, was selected from Anhui Agricultural University. Seedlings were germinated and cultured in a growth chamber at 24 ± 1 °C (day) and 16 ± 1 °C (night), with a relative humidity of 75–80% and a light intensity of 300 µmol·m− 2·s− 1 under a 14-h/10-h (light/dark) photoperiod. After 20 days, unfolded functional leaves were collected for RNA extraction by total RNA kit (Takara Biomedical Technology Co., Beijing, China). Primer software v6.0 (Premier Biosoft International, Palo Alto, CA, USA) was used to design specific gene primers (Table S3). The gene encoding actin was used as the control. RT-qPCR was performed using SYBR GREEN Master Mix (Vazyme Biotechnology Co., Ltd., Nanjing, China). Relative gene expression levels were calculated using the 2−ΔΔCT method .
Arabidopsis (Columbia-0) plants were placed in a growth chamber at 25 °C (light for 18 h) and 15 °C (dark for 6 h), with a light intensity of 300 µmol·m− 2·s− 1. After 30 days, leaves of wild type and transgenic Arabidopsis were collected for analysis at the gene expression level. Transgenic T3 generation and wild type seedings were treated with 1/2 Murashige and Skoog (MS) liquid medium containing 150 mM NaCl  and 250 mmol·L− 1 mannito , respectively. Seedings in 1/2 MS liquid medium were subjected to 40 °C to mimic heat stress . Seedlings under normal conditions were used as the control.
Generation of BrDREB2B‑overexpressing Arabidopsis
BrDREB2B cDNA was cloned into the pMD19-T simple vector. BrDREB2B was inserted into plasmid pCambia1305-35 S-nFLAG-cMYC, and transgenic Arabidopsis plants were cultivated by the inflorescence infection method. The regenerated plants were placed in a medium containing 50% hygromycin B and grown in an incubator. Four weeks later, the copy number of BrDREB2B of each transgenic line was detected by Quantitative Real-time PCR (qRT-PCR), and three overexpression positive transgenic lines were identified (OE#1, OE#7, and OE#8).
Physiological index measurement of transgenic plants
The transgenic T3 generation seeds and wild-type seeds were planted on ½ MS medium at 40 °C, ½ MS medium with 150 mmol·L− 1 NaCl, and ½ MS medium with 250 mmol·L− 1 mannitol. Normal ½ MS medium was used as the control group, and the germination rate was observed after one week. The 5-day-old transgenic and wild-type seedlings were carefully transferred to ½ MS medium containing 150 mmol·L− 1 NaCl. Normal ½ MS medium was used as the control group. After 2 weeks, 20 transgenic and wild-type plants were randomly selected for root length measurement. The 3-week-old transgenic and wild-type plants were placed in an artificial climate incubator at 40 °C/30°C (day for 18 h/night for 6 h) for 3 d, compared with 25 °C/15°C (day for 18 h/night for 6 h) as a control. The growth of transgenic and wild-type plants was observed. Transgenic and wild-type plants were deprived of water for 14 d and then re-watered for 3 d. We observed the growth of the plants and determined the survival rate.
Relative electrical conductivity (REC) was measured according to the method described by Baziramakenga et al. with some modifications . The H2O2 and O2− content, as well as the total antioxidant capacity, were measured using the Solarbio reagent kit (Cat #BC3595, 1290 and 1315, Beijing Solarbio Science & Technology Co., Ltd., Beijing, China).
Identification and characterization of the BrDREBs
A total of 65 BrDREB genes were identified from the Brassica rapa genome (Table S1), which are unevenly distributed on 10 chromosomes. The ORFs of the BrDREB genes ranged from 471 to 1158 bp in length, encoding polypeptides 156 aa to 385 aa in length. The predicted molecular mass was between 17.46 and 43.05 kDa, and the isoelectric point also varied widely from 4.48 to 9.73, suggesting that different DREB proteins might operate in different microenvironments. Predictions from the WoLF PSORT website showed that seven genes were subcellularly located on the chloroplast, and the other 58 genes were all located on the nucleus.
Gene structure and conserved motif analysis of the BrDREB gene family
To further study the structural characteristics of the BrDREB gene family, we analyzed the gene structure and conserved motifs according to the full-length phylogenetic relationships (Fig. 1). Through cluster analysis, it was found that the closer the genetic relationship of genes, the more similar the number and position of motif distribution, suggesting functional similarity within subgroups. A total of ten conserved motifs were found in the DREB family of Brassica rapa. Motif 1 was found in all BrDREB proteins, while Motif 9 was found only in BrDREB2G, BrDREB2A3, and BrDREB2F2. In addition, structure analysis revealed that all BrDREB genes had no introns and only contained one exon, from which it can be inferred that the gene structure is highly conserved.
In order to understand the pedigree and functional characteristics of the BrDREB gene family, a multiple sequence alignment of the DREB family genes in Chinese cabbage and Arabidopsis was carried out. Thus, a phylogenetic tree was constructed (Fig. 2), and BrDREB family genes were divided into 6 subgroups through cluster analysis (Groups 1–6). The BrDREB gene was unevenly distributed among 5 subgroups. Among them, Group 4 had the most BrDREB genes, Group 3 had only two BrDREB genes, and Group 6 had no BrDREB genes. Through cluster analysis, it was found that the DREB genes of Chinese cabbage and Arabidopsis have high homology. It was speculated that the orthologous genes may have the same function.
Analysis of cis-acting elements of the BrDREB gene family
To further investigate the function of BrDREB genes, we predicted the cis-acting element of the putative promoter region of BrDREB genes using the PlantCARE database. The cis-acting element analysis was carried out on the promoters of the DREB gene family in Brassica rapa (Fig. 3). As shown in Fig. 3, we identified nine cis-acting elements according to their functional annotations, of which five elements were related to hormone response and three elements were related to adversity stress. Some diverse distribution patterns of cis-acting elements were observed in the promoter region of BrDREB genes, indicating that BrDREB was particular to various different biological processes. It was notable that all BrDREB genes contained the cis-acting related to hormone regulation, such as gibberellin, methyl jasmonate (MeJA), abscisic acid, auxin, and salicylic acid responsiveness elements. These results suggested that the expression of DREB genes in Brassica rapa is regulated by cis-elements associated with hormone signal transduction and abiotic stress tolerance.
Chromosome distribution and collinearity analysis of the BrDREB gene family
We analyzed the distribution of the DREB family genes on chromosomes, and the results showed that BrDREB the gene family members were unevenly distributed on 10 chromosomes (Fig. 4). Among them, chromosome A03 was the most populated, with 12 genes, while the A06 chromosome was the least populated, with only 2 genes. In addition, some DREB genes had tandem duplications, resulting in the presence of homologous genes in BrDREB, which enlarged the DREB genome of Brassica rapa.
Collinearity analysis of DREB genes between Brassica rapa and Arabidopsis found that a large number of members in the DREB gene family had a collinear relationship: 39 Arabidopsis DREB genes and 65 DREB genes (Fig. 5). There were 152 gene pairs, indicating high homology between AtDREBs and BrDREBs (Tab. S2). The ratio of nonsynonymous substitution rate (Ka) to synonymous substitution rate (Ks) of all the 152 gene pairs were below 1. These genes underwent strong purification and selection in evolution and have relatively consistent functions and effects.
Genetic transformation of Arabidopsis
To further study the function of DREB genes, we analyzed the expression level of DREB genes under heat and cold stress (Fig. 6a). We found that the expression level of the BrDREB2B gene increased under heat stress, so we generated transgenic Arabidopsis plants overexpressing BrDREB2B under the control of the CaMV35S promoter (Fig. 6b). Quantitative RT-PCR results showed that the transcript level of BrDREB2B was highly induced in selected transgenic plants, confirming that the plants were successfully transformed with pCambia1305-35 S-BrDREB2B-nFLAG-cMYC (Fig. 6c).
Effects of different stresses on the germination of transgenic plants
The seeds of the T3 BrDREB2B overexpression lines (OE#1, OE#7 and OE#8) and wild-type were subjected to high temperature (40 °C), salt (NaCl, 150 mM), and drought (mannitol, 250 mM) stress treatments. The germination status was observed (Fig. 7a). In the germination stage, seed germination of the wild-type and transgenic plants did not differ under normal conditions (Fig. 7). However, in the 150 mM mannitol treatment, the seed of transgenic lines displayed a significantly higher germination rate than that of the wild type. In the 150 mM NaCl and 40 °C heat stress treatments, seed germination of the wild type was seriously inhibited. In the presence of 150 mM NaCl, 50% of the wild-type seeds germinated, whereas the germination rates of the OE#1, OE#7 and OE#8 transgenic lines were 90, 70, and 60%. At 40 °C, we found 30% of the wild-type seeds germinated, whereas the germination rates of the OE#1, OE#7, and OE#8 transgenic lines were 90, 80, and 40%, respectively.
Effects of salt stress on the root growth of transgenic plants
Transgenic and WT Arabidopsis seeds were grown on 1/2 MS medium for 5 days and then transferred to 1/2 MS medium containing 150 mM NaCl. Under the normal condition, the root length of transgenic seedlings was not significantly different from that of the wild type (Fig. 8). However, under 150 mM NaCl treatment for 14 d, the average root length of transgenic lines was 4.6 cm, which was 2.42 times that of the wild type plants. The results showed the enhanced salt tolerance of the BrDREB2B transgenic lines.
Effects of different stresses on the growth of transgenic plants
As shown in Fig. 9a, under heat stress, transgenic lines overexpressing BrDREB2B grew well with green leaves, while the leaves of the wild type grew weakly. Under drought stress for 7 d, some leaves of the wild-type plants were dry and curled, while leaves of the transgenic Arabidopsis overexpressing BrDREB2B were still green (Fig. 9b). Two weeks after watering was stopped, the wild-type plants showed a state of dryness and wilting, and a large number of them even died, whereas the OE#1 transgenic plants showed mild wilting, and the other two lines showed a state of leaf curling and a lack of water. After 3 days of recovery, most of the transgenic plants returned to normal growth, while the wild-type plants were weak and wilted with curled leaves. The survival rate of the OE#1 transgenic plants exceeded 90% after 3 days of rewatering, while the survival rate of wild-type plants was less than 10%, which was significantly lower than that of transgenic plants. Over-expression of BrDREB2B significantly enhanced the survival rate of transgenic Arabidopsis under heat and drought stresses.
The effects of different stresses on the physiological indices of transgenic plants
In addition, the physiological characteristics including the accumulation of H2O2 and O2·−, REC, and total antioxidant capacity also reflected the stress tolerance. We detected these indices in the wild type and three lines of transgenic Arabidopsis over-expressing BrDREB2B (OE#1, OE#7, and OE#8). Under normal conditions, the leaves of WT and transgenic Arabidopsis showed no difference. Compared with the wild-type, the REC, H2O2 content and O2·− production rate of the transgenic plants under heat and salt stresses decreased, while the total antioxidant capacity increased (Fig. 10). This indicated that compared with wild-type plants, transgenic plants had greater reactive oxygen species scavenging ability, thereby maintaining cell homeostasis. These results further indicated that overexpression of BrDREB2B improved the tolerance of transgenic Arabidopsis to salt and heat stresses.
According to the number of AP2 domains and the existence of other DNA binding domains, AP2/ERF superfamily can be divided into AP2, ERF, DREB, RAV and Soloist subfamilies . Although the sequence of AP2/ERF domain is highly conserved, each subfamily shows different DNA binding characteristics . Generally, the ERF subfamily binds to an AGCCGCC sequence . AP2 subfamily binds to the GCAC(A/G)N(A/T) TCCC(A/G)ANG(C/T) element [42, 43] and is regulated by microRNA (mir172) [44, 45]. The RAV family binds to the CAACA and CACCTG sequences . DREB subfamily usually interacts with CCGA/CC core sequence .
AP2/ERF TFs have been successfully identified in many plants, such as Arabidopsis, rice [8, 48], grape , poplar (Populus tomentosa) , corn , wheat (Triticum wheat) , soybean (Glycine max)  and Chinese cabbage (Brassica rapa ssp. pekinensis) . AP2/ERF TFs can regulate various biological processes in plant function and development [53, 54], such as drought (SHN1, SHN2 and SHN3), salt (AP37, CaRAV1), freezing (TaERF1) [55,56,57], cell dedifferentiation (WIND1) , ABA induced crown root sprouting (CRL5) , plant height (NSAP2) and leaf shape [60,61,62].
As an important branch of AP2/ERF, DREB subfamily is used as a feasible candidate gene to improve abiotic stress tolerance of crops . The DREB-type transcription factors have recently been recognized in various plants, for example, Arabidopsis , rice , tomato , tobacco , barley , sorghum , and maize . The DREB genes also play a key role in plant responses to multiple abiotic stresses, including osmotic (CkDREB) , drought (OsDREB1) [70, 71], cold (AtCBF1) , heat (ZmDREB2A, AtDREB1A) [73, 74], and high-salt stress (CaDREBLP1) . However, the DREB genes have not been systematically studied in Chinese cabbage.
In this study, a total of 65 DREB family members were identified in Chinese cabbage (Table S1). Combined with the previous identification of 56 DREB genes in Arabidopsis, 56 in rice, and 38 in grape, the BrDREB family members are more abundant. It could be due to the triploidization of the whole genome of cabbage. The duplication of chromosome fragments or the entire genome is the main source of evolution, including the generation of new gene functions and expression patterns . In this study, some DREB genes had tandem duplications, resulting in the presence of homologous genes in BrDREB, which enlarged the DREB genome of Brassica rapa (Fig. 4). Furthermore, the genome size of each species also has an impact on the number of DREB family members. A large number of members in the DREB gene family had a collinear relationship between AtDREBs and BrDREBs (Fig. 5, Tab. S2). And the ratio of Ka to Ks of all the 152 gene pairs were below 1. These genes had high homology and underwent strong purification and selection in evolution and have relatively consistent functions and effects.
To investigate the evolutionary relatedness of the identified sequences together with DREB genes encoded by the other fully sequenced B. rapa and Arabidopsis, we performed phylogenetic reconstruction using the conserved DREB transcription factor domain. Based on phylogenetic tree analysis from (A) thaliana and (B) rapa, BrDREBs are classified into 6 subgroups, which is coherent with prior statements [18, 77]. However, the BrDREB gene was not present in group 6, so it can be assumed that a loss event occurred in the genetic evolution of B. rapa. The gene structure analysis revealed that BrDREBs contain no introns and only one exon, which is consistent with previous results for Arabidopsis . To characterize the distribution of genes, we found that these genes are unevenly distributed on 10 chromosomes, and some genes exhibit tandem duplication. Some gene duplication events of CBF/DREB family TFs in B. oleracea were found . This indicates that BrDREB could be generated through chromosome duplication and gene duplication events.
Furthermore, we analyzed the cis-acting elements on the promoters of the BrDREB gene. These elements are mainly related to hormone regulation (such as response to abscisic acid, salicylic acid, and auxin) and abiotic stress (such as response to drought and low-temperature stress). The cis-acting elements analysis predicts that the function of the BrDREBs could play a role in hormonal response and abiotic stress responses, which is similar to the results of the study of MnDREBs in mulberry . Mizoi et al. found DREB transcription factors can regulate the expression of stress resistance genes through ABA-dependent and ABA-independent pathways . In tomato, the ectopic expression of SlDREB3 can increase the growth of plant roots and improve the photosynthetic rate by reducing the level of ABA . Taken together, the BrDREBs may be involved in the regulation of hormone metabolism pathways to improve plant resistance to environmental stresses such as drought, low temperature, and salinization.
In order to further explore the role of BrDREB2B in abiotic stress, we constructed a pCambia1305-35 S-BrDREB2B-nFLAG-cMYC overexpression vector and obtained transgenic Arabidopsis plants. The results indicate that transgenic plants overexpressing BrDREB2B showed better growth status and higher survival rates after drought, high-salt and high-temperature stresses than wild-type plants. Similarly, overexpression of OsDREB2A , TaDREB3 , and GmDREB2A  in Arabidopsis enhanced tolerance to abiotic stress. SsDREB2D and SsDREB2F play a key role in sugarcane dehydration reaction, which helps sugarcane plants recover from drought stress and low temperature stress . OsDREB2A and OsDREB2B were reported to involve in the embryo and endosperm development in rice . Therefore, BrDREB2B may play a key role in response and resistance to abiotic stress in wucai and can be used as a candidate gene for molecular breeding to increase the yield of wucai.
The accumulation of reactive oxygen species (ROS) is closely related to the environmental stress of plants, leading to oxidative damage and lipid peroxidation of plant cell membranes . In this study, the relative electrolytic leakage content, H2O2 content, and superoxide anion production rate of the three overexpression transgenic lines after different stress treatments decreased. At the same time, the total antioxidant capacity increased compared with the wild-type. These results verify that BrDREB2B improves its tolerance to abiotic stress by maintaining the homeostasis of intracellular reactive oxygen species, which is consistent with a report for soybean . In carrot, DcDREB1A overexpression in transgenic plants resulted in increased superoxide dismutase (SOD) and peroxidase (POD) activities under drought stress .
In this study, 65 BrDREBs were identified in Brassica rapa and the gene structure, protein conserved domains, evolutionary relationship, and cis-acting elements were analyzed. There is high homology in the DREB gene of Chinese cabbage and Arabidopsis. In addition, many abiotic stress response cis-elements were found in the promoter regions of BrDREB genes. Finally, the heterologous expression of BrDREB2B improved the salt, heat, and drought stress tolerance of transgenic Arabidopsis. This work has provided a solid foundation for the investigation of the function and the mechanism of response to abiotic stress of the DREB TF family in wucai.
Availability of data and materials
The raw RNA-Seq data used in this study have been deposited in the Nation Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA) database under the accession number PRJNA694542 (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA694542).
Boyer JS. Plant productivity and environment, Science (New York, N.Y.). 1982; 218:443–448.
Mittler R. Abiotic stress, the field environment and stress combination. Trends Plant Sci. 2006;11:15–9.
Goyal RK, Tulpan D, Chomistek N, Fundora DG-P, West C, Ellis BE, Frick M, Laroche A, Foroud NA. Analysis of MAPK and MAPKK gene families in wheat and related Triticeae species. BMC Genomics. 2018;19(1):178–203.
Wan Y, Mao M, Wan D, Yang Q, Yang F, Mandlaa, Li G, Wang R. Identification of the WRKY gene family and functional analysis of two genes in Caragana intermedia. BMC Plant Biol. 2018;18:31–46.
Wessler SR. Homing into the origin of the AP2 DNA binding domain. Trends Plant Sci. 2005;10:54–6.
Mizoi J, Shinozaki K, Yamaguchi-Shinozaki K. AP2/ERF family transcription factors in plant abiotic stress responses. Biochim Biophys Acta. 2011;1819(2):86–96.
Licausi F, Giorgi FM, Zenoni S, Osti F, Pezzotti M, Perata P. Genomic and transcriptomic analysis of the AP2/ERF superfamily in Vitis vinifera. BMC Genomics. 2010;11:719–34.
Nakano T, Suzuki K, Fujimura T, Shinshi H. Genome-wide analysis of the ERF gene family in Arabidopsis and rice. Plant Physiol. 2006;140(2):411–32.
Hu L, Liu S. AP2/ERF family transcription factors in plant abiotic stress responses. Genet Mol Biol. 2011;34(4):624–33.
Ali N, Hadi F. CBF/DREB transcription factor genes play role in cadmium tolerance and phytoaccumulation in Ricinus communis under molybdenum treatments. Chemosphere. 2018;208: 425–432.
Eckardt NA. DREB Duo Defines Distinct Drought and Cold Response Pathways. Plant cell. 2019;31:1196–7.
Niu X, Luo TL, Zhao HY, Su YL, Ji WQ, Li HF. Identification of wheat DREB genes and functional characterization of TaDREB3 in response to abiotic stresses. Gene. 2020;740: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;94:1035–40.
Liu Q, Kasuga M, Sakuma Y, Abe H, Miura S, Yamaguchi-Shinozaki K, Shinozaki K. Two transcription factors, DREB1 and DREB2, with an EREBP/AP2 DNA binding domain separate two cellular signal transduction pathways in drought- and low-temperature-responsive gene expression, respectively, in Arabidopsis. Plant Cell. 1998;10:1391–1406.
Dubouzet JG, Sakuma Y, Ito Y, Kasuga M, Dubouzet EG, 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.
Kume S, Kobayashi F, Ishibashi M, Ohno R, Nakamura C, Takumi S. Differential and coordinated expression of Cbf and Cor/Lea genes during long-term cold acclimation in two wheat cultivars showing distinct levels of freezing tolerance. Genes Genet Syst. 2005;80:185–97.
Li XP, Tian AG, Luo GZ, Gong ZZ, Zhang JS, Chen SY. Soybean DRE-binding transcription factors that are responsive to abiotic stresses. Theoretical Appl Genet. 2005;110:1355–62.
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;290:998–1009.
Agarwal PK, Agarwal P, Reddy MK, Sopory SK. Role of DREB transcription factors in abiotic and biotic stress tolerance in plants. Plant Cell Reports. 2006;25:1263–1274.
Nakashima K, Shinwari ZK, Sakuma Y, Seki M, Miura S, Shinozaki K, Yamaguchi-Shinozaki K. Organization and expression of two Arabidopsis DREB2 genes encoding DRE-binding proteins involved in dehydration- and high-salinity-responsive gene expression. Plant Mol Biol. 2000;42:657–65.
Chen H, Liu L, Wang L, Wang S, Cheng X. VrDREB2A, a DREB-binding transcription factor from Vigna radiata, increased drought and high-salt tolerance in transgenic Arabidopsis thaliana. J Plant Res. 2016;129:263–73.
Sb A, Sp B, Jkt C, et al. Antioxidant, physiological and biochemical responses of drought susceptible and drought tolerant mustard (Brassica juncea L) genotypes to rhizobacterial inoculation under water deficit stress. Plant Physiol Biochem. 2019;143:19–28.
Yuan L, Xie S, Nie L, Zheng Y, Wang J, Huang J, Zhao M, Zhu S, Hou J, Chen G, Wang C. Comparative Proteomics Reveals Cold Acclimation Machinery Through Enhanced Carbohydrate and Amino Acid Metabolism in Wucai (Brassica Campestris L.). Plants-Basel. 2019;8:474–474.
Nie L, Zheng Y, Zhang L, Wu Y, Zhu S, Hou J, Chen G, Tang X, Wang C, Yuan L. Characterization and transcriptomic analysis of a novel yellow-green leaf wucai (Brassica campestris L.) germplasm. BMC Genomics. 2021;22:258–258.
Zeng RF, Zhou JJ, Liu SR, et al. Genome-Wide Identification and Characterization of SQUAMOSA—Promoter-Binding Protein (SBP) Genes Involved in the Flowering Development of Citrus Clementina. Biomolecules. 2019;9(2):66–87.
El-Gebali S, Mistry J, Bateman A, Eddy SR, Luciani A, Potter SC, Qureshi M, Richardson LJ, Salazar GA, Smart A, Sonnhammer ELL, Hirsh L, Paladin L, Piovesan D, Tosatto SCE, Finn RD. The Pfam protein families database in 2019. Nucleic Acids Research. 2019;47:D427-D432.
Cheng F, Liu SY, Wu J, Fang L, Sun SL, Liu B, Li PX, Hua W, Wang XW. BRAD, the genetics and genomics database for Brassica plants. Bmc Plant Biology. 2011;11:136–136.
Cheng X, Li G, Manzoor MA, et al. In Silico Genome-Wide Analysis of Respiratory Burst Oxidase Homolog (RBOH) Family Genes in Five Fruit-Producing Trees, and Potential Functional Analysis on Lignification of Stone Cells in Chinese White Pear. Cells. 2019;8(6):520–40.
Gasteiger E, Gattiker A, Hoogland C, Ivanyi I, Appel RD, Bairoch A. ExPASy: the proteomics server for in-depth protein knowledge and analysis. Nucleic Acids Research. 2003;31:3784–3788.
Horton P, Park KJ, Obayashi T, Fujita N, Harada H, Adams-Collier CJ, Nakai K. WoLF PSORT: protein localization predictor. Nucleic Acids Research. 2007;35:W585-W587.
Kumar S, Stecher G, Tamura K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol Biol Evol. 2016;33:1870–4.
Lescot M, Dehais P, Thijs G, Marchal K, Moreau Y, Van de Peer Y, Rouze P, Rombauts S. Plant CARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences. Nucleic Acids Research. 2002;30:325–327.
Chen CJ, Chen H, Zhang Y, Thomas HR, Frank MH, He YH, Xia R. TBtools: An Integrative Toolkit Developed for Interactive Analyses of Big Biological Data. Molecular Plant. 2020;13:1194–1202.
Wang YP, Tang HB, DeBarry JD, Tan X, Li JP, Wang XY, Lee TH, Jin HZ, Marler B, Guo H, Kissinger JC, Paterson AH. MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Research. 2012;40:14.
Schmittgen TD, Livak KJ. Analyzing real-time PCR data by the comparative C(T) method. Nat Protoc. 2008; 3(6):1101–1108.
Liu JX, Hao-Bin HU, Wang X, et al. Effect of Nitric Oxide on Proline Accumulation in Ryegrass Seedlings Subjected to Salt Stress. Acta Agrestia Sinica. 2010;18(6):786–91.
Fan X, Zhao L L, Zhai H H, et al. Functional Characterization of AtNEK6 Overexpression in Cotton Under Drought and Salt Stress. Scientia Agricultura Sinica, 2018.
Yuan Lingyun Zheng Yushan, Nie, Libing, et al. Transcriptional profiling reveals changes in gene regulation and signaling transduction pathways during temperature stress in wucai (Brassica campestris L.). BMC Genomics. 2021;22(1):687–687.
Baziramakenga R, Leroux GD, Simard RRJJoCE. Effects of benzoic and cinnamic acids on membrane permeability of soybean roots. J Chem Ecolog. 1995;21:1271–85.
Thamilarasan Senthil Kumar, Park Jong-In. Hee-Jeong Jung and Ill-Sup Nou;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;15:422–34.
Ohme-Takagi M, Shinshi H. Ethylene-inducible DNA binding proteins that interact with an ethylene-responsive element. Plant Cell. 1995;7(2):173–82.
Nole-Wilson S, Krizek BA. DNA binding properties of the Arabidopsis floral development protein AINTEGUMENTA. Nucleic Acids Res. 2000;28(21):4076–82.
Gong W, He K, Covington M, Dinesh-Kumar SP, Snyder M, Harmer SL, Zhu YX, Deng XW. The development of protein microarrays and their applications in DNA-protein and protein-protein interaction analyses of Arabidopsis transcription factors. Mol Plant. 2008;1(1):27–41.
Shigyo M, Hasebe M, Ito M: Molecular evolution of the AP2 subfamily. Gene. 2006;366(2):256–265.
Kim S, Soltis PS, Wall K, Soltis DE: Phylogeny and domain evolution in the APETALA2-like gene family. Mol Biol Evol 2006, 23(1):107–120.
Sohn KH, Lee SC, Jung HW, Hong JK, Hwang BK. Expression and functional roles of the pepper pathogen-induced transcription factor RAV1 in bacterial disease resistance, and drought and salt stress tolerance. Plant Mol Biol. 2006;61(6):897–915.
Jiang C, Iu B, Singh J. Requirement of a CCGAC cis-acting element for cold induction of the BN115 gene from winter Brassica napus. Plant Mol Biol. 1996;30(3):679–84.
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;245(1–4):3–14
Zhuang J, Cai B, Peng RH, Zhu B, Jin XF, Xue Y, Gao F, Fu XY, Tian YS, Zhao W, et al. Genome-wide analysis of the AP2/ERF gene family in Populus trichocarpa. Biochem Biophys Res Commun. 2008;371(3):468–74.
Hu L, Liu S. Genome-wide identification and phylogenetic analysis of the ERF gene family in cucumbers. Genet Mol Biol. 2011;34(4):624–33.
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;38(2):745–53.
Zhang G, Chen M, Chen X, Xu Z, Guan S, Li LC, Li A, Guo J, Mao L, Ma Y. Phylogeny, gene structures, and expression patterns of the ERF gene family in soybean (Glycine max L.). J Exp Bot. 2008;59(15):4095–107.
Xu ZS, Chen M, Li LC, Ma YZ. Functions and application of the AP2/ERF transcription factor family in crop improvement. J Integr Plant Biol. 2011;53(7):570–85.
Gutterson N, Reuber TL. Regulation of disease resistance pathways by AP2/ERF transcription factors. Curr Opin Plant Biol. 2004;7(4):465–71.
Aharoni A, Dixit S, Jetter R, Thoenes E, van Arkel G, Pereira A. The SHINE clade of AP2 domain transcription factors activates wax biosynthesis, alters cuticle properties, and confers drought tolerance when overexpressed in Arabidopsis. Plant Cell. 2004;16(9):2463–80.
Oh SJ, Kim YS, Kwon CW, Park HK, Jeong JS, Kim JK. Overexpression of the transcription factor AP37 in rice improves grain yield under drought conditions. Plant Physiol. 2009;150(3):1368–79.
Xu ZS, Xia LQ, Chen M, Cheng XG, Zhang RY, Li LC, Zhao YX, Lu Y, Ni ZY, Liu L, et al. Isolation and molecular characterization of the Triticum aestivum L. ethylene-responsive factor 1 (TaERF1) that increases multiple stress tolerance. Plant Mol Biol. 2007;65(6):719–32.
Iwase A, Mitsuda N, Koyama T, Hiratsu K, Kojima M, Arai T, Inoue Y, Seki M, Sakakibara H, Sugimoto K, et al. The AP2/ERF transcription factor WIND1 controls cell dedifferentiation in Arabidopsis. Curr Biol. 2011;21(6):508–14.
Kitomi Y, Ito H, Hobo T, Aya K, Kitano H, Inukai Y. The auxin responsive AP2/ERF transcription factor CROWN ROOTLESS5 is involved in crown root initiation in rice through the induction of OsRR1, a type-A response regulator of cytokinin signaling. Plant J. 2011;67(3):472–84.
Yan X, Zhang L, Chen B, Xiong Z, Chen C, Wang L, Yu J, Lu C, Wei W. Functional identification and characterization of the Brassica napus transcription factor gene BnAP2, the ortholog of Arabidopsis thaliana APETALA2. PLoS One. 2012;7(3):e33890.
Luo H, Chen S, Jiang J, Teng N, Chen Y, Chen F. The AP2-like gene NsAP2 from water lily is involved in floral organogenesis and plant height. J Plant Physiol. 2012;169(10):992–8.
Jiang F, Guo M, Yang F, Duncan K, Jackson D, Rafalski A, Wang S, Li B. Mutations in an AP2 transcription factor-like gene affect internode length and leaf shape in maize. PLoS One. 2012;7(5):e37040.
Okamuro JK, Caster B, Villarroel R, Montagu MV, Jofuku KD. The AP2 domain of APETALA2 defines a large new family of DNA binding proteins in Arabidopsis. Plant Biol. 1997;94:7076–81.
Zhou J, Tang X, Martin GBJEJ. The Pto kinase conferring resistance to tomato bacterial speck disease interacts with proteins that bind a cis-element of pathogenesis-related genes. EMBO J. 2014;16:3207–18.
Shinshi OTJTPCO. Ethylene-inducible DNA binding proteins that interact with an ethylene-responsive element. Plant Cell. 1995;7:173–82.
Guo B, Wei Y, Xu R, Shen L, Xu RJPO. Genome-Wide Analysis of APETALA2/Ethylene-Responsive Factor (AP2/ERF) Gene Family in Barley (Hordeum vulgare L.). Plos One. 2016;11:e0161322.
Akbudak MA, Filiz E, Kontbay KJB. DREB2 (dehydration-responsive element-binding protein 2) type transcription factor in sorghum (Sorghum bicolor): genome-wide identification, characterization and expression profiles under cadmium and salt stresses. 3 Biotech. 2018;8(10):426–41.
Qin F, Kakimoto M, Sakuma Y, Maruyama K, Osakabe Y, Tran LP, Shinozaki K, Yamaguchi-Shinozaki K. Regulation and functional analysis of ZmDREB2A in response to drought and heat stresses in Zea mays L. Plant Cell Physiol. 2006; 2006,47(-):S225.
Wang X, Chen X, Liu Y, Gao H, Wang Z, Sun G. CkDREB gene in Caragana korshinskii is involved in the regulation of stress response to multiple abiotic stresses as an AP2/EREBP transcription factor. Mol Biol Rep. 2011;38(4):2801–11.
Ito Y, Katsura K, Maruyama K, Taji T, Kobayashi M, Seki M, Shinozaki K, 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;47(1):141–53.
Liang CL, Li YN, Zhang XP, Song Y, Wang W, Fang J, Cui WM, Jia XD. Immunotoxicologic assessment of genetically modified drought-resistant wheat T349 with GmDREB1. Zhonghua Yu Fang Yi Xue Za Zhi. 2012;46(6):556–60.
Jaglo-Ottosen KR, Gilmour SJ, Zarka DG, Schabenberger O, Thomashow MF. Arabidopsis CBF1 overexpression induces COR genes and enhances freezing tolerance. Science. 1998;280(5360):104–6.
Qin F, Kakimoto M, Sakuma Y, Maruyama K, Osakabe Y, Tran LS, Shinozaki K, Yamaguchi-Shinozaki K. Regulation and functional analysis of ZmDREB2A in response to drought and heat stresses in Zea mays L. Plant J. 2007;50(1):54–69.
Hong B, Ma C, Yang Y, Wang T, Yamaguchi-Shinozaki K, Gao J. Overexpression of AtDREB1A in chrysanthemum enhances tolerance to heat stress. Plant Mol Biol. 2009;70(3):231–40.
Hong JP, Kim WT. Isolation and functional characterization of the Ca-DREBLP1 gene encoding a dehydration-responsive element bindingfactor-like protein 1 in hot pepper (Capsicum annuum L. cv Pukang). Planta. 2005;220(6):875–88.
Vakirlis N, Sarilar V, Drillon G, Fleiss A, Agier N, Meyniel JP, Blanpain L, Carbone A, Devillers H, Dubois KJGR. Reconstruction of ancestral chromosome architecture and gene repertoire reveals principles of genome evolution in a model yeast genus. Genome research. 2016;26(7):918–32.
Akhtar M, Jaiswal A, Taj G, Jaiswal JP, Qureshi MI, Singh NK. DREB1/CBF transcription factors: their structure, function and role in abiotic stress tolerance in plants. J Genet. 2012;91:385–95.
Liu XQ, Zhu JJ, Wei CJ, Guo Q, Bian CK, Xiang ZH, Zhao AC. Genome-wide identification and characterization of the DREB transcription factor gene family in mulberry. Biologia Plantarum. 2015;59:253–265.
Mizoi J, Shinozaki K, Yamaguchi-Shinozaki K. AP2/ERF family transcription factors in plant abiotic stress responses. Biochim. Biophys. Acta-Gene Regul. Mech. 2012;1819:86–96.
Upadhyay RK, Gupta A, Soni D, Garg R, Pathre UV, Nath P, Sane AP. Ectopic expression of a tomato DREB gene affects several ABA processes and influences plant growth and root architecture in an age-dependent manner. J Plant Physiol. 2017;214:97–107.
Mallikarjuna G, Mallikarjuna K, Reddy MK, Kaul T. Expression of OsDREB2A transcription factor confers enhanced dehydration and salt stress tolerance in rice (Oryza sativa L.). Biotechnol. Lett. 2011;33:1689–1697.
Mizoi J, Ohori T, Moriwaki T, Kidokoro S, Todaka D, Maruyama K, Kusakabe K, Osakabe Y, Shinozaki K, Yamaguchi-Shinozaki K. GmDREB2A;2, a Canonical Dehydration-Responsive element-binding protein2-Type Transcription Factor in Soybean, Is Posttranslationally Regulated and Mediates Dehydration-Responsive Element-Dependent Gene Expression. Plant Physiol. 2013;161:346–61.
Li, Z., Wang, G., Liu, X. et al. Genome-wide identification and expression profiling of DREB genes in Saccharum spontaneum. BMC Genomics. 2021;22:456–470.
Herath V. Small family, big impact: in silico analysis of DREB2 transcriptionfactor family in rice. Comput Biol Chem. 2016;65:128–139.
Nejadsadeghi L, Maali-Amiri R, Zeinali H, Ramezanpour S, Sadeghzade B. Membrane fatty acid compositions and cold-induced responses in tetraploid and hexaploid wheats. Mol Biol Rep. 2015;42:363–72.
Zhou YB, Chen M, Guo JK, Wang YX, Min DH, Jiang QY, Ji HT, Huang W, Wei CY, Xu HJ, Chen X, Li LC, Xu ZS, Cheng XG, Wang CX, Wang CS, Ma YZ. Overexpression of soybean DREB1 enhances drought stress tolerance of transgenic wheat in the field. J Experiment Botany. 2020;71:1842–57.
Li T, Huang Y, Khadr A, Wang YH, Xu ZS, Xiong AS. DcDREB1A, a DREB-binding transcription factor from Daucus carota, enhances drought tolerance in transgenic Arabidopsis thaliana and modulates lignin levels by regulating lignin-biosynthesis-related genes. Environ. Exp. Bot. 2020;169:12.
We thank LetPub (www.letpub.com) for its linguistic assistance during the preparation of this manuscript.
This research study was supported by Performance Award of Anhui Pakchoi Germplasm Resource Bank (Nursery) Project, University Natural Science Research Project of Anhui Province (KJ2020ZD11), the Graduate Innovation Fund of Anhui Agricultural University (2021yjs-6), and Collaborative Innovation Project of Universities in Anhui Province (GXXT-2021-054). We thank LetPub (www.letpub.com) for its linguistic assistance during the preparation of this manuscript.
Ethics approval and consent to participate
Line WS-1 was screened by the corresponding author, WS-1 was identified by the corresponding author, and all seed materials were kept in the Vegetable Breeding Laboratory of the College of Horticulture, Anhui Agricultural University. No specific permits were required for the described field studies. The location is not privately-owned or protected in any way, and the field studies did not involve endangered or protected species. We complied with the IUCN Policy Statement on Research, Involving Species at Risk of Extinction and the Convention on the Trade in Endangered Species of Wild Fauna and Flora.
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.
About this article
Cite this article
Wu, Y., Zhang, L., Nie, L. et al. Genome-wide analysis of the DREB family genes and functional identification of the involvement of BrDREB2B in abiotic stress in wucai (Brassica campestris L.). BMC Genomics 23, 598 (2022). https://doi.org/10.1186/s12864-022-08812-1
- Brassica campestris
- DREB gene family
- Genome-wide identification
- Functional verification