Next Article in Journal
Obstacles against the Marketing of Curcumin as a Drug
Next Article in Special Issue
How Machine Learning Methods Helped Find Putative Rye Wax Genes Among GBS Data
Previous Article in Journal
Quaternary Ammonium Chitosans: The Importance of the Positive Fixed Charge of the Drug Delivery Systems
Previous Article in Special Issue
Combination of Linkage Mapping, GWAS, and GP to Dissect the Genetic Basis of Common Rust Resistance in Tropical Maize Germplasm
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Silencing of HvGSK1.1—A GSK3/SHAGGY-Like Kinase–Enhances Barley (Hordeum vulgare L.) Growth in Normal and in Salt Stress Conditions

by
Yuliya Kloc
1,
Marta Dmochowska-Boguta
1,
Andrzej Zielezinski
2,
Anna Nadolska-Orczyk
3,
Wojciech M. Karlowski
2 and
Waclaw Orczyk
1,*
1
Department of Genetic Engineering, Plant Breeding and Acclimatization, Institute–National Research Institute, Radzikow, 05-870 Blonie, Poland
2
Department of Computational Biology, Institute of Molecular Biology and Biotechnology, Faculty of Biology, Adam Mickiewicz University, Uniwersytetu Poznańskiego 6, 61-614 Poznan, Poland
3
Department of Functional Genomics, Plant Breeding and Acclimatization, Institute–National Research Institute, Radzikow, 05-870 Blonie, Poland
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2020, 21(18), 6616; https://doi.org/10.3390/ijms21186616
Submission received: 21 August 2020 / Revised: 4 September 2020 / Accepted: 8 September 2020 / Published: 10 September 2020
(This article belongs to the Special Issue Functional Genomics for Plant Breeding)

Abstract

:
Glycogen synthase kinase 3 (GSK3) is a highly conserved kinase present in all eukaryotes and functions as a key regulator of a wide range of physiological and developmental processes. The kinase, known in land plants as GSK3/SHAGGY-like kinase (GSK), is a key player in the brassinosteroid (BR) signaling pathway. The GSK genes, through the BRs, affect diverse developmental processes and modulate responses to environmental factors. In this work, we describe functional analysis of HvGSK1.1, which is one of the GSK3/SHAGGY-like orthologs in barley. The RNAi-mediated silencing of the target HvGSK1.1 gene was associated with modified expression of its paralogs HvGSK1.2, HvGSK2.1, HvGSK3.1, and HvGSK4.1 in plants grown in normal and in salt stress conditions. Low nucleotide similarity between the silencing fragment and barley GSK genes and the presence of BR-dependent transcription factors’ binding sites in promoter regions of barley and rice GSK genes imply an innate mechanism responsible for co-regulation of the genes. The results of the leaf inclination assay indicated that silencing of HvGSK1.1 and the changes of GSK paralogs enhanced the BR-dependent signaling in the plants. The strongest phenotype of transgenic lines with downregulated HvGSK1.1 and GSK paralogs had greater biomass of the seedlings grown in normal conditions and salt stress as well as elevated kernel weight of plants grown in normal conditions. Both traits showed a strong negative correlation with the transcript level of the target gene and the paralogs. The characteristics of barley lines with silenced expression of HvGSK1.1 are compatible with the expected phenotypes of plants with enhanced BR signaling. The results show that manipulation of the GSK-encoding genes provides data to explore their biological functions and confirm it as a feasible strategy to generate plants with improved agricultural traits.

1. Introduction

Glycogen synthase kinase 3 (GSK3) is a highly conserved kinase present in all eukaryotes and functioning as a key regulator of a wide range of physiological and developmental processes. The kinase, known in land plants as GSK3/SHAGGY-like kinase (GSK), is a key player in the brassinosteroid (BR) signaling pathway and, through the BRs, the GSKs affect diverse developmental processes and modulate responses to environmental factors [1]. The regulatory circuits connecting BR-dependent receptors at the cell surface with genetic effectors regulating more than 1000 genes have been defined in the last decade (for comprehensive reviews, see Zhu et al. [2], Belkhadir, and Jaillais [3]). Briefly, BR-dependent signaling is initiated by binding of the brassinolide (BL) molecule by BRI1/BRL1-BRL3 kinases. This step triggers the phosphorylation cascade of several components including BSK kinase and BSU1 phosphatase. The former specifically dephosphorylates and inactivates BIN2 kinase [4], which is a homologous protein showing high sequence identity to GSKs. The BIN2/GSK3 phosphorylates the two BR-dependent transcription factors (TF) BZR1 and BES1/BZR2 [2,5,6]. This consecutively inhibits their nuclear localization and DNA-binding activity [7,8] and, consequently, restrains BR-dependent gene regulation. BZR1 and BZR2 regulate expression of thousands of genes by interacting with other transcription regulators and binding to BR-responsive promoters [9,10]. Thus, BIN2/GSK3 functions as a negative regulator of BR signaling [6,11]. Initiation of the signaling cascade by BRs leads to inhibition of BIN2/GSK3 kinase, accumulation of TFs BZR1/BZR, and promotion of BR-dependent genes and traits.
Zhu et al. [12] added KIB1 as another component of the signaling cascade. This protein interacts in a BR-dependent mode with BIN2 and directs it to ubiquitination. Thus, KIB1 functions as an important positive regulator of BR signaling, which also acts under the BRs’ control. The activity of BIN2/GSK3, accumulation of BZR1, and BR-dependent growth stimulation are regulated by the Target of Rapamycin1 (TOR1). The TOR1 protein, functioning in the sugar signaling pathway, integrates nutrient and energy signaling with growth homeostasis. This step intertwines sugar-dependent and BR-dependent signaling and provides balanced regulation of steroid-promoted growth with carbon and energy availability [13,14]. GSK3/SHAGGY-like kinases (GSK) are ancient non-receptor Ser/Thr protein kinases found in all eukaryotes. Mammalian GSK3 exists in two isoforms, GSK3α and GSK3β, differing in the localization of their phosphate-binding pocket [15]. This domain in land plants’ GSKs is identical to mammalian GSK3β [16]. Additionally, in plants, the GSK proteins are encoded by multigene families.
There are 10 Arabidopsis GSK3/SHAGGY-like kinases (AtSK or ASK) assigned to four functional/phylogenetic groups [16]. They are involved in regulating diverse cellular and developmental processes including plant growth [6], flower [17], stomatal [18,19] and root development [20], seed size [21], and responses to environmental stresses [22,23]. At least 7 of the 10 Arabidopsis AtSKs (i.e., AtSK11, AtSK12, AtSK13 from group I, AtSK21, AtSK22, AtSK23 from group II and AtSK31 from group III) function in BR-dependent regulation, as it was shown in genetic screens and application of GSK3-specific inhibitor bikinin [24,25]. Since no data indicate similar roles for AtSK32, AtSK41, and AtSK42, their detailed functions remain to be elucidated. The structure of the ATP-binding pocket and sensitivity to bikinin indicate involvement in BR-dependent signaling. In contrast, insensitivity to bikinin, observed in AtSK42, argues against such a function [24]. The function related to carbohydrate metabolism and salt tolerance of group IV GSKs may be anticipated based on MsK4, which is the Medicago sativa ortholog of AtSK41. The gene was reported to regulate salt tolerance by adjusting carbohydrate metabolism in response to environmental stress [26,27].
Screening of genetic mutants and the results of exogenous application of BRs indicate that this group of phytohormones plays an important role in plant tolerance to abiotic stresses. BR application enhanced tolerance to freezing or chilling stress, a water deficit, and salt stress [28,29,30,31,32]. Treatments with BRs alleviated injuries caused by herbicide or insecticide treatments [33,34]. Zhu et al. [35] reported that enhanced salt tolerance observed in Nicotiana benthamiana, as a result of treatment with 0.1 µM brassinolide (BL), was induced via the nitric oxide signaling pathway. Exogenous BL application to rice improved plant growth and salinity tolerance by better osmoprotection, which is an increased level of proline and higher activities of antioxidant enzymes [36]. Compatible results, presented by Yusuf et al. [37], indicated that 24-epibrassinolide (EBL) treatment of wheat alleviated aluminum and salt stress injuries by modifying proline metabolism and improving Reactive Oxygen Species (ROS) scavenging. The ubiquitin proteasome system, known to participate in plant responses to environmental stresses was reported by Cui et al. [38] to be tightly connected with BR-dependent mechanisms of growth and salt tolerance. A number of reports on BR-dependent traits of cultivated plants have addressed the molecular mechanisms of BR-dependent regulation with the focus on GSK participation. Grain length locus GL3.3 in rice was annotated as OsGSK41 (OsGSK5). Plants with the gl3.3 loss-of-function mutation had increased cell length in spikelet hulls and elevated grain weight, which indicates the negative regulation of the trait by OsGSK41. Another mechanism of the effect of OsGSK41 on grain weight was reported by Hu et al. [39]. Investigating the qTGW3 locus and the qTGW3-encoded OsGSK41, they observed an interaction of this kinase with AUXIN RESPONSE FACTOR 4 and auxin-dependent signaling. The grain weight, although BR-dependent, requires tight interaction of BRs with auxin. Several GS loci (such as gs3 and GS9) further enhanced the trait by epistatic interaction with gl3.3 by producing extra-long grains [40]. The opposite mode of the OsGSK5 function was reported by Thitisaksakul et al. [41], where overexpression of the gene correlated with higher starch accumulation and better growth under severe salinity stress. As the authors highlighted, the results were compatible with another possible role of the GSKs assigned to group IV, namely, regulation of carbohydrate metabolism and the osmotic stress response.
Developmental processes and plant architecture pose another important target for modifications. The proper size of rice mesocotyls, which is a key trait facilitating rice cultivation, was found by Sun et al. [42] to be under strong selection during rice domestication. The authors revealed that OsGSK2 determine mesocotyl length variation and certain variants were selected as a result of rice domestication. BRs regulate leaf erectness, which is a trait with a potential strong impact on rice productivity. In this case, the signaling involves GSK3-like kinases and U-type cyclin [43]. Coordinated, organ-specific expression of cotton GSK3-like kinases and their involvement in the regulation of cotton fiber development and stress responses were reported by Wang et al. [44]. In barley, based on encoded amino acid sequences and protein kinase domains, seven transcriptionally active GSK-encoding genes were assigned to four groups, which is similar to AtSKs in Arabidopsis [45].
In this work, we describe the functional analysis of HvGSK1.1, which is one of the GSK3/SHAGGY-like orthologs in barley. RNAi-mediated silencing of the target HvGSK1.1 gene was associated with modified expression of its paralogs HvGSK1.2, HvGSK2.1, HvGSK3.1, and HvGSK4.1 in plants grown in normal conditions and in salt stress conditions. The expression changes of the target gene and its paralogs is discussed as the mutual co-regulation of genes belonging to one gene family. The conclusion is supported by the analysis of transcription factor binding sites identified in the promoter regions of the genes. The strongest phenotype of transgenic lines with downregulated HvGSK1.1 was greater biomass of the seedlings and greater weight of 1000 kernels. Both traits showed a strong negative correlation with the transcript level of the target gene and the selected paralogs. The characteristics of the plants are in agreement with the traits known to be regulated by BRs.

2. Results

The target HvGSK1.1 gene was annotated in our earlier study [45]. It was selected for the silencing experiments because of its similarity to rice OsGSK1, which is a knock-out mutation that enhanced plant tolerance to abiotic stresses [46]. The barley HvGSK1.1 has 12 exons and a similar exon-intron structure to the rice OsGSK1. The transcripts of OsGSK1 and HvGSK1.1 showed 73% nucleotide similarity. Kinase domains of the peptides encoded by the genes showed 84.2% identity and the identical pattern of the Lys-containing ATP-binding site and Asp-containing kinase active sites (Figure S1). The proteins encoded by the genes HvGSK1.1 (HORVU3Hr1G034440) and OsGSK1 (Os01g10840) showed a high level of conservation, e.g., very similar sequence length (406 and 403 aa residues, respectively), 84.2% identity, and 93% similarity. Very similar location of the protein kinase domain (Pfam AC: PF00069) (71–356 and 65–349 respectively) and the same pattern of ATP binding and kinase active site additionally confirm the close relation of both proteins (Figure S1).

2.1. Silencing Vector Construction and Agro-Based Transformation of Immature Barley Embryos

The 262 bp long open reading frame (ORF) fragment of HvGSK1.1 used for silencing had a relatively low nucleotide similarity with other barley GSK paralogs. The longest identical stretches were 17 nt, 16 nt, 14 nt, and 12 nt long with HvGSK1.2, HvGSK2.1, HvGSK3.1, and HvGSK4.1, respectively. In all cases, the identical stretches were much shorter than the 21 nt necessary for the efficient activation of cross-silencing.
Inoculation of 1975 scutella with Agrobacterium AGL1: pSoup and pBract207 with silencing cassettes and culture on hygromycin-containing media resulted in regeneration of 14 plants from independent explants. The presence of T-DNA with the silencing cassette was confirmed in 11 T0 plants. Transgene segregation in T1 seedlings indicated that T-DNA was integrated in a single locus in 10 lines and multiple loci in 1 line. Individual plants representing 10 lines with single transgene integration were analyzed for hemi-/homozygosity using qPCR and gDNA as a template. This approach allowed for fast selection of homozygous plants in the T1–T2 generations (Table 1).

2.2. Silencing of HvGSK1.1 Modifies Expression of GSK Paralogs in Plants Grown in Normal and Salt Stress Conditions

The seeds of the T4 generation of selected homozygous plants were used for all further molecular and phenotype analysis. The relative transcript level of the targeted HvGSK1.1 gene in plants grown in normal conditions was significantly lower at p < 0.001 in five lines and p < 0.01 in one line. Relative expression was 0.15 (line #1 and #6), 0.31 (line #5), 0.42 (line #8), 0.51 (line #9), 0.72 (line #10), and 0.82 (line #4) (Figure 1).
Expression of HvGSK1.1 in salt stress was lower when compared to plants in normal growth conditions. Relative expression in the non-transgenic control in salt stress was 0.47 of expression in normal conditions. In 5 lines #1 (0.06), #6 (0.07), #5 (0.22), #8 (0.19), and line #9 (0.20)-relative expression in salt stress conditions was significantly (p < 0.001) lower (Figure 1).
Silencing of the target HvGSK1.1 modified expression of four GSK paralogs HvGSK1.2, HvGSK2.1, HvGSK3.1, and HvGSK4.1. Relative expression of the paralogs varied depending on the line and the growth conditions. Relative expression of HvGSK1.2 in the five lines: #1 (0.33), #6 (0.29), #5 (0.68), #10 (0.63), and #4 (0.64) grown in normal conditions was significantly (p < 0.001 and p < 0.01) lower than in the non-transgenic control (Figure 2A). Expression in salt stress conditions in the five lines #1 (0.35), #6 (0.23), #8 (0.60), #9 (0.37), and # 4 (0.81) was significantly (p < 0.001 and p < 0.05) lower than in the non-transgenic control (Figure 2A). Expression of HvGSK2.1 was slightly downregulated in transgenic lines grown in normal conditions (Figure 2B). In salt stress, however, relative expression of the four lines #1 (0.62), #6 (0.69), #8 (0.57), and #9 (0.61) was significantly (p < 0.01 and p < 0.001) lower when compared to the non-transgenic control (Figure 2B).
Expression of HvGSK3.1 was significantly downregulated in most of the transgenic lines grown in either of the tested conditions. In normal conditions, expression in lines #1 (0.64), #6 (0.48), #5 (0.54), #10 (0.62), and #4 (0.45) was significantly (p < 0.001 and p < 0.01) lower when compared with the non-transgenic control (Figure 2C). In salt stress conditions, expression was downregulated in three lines: #8 (0.37), #9 (0.44), and #10 (0.55) (Figure 2C). Expression of HvGSK4.1 in normal conditions was significantly (p < 0.05, p < 0.01 and p < 0.001) downregulated in four transgenic lines #6 (0.71), #5 (0.6), #10 (0.44), and #4 (0.6) (Figure 2D). In salt stress conditions, the three lines #8 (0.56), #9 (0.74), and #4 (0.81) were significantly downregulated (Figure 2D). Northern blot hybridization confirmed the presence of the siRNA molecules complementary to the ORF fragment present in the silencing cassette. The small RNAs were detected in transgenic plants of selected lines in the T2 generation. They were not detectable in the non-transgenic control (Figure S2).

2.3. Leaf Inclination Biotest

The leaf inclination biotest was developed as a semi-quantitative assay of BRs [47] and BR-dependent signaling [48]. The test was adapted for barley and showed that higher concentrations of EBL and of bikinin were associated with bigger lamina joint angles (Figure S3A,B). The inclination angle in the control was 36.8° (SD 3°). After treatment with EBL (0.01 µM), it was 40.7° (SD 4.9°), and, after bikinin (5 µM) treatment, it was 37° (SD 7.6°). After simultaneous treatment with both compounds EBL (0.01 µM) and bikinin (5 µM), the angle was 59.2° (SD 14°) (Figure S3C). The concentrations 0.01 µM EBL and 5µM bikinin were used to test the response of the non-transgenic control (Figure 3A) and transgenic lines #1, #8, and #4, which represented strong, medium, and weak silencing of HvGSK1.1 (Figure 3B–D, Figure S4). The inclination angles in water-incubated line #1 was 61°. After 0.01 µM EBL, the angle was 74.7°. After 5 µM bikinin, the angle was 70° and the combined treatment with both compounds increased the angle to 76.7° (Figure 3B, Figure S4). The inclination angle in water-incubated line #8 was 39.3°. After 0.01 µM EBL, the angle was 42°. After 5 µM bikinin, the angle was 54° and the combined treatment with both compounds gave a 36° angle (Figure 3C, Figure S4). The inclination angles in the water-incubated line #4 was 42°. After 0.01 µM EBL, the angle was 37.3°. After 5 µM bikinin, the angle was 43.3° and combined treatment with both compounds elevated the angle to 44.5° (Figure 3D, Figure S4). The inclination angles of line #1 in the reaction for either a single or both compounds were greater than in the non-transgenic control. The reaction of lines #4 and #8 with the weaker silencing of HvGSK1.1 was similar to the non-transgenic control.

2.4. Silencing of HvGSK1.1 Increases Biomass of Plants Grown in Normal and Salt Stress Conditions and Kernel Weight

Relative biomass of all lines was greater than biomass of the control from the same conditions and, for six lines, the values were significantly different (Figure 4 and Figure 5). Lines #1 (1.21), #6 (1.22), and #8 (1.25) exceeded the biomass of the control by 21% to 25% with high statistical significance p < 0.001. Biomass of lines #5 (1.21) and #10 (1.18) was greater than the control with p < 0.01 and, for line #9 (1.15), the significance was p < 0.05 (Figure 4).
Biomass of the plants grown in salt stress was smaller and, for the non-transgenic control, it was 0.46 of the biomass in normal conditions. However, relative biomass of all lines #1 (0.62), #6 (0.57), #5 (0.55), #8 (0.57), #9 (0.55), #10 (0.54), and #4 (0.51) was greater than the control from the same conditions. The values shown by three lines #1, #6, and #8 exceeded the control by 24% to 35% with significance at p < 0.01 and p < 0.001. The values of another lines #5, #9, and #10 exceeded the control by 17% to 20% with a significance at p < 0.05 (Figure 4).
Thousand kernel weight (TKW) of the control plants grown in normal conditions was 31 g (Figure 6). In transgenic lines, the values were: #1 (43), #6 (47.9), #5 (37.4), #8 (35.8), #9 (31.6), #10 (35.3), and #4 (32). The TKW values of #1 and #6 were significantly (p < 0.05) higher than those of the control (Figure 6 and Figure 7).

2.5. Silencing of HvGSK1.1 Correlates with Expression of GSK Paralogs, Plant Biomass, and Thousand Kernel Weight

Expression of HvGSK1.1 in plants grown in normal conditions correlated with expression of three paralogs: HvGSK1.2 (r = 0.59), HvGSK3.1 (r = 0.71), and HvGSK4.1 (r = 0.49). In salt stress, it correlated with HvGSK1.2 (r = 0.77) and HvGSK2.1 (r = 0.63) (Figure 8). Expression of the four paralogs showed a net of mutual positive correlations, which were different in normal and salt stress conditions (Figure S5). The most clear-cut characteristics of all transgenic lines had greater biomass of the seedlings grown in normal and in salt stress conditions (Figure 4 and Figure 5). Statistical analysis revealed a strong negative correlation between the biomass and the expression of HvGSK1.1 in plants grown in normal conditions (r = −0.77) and in salt stress (r = −0.83) (Figure 8). The biomass was also negatively correlated with the expression of HvGSK4.1 (r = −0.60) in normal conditions and with expression of HvGSK1.2 (r = −0.69) and HvGSK2.1 (r = −0.70) in salt stress conditions (Figure 8). This was consistent with the positive mutual correlations between the paralogs’ expressions (Figure S5). Kernel weight, checked in plants in normal conditions and shown as thousand kernel weight (TKW), showed a strong negative correlation with HvGSK1.1 (r = −0.82), HvGSK1.2 (r = −0.87), and HvGSK3.1 (r = −0.88) (Figure 8, Figure S5).

3. Discussion

The HvGSK1.1 gene was the first barley GSK-encoding gene identified by us. It was selected for RNAi-based silencing because of its high sequence similarity to rice OsGSK1, whose knockout mutation conferred elevated tolerance to environmental stresses [46]. The ORF fragment of HvGSK1.1 selected for construction of the silencing cassette had a low-level similarity to the cDNA sequences of the GSK paralogs. The longest identical stretches between the silencing cassette and the paralogs, from 12 to 17 nt, excluded efficient activation of cross-silencing. The Gateway strategy proved to be very useful for constructing pBract207-based silencing cassette using the fragment of the target gene. As a result of Agro-based transformation, the T-DNA was integrated as a single locus in 10 plants. Only a single plant contained the T-DNA integrated into more than one locus. The 1:10 rate of multiple vs. single integration (Table 1) is lower when compared with the average 46% rate of multiple integrations reported by Bartlett et al. [49] in a large group of transgenic barley plants. The relatively high rate of lines with single integration is an advantage considering the expected higher stability of the transgene expression in the subsequent generations. As found by Zalewski et al. [50], this is particularly important when the targeted gene is developmentally regulated and it is involved in hormone regulation. Results of the Northern blot (Figure S2) confirmed that short (22–25 nt) RNA molecules, complementary to the silencing fragment of the target gene, were present in the selected transgenic lines and they were undetectable in the non-transgenic control plants. The presence of this class of short RNA molecules, complementary to the silencing cassette, confirms the synthesis of the hpRNA transcript and its processing into siRNAs. The analysis of target gene expression in T2 plants showed that the silencing was strong and ranged from 0.1 to 0.6 in all transgenic lines Figure S6). Seed multiplication necessary for all experiments required two additional generations of self-pollinated plants. In T4 generation plants, the target gene silencing remained in seven lines. These lines represented a gradient of HvGSK1.1 silencing ranging from 0.15 to 0.82 (Figure 1) and a set of molecular and phenotype changes.
The HvGSK1.1 silencing was associated with modified expression of the GSK-encoding paralogs. The original silencing of the target gene associated with this modified expression of the paralogs showed a net of mutually correlated changes present in both growth conditions (Figure 8, Figure S5). A low level of nucleotide similarity between the silencing cassette and the paralogs excluded cross silencing as the process responsible for the co-regulation. We propose that such an extensive net of correlations, which was further strengthened by the salt stress, implies an innate mechanism of GSK gene regulation. Although, at this point, there are no experimental data to explain it, we think that the pattern of transcription factor binding sites (TFBS) found in promoter regions of the genes indicates such a mechanism. The promoter regions analysis of several barley GSKs showed a presence in multiple binding sites of the two BR-dependent transcription factors BZR1 and BZR2 (Table S1). The two BR-dependent TFs are the targets of the GSK-dependent phosphorylation and inactivation [2,51]. This combination of positive transcriptional regulation and a negative post-transcriptional regulation indicates a possible feedback mechanism coordinating expression of the GSK genes. We think that the net of mutual co-regulation could be the result of such a mechanism. This notion is further confirmed by (i) the presence of the BZR1/BZR2 binding sites in promoter regions of OsSK11 and OsGSK1, the two rice orthologs of HvGSK1.1 (Table S1), and (ii) a large number of common TFBS present in the promoter regions of the selected barley, rice, and Arabidopsis GSK genes (Figure S7, Table S2). These observations are in agreement with extensive and strong co-regulation of the GSK genes found in generative organs in barley [45]. They are also compatible with the results of chromatin immuno-precipitation showing that several GSK-encoding (AtSK12, AtSK13, AtSK21, and AtSK41) and BZR2-encoding genes are among the specific targets of the BZR1 transcription factor [9]. The presence of BR-specific TFBS in at least four GSK-encoding genes in Arabidopsis indicates a complex network of co-regulation of the GSK paralogs. Our results showing coordinated regulation of the GSKs in barley provide novel data on the mechanisms of BR homeostasis in plants. A similar data, concerning RNAi-based silencing of TaCKX1 and associated with this co-regulation of CKX paralogs, were reported by Jablonski et al. [52]. The authors found that diverse levels of target gene silencing were associated with different models of co-expression with other TaCKX genes and parameters of yield-related traits in wheat.
Bikinin was identified as a GSK3 inhibiting molecule. In Arabidopsis, 7 out of the 10 ASKs were explicitly inhibited by this compound when enhancing BR-dependent signaling [24,53]. We hypothesized that HvGSK1.1 silencing would enhance the BR-dependent response and mimic the effect of the bikinin treatment. To verify this, we adapted the procedure of leaf inclination for barley based on the original test developed for rice [47] and confirmed by the findings of Yamamuro et al. [54] that leaf inclination angle and internode elongation were specifically regulated by BRs. The test adapted for barley showed that higher concentrations of EBL were associated with greater angles of the lamina joints. Taking into account the function, bikinin could be considered as the physiological equivalent of the silencing of the GSK-encoding genes. This effect was clearly visible in seedling fragments treated with the tested concentrations of EBL and bikinin (Figure S3A,B). The simultaneous treatment with EBL and bikinin showed a stronger effect than the impact of either of the compounds applied separately (Figure S3C). The changes observed after EBL and bikinin treatment in leaves of the non-transgenic control and the three lines #1, #8, and #4 indicated that silencing of HvGSK1.1 enhanced the reaction for exogenous BL in a similar way as the bikinin treatment. Based on these results, we concluded that silencing of HvGSK1.1 enhanced the BR-dependent signaling in the plants. The final effect was the result of the experimental silencing of HvGSK1.1 and correlated with this down-regulation of other GSK paralogs. It agrees with findings that bikinin inhibiting the GSKs acts as an activator of BR signaling and mimics the exogenous application of BRs [24]. We propose that both the original silencing of HvGSK1.1 correlated with this down-regulation of GSK paralogs enhanced the BR-dependent signaling and, in that way, promoted certain traits regulated by BRs.
The most notable feature of all transgenic lines found in normal and in salt stress conditions was their greater biomass compared with the non-transgenic control. The trait showed a strong negative correlation with the expression of the HvGSK1.1 and the GSK paralogs. The observation agrees with the reported greater weight of Arabidopsis plants grown in normal and in salt stress growth conditions after EBL treatment [28]. Similar results of growth promotion and alleviation of salt stress have been presented in several other studies [36,55,56], and this may involve auxin-dependent regulation [57]. Another trait associated with the experimental silencing of HvGSK1.1 and correlated with this down-regulation of the GSKs was increased weight of the kernels. The trait, shown as TKW, was significantly bigger in lines #1 and #6 with the strongest HvGSK1.1 silencing (Figure 6 and Figure 7). It showed a strong negative correlation with HvGSK1.1 (r = −0.82), HvGSK1.2 (r = −0.87), and HvGSK3.1 (r = −0.88) (Figure 8, Figure S5). This set of results fully agrees with articles showing that productivity-related traits are stimulated by BRs and BR-dependent signaling. These include significantly increased seed number in a BR-enhanced mutant [58], smaller seeds in BR-deficient mutants, and confirmation that seed size, mass, and shape in Arabidopsis were regulated by BRs [59]. It agrees with a finding by Che et al. [21] that mutation of grain lenths2 loci in rice activated BR-dependent responses and, as a consequence, increased the weight of rice grains by 27%. The former authors concluded that modulating specific brassinosteroid responses could improve plant productivity. Wu et al. reported compatible results of significantly increased grain yield and weight in rice with enhanced production of BRs [60]. Enhanced barley growth observed after silencing of HvGSK1.1 is also compatible with results of Wan et al. [61]. They reported that PeGSK1 acts as a negative regulator of cell growth in Moso bamboo, which is one of the fastest growing monocots.
The characteristics of barley lines with silenced expression of HvGSK1.1 are compatible with the expected phenotypes of plants with enhanced BR signaling. The results show that manipulation of the GSK3-encoding genes provides data to explore their biological functions and, at the same time, indicate it as a feasible strategy to generate plants with improved agricultural traits.

4. Materials and Methods

4.1. Bioinformatics of HvGSK1.1 and Rice Orthologs

The pairwise alignment between protein sequences of HvGSK1.1 (HORVU3Hr1G034440) and OsGSK1 (Os01g10840) was performed in the needle program from the EMBOSS package [PMID: 30976793]. The identification of protein domains in HvGSK1.1 and OsGSK1 was performed in InterProScan [PMID: 30398656]. The prediction of transcription factors’ binding sites (TFBS) in the promoter regions (up to 4000 nucleotides) of GSK-encoding genes in barley, rice, and Arabidopsis was carried out by the Match™ program [PMID: 12824369] and the TRANSFAC® database (Release 2020.1) [PMID: 18436575]. Match™ was set to use only high-quality TFBS matrices to minimize false positives. The predicted TFBS were also confirmed by the CiiiDER tool (with default parameters) [PMID: 31483836] and non-redundant transcription factor plant models from the JASPAR database (Release JASPAR2020) [62].

4.2. Plant Material and Growth in Normal and Salt Stress Conditions

The barley (Hordeum vulgare L.) cultivar Golden Promise was used as a source of plant material in all experiments. Barley kernels, after 48 h of imbibition, were planted in pots (14-cm diameter) filled with peat substrate mixed with sand in a 10:3 v/v ratio. The seedlings were cultivated in a growth chamber with a 16-h photoperiod at 21 °C in the day and 18 °C in the night. The relative humidity was in the range of 60–80%, and the light intensity was 130–230 μmol photons m−2 s−1. Plants were irrigated twice a week and fertilized once a week with multi-component soil fertilizer Florovit (http://florovit.pl/), according to the manufacturer’s instructions. All experiments of transcript quantification and phenotype characterization were done using homogenous seed material of T4 generation obtained after self-pollination of selected T2 homozygotes grown in the same conditions as the donor plants. Germinating seeds of the tested lines and the non-transgenic control, after 48 h imbibition in Hoagland solution [63] or Hoagland supplemented with 25 mM NaCl, were placed on filter paper, which was tightly rolled, saturated with Hoagland solution, or Hoagland supplemented with 200 mM NaCl and placed in a growth chamber in the same conditions as the donor plants. Cultivation in Hoagland or Hoagland supplemented with NaCl is referred as growth in normal or in salt stress conditions, respectively. The 14-day old seedlings grown in normal and in salt stress conditions were used for biomass quantification and for leaf collection, RNA extraction, and transcript quantification. The plants cultivated in soil in the same conditions as the donor plants were used for kernel characterization.

4.3. Construction of Silencing Vector and Agrobacterium-Mediated Transformation of Barley

A barley clone representing a cDNA of HvGSK1.1 transcript (HORVU3Hr1G034440.2) was amplified using AmpliTaq Gold DNA Polymerase (Thermo-Fisher Scientific, Waltham, MA, USA), the primers 1FRGSK1 and 4RevGSK1 (Table 2), and cDNA of barley cv. Golden Promise as a template. The amplified fragment was cloned to pCR2.1TOPO (Invitrogen by Thermo Fisher Scientific, Waltham, MA, USA), sequenced (Genomed, Warsaw, Poland), and assigned as HvGSK1.1 (Figure 9A).
Gateway recombination was performed in 5 μL of a reaction mixture containing entry and destination molecules in a 2:1 molar ratio, total DNA 75–100 ng, and 1 μL of LR-Clonase 5x enzyme mix (Life Technologies by Thermo Fisher Scientific, Waltham, MA, USA). The 1 μL of 5x water-diluted reaction mix was electroporated (2.5 kV, 25 μF, 200 Ω, Gene Pulser II, Bio-Rad, Hercules, CA, USA) into 50 μL of electro-competent E. coli OmniMAX2-T1 strain (Life Technologies) and cultured on LB with kanamycin (50 mg/L). Silencing cassettes were confirmed by sequencing using primers anchored in the Ubi promoter, intron, and nos terminator of pBRACT207 (Figure 9B). The plasmid with an expected sequence and structure of the silencing cassette was electroporated into Agrobacterium tumefaciens AGL1 strain containing a pSoup helper plasmid [65]. Immature embryos, isolated from barley cv. Golden Promise, were inoculated with Agrobacterium AGL1 strain containing the pBRACT207-based construct and were in vitro cultured, according to the modified procedure [49,66]. Regenerated plants were transferred to soil and grown in conditions identical to barley donor plants. In vitro regenerated, non-transgenic plants were used as controls for all phenotypic observations and all molecular analyses.

4.4. Nucleic Acid Isolation, Transcript Quantification Leaf Inclination Test

Genomic DNA (gDNA) was isolated from plant samples, according to the modified CTAB procedure [67]. Total RNA was extracted from leaves of 14-day-old plants grown in normal conditions (Hoagland solution) and in salt stress conditions (Hoagland with 200 mM NaCl) using RNA 3-zone reagent solution (Novazym, Poznan, Poland). RNA concentration and A260/280 ratio (always higher than 1.8) were measured using a NanoDrop spectrophotometer (NanoDrop Technologies, Wilmington, DE, USA). The RNA quality was further determined in agarose gel. Isolated total RNA was treated with 2 U of DNase RNase-free (Roche Diagnostics, Mannheim, Germany) and 2 U of Protector RNase inhibitor followed by DNase inactivation, according to the manufacturer’s protocol. Genomic DNA impurities were checked by polymerase chain reaction (PCR) with primers qAct1 and qAct2a (Table 2) specific to the β-actin encoding gene, and 100 ng of DNase treated RNA as a template. Complete removal of gDNA was confirmed by a lack of detectable amplicon after 36 cycles of amplification. Two micrograms of RNA (non-degraded, DNase-treated with undetectable gDNA impurities) were used as a template for the reverse transcription reaction with oligo d(T) primers and RevertAid First Strand cDNA Synthesis Kit (ThermoFisher Scientific, USA). Obtained cDNA was diluted five-fold and used directly as a template for quantitative PCR (qPCR). The standard qPCR reaction mix was composed of 2.2 µL of 5x HOT FIREPol EvaGreen qPCR Mix Plus (ROX) (Solis Biodyne, Tartu, Estonia), 0.25 µL of primer F (10 µM), 0.25 µL of primer R (10 µM), 1 µL of cDNA, and water to 11 µL. The reaction was performed in a Rotor-Gene 6000 model 5-plex thermocycler (Corbett Life Science, Chadstone, Australia). The efficiencies of amplification for all primer pairs were in the range of 0.85 to 1.00 and R2 0.99875–1.00000. The specificity of amplification was verified by melting curve analysis. Template concentrations ranging from 104 to 108 copies of analyzed amplicon and 104–108 copies of the reference gene ADP-ribosylation factor (ADP RF) AJ508228 per reaction were used as the standards for qPCR. Primers for transcript quantification were designed based on the sequences deposited in NCBI and Ensemble Plants (Table 2). Threshold line, Ct values, standard curves, and relative quantifications were determined using the proprietary Rotor-Gene 6000 software v 1.7 (Corbett Life Science, Chadstone, Australia). The transcripts’ quantification data represent a medium of at least three biological replicates with three technical repetitions each. The whole procedure of RNA isolation, reverse transcription, qPCR conditions, and data analysis met the Minimum Information for Publication of Quantitative Real-Time PCR Experiments (MIQE) criteria outlined by Bustin et al. [68].
Identification of transgene hemi-zygotes and homo-zygotes in the T1 plants was made using a qPCR reaction with primers designed for hptII gene, primers designed for β-actin-encoding gene, which served as a reference, and using barley genomic DNA (gDNA) isolated from the tested T1 plants as a template. The reaction was carried out in a 11 μL mixture containing 2.2 µL of HOT EvaGreen qPCR Mix Plus and 0.25 µL (10 µM) of each primer for hptII (pBr_Hyg_F2/pBr_Hyg_R2) or β-actin (qAct1/qAct2a) amplification (Table 2), 7.3 µL of H2O, and 1 µL containing 50 ng of gDNA as a template. Each reaction was run in triplicate. The homozygotes were discriminated from hemizygotes based on the values of hptII to actin ratio calculated for at least six plants representing the T1 progeny of each transgenic line.
The fraction of small RNA was obtained from total RNA after precipitating high molecular weight RNA using 10% PEG (MW 8000) and 1 M NaCl and collecting small RNAs using a standard isopropanol method [69]. 10 µg of samples of fractionated RNA were separated by PAGE (15% acrylamide:bis acrylamide 19:1 v/v) in TBE denaturing buffer, transferred to nylon membrane positively charged using a semidry electroblotter (Bio-Rad, Hercules, CA, USA), and immobilized by UV cross-linking. The ORF-derived silencing fragment was used as a template for probe labeling using the PCR DIG Probe synthesis kit. Hybridization and chemiluminescent detection were done using a DIG detection kit. Labeling, hybridization, and detection were done, according to the manufacturer’s protocols (Roche Diagnostics, Mannheim, Germany).
The leaf inclination biotest was adapted from the protocol of Wada et al. [47]. The excised upper part of the seedlings with the first leaf developed and the second leaf emerging were incubated for 72 h in the tested concentrations of 24-epibrassinolide (EBL), (OlChemIm, Olomouc, Czech Republic) and bikinin (Sigma-Aldrich, Steinheim, Germany, a kind gift from prof. R. Malinowski). After incubation, the leaf fragments were used to assess leaf inclination angles. At least five plants of lines #1, #8, and #4 representing strong, medium, and weak HvGSK1.1 silencing, respectively, were used. At least six seedlings representing transgenic lines and the non-transgenic control were used to assess the biomass of the plants grown in normal conditions and salt stress conditions. The kernels collected from at least three plants of transgenic lines and non-transgenic control grown in soil in the same conditions as donor plants were used to calculate the thousand kernel weight (TKW).

4.5. Statistical Analyses

The Pearson correlation coefficient was calculated (Statistica version 13, Tibco, Palo Alto, CA, USA) between the values of the relative expressions of the analyzed genes, the biomass of the seedlings, and the values of thousand kernel weight (TKW). The correlation coefficient was significant at p < 0.05. The values between 0 and ǀ0.3ǀ were defined as a weak (low) correlation, the values between ǀ0.3ǀ and ǀ0.7ǀ were defined as a medium (moderate), and the values between ǀ0.7ǀ and ǀ1ǀ were defined as a strong (high) correlation. Data are presented as the mean value with the standard deviation. Quantitative RT-PCR data represented the medium of at least three independent biological replications with each done using the three technical repetitions. Standard deviation represents the variation between the biological repetitions. The statistical significance of the differences was analyzed using the ANOVA test, which was followed by the LSD post-hoc test (STATISTICA 10, StatSoft). The differences were considered statistically significant at p < 0.05.

Supplementary Materials

Supplementary Materials can be found at https://www.mdpi.com/1422-0067/21/18/6616/s1. Figure S1: Protein domain organization of GSK orthologs encoded by: HvGSK1.1 (HORVU3Hr1G034440.2) and OsGSK1 (Os01g10840). The amino acid defining features of the kinase domains i.e. protein kinase domain (Pfam AC: PF00069), ATP binding site and kinase active site are indicated. Figure S2: Northern blot detection of low molecular weight RNA in barley leaves in T2 generation using DIG labeled probe complementary to HvGSK1.1 ORF fragment in the silencing cassette. Indicated are U6 and 25 nt marker (lane 1), RNA from non-transgenic control plants (lanes 5, 9), RNA from lines #1, #2, #4 and #6. Lanes 2, 3, 4 and 6 were not loaded with the RNA to avoid cross contamination of neighboring lanes with the RNA samples. Figure S3: Representative pictures of leaf fragments from wild type barley cv. Golden Promise treated with indicated concentrations of 24-epibrassinolide (EBL) (A), bikinin (B). Pictures of leaf fragments treated with threshold concentrations of 0.01μM EBL and 5μM bikinin applied separately and together (C). Indicated are medium inclination angles and standard deviation (SD) of at least five plants. Bar represents 5 cm. Figure S4: Leaf inclination angles of nontransgenic control and transgenic lines #1, #8 and #4, which represented strong, medium and weak silencing of HvGSK1.1. Leaf fragments were treated with 24-epibrassinolide (EBL) 0.01 μM, bikinin 5 μM and EBL (0.01 μM) and bikinin (5 μM) applied together. Data represent medium inclination angles and standard deviation (SD) of at least five plants. Figure S5: Graphical representation of Pearson correlation net between expression of HvGSK1.1 and HvGSK paralogs, and the values of thousand kernel weight (TKW) measured in the non-transgenic control and the transgenic lines grown in normal conditions (A) and in salt stress conditions (Hoagland with NaCl 200 mM) (B). Positive and negative correlations are indicated with red and blue respectively. Bar width denotes the r values according to the scale. Figure S6: Relative expression of HvGSK1.1 in leaves of T2 plants grown in normal (Hoagland medium) and in salt stress conditions (Hoagland medium supplemented with NaCl 200mM). Relative expression of HvGSK1.1 in leaves of non-transgenic control plants grown in normal conditions was assumed as 1.0. Figure S7: The Upset plot shows the number of shared transcription factor binding sites (TFBS) in the promoter regions of GSK-encoding genes from barley (HvGSK1.1, HvGSK1.2, HvGSK1.3, HvGSK2.1, HvGSK3.1, HvGSK4.1) and rice (OsGSK1 and OsGSK1.1). Promoter region of HvGSK1.1 showed the highest number of shared TFBS with (n = 280) and OsGSK1 (n = 249). Table S1: The number of binding sites for brassinosteroid-dependent transcription factors BZR1 and BZR2 identified in the GSK-encoding genes based on TRANSFAC. Table S2: The list of all transcription factors binding sites (TFBS) in promoter regions of GSK-encoding genes in barley, rice and Arabidopsis.

Author Contributions

Y.K.—carried the experiments, vector construction, barley genetic transformation, selection and molecular analyses of transgenic lines, phenotyping, and expression profiling. M.D.-B.—phenotyping and statistical analyses. A.Z.—computational genomic and phylogenetic analyses. A.N.-O.—contribution to conception and discussion of the project and participation in manuscript writing. W.M.K.—concept and supervision of the computational part of the study, interpretation of the result, and participation in manuscript writing. W.O.—concept of the project, analysis and interpretation of the results, and writing the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

The National Science Center grant No: UMO-2014/13/B/NZ9/02381 (WO) funded the research.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study, in the collection, analyses, or interpretation of data, in the writing of the manuscript, or in the decision to publish the results.

References

  1. Saini, S.; Sharma, I.; Pati, P.K. Versatile roles of brassinosteroid in plants in the context of its homoeostasis, signaling and crosstalks. Front. Plant Sci. 2015, 6, 950. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Zhu, J.Y.; Sae-Seaw, J.; Wang, Z.Y. Brassinosteroid signalling. Development 2013, 140, 1615–1620. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Belkhadir, Y.; Jaillais, Y. The molecular circuitry of brassinosteroid signaling. New Phytol. 2015, 206, 522–540. [Google Scholar] [CrossRef]
  4. Kim, T.W.; Guan, S.; Sun, Y.; Deng, Z.; Tang, W.; Shang, J.X.; Sun, Y.; Burlingame, A.L.; Wang, Z.Y. Brassinosteroid signal transduction from cell-surface receptor kinases to nuclear transcription factors. Nat. Cell Biol. 2009, 11, 1254–1260. [Google Scholar] [CrossRef] [PubMed]
  5. Wang, Z.Y.; Nakano, T.; Gendron, J.; He, J.X.; Chen, M.; Vafeados, D.; Yang, Y.L.; Fujioka, S.; Yoshida, S.; Asami, T.; et al. Nuclear-localized BZR1 mediates brassinosteroid-induced growth and feedback suppression of brassinosteroid biosynthesis. Dev. Cell 2002, 2, 505–513. [Google Scholar] [CrossRef] [Green Version]
  6. Yin, Y.H.; Wang, Z.Y.; Mora-Garcia, S.; Li, J.M.; Yoshida, S.; Asami, T.; Chory, J. BES1 accumulates in the nucleus in response to brassinosteroids to regulate gene expression and promote stem elongation. Cell 2002, 109, 181–191. [Google Scholar] [CrossRef] [Green Version]
  7. Gampala, S.S.; Kim, T.W.; He, J.X.; Tang, W.Q.; Deng, Z.P.; Bai, M.Y.; Guan, S.H.; Lalonde, S.; Sun, Y.; Gendron, J.M.; et al. An essential role for 14-3-3 proteins in brassinosteroid signal transduction in Arabidopsis. Dev. Cell 2007, 13, 177–189. [Google Scholar] [CrossRef] [Green Version]
  8. Vert, G.; Chory, J. Downstream nuclear events in brassinosteroid signalling. Nature 2006, 441, 96–100. [Google Scholar] [CrossRef]
  9. Sun, Y.; Fan, X.Y.; Cao, D.M.; Tang, W.Q.; He, K.; Zhu, J.Y.; He, J.X.; Bai, M.Y.; Zhu, S.W.; Oh, E.; et al. Integration of Brassinosteroid Signal Transduction with the Transcription Network for Plant Growth Regulation in Arabidopsis. Dev. Cell 2010, 19, 765–777. [Google Scholar] [CrossRef] [Green Version]
  10. Oh, E.; Zhu, J.Y.; Bai, M.Y.; Arenhart, R.A.; Sun, Y.; Wang, Z.Y. Cell elongation is regulated through a central circuit of interacting transcription factors in the Arabidopsis hypocotyl. Elife 2014, 3, e03031. [Google Scholar] [CrossRef]
  11. He, J.X.; Gendron, J.M.; Yang, Y.L.; Li, J.M.; Wang, Z.Y. The GSK3-like kinase BIN2 phosphorylates and destabilizes BZR1, a positive regulator of the brassinosteroid signaling pathway in Arabidopsis. Proc. Natl. Acad. Sci. USA 2002, 99, 10185–10190. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Zhu, J.Y.; Li, Y.; Cao, D.M.; Yang, H.; Oh, E.; Bi, Y.; Zhu, S.; Wang, Z.Y. The F-box Protein KIB1 Mediates Brassinosteroid-Induced Inactivation and Degradation of GSK3-like Kinases in Arabidopsis. Mol. Cell 2017, 66, 648–657 e4. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Zhang, Z.Z.; Zhu, J.Y.; Roh, J.; Marchive, C.; Kim, S.K.; Meyer, C.; Sun, Y.; Wang, W.F.; Wang, Z.Y. TOR Signaling Promotes Accumulation of BZR1 to Balance Growth with Carbon Availability in Arabidopsis. Curr. Biol. 2016, 26, 1854–1860. [Google Scholar] [CrossRef] [PubMed]
  14. Xiong, F.J.; Zhang, R.; Meng, Z.G.; Deng, K.X.; Que, Y.M.; Zhuo, F.P.; Feng, L.; Guo, S.D.; Datla, R.; Ren, M.Z. Brassinosteriod Insensitive 2 (BIN2) acts as a downstream effector of the Target of Rapamycin (TOR) signaling pathway to regulate photoautotrophic growth in Arabidopsis. New Phytol. 2017, 213, 233–249. [Google Scholar] [CrossRef] [Green Version]
  15. Kim, L.; Kimmel, A.R. GSK3 at the edge: Regulation of developmental specification and cell polarization. Curr. Drug Targets 2006, 7, 1411–1419. [Google Scholar] [CrossRef]
  16. Jonak, C.; Hirt, H. Glycogen synthase kinase 3/SHAGGY-like kinases in plants: An emerging family with novel functions. Trends Plant Sci. 2002, 7, 457–461. [Google Scholar] [CrossRef]
  17. Dornelas, M.C.; van Lammeren, A.A.M.; Kreis, M. Arabidopsis thaliana SHAGGY-related protein kinases (AtSK11 and 12) function in perianth and gynoecium development. Plant J. 2000, 21, 419–429. [Google Scholar] [CrossRef]
  18. Kim, T.W.; Michniewicz, M.; Bergmann, D.C.; Wang, Z.Y. Brassinosteroid regulates stomatal development by GSK3-mediated inhibition of a MAPK pathway. Nature 2012, 482, 419–422. [Google Scholar] [CrossRef] [Green Version]
  19. Khan, M.; Rozhon, W.; Bigeard, J.; Pflieger, D.; Husar, S.; Pitzschke, A.; Teige, M.; Jonak, C.; Hirt, H.; Poppenberger, B. Brassinosteroid-regulated GSK3/Shaggy-like Kinases Phosphorylate Mitogen-activated Protein (MAP) Kinase Kinases, Which Control Stomata Development in Arabidopsis thaliana. J. Biol. Chem. 2013, 288, 7519–7527. [Google Scholar] [CrossRef] [Green Version]
  20. Cho, H.; Ryu, H.; Rho, S.; Hill, K.; Smith, S.; Audenaert, D.; Park, J.; Han, S.; Beeckman, T.; Bennett, M.J.; et al. A secreted peptide acts on BIN2-mediated phosphorylation of ARFs to potentiate auxin response during lateral root development. Nat. Cell Biol. 2014, 16, 66–76. [Google Scholar] [CrossRef]
  21. Che, R.; Tong, H.; Shi, B.; Liu, Y.; Fang, S.; Liu, D.; Xiao, Y.; Hu, B.; Liu, L.; Wang, H.; et al. Control of grain size and rice yield by GL2-mediated brassinosteroid responses. Nat. Plants 2015, 2, 15195. [Google Scholar] [CrossRef] [PubMed]
  22. Eremina, M.; Unterholzner, S.J.; Rathnayake, A.I.; Castellanos, M.; Khan, M.; Kugler, K.G.; May, S.T.; Mayer, K.F.X.; Rozhon, W.; Poppenberger, B. Brassinosteroids participate in the control of basal and acquired freezing tolerance of plants. Proc. Natl. Acad. Sci. USA 2016, 201611477. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Dal Santo, S.; Stampfl, H.; Krasensky, J.; Kempa, S.; Gibon, Y.; Petutschnig, E.; Rozhon, W.; Heuck, A.; Clausen, T.; Jonak, C. Stress-Induced GSK3 Regulates the Redox Stress Response by Phosphorylating Glucose-6-Phosphate Dehydrogenase in Arabidopsis. Plant Cell 2012, 24, 3380–3392. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. De Rybel, B.; Audenaert, D.; Vert, G.; Rozhon, W.; Mayerhofer, J.; Peelman, F.; Coutuer, S.; Denayer, T.; Jansen, L.; Nguyen, L.; et al. Chemical Inhibition of a Subset of Arabidopsis thaliana GSK3-like Kinases Activates Brassinosteroid Signaling. Chem. Biol. 2009, 16, 594–604. [Google Scholar] [CrossRef] [PubMed]
  25. Saidi, Y.; Hearn, T.J.; Coates, J.C. Function and evolution of ‘green’ GSK3/Shaggy-like kinases. Trends Plant Sci. 2012, 17, 39–46. [Google Scholar] [CrossRef]
  26. Dornelas, M.C.; Wittich, P.; von Recklinghausen, I.; Lammeren, A.; Kreis, M. Characterization of three novel members of the Arabidopsis SHAGGY-related protein kinase (ASK) multigene family. Plant Mol. Biol. Report. 1999, 39, 137–147. [Google Scholar] [CrossRef]
  27. Kempa, S.; Rozhon, W.; Samaj, J.; Erban, A.; Baluska, F.; Becker, T.; Haselmayer, J.; Schleiff, E.; Kopka, J.; Hirt, H.; et al. A plastid-localized glycogen synthase kinase 3 modulates stress tolerance and carbohydrate metabolism. Plant J. 2007, 49, 1076–1090. [Google Scholar] [CrossRef] [Green Version]
  28. Anuradha, S.; Rao, S.S.R. Application of brassinosteroids to rice seeds (Oryza sativa L.) reduced the impact of salt stress on growth, prevented photosynthetic pigment loss and increased nitrate reductase activity. Plant Growth Regul. 2003, 40, 29–32. [Google Scholar] [CrossRef]
  29. Krishna, P. Brassinosteroid-mediated stress responses. J. Plant Growth Regul. 2003, 22, 289–297. [Google Scholar] [CrossRef]
  30. Kagale, S.; Divi, U.K.; Krochko, J.E.; Keller, W.A.; Krishna, P. Brassinosteroid confers tolerance in Arabidopsis thaliana and Brassica napus to a range of abiotic stresses. Planta 2007, 225, 353–364. [Google Scholar] [CrossRef]
  31. Singh, I.; Shono, M. Physiological and molecular effects of 24-epibrassinolide, a brassinosteroid on thermotolerance of tomato. Plant Growth Regul. 2005, 47, 111–119. [Google Scholar] [CrossRef]
  32. Dhaubhadel, S.; Chaudhary, S.; Dobinson, K.F.; Krishna, P. Treatment with 24-epibrassinolide, a brassinosteroid, increases the basic thermotolerance of Brassica napus and tomato seedlings. Plant Mol. Biol. 1999, 40, 333–342. [Google Scholar] [CrossRef] [PubMed]
  33. Pinol, R.; Simon, E. Effect of 24-Epibrassinolide on Chlorophyll Fluorescence and Photosynthetic CO2 Assimilation in Vicia faba Plants Treated with the Photosynthesis-Inhibiting Herbicide Terbutryn. J. Plant Growth Regul. 2009, 28, 97–105. [Google Scholar] [CrossRef]
  34. Sharma, I.; Bhardwaj, R.; Pati, P.K. Stress modulation response of 24-epibrassinolide against imidacloprid in an elite indica rice variety Pusa Basmati-1. Pestic. Biochem. Physiol. 2013, 105, 144–153. [Google Scholar] [CrossRef]
  35. Zhu, T.; Deng, X.G.; Tan, W.R.; Zhou, X.; Luo, S.S.; Han, X.Y.; Zhang, D.W.; Lin, H.H. Nitric oxide is involved in brassinosteroid-induced alternative respiratory pathway in Nicotiana benthamiana seedlings’ response to salt stress. Physiol. Plant. 2016, 156, 150–163. [Google Scholar] [CrossRef]
  36. Sharma, I.; Ching, E.; Saini, S.; Bhardwaj, R.; Pati, P.K. Exogenous application of brassinosteroid offers tolerance to salinity by altering stress responses in rice variety Pusa Basmati-1. Plant Physiol. Biochem. Ppb/Soc. Fr. Physiol. Veg. 2013, 69, 17–26. [Google Scholar] [CrossRef] [PubMed]
  37. Yusuf, M.; Fariduddin, Q.; Khan, T.A.; Hayat, S. Epibrassinolide reverses the stress generated by combination of excess aluminum and salt in two wheat cultivars through altered proline metabolism and antioxidants. S. Afr. J. Bot. 2017, 112, 391–398. [Google Scholar] [CrossRef]
  38. Cui, F.; Liu, L.J.; Zhao, Q.Z.; Zhang, Z.H.; Li, Q.L.; Lin, B.Y.; Wu, Y.R.; Tang, S.Y.; Xie, Q. Arabidopsis Ubiquitin Conjugase UBC32 Is an ERAD Component That Functions in Brassinosteroid-Mediated Salt Stress Tolerance. Plant Cell 2012, 24, 233–244. [Google Scholar] [CrossRef] [Green Version]
  39. Hu, Z.J.; Lu, S.J.; Wang, M.J.; He, H.H.; Sun, L.; Wang, H.R.; Liu, X.H.; Jiang, L.; Sun, J.L.; Xin, X.Y.; et al. A Novel QTL qTGW3 Encodes the GSK3/SHAGGY-Like Kinase OsGSK5/OsSK41 that Interacts with OsARF4 to Negatively Regulate Grain Size and Weight in Rice. Mol. Plant 2018, 11, 736–749. [Google Scholar] [CrossRef] [Green Version]
  40. Xia, D.; Zhou, H.; Liu, R.J.; Dan, W.H.; Li, P.B.; Wu, B.; Chen, J.X.; Wang, L.Q.; Gao, G.J.; Zhang, Q.L.; et al. GL3.3, a Novel QTL Encoding a GSK3/SHAGGY-like Kinase, Epistatically Interacts with GS3 to Produce Extra-long Grains in Rice. Mol. Plant 2018, 11, 754–756. [Google Scholar] [CrossRef] [Green Version]
  41. Thitisaksakul, M.; Arias, M.C.; Dong, S.; Beckles, D.M. Overexpression of GSK3-like Kinase 5 (OsGSK5) in rice (Oryza sativa) enhances salinity tolerance in part via preferential carbon allocation to root starch. Funct. Plant Biol. 2017, 44, 705. [Google Scholar] [CrossRef] [PubMed]
  42. Sun, S.Y.; Wang, T.; Wang, L.L.; Li, X.M.; Jia, Y.C.; Liu, C.; Huang, X.H.; Xie, W.B.; Wang, X.L. Natural selection of a GSK3 determines rice mesocotyl domestication by coordinating strigolactone and brassinosteroid signaling. Nat Commun. 2018, 9, 1–13. [Google Scholar] [CrossRef]
  43. Sun, S.Y.; Chen, D.H.; Li, X.M.; Qiao, S.L.; Shi, C.N.; Li, C.X.; Shen, H.Y.; Wang, X.L. Brassinosteroid Signaling Regulates Leaf Erectness in Oryza sativa via the Control of a Specific U-Type Cyclin and Cell Proliferation. Dev. Cell 2015, 34, 220–228. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Wang, L.L.; Yang, Z.E.; Zhang, B.; Yu, D.Q.; Liu, J.; Gong, Q.; Qanmber, G.; Li, Y.; Lu, L.L.; Lin, Y.J.; et al. Genome-wide characterization and phylogenetic analysis of GSK gene family in three species of cotton: Evidence for a role of some GSKs in fiber development and responses to stress. BMC Plant Biol. 2018, 18, 330. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Groszyk, J.; Yanushevska, Y.; Zielezinski, A.; Nadolska-Orczyk, A.; Karlowski, W.M.; Orczyk, W. Annotation and profiling of barley GLYCOGEN SYNTHASE3/Shaggy-like genes indicated shift in organ-preferential expression. PLoS ONE 2018, 13, e0199364. [Google Scholar] [CrossRef] [PubMed]
  46. Koh, S.; Lee, S.C.; Kim, M.K.; Koh, J.H.; Lee, S.; An, G.; Choe, S.; Kim, S.R. T-DNA tagged knockout mutation of rice OsGSK1, an orthologue of Arabidopsis BIN2, with enhanced tolerance to various abiotic stresses. Plant Mol. Biol. 2007, 65, 453–466. [Google Scholar] [CrossRef]
  47. Wada, K.; Marumo, S.; Abe, H.; Morishita, T.; Nakamura, K.; Uchiyama, M.; Mori, K. A rice lamina inclination test—A micro quantitative bioassay for brassinosteroids. Agric. Biol. Chem. 1984, 48, 719–726. [Google Scholar]
  48. Gan, L.J.; Wu, H.; Wu, D.P.; Zhang, Z.F.; Guo, Z.F.; Yang, N.; Xia, K.; Zhou, X.; Oh, K.; Matsuoka, M.; et al. Methyl jasmonate inhibits lamina joint inclination by repressing brassinosteroid biosynthesis and signaling in rice. Plant Sci. 2015, 241, 238–245. [Google Scholar] [CrossRef]
  49. Bartlett, J.G.; Alves, S.C.; Smedley, M.; Snape, J.W.; Harwood, W.A. High-throughput Agrobacterium-mediated barley transformation. Plant Methods 2008, 4, 22. [Google Scholar] [CrossRef] [Green Version]
  50. Zalewski, W.; Orczyk, W.; Gasparis, S.; Nadolska-Orczyk, A. HvCKX2 gene silencing by biolistic or Agrobacterium-mediated transformation in barley leads to different phenotypes. BMC Plant Biol. 2012, 12, 206. [Google Scholar] [CrossRef] [Green Version]
  51. Wang, Z.Y.; Bai, M.Y.; Oh, E.; Zhu, J.Y. Brassinosteroid signaling network and regulation of photomorphogenesis. Annu. Rev. Genet. 2012, 46, 701–724. [Google Scholar] [CrossRef] [PubMed]
  52. Jablonski, B.; Ogonowska, H.; Szala, K.; Bajguz, A.; Orczyk, W.; Nadolska-Orczyk, A. Silencing of TaCKX1 Mediates Expression of Other TaCKX Genes to Increase Yield Parameters in Wheat. Int. J. Mol. Sci. 2020, 21, 4809. [Google Scholar] [CrossRef] [PubMed]
  53. Rozhon, W.; Wang, W.; Berthiller, F.; Mayerhofer, J.; Chen, T.; Petutschnig, E.; Sieberer, T.; Poppenberger, B.; Jonak, C. Bikinin-like inhibitors targeting GSK3/Shaggy-like kinases: Characterisation of novel compounds and elucidation of their catabolism in planta. BMC Plant Biol. 2014, 14, 172. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Yamamuro, C.; Ihara, Y.; Wu, X.; Noguchi, T.; Fujioka, S.; Takatsuto, S.; Ashikari, M.; Kitano, H.; Matsuoka, M. Loss of function of a rice brassinosteroid insensitive1 homolog prevents internode elongation and bending of the lamina joint. Plant Cell 2000, 12, 1591–1605. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Bajguz, A.; Hayat, S. Effects of brassinosteroids on the plant responses to environmental stresses. Plant Physiol. Biochem. Ppb/Soc. Fr. Physiol. Veg. 2009, 47, 1–8. [Google Scholar] [CrossRef] [PubMed]
  56. Vriet, C.; Russinova, E.; Reuzeau, C. Boosting crop yields with plant steroids. Plant Cell 2012, 24, 842–857. [Google Scholar] [CrossRef] [Green Version]
  57. Perez-Perez, J.M.; Ponce, M.R.; Micol, J.L. The UCU1 Arabidopsis gene encodes a SHAGGY/GSK3-like kinase required for cell expansion along the proximodistal axis. Dev. Biol. 2002, 242, 161–173. [Google Scholar] [CrossRef] [Green Version]
  58. Huang, H.Y.; Jiang, W.B.; Hu, Y.W.; Wu, P.; Zhu, J.Y.; Liang, W.Q.; Wang, Z.Y.; Lin, W.H. BR signal influences Arabidopsis ovule and seed number through regulating related genes expression by BZR1. Mol. Plant 2013, 6, 456–469. [Google Scholar] [CrossRef] [Green Version]
  59. Jiang, W.B.; Huang, H.Y.; Hu, Y.W.; Zhu, S.W.; Wang, Z.Y.; Lin, W.H. Brassinosteroid regulates seed size and shape in Arabidopsis. Plant Physiol. 2013, 162, 1965–1977. [Google Scholar] [CrossRef]
  60. Wu, C.Y.; Trieu, A.; Radhakrishnan, P.; Kwok, S.F.; Harris, S.; Zhang, K.; Wang, J.; Wan, J.; Zhai, H.; Takatsuto, S.; et al. Brassinosteroids regulate grain filling in rice. Plant Cell 2008, 20, 2130–2145. [Google Scholar] [CrossRef] [Green Version]
  61. Wan, T.T.; Li, Q.Z.; Lou, S.T.; Yang, Y.; Peng, L.F.; Lin, Z.Z.; Hu, Q.; Ma, L.Y. GSK3/shaggy-like kinase 1 ubiquitously regulates cell growth from Arabidopsis to Moso bamboo (Phyllostachys edulis). Plant Sci. 2019, 283, 290–300. [Google Scholar]
  62. Gearing, L.J.; Cumming, H.E.; Chapman, R.; Finkel, A.M.; Woodhouse, I.B.; Luu, K.; Gould, J.A.; Forster, S.C.; Hertzog, P.J. CiiiDER: A tool for predicting and analysing transcription factor binding sites. PLoS ONE 2019, 14, e0215495. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Hoagland, D.R.; Arnon, D.I. The water-culture method for growing plants without soil. Circ. Calif. Agric. Exp. Stn. 1950, 347, 1–32. [Google Scholar]
  64. Himmelbach, A.; Zierold, U.; Hensel, G.; Riechen, J.; Douchkov, D.; Schweizer, P.; Kumlehn, J. A set of modular binary vectors for transformation of cereals. Plant Physiol. 2007, 145, 1192–1200. [Google Scholar] [CrossRef] [Green Version]
  65. Hellens, R.P.; Edwards, E.A.; Leyland, N.R.; Bean, S.; Mullineaux, P.M. pGreen: A versatile and flexible binary Ti vector for Agrobacterium-mediated plant transformation. Plant Mol. Biol. 2000, 42, 819–832. [Google Scholar] [CrossRef]
  66. Zalewski, W.; Galuszka, P.; Gasparis, S.; Orczyk, W.; Nadolska-Orczyk, A. Silencing of the HvCKX1 gene decreases the cytokinin oxidase/dehydrogenase level in barley and leads to higher plant productivity. J. Exp. Bot. 2010, 61, 1839–1851. [Google Scholar] [CrossRef] [Green Version]
  67. Murray, M.G.; Thompson, W.F. Rapid Isolation of High Molecular-Weight Plant DNA. Nucleic Acids Res. 1980, 8, 4321–4325. [Google Scholar] [CrossRef] [Green Version]
  68. Bustin, S.A.; Benes, V.; Garson, J.A.; Hellemans, J.; Huggett, J.; Kubista, M.; Mueller, R.; Nolan, T.; Pfaffl, M.W.; Shipley, G.L.; et al. The MIQE Guidelines: Minimum Information for Publication of Quantitative Real-Time PCR Experiments. Clin. Chem. 2009, 55, 611–622. [Google Scholar] [CrossRef] [Green Version]
  69. Goto, K.; Kanazawa, A.; Kusaba, M.; Masuta, C. A Simple and Rapid Method to Detect Plant siRNAs Using Nonradioactive Probes. Plant Mol. Biol. Report. 2003, 21, 51–58. [Google Scholar] [CrossRef]
Figure 1. Relative expression of HvGSK1.1 (HORVU3Hr1G034440.2) in leaves of T4 plants grown in normal (Hoagland medium) and in salt stress conditions (Hoagland medium supplemented with NaCl 200 mM). Relative expression of HvGSK1.1 in leaves of the non-transgenic control plants grown in normal conditions was assumed as 1.0. The ADP-ribosylation factor (ARF) AJ508228 was used as the reference gene. Data represent mean values and standard deviation of at least three biological replicates with three technical repetitions. Indicated are significant differences between the control and the line grown in the same conditions: ** p < 0.01, *** p < 0.001.
Figure 1. Relative expression of HvGSK1.1 (HORVU3Hr1G034440.2) in leaves of T4 plants grown in normal (Hoagland medium) and in salt stress conditions (Hoagland medium supplemented with NaCl 200 mM). Relative expression of HvGSK1.1 in leaves of the non-transgenic control plants grown in normal conditions was assumed as 1.0. The ADP-ribosylation factor (ARF) AJ508228 was used as the reference gene. Data represent mean values and standard deviation of at least three biological replicates with three technical repetitions. Indicated are significant differences between the control and the line grown in the same conditions: ** p < 0.01, *** p < 0.001.
Ijms 21 06616 g001
Figure 2. Relative expression of HvGSK paralogs HvGSK1.2 (A), HvGSK2.1 (B), HvGSK3.1 (C), and HvGSK4.1 (D) in leaves of T4 plants of transgenic lines grown in normal (Hoagland medium, black plots) and in salt stress conditions (Hoagland medium supplemented with NaCl 200 mM, grey plots). Relative expression of each of the tested paralogs in non-transgenic control plants grown in normal conditions was assumed as 1.0. The ADP-ribosylation factor (ARF) AJ508228 was used as the reference gene. Data represent mean values and standard deviation of at least three biological replicates with three technical repetitions. Indicated are significant differences between the control and the line grown in the same conditions: * p < 0.05, ** p < 0.01, *** p < 0.001.
Figure 2. Relative expression of HvGSK paralogs HvGSK1.2 (A), HvGSK2.1 (B), HvGSK3.1 (C), and HvGSK4.1 (D) in leaves of T4 plants of transgenic lines grown in normal (Hoagland medium, black plots) and in salt stress conditions (Hoagland medium supplemented with NaCl 200 mM, grey plots). Relative expression of each of the tested paralogs in non-transgenic control plants grown in normal conditions was assumed as 1.0. The ADP-ribosylation factor (ARF) AJ508228 was used as the reference gene. Data represent mean values and standard deviation of at least three biological replicates with three technical repetitions. Indicated are significant differences between the control and the line grown in the same conditions: * p < 0.05, ** p < 0.01, *** p < 0.001.
Ijms 21 06616 g002
Figure 3. Representative pictures of leaf fragments treated with 24-epibrassinolide (EBL) 0.01 µM, bikinin 5 µM, and EBL (0.01 µM) and bikinin (5 µM) applied together. (A)–the non-transgenic control plants, (B)–line #1 with strong silencing of HvGSK1.1, (C)–line #8 with intermediate silencing of HvGSK1.1, and (D)–line #4 with alleviated silencing of HvGSK1.1. Indicated are medium inclination angles and standard deviation (SD) of at least five plants. Bar represents 5 cm.
Figure 3. Representative pictures of leaf fragments treated with 24-epibrassinolide (EBL) 0.01 µM, bikinin 5 µM, and EBL (0.01 µM) and bikinin (5 µM) applied together. (A)–the non-transgenic control plants, (B)–line #1 with strong silencing of HvGSK1.1, (C)–line #8 with intermediate silencing of HvGSK1.1, and (D)–line #4 with alleviated silencing of HvGSK1.1. Indicated are medium inclination angles and standard deviation (SD) of at least five plants. Bar represents 5 cm.
Ijms 21 06616 g003
Figure 4. Relative biomass of the non-transgenic control and the transgenic plants of T4 generation grown in normal (Hoagland medium) and in salt stress conditions (Hoagland medium supplemented with NaCl 200 mM). Data represent mean values and standard deviation of at least six plants grown in normal and salt stress conditions. Indicated are significant differences between the control and the line grown in the same conditions: * p < 0.05, ** p < 0.01, *** p < 0.001.
Figure 4. Relative biomass of the non-transgenic control and the transgenic plants of T4 generation grown in normal (Hoagland medium) and in salt stress conditions (Hoagland medium supplemented with NaCl 200 mM). Data represent mean values and standard deviation of at least six plants grown in normal and salt stress conditions. Indicated are significant differences between the control and the line grown in the same conditions: * p < 0.05, ** p < 0.01, *** p < 0.001.
Ijms 21 06616 g004
Figure 5. Representative pictures of 14-days-old seedlings of the non-transgenic control and the line #1. (A) 14-days-old seedlings of the non-transgenic control and the line #1 grown in normal conditions (Hoagland medium). (B) 14-days-old seedlings of the non-transgenic control and the line #1 grown in salt stress conditions (Hoagland medium supplemented with NaCl 200 mM).
Figure 5. Representative pictures of 14-days-old seedlings of the non-transgenic control and the line #1. (A) 14-days-old seedlings of the non-transgenic control and the line #1 grown in normal conditions (Hoagland medium). (B) 14-days-old seedlings of the non-transgenic control and the line #1 grown in salt stress conditions (Hoagland medium supplemented with NaCl 200 mM).
Ijms 21 06616 g005
Figure 6. Thousand kernels weight (TKW) in the non-transgenic control and the transgenic plants grown in soil in normal conditions. Data represent mean values and standard deviation. Indicated are significant differences (* p < 0.05) between the control and the tested lines.
Figure 6. Thousand kernels weight (TKW) in the non-transgenic control and the transgenic plants grown in soil in normal conditions. Data represent mean values and standard deviation. Indicated are significant differences (* p < 0.05) between the control and the tested lines.
Ijms 21 06616 g006
Figure 7. Representative pictures of kernels of the non-transgenic control and the transgenic lines #1 and #6.
Figure 7. Representative pictures of kernels of the non-transgenic control and the transgenic lines #1 and #6.
Ijms 21 06616 g007
Figure 8. Heat map of the Pearson correlations between expression of HvGSK1.1, expression of HvGSK paralogs, and values of thousand kernel weight (TKW) measured in the non-transgenic control and the transgenic lines grown in normal conditions (A) and in salt stress conditions (Hoagland with NaCl 200 mM) (B). Indicated are r values and p values in brackets. Black font indicates significant correlations for p < 0.05. Red blocks indicate a positive correlation and blue blocks indicate a negative correlation.
Figure 8. Heat map of the Pearson correlations between expression of HvGSK1.1, expression of HvGSK paralogs, and values of thousand kernel weight (TKW) measured in the non-transgenic control and the transgenic lines grown in normal conditions (A) and in salt stress conditions (Hoagland with NaCl 200 mM) (B). Indicated are r values and p values in brackets. Black font indicates significant correlations for p < 0.05. Red blocks indicate a positive correlation and blue blocks indicate a negative correlation.
Ijms 21 06616 g008
Figure 9. Schematic representation of cDNA of HvGSK1.1 gene (A) and pBract207 vector with HvGSK1.1 fragment built into the silencing cassette (B). Open reading frame (ORF) fragment used for construction of the silencing cassette in the pBract207 vector is indicated (A). The most important functional regions of the pBRACT207 vector are indicated (B). The open reading frame (ORF) fragment of HvGSK1.1 used for gene silencing was amplified with AmpliTaq Gold DNA Polymerase (ThermoFisher Scientific) with the primers ORFgskFR/ORFgskRe (Table 2), and the vector pCR2.1TOPO with cloned HvGSK1.1 cDNA as a template. Gel purified amplicon was blunt-end cloned into the pIPKTA38 [64]. The obtained plasmid pIPK:ORF was used as the entry clone for Gateway recombination (Invitrogen) with destination vector pBRACT207 (John Innes Center, Norwich, UK, www.bract.org).
Figure 9. Schematic representation of cDNA of HvGSK1.1 gene (A) and pBract207 vector with HvGSK1.1 fragment built into the silencing cassette (B). Open reading frame (ORF) fragment used for construction of the silencing cassette in the pBract207 vector is indicated (A). The most important functional regions of the pBRACT207 vector are indicated (B). The open reading frame (ORF) fragment of HvGSK1.1 used for gene silencing was amplified with AmpliTaq Gold DNA Polymerase (ThermoFisher Scientific) with the primers ORFgskFR/ORFgskRe (Table 2), and the vector pCR2.1TOPO with cloned HvGSK1.1 cDNA as a template. Gel purified amplicon was blunt-end cloned into the pIPKTA38 [64]. The obtained plasmid pIPK:ORF was used as the entry clone for Gateway recombination (Invitrogen) with destination vector pBRACT207 (John Innes Center, Norwich, UK, www.bract.org).
Ijms 21 06616 g009
Table 1. Agrobacterium-mediated transformation and in vitro culture of immature barley embryos.
Table 1. Agrobacterium-mediated transformation and in vitro culture of immature barley embryos.
VectorCultured Immature EmbryosRegenerated PlantsPlants with Confirmed t-DNAPlants with Single T-DNA Integration LocusPlants with Multiple T-DNA Integration Loci
T-DNA of silencing vector19751411101
Control of in vitro culture and plant regeneration107Good plant regeneration on hygromycin-free medium, 20–25 shoots from one immature embryo.
Control of in vitro culture and hygromycin selection107No shoots were regenerated on hygromycin-containing medium.
Table 2. List of primers and reaction conditions used in this study.
Table 2. List of primers and reaction conditions used in this study.
Gene Name Ensemble Plants ID NCBI IDPrimersAmplicon Length and Reaction Conditions
SymbolSequence
HvGSK1.1 HORVU3Hr1G034440.2 AK2512871FRGSK1GGTGGTGGTAGAGTAGGAGTA1719 bp, 95 °C 5 min, 35 cycles (95 °C 30 s, 60 °C 30 s, 72 °C 30 s) 72 °C 5 min.
4RevGSK1TTACAGATCAGCTATGGCAAT
ORFgskFRTGGTGAAAGTGGTGTGGACC262 bp, 95 °C 5 min, 35 cycles (95 °C 30 s, 56 °C 30 s, 72 °C 30 s) 72 °C 5 min.
ORFgskReGGTCCCGAAGCTCATCAAAG
437AK251_FrAGGGAACAGAGACTGGTCACAT108 bp, 95 °C 15 min, 45 cycles (95 °C 25 s, 60 °C 25 s, 72 °C 25 s) 72 °C 1 min.
545AK251_RevAATGAACCTTGACCAACAATCC
HvGSK1.2 HORVU5Hr1G117030.1 AK368391AK368_FRTCTGGGCACACCTACAAGGG139 bp, 95 °C 15 min, 45 cycles (95 °C 25 s, 60 °C 25 s, 72 °C 25 s) 72 °C 1 min.
AK368_RevTGGAGACCAGGTCCACTGCT
HvGSK2.1 HORVU3Hr1G026020.1 AK364823403AL364_FRAGTGCTTGGAGACTGGAGAGAC122 bp, 95 °C 15 min, 45 cycles (95 °C 25 s, 60 °C 25 s, 72 °C 25 s) 72 °C 1 min.
524AK364_ReGTGCTTCAGAGAGACGACATTG
HvGSK3.1 HORVU1Hr1G048580.7AK3625471270AK362_FrAAAGTGGCGTTGATCAGTTGG123 bp, 95 °C 15 min, 45 cycles (95 °C 25 s, 60 °C 25 s, 72 °C 25 s) 72 °C 1 min.
1393AK362_RevCAGGGATGAGCTTTTATCTGAGG
HvGSK4.1 HORVU5Hr1G119790.18AK358344AK358_FRGCGAGAAGGCAGAACCTGTT133 bp, 95 °C 15 min, 45 cycles (95 °C 25 s, 60 °C 25 s, 72 °C 25 s) 72 °C 1 min.
AK358_RevTGTCACCCACCCACACAAAG
ARF AJ508228Ref2_FRGCTCTCCAACAACATTGCCAAC162 bp, 95 °C 15 min, 45 cycles (95 °C 25 s, 60 °C 25 s, 72 °C 25 s) 72 °C 1 min.
Ref2_RevGCTTCTGCCTGTCACATACGC
ACTB (actin) AY145451qAct1AGCAACTGGGATGACATGGAG172 bp, 95 °C 15 min, 45 cycles (95 °C 25 s, 60 °C 25 s, 72 °C 25 s) 72 °C 1 min.
qAct2aCGTACATGGCAGGAACATTG
hptIIpBr_Hyg_F2GACGGCAATTTCGATGATG205 bp, 95 °C 15 min, 45 cycles (95 °C 25 s, 60 °C 25 s, 72 °C 25 s) 72 °C 1 min.
pBr_Hyg_R2CCGGTCGGCATCTACTCTAT

Share and Cite

MDPI and ACS Style

Kloc, Y.; Dmochowska-Boguta, M.; Zielezinski, A.; Nadolska-Orczyk, A.; Karlowski, W.M.; Orczyk, W. Silencing of HvGSK1.1—A GSK3/SHAGGY-Like Kinase–Enhances Barley (Hordeum vulgare L.) Growth in Normal and in Salt Stress Conditions. Int. J. Mol. Sci. 2020, 21, 6616. https://doi.org/10.3390/ijms21186616

AMA Style

Kloc Y, Dmochowska-Boguta M, Zielezinski A, Nadolska-Orczyk A, Karlowski WM, Orczyk W. Silencing of HvGSK1.1—A GSK3/SHAGGY-Like Kinase–Enhances Barley (Hordeum vulgare L.) Growth in Normal and in Salt Stress Conditions. International Journal of Molecular Sciences. 2020; 21(18):6616. https://doi.org/10.3390/ijms21186616

Chicago/Turabian Style

Kloc, Yuliya, Marta Dmochowska-Boguta, Andrzej Zielezinski, Anna Nadolska-Orczyk, Wojciech M. Karlowski, and Waclaw Orczyk. 2020. "Silencing of HvGSK1.1—A GSK3/SHAGGY-Like Kinase–Enhances Barley (Hordeum vulgare L.) Growth in Normal and in Salt Stress Conditions" International Journal of Molecular Sciences 21, no. 18: 6616. https://doi.org/10.3390/ijms21186616

APA Style

Kloc, Y., Dmochowska-Boguta, M., Zielezinski, A., Nadolska-Orczyk, A., Karlowski, W. M., & Orczyk, W. (2020). Silencing of HvGSK1.1—A GSK3/SHAGGY-Like Kinase–Enhances Barley (Hordeum vulgare L.) Growth in Normal and in Salt Stress Conditions. International Journal of Molecular Sciences, 21(18), 6616. https://doi.org/10.3390/ijms21186616

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop