Skip to main content

Molecular characterization reveals that OsSAPK3 improves drought tolerance and grain yield in rice

Abstract

Background

Many data suggest that the sucrose non-fermenting 1-related kinases 2 (SnRK2s) are very important to abiotic stress for plants. In rice, these kinases are known as osmotic stress/ABA–activated protein kinases (SAPKs). Osmotic stress/ABA–activated protein kinase 3 (OsSAPK3) is a member of SnRK2II in rice, but its function is still unclear.

Results

The expression of OsSAPK3 was up regulated by drought, NaCl, PEG and ABA. OsSAPK3 mutated seedings (sapk3-1 and sapk3-2) showed reduced hypersensitivity to exogenous ABA. In addition, under drought conditions, sapk3-1 and sapk3-2 showed more intolerance to drought, including decreased survival rate, increased water loss rate, increased stomatal conductance and significantly decreased expression levels of SLAC1 and SLAC7. Physiological and metabolic analyses showed that OsSAPK3 might play an important role in drought stress signaling pathway by affecting osmotic adjustment and osmolytes, ROS detoxification and expression of ABA dependent and independent dehydration-responsive genes. All gronomic traits analyses demonstrated that OsSAPK3 could improve rice yield by affecting the regulation of tiller numbers and grain size.

Conclusion

OsSAPK3 plays an important role in both ABA-dependent and ABA-independent drought stress responses. More interestingly, OsSAPK3 could improve rice yield by indirectly regulating tiller number and grain size. These findings provide new insight for the development of drought-resistant rice.

Peer Review reports

Background

Rice (Oryza sativa L.) is the most important staple crop worldwide. Improving rice yield is an urgent need to solve the crisis of food shortage. Rice yield mainly depends on growth environment, plant height, tillering and other important agronomic traits [1].

Under optimal environmental conditions, rice yield is mainly controlled by tiller numbers per plant, grain yield per plant and thousand-grain weight [1]. Starting with shoot branching, rice tiller experience two stages, axillary meristem formation and tiller bud outgrowth [2]. Therefore, the number and outgrowth rate of tiller buds determine final number of tillers [3]. Previous studies have reported that effective tillers per plant are essential determinant of rice yield, which are affected by genetic and environmental factors [2]. Many genes related to tiller bud formation and outgrowth have been identified in rice, such as MOC1 [4] and MOC2 [5]. Earlier studies have also reported that due to the defect of tiller bud formation, moc1 mutant plants have only a main culm without any tillers [4]. Due to insufficient supply of sucrose, moc2 mutant plants show significantly reduced tiller numbers, a reduced outgrowth rate and a dwarf phenotype [5]. More interestingly, previous study indicates that tiller enhancer (TE) controls rice tillering by mediating the degradation of MOC1 protein, which encoded an activator of the APC/CTE E3 ubiquitin ligase in rice [6]. In addition, APC/CTE activity is inhibited via phosphorylation of TE by SNF1-related protein kinases (SnRK2s) [7]. These results show that protein kinase plays an important role in regulating rice tillering.

Among environmental stresses, drought is the main threat that affects the development and growth of rice [8]. Drought decreases stomatal conductance, transpiration rate, water use efficiency, relative water content and photosynthesis rate [9]. Under drought conditions, leaf expansion, root growth, plant height and tillering are severely inhibited [10]. All of these morphological and physiological changes are responsible for a reduction in grain yield under drought condition. Drought stress occurring at booting stage [11] (Shao et al., 2014), flowering stage [12] (Liu et al., 2006) and filling stage [13] (Zhang et al., 2018) has larger detrimental influences on rice yield. Water deficit also increases the production of reactive oxygen species (ROS), which leads to peroxidation of lipids, denaturation of proteins, mutation of DNA and various types of cellular oxidative damage [10]. To overcome the damage caused by higher ROS levels, plants produce antioxidants as a tolerance mechanism [10]. Therefore, understanding the mechanism of drought stress response, especially the antioxidant mechanism under drought stress, will help to improve crop productivity under drought stress.

Numerous reports show that protein kinases involve in plant responses to biotic and abiotic stresses. With the ability to phosphorylate specific substrates, protein kinases are key components in plant drought stress response [14]. As plant specific protein kinases, SnRK2 kinases (SnRK2s) have been found in many seed plants, such as rice, maize, tobacco, tomato, wheat and soybean [15]. According to the conservation of sequences and active domains, SnRK2s are divided into three subclasses [11]. In Arabidopsis thaliana, SnRK2.1, 2.4, 2.5, 2.9, and 2.10 belong to groupI, SnRK2.7 and 2.8 to group II, and SnRK2.2, 2.3, and 2.6 to group III. In rice, these kinases are designated as osmotic stress/ABA–activated protein kinases (SAPKs). SnRK2I includes SAPK4-SAPK7, SnRK2II includes SAPK1-SAPK3, and SnRK2III includes SAPK8-SAPK10 [16]. Among SnRK2s, subclass III members act as key positive regulators in ABA signaling, which regulate the expression of stress-responsive genes in an ABA-dependent manner [17, 18]. In A. thaliana, three subclass III SnRK2s (SnRK2.2, 2.3, and 2.6) are mainly involved in ABA-dependent drought stress signaling. The snrk2.2/snrk2.6/snrk2.3 triple mutant shows severe ABA-insensitive viviparity and drought-sensitive phenotype [19,20,21,22]. In rice, members of SnRK2 subclass III (SAPK8–10) are involved in osmotic stress response by regulating ABA-dependent gene expression [23,24,25]. The protein phosphorylation network of subclass III SnRK2 has also been elucidated, which involves mitogen-activated protein (MAP) kinase (MAPK), epigenetic regulation and RNA processing [18, 26]. Subclass I SnRK2s have been identified as playing important roles in early processes of drought stress [27,28,29,30]. For example, SnRK2.4 and SnRK2.10 are activated by ABI1 and PP2CA under osmotic stress in A. thaliana [27, 28]. OsSAPK6 can phosphorylate and activate OsbZIP10 and OsbZIP46 in the ABA signaling pathway [29, 30]. For subclass II SnRK2s, SnRK2.8 positively regulates drought resistance response and participates in root growth [31]. A recent report showed that SnRK2.7 and SnRK2.8 regulated drought response genes and transcription factors bound to ABA response elements in A. thaliana [32]. Recent report indicated that SAPK2 is involved in the osmotic stress response by phosphorylating OsbZIP23 and OsbZIP46 and promoting the transcription of stress response genes [33, 34]. These data indicate that SnRK2s are crucial for abiotic stress responses in rice. Therefore, it is necessary to clarify the special function of SnRK2s.

Up till now, the function of OsSAPK3 has not been reported yet. In this study, sapk3 mutants (sapk3-1 and sapk3-2) generated by CRISPR/Cas9 system were used to characterize the role of OsSAPK3 in responding to drought stress and regulating rice yield. These results provide new insight for functional analysis of SnRK2s and engineering of drought resistant rice.

Results

Cloning and sequence analysis of OsSAPK3

The open reading frame of OsSAPK3 consists of 1,005 bp, which encoded 334 amino acids (Annotation identified in Introduction to the Rice Genome Annotation Project. database:LOC_Os10g41490 http://rice.plantbiology.msu.edu/index.shtml). SnRK2s contain a Ser/Thr enzyme activity conservative domain DFGYSKSSVLHSQPKSTVGTPAYIAPE and an ATP binding site GXGXXGX. The C-terminal sequence has high polymorphism, mainly including response to osmotic and ABA stress [16]. OsSAPK1, OsSAPK2 and OsSAPK3 are all rice SnRK2 II subfamily members. Similar to SnRK2 III members that play an important role in ABA pathway, OsSAPK3 contained several conservative domains, such as an ATP binding domain (ATP) (marked with black), the activation loop (marked with red lines) and SnRK2s conserved motif (marked with blue lines) (Fig. 1 A, Supplementary Table 1). But there is a significant difference in the ABA-activated motif (ABA box) (marked with yellow lines) comparing with ABA-pathway SnRK2s (Fig. 1 A, Supplementary Table 1).

Fig. 1
figure 1

Sequence analysis and phylogenetic relationship of SnRK2II and SnRK2III subfamily genes. A Amino acids comparison and B Phylogenetic relationship of selected SnRK2 members in Arabidopsis thaliana (At), Oryza sativa (Os) and Zea mays (Zm). The ATP-binding domain (ATP) and the activation loop are marked with black and red lines respectively. The SnRK2s conserved motif and the ABA-activated SnRK2s motif (ABA box) are marked with blue and yellow lines respectively

To analyze the closeness of OsSAPK3 to ABA dependent SnRK2 III members, a phylogenetic tree was constructed based on their amino acids sequence of selected SnRK2 II and SnRK2 III members from A. thaliana (At), Oryza sativa (Os) and Zea mays (Zm) (Supplementary Table 1). These results showed that OsSAPK3 was similar to ZmSAPK3 (Fig. 1 B). In addition, it is highly similar to AtSnRK2.8 in A. thaliana and OsSAPK2 in rice (Fig. 1 B). These results indicated that OsSAPK3 might have functions different from subclass III SnRK2s but similar to OsSAPK2 and AtSnRK2.8.

Expression profile and localization of OsSAPK3

To analyze the responsiveness of OsSAPK3 under different abiotic stresses in rice, OsSAPK3 expression profiles were conducted via quantitative RT-PCR (qRT-PCR). Under drought, NaCl and PEG treatments, the expression of OsSAPK3 remarkably increased (Fig. 2 A). In ABA-treated plants, OsSAPK3 expression was weakly induced (Fig. 2 A).These results indicated that OsSAPK3 played a role in osmotic stress response, but had moderate role in ABA tolerance.

Fig. 2
figure 2

The expression pattern of OsSAPK3 A Transcription analyses of OsSAPK3 under various stress, such as 150 mM NaCl, 20% PEG6000, exposed to the air or 100 μM ABA. B Transcription analyses of OsSAPK3 in different rice tissues by qRT-PCR. C GUS staining analyses of ProSAPK3-GUS expression in different rice tissues. a, root; b, stem; c, leaf. (D) Subcellular localization of OsSAPK3 in rice protoplasts. Scale bar = 10 μm. Data represent means ± SD of three biological replicates (5 plants for each replicate)

To investigate the biological role of OsSAPK3, the spatial and temporal expression of OsSAPK3 in different tissues were analyzed. The qRT-PCR results displayed that OsSAPK3 expression level was highest in roots, followed by leaves (Fig. 2 B). To further verify this result, ProOsSAPK3-Gus plants were used for detecting tissue expression patterns. The GUS activity was relatively higher in roots than in leaves (Fig. 2 C), which was consistent with qRT-PCR analyses. The subcellular localization of OsSAPK3 was assessed through transient expression of the OsSAPK3-GFP fusion protein in rice protoplasts. Fluorescence microscopy analyses showed that OsSAPK3 was localized in the nucleus and cytoplasm (Fig. 2 D).

OsSAPK3 mutated plants show decreased ABA sensitivity

To explore the function of OsSAPK3 in responses to ABA stress, OsSAPK3 mutated plants were generated using CRISPR/Cas9 gene editing technology. Two homozygous lines were chosen for ABA tolerance evaluation, which were named sapk3-1 and sapk3-2.

The sapk3-1 plants contained a 7-bp deletion in the first exon of OsSAPK3, while sapk3-2 plants carried a 14-bp deletion in second exon (Fig. 3 A, Supplementary Table 2). Amino acid sequence analyses showed that OsSAPK3 protein from sapk3-1 and sapk3-2 mutant plants had only 5 or 48 amino acids due to premature termination of translation. (Fig. 3 B, Supplementary Table 2). In addition, OsSAPK3 expression level was much lower in sapk3-1 and sapk3-2 plants than in wild-type plants (Fig. 3 C). These results indicated that sapk3-1 and sapk3-2 mutant lines were OsSAPK3 loss-of-function mutants.

Fig. 3
figure 3

OsSAPK3 mutated plants show decreased ABA sensitivity. A Sequence comparison of wild-type plants and two sapk3 mutant lines. B Protein structures of wild-type plants and two sapk3 mutant lines. C Transcript accumulation of OsSAPK3 in wild-type plants and two sapk3 mutant lines. D Phenotypes of wild-type plants and two sapk3 mutant lines before and after 100 μM ABA treatment for10 days. Scale bars = 2 cm. E Shoot and F Root length corresponding to D. Data represent means ± SD of three biological replicates (16 plants for each replicate). **P < 0.01 (Student’s t-test)

In order to determine whether OsSAPK3 participates in ABA responses, shoot and root length of the wild-typeplants and two sapk3 mutant lines were examined after treatment with 100 μM exogenous ABA for 10 days. Compared with wild-type plants, the shoot and root growth inhibition significantly reduced in two sapk3 mutant lines after ABA treatment (Fig. 3 D-F), which indicated that sapk3-1 and sapk3-2 plants were not sensitive to ABA. These results demonstrated that OsSAPK3 played a role in ABA dependent signaling pathway.

OsSAPK3 mutated plants show reduced drought stress tolerance

To confirm the function of OsSAPK3 in drought tolerance in rice, the performance of wild-type plants and two sapk3 mutant lines under drought stress was examined by withholding water for 7 days. The survival rate of wild-type plants (46.4%) was significantly higher than that in two sapk3 mutant lines (i.e., sapk3-1 and sapk3-2 plants survival rates were only 20.4% and 22.5%, respectively) (Fig. 4 A and B). This result demonstrated that OsSAPK3 played a positive role in drought stress response.

Fig. 4
figure 4

OsSAPK3 mutated plants show reduced drought stress tolerance A Phenotype of wild-type plants and two sapk3 mutant lines before and after drought stress treament for 7 days. Scale bars = 2 cm. B Survival rates corresponding to A. C Relative water loss rate of leaves from 2-week-old wild-type plants and two sapk3 mutant lines under drought stress. D Stomatal density and E Stomatal conductance from wild-type plants and two sapk3 mutant lines before and after drought stress. F Transcript level of stomata genes OsSLAC1 and G OsSLAC7 in wild-type plants and two sapk3 mutant lines before and after drought stress. Values represent the means ± SD of three biological replicates (5 plants for each replicate). **P < 0.01 (Student’s t-test)

Since the drought tolerance is closely related to relative water loss rate, the relative water loss rate was observed. As shown in Fig. 4, the relative water loss rate of wild-type leaves was higher than that in two sapk3 mutant lines (Fig. 4 C). Stomatal status is important for drought tolerance, then the stomatal density and stomatal conductance at 6-leaves stage were examined. There was no obvious difference for stomatal density between the wild-type plants and two sapk3 mutant lines (Fig. 4 D). Under normal conditions, the stomatal conductance of all plants had no significant difference between wild-type plants and sapk3 mutants. However, stomatal conductance of wild-type plants was significantly lower than that in two sapk3 mutant lines under drought conditions (Fig. 4 E). Based on these results, OsSAPK3 may change the stomatal conductance by participating in ABA dependent stomatal regulation. To test this assumption, the expression levels of slow anion channel-associated 1 (OsSLAC1) and slow anion channel-associated 7 (OsSLAC7) were verified. The expression levels of these two genes in two sapk3 mutant lines were significantly lower than that in wild-type plants under drought conditions (Fig. 4 F and D).

These results indicated that OsSAPK3 improved the tolerance to drought stress by regulating stomatal conductance through regulating the expression of OsSLAC1 and OsSLAC7.

Metabolic regulation of sapk3 mutants in response to drought stress

Plants adapt to drought stress through physiological and biochemical changes, such as the accumulation of proline and soluble sugar [23]. Thus, the content of proline and soluble sugar was examined. Under normal conditions, there is no significant difference in proline and soluble sugar content among all plants (Fig. 5 A and B). After drought treatment, the proline and soluble sugar content in two sapk3 mutant lines were significantly lower than that in wild-type plants (Fig. 5 A and B). Proline biosynthesis is catalyzed by 1-pyrroline-5- carboxylate synthetase (OsP5CS) [35]. According to qRT-PCR data, the OsP5CS expression level was much lower in sapk3 mutants than in wild-type plants under drought conditions (Fig. 5 E), which was consistent with proline content. These results suggested that OsSAPK3 might positively regulate the resistance to drought stress via changing the accumulation of osmolyte, such as proline and soluble sugars.

Fig. 5
figure 5

Metabolic regulation of sapk3 mutants in response to drought stress. A Proline B Soluble sugar C MDA contents in the seedlings of wild-type plants and two sapk3 mutant lines before and after drought stress treament for 7 days. D Relative ion leakage in wild-type plants and two sapk3 mutant lines leaves before and after drought stress treament for 7 days. E Relative expression levels of OsP5CS in the seedlings of wild-type plants and two sapk3 mutant lines before and after drought stress treament for 7 days. F AsA content in the seedlings of wild-type plants and two sapk3 mutant lines before and after drought stress treament for 7 days. G CAT H SOD I POD activities in the seedlings of wild-type plants and two sapk3 mutant lines before and after drought stress treament for 7 days. Values represent the means ± SD of three biological replicates (5 plants for each replicate). **P < 0.01 (Student’s t-test). FW: Fresh weight

OsSAPK3 affects ROS detoxification via regulating expression of antioxidant genes under drought conditions

Plants under abiotic stress generate various kinds of ROS molecules [36]. ROS concentrations is frequently quantified by measuring the amount of Malondialdehyde (MDA) and relative electrolytic leakage [37]. MDA content and relative electrolytic leakage of the leaves from wild-type plants and two sapk3 mutant lines were compared. After drought treatment, two sapk3 mutant lines accumulated more MDA than the wild-type plants (Fig. 5 C). Moreover, relative electrolytic leakage of two sapk3 mutant lines was higher than that in wild-type plants (Fig. 5 D). Given the increased ROS-induced damage of sapk3 mutants, the activities of catalase (CAT), superoxide dismutase (SOD), peroxidase (POD) and the content ascorbic acid (AsA) from wild-type plants and two sapk3 mutant lines were also measured. Compared with wild-type plants, two sapk3 mutant lines showed lower CAT, SOD, POD activities and AsA content than wild-type plants under drought conditions (Fig. 5 F-I).

To determine whether OsSAPK3 is involved in ROS-scavenging, the expression levels of several antioxidant genes including OsCAT, OsSOD1, OsSOD2 and ascorbic acid peroxidase 2 (OsAPX2) were assessed. The expression levels of these genes were significantly lower in two sapk3 mutant lines than that in wild-type plants (Fig. 6A-D). These results demonstrated that OsSAPK3 was also important for improving ROS tolerance by regulating ROS detoxification.

Fig. 6
figure 6

The expression of abiotic stress-responsive genes. A-I Relative expression levels of OsCAT, OsSOD1, OsSOD2, OsAPX2, OsTRAB1, OsbZIP23, OsbZIP46, OsABI5 and OsDREB2A from the seedings of 2-week-old wild-type plants and two sapk3 mutant lines before and after drought stress. Values represent the means ± SD of three biological replicates (5 plants for each replicate). **P < 0.01 (Student’s t-test)

OsSAPK3 involved in drought stress response in ABA dependent and independent manner

In order to explore the possible mechanism of OsSAPK3 in the regulation of drought stress, the expression of 5 well-characterized drought resistance-related genes (i.e., OsTRAB1 [38], OsABI5 [39], OsbZIP23 [34], OsbZIP46 [33] and OsDREB2A [40]) were analyzed. Among these genes, OsTRAB1, OsABI5, OsbZIP23 and OsbZIP46 are involved in ABA dependent signaling, OsDREB2A is involved in ABA independent regulatory systems. All their expression levels were significantly lower in two sapk3 mutant lines than that in wild-type plants after drought treatment (F ig. 6E-I). These results suggested that OsSAPK3 upregulated the expression of some stress-responsive genes both in ABA-dependent and ABA-independent manner under drought stress.

OsSAPK3 improved rice yield by altering tiller numbers and grain size

In order to assess the effects of OsSAPK3 on rice growth and productivity, the wild-type plants and two sapk3 mutant lines were grown in the field. At maturity, compared to wild-type plants, plant height of two sapk3 mutant lines decreased significantly (Fig. 7 A and C). In addition to dwarf phenotype, two sapk3 mutant lines showed significantly fewer tiller numbers than wild-type plants (Fig. 7 E). However, effective tiller numbers of sapk3 mutant plants were slightly fewer than wild-type plants (Fig. 7 F). The effective tiller rate of sapk3 mutants was significantly higher than wild-type plants (Fig. 7 G). These results indicated that OsSAPK3 participated in the regulation of rice tillering.

Fig. 7
figure 7

OsSAPK3 affected growth phenotypes in rice A Growth phenotypes B Grain phenotypes C Plant height D Grain length E Tiller numbers F Effective tiller numbers G Effective tiller rate H Grain yield per plant of wild-type plants and two sapk3 mutant lines.Values represent the means ± SD of three biological replicates (25 plants for each replicate). *P < 0.05 and **P < 0.01 (Student’s t-test)

The tiller number is one key factor determining rice grain yield. Efficient utilization of nitrogen (N) is one important factor affecting rice tillering [41]. So far, many nitrate transporter genes, such as nitrate and di/tripeptide transporters (NPFs), have been found in rice [42]. OsNPF7.1, OsNPF7.2 and OsNPF7.4 was reported to positively regulate the tiller number [41, 42]. Additionally, OsSAPK2, which is homologous to OsSAPK3, could enhance grain production by regulating nitrogen utilization efficiency [43]. The function of OsSAPK3 on rice tillering was investigated by analysing tiller buds development of wild-type plants and two sapk3 mutant lines. Tiller buds grew slower in two sapk3 mutant lines than that in wild-type plants, and this phenomenon was observed continuously for 34 days after germination (DAG) (Fig. 8 A and B). Based on these results, it could be concluded that OsSAPK3 promoted rice tiller bud outgrowth between 20 and 34 DAG. To investigate how OsSAPK3 altered the tiller buds, the expression levels of OsNPF7.1, OsNPF7.2 and OsNPF7.4 were measured. These results showed that the expression level of OsNPF7.2 in two sapk3 mutant lines was significantly lower than that in wild-type plants (Fig. 8 D). However, expression levels of OsNPF7.1 and OsNPF7.4 exhibited no significant difference between wild-type plants and two sapk3 mutant lines (Fig. 8 C and E). These results suggested that OsSAPK3 positively regulated rice tillering by promoting rice tiller bud outgrowth through regulating expression of nitrate transporter genes.

Fig. 8
figure 8

OsSAPK3 altered tiller bud outgrowth in rice A Comparison of tiller buds at 20 days after germination among wild-type plants and two sapk3 mutant lines. bar = 1 cm. B Statistical analysis at different days after germination among wild-type plants and two sapk3 mutant lines. Values represent the means ± SD of 3 biological replicates (25 plants for each replicate). C Expression levels of OsNPF7.1 D OsNPF7.2 E OsNPF7.4 F OsMIS2, G OsGSN1 and H OsGIF1 among wild-type plants and two sapk3 mutant lines. Values represent the means ± SD of three biological replicates (5 plants for each replicate). **P < 0.01 (Student’s t-test)

On the other hand, rice yield is also estimated via panicle length, grain number, seed setting rate per panicle, grain size and grain yield per plant [1]. Analyses of panicle presented that there was no significant difference in panicle length, grain number and seed setting rate per panicle between two sapk3 mutant lines and wild-type plants (Fig. S1 A and C-E). This result suggested that OsSAPK3 might not participate in the regulation of pollen fertility. In addition, there was no significant difference in grain width and 1000-seed weight (Fig. S1 B and F). Further analyses of grain size and yield revealed that grain length and grain yield per plant of two sapk3 mutant lines were significantly lower than that of wild-type (Fig. 7 B, D and H).

Several genes related to grain size in rice have been identified, such as GRAIN SIZE MINI SEED 2 (MIS2) [44], AND NUMBER1(GSN1) [45] and GRAIN INCOMPLETE FILLING 1 (GIF1) [46]. In order to investigate how OsSAPK3 changed grain length, the expression levels of OsMIS2, OsGSN1 and OsGIF1 were measured. These results demonstrated that the expression levels of OsMIS and OsGIF1 in two sapk3 mutant lines were significantly lower than that in wild-type plants (Fig. 8 F and H). But the expression level of OsGSN1 showed no significant difference between two sapk3 mutant lines and wild-type plants (Fig. 8 G). These results suggested that OsSAPK3 regulated grain length by promoting expression levels of seed size related genes.

Discussion

OsSAPK3 confers drought tolerance in an ABA-dependent manner in rice

Plant-specific kinases, SnRK2s are considered to be major regulators in the plant responses to challenging environmental conditions [17]. OsSAPK2, member of subclass II SnRK2s, is involved in the osmotic stress response by regulating the expression of ABA dependent genes [47, 48]. These reports suggest that subclass II SnRK2s play important roles in osmotic stress signaling and promotes us to explore the potential role of OsSAPK3. However, direct experimental evidence of the role of OsSAPK3 in abiotic stress is largely lacking.

In this study, sequence alignment showed that OsSAPK3 differed dramatically from subclass III SnRK2s in terms of ABA-activated motif (marked with yellow lines) (Fig. 1 A). But it was highly similar to AtSnRK2.8 in A. thaliana and OsSAPK2 in rice (Fig. 1 B). SnRK2.8 is strongly activated by salt and mannitol, slightly activated by ABA [31]. OsSAPK2 is strongly induced by drought, NaCl and PEG, but not by ABA in rice [49]. It is found that OsSAPK3 was significantly induced by drought, NaCl and PEG stress, and weakly induced by ABA (Fig. 2 A). In addition, OsSAPK3 was highly expressed in roots and weakly in leaves (Fig. 2 B and C) which was consistent with the expression patterns of AtSnRK2.8 [31], but different from OsSAPK2 [49]. These findings suggested that OsSAPK3 might acquire functions different from subclass III SnRK2s but similar to OsSAPK2 and AtSnRK2.8 during long-term evolution.

Mizoguchi et al. report that AtSnRK2.8 play important roles in drought stress signaling in A. thaliana [32]. Here, the results demonstrated that OsSAPK3 was localized in the nucleus and cytoplasm (Fig. 2 D), which was quite similar to AtSnRK2.8 [32]. In addition, the hypersensitivity to exogenous ABA reduced in sapk3 mutants lines (Fig. 3 D-F). Further research showed that sapk3 mutants exhibited decreased drought tolerance (Fig. 4 A), as charactered by lower survival rate (Fig. 4 B) and higher water loss rate under drought conditions (Fig. 4 C). Under environmental stresses, ABA induces stomatal closure through the activation and inactivation of ion channels, such as SLAC1 and SLAC7 in A. thaliana [50, 51]. In rice, SAPK8 has been found to activate the function of OsSLAC1 [52]. SAPK2 [29] and SAPK9 [53] positively regulated the stomatal closure by up regulating the expression of OsSLAC1 and OsSLAC7. Under drought conditions, stomatal conductance in sapk3 mutant plants leaves was remarkably increased (Fig. 4E), and expression patterns of SLAC1 and SLAC7 in sapk3 mutant plants leaves were remarkably decreased compared with wide type plants (Fig. 4 F and G). These results suggested that OsSAPK3 functioned positively in drought stress response by ABA-dependent manner.

OsSAPK3 affects physiological and metabolic regulation in response to drought stress

In rice, understanding the physiological adaptation of drought resistance is very useful for developing high-yield varieties under drought conditions. Previous studies have shown that plant growth is severely restricted by abiotic stresses, mainly due to protein denaturation, lipid peroxidation, cell homeostasis destruction, cell oxidative damage and DNA mutations caused by the accumulation of ROS [54, 55]. Lipid peroxidation induced by ROS accumulation can damage cell membranes [55]. Therefore, relative ion leakage is a good indicator to evaluate membrane stability under cold, drought or salt stress conditions in rice and wheat [56, 57]. The results in this study revealed that MDA content and relative ion leakage in sapk3 mutants were significantly higher than that in wild-type plants (Fig. 5 C and D), suggesting that OsSAPK3 played a crucial role in alleviating oxidative damage in rice.

According to previous reports, ROS detoxification plays a key role in drought and heat stress response [58, 59]. The antioxidant system that protects plants against adverse effects of ROS contains enzymatic antioxidants and non enzymatic molecules. Enzymatic antioxidants include SOD, CAT and APX [60]. AsA serves as non enzymatic antioxidant [60]. Therefore, CAT, SOD, POD activities and AsA content can be used as important indicators for drought tolerance in plants. In this study, activities of CAT, SOD, POD and AsA content of two sapk3 mutant lines were significantly lower than that in wild-type plants (Fig. 5 F-I). To understand the molecular mechanism of decreased ROS levels and higher antioxidant activities in two sapk3 mutant lines, the expression levels of OsCAT, OsSOD1, OsSOD2 and OsAPX2 were tested. In this study, the expression levels of OsCAT, OsSOD1, OsSOD2 and OsAPX2 exhibited significant downregulation in two sapk3 mutant lines (Fig. 6 A-D), which was consistent with lower CAT, SOD, POD activities and AsA content in two sapk3 mutant lines (Fig. 5 F-I). These results supported that OsSAPK3 regulated the accumulation of ROS through altering the expression of ROS-scavenging genes under drought stress. ABA not-activated SnRK2s, such as SnRK2.4 and SnRK2.10, take part in ROS production/accumulation under salt or high osmotic stress conditions in A. thaliana [61]. In rice, OsSAPK2 improved ROS detoxification by promoting the generation of ROS scavengers under salt or drought stress [33, 42]. SAPK9 decreased cellular oxidative damage by reducing ROS accumulation [53]. These results suggested that OsSAPK3 could be involved in the SnRK2-dependent ROS production and scavenging.

Early studies showed that under dehydration conditions, the accumulation of compatible osmoregulation substances such as proline and soluble sugar could certify the normal metabolic pathway by providing osmotic balance [54]. Under drought stress, plants can repair damage by accumulating proline to increase antioxidant activity [34]. Therefore, proline content can be used as an important reference for drought resistance. In transgenic rice, overexpression of OsP5CS showed a great improvement for drought tolerance [55]. SAPK2 [29] and SAPK9 [53] positively regulates drought or salt stress tolerance by increasing the osmotic adjustment and stomatal closure in rice. To test if OsSAPK3 plays a role in proline and soluble sugar accumulation under drought stress, proline and soluble sugar contents were examined. Under drought condition, contents of proline and soluble sugar in two sapk3 mutant lines were significantly reduced compared with wild-type plants (Fig. 5 A and B). Along with this, relative expression of OsP5CS was significantly reduced in two sapk3 mutant lines (Fig. 5 E). These results implied that OsSAPK3 functioned positively by influencing accumulation of compatible osmolytes in drought stress response.

Subclass III SnRK2s phosphorylate and activate several transcription factors, such as ABA-responsive element binding factor (ABF/AREB) transcription factors OsTRAB1 (OsbZIP66) [38], OsABI5 [39], OsbZIP23 [34] and OsbZIP46 [33]. Dehydration responsive element-binding protein 2A (OsDREB2A) mediates transcriptional changes to acquire stress resistance in ABA-independent pathways [40]. In this study, the transcript levels of OsTRAB1, OsABI5, OsbZIP23 and OsbZIP46 were all significantly lower in two sapk3 mutant lines than in wild-type plants under drought conditions (Fig. 6 E–H). In addition, the analyses of OsDREB2A expression level showed that there was also significantly decreased in two sapk3 mutant lines than in wild-type plants under drought conditions (Fig. 6I). As ABA-responsive bZIP transcription factors, OsbZIP23 and OsbZIP46 can be activated by ABA-dependent SnRKs including SAPK2 [29], SAPK6 [29] and SAPK9 [53]. These results suggested that OsSAPK3 had several substrates in both ABA-dependent and ABA-independent signaling.

In total, OsSAPK3 might play an important role in both ABA-dependent and ABA-independent drought stress signaling pathways.

OsSAPK3 improved rice yield by altering tiller numbers and grain length

It has reported that overexpression of SRK2C/SnRK2.8 enhances plant growth in A. thaliana, which may due to the phosphorylation of enzymes involved in metabolic processes [31]. SAPK9 increases rice yield by influencing panicle weight and spikelet fertility [53]. Under drought stress in reproductive period, sapk2 mutant showed dwarf phenotype, with decreased grain number per panicle and grain yield per plant [40]. To gain more comprehensively understanding of OsSAPK3, plants height between wild-type plants and two sapk3 mutant lines was evaluated. The plant height of two sapk3 mutant lines decreased significantly (Fig. 7 A and C). These results implied that OsSAPK3 had a similar function to AtSnRK2.8 or OsSAPK2 in regulating plant growth.

Previous studies have reported that effective tillers per plant is an essential determinant of rice yield, which is affected by genetic and environmental factors [2]. TE controls tillering and shoot branching [6]. In addition, TE-OE lines have much fewer tillers, whereas the loss-of-function te mutants have much more tillers [6]. More interesting, subclass III SAPKs (SAPK8, SAPK9 and SAPK10) are involved in shoot branching and tillering by inhibiting APC/CTE activity through phosphorylating TE [7]. In addition to dwarf phenotype, tiller numbers and effective tillers reduced significantly in two sapk3 mutant lines (Fig. 7 E and F). These results indicated that, similar to subclass III SAPKs, OsSAPK3 might participated in the regulation of rice tillering.

In previous research, OsSAPK2 can improve grain yield by adjusting nitrogen utilization efficiency [40]. Rice contains a large gene family, which is used for the uptake and transport of nitrate and small peptides. For example, OsNPF7.1, OsNPF7.1 or OsNPF7.4 participate in the regulation of crop yield by promoting the growth of axillary buds and by increasing the number of tillers [38]. Overexpression of OsNPF7.1, OsNPF7.1 and OsNPF7.4 could enhance grain yield by improving the uptake of nitrate [39]. In order to further explore its possible regulatory mechanism, the development of tiller buds and the expression levels of OsNPF7.1, OsNPF7.2 and OsNPF7.4 between two sapk3 mutant lines and wild-type plants were investigated. These results showed that tiller buds grew shorter in two sapk3 mutant lines (Fig. 8 A and B). In addition, the expression level of OsNPF7.2 in two sapk3 mutant lines was significantly lower than that in wild-type plants (Fig. 8 D). These results suggested that OsSAPK3 regulated rice tiller numbers by affecting the growth of axillary buds via promoting expression of nitrate transporter genes.

On the other hand, results from analysis of other agronomic traits showed that, compared with wild-type plants, grain length and grain yield per plant decreased significantly in two sapk3 mutant lines (Fig. 7 D and H). However, there is no significant change in other traits, such as grain width, panicle length, grain number per panicle, seed setting rate per panicle and 1000-seed weight (Fig. S1 A-F). Grain size is directly associated with grain yield in rice. In rice, loss-of-function of GSN1 caused larger grains by decreasing cell number due to reduced cell division during spikelet development [45]. MIS2 controls grain size by regulating epidermal cell size and cell number. In mis2 mutant, the grain showed reduced length, width and thickness [44]. Overexpression of GIF1, which encodes a cell-wall invertase, increased grain size [46]. In this study, the expression levels of OsMIS and OsGIF1 in two sapk3 mutant lines was significantly lower than that in wild-type plants (Fig. 8 F and H). These results suggested that OsSAPK3 regulated grain length by promoting expression levels of seed size related genes. In conclusion, these results demonstrated that OsSAPK3 improved rice yield by participating in the regulation of tiller numbers and grain size indirectly.

Conclusions

In this study, OsSAPK3 was significantly induced by drought, NaCl and PEG stress, and weakly induced by ABA. OsSAPK3 was localized in the nucleus and cytoplasm. In order to characterize the role of OsSAPK3 in drought stress response, OsSAPK3 mutated plants (sapk3-1 and sapk3-2) was generated using CRISPR/Cas9 system. sapk3-1 and sapk3-2 seedings reduced the hypersensitivity to exogenous ABA. In addition, sapk3-1 and sapk3-2 showed decreased stress tolerance, accompanied by lower survival rates, higher water loss rate, increased stomatal conductance and remarkably decreased expression of SLAC1 and SLAC7 under drought stress. These results suggested that OsSAPK3 functioned positively in drought stress response by ABA-dependent manner. In addition, physiological and metabolic analyses showed that OsSAPK3 played an important role in both ABA dependent and independent abiotic stress signaling pathway, including influencing accumulation of compatible osmolytes, ROS detoxification, expression of ABA dependent and independent dehydration-responsive genes. Results of agronomic traits demonstrated that OsSAPK3 improved rice yield by regulating tiller numbers and grain size. The experimental of tiller bud development suggested that OsSAPK3 regulated rice tiller numbers by affecting axillary bud growth via promoting the expression of nitrate transporter genes. Further studies suggested that OsSAPK3 regulated grain length by promoting expression levels of seed size related genes. These findings can be used to improve rice yield under drought stress. However, the in-depth mechanism needs further study.

Materials and methods

Plant cultivation and agronomic traits analyses

Wild-type rice plants (Oryza sativa L.) used in this paper, were from the laboratory of XiaoYan Yang, Faculty of Life Science and Technology, Kunming University of Science and Technology.

All seeds of different genotypes were harvested at the same time and stored in the same environment. Seeds of wild-type plants and two sapk3 mutant lines germinated on 1/2 Murashige and Skoog (MS) medium simultaneously for one week.

For ABA sensitivity assay, one week seedlings were transplanted into 1/2 MS medium containing 100 μM ABA. After 10 days treatment, shoot and root height were measured.

For abiotic stress experiment, one week seedlings were transplanted into the pot with same amount of soil, The drought phenotype of two weeks seedlings was identified after withdrawing water for 7 days and re-watering for 7 days.

For basic agronomic traits analysis, rice plants were grown in experimental plot of rice in Yuxi Normal University from March to August. Ten plants at a spacing of 16.5 cm × 26.5 cm were planted in a row and 5 rows of each line were planted. At reproductive stage, 25 plants of wild-type plants and two sapk3 mutant lines were randomly chosen to detect agronomic traits.

Generation and screening of transgenic plants

OsSAPK3 mutants were generated by the CRISPR/Cas9 system. The CRISPR/Cas9 plasmid was designed according to the protocol described previously [62]. Concisely, the first and second exons of SAPK2 were selected for guide RNA design. Double-strand DNA generated by annealing the oligo pairs, and then was cloned into the pYLCRISPR/Cas9Pubi-H vector. Wild-type rice plants (Oryza sativa L.) was used for transformation. Then transgenetic seedlings were kept in growth chamber at 28C under long-day conditions (14 h light/10 h dark cycles). For mutation detection, genomic DNA extracted from mutant seedlings (all plant) were used for PCR. Then PCR products (sequence is in Supplementary Table 3) were identified by comparing the 19-bp gRNA target sequences (gtttcgagggggccaatgg, gattgatgaaaatgtgcag) to the rice reference genome (sequence is in Supplementary Table 3).

Based on mutation detection results, we selected two independent homozygous mutant lines in the T1 generation, which are named sapk3-1 and sapk3-2. The primers used for CRISPR/Cas9 (U6a-SAPK3-F, U6a-SAPK3-R, U3-SAPK3-F, U3-SAPK3-R) and mutation detection (Cas9-SAPK3-F and Cas9-SAPK3-R) were listed in the Supplementary Table 3. Nucleotide sequence of SAPK3 CDS and predicted amino acid sequences in different genotypes were listed in the Supplementary Table 2.

Phylogenetic tree and alignment

Sequences of selected SnRK2II and SnRK2III members in A. thaliana (At), Oryza sativa (Os) and Zea mays (Zm) were aligned with CLUSTALX (Sequences were listed in the Supplementary Table 1) [16, 63]. The phylogenic tree was constructed by MEGA 5.0. The reliability for the internal branch was evaluated by the bootstrap with 1000 bootstrap replicates and marked above the nodes.

RNA extraction and qRT-PCR analyses

To detect the transcript level of target genes under different stresses, wild-type plants were kept in pot with same soil in growth chamber at 28C under long-day conditions (14 h light/10 h dark cycles) for 2 weeks. Two weeks old seedlings were treated with 100 μM ABA, 175 mM NaCl or 20% PEG 6000 (m/v) for 24 h. For drought stress, two weeks old seedlings were exposed to the air without water supply for 24 h. Leaves were sampled every 3 h.

For the qRT-PCR analysis, total RNA was isolated from rice seedlings using the TriZol reagent (Invitrogen). The cDNAs were obtained by using Superscript II in accordance with manufacturer’s instructions (Invitrogen). The qRT-PCR analysis was performed using SYBR Premix Ex Taq kit (Takara). The primers used in qRT-PCR analysis were listed in Supplementary Table 3.

GUS Staining of OsSAPK3

For GUS reporter analysis, OsSAPK3 promoter was amplified using primers Promoter-SAPK3-F and Promoter-SAPK3-R (Supplementary Table 3). The fused ProSAPK3-GUS was cloned into the P1300 vector. For the GUS staining analysis, leaves, stems and roots were destained with pure ethanol and then examined. The GUS activity was detected using the GUS blue kit (Huayueyang, Beijing, China) according to the manufacturer’s manual.

Subcellular localization

For the transient expression assay in rice protoplasts, the full-length OsSAPK3 CDS sequences were inserted into the 35S: GFP (PAN580) empty vector to generate the green fluorescent protein (GFP) recombinant vector 35S-SAPK3-GFP (pAN580). The 35S-GFP empty vector and the 35S-SAPK3-GFP recombinant vector were inserted into separate rice protoplasts as previously described [64]. Fluorescent signals were detected using the Zeiss LSM 710 laser scanning confocal microscope.

Determination of metabolites

Water loss rate, relative ion leakage, the content of proline, soluble sugar and malondialdehyde (MDA) were checked following the method as previously described [49].

For water loss rate, flag leaves were detached and left on the laboratory bench at room temperature. Then weighed every 20 min. Water loss rate was expressed as a percentage of initial fresh weight.

For relative ion leakage, About 1 g flag leaves separated from different lines were cut into 5 mm length and placed in test tubes containing 10 ml deionized water. The tubes were covered with plastic caps and placed in a water bath maintained at the constant temperature of 22 °C for 2 h. The conductance of H2O was measured by conductivity meter (HORIBA TWIN COND B-173).

For the content of proline, about 0.5 g dried leaf segments from different lines were ground into powder with liquid nitrogen, and then homogenized with 10 ml of 3% sulphosalicylic acid in tube. Collected the supernatant after centrifuging for 20 min (3000 × g). 2 ml of supernatant was reacted with 2 ml acid ninhydrin and 2 ml glacialacetic acid in a test tube at 100C for 1 h. Then cooled on ice. The absorbance at 520 nm was measured by spectrophotometer.

Total soluble sugar content in leaves was determined using anthrone reagent. Approximately 0.5 g dried leaf segments from different lines were ground into powder with liquid nitrogen, and then homogenized with 2 ml 80% ethanol in shaker at 200 rpm for 1 h. Following centrifuge at 6,000 × g for 10 min, and then collected as much supernatant as possible. Added equal volume of chloroform, completely mix, and then centrifuged at 12,000 × g for 10 min. The aqueous part was transferred to a new tube, 50 μl of each was mixed with 4.95 ml anthrone reagent and then boiled for 15 min. Measured the optical density of glucose standards at 620 nm by spectrophotometer.

For the content of MDA, about 1 g leaf segments from different lines were homogenized in 10 ml of 10% trichloroacetic (v/v) and centrifuged at 5,000 × g for 10 min. Following 2 ml of supernatant was reacted with 2 ml thiobarbituric acid in a test tube at 100C for 15 min, quickly cooled on ice, and the absorbance at 532 nm was measured by spectrophotometer.

The activities of SOD, POD, CAT and AsA content were detected using copper-zinc superoxide dismutase (CuZn-SOD) assay kit (A001-4–1), peroxidase assay kit (A084-3–1) catalase (CAT) assay kit (A007-1–1) and Vitamin C assay kit (A009-1–1) (http://www.njjcbio.com/, Nanjing Jiancheng Bioengineering Institute) according to the manufacturer’s manual.

Two-week-old seedlings of two sapk3 mutant lines and wide type plants withheld water for 7 days. According to the manual, leaf segments (0.5 g) from control and dehydration lines were homogenized using a chilled mortar.

For total SOD activity, 0.5 g leaf segments were ground into homogenate with 0.2 ml phosphate buffer solution (0.1 mol/L pH7-7.4), and then homogenized with 0.2 ml reagent No. 7. Centrifuged at 4,000 × g for 15 min. Then collected supernatant for Cu/Zn SOD determination. 2 ml of supernatant was reacted with 0.1 ml reagent No. 1, No. 2, No. 3 and No. 4 in a test tube at 37C for 40 min. Then reacted with 2 ml chromogenic agent for 10 min at room temperature. Measured the absorbance at 550 nm.

For POD activity, 0.5 g leaf segments were ground into homogenate with 0.2 ml phosphate buffer solution (0.1 mol/L pH7-7.4). Then collected the supernatant after centrifuging for 10 min (3500 × g). 0.1 ml of supernatant was reacted with 2.4 ml reagent No. 1, 0.3 ml reagent No. 2 and 0.2 ml reagent No. 3 in a test tube at 37C for 30 min. Then reacted with reagent No. 4 for 10 min at 37C. Collected the supernatant after centrifuging for 10 min (3500 × g). Measure the absorbance at 420 nm.

For CAT activity, 0.5 g leaf segments were ground into homogenate with 4.5 ml phosphate buffer solution (0.1 mol/L pH7-7.4). Then collected the supernatant after centrifuging for 10 min (2000 × g). 0.05 ml supernatant was reacted with 1 ml reagent No. 1 and 0.1 ml reagent No. 2 in a test tube at 37C for 1 min. Then mixed with 1 ml reagent No. 3 and 0.1 ml reagent No. 4. Measured the absorbance at 405 nm.

For the content of ascorbic acid (AsA), 0.5 g leaf segments were ground into homogenate with 5 mL of 5% (w/v) m-phosphoric acid. Following centrifugation at 10,000 × g for 15 min at 4 °C, the supernatant was used for the determination of AsA. Following 2 ml of supernatant was reacted with specified dosage of R1, R2, R3 and R4 (According to the manual) in a test tube at 22 ◦C for 30 min, quickly cooled on ice, and the absorbance at 536 nm was measured.

Stomatal number and stomatal conductance analysis

Leaves from two weeks wild-type plants and two sapk3 mutant lines were fixed with 2.5% glutaraldehyde and observed by scanning electron microscopy (S-3400 N, Hitachi, Japan). The second leaves of two weeks rice plants were selected for the stomatal conductance analysis, which was performed using the LI-6400XT Portable Photosynthesis System (LI-COR, USA) before and after 7 days exposure to drought stress.

Statistical analyses

Three independent experiments were conducted for drought stress analysis and gene expression analyses. Three biological replicates were performed for each experiment. For agronomic traits analyses, 25 plants of wild-type plants or two sapk3 mutant lines were randomly chosen to detect agronomic traits. Excel 2010 was used for making charts. All results are presented as means ± standard derivation (SD) of three biological replicates. Statistically significant difference analysis was conducted by “Student’s t-test” using SPSS statistics software.

Availability of data and materials

All data generated or analyzed during this study are included in this published article and its supplementary information files. The datasets and materials used or analysed during the current study available from the corresponding author on reasonable request.

References

  1. Gouda G, Gupta MK, Donde R, Mohapatra T, Vadde R, Behera L. Marker-assisted selection for grain number and yield-related traits of rice (Oryza sativa L.). Physiol Mol Biol Plants. 2020;26(5):885–98.

    Article  CAS  Google Scholar 

  2. Xing Y, Zhang Q. Genetic and molecular bases of rice yield. Annu Rev Plant Biol. 2010;61:421–42.

    Article  CAS  Google Scholar 

  3. Wang YH, Li JY. Branching in rice. Curr Opin Plant Biol. 2011;14(1):94–9.

    Article  CAS  Google Scholar 

  4. Li X, Qian Q, Fu Z, Wang Y, Xiong G, Zeng D, Wang X, Liu X, Teng S, Hiroshi F, Yuan M, Luo D, Han B, Li J. Control of tillering in rice. Nature. 2003;422:618–21.

    Article  CAS  Google Scholar 

  5. Koumoto T, Shimada H, Kusano H, She KC, Iwamoto M, Takano M. Ricemonoculm mutation moc2, which inhibits outgrowth of the second tillers, is ascribed to lack of a fructose-1, 6-bisphosphatase. Plant Biotech. 2013;30:47–56.

    Article  CAS  Google Scholar 

  6. Lin Q, Wang D, Dong H, Gu S, Cheng Z, Gong J, Wan J. Rice APC/CTE controls tillering by mediating the degradation of MONOCULM 1. Nat Commun. 2012;3:752.

    Article  Google Scholar 

  7. Lin Q, Wu F, Sheng P, Zhang Z, Zhang X, Guo X, Wan J. The SnRK2-APC/CTE regulatory module mediates the antagonistic action of gibberellic acid and abscisic acid pathways. Nat Commun. 2015;6:7981.

    Article  CAS  Google Scholar 

  8. Bahuguna RN, Gupta P, Bagri J, Singh D, Azri Kusuma D, Tao L. Forward and reverse genetics approaches for combined stress tolerance in rice. Ind J Plant Physiol. 2018;23:630–46.

    Article  Google Scholar 

  9. Yang PM, Huang QC, Qin GY, Zhao S, Zhou JG. Different drought stress responses in photosynthesis and reactive oxygen metabolism between autotetraploid and diploid rice. Photosynthetica. 2014;52(2):193–202.

    Article  CAS  Google Scholar 

  10. Ji K, Wang Y, Sun W, Lou Q, Mei H, Shen S, Chen H. Drought-responsive mechanisms in rice genotypes with contrasting drought tolerance during reproductive stage. J Plant Physiol. 2012;169(4):336–44.

    Article  CAS  Google Scholar 

  11. Shao GC, Deng S, Liu N, Yu SE, Wang MH, She DL. Effects of controlled irrigation and drainage on growth, grain yield and water use in paddy rice. Eur J Agron. 2014;53:1–9.

    Article  Google Scholar 

  12. Liu JX, Liao DQ, Oane R, Estenor L, Yang XE, Li ZC, Bennett J. Genetic variation in the sensitivity of anther dehiscence to drought stress in rice. Field Crops Res. 2006;97(1):87–100.

    Article  Google Scholar 

  13. Zhang J, Zhang S, Cheng M, Jiang H, Zhang X, Peng C, Jin J. Effect of drought on agronomic traits of rice and wheat: a meta-analysis. Int J Environ Res Public Health. 2018;15(5):839.

    Article  Google Scholar 

  14. Yoshida T, Mogami J, Yamaguchi-Shinozaki K. ABA-dependent and ABA-independent signaling in response to osmotic stress in plants. Curr Opin Plant Biol. 2014;21:133–9.

    Article  CAS  Google Scholar 

  15. Shinozawa A, Otake R, Takezawa D, Umezawa T, Komatsu K, Tanaka K, Amagai A, Ishikawa S, Hara Y, Kamisugi Y. SnRK2 protein kinases represent an ancient system in plants for adaptation to a terrestrial environment. Commun Biol. 2019;2:30.

    Article  Google Scholar 

  16. Kobayashi Y, Yamamoto S, Minami H, Kagaya Y, Hattori T. Differential activation of the rice sucrose nonfermenting1–related protein kinase2 family by hyperosmotic stress and abscisic acid. Plant Cell. 2004;16:1163–77.

    Article  CAS  Google Scholar 

  17. Takahashi F, Kuromori T, Sato H, Shinozaki K. Regulatory gene networks in drought stress responses and resistance in plants. Adv Exp Med Biol. 2018;1081:189–214.

    Article  CAS  Google Scholar 

  18. Umezawa T, Nakashima K, Miyakawa T, Kuromori T, Tanokura M, Shinozaki K, Yamaguchi-Shinozaki K. Molecular basis of the core regulatory network in ABA responses: sensing, signaling and transport. Plant Cell Physiol. 2010;51:1821–39.

    Article  CAS  Google Scholar 

  19. Fujii H, Chinnusamy V, Rodrigues A, Rubio S, Antoni R, Park SY, Cutler SR, Sheen J, Rodriguez PL, Zhu JK. In vitro reconstitution of an abscisic acid signalling pathway. Nature. 2009;462:660–4.

    Article  CAS  Google Scholar 

  20. Fujita Y, Nakashima K, Yoshida T, Katagiri T, Kidokoro S, Kanamori N, Umezawa T, Fujita M, Maruyama K, Ishiyama K. Three SnRK2 protein kinases are the main positive regulators of abscisic acid signaling in response to water stress in Arabidopsis. Plant Cell Physiol. 2009;50:2123–32.

    Article  CAS  Google Scholar 

  21. Nakashima K, Fujita Y, Kanamori N, Katagiri T, Umezawa T, Kidokoro S, Maruyama K, Yoshida T, Ishiyama K, Kobayashi M. Three Arabidopsis SnRK2 protein kinases, SRK2D/SnRK2.2, SRK2E/SnRK2.6/OST1 and SRK2I/SnRK2.3, involved in ABA signaling are essential for the control of seed development and dormancy. Plant Cell Physiol. 2009;50:1345–63.

    Article  CAS  Google Scholar 

  22. Umezawa T, Sugiyama N, Mizoguchi M, Hayashi S, Myouga F, Yamaguchi-Shinozaki K, Ishihama Y, Hirayama T, Shinozaki K. Type 2C protein phosphatases directly regulate abscisic acid-activated protein kinases in Arabidopsis. Proc Natl Acad Sci USA. 2009;106:17588–93.

    Article  CAS  Google Scholar 

  23. Liu J, Sun X, Liao W, Zhang J, Liang J, Xu W. Involvement of OsGF14b adaptation in the drought resistance of rice plants. Rice. 2019;12(1):82.

    Article  Google Scholar 

  24. Li X, Yu B, Wu Q, Min Q, Zeng R, Xie Z. Osmads23 phosphorylated by SAPK9 confers drought and salt tolerance by regulating aba biosynthesis in rice. PLoS Genet. 2021;17(8):e1009699.

    Article  CAS  Google Scholar 

  25. Wu Q, Liu Y, Xie Z, Yu B, Sun Y, Huang J. OsNAC016 regulates plant architecture and drought tolerance by interacting with the kinases GSK2 and SAPK8. Plant Physiol. 2022;189(3):1296–313.

    Article  CAS  Google Scholar 

  26. Wang P, Hsu CC, Du Y, Zhu P, Zhao C, Fu X, Zhang C, Paez JS, Macho AP, Tao WA. Mapping proteome-wide targets of protein kinases in plant stress responses. Proc Natl Acad Sci USA. 2020;117(6):3270–80.

    Article  CAS  Google Scholar 

  27. McLoughlin F, Galvan-Ampudia CS, Julkowska MM, Caarls L, Lauriere C. The Snf1-related protein kinases SnRK2.4 and SnRK2.10 are involved in maintenance of root system architecture during salt stress. Plant J. 2012;72:436–49.

    Article  CAS  Google Scholar 

  28. Julkowska MM, McLoughlin F, Galvan-Ampudia CS, Rankenberg JM, Kawa D, Klimecka M, Haring MA, Munnik T, Kooijman EE, Testerink C. Identification and functional characterization of the Arabidopsis Snf1-related protein kinase SnRK2.4 phosphatidic acid-binding domain. Plant Cell Environ. 2015;38:614–24.

    Article  CAS  Google Scholar 

  29. Chae MJ, Lee JS, Nam MH, Cho K, Hong JY, Yi SA. A rice dehydration-inducible SNF1-related protein kinase 2 phosphorylates an abscisic acid responsive element-binding factor and associates with ABA signaling. Plant Mol Biol. 2007;63:151–69.

    Article  CAS  Google Scholar 

  30. Chang Y, Nguyen BH, Xie Y, Xiao B, Tang N, Zhu W, Xiong L. Co-overexpression of the constitutively active form of OsbZIP46 and ABA-activated protein kinase SAPK6 improves drought and temperature stress resistance in rice. Front Plant Sci. 2017;8:1102.

    Article  Google Scholar 

  31. Shin R, Alvarez S, Burch AY, Jez JM, Schachtman DP. Phosphoproteomic identification oftaigets of the Arabidopsis sucrose nonfermenting-like kinase SnRK2.8 reveals a connection to metabolic processes. Proc Natl Acad Sci U S A. 2007;104(15):6460–5.

    Article  CAS  Google Scholar 

  32. Mizoguchi M, Umezawa T, Nakashima K, Kidokoro S, Takasaki H, Fujita Y. Two closely related subclassII SnRK2 protein kinases cooperatively regulate drought-inducible gene expression. Plant Cell Physiol. 2010;51:842–7.

    Article  CAS  Google Scholar 

  33. Tang N, Zhang H, Li X, Xiao J, Xiong L. Constitutive activation of transcription factor OsbZIP46 improves drought tolerance in rice. Plant Physiol. 2012;158:1755–68.

    Article  CAS  Google Scholar 

  34. Zong W, Tang N, Yang J, Peng L, Ma S, Xu Y. Feedback regulation of ABA signaling and biosynthesis by a bZIP transcription factor targets drought resistance related genes. Plant Physiol. 2016;171(4):2810–25.

    Article  CAS  Google Scholar 

  35. Verbruggen N, Hermans C. Proline accumulation in plants: a review. Amino Acids. 2008;35(4):753–9.

    Article  CAS  Google Scholar 

  36. Ouyang SQ, Liu YF, Liu P, Lei G, He SJ, Ma B, Zhang WK, Zhang JS, Chen SY. Receptor-like kinase OsSIK1 improves drought and salt stress tolerance in rice (Oryza sativa) plants. Plant J. 2010;62:316–29.

    Article  CAS  Google Scholar 

  37. Jambunathan N. Determination and detection of reactive oxygen species (ROS), lipid peroxidation, and electrolyte leakage in plants. Methods Mol Biol. 2010;639:292–8.

    Google Scholar 

  38. Wang Y, Hou Y, Qiu J, Wang H, Wang S, Tang L, Zhang J. Abscisic acid promotes jasmonic acid biosynthesis via a “SAPK10-bZIP72-AOC” pathway to synergistically inhibit seed germination in rice ( Oryza sativa ). New Phytol. 2020;228(4):1336–53.

    Article  CAS  Google Scholar 

  39. Zou M, Guan Y, Ren H, Zhang F, Chen F. A bZIP transcription factor, OsABI5, is involved in rice fertility and stress tolerance. PlantMol Biol. 2008;66:675–83.

    CAS  Google Scholar 

  40. Zhang X, Tang Y, Ma Q, Yang C, Mu Y, Suo H, Nian H. OsDREB2A, a rice transcription factor, significantly affects salt tolerance in transgenic soybean. PLoS ONE. 2013;8(12):e83011.

    Article  Google Scholar 

  41. Wang J, Lu K, Nie H, Zeng Q, Wu B, Qian J, Fang Z. Rice nitrate transporter OsNPF7.2 positively regulates tiller number and grain yield. Rice. 2018;11(1):12.

    Article  Google Scholar 

  42. Huang W, Nie H, Feng F, Wang J, Lu K, Fang Z. Altered expression of OsNPF7.1 and OsNPF7.4 differentially regulates tillering and grain yield in rice. Plant Sci. 2019;283:23–31.

    Article  CAS  Google Scholar 

  43. Lou DJ, Chen Z, Yu DQ, Yang XY. SAPK2 contributes to rice yield by modulating nitrogen metabolic processes under reproductive stage drought stress. Rice. 2020;13(1):35.

    Article  Google Scholar 

  44. Chun Y, Fang J, Zafar SA, Shang J, Zhao J, Yuan S, Li X. MINI SEED 2 (MIS2) Encodes a receptor-like kinase that controls grain size and shape in rice. Rice. 2020;13(1):7.

    Article  Google Scholar 

  45. Guo T, Chen K, Dong NQ, Shi CL, Ye WW, Gao JP, Lin HX. GRAIN SIZE AND NUMBER1 negatively regulates the OsMKKK10-OsMKK4-OsMPK6 cascade to coordinate the trade-off between grain number per panicle and grain size in rice. Plant Cell. 2018;30(4):871–88.

    Article  CAS  Google Scholar 

  46. Wang E, Wang J, Zhu X, Hao W, Wang L, Li Q, Zhang L, He W, Lu B, Lin H, Ma H, Zhang G, He Z. Control of rice grain-filling and yield by a gene with a potential signature of domestication. Nat Genet. 2008;40(11):1370–4.

    Article  CAS  Google Scholar 

  47. Liu X, Li Z, Hou Y, Wang Y, Wang H, Tong X, Zhang J. Protein interactomic analysis of SAPKs and ABA-Inducible bZIPs revealed key roles of SAPK10 in rice flowering. Int J Mol Sci. 2019;20(6):1427.

    Article  CAS  Google Scholar 

  48. Lou DJ, Wang HP, Yu DQ. The sucrose non-fermenting-1-related protein kinases SAPK1 and SAPK2 function collaboratively as positive regulators of salt stress tolerance in rice. BMC Plant Biol. 2018;18(1):203.

    Article  CAS  Google Scholar 

  49. Lou DJ, Wang HP, Liang G, Yu DQ. OsSAPK2 confers abscisic acid sensitivity and tolerance to drought stress in rice. Front Plant Sci. 2017;8:993.

    Article  Google Scholar 

  50. Vahisalu T, Puzorjova I, Brosche M, Valk E, Lepiku M, Moldau H, Pechter P, Wang YS, Lindgren O, Salojarvi J. Ozone-triggered rapid stomatal response involves the production of reactive oxygen species, and is controlled by SLAC1 and OST1. Plant J. 2010;62:442–53.

    Article  CAS  Google Scholar 

  51. Fan X, Wu J, Chen T, Tie W, Chen H, Zhou F, Lin Y. Loss-of-function mutation of rice SLAC7 decreases chloroplast stability and induces a photoprotection mechanism in rice. J Integr Plant Biol. 2015;57(12):1063–77.

    Article  CAS  Google Scholar 

  52. Sun SJ, Qi GN, Gao QF, Wang HQ, Yao FY, Hussain J, Wang YF. Protein kinase OsSAPK8 functions as an essential activator of S-type anion channel OsSLAC1, which is nitrate-selective in rice. Planta. 2016;243(2):489–500.

    Article  CAS  Google Scholar 

  53. Dey A, Samanta MK, Gayen S, Maiti MK. The sucrose non-fermenting 1-related kinase 2 gene SAPK9 improves drought tolerance and grain yield in rice by modulating cellular osmotic potential, stomatal closure and stress-responsive gene expression. BMC Plant Biol. 2016;16(1):158.

    Article  Google Scholar 

  54. Nguyen GN, Hailstones DL, Wilkes M, Sutton BG. DROUGHT STRESS: role of carbohydrate metabolism in drought-induced male sterility in rice anthers. J Agron Crop Sci. 2010;196(5):346–57.

    Article  CAS  Google Scholar 

  55. Oladosu Y, Rafii MY, Samuel C, Fatai A, Magaji U, Kareem I, Kolapo K. Drought resistance in rice from conventional to molecular breeding: a review. Int J Mol Sci. 2019;20(14):3519.

    Article  CAS  Google Scholar 

  56. Raza H, Khan MR, Zafar SA, Kirch HH, Bartles D. Aldehyde dehydrogenase 3I1 gene is recruited in conferring multiple abiotic stress tolerance in plants. Plant Biol. 2022;24(1):85–94.

    Article  CAS  Google Scholar 

  57. Anwar A, Wang K, Wang J, Shi L, Du L, Ye X. Expression of arabidopsis ornithine aminotransferase (AtOAT) encoded gene enhances multiple abiotic stress tolerances in wheat. Plant Cell Rep. 2021;40(7):1155–70.

    Article  CAS  Google Scholar 

  58. Zafar SA, Hameed A, Ashraf M, Khan AS, Qamar Z, Li X, Siddique KHM. Agronomic, physiological and molecular characterisation of rice mutants revealed the key role of reactive oxygen species and catalase in high-temperature stress tolerance. Funct Plant Biol. 2020;47(5):440–53.

    Article  CAS  Google Scholar 

  59. Ahmad H, Zafar SA, Naeem MK, Shokat S, Inam S, Naveed AS, Xu JL, Li ZK, Ali GM, Khan MR. Impact of pre-anthesis drought stress on physiology, yield-related traits and drought responsive genes in green super rice. Front Genet. 2022;13:832542.

    Article  CAS  Google Scholar 

  60. Sahebi M, Hanafi MM, Rafii MY, Mahmud TMM, Azizi P, Osman M, Miah G. Improvement of drought tolerance in rice (Oryza sativa L.): Genetics, genomic tools, and the WRKY gene family. BioMed Res Int. 2018;2018:3158474.

    Article  Google Scholar 

  61. Mazur R, Maszkowska J, Anielska-Mazur A, Garstka M, Polkowska-Kowalczyk L, Czajkowska A. The SnRK21.0 kinase mitigates the adverse effects of salinity by protecting photosynthetic machinery. Plant Physiol. 2021;187(4):2785–802.

    Article  CAS  Google Scholar 

  62. Liu W, Xie X, Ma X, Li J, Chen J, Liu YG. DSDecode: A web-based tool for decoding of sequencing chromatograms for genotyping of targeted mutations. Mol Plant. 2015;8(9):1431–3.

    Article  CAS  Google Scholar 

  63. Huai J, Wang M, He J, Zheng J, Dong Z, Lv H, Wang G. Cloning and characterization of the SnRK2 gene family from Zea mays. Plant Cell Rep. 2008;27(12):1861–8.

    Article  CAS  Google Scholar 

  64. Chen S, Tao L, Zeng L, Vega-Sanchez ME, Umemura K, Wang GL. A highly efficient transient protoplast system for analyzing defence gene expression and protein-protein interactions in rice. Mol Plant Pathol. 2006;7(5):417–27.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This research was supported by National Natural Science Foundation of China (31960371, 21562029).

Author information

Authors and Affiliations

Authors

Contributions

DJL., XYY. designed the experiments, DJL., SPL., YL. performed the experiments and wrote the manuscript, ZC., DQY., XYY. analyzed the data and edited the article. Both authors read and approved the final article.

Corresponding author

Correspondence to Xiaoyan Yang.

Ethics declarations

Ethics approval and consent to participate

Experimental research and field studies on plants (either cultivated or wild) in this paper, including the collection of plant material, comply with relevant institutional, national, and international guidelines and legislation. Wild-type rice plants (Oryza sativa L.) used in this paper, were from the laboratory of XYY, Faculty of Life Science and Technology, Kunming University of Science and Technology. Permissions had been obtained.

Consent for publication

Not applicable.

Competing interests

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Supplementary Table 1. 

Amino acid sequences of selectedSnRK2s inthe phylogenetic tree.

Additional file 2: Supplementary Table 2.

Nucleotide sequence of SAPK3 CDS and predicted amino acid sequences in different genotypes.

Additional file 3: Supplementary Table 3.

Primers and oligos used in this study.

Additional file 4: Supplementary Figure 1. 

Agronomic traits of sapk3mutant lines. (A) Panicle phenotypes of wild-type plants and two sapk3 mutant lines. (B) Grain width (C) Panicle length (D) Grain number per panicle (E) Setting rate per panicle (F) 1000-seed weight of wild-type plants and two sapk3 mutant lines. Values represent the means ± SD of three biological replicates (25 plants for each replicate).

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Lou, D., Lu, S., Chen, Z. et al. Molecular characterization reveals that OsSAPK3 improves drought tolerance and grain yield in rice. BMC Plant Biol 23, 53 (2023). https://doi.org/10.1186/s12870-023-04071-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12870-023-04071-8

Keywords