Brassinazole Resistant 1 Activity Is Organ-Specific and Genotype-Dependent in Barley Seedlings

Brassinosteroids (BRs) control many plant developmental processes by regulating different groups of transcription factors, and consequently gene expressions. The most known is BZR1, the main member of the BES1 family. However, to date, it is poorly characterized in crop species. The main goal of the presented study was to identify HvBZR1 and determine its activity in 5-day-old barley (the stage is related to one leaf on the main shoot and a few seminal roots) using two cultivars with different sensitivities to BRs. Using the anti-OsBZR1 antibody, we identified the forms of HvBZR1 transcription factor with different molecular weights, which can be related to different phosphorylated forms of serine/threonine residues. Two phosphorylated forms in the shoots and one dephosphorylated form in the roots were determined. A minor amount of the dephosphorylated form of the HvBZR1 in the Haruna Nijo shoots was also found. The phosphorylated forms gave a higher band intensity for Golden Promise than Haruna Nijo. The bands were similar in their intensity, when two different phosphorylated forms were compared in Golden Promise, while a reduced intensity was detected for the phosphorylated form with a lower molecular weight for Haruna Nijo. Degradation of the phosphorylated forms in the shoots (complete degradation in Golden Promise and significant but not complete in Haruna Nijo) and the presence of the dephosphorylated form in the roots were proven for the etiolated barley. In the case of Haruna Nijo, a wider range of the regulators of the BR biosynthesis and signaling pathways induced the expected effects, 24-EBL (0.001 µM) and bikinin (10 and 50 µM) caused low amount of the phosphorylated forms, and at the same time, a tiny band of dephosphorylated form was detected. However, the expression of genes related to the BR biosynthesis and signaling pathways was not a determinant for the protein amount.

In arabidopsis (Arabidopsis thaliana L.), the last step of the BR signaling pathway leads to the regulation of Brassinazole Resistant 1 (BZR1) activity. The BZR1 is a DNA-binding protein, and regulates thousands of the BR-responsive genes [23][24][25]. A reduced transcript level of OsBZR1 led to the dwarf rice phenotype, its erected leaves, and reduced sensitivity to BRs [26]. TaBZR2 overexpression improved the drought tolerance (survival rate) of wheat as a reduced electrolyte leakage and malondialdehyde content. The higher proline content was observed during drought conditions [16].
A sequence identity search with the Basic Local Alignment Search Tool (BLAST) program in the UniProt Knowledgebase [54] for OsBZR1 (Q7XI96) [26,30,31,33,50] was performed. The F2CRW8 of Hordeum vulgare subsp. vulgare with an 83.6% identity of the sequence was found ( Figure S1). The bioinformatic analysis showed a high degree of identity between OsBZR1 and putative BZR1 orthologs in the other crop species, i.e., common wheat, 86.3-85.6%; foxtail millet (Setaria italica L.), 82.6%; broomcorn millet (Panicum miliaceum L.), 82.8%; sorghum (Sorghum bicolor L.), 78.7%; and brachypodium, 76.9%. The F2CRW8 consists of 317 amino acid residues (aa) and has BZR1 typical domains, i.e., DNA binding domain, 14-3-3 domain, PEST domain, and EAR motif ( Figure S1). The molecular weight (MW) of F2CRW8 protein is 33.9 kDa, and the isoelectric point (IP) is 8.0537; while the MW of OsBZR1 (298 aa) is 31.9 kDa, and IP is 8.4966 (Ensemble Plants). The F2CRW8 protein sequence is deposited with an accession number BAJ85589.1, and a coding sequence with accession number AK354370 in the National Center for Biotechnology Information (NCBI) or HORVU.MOREX.r2.2HG0108340 Ensemble Plants databases [55]. Using the BLASTN program, the nucleotide database was searched for the highly similar sequences (megablast) and the results indicated an 86.9% identity between H. vulgare subsp. vulgare AK354370 sequence and rice Os07t0580500-01 (Ensemble Plants, O. sativa Japonica Group, IGRSP-1.0 [56,57], older version Os07g39220.1 [26,34]), encoding OsBZR1 [26,30,31,33,50] ( Figure S2). The AK354370 sequence was used for the construction of the expression vector and F2CRW8 recombinant protein. The F2CRW8 recombinant protein was detected at the MW of~30 kDa ( Figure S3). This band corresponds to the dephosphorylated form of OsBZR1 [51,58,59] and SiBZR1 [60]. The Western blot analysis performed for the proteins extracted from shoots and roots of 5-day-old barley seedlings (Figure 1a) allowed for the identification of the presence of HvBZR1 in three different forms (Figure 1b). Two bands, thick one and thin one, located very close to each other, were obtained after analysis of the protein extract from the shoots of both cultivars, Golden Promise and Haruna Nijo. Another form represented by a thick band with a lower MW was found in the roots of both cultivars, and in a low amount in Haruna Nijo shoots. Small differences in the MW of the HvBZR1 phosphorylated forms can result from a diverse level of phosphorylation on serine and threonine residues, which in turn, is a major mechanism regulating the activity of cell proteins and playing a central role in signal transduction pathways [26,50,61]. The number of phosphorylation sites of serine and threonine residues is dependent on the species; there are 28 phosphorylation sites of serine and threonine residues in barley and rice ( Figure S1), and 21 in arabidopsis [26].
A few works were focused on the detection of the OsBZR1 protein in rice, while the HvBZR1 has not been previously studied in barley at all. Inactivation of the phosphorylated forms of the BZR1 by Constitutive Photomorphogenic 1 (COP1) in dark conditions and degradation of the dephosphorylated form of the BZR1 by SINA E3 ubiquitin ligases in the light were previously proven [20,62,63]. Therefore, we compared plants grown in the light with those grown in dark conditions (5-day-old etiolated plants) to prove F2CRW8 as HvBZR1 (Figure 2b). For the first time, degradation of the phosphorylated forms of the HvBZR1 in the shoots (complete degradation in Golden Promise, and significant, but not completed degradation in Haruna Nijo), and the presence of the dephosphorylated form of the HvBZR1 in the roots were proven for etiolated barley seedlings ( Figure 2b). The results indicated, that the HvBZR1 activity was organ-specific. Similarly, the identification of phosphorylated HvBZR1 forms corresponded to the HvGSK2.1 appearance in barley shoots (Figure 1b,o), while the root-specific dephosphorylated HvBZR1 form was related to HvGSK2.1 absence. This can be considered as a specific feature for Poaceae or crop species.
To link the influence of HvBZR1 activity with barley phenotype, the morphology of 5-day-old plants was assessed. All of the parameters describing the growth of shoots and roots showed a similar cultivar-dependent pattern (Figure 1a,c-n). The shoots of the Haruna Nijo were longer, and their fresh weight and dry weight were improved, when compared to the Golden Promise (Figure 1c-e). The opposite effect was found for the roots; their length and weight were significantly higher in the Golden Promise (Figure 1g-i). In all cases, the relative water content was similar (Figure 1f,j). The result indicates that despite the reduced root growth, which may directly affect the efficiency of water and nutrient uptake, the Haruna Nijo was able to develop much stronger shoots (Figure 1k-n). These traits differentiated both genotypes and suggested that growth changes could be related to the difference in the dephosphorylated HvBZR1 form in the shoots of both cultivars (a small amount in the Haruna Nijo, and its absence in the Golden Promise, Figure 1b). (c-f), root growth (g-j), and shoot-to-root ratio (k-n); (o) immunodetection of HvGSK2.1 kinase and α-tubulin as a control; (p,q) expression profile of the genes related to the BR biosynthesis and signaling pathways. The results present the mean with a standard error of the mean (phenotypic traits n = 6, (c-n); expression profile n = 3, (p,q). The asterisks indicate significant differences, revealed by the Student's t-test for p < 0.05 (*), p < 0.01 (**), and p < 0.001 (***). (c-f) growth parameters of shoots, (g-j) roots, and (k-n) shoot-to-root ratio. The results present the mean with a standard error (n = 6). The asterisks indicate significant differences, revealed by the Student's t-test for p < 0.05 (*), p < 0.01 (**), and p < 0.001 (***).
The BZR1 activity is regulated by the BIN2 as well as other GSKs in the BR signaling pathway [26]. Therefore, we determined the expression of key genes in this pathway (Figure 1p,q). The BL molecules are used to initiate the BR signaling pathway. Therefore, we also determined the expression of the first gene controlling the BR biosynthesis, i.e., HvDWF4. The DWF4 expression is regulated by BZR1, and the DWF4 activity is inhibited by Brz [64]. HvBRI1 and HvBAK1 encode a transmembrane receptor BRI1-BAK1, which captures the BR molecules with the BL structure (i.e., 24-EBL). HvBSU1 encodes the phosphatase BSU1, which is regulated by the BR signaling pathway and inhibits proteins from the GSK3s family. Moreover, we analyzed the expression of genes encoded GSK3s family members [46], and finally the HvBZR1 expression. The genes encoding enzymes that control the subsequent steps of the BR signaling pathway were expressed in cultivar-and organ-dependent mode (Figure 1p,q). The Real-Time PCR analysis allowed to distinguish four types of gene responses. The first group consisted of the genes in which the expression was organ-dependent in the Golden Promise, but at the same time, they did not change in the shoots and roots of the Haruna Nijo. These were the genes specifically upregulated in the roots (HvDWF4, HvBRI1, HvGSK1.1) and in the shoots (HvBSU1, HvGSK2.1, HvGSK3.1, HvGSK4.1) of the Golden Promise. The second group consisted of genes upregulated in Haruna Nijo roots (HvBAK1 and HvGSK2.2), while organ-independent in the Golden Promise. The third group included two genes with a similar expression in the organs of both cultivars (HvGSK1.2 and HvGSK1.3). However, one of them (HvGSK1.3) was clearly cultivar-dependent, as a significantly higher expression was detected in both the shoots and roots of the Golden Promise. The HvBZR1 represented the fourth type of response, as the upregulation in the roots, but not the shoots of both cultivars, was found. The result for HvBZR1 was consistent with data for other plant species [61].
Despite the relatively low expression of HvBZR1 (compared to the other genes, Figure 1p

Barley Response after 24-EBL, Bikinin, and Brz Treatment
To determine the contribution of the phosphorylated and dephosphorylated forms of HvBZR1 in the regulation of plant growth, the Golden Promise and Haruna Nijo cultivars were treated with 24-EBL, bikinin, and Brz at the two concentrations. This differentiate the shoot and root phenotypic traits between the cultivars in our previous study [18]. The 24-EBL and bikinin lead to the inhibition of GSK3s activity, and consequently promote the dephosphorylation of BZR1 and its nuclear localization [50,61]. In addition, 24-EBL inhibited the GSK3s activity through the activation of the BR signaling pathway, while bikinin through the inhibition of LEYV (specific for the second group of the GSK3s family members) and MEYV motif (specific for the first and third groups of GSK3s family members) [48]. On the contrary, the Brz treatment leads to an inhibition of the BR biosynthesis pathway and impaired BR signaling pathway, thus allowing for the GSK3s to be active, and consequently promotes the phosphorylation of BZR1 and its cytoplasmic localization. In each case, the concentration of the active substance was selected based on the earlier studies, which focused on their dose-dependent regulation of barley growth [18]. The 0.001 µM 24-EBL and 10 µM bikinin stimulated root elongation, while 1 µM 24-EBL, 50 µM bikinin, 10 µM Brz, and 50 µM Brz inhibited barley shoot and root development (Figure 2a,c-n).
In the case of each treatment (24-EBL, bikinin, and Brz), the bands typical for the untreated control plants were detected, i.e., two phosphorylated forms of the HvBZR1 in a large amount in shoots of both cultivars, a major amount of dephosphorylated form in roots of both cultivars, and a minor amount of dephosphorylated form in shoots of the Haruna Nijo ( Figure 2b). A more detailed analysis of bands showed differences in the amount of HvBZR1 between the cultivars and after the treatment with active substances (Figure 2b). In shoots, phosphorylated forms had higher band intensities for Golden Promise than Haruna Nijo. The bands were similar in their intensity when two phosphorylated forms were compared in Golden Promise, while a reduced intensity was detected for the phosphorylated form with a lower MW for Haruna Nijo. Similarly, in roots, the dephosphorylated form with a higher content for Golden Promise than Haruna Nijo was found. The expected (similar as presented in the literature for other species) response in shoots was only observed after the treatment of Golden Promise plants with bikinin used in the low concentration (10 µM). A low amount of phosphorylated forms was noted, even if dephosphorylation was not observed (absence of bands related to the dephosphorylated form). In the case of Haruna Nijo, the wider range of treatments induced the expected effects. 24-EBL in a low concentration (0.001 µM) and bikinin at both concentrations (10 and 50 µM) caused the low amount of phosphorylated forms, and at the same time, the tiny band of dephosphorylated form was detected. In the same way, Brz at both concentrations (10 and 50 µM) induced an accumulation of phosphorylated forms and promoted the absence of dephosphorylated form in Haruna Nijo shoots. In Golden Promise roots, the expected response was found after the 1 µM 24-EBL treatment, i.e., a high level of dephosphorylated form. For the other treatments, the effects were relatively comparable to the untreated control, particularly in Haruna Nijo root samples. The shoot-to-root ratio of both cultivars increased after the treatment with 1 µM 24-EBL and 50 µM bikinin, but decreased when a low concentration of the active substances was used. Brz led to genotype-dependent changes in the shootto-root ratio, i.e., increased parameters in Golden Promise and decreased in Haruna Nijo. Altogether, the results proved the organ-specific (occurrence of phosphorylated forms of the HvBZR1 in the shoots, and dephosphorylated form of the HvBZR1 in the roots) and genotype-dependent (major changes in HvBZR1 for Haruna Nijo shoots, minor changes in HvBZR1 for Golden Promise shoots) activity of HvBZR1.
Following the changes in HvBZR1 activity after the different treatments, the plant growth changes induced by active substances were also organ-specific and genotypedependent. In the first case, roots were more sensitive to different concentrations of active substances than shoots, particularly to 24-EBL and bikinin. The pattern of changes was similar in both cultivars, in which a low concentration increased the growth parameters and a high concentration decreased them. In the second case, shoot growth was more dependent on the genotype. The Golden Promise had decreased the shoot growth after the 24-EBL and bikinin treatment, in contrast to the improved growth of the shoots of Haruna Nijo.
In addition, a detailed statistical analysis proved that all of the measured traits in shoots and only dry weight in roots were genotype-dependent (Table 1). Following the results, we performed a quantitative analysis of the genes controlling consecutive stages of the BR signaling pathway (Figure 3). 24-EBL and bikinin at the high concentration reduced expression of the genes in barley roots, e.g., HvBSU1, HvGSK2.1, and HvGSK2.2. Similarly, the 10 µM bikinin resulted in a significant reduction of the expression of some genes in the roots of Haruna Nijo (HvBAK1, HvBSU1, and HvGSK1.2) and in Golden Promise (HvGSK2.1). In the shoots, the 50 µM 24-EBL reduced expression of the genes in Golden Promise, while most of the genes were not altered in Haruna Nijo. In Golden Promise, a few genes were downregulated by Brz, which was applied in the high concentration. Two-way ANOVA indicated that most of the genes (with the exception of the HvBRI1) in shoots and only HvGSK1.1 and HvGSK2.1 in the roots are genotype-dependent ( Table 1). The results indicated that the following steps of the BR signaling pathway can be organ-specific and genotype-dependent in 5-day-old barley seedlings. Table 1. Results of two-way ANOVA for the genotype (Golden Promise and Haruna Nijo) or treatment (0.001 and 1 µM 24-EBL, 10 and 50 µM bikinin, 10 and 50 µM Brz) calculated for different growth parameters (shoot and root length, their fresh and dry weight, and related shoot-to-root ratio) and genes (HvDWF4, HvBRI1, HvBAK1, HvBSU1, HvGSK3, HvGSK1.1,  HvGSK1.2, HvGSK1.3, HvGSK2.1, HvGSK2.2, HvGSK3.1, HvGSK4.1, HvBZR1). The asterisks indicate significant dependence for p < 0.05 (*), p < 0.01 (**), and p < 0.001 (***).

Discussion
In the last few years, the role of GSK3s and BZR1 was extensively studied in rice architecture and grain formation [12,51,60,65]. However, the proteins were rarely studied at the early developmental stages and during stress responses. Similarly, few articles have provided information regarding GSK3s in barley [18,46], and there is no available data concerning HvBZR1. In our study, we present the organ-specific HvBZR1 activity and genotype-dependent responses to the active substances, which modify the BRs biosynthesis and signaling pathways, i.e., 24-EBL, bikinin, and Brz.
The barley cultivars Golden Promise and Haruna Nijo displayed a different sensitivity to 24-EBL, bikinin, and Brz in the lamina joint inclination test [18]. They had different phenotypes, most probably as a result of differences in the BR-related pathways. Some works indicated that the lamina joint depends on the endogenous BRs level [51,60]. Our experiments proved that different phenotype of both cultivars corelate with different levels of HvBZR1 and HvGSK2.1 proteins. HvGSK2.1 is a homolog of arabidopsis BIN2 and rice OsGSK2 protein kinases [46]. Here, the most important result consists of a different pattern of the HvBZR1 forms between both cultivars. Similar to OsBZR1 that was detected in the protoplasts of wild-type rice plants [50], we detected two phosphorylated forms and one dephosphorylated form of HvBZR1 in barley with different MW levels. Differences in MW correspond to various levels of phosphorylation in serine and threonine residues [26,51]. Furthermore, in the dark condition, phosphorylated forms of the HvBZR1 degraded completely in the shoots of the Golden Promise, while a significant but not completed degradation was found in the Haruna Nijo.
The treatment of a 5-day-old barley plant with 24-EBL and bikinin was expected to promote the dephosphorylated HvBZR1 form, and the Brz treatment to promote the phosphorylated HvBZR1 forms. Unexpectedly, we found that the HvBZR1 activity was strongly organ-dependent, and the HvBZR1 divergent pattern was found between the shoot and root tissues (Figure 2b). The phosphorylated forms were detected in the shoots and the dephosphorylated form in the roots. Their amount was only slightly modified by the active substances. The result was strengthened by the HvGSK2.1 analysis. HvGSK2.1 was present in the shoots, but it was undetected in the roots of both barley cultivars. The activity of PUB40, the U-box E3 ubiquitin ligase, led to the degradation of both forms of the BZR1 in arabidopsis roots, but not in the shoots, when BRs were present in the low amount [61], and OsPUB24 induced OsBZR1 proteasomal degradation by ubiquitination of OsBZR1 in rice [50]. Both ligases required phosphorylation by OsGSK2/BIN2 to induce their activity. Earlier studies with protoplasts allowed for the study of OsBZR1 changes only in the shoots [50]. For the first time, our study proved the organ-specific location and response of HvBZR1 in barley shoots and roots at the early stage of plant development.
A higher amount of the phosphorylated HvBZR1 forms in shoots than roots and the difference between two phosphorylated forms in Golden Promise (an equal amount of both forms) and Haruna Nijo (a lower amount of the form with lower MW), may be related to the HvBZR1 specific role in the regulation of plant aboveground development. In the case of Haruna Nijo, the wider range of treatments induced the expected effects, 24-EBL and bikinin caused the low amount of phosphorylated forms and at the same time, the tiny band of dephosphorylated form was detected. The pleiotropic effect was suggested in the case of bikinin [48]. Indeed, results consisting of the significant degradation of phosphorylated forms of the HvBZR1 in Haruna Nijo, its stability in the Golden Promise plant, and GSK3s inhibition. Altogether, it raises the following questions: 1. Is the HvBZR1 activity regulated by GSK3s?; 2. and/or is the regulation determined by the affinity of bikinin for the bikinin-binding domains, i.e., LEYV in GSK group II (GSK2.1 and GSK2.2), MEYV in group I (GSK1.1, GSK1.2, and GSK1.3), and III (GSK3.1) ? The hypothesis was partially supported by the results of gene expression. We observed a lower expression of the HvBZR1 gene in the shoots than the roots of barley. Similarly, a higher expression of BZR1 was observed in arabidopsis and rice roots [50,61]. However, the expression of genes related to BR biosynthesis and signaling pathways, in our case, was not a determinant for the protein amount. For instance, the upregulation of HvGSK2.1 was not followed by the improved amount of HvGSK2.1 protein in barley roots.
Recent studies with arabidopsis mutants related to individual steps of the BR biosynthesis pathway indicated a highly site-specific BRs function and short-distance signal transduction in roots [66]. The nuclear localization of BZR1 was specific for the root meristem zone, and its role was assigned to cell division, particularly the Quiescent Center (QC), Columella Cell (CC), and Columella Stem Cell (CSC) [67,68]. The QC of the root meristem is under control of the Brassinosteroids at the Vascular and Organizing Centre (BRAVO) TF, inhibiting the QC division [69]. In addition, 24-EBL or bikinin treatment can lead to GSK3 deactivation, thus the activation of BZR1, and inhibition of BRAVO. As a result, this can induce cell division and root elongation. Indeed, both barley cultivars had elongated roots after the treatment with low concentrations of 24-EBL and bikinin, and the elongation resulted from the cell division rather than from the cell elongation (improved fresh and dry weight, stable RWC). On the contrary, roots phenotype changes after the plant treatment with a high concentration of 24-EBL and bikinin were consistent with the GSK3 role in the regulation of auxin signaling pathway [68,70,71], and changes in the PIN-Formed 2 (PIN2) level [72,73]. The changes in Golden Promise and Haruna Nijo root phenotypes suggest that the application of 24-EBL and bikinin, but not Brz, led to a modification in auxin signaling, but the hypothesis should be verified in the future.
The presented results confirmed that effect of 24-EBL, bikinin, and Brz on barley development is dose-and genotype-dependent as well as organ-and developmental stage-specific ( Figure 2) [74,75]. However, considering our results and the other studies, the important differences in both, the role of HvBZR1 in plant growth and the effect of active substances in the modification of the BR biosynthesis and signaling pathways, were recognized when barley (monocotyledon) was compared to arabidopsis (dicotyledon). In conclusion, similar phenotypical changes after the bikinin and 24-EBL treatments suggest a role of GSK3s in the regulation of HvBZR1 activity. HvBZR1 can be localized in the nucleus of root cells, and its activity can lead to the BR-dependent regulation of gene expression. In barley shoots, HvBZR1 can be localized mainly in the cytoplasm. The phosphorylation or dephosphorylation of serine and threonine allow for the interaction of HvBZR1 with other proteins. Finally, phenotype changes are also dependent on the different plant sensitivity to BR, as we observed for Golden Promise and Haruna Nijo. However, further work is required to determine the factors influencing the plant sensitivity to BRs, and the role of dephosphorylated and phosphorylated forms of the BZR1.

Treatments and Growth Traits
Grains were germinated in the presence of 0.001 and 1 µM 24-EBL, 10 and 50 µM bikinin, 10 and 50 µM Brz, 0.7% EtOH (control for 24-EBL), and 0.7% DMSO (control for bikinin and Brz). Plants were cultivated as previously reported [18]. The 5-day-old barley shoots and roots (the stage is related to one leaf on the main shoot and a few seminal roots) were collected between 10:00 and 12:00. Plants were photographed, and the shoot and root length were measured with the ImageJ software (v. 1.51k). The shoot and root length were measured between the mesocotyl and the 1st leaf tip and between the mesocotyl and the tip of the longest primary root, respectively. Plants were weighed, the fresh (FW) and dry weight (DW) were used to calculate the relative water content (RWC) according to the formula: RWC = (Fw − Dw/Fw) × 100%.

RNA Extraction, cDNA Synthesis, and Real-Time PCR Analysis
Total RNA was extracted using the TRI Reagent Solution (Life Sciences Solutions, Carlsbad, CA, USA). The cDNA was synthesized using the Maxima H Minus First Strand cDNA Synthesis Kit with dsDNase (Thermo Fisher Scientific, Vilnius, Lithuania) with oligo(dT) as primers. Real-Time PCR was carried out using the 5 × HOT FIREPol EvaGreen qPCR Mix Plus, no ROX (Solis BioDyne, Tartu, Estonia) and Rotor-Gene 6000q series (Corbett Life Science, Mortlake, Australia) according to the manufacturer's protocol. The ADP-ribosylation factor (AFR) and Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) were used as internal controls. The gene-specific primers used for the Real-Time PCR are listed in Table S1. The relative expression of the genes was calculated from the number of copies of the investigated genes and the geometric mean of the reference genes, i.e., ARF and GAPDH, which were determined using the Rotor-Gene 6000 series software 1.7 (Corbett Life Science, Chadstone, Australia).

Data Analysis
Statistical analyses (Student's t-test) were performed using the Microsoft Excel 2016 (Microsoft Office, Warszawa, Poland) and the statistical software for the Excel, i.e., XLSTST (Addinsoft, Paris, France). Two-way ANOVA were performed using the Statistica 13.0 (StatSoft, Kraków, Poland). All of the phenotypical traits are shown as the results for six biological replicates. Three biological replicates were performed for each gene and three technical repeats for each sample. The mean values with a standard deviation were presented. Graphs were generated using the Microsoft Excel 2016 and the Microsoft PowerPoint 2016 (Microsoft Office, Warszawa, Poland). Results of the experiments are presented in Table S2.

Conclusions
The BZR1 in plants was detected in three forms with different MWs corresponding to the phosphorylation level of serine/threonine residues. Additionally, the activity of HvBZR1 was organ-specific. Phosphorylated forms of the HvBZR1 with the cytoplasmic localization were characteristic for the shoots. In contrast, HvBZR1 in roots was present as a dephosphorylated form with nuclear localization. The changes in the accumulation of dephosphorylated and phosphorylated forms of HvBZR1 were observed between the Golden Promise and Haruna Nijo cultivars.