1. Introduction
Brassinosteroids (BRs) are a class of plant polyhydroxylated steroid hormones, which are involved in many developmental processes and stress responses. Brassinolide (BL), an end-product of the BR biosynthetic pathway, is a signaling molecule recognized by Brassinosteroid Insensitive 1 (BRI1) [
1]. The binding of BL by BRI1 leads to conformational changes in this receptor and the disconnection of BRI1 Kinase Inhibitor 1 (BKI1) [
2], resulting in its association with BRI1 associated receptor kinase 1 (BAK1) [
3,
4]. The BRI1-BAK1 transmembrane receptor initiates a cascade of phosphorylation and dephosphorylation of cytoplasmic relay proteins, leading to dephosphorylation and the inactivation of kinases from the Glycogen Synthase Kinase 3 (GSK3) family [
5,
6,
7]. The GSK3 family controls many transcription factors (TFs), e.g., those regulating cell elongation and cell division, root meristem and root development, lateral root development, stomatal development, xylem differentiation, vascular cambial activity, chloroplast development, photomorphogenesis, hypocotyl elongation, floral organ development, flowering, seed development, fruit ripening, sugar content in fruits, and responses to stress (salt, drought, cold, and biotic) [
8]. The GSK3s family is represented by four kinase groups [
9]; nine are known to date in rice (
Oryza sativa L.) [
10], seven in barley (
Hordeum vulgare L.) [
11], eleven in maize (
Zea mays L.) [
12], and ten in Arabidopsis (
Arabidopsis thaliana L.) [
13]. At low BR concentrations, GSK3 phosphorylates Brassinzaole Resistant 1 (BZR1), one of the major TFs, leading to its proteasomal degradation [
14]. The activity of BZR1 is regulated by 14-3-3 protein [
15,
16] and Protein Phosphatase 2A [
17].
The 4-[(5-bromo-2-piridynyl)amino]-4-oxobutanoic acid (bikinin) compound has been identified among the library of 10,000 compounds as one that induces constitutive BR responses in Arabidopsis, i.e., a significant increase in hypocotyl length; long and bending petioles; and blade-shaped, pale-green leaves, as well as effects comparable to those obtained with treatment with BL at micromolar concentrations, in which lateral root density was reduced [
18]. Consecutive analyses with bikinin and its inactive variant showed that treatment induced petiole and hypocotyl elongation under light and dark growth conditions, allowing researchers to conclude that Brassinosteroid Insensitive 2 (BIN2–the best characterized kinase from the GSK3 family) is a direct target of this compound. Due to its role in inactivating BIN2, the name of bikinin was derived from its role in the BR signaling pathway, i.e., BIN2 kinase inhibitor. Bikinin folds the lamina joint of barley, similarly to 24-epibrassinolide (24-EBL), and has resulted in a similar phenotype under optimal and salt stress conditions [
19].
Previously, BRs have been shown to regulate the fruit ripening of tomato (
Solanum lycopersicon L.) [
20,
21], persimmon (
Diospyros kaki L.) [
22], grapevine (
Vitis vinifera L.) [
23], pear (
Pyrus ussuriensis Maxim.), and apple (
Malus domestica Borkh.) [
24] plants, as well as affecting the yield of crops by regulating grain formation, which has been the best characterized in rice [
25,
26,
27,
28,
29]. Moreover, BRs have been shown to positively regulate salt stress tolerance in many species, e.g., canola (
Brassica napus L.) [
30], brown mustard (
Brassica juncea L.) [
31], pepper (
Capsicum annuum L.) [
32], pea (
Pisum sativum L.) [
33], tomato [
34], potato (
Solanum tuberosum L.) [
35], barley [
19], rice [
36], and wheat (
Triticum aestivum L.) [
37]. The regulatory role of BRs has been investigated by applying bioactive 24-EBL [
38], BL, and 28-homobrassinolide [
39] or by using mutants with the genetic dissection of the BR biosynthetic pathway [
40,
41,
42,
43,
44]. However, the role of GSK3s in barley development using bikinin has not yet been reported. Here we present the results of the application bikinin to barley plants and the long-term effects of GSK3 family inhibition on barley growth under high-salinity conditions. According to the results presented by Honda et al. [
44] and Groszyk and Szechyńska-Hebda [
45], we considered that the best approach to analyze GSK3s function in barley development was to use a GSK3s inhibitor, i.e., bikinin. We assumed that because of its specificity for GSK3s [
46] and its lack of toxic effects, as with another GSK3s inhibitor, i.e., lithium [
47,
48,
49], we would be able to characterize the initial molecular and physiological processes in which GSK3s are involved. The analyses were performed during the outset of the formation of side shoots phase, when tillers were still invisible, i.e., the BBCH20 growth stage according to the BBCH-scale, which is used to identify the phenological development stages of plants [
50]. Therefore, the aim of this study was to understand the long-term response of barley plants exposed to salinity to bikinin treatment.
3. Discussion
Crops grown under stress conditions show restricted growth and, eventually, reduced yields. Among other treatments, brassinosteroids (BRs) mitigate the effects of stress and improve plant growth. In this study, the following hypotheses were verified: (1) barley’s (Hordeum vulgare L.) tolerance of salinity depends on its sensitivity to BRs and the activity of the transcription factor (TF) HvBZR1; (2) treatment of barley with bikinin lowers kinase suppression of HvBZR1, which can be observed through a reduction in the amount of the phosphorylated form of HvBZR1; and (3) the use of two barley varieties with differing responses to salinity allowed the identification of regulatory elements of the BR biosynthesis pathway.
To test our hypotheses, we used two barley cultivars, i.e., Golden Promise and Haruna Nijo, with differing sensitivity to BRs, which was determined by means of the lamina joint inclination test, in which Haruna Nijo presented greater sensitivity [
19]. Barley plants with a 2nd unfolded leaf of the same length or ~1 cm longer than the 1st one were sprayed with a diluted series of bikinin, in parallel with the start of NaCl treatment. Barley was grown under controlled conditions until the beginning of the 5th leaf emergence stage. The 3rd and the 4th leaves were used for molecular analyses to determine the effects of the applied inhibitor. Plants treated with 0.11% DMSO were used as a baseline for plants treated with bikinin solutions. Non-treated plants were used as a benchmark for salinity. Similarly to previous studies [
19], barley length and fresh weight depended on the genotype’s response to growth conditions, and early leaf length and area were characteristic of the genotype (
Table 1). In contrast, the consecutive leaves’ traits depended on the growth conditions (
Table 1). Barley’s response to bikinin treatment was genotype-dependent and the greater (albeit insignificant) changes were observed in Haruna Nijo treated with higher concentrations of bikinin (
Figure 1 and
Figure 2). This was also observed in the growth of the 4th leaf (
Figure 1g,h). In contrast to Haruna Nijo, Golden Promise had a faster phenotypic response to bikinin treatment, as observed in the increased length of the 3rd leaf (
Figure 1g). These traits may depend on the varieties’ sensitivity to BRs. Despite the greater sensitivity to 24-EBL and bikinin observed in Haruna Nijo, greater changes in lamina inclination were observed at higher concentrations in Golden Promise [
19]. BR application promoted hypocotyl and epicotyl elongation of soybeans (
Glycine max L.) [
52], petioles of carrot (
Daucus carota L.) [
53] shoots and roots of barley [
19].
Bikinin at the highest concentration also increased the fresh and dry weights of roots in Haruna Nijo under saline conditions, resulting in root lengths similar to plants under optimal conditions (
Figure 1). When 5-day-old plants were analyzed, the roots of Haruna Nijo under conditions of salinity were found to be longer than those of Golden Promise [
19]. In another study, these changes were associated with the characteristics of root BZR1 transcription factor activity, which has been shown in many studies to be a factor regulating cell division in the Quiescent Center (QC), Columella Cells, and Columella Stem Cell [
45,
54,
55,
56]. The QC of the root meristem is controlled by Brassinosteroids at the Vascular and Organizing Center (BRAVO) [
57]. The increased cell division levels observed in Haruna Nijo roots after bikinin treatment, resulting in higher fresh and dry weights, may be associated with a higher cell number, but this should be verified in the future. Root phenotype may be connected with the role of Glycogen Synthase Kinase 3 (GSK3) family in the regulation of the auxin signaling pathway [
56,
58,
59] and signal transduction due to the fact that the activity of the
HvGSK2.1 kinase from the GSK3 family has not been detected in barley roots [
45]. However,
AtSK11 and
AtSK12 have been characterized as a genes that induce root growth under osmotic stress [
60], and the best-characterized kinase from GSK3 family, i.e., Brassinosteroid Insensitive 2 (BIN2), has been detected in Arabidopsis roots [
61,
62]. Orthologs of these genes have been identified in barley [
11] and the expression of
HvGSK1.1 in Golden Promise roots was the highest compared to that of other genes from this family [
45]. The first functional analysis of
HvGSK1.1 showed that reduced expression by RNAi resulted in higher weights of transgenic plants under conditions of salinity [
63]. However, the figure presented in the article shows that the authors of that study compared barley in different developmental stages [
63] and that the inhibition of kinase from GSK3 family led to faster growth. Analysis of
OsGSK1 (also known as OsGSK21, OSKζ) with the Os01g10840 locus [
10,
64,
65], an ortholog of BIN2 belonging to a class represented by four orthologs in rice [
10] and two in barley [
11], showed that knockout mutants had enhanced tolerance to cold, heat, salt, and drought stresses [
65]. In contrast, the overexpression of full-length
OsGSK1 led to a stunted growth phenotype, similar to that of the gain-of-function
bin2 mutant [
65]. Rice with knock-out of
OsGSK1 exhibited a lower wilting ratio and improved Fv/Fm compared to controls under the same conditions [
65]. Similarly, with increasing bikinin concentrations, both barley genotypes had higher relative water content and relative turgidity and lower water deficits in consecutive concentrations than in the controls grown under saline conditions (
Figure 3). Despite the changes in both genotypes, the water content was genotype-dependent and depended on the control conditions and the genotypes’ responses to bikinin treatment (
Table 3 and
Table 4). On the other hand, photosystem II (PSII) activity was genotype-dependent and depended on growth conditions and barley’s response to salinity (
Table 5) but changes in consecutive parameters depended on the genotypes’ responses to bikinin treatment (
Table 6). Salinity-induced changes were alleviated by bikinin (e.g., PI ABS, DIo/RC in Golden Promise; Fv/Fo, φDo in Haruna Nijo). A positive effect of BR treatment on PSII activity has been observed in cucumber (
Cucumis sativus L.) [
66,
67], wheat (
Triticum aestivum L.) [
68,
69], soybean (
Glycine max L.) [
70], maize (
Zea mays L.), spinach (
Spinacia oleracea L.) [
71], and mung beans (
Vigna radiata) [
72]. Differentially expressed gene (DEG) analysis revealed an opposite and genotype-dependent response. 9 out of 117 identified DEGs in Golden Promise and 24 out of 150 DEGs identified in Haruna Nijo played a role in photosynthesis (
Figure 8). Previously, the GSK3 family was identified to contain kinases controlling stomata development via the regulation of YDA, MKK4/5, and Speechless [
64,
73,
74,
75,
76].
The short-term response (30 min or 2 h) to BL or bikinin led to the identification of 272 genes involved in BR metabolism, BR biosynthesis, hormone-mediated signaling, auxin, and the response to abiotic stimuli [
18]. The long-term response (approximately 17 days) resulted in the regulation of genes involved in photosynthesis (
Figure 7 and
Figure 8). The differences indicate that the bikinin response led to rapid changes in phytohormonal regulation that regulated plant phenotype, mainly affecting photosynthetic efficiency and water content in the days following plant growth (
Figure 3 and
Figure 7).
The most distinctive results were obtained in the response to two concentrations of bikinin. Although the transcript level in Golden Promise increased with higher bikinin concentrations, it decreased in Haruna Nijo and vice versa (
Figure 8). However, photosynthetic activity parameters and transcriptome analysis data confirmed the genotype-dependent response to bikinin treatment. According to a previous study [
45], we identified a phosphorylated form of
HvBZR1 in the 3rd and the 4th barley leaves that accumulated more under conditions of salinity and the amount of which was reduced during treatment with 50 µM and 100 µM bikinin in Golden Promise and after each bikinin concentration in Haruna Nijo (
Figure 9). As before, the results indicated two levels of phosphorylation with a lower amount of
HvBZR1 over a lower molecular weight. Both proteins,
HvBZR1 and
HvGSK2.1, detected in leaves showed greater amounts than in 5-day-old barley shoots [
45]. In contrast, the expression profiles of these genes showed an opposite response. The low expression profile of
HvBZR1, accompanied by the high
HvBZR1 protein level, suggests that this TF is stabilized through the phosphorylation process and stored in the cytoplasm. To date, many security proteins for BZR1 have been identified [
16,
77,
78,
79]. Contrary to the highest expression level of the GSK3 family, protein accumulation in leaves was lower than for
HvBZR1. This suggests that
HvGSK2.1 may be an unstable protein with short-term activity. Similarly to
HvGSK2.1 and
HvBZR1, the expression levels of genes controlling the consecutive stages of the BR signaling pathway were reduced by salt stress, but their expression was stable and unchanged after bikinin treatment.
The phenotypic changes observed in the two barley genotypes grown under saline stress confirmed previous results obtained for the same genotypes in a different experimental design. Despite the presence of shared traits, many data suggested a genotype-dependent, exogenous BR-linked response by barley plants to bikinin and the consequent inhibition of GSK3 activities, affecting plant development. Nevertheless, many aspects remain unexplained and further studies are needed to further characterize barley. In summary, our results indicate that barley’s salt tolerance depends on its sensitivity to BRs, but the activity of the HvBZR1 TF should be verified. We confirmed that barley treatment with bikinin reduced the amount of the phosphorylated form of HvBZR1. The use of two barley varieties with different responses to salinity led to the identification of the role of BR signaling in photosynthesic activity. The results suggested that salinity reduces the expression of the genes controlling the BR signaling pathway. Moreover, the results also suggested that the functional analysis of GSK3 in stress responses can be a tool for the breeding of crops to improve plant resistance to salinity or other stresses.
4. Materials and Methods
4.1. Chemicals
Bikinin (CAS 188011-69-0, purity ≥98%) was purchased from Sigma-Aldrich (Schnelldorf, Germany). Solutions of bikinin (10 µM, 50 µM, and 100 µM) were prepared from 91.5 mM stock, dissolved in 100% dimethyl sulfoxide (DMSO) (Sigma-Aldrich, Schnelldorf, Germany). The controls used in the experiments contained the same concentration of the solvent solution, and these were used as a background for the dilution of bikinin.
4.2. Plant Material
The barley (
Hordeum vulgare L.) cultivars Golden Promise (United States Department of Agriculture, GRAIN-Global, USA, accession number 343079) and Haruna Nijo (Gene Bank Dept., CRI Prague-Ruzyně, accession number 03C0602168) were used in the experiments [
45]. Grains were provided from both Gene Bank and imbibed in Petri dishes with three layers of filter paper soaked with spring water (Żywiec-Zdrój S.A., Węgierska Górka, Poland) for 48 h in a refrigerator at 4 °C, then germinated in darkness for 72 h in an incubator at 23 °C, then 8 plants were planted in 12 L buckets (23 cm × 33 cm × 19 cm) filled with solid substrate (Hollas, Pasłęk, Poland) and sand (ratio 4:1), and grown in a greenhouse (October 2019–January 2020) for ~3 months in the case of Haruna Nijo and ~4 months in the case of Golden Promise until complete harvest. Grains from each plant were collected separately and studied as single-seed descent (SSD) lines.
4.3. Growth Conditions
SSD lines were used for the experiment. Grains were imbibed and germinated as described above. Then, eight seedlings were planted in 1 L pots filled with soil substrate (Hollas, Pasłęk, Poland) and sand (2:1). Plants were cultivated in a phytotron chamber at a 16 h photoperiod at 20 °C during the day and 18 °C at night, with a daylight intensity of 200 µmol photons m2 s−1, and a humidity of 70%. Plants were watered to a soil humidity of 70%, with the growth substrate fully watered (as 100%) and fully desiccated (as 0%). At the stage when the 2nd leaf had a similar length to the 1st leaf, sodium chloride application and bikinin treatment were performed. Sodium chloride at a final concentration of 150 mM was administered every second day 3 times with 50 mM NaCl once. Then, pots were watered with tap water to 70% humidity and measured using technical scales and weighed up to 1051 g (fully watered, 1152 g; fully desiccated, 816 g; pot weight, 31 g). Plants watered with tap water only were used as controls and grown under optimal conditions.
4.4. Barley Treatment
Barley at the stage described above (1st and 2nd leaf of approximately equal length) was treated with 10 µM (BK10), 50 µM (BK50), or 100 µM bikinin (BK100) or 0.11% DMSO as a control solvent solution (CK). All solutions were prepared as follows:
BK10: 0.98 µL 91.5 mM bikinin, 8.85 µL 100% DMSO, 25 µL Tween 20 (Sigma-Aldrich, Schnelldorf, Germany), 9 mL deionized water;
BK50: 4.92 µL 91.5 mM bikinin, 4.91 µL 100% DMSO, 25 µL Tween 20, 9 mL deionized water;
BK100: 9.83 µL 91.5 mM bikinin, 25 µL Tween 20, 9 mL deionized water; or
CK (0.11% DMSO): 9.83 µL 100% DMSO, 25 µL Tween 20, 9 mL deionized water.
Barley at the 5th leaf development stage (BBCH20) was used for assessments of molecular and phenotypical characteristics. The experiments were repeated twice.
4.5. Physiological Trait Measurements
Barley at the 5th leaf development stage was used for analysis. First, plants were photographed. Second, 5 cm fragments from the 3rd leaf were used for RWC analysis. Third, the 3rd and the 4th leaves were collected and stored at −80 °C for total RNA and protein extractions. Analyses were performed in five and three biological replicates, respectively. In six plants for each growing condition, shoot and root length were measured, and the fresh and dry biomasses of shoots and roots were weighed. Chlorophyll a fluorescence was measured using a fluorometer (FluorPen FP 100, ICT International, Australia) and leaf-clips with a window diameter of approximately 3 mm. Chlorophyll a parameters were measured according to manufacturer’s protocols with an adaptation to darkness of about 30 min.
4.6. Total RNA Extraction, cDNA Synthesis, and Real-Time PCR Analysis
Analyses were performed as described previously. Total RNA was extracted from the 3rd and the 4th leaves using TRI Reagent solutions (Invitrogen, Waltham, MA, USA). Genomic DNA was removed using DNase I, RNase-free (Thermo Fisher Scientific, Waltham, MA, USA). cDNA was synthetized using a Revert Aid cDNA Synthesis Kit (Thermo Fisher Scientific, Waltham, MA, USA). Real-Time PCR was carried out using the 5 x HOT FIREPol EvaGreen qPCR Mix Plus (noROX) (Solis BioDyne, Tartu, Estonia) kit and Rotor Gene 6000q Series (Corbett Life Science, Mortlake, Australia) thermalcycler according to the manufacturer’s protocols. The barley
ADP-rybosilation factor and
Glyceraldehyde-3-phosphate dehydrogenase (
GAPDH) genes were used as the internal controls. For each gene, three biological replicates were performed in three technical repeats, and the average value of the standard curve and standard error was shown. Gene-specific primers used for real-time PCR were published by Groszyk et al. [
11] and Groszyk and Szechyńska-Hebda [
45].
4.7. RNA-Seq Analysis
The mRNA sequencing service was outsourced to the commercial service laboratory of BGI Genomics Co., Ltd. (Hong Kong, China). Analyses were performed using Interdisciplinary Centre for Mathematical and Computational Modelling Warsaw University (Poland) and Galaxy software [
80]. Bioinformatics analysis was carried out by filtering out low-quality readings, containing unknown bases (N) and/or low-certainty readings (Q ≤ 20), and removing adapters using trimmomatic v. 0.39 [
81]. The next step was mapping the readings to the reference genome of
Hordeum vulgare (Hordeum_vulgare.IBSC_v2) [
82] with the HISAT2 v. 2.2.1 program [
83] and calculating the expression levels of the mapped genes using the StringTie v. 2.1.7 program [
84]. Expression normalization between the analyzed samples was performed using the DESeq2 tool [
85]. For the functional enrichment analysis of DEGs, we used the online Mercator4 v 5.0 tool [
86] and the PHANTER 7.0 database [
87].
4.8. Western Blot Detection of HvGSK2 and HvBZR1
Commercial anti-OsGSK2 (AbP80050-A-SE) and anti-OsBZR1 (AS16 3219) polyclonal antibodies were purchased from Beijing Protein Innovation Co., (Beijing, China) and Agrisera (Vännäs, Sweden), and used to detect HvGSK2.1 and HvBZR1, respectively. Total protein was extracted from the 3rd and the 4th leaves of barley. The plant material was ground to powder in liquid nitrogen and treated with a 1× SDS sample buffer (5 µL per 1 mg ground leaves). Supernatants were denatured at 70 °C for 10 min, centrifuged, and used for SDS-PAGE and immunoblot analysis. Anti-GSK2 and anti-OsBZR1 antibodies were used at dilutions of 1:5000 and 1:10000, respectively. Detection was performed using a PVDF membrane (Bio-Rad Laboratories, Hercules, CA, USA) and AgriseraECL SuperBright, AS16 ECL-S solutions (Agrisera, Vännäs, Sweden).
4.9. Data Analysis
Statistical analysis was performed using Microsoft Excel Professional Plus 2016 (Microsoft Office, Warszawa, Poland) and Statistica 13.0 (StatSoft, Kraków, Poland). Graphs were generated using Microsoft Excel Professional Plus 2016 and Microsoft PowerPoint Professional Plus 2016 (Microsoft Office, Warszawa, Poland).