Genome-Wide Identification and Analysis of the Cytochrome B5 Protein Family in Chinese Cabbage (Brassica rapa L. ssp. Pekinensis)

Cytochrome B5 (CB5) family proteins play an important role in various oxidation/reduction reactions in cells as the electron donor and are involved in a variety of biotic and abiotic stress processes. However, the function of the CB5s in Brassica rapa is still unclear. In this study, we carried out genome-wide identification, characterization, and expression analysis of BrCB5s in different tissues under adversities and stresses. It was identified that fifteen BrCB5s were distributed on different chromosomes, which were classified into seven groups (A-G) according to its phylogenetic relationship. Phylogenetic analysis of the CB5 protein sequences from six species showed that the BrCB5s conduct a close evolutionary process with the CB5s of Arabidopsis thaliana and far from those of Oryza sativa. Protein interaction analysis showed that 40 interaction patterns were predicted including two Sucrose Transporter 4 subfamily proteins (SUT 4) and Fatty Acid Hydroxylase 2 protein (FAH 2) can interact with most members of BrCB5s. The expression profile analysis indicated that BrCB5s were differentially expressed in different tissues, and the transcript abundances were significantly different under various abiotic stresses and plant hormone treatments. Our study provides a basis for a better understanding of the characteristics and biological functions of the CB5 family genes in Chinese cabbage during plant development, especially in plant responses to multiple stresses.


Introduction
Cytochrome P450 (P450) belongs to a family of hemebinding proteins, which catalyzes multiple functional monooxygenase reactions involved in oxidative metabolism. In plants, P450s play a role in the generation of secondary metabolites [1,2], some of which are synthesized to organize and integrate vital biological processes, and the others are accumulated as defense responders to biotic or abiotic stresses. P450 proteins are important to plants in processes from biosynthesis and metabolism to growth regulation.
Cyt b5 proteins (CB5s), which enhance the turnover of related catabolic enzymes, are important family members of P450 [3]. CB5, anchored to the endoplasmic reticulum (ER), was firstly observed in the larvae of the silkworm Platysumia cecropia by Sanborn and Williams in 1950 [4]. CB5s are small (~15kD) heme-binging proteins ubiquitously expressed in animals, plants, fungi, and purple photosynthetic bacteria [5] and function as electron transporters.
Due to their roles in cell detoxification and drug metabolism, CB5s have been studied extensively in animal [6][7][8][9]. However, the functions of these proteins have yet to be understood in plants. With the development of molecular biology and sequencing technology, the whole genomes of many species have been known. Multiple CB5 isoforms have been discovered in higher plants. For example, seven CB5s have been found in the model plant Arabidopsis thaliana and seventeen in the plant Oryza sativa; in contrast, only a signal copy of CB5 has been discovered in mammals [10]. A hypothesis has been proposed that a large number of CB5 isoforms are needed in response to the increasing number of P450 proteins, as CB5s enhance the activities of P450 proteins by supplying electrons or by physical interaction independent of electron donation [11,12]. A classical experiment has been designed to understand the relationship between CB5s and P450s via observing the discoloration of flower petals in a CB5 knockdown mutant of petunia [13]. The result showed that the change of the flower color is accompanied by the decrease in the activity of a biosynthetic P450 enzyme 3‚5 ′ -hydroxylase. Recently, the glucosinolate (GLS) levels of two T-DNA insertion mutants of CB5 isoform C (CB5C) from Arabidopsis thaliana were characterized. The first one (cb5c-1) was a knockdown mutant with an insertion in the coding region, while the other one (cb5c-2) was a "knockabout" mutant with an insertion in the 3 ′ untranslated region of the gene [14]. GLS relates to plant defense [15,16], and P450 enzymes carry out central catalytic steps in the GLS biosynthetic pathway [17]. These two mutants lead to a subtle and distinct decrease in the levels of GLS under methyl jasmonate treatments. These findings suggest CB5 is required by P450 to be fully functional.
In this study, 12 putative BrCB5s with 3 BrCB5 NADPHdependent reductases of Chinese cabbage were selected from the Brassica database (http://brassicadb.org/brad/) [18]. We performed a genome-wide bioinformatics analysis of the BrCB5s, including genome location, gene structure, and evolutionary divergence. We identified the expression patterns of these BrCB5s by quantitative real-time PCR (qRT-PCR) in different tissues and in response to various treatments. In addition, we conducted further experiments on the functional characterization of the BrCB5s in Chinese cabbage.
Chromosomes 02, 04, and 05 have two BrCB5 genes and the chromosomes 03 and 09 have three, and the last three chromosomes have only one gene.

Gene Structures and Conserved Motifs of BrCB5s.
Intron/exon regions of the BrCB5s were identified by aligning the CDSs to the genomic sequence. The results showed that all the BrCB5 gene sequences contained introns except BrCB5a. Although the number of introns varied from zero to eight, most of the genes (10 out of 15) contained two introns. Two genes contained eight introns, and three genes contained seven, one, and zero introns, respectively ( Figure 1(b)). Additionally, the number of introns in the genes of the same subfamily of BrCB5s was not always the same. For instance, BrCB5a and BrCB5m, which belong to the same subfamily F, contained zero and one introns, respectively ( Figure 1).
Simple sequence repeat (SSR) markers are extensively used in plant genetic mapping and molecular breeding due to genetic codominance abundance, wide distribution in genomes, multiallelic variation, high reproducibility, and high level of polymorphisms [19]. In this study, 17 SSR markers, including eleven di-, five tri-, and one tetranucleotide motifs, were detected in the 15 BrCB5s using the online SSR identification tool SSRIT (Table 3). BrCB5b, BrCB5e, BrCB5j, BrCB5k, and BrCB5n had one SSR marker. BrCB5g, BrCB5l, and BrCB5o had two SSR markers. BrCB5d and BrCB5m had three SSR markers, while BrCB5a, BrCB5c, BrCB5f, BrCB5h, and BrCB5i had no SSR markers. Among these SSR markers, ten were found in introns and seven were found in exons.
To better understand the function of BrCB5s, we used the MEME web server (http://meme.nbcr.net/meme/cgi-bin/ meme.cgi) to analyze the domain distribution in BrCB5s (Supplementary file 2. Figure S2). Motif 1, specified as the N-terminal hydrophilic haem-binding domain, was found in 12 of the 15 BrCB5s. The other three CB5s without this motif were found to belong to the NADH-dependent reductase subfamily (Figure 1(a)). Motif 3, specified as the C-terminal hydrophobic region, was predicted to be present in 8 of the 15 BrCB5s, which is less conserved compared with the haem-binding domain, as it is related to the transmembrane function [20]. By MEME analysis, a motif, namely motif 2, was found in 10 of 15 BrCB5s, including BrCB5c, BrCB5 e-l, and BrCB5n. Multiple sequence alignment was also conducted on the BrCB5 proteins (Supplementary file 3. Figure S3). The result confirmed that all the BrCB5 proteins had a "HPGG" haem-binding domain except these three NADH-dependent reductases (BrCB5b, BrCB5d, and BrCB5o). The BrCB5 family was found to contain a nonconserved C-terminal binding domain. All the results were similar to that by the MEME analysis (Supplementary file 2. Figure S2).

Phylogenetic Analysis and Duplication Events of BrCB5
Genes. Classifying genes and phylogenetic analysis are imperative for identifying the functions of a gene family. As Chinese cabbage is one of the most important leafy head vegetables, we analyze its phylogenetic relationship with several model crops: Raphanus sativus, Oryza sativa, Solanum lycopersicum, and Glycine max. CB5s involved in these species were obtained using BLAST online tool on NCBI (https://www.ncbi.nlm.nih.gov/). MEGA5 with the bootstrap neighbor-joining method was used to construct the phylogenetic tree (Figure 2), including 15 predicted BrCB5s, 7 Arabidopsis CB5s (AtCB5), 13 Raphanus sativus CB5s (RsCB5s), 17 Oryza sativa CB5s (OsCB5), 12 Solanum lycopersicum CB5s (SlCB5s), and 22 Glycine max CB5s (GlyCB5s) (including their NADPH-dependent reductases). The result showed that, as a member of the Cruciferae family, CB5s of Chinese cabbage is closely related to Raphanus sativus and Arabidopsis, which clustered together with Chinese cabbage. But for Oryza sativa, lots of its protein members branched alone, indicating the furthest relationship ( Figure 2). BrCB5s   Figure S2). The syntenic relationship of CB5s between Arabidopsis thaliana and Brassica rapa is shown in Table 1 with the data downloaded from the Brassica database (http:// brassicadb.org/brad/); at the same time, as the great value to the decryption of the evolutionary mechanism, we also analyzed the genome duplication events in the evolutionary process of Chinese cabbage. The result indicated that 15 BrCB5s were derived from 7 blocks of 4 tPCK (translocation Proto-Calepineae Karyotype) chromosomes. Coincidently, they were evenly distributed on three subgenomes (LF, MF1, and MF2). There were one to three copies of BrCB5s syntenically corresponding to an AtCB5. For example, BrCB5j, BrCB5g, and BrCB5f corresponded to AtCB5d in the J block; BrCB5o, d, and b corresponded to AtCB5g in the R block, while the other BrCB5s were either duplicated or singletons.

Protein
Interaction Analysis of BrCB5s. The knowledge of the functional interactions of proteins is indispensable to a widely understanding of the molecular machinery. The STRING database (http://string-db.org/) is known as an online tool on which we can get predicted protein-protein association information related to our target proteins. In this study, 40 predicted interactive candidates of BrCB5 proteins were identified (Table 4).
Among these candidates, two encoded by Br031692 and Bra019972 which belonged to the Sucrose Transporter (SUT 4) subfamily interact with 10 of the total 15 BrCB5s. It was reported that the Arabidopsis SUT4 (AtSUT4) was a functional interactive protein to 5 members of the AtCB5 family which include a total of 7 members [21]. A previous study also showed that an Apple Sucrose Transporter, MdSUT1, interacted physically with an apple cytochrome b5 (MdCYB5) in vitro and vivo, which was verified by the yeast two-hybrid, immunocoprecipitation, and bimolecular fluorescence complementation assays [22].
Additionally, there are 2 predicted proteins that can interact with 9 and 7 members of BrCB5s, respectively. One is Bra006538, a NADH-cytochrome b5 reductase, and another one is a FAH 2 protein (Fatty Acid Hydroxylase 2) encoded by Bra013479. It was reported that NADH-Cyt b5 reductase has the ability to specifically reduce the content of CB5 in the presence of NADH [23]. It was also reported that AtFAH1 and AtFAH2 physically interact with AtCb5s from where they obtain electrons [24].
Moreover, 10, 6, and 20 predicted functional proteins interact with 3, 2, and 1 members of BrCB5s, respectively, and BrCB5e and m have no predicted functional proteins with a score ≥ 0:6 which was given by the STRING online tool.

Expression Patterns of BrCB5s in Various Tissues, in
Response to Abiotic Stresses and Hormone Treatments. To explore the functions of the BrCB5s in Chinese cabbage growth and development, mRNA expression analysis was performed in six tissues, including root (R), stem (S), young leaf (YL), old leaf (OL), flower bud (FB), and immature silique (IS). According to the results, all the 15 BrCB5s had a higher expression in the old leaves than in the young leaves; all the BrCB5s had a higher expression level in the immature silique than that in the stem; all the BrCB5s had a higher expression in the roots than in the stems, except BrCB5a and BrCB5d, which had a similar expression level in both tissues ( Figure 3). In addition, eight of fifteen BrCB5s (BrCB5a, BrCB5b, BrCB5d, BrCB5f, BrCB5g, BrCB5h, and BrCB5j) were expressed mainly in immature siliques. There were five BrCB5s (BrCB5c, BrCB5i, BrCB5k, BrCB5l, and BrCB5n), and two BrCB5s (BrCB5e and BrCB5m) were expressed mainly in old leaves and in roots, respectively. Some paralogs showed a similar expression pattern in different tissues such as BrCB5g/j and BrCB5d/b, while some exhibited different expression tendencies in different tissues, for example, BrCB5c/e, BrCB5h/i, and BrCB5k/m. The varied expression patterns of the BrCB5s indicated that they might play different roles in the growth and development of those tissues.
The responses of the BrCB5s to abiotic stresses, such as the common osmotic reagent (polyethylene glycol, PEG 6000 ), salt (NaCl), heat (35°C), and cold (4°C) stress, were investigated to understand the potential roles of those BrCB5 genes. The expression of these BrCB5s displayed different abundances under different stresses (Figures 4 and 5 Figure 2: Phylogenetic analysis of CB5 proteins in Arabidopsis thaliana, Brassica rapa, Raphanus sativus, Oryza sativa, Solanum lycopersicum, and Glycine max, whose members were denoted by hollow circles, green circles, blue inverted triangles, cyan triangles, red diamonds, and yellow squares, respectively. The tree, based on the core CB5 domains in the six species, was constructed using the neighbor-joining (NJ) method (bootstrap 1000 replicates) by MEGA5 software. Table  4: Predicted functional partners of BrCB5s using the STRING Database (https://string-db.org/), which is known of predicting protein-protein interactions. The scores and description of the predicted functional partners were also downloaded from the database. Here, we selected predicted interacting proteins with a score more than 0.6.
Predicted functional partners BrCB5 proteins (score ≥ 0:6) Description of predicted functional partners      It is known that plant hormone plays crucial roles in plant growth and defense signaling. In order to explore the expression pattern of BrCB5s under plant hormone stress, Gibberellin A3 (GA3), Abscisic Acid (ABA), and Salicylic Acid (SA) were used to treat Chinese cabbage plants ( Figure 6).
For the SA treatment, the expression levels of 13 BrCB5s were higher than that of the CK group at 3 h and 24 Figure 4: Expression analysis of the BrCB5 genes under salt stresses. Three-week-old plants were treated with 20% (w/v) PEG 6000 and 200 mM NaCl for 0, 3, and 24 h before the mature leaves were harvested. CK plants were treated with the same quantity water compared with salt. Expression of the BrCB5 genes was normalized to those of BrACT1 and shown relative to the expression of CK at 0 h. The 2 -ΔΔCt method was used to calculate the expression levels of target genes in different tissues. * indicated that the expression level is significantly different from the value of the control ( * p < 0:05, * * p < 0:01).
Under the ABA treatment, compared with the CK group, the expression levels of most BrCB5s were upregulated, while the expression of BrCB5f and BrCB5i was almost the same at 3 h with that of the CK group. BrCB5m showed a specific expression pattern with a basically unchanged expression level compared with that in the CK group at both 3 h and 24 h. There were four genes (BrCB5b, BrCB5d, BrCB5j, and BrCB5m) which were downgraduated by about 3-fold, 9-fold, 100-fold, and 20-fold from 0 h to 24 h, respectively, while the expression level of BrCB5f increased by 6  were downregulated at 3 h, but upregulated at 24 h; BrCB5d maintained the same expression level with that in the CK group at 3 h, and then the expression decreased at 24 h. In the experimental group, under the GA3 treatment, the expression level of two genes (BrCB5f and BrCB5o) increased continuously and that of three genes (BrCB5d, BrCB5e, and BrCB5l) declined. The others had a significantly decreased expression level at 3 h, and then the expression increased at 24 h. For instance, the expression level of BrCB5i was declined by more than 20-fold at 3 h and then increased by almost 22-fold at 24 h. Generally, the expression levels of CB5s were higher compared with that of the CK group at the same time point after treatment, respectively, but a conspicuously declined trend was found in the experimental group. For instance, BrCB5i and BrCB5j showed their lowest expression levels 3 hours after the four phytohormone treatments, while BrCB5e and BrCB5l had their lowest levels 24 hours after the treatments.
BrCB5m was not as unique as it used to be, for BrCB5b, BrCB5d, and BrCB5j showed a similar downward trend under all these stresses.

Discussion
The CB5 genes may play an important role in Chinese cabbage. However, there is little information about the function of the BrCB5s. It is unknown how the BrCB5s regulate growth and development and how they respond to a variety of biotic and abiotic stresses. To better and accurately explore the function of BrCB5s, we carried out bioinformatics analyses according to our previous studies [25,26], which have been successfully applied to genome-wide analysis of GRF and VQ family genes in Chinese cabbage. Additionally, expression patterns of BrCB5s under some conventional abiotic and hormone treatments were carried out to explore the stress tolerance mechanisms of Chinese cabbage. Genome duplication, which results in nonfunctionalization, subfunctionalization, and neofunctionalization, plays a vital role in expanding genome content and diversifying gene function [27][28][29]. Although the whole genome triplication (WGT) event was undergone by Brassica rapa [30] after its divergence from Arabidopsis thaliana, the B. rapa genome is approximately 4-fold larger than the Arabidopsis genome, and its gene number is twice than that of Arabidopsis [31,32], suggesting that a large number of genes were lost during the genome duplication. In our study, we found 15 BrCB5s had a syntenic relationship with 7 AtCB5s (Table 1). Moreover, the WGT event had greatly expanded the gene family members of B. rapa. The phylogenetic and gene duplication analyses showed that the BrCB5s without tandem duplication contained two triplets, four duplicates, and one singleton. These results are similar with the A. thaliana genome, in which a large proportion of gene families are divided into low tandem and high segmental duplication class [33]. In our study, there were also some paralogs showing different expression patterns, for example, BrCB5h/BrCB5i and BrCB5c/BrCB5e in different tissues and BrCB5a/BrCB5m and BrCB5b/BrCB5o under various abiotic and hormone stresses, revealing that BrCB5 paralogs might be maintained by subfunctionalization. Our results also confirmed the previously reported finding that duplicate genes can develop divergent patterns of gene expression for stably maintaining by subfunctionalization [34].
The expression patterns of BrCB5s in different tissues indicated that they might play important roles in the growth and development of specific tissues or organs. For instance, BrCB5a, BrCB5b, BrCB5d, BrCB5h, and BrCB5o had the highest expression level in immature siliques, suggesting that these genes are involved mainly in the growth and development of immature siliques in Chinese cabbage. Thus, we might improve Chinese cabbage seed harvest through regulating the expression of these genes.
Previous studies reported that CB5s interact with P450 enzymes either by supplying electrons or by modulating their activity via physical interaction independent of electron donation [11,13]. It is known that P450s carry out essential enzymatic steps in the glucosinolate (GLS) biosynthetic pathway. A. thaliana produces GLS as their major class of specialized metabolites involved in plant defense [16,35]. Additionally, it was confirmed that a CB5 protein, namely CB5C, can improve the efficiency of this pathway [15]. The GLS levels have been investigated using 2 mutants of the CB5C gene, and the results show that both mutants lead to subtle but distinct alterations in the levels of individual GLS. We explored the expression levels of some BrGLS genes under the NaCl treatment (Supplementary file 4. Figure S4 and Supplementary file 6. Table S2). The results showed that there existed some correlations between the two gene families under NaCl stress. Based on these results, GLS was used as an indicator to explore the role of CB5s in plants under stresses.
Salinity is considered the major abiotic stress affecting plant physiology and development [36,37]. In our study, we found that, under the salt treatment, most of BrCB5s were upregulated at both time points we selected ( Figure 4). One study showed that the GLS content is increased after 5 days under NaCl stress in B. rapa [38], which is consistent with our results. However, little is known about the expression of CB5 under PEG 6000 stress. We found that most BrCB5s were greatly upregulated at 3 h and then downregulated at 24 h upon the PEG 6000 treatment. Water stress increases the glucosinolate accumulation in Brassica species-Nasturtium officinale L. [39], Brassica oleracea L. var. capitata [40], Brassica oleracea L. var. italica [41,42], Brassica napus L. [43], Brassica rapa ssp. rapifera L. [44], and Brassica carinata L. [45]. It suggests that the response of BrCB5s to PEG 6000 was an instantaneous and rapid process.
It is reported that elevated temperatures (21-34°C) can increase the glucosinolate levels in Brassica rapa and lowmedium temperatures (15-27°C) can decrease the glucosinolate levels [46]. However, in our study, we found that the expression of almost all BrCB5s was downregulated at 35°C, while the expression reached the highest at 3 h under the 4°C treatment ( Figure 5). The discrepancy has been observed by Justen et al. [47], and they refer these differences to different growth conditions and distinct genotypes.
It is known that ABA and SA play critical roles in mediating plant defense responses against pathogens and abiotic stresses [48,49]. ABA is mainly related to plant defense against abiotic stresses. Environment factors such as drought, salinity, cold, heat stress, and wounding have been reported to trigger the increase of the ABA level [50,51]. SA plays an important role in response to biotic stresses as evidenced by the fact that pathogen infection can lead to an elevated level of SA [52]. Not all the stress-response signaling pathways are specific, and plenty of evidences had been provided for the cross-talk of ABA, SA with GAs in regulating plant defense responses [52,53]. GAs act throughout the whole plant life cycle to promote the growth of organs via enhancing cell division and elongation and to promote developmental phase transitions containing seed dormancy and germination, juvenile and adult growth phases, and vegetative and reproductive development [54,55].We also examined the expression profiles of the BrCB5s in response to exogenous GA3 ( Figure 6). In our study, we found that the majority of the BrCB5s were dramatically either up-or downregulated under the SA and ABA treatments and the expression of most BrCB5s was downregulated under the GA3 treatment.
BrCB5s can be used as electron donors in various oxidation/reduction reactions in cells. It has also been found that BrCB5s regulate the balance of active oxygen (ROS) in plants and are involved in response to a variety of biotic and abiotic stresses. Our results revealed that BrCB5s were more sensitive to NaCl and PEG 6000 treatments than heat and cold stresses. However, for BrVQs, the response seems to be the opposite [26]. We also found that the expression patterns of BrCB5s under NaCl and PEG 6000 treatments were significantly different. The expression levels of BrCB5s peaked in a few hours under the PEG 6000 treatment and then started to decrease, while the expression levels of BrCB5s continued to increase at 24 h when treated by NaCl. For the phytohormone treatment, we found that the expression of BrCB5s responded to SA, ABA, and GA3 stresses with different tendencies, indicating that BrCB5s might function in more than one signaling pathways and even mediate the cross-link of different pathways.
Additionally, some of BrCB5s under certain stresses were expressed at low levels, such as BrCB5e and BrCB5m under high-temperature treatment; it is also probably that they are not functional in Chinese cabbage without or with stress.

Conclusions
We identified 15 members of the Chinese cabbage CB5 gene family, which were classified into seven subfamilies. The phylogenetic relationships of the CB5s among Chinese cabbage, Rice, and Arabidopsis suggested that BrCB5s were more closely related to AtCB5s than OsCB5s. Phylogenetic and duplication event analysis suggested that whole genome duplication might be the main contributor to the expansion of the BrCB5s. The number of CB5 genes in Chinese cabbage was more than twice than that of AtCB5s. When treated with various abiotic stresses and hormone stimuli, differential expression of BrCB5s was observed. These data indicate that BrCB5s play an important role in plant growth and development. BrCB5s might mediate the cross-link between abiotic stresses and hormone signaling. Our work provides a basis for a further understanding of the characteristics and functions of the CB5 family in Chinese cabbage.

Identification and Comparison of CB5 Gene Family
Members in Chinese Cabbage. Bioinformatics methods were carried out according to our previous studies [25,26], which have been successfully applied to genome-wide analysis of GRF and VQ family genes in Chinese cabbage. Briefly, the nucleotide and protein sequences of BrCB5s were identified based on the B. rapa line Chiifu genome sequence (http:// brassicadb.org) [31]. The amino acid sequences were aligned by the software DNAMAN 6.0.40 (Lynnon Biosoft, Quebec, QC, Canada). Intron/exon structure analysis was performed by the Gene Structure display Server (GSDS) (http://gsds .cbi.pku.edu.cn/). The GC content was calculated by DNAS-TAR (Madison, WI, USA). The number of amino acids, molecular weight (MW), and theoretical isoelectric point (pI) were computed by the ProtParam tool (http://web. expasy.org/protparam/). Phylogenetic trees were constructed by the MEGA 5 software using the neighbor-joining method with 1000 bootstrap replicates.
The SSR markers were detected by the SSRIT software (http://archive.gramene.org/db/markers/ssrtool) with the parameters adjusted for identification of perfect di-, tri-, tetra-, penta-, and hexanucleotide motifs with a minimum of 6, 5, 4, and 4 repeats, respectively. The distribution of the conserved motifs and domains was detected by the MEME suite (http://meme-suite.org/tools/meme). The Arabidopsis and rice GRF protein sequences were downloaded from the Arabidopsis Information Resource (TAIR: http://www. arabidopsis.org/) and the Institute for Genomic Research Rice Genome Annotation project (TIGR: http://www. tigr.org/), respectively.

Plant
Growth and Treatments. Chinese cabbage cultivar "Zaohuangbai" plants were used in our study. The seeds were germinated in a glass petri dish with clean water at 20 ± 2°C for 24 h; then the glass petri dish was placed in a 4°C fridge for 15 days. After the vernalization, the seedlings were transferred into pots with soil and grown in a greenhouse at 20 ± 2°C with a 16 h light/8 h dark photoperiod. Artificial pollination was carried out after 7 weeks. 15 days after fertilization, samples were collected from root (R), stem (S), old leaf (OL, rosette leaves), young leaves (YL, cauline young leaves), flower bud (FB), and immature silique (IS) in three biological replicates for analyzing the expression of BrCB5s in different tissues.
For salinity, osmotic and hormone treatments, the seeds were germinated in a glass petri dish with clean water at 20 ± 2°C for 24 h. After germination, seedlings were transferred into pots with growth medium and grown in a greenhouse at 20 ± 2°C with a photoperiod of 16 h light and 8 h dark. Three-week-old seedlings of similar size were selected for the abiotic and hormone treatments. The plants were irrigated with 200 mM NaCl or 20% (w/v) polyethylene glycol (PEG 6000 ) while the CK group was watered with distilled water until the solution flowed out from the bottom of the pot. For temperature stresses, seedlings were transferred to incubators at 35°C and 4°C, respectively. For hormone treatments, the plant leaves were sprayed doubled-sided with 200 μM Gibberellin A3 (GA3), 100 μM Abscisic Acid (ABA), and 200 μM Salicylic Acid (SA) solutions, respectively, while leaves of the CK group were sprayed with distilled water until drops began to fall from the leaves' surface. The leaves of the treated seedlings were harvested after 0, 3, and 24 h of the above abiotic and hormone treatments.
All materials were immediately frozen in liquid nitrogen and stored at -80°C until RNA isolation.

RT qPCR Expression
Analysis of BrCB5s. Total RNA was extracted from each sample using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) and treated with RNase-free DNase I (Takara, Dalian, China) for 45 min according to the manufacturer's protocol. First-strand cDNA was synthesized using PrimeScript 1 st Strand cDNA synthesis Kit (Takara). Quantitative real-time PCR (qRT-PCR) was carried out using a SYBR Green Master mix (Takara, Dalian, China) on an IQ5 Real-Time PCR Detection System (Bio-Rad, Hercules, CA, USA). The qRT-PCR primers designed for the BrCB5s and actin gene (BrACT1) are listed in Supplementary file 6. Table S2, and the BrACT1 was proved to be successful in being the internal control in the research published by Lee et al. [56] in 2103. So, BrACT1 was used as a constitutive expression control in the qRT-PCR experiments. The PCR cycling conditions comprised an initial polymerase activation step of 95°C for 1 min, followed by 40 cycles of 95°C for 10 s and 60°C for 30 s. After each PCR run, a dissociation curve was designed to confirm the specificity of the product and to avoid the production of primer dimers. Three replicates of each sample were conducted to calculate the average Ct values. The relative expression level was calculated by the comparative 2 -ΔΔCt method.

Data Availability
The data used to support the findings of this study are included within the article.