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 Biology volume 16, Article number: 158 (2016)
Family members of sucrose non-fermenting 1-related kinase 2 (SnRK2), being plant-specific serine/threonine protein kinases, constitute the central core of abscisic acid (ABA)-dependent and ABA-independent signaling pathways, and are key regulators of abiotic stress adaptation in plants. We report here the functional characterization of SAPK9 gene, one of the 10 SnRK2s of rice, through developing gain-of-function and loss-of-function phenotypes by transgenesis.
The gene expression profiling revealed that the abundance of single gene-derived SAPK9 transcript was significantly higher in drought-tolerant rice genotypes than the drought-sensitive ones, and its expression was comparatively greater in reproductive stage than the vegetative stage. The highest expression of SAPK9 gene in drought-tolerant Oryza rufipogon prompted us to clone and characterise the CDS of this allele in details. The SAPK9 transcript expression was found to be highest in leaf and upregulated during drought stress and ABA treatment. In silico homology modelling of SAPK9 with Arabidopsis OST1 protein showed the bilobal kinase fold structure of SAPK9, which upon bacterial expression was able to phosphorylate itself, histone III and OsbZIP23 as substrates in vitro. Transgenic overexpression (OE) of SAPK9 CDS from O. rufipogon in a drought-sensitive indica rice genotype exhibited significantly improved drought tolerance in comparison to transgenic silencing (RNAi) lines and non-transgenic (NT) plants. In contrast to RNAi and NT plants, the enhanced drought tolerance of OE lines was concurrently supported by the upgraded physiological indices with respect to water retention capacity, soluble sugar and proline content, stomatal closure, membrane stability, and cellular detoxification. Upregulated transcript expressions of six ABA-dependent stress-responsive genes and increased sensitivity to exogenous ABA of OE lines indicate that the SAPK9 is a positive regulator of ABA-mediated stress signaling pathways in rice. The yield-related traits of OE lines were augmented significantly, which resulted from the highest percentage of fertile pollens in OE lines when compared with RNAi and NT plants.
The present study establishes the functional role of SAPK9 as transactivating kinase and potential transcriptional activator in drought stress adaptation of rice plant. The SAPK9 gene has potential usefulness in transgenic breeding for improving drought tolerance and grain yield in crop plants.
Plants, being sessile are incessantly confronted by different environmental stresses, which include drought, high salinity, and extreme temperature, affecting both biomass productivity and grain yield of crops. To cope with such adverse multiple stresses, diverse molecular and physiological mechanisms have been evolved by the plant kingdom in general and a plant species in particular. Therefore, it is indispensable to comprehend the specific mechanism of responses by any crop species to such stresses, with the ultimate aim of improving crop performance. Rice (Oryza sativa L.), being one of the most important cereal crops, feeds more than half of the world population and is adversely affected by drought at the morphological, physiological and molecular level. The phytohormone abscisic acid (ABA), which is produced under the drought stress functions to regulate several developmental and physiological processes including seed maturation, germination, seedling growth and transpiration. The ABA level increases under water deficit condition in plants triggering stomatal closure and responses to stress tolerance . Stress signals are recognized by specialized signaling pathways which transmit them to different cellular compartments, and the numerous evidences demonstrated that protein kinases play vital roles in the responses to such environmental stimuli . The sucrose non-fermenting 1-related kinase 2 (SnRK2) family members, which function in diverse developmental processes in plants [3, 4], have been shown to be the positive regulators of plant response to abiotic stresses [5–8]. The members of SnRK2 family work at the merging point of the ABA-dependent and ABA-independent stress signaling pathways.
In Arabidopsis, ten members of the SnRK2 family have been identified and divided into three subclasses based on amino acid sequence similarity . The subclass I comprises of kinases that are not activated by ABA, subclass II are either not activated or very weakly activated by ABA (depending on plant species). However, the subclass III is strongly activated by ABA. The amino acid sequences of all the SnRK2s can be divided into two regions, the highly conserved N-terminal kinase domain, and the C-terminal regulatory domain. The C-terminal domain contains stretches of acidic amino acids, either glutamic acid (E) (subclass I) or aspartic acid (D) (subclass II and III). Further, the C-terminal regulatory domain consists of two subdomains, subdomain I and subdomain II. The subdomain I which is required for activation by osmotic stress, independent of ABA, is present in all SnRK2 family members; whereas subdomain II is necessary for the ABA response and is specific to the ABA-dependent SnRK2s only. All the SnRK2 members, except SnRK2.9, were found to be rapidly induced by different osmolyte treatments, such as mannitol, sorbitol, sucrose or NaCl and few of them also by ABA . The Arabidopsis OST1/SnRK2.6/SRK2E protein kinase has been shown to regulate the ABA-mediated stomatal closure and act upstream of reactive oxygen species production . However, the srk2e (ost1) mutant in Arabidopsis has been found to be defective in ABA-induced stomatal closure and it showed a wilty phenotype . The srk2d/e/i triple mutant (for SnRK2.2/SRK2D, SnRK2.6/SRK2E and SnRK2.3/SRK2I genes) of Arabidopsis displayed dramatically decreased drought tolerance and extreme insensitivity to ABA, as documented by defects in seed germination and seedling growth as well as decreased expression of ABA- and stress-inducible genes . The knockout of these three genes, which belong to the subclass III of SnRK2 family, therefore, almost completely blocks ABA responses, demonstrating that they are the essential components of ABA-stress signaling pathway in Arabidopsis .
Eleven SnRK2 members designated as ZmSnRK2s have been identified in maize, and most of them are inducible by one or more abiotic stresses . The ZmSnRK2.8 protein, which is highly homologous to Arabidopsis OST1, has been found to be involved in diverse stress signaling pathways, particularly in salt stress tolerance [12, 13]. In wheat, the first characterized SnRK2 member, PKABA1 was reported to be induced by hyperosmotic stress, ABA, and multiple other environmental factors [14, 15]. Afterwards, three more wheat SnRK2 members viz., TaSnRK2.3, TaSnRK2.4 and TaSnRK2.8 have been characterized and found to be involved in development and abiotic stress tolerance [16–18]. Thus, substantial evidence showed that the SnRK2 protein kinases are involved in multiple environmental stress responses and all have potential biotechnological utility for generation of high yielding abiotic stress tolerant crops . In vitro studies have documented that the ABA-activated SnRK2s phosphorylate the downstream target proteins in different plant species; and this phosphorylation is needed for the transcriptional activity of the individual target proteins, which in turn induce the expression of hierarchically downstream genes to mitigate the stress condition. This includes bZIP transcription factors, such as TRAB1 from rice , TaABF from wheat , and AREB1 from Arabidopsis , or RNA-binding VfAKIP1 proteins from Vicia faba .
In rice, the ten SnRK2 members, designated as SAPK1-10 (osmotic stress/ABA-activated protein kinase 1-10) have been identified. All of them were found to be activated by hyperosmotic stress, and SAPK8-10 were also induced by ABA . The overexpression of SAPK4 showed improved salt tolerance in rice . In the domain swapping experiments using rice SAPKs, it has been observed that the grafting of the non-catalytic C-terminal region from SAPK8 (Glu-254 to Met-372), which contains the domain involved in ABA-responsive activation, onto the catalytic domain of SAPK2, which doesn’t contain the domain for ABA-responsive activation, is sufficient to confer ABA responsiveness in SAPK2 . None of the three ABA-dependent and osmotic stress-activated SAPKs, i.e. SAPK8, SAPK9, and SAPK10, belonging to subclass III of SnRK2 family, has been functionally characterized by transgenic approach till date. Therefore, understanding the molecular basis of subclass III SAPK gene function in rice is necessary for the development of drought-tolerant transgenic rice or to design SAPK gene-based marker for molecular breeding of drought-tolerant rice cultivar. We were interested in the detailed characterization of the SAPK9 gene and its functional role in drought tolerance through (i) expression profiling of the gene in selected drought-sensitive and drought-tolerant rice genotypes, (ii) determining the allelic polymorphism in the coding DNA sequence (CDS) of the gene in these genotypes, (iii) cloning the highly expressed SAPK9 CDS from the drought-tolerant rice genotype Oryza rufipogon followed by assaying the kinase activity of the recombinant protein, and (iv) developing transgenic rice lines for overexpression and endogenous silencing of this particular gene. Transgenic experiments demonstrated that in contrast to silencing, the overexpression of SAPK9 increased drought tolerance and grain yield by adjusting osmotic potential and stomatal closure, thereby decreasing the cellular membrane damage and reactive oxygen species activity. In addition, the increased and decreased ABA sensitivity and ABA-responsive gene expression in transgenically overexpressed and silenced lines, respectively, confirmed that the SAPK9 is a positive regulator of ABA-dependent stress-responsive signaling pathway in rice.
Plant materials, growth condition and stress treatments
Nine cultivated indica rice (Oryza sativa L.) genotypes and two wild rice progenitors were selected for the experimental works. The cultivated rice genotypes comprised of drought-tolerant- Manipuri, Nagina22, Vandana and drought-sensitive- Swarna, HRC300, IR20, IR36, IR64, and IR72; while the two wild rice genotypes - O. rufipogon and O. nivara are drought-tolerant. The seed samples of all rice genotypes were provided by Indian Agricultural Research Institute, New Delhi. All the experimental works were performed using the seeds of same harvest and storage conditions, and rice genotypes were usually grown for 30 days inside the glasshouse, which was maintained at 25/28 °C temperature with 70 % relative humidity and 16/8 h light/dark photoperiod. Thereafter, the potted plants were transferred into the net-house. To analyse the effects of drought stress and ABA on the expression of SAPK9 gene, rice seedlings were subjected to dehydration and ABA (100 μM) treatment for a period of 6–48 h and an untreated sample at 0 h was used as control. The transcript level of SAPK9 gene was analysed from leaf tissues of rice genotypes grown under: (i) normal growth condition- designated as before stress (BS), (ii) withholding water supply for 8 days- referred as after stress (AS), and (iii) resuming water supply for 3 days after the drought stress period- referred as after recovery (AR). The experiments were performed at both vegetative and reproductive (or grain filling) stages. For T1 transgenic lines and non-transgenic (NT) plants drought stress treatment at the vegetative and early reproductive (panicle initiation) stages, and analysis of grain yield under drought condition were performed following the reported protocol .
Cellular RNA isolation, first strand cDNA synthesis and real‑time PCR analysis
The total cellular RNA was isolated from rice leaves using RNeasy Mini Kit (Qiagen) following the vendor’s protocol. Total RNA (2 μg) was used to synthesize the first strand cDNA with the help of gene-specific 3'-primer and Transcriptor First Strand cDNA synthesis kit version 6.0 (Roche Diagnostics India Pvt. Ltd.) following the instructions of the manufacturer. For transcript expression profiling, real-time PCR was performed by means of SYBR green-based relative quantification method using RealMasterMix SYBR ROX (5 PRIME) kit in Eppendorf Realplex2 Master Cycler. The relative gene expression level was determined following the reported literature . For each sample, three replicates were taken. In each case, rice polyubiquitin1 gene (OsUbi1) was used as internal reference. The different primers used for the real-time PCR are provided in Additional file 1: Table S1.
Standard protocols were followed to isolate the genomic DNA from rice leaves and to perform Southern blot hybridization . In brief, the HindIII-digested genomic DNA (12 μg) upon electrophoresis on agarose gel overnight, was transferred onto nylon membrane (Hybond-N+). A 470 bp DNA fragment from the middle of SAPK9 coding DNA sequence (CDS) from O. rufipogon was PCR-amplified using gene-specific primers SA9SF-SA9SR (Additional file 1: Table S1) and radiolabelled with P32-dCTP (3500 Ci/mmol) by random priming using rediprime II DNA labelling system (GE Healthcare, USA) following the vendor’s instructions. For autoradiography, the Cylone® Plus phosphor system (Perkin Elmer) was used to scan the multi-sensitive X-ray film (Perkin Elmer).
Cloning and sequence analysis of the SAPK9 CDS
The full-length CDS of SAPK9 gene was PCR amplified using 1st strand cDNA sample as a template from the drought-sensitive and drought-tolerant rice genotypes mentioned in plant materials section. The gene-specific primers SK9F and SK9R (Additional file 1: Table S1) were used for the PCR amplification in a thermocycler (Applied Biosystems) with the following thermal profile: an initial denaturation at 98 °C for 2 min, followed by 30 cycles of 98 °C-for 10 s, 62 °C for 15 s, 72 °C for 1 min and a final extension at 72 °C for 8 min. All the PCR amplicons (~1086 bp size) were individually cloned by blunt-end ligation into pUC18 plasmid digested with SmaI restriction enzyme. Positive clones were sequenced, and nucleotide polymorphism of all the CDSs were analysed. The multiple sequence alignment of the CDS-derived polypeptides obtained from all the genotypes was performed using the Jalview 2 software (http://www.jalview.org/Download).
Multiple amino acid sequence alignment of SnRK2 family proteins from Arabidopsis, rice, and a few other plant species were performed using the ClustalX2 (http://www.clusal.org/clustal2/). The phylogenetic tree was constructed using the neighbor-joining method in MEGA6 (http://www.megasoftware.net/mega.php). To evaluate the reliability of the tree, a bootstrap analysis was performed using 1000 replicates in MEGA6.
The homology model of SAPK9 was made using the MODELLER 9.15 (https://salilab.org/modeller/download_installation.html). The solved crystal structure of the OST1/SnRK2.6 protein from Arabidopsis (PDB ID: 3UC4) was taken as a template for constructing the in silico model of SAPK9. The resulting model was energy minimized using Insight II (2000.1, Accelrys Inc.) followed by stereo-chemical evaluation using MolProbity (http://molprobity.biochem.duke.edu/). The pictorial representations were prepared in PyMol (http://www.pymol.org).
In vitro kinase assay
The recombinant pRSET plasmid (Invitrogen) containing the 1086 bp CDS of SAPK9 in XhoI and EcoRI restriction sites allowed expression of 6xHis N-terminal tagged SAPK9 protein in E. coli BL21 (DE3) pLysE strain (Invitrogen) upon 1 mM IPTG induction. Similarly, first 900 bp (encoding 300-amino acid) from the OsbZIP23 CDS of O. rufipogon  was cloned into a pRSET vector in BamHI and EcoRI restriction sites, transformed and expressed in pLysE cells with 0.5 mM IPTG induction. The expressed proteins were purified in native condition and used for in vitro kinase assay. In vitro phosphorylation of generic substrate histone III (Sigma) was performed as described previously . In vitro phosphorylation of OsbZIP23 was performed by incubating the individual reaction mixture for 5, 25 and 40 min at 25 °C following the above-mentioned protocol. The products were fractionated in 12 % SDS-PAGE and visualized by autoradiography.
Preparation of constructs for overexpression and silencing of SAPK9 gene, and development of transgenic rice plants
For the preparation of SAPK9 overexpression (OE) construct, the full-length 1086 bp CDS of SAPK9 was cloned (as SalI-KpnI fragment) from O. rufipogon under the OsUbi1 promoter (as HindIII-SalI fragment) and NOS transcription terminator (as KpnI-EcoRI fragment) in the pCAMBIA1301 binary plasmid. Thus, the SAPK9 OE genetic construct was prepared (Additional file 2: Figure S5a). On the other hand, the RNAi-mediated gene silencing (RNAi) construct of SAPK9 was prepared (Additional file 2: Figure S5b) by cloning a 605 bp DNA fragment from 5′-part of SAPK9 CDS in sense (as SalI-BamHI fragment) and antisense (as BamHI-KpnI fragment) direction flanked by an arbitrary linker DNA (having BamHI site at both ends) under the OsUbi1 promoter (as HindIII-SalI fragment) and NOS transcription terminator (as KpnI-EcoRI fragment) in pCAMBIA1301 binary plasmid. Drought-sensitive IR20 rice cultivar was genetically transformed separately with these OE and RNAi constructs by Agrobacterium-mediated transfer technique, as reported earlier . The transgenic rice lines developed with the OE and RNAi constructs were designated as SAOE#1, 2, 3 etc. and RNAi#1, 2, 3 etc, respectively. The list of primers used for preparing genetic constructs is given in Additional file 1: Table S1.
Western blot analysis of SAPK9 protein in transgenic lines and non-transgenic rice plants
Protein samples (40 μg each) isolated from the leaf tissues of transgenic lines and NT rice plants were used for western blotting experiment as per the method reported earlier . The custom made 15-amino acid peptide (MERNAAGPLGMEMNC) from the N-terminal end of SAPK9 was used to raise the polyclonal antibody in rabbit (IMGENEX India). Affinity purified polyclonal antibody of SAPK9 was taken as primary antibody (with 1: 1000 dilutions). The monoclonal antibody of plant β-actin (Sigma, catalog no.-A0480) was used (with 1:500 dilutions) to detect actin protein in the samples as a loading control. Immuno-detection was performed using the Lumi-LightPLUS Western Blotting Kit (Roche Diagnostics India Pvt. Ltd.), as per the vendor’s protocol.
Estimation of water loss rate, relative water content, proline content and soluble sugar content in leaf samples of transgenic and non-transgenic rice plants
Water loss rate (WLR) and relative water content (RWC) in the leaves of transgenic and NT rice plants were determined on the basis of the published protocols [30, 31]. From the leaves (50 mg) of transgenic lines and NT rice plants free proline content was estimated as per the published literature . Similarly, total soluble sugar content in rice leaves (100 mg) was determined using anthrone reagent . Consistent results were obtained from two sets of the experiment, and the result from one experimental set is documented here. Average of three replicates was calculated to represent each data point.
Analysis of rice stomata by scanning electron microscopy
Scanning electron microscopy analysis of leaf samples from transgenic and NT rice plants (collected before and after drought stress) were performed following the described method . The samples were mounted onto the sample stages and coating treatment was performed. The stomatal pictures were obtained using a scanning electron microscope (ZEISS, Germany).
Estimation of malondialdehyde (MDA) content and detection of reactive oxygen species (ROS) in leaf samples of rice plants
Rice leaf (100 mg) of transgenic lines and NT plants was used to determine the MDA content following the reported method . In vivo localization of ROS in intact leaves of the transgenic lines and NT plants were carried out following the method described earlier . Consistent results were obtained from two sets of experiment and the result from one experimental set is documented here. Average of three replicates was calculated to represent each data point for MDA content.
Subcellular localization of the SAPK9 protein
The SAPK9 CDS was cloned in frame at the 5′-end of the green fluorescence protein (GFP) gene at BglII and SpeI restriction sites under the control of constitutive CaMV35S promoter in the binary plasmid pCAMBIA1302. The list of primers used is given in Additional file 1: Table S1. Agrobacterium-mediated transformation of onion epidermal cells was carried out following the previously published protocol .
Assay of ABA sensitivity during seed germination and post-germination stages
Sensitivity towards ABA during seed germination and post-germination (seedling) stages was analysed using the different concentrations of ABA (0, 1, 3, 6 μM) in the growth media following the reported protocol .
Estimation of spikelet fertility and pollen maturation
The transgenic lines and NT plants were grown normally in PVC pipes. Upon attaining flowering stage, plants were subjected to drought stress for a period of 10 days followed by recovery until seed maturation occurred. The panicle weight and spikelet fertility were calculated in replicates. For pollen staining, 1 % I2-KI solution based method was followed as described previously .
Sequence information of the genes mentioned in this study is available in the GenBank database having their respective accession numbers: Rice polyubiquitin 1 (OsUbi1) gene (AF184279), OsUbi1 promoter (AY785814), SAPK9 CDS from wild rice O. rufipogon (KT387673), OsRab16B gene (AF333275), OsRab21 gene (Y00842) OsLEA3-1 gene (DQ789359), OsbZIP23 gene (KP779638), TRAB1 gene (BAD09357), OsbZIP46 gene (ADK60888), SAPK8 gene (AB125309), SAPK10 gene (AB125311), OsSLAC1 gene (LOC_Os04g48530.1) and OsSLAC7 gene (LOC_Os01g28840).
Results and discussion
The single copy SAPK9 gene displays differential transcript expression profile and allelic polymorphism in drought-tolerant and drought-sensitive rice genotypes
The ability of a plant species to withstand drought stress condition is a genotype-dependent functional activity of the particular genes involved in the stress-adaptive pathways. Hypothesizing that the SAPK9 is a key player in drought tolerance, the relative expression level of SAPK9 gene was analysed in nine selected drought-sensitive and drought-tolerant indica rice genotypes along with two wild rice progenitors during the vegetative and reproductive stages of growth through real-time PCR. For the analysis, first strand cDNA was synthesized from leaf total RNA of all the genotypes under the condition of before stress (BS), after stress (AS) and after recovery (AR). The real-time PCR analysis of BS, AS and AR samples revealed that the expression level of SAPK9 transcript in O. rufipogon is substantially greater than the other drought-tolerant genotypes viz., O. nivara, Nagina22, Manipuri and Vandana in both vegetative stage (Fig. 1a) and reproductive stage (Fig. 1b). On the other hand, the SAPK9 expressions in drought-sensitive genotypes viz., HRC300, IR20, IR36, IR64, IR72, and Swarna were comparatively much lower than the drought-tolerant genotypes in both the growth stages (Fig. 1a and b). It was also observed that the expression of SAPK9 transcript is elevated in reproductive stage (Fig. 1b) compared to the vegetative stage (Fig. 1a). To check the influence of gene copy number if any, on the differential expression of SAPK9 transcript in the rice genotypes under the same growth conditions, we performed Southern hybridization using SAPK9 gene-specific probe. A single band observed in Southern blot indicated the presence of one copy of the endogenous SAPK9 gene in all the rice genotypes investigated (Fig. 1c). The occurrence of a single copy endogenous SAPK9 gene in the selected rice genotypes implies that the gene copy number does not influence the differential expression of SAPK9 gene in drought-sensitive and drought-tolerant rice genotypes.
Single nucleotide polymorphism (SNP) in the CDS is of great importance in the crop improvement programme as the functions of several genes are known to be modulated by the associated SNPs leading to the differences in plant performance. Assessing genetic diversity of candidate gene sequences involved in stress-responsive pathways lead to the identification of specific alleles which are linked to particular agronomic traits of abiotic stress responses. Trait enhancing superior alleles, which are dominant and express better in wild species than cultivated varieties and are left behind during evolution and domestication of present-day commercial cultivars, can be transferred to elite genetic backgrounds for improvement of desirable traits of a particular crop species . Therefore, it was of interest to investigate the presence of naturally occurring polymorphisms in the CDS of SAPK9 gene in all the chosen rice genotypes.
For the analysis of natural allelic polymorphism, the full-length CDS of SAPK9 gene was amplified through RT-PCR from total RNA of leaf tissues using two gene-specific primers. All the genotypes showed positive PCR amplification with an expected product size of ~1086 bp (Additional file 3: Figure S1). The sequence of each PCR amplicon was determined after cloning them individually. Three non-synonymous SNPs were detected among the SAPK9 CDSs obtained from eleven rice genotypes upon multiple sequence alignment. The first non-synonymous SNP was found in the CDS of Swarna, creating the nucleotide modification (T to A) at the 593rd position from the start codon, which led to the 198th amino acid change from isoleucine (I) to asparagine (N) (Additional file 4: Figure S2). The resulted amino acid change was located in the αF-helix of SAPK9 protein (Fig. 1d). The 2nd non-synonymous SNP was noticed in O. rufipogon, causing the nucleotide modification (C to T) at the 746th position from the start codon with the 249th amino acid change from proline (P) to leucine (L) (Additional file 4: Figure S2), and this was observed in the αH-helix of the SAPK9 protein (Fig. 1d). The 3rd non-synonymous SNP, being observed in a few genotypes, resulted in the nucleotide modification (A to C) at the 970th position from the start codon, which changed the 324th amino acid from serine (S) to arginine (R). Serine is present at 324th amino acid position in Swarna, IR20, IR36, IR64, HRC300, and Vandana; whereas arginine occurs at the same position in O. rufipogon, O. nivara, IR72, Nagina22 and Manipuri (Additional file 4: Figure S2). This last amino acid change was found to be located in between the SnRK2 box and ABA-dependent activation box (ABA-box) of SAPK9 protein (Fig. 1d). Various reports of non-synonymous SNPs leading to amino acid substitution in the CDS of drought responsive genes and their correlation with the particular trait have been published earlier . However, the present study showed that the natural occurrence of three non-synonymous SNPs did not affect any important functional domains and also could not be correlated with the differential expression level of the SAPK9 gene in the chosen drought-sensitive and drought-tolerant rice genotypes. Therefore, the highly expressed SAPK9 CDS of wild progenitor O. rufipogon was selected for further functional analyses.
The SAPK9 gene from drought-tolerant O. rufipogon encodes a protein belonging to subclass III of SnRK2 family, exhibits higher expression in leaf tissues and responses to both drought stress and ABA treatment
The full-length CDS of SAPK9 gene was cloned after isolation from the drought-tolerant rice genotype O. rufipogon by RT-PCR (Fig. 2a). Comparison of the cloned gene sequence (accession no. KT387673) with a previously reported sequence (accession no. AB125310) of SAPK9 gene from a japonica rice cultivar revealed alternation of the nucleotide at two positions. A nucleotide modification (C to T) at the 746th position from the start codon created the 249th amino acid change from proline (P) to leucine (L); and another nucleotide modification (A to C) at the 970th position from the start codon altered the 324th amino acid from serine (S) to arginine (R) (Additional file 5: Figure S3). The phylogenetic analysis of rice SnRK2 family members (designated as SAPKs) with Arabidopsis SnRK2s and related proteins from other plants revealed that the ten SAPK proteins of rice are clustered into three subclasses (Fig. 2b). Each subclass contains members from both Arabidopsis and rice SnRK2, demonstrating that the distinction of these subclasses originated before the divergence of dicots and monocots. The SAPK9 was found to be clustered into subclass III of SnRK2 family along with the well characterized ABA-activated SnRK2 genes, i.e. AAPK of Vicia faba  and OST1 (SnRK2.6/SRK2E) of Arabidopsis [6, 7], which have been reported to be involved in the ABA-regulated stomatal closure. Apart from these two genes, the SAPK9 was observed to be grouped with SnRK2.2 and SnRK2.3 of Arabidopsis , and SAPK8 and SAPK10 of rice  in the subclass III of SnRK2 family. Growing evidence indicate that subclass III SnRK2s are the global regulators of multiple stress signaling pathways, and these proteins have certain characteristic functional domains (Fig. 1d).
The tissue specificity of SAPK9 gene was monitored through real-time PCR analyses using the total RNA sample isolated from different tissues, namely root, shoot, stem, leaf, leaf sheath and panicle. The strongest expression of SAPK9 transcript was found in the leaf followed by root (Fig. 2c). Further, to analyse the transcriptional responses of SAPK9 gene under drought stress and ABA treatment, real-time PCR analysis was performed. The SAPK9 expression was observed to be significantly induced by exogenous ABA within 6 h of treatment in comparison to drought stress. However, the SAPK9 expression upon drought stress increased significantly after 24 h (Fig. 2d). The results indicate that the SAPK9 is involved in responses to both drought stress and ABA treatment, confirming the subclass III nature of this SnRK2 member .
The SAPK9 protein has a characteristic kinase fold structure and the recombinant protein possesses autophosphorylation and transphosphorylation activities
Analysis of the amino acid sequence (362 amino acids) derived from the cloned SAPK9 CDS revealed that the predicted molecular mass of SAPK9 is ∼ 40 kDa. The SAPK9 shares 80 and 76 % sequence similarity with the two reported Snf1-related kinase proteins viz., SnRK2.6 and SnRK2.3 of Arabidopsis thaliana, respectively (Additional file 6: Figure S4). Using Arabidopsis SnRK2.6 (PDB ID: 3UC4) as a template, the homology modelling revealed the bilobal (N-lobe and C-lobe) kinase fold structure of SAPK9 (Fig. 3a). The quality of the model displayed 94.5 % residues in the core region and 2.5 % residues in the allowed region. Interestingly, the in silico structure of catalytic SAPK9 displayed the canonical Ser/Thr kinase fold, identical to other Snf1 kinase domains. The well-ordered characteristic SnRK2 box (Fig. 3a) forming a single α-helix was found to be present in the N-terminal lobe packed parallel against the αC helix. In the N-terminal lobe, five β-sheets plus one α-helix called helix αC was found. The C-terminal lobe was found to be larger and mainly helical (Fig. 3a). The activation loop was found to be present at amino acid position 161 to 187. In the activation loop, S176 and T177 were predicted to be the phosphorylation site for SAPK9 activation (Fig. 3a). The catalytic cleft resides at the junction of the larger C-terminal lobe, and this cleft contains the substrate-binding site, the ATP binding G-loop and the Mg2+ binding site (Fig. 3b). It has been established before that the αC-helix in the N-terminal lobe is required for correct folding of the catalytic centre and kinase activation loop, and it is stabilized by the SnRK2 box (Fig. 3a and b) .
To study the kinase activity of recombinant SAPK9 in vitro, the bacterially expressed 6xHis-tagged SAPK9 protein was purified by Ni-NTA chromatography under native condition (Fig. 3c) and kinase assay was performed. It was observed that the generic substrate histone III was phosphorylated by the recombinant kinase (Fig. 3d). Moreover, the occurrence of an extra band corresponding to the size of recombinant SAPK9 (~42 kDa) suggested that it is able to autophosphorylate. This was confirmed by performing the kinase assay in the absence of substrate (Fig. 3d). These results are in accordance with the previous reports on the kinase activity of subclass III SnRK2s, OST1/SnRK2.6 and ZmSnRK2.8 [43, 44]. It is known that the AREB/ABF proteins require the phosphorylation of their multiple conserved sites by SnRK2 protein kinases for ABA-dependent activation [21, 41, 45]. Previous studies have shown that the SAPK9 can phosphorylate the downstream AREB/ABF type transcription factors such as TRAB1 and OsbZIP46 [46, 47]. As OsbZIP23 is phylogenetically very close to OsbZIP46, therefore we investigated the potential phosphorylation activity of SAPK9 on OsbZIP23 substrate. For this, the N-terminal 6xHis-tagged OsbZIP23 was expressed in the bacterial system and the recombinant OsbZIP23 protein was purified by Ni-NTA chromatography under native condition (Fig. 3e) and subsequently, kinase assay was performed for different time intervals. It was found that the SAPK9 phosphorylates OsbZIP23, and the phosphorylation efficiency is increased with the increase in reaction time (Fig. 3f). These results suggest that the SAPK9 protein has both autophosphorylation and transphosphorylation activities in vitro, and has the potentiality of transactivation of the OsbZIP transcription factors in vivo.
Generation of transgenic rice lines for SAPK9 overexpression and endogenous gene silencing
To investigate the gain-of-function and loss-of-function phenotypes of SAPK9 gene in rice through transgenesis, gene overexpression (OE) and RNAi-mediated gene silencing (RNAi) constructs were prepared using OsUbi1 promoter and nopaline synthase (NOS) transcription terminator in the pCAMBIA1301 binary plasmid (Additional file 2: Figure S5A and B). Both the constructs were then separately incorporated into drought-sensitive indica rice cultivar IR20 by Agrobacterium-mediated transformation. The putative OE and RNAi transformants were subjected to hygromycin selection during in vitro culture. The hygromycin selected transformants were verified by means of Southern hybridization using SAPK9 gene-specific probe to detect transgene integration patterns in T0 plants. Southern blot analysis detected single and multiple integration lines in both OE and RNAi (Fig. 4a and b). The OE lines of SAPK9 were designated as SAOE#1, 2, 3 etc., and RNAi lines were designated as RNAi#1, 2, 3 etc. The single integration T1 plants of OE and RNAi lines were examined for the expression level of SAPK9 transcript by real-time PCR. The real-time PCR analysis revealed that the relative expression level of SAPK9 transcript increased significantly (P < 0.01) in OE lines, but decreased significantly (P < 0.01) in RNAi lines compared to non-transgenic (NT) plants (Fig. 4c). Similarly, western blot analysis disclosed the increased expression of SAPK9 protein in OE lines in comparison to RNAi lines and NT plants (Fig. 4d), corroborating the real-time PCR result. Further, we have also checked the expression of SAPK9 gene at the reproductive stage of the transgenic and NT plants (Additional file 7: Figure S7), and the analysis revealed the similar level and trend of expression as observed in the vegetative stage. The high sequence similarity of SAPK9 with the other two members of SnRK2 subclass III in rice viz., SAPK8 and SAPK10, led us to analyse the relative expression levels of these two genes in transgenic lines and NT plants. It was observed that in OE and RNAi lines the relative expression levels of both the genes were not much altered in comparison to NT plants (Additional file 8: Figure S6).
Overexpression of SAPK9 improves drought tolerance in rice by increasing the water retention capacity in transgenic plants through osmotic adjustment and stomatal closure
To examine the drought stress tolerance capacity of SAPK9 OE and RNAi lines at the vegetative and early reproductive (panicle initiation) stages, hygromycin selected T1 plants were transferred into the pots and allowed to grow under glasshouse conditions for 30 days along with non-transgenic (NT) plants. After one week of transferring the pots in net-house, each set of plants from vegetative and panicle initiation stage were subjected to drought stress by withholding water supply until the visible effect of dehydration was noticed in RNAi lines and NT plants. After the duration of drought stress, plants were watered for 3 days in each case, and the survival rate (%) was determined. A significantly (P < 0.01) greater survival rate was observed in OE lines when compared with RNAi lines and NT plants during both vegetative stage (Fig. 5a) and reproductive stage (Fig. 5b). The observed increase and decrease in the drought tolerance ability of OE and RNAi lines, respectively, clearly indicated that the SAPK9 plays a strong role in drought tolerance in rice. Previous studies on either overexpression or silencing of the ABA-responsive kinase genes in rice and other plant species have documented similar results of drought tolerance or drought sensitivity, respectively [5, 9, 16, 31, 48].
Globally, water availability dictates the yield of crop plants including rice . Plant with high water holding capacity can sustain the severity of drought stress better than the others. The water status in a plant body is principally determined by two critical parameters- water loss rate (WLR) and relative water content (RWC) [31, 49]. The WLR and RWC were measured in leaf tissues of OE, RNAi lines, and NT plants. It was observed that the OE lines exhibited lower WLR and higher RWC in contrast to the RNAi lines, which showed higher WLR and lower RWC in comparison with NT plants (Fig. 5c). These results indicate that the SAPK9 plays an important functional role in enhancing the water retention capacity of the plants under water deficit condition.
Additionally, we investigated the mechanism by which the SAPK9 endows the OE rice plant for its improved water holding capacity. The two major physiological mechanisms that reduce water loss in plants under drought stress condition are the osmotic adjustment and stomatal closure. It has been established that plants under drought stress condition accumulate compatible osmolytes such as free proline and soluble sugars in order to keep the cellular structure intact by adjusting the intracellular osmotic potential . The physiological mechanism by which the overexpression of SAPK9 improves the water retention capacity of transgenic rice lines under drought stress was elucidated by determining the contents of proline and soluble sugar in the transgenic and NT plants. Analysis of the experimental data revealed not much difference in proline and soluble sugar contents of the OE, RNAi lines and NT plants while growing them under normal condition. However, upon drought stress imposition, higher contents of proline and soluble sugar were accumulated in the OE lines compared with RNAi and NT plants (Fig. 5d). The significant increase in the cellular osmolyte contents in OE lines during drought stress strongly suggests that the accumulation of these osmolytes occurs due to the enhanced activity of SAPK9, the possible molecular mechanism of which has been discussed in the later section. On the contrary, the reduced activity of SAPK9 in RNAi lines results in decreased accumulation of osmolytes, consequently leading to drought sensitivity. Similar phenotypic observations associated with cellular osmolyte concentration have been documented previously upon up-regulation or down-regulation of candidate genes [25, 31, 50, 51]. It is to be noted that in our study the analysis has been performed after a drought stress treatment over a period of 8 days under net-house condition (when the clear phenotypic difference was observed among OE, RNAi and NT plants). However, for a comprehensive physiological evaluation of OE and RNAi lines to compare with the NT plants, the analyses need to be carried out over a longer duration in the field condition, as has been suggested before . Previous reports have shown that the members of subclass III SnRK2, i.e., Arabidopsis OST1/SnRK2.6 and fava bean AAPK are involved in the ABA-regulated stomatal movement [5–7]. Hence in the present study, the stomatal status was observed and counted in OE lines, RNAi lines, and NT plants. Under normal growth conditions, the differences in the stomatal status of transgenic lines and NT plants were minimal. After the drought stress treatment, ~20–23 and ~27 % of stomata were completely closed in RNAi and NT plants, respectively; while greater proportions (~40–50 %) of closed stomata were present in OE lines (Fig. 6a). Accordingly, the proportions of completely opened stomata were significantly lowest in OE lines among the three sets of plants (Fig. 6a). Interestingly, the proportions of partially opened stomata in OE lines were comparatively highest. The variation in the stomatal movement of transgenic lines and NT plants led us to verify the expression level of two slow anion channel-associated (SLAC) genes viz., OsSLAC1 and OsSLAC7. The results showed that the transcript expression levels of both the genes were significantly (P < 0.01) higher in OE lines, compared with the RNAi lines and NT plants (Fig. 6b). Recent investigations have documented that these two stomatal proteins are involved in anion transport [53, 54], and the SAPK8 has been found to activate the function of OsSLAC1 . Our results suggest that the stomatal closure is positively regulated by the functional activity of SAPK9 upon OsSLAC1 and OsSLAC7 gene expression. It is anticipated that the phosphorylation ability of SAPK9 might activate the transcription factor(s) specifically required for the expression of these two genes. However, the detailed molecular mechanism of SAPK9-regulated stomatal closure requires further study.
Transgenic rice lines overexpressing SAPK9 show reduced lipid peroxidation and increased antioxidant activity
Stresses usually cause harm in plants via oxidative damage by generating reactive oxygen species (ROS), such as H2O2 and O2- ions. The ROS are toxic and highly reactive oxygen derivatives that damage macromolecules such as DNA, protein and carbohydrate, resulting in cell death . ROS causes lipid peroxidation, which results in the production of malondialdehyde (MDA) . The MDA is a stress-specific molecular marker that measures the extent of cellular damage due to the imposition stress in plants. Since the overexpression of SAPK9 improved drought tolerance in rice plants through osmotic adjustment and stomatal closure, hence we examined whether the SAPK9 functions in stress tolerance through ROS detoxification. The content of MDA in OE lines, RNAi lines, and NT plants was measured and analysed. There were marginal variations in the content of MDA in OE, RNAi lines, and NT plants during normal growth condition. However, under drought stress, the MDA content in OE lines was significantly (P < 0.01) lowest among the three sets of plants, although the lipid peroxidation levels increased in all plants compared to the normal condition (Fig. 7a). The reduced MDA level under drought stress in the OE lines implied that the cellular oxidative damage is less severe in these lines in comparison with RNAi lines and NT plants. Therefore, we were interested in documenting the ROS accumulation by staining the leaves of OE, RNAi lines, and NT plants with 3,3′-diaminobenzidine (DAB) and nitro blue tetrazolium (NBT). Under normal growth conditions, there were no differences in the localized peroxide and free radicals in leaf tissues amongst transgenic lines and NT plants. However, upon exposure to drought condition, the leaf tissues of RNAi and NT plant displayed higher levels of H2O2 and O2- accumulation in comparison with the OE lines (Fig. 7b and c). These results confirm that the SAPK9 overexpression in planta has enhanced the drought tolerance phenotype of OE plants by ROS detoxification, which resulted in less membrane damage.
Overexpression of SAPK9 causes elevated transcript expression of ABA- and stress-responsive genes resulting in increased ABA sensitivity at germination and post-germination stages of transgenic rice
To understand the comprehensive function of SAPK9 under drought stress condition, we analysed the expression profiles of some drought-inducible genes, such as TRAB1, OsbZIP23, OsbZIP46, OsLEA3-1, OsRab16B, and OsRab21, which work hierarchically downstream of SAPK9 in the ABA-signaling pathway. The TRAB1, OsbZIP23, and OsbZIP46 are the rice bZIP transcription factors, known to be induced by ABA and drought stress [25, 47, 56, 57]. Two dehydrin proteins are encoded by OsRab16B and OsRab21, which are transcriptionally responsive to abiotic stress [58, 59]. OsLEA3-1 is a late embryogenesis abundant (LEA) protein proposed to protect cell membrane structure stable by minimizing the damage caused by increasing ion concentration in dehydrating cells  and improves drought tolerance . We found that the relative expression level of TRAB1, OsbZIP23, OsbZIP46, OsLEA3-1, OsRab16B, and OsRab21 increased significantly (P < 0.01) in OE lines, compared with the NT plants (Fig. 8a). However, the levels of expression of these genes in RNAi lines were significantly (P < 0.01) lower in comparison with NT plants (Fig. 8a). These results suggest that the SAPK9 might be a transcriptional activator of these drought-stress associated genes in the ABA signaling pathway in rice, and overexpression or downregulation of SAPK9 alters the expression of these set of genes. Since most of the transcriptional activator proteins are usually localized in the nucleus, therefore we examined the subcellular localization of SAPK9 protein in onion epidermal cells. We found increased fluorescence in cytoplasm and nucleus due the expression of SAPK9–GFP fusion protein compared to that of only GFP, indicating the SAPK9 protein was localized in these cellular compartments (Fig. 8b). The subcellular localization of SAPK9 is similar to the previous reports on other SnRK2 family members [16–18]. The exact molecular mechanism of SAPK9-mediated transcriptional activation is yet to be investigated.
The role of phytohormone ABA in the maintenance of seed dormancy, hindrance of seed germination and inhibition of seedling growth is well documented . Since the expression of SAPK9 gene is profoundly induced by ABA (Fig. 2d), and ABA-responsive gene expressions are increased in OE lines (Fig. 8a); hence we have investigated the ABA sensitivity of transgenic lines at germination and post-germination stages. The seeds from two independent OE lines (SAOE#1 and SAOE#2), two RNAi lines (RNAi#4 and RNAi#5) and non-transgenic (NT) plants were germinated on MS agar media with variable concentrations of ABA (0, 1, 3 and 6 μM). There were no significant differences in the germination rate of OE lines at 0 and 1 μM ABA, but at 3 and 6 μM ABA concentration germination rate of OE lines was significantly lowest in comparison with RNAi lines and NT plants (Fig. 9a and b), implying that the SAPK9 OE lines are most sensitive to ABA at the germination stage. We also investigated the ABA sensitivity of OE and RNAi lines at the post-germination stage. The results showed no significant differences in the shoots and roots length of OE, RNAi and NT seedlings grown at 0 μM concentration of ABA for 14 days. However, at 1, 3 and 6 μM ABA concentrations, the length of shoots and roots of OE lines decreased significantly (P < 0.01) in comparison with RNAi lines and NT plants (Fig. 9c and d). These findings indicate that the overexpression of SAPK9 increases ABA sensitivity at the post-germination stage as well. Together with the finding on decreased ABA sensitivity of SAPK9 down-regulated RNAi lines, we infer that the SAPK9 is a positive regulator of ABA-dependent stress signaling pathway in rice. Our observations are consistent with the previous reports that overexpression and silencing of ABA-responsive gene results in differential ABA sensitivity at both germination and post-germination stages in rice [25, 31, 63].
Overexpression of SAPK9 increases grain yield in transgenic rice under drought stress by improving pollen maturation, spikelet fertility, and panicle weight
Since there was a striking difference between the transgenic OE and RNAi plants with respect to the duration of drought tolerance, therefore we were encouraged to examine the grain yield attributing traits of these plants grown under drought stress. For this, two sets of OE lines, RNAi lines and NT plants were grown in PVC pipes under glasshouse condition till they reach the panicle imitation stage (~60 days of seeding), and subsequently transferred to net-house. Upon attaining the flowering stage, plants were subjected to drought stress condition for a period of 10 days. After completion of the stress period plants were watered normally, and allowed to recover through the flowering stage until the seed maturation occurs (Fig. 10a). Under drought stress condition, the panicle weight and spikelet fertility of OE lines were found to be significantly (P < 0.01) highest among the three sets of plants (Fig. 10b, c and d). There was no obvious difference in grain length and width amongst transgenic lines and NT plants (Fig. 10b). These results indicate that the SAPK9 gene plays a significant role in increasing the grain yield in rice by influencing panicle weight and spikelet fertility. Since panicle weight and spikelet fertility are correlated with the pollen maturation and viability, therefore the number of matured and viable pollens in OE, RNAi, and NT plants were investigated using I2-KI staining. The analysis revealed that OE lines had a higher proportion of mature pollen than those of RNAi lines and NT plants (Fig. 10e). The proportion of mature pollen reflects the pollen viability. The more number of mature pollens correlating with the overexpression of SAPK9 suggested that the SAPK9 functions in increasing pollen viability. Additionally, we have already found that the relative expression level of SAPK9 is significantly high in panicle under drought stress condition (Fig. 2c). Collectively, these results document that the SAPK9 enhances spikelet fertility by regulating pollen maturation and thereby improves grain yield. Interestingly, previous reports have documented that the SAPK9 plays some functional role in the development of panicle  as well in the pollen  via ABA-signaling pathway. However, the exact molecular mechanism by which the SAPK9 gene regulates the spikelet fertility and grain yield in rice remain to be elucidated.
Possible molecular mechanism of SAPK9 action
Based on the findings of the present study and the information available in the literature, a molecular interactive pathway for SAPK9 has been proposed (Additional file 9: Figure S8). Our work clearly established that the SAPK9 acts as a transactivating kinase and a potential transcriptional activator for modulating the activity of several key components in the ABA-dependent stress adaptation pathway in rice plant. We hypothesize that an unknown factor(s) might be involved in mediating the transcriptional activity of SAPK9, which requires further investigation.
We have characterized the structural, biochemical and physiological functions of SAPK9, a subclass III SnRK2 gene of rice through developing gain-of-function and loss-of-function mutants by transgenesis. The SAPK9 transcript expression has been found to be differentially regulated in the selected drought-tolerant and drought-sensitive rice genotypes, and its expression is comparatively more elevated in reproductive stage than the vegetative stage. The SAPK9 transcript has been observed to be upregulated by drought stress and ABA treatments. The SAPK9 protein has characteristic kinase fold structure, and is able to transactivate its target substrates, including histone III and OsbZIP23; and itself is autophosphorylated. Constitutive overexpression of SAPK9 gene from a drought-tolerant wild rice genotype and RNAi-mediated silencing of endogenous SAPK9 gene in a drought-sensitive cultivar have revealed that the SAPK9 positively regulates drought stress tolerance by boosting the osmotic adjustment and stomatal closure of the plant. The enhanced drought-tolerant transgenic rice lines overexpressing SAPK9 display less cellular oxidative damage resulting from the reduced reactive oxygen species accumulation. The SAPK9 overexpressed lines exhibit increased sensitivity to exogenous ABA and enhanced transcription of other hierarchically downstream ABA-responsive genes, indicating that the SAPK9, indeed, is a positive regulator of ABA-mediated stress signaling pathway in rice plant. The increased functional activity of SAPK9 overexpression in transgenic lines improves grain yield-related traits such as panicle weight and spikelet fertility by increasing the pollen viability. Together, the present findings strengthen our knowledge about the functional role of SAPK9 as transactivating kinase and a probable transcriptional activator, which can be utilized as a promising gene-based molecular marker in transgenic breeding for generating crop plants with improved drought tolerance and grain yield.
ABA, Abscisic acid; CDS, Coding DNA sequence; DAB, 3,3′-diaminobenzidine; MDA, Malondialdehyde; NBT, Nitro blue tetrazolium; NT, Non-transgenic; OE, Overexpression; RNAi, gene silencing by RNA interference; ROS, Reactive oxygen species; RWC, Relative water content; SAPK, Osmotic stress/ABA-activated protein kinase; SLAC, Slow anion channel-associated; SNP, Single nucleotide polymorphism; SnRK2, Sucrose non-fermenting 1-related kinase 2; WLR, Water loss rate.
Cutler SR, Rodriguez PL, Finkelstein RR, Abrams SR. Abscisic acid: emergence of a core signaling network. Annu Rev Plant Biol. 2010;61:651–79.
Hong SW, Jon JH, Kwak JM, Nam HG. Identification of a receptor-like protein kinase gene rapidly induced by abscisic acid, dehydration, high salt, and cold treatments in Arabidopsis thaliana. Plant Physiol. 1997;113(4):1203–12.
Fujii H, Verslues PE, Zhu JK. Arabidopsis decuple mutant reveals the importance of SnRK2 kinases in osmotic stress responses in vivo. Proc Natl Acad Sci. 2011;108(4):1717–22.
Kulik A, Wawer I, Krzywińska E, Bucholc M, Dobrowolska G. SnRK2 protein kinases—key regulators of plant response to abiotic stresses. OMICS. 2011;15(12):859–72.
Li J, Wang XQ, Watson MB, Assmann SM. Regulation of abscisic acid-induced stomatal closure and anion channels by guard cell AAPK kinase. Science. 2000;287(5451):300–3.
Mustilli AC, Merlot S, Vavasseur A, Fenzi F, Giraudat J. Arabidopsis OST1 protein kinase mediates the regulation of stomatal aperture by abscisic acid and acts upstream of reactive oxygen species production. Plant Cell. 2002;14(12):3089–99.
Yoshida R, Hobo T, Ichimura K, Mizoguchi T, Takahashi F, Aronso J, Shinozaki K. ABA-activated SnRK2 protein kinase is required for dehydration stress signaling in Arabidopsis. Plant Cell Physiol. 2002;43(12):1473–83.
Umezawa T, Yoshida R, Maruyama K, Yamaguchi-Shinozaki K, Shinozaki K. SRK2C, a SNF1-related protein kinase 2, improves drought tolerance by controlling stress-responsive gene expression in Arabidopsis thaliana. Proc Natl Acad Sci U S A. 2004;101(49):17306–11.
Boudsocq M, Barbier-Brygoo H, Laurière C. Identification of nine sucrose nonfermenting 1-related protein kinases 2 activated by hyperosmotic and saline stresses in Arabidopsis thaliana. J Biol Chem. 2004;279(40):41758–66.
Boudsocq M, Laurière C. Osmotic signaling in plants. Multiple pathways mediated by emerging kinase families. Plant Physiol. 2005;138(3):1185–94.
Fujita Y, Yoshida T, Yamaguchi‐Shinozaki K. Pivotal role of the AREB/ABF‐SnRK2 pathway in ABRE‐mediated transcription in response to osmotic stress in plants. Physiologiaplantarum. 2013;147(1):15–27.
Huai J, Wang M, He J, Zheng J, Dong Z, Lv H, et al. Cloning and characterization of the SnRK2 gene family from Zea mays. Plant Cell Rep. 2008;27(12):1861–8.
Ying S, Zhang DF, Li HY, Liu YH, Shi YS, Song YC, et al. Cloning and characterization of a maize SnRK2 protein kinase gene confers enhanced salt tolerance in transgenic Arabidopsis. Plant Cell Rep. 2011;30(9):1683–99.
Anderberg RJ, Walker-Simmons MK. Isolation of a wheat cDNA clone for an abscisic acid-inducible transcript with homology to protein kinases. Proc Natl Acad Sci. 1992;89(21):10183–7.
Xu ZS, Liu L, Ni ZY, Liu P, Chen M, Li LC, et al. W55a encodes a novel protein kinase that is involved in multiple stress responses. J Integr Plant Biol. 2009;51(1):58–66.
Tian S, Mao X, Zhang H, Chen S, Zhai C, Yang S, et al. Cloning and characterization of TaSnRK2. 3, a novel SnRK2 gene in common wheat. J Exp Bot. 2013;64(7):2063–80.
Mao X, Zhang H, Tian S, Chang X, Jing R. TaSnRK2. 4, an SNF1-type serine/threonine protein kinase of wheat (Triticum aestivum L.), confers enhanced multistress tolerance in Arabidopsis. J Exp Bot. 2010;61(3):683–96.
Zhang H, Mao X, Wang C, Jing R. Overexpression of a Common Wheat Gene TaSnRK2.8 Enhances Tolerance to Drought, Salt and Low Temperature in Arabidopsis. PLoS ONE. 2010;5(12):e16041.
Kagaya Y, Hobo T, Murata M, Ban A, Hattori T. Abscisic acid–induced transcription is mediated by phosphorylation of an abscisic acid response element binding factor, TRAB1. Plant Cell. 2002;14(12):3177–89.
Johnson RR, Wagner RL, Verhey SD, Walker-Simmons MK. The abscisic acid-responsive kinase PKABA1 interacts with a seed-specific abscisic acid response element-binding factor, TaABF, and phosphorylates TaABF peptide sequences. Plant Physiol. 2002;130(2):837–46.
Furihata T, Maruyama K, Fujita Y, Umezawa T, Yoshida R, Shinozaki K, et al. Abscisic acid-dependent multisite phosphorylation regulates the activity of a transcription activator AREB1. Proc Natl Acad Sci U S A. 2006;103(6):1988–93.
Li J, Kinoshita T, Pandey S, Ng CKY, Gygi SP, Shimazaki KI, et al. Modulation of an RNA-binding protein by abscisic-acid-activated protein kinase. Nature. 2002;418(6899):793–7.
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(5):1163–77.
Diédhiou CJ, Popova OV, Dietz KJ, Golldack D. The SNF1-type serine-threonine protein kinase SAPK4 regulates stress-responsive gene expression in rice. BMC Plant Biol. 2008;8(1):49.
Dey A, Samanta MK, Gayen S, Sen SK, Maiti MK. Enhanced Gene Expression Rather than Natural Polymorphism in Coding Sequence of the OsbZIP23 Determines Drought Tolerance and Yield Improvement in Rice Genotypes. PLoS ONE. 2016;11(3):e0150763. doi:10.1371/journal.pone.0150763.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2− ΔΔCT method. Methods. 2001;25(4):402–8.
Sambrook J, Fritsch EF, Maniatis T. Molecular Cloning: A Laboratory Manual. Cold Spring Harbor: Cold Spring Harbor Laboratory; 1989.
Feng CZ, Chen Y, Wang C, Kong YH, Wu WH, Chen YF. Arabidopsis RAV1 transcription factor, phosphorylated by SnRK2 kinases, regulates the expression of ABI3, ABI4, and ABI5 during seed germination and early seedling development. Plant J. 2014;80(4):654–68.
Kumar KK, Maruthasalam S, Loganathan M, Sudhakar D, Balasubryamaniam P. An improved Agrobacterium-mediated transformation protocol for recalcitrant elite indica rice cultivars. Plant MolBiol Rep. 2005;23:67–73.
Hu W, Huang C, Deng XM, Zhou SY, Chen LH, Li Y, et al. TaASR1, a transcription factor gene in wheat, confers drought stress tolerance in transgenic tobacco. Plant Cell Environ. 2013;36:1449–64.
Wei S, Hu W, Deng X, Zhang Y, Liu X, Zhao X, et al. A rice calcium-dependent protein kinase OsCPK9 positively regulates drought stress tolerance and spikelet fertility. BMC Plant Biol. 2014;14(1):133.
Ye H, Du H, Tang N, Li X, Xiong L. Identification and expression profiling analysis of TIFY family genes involved in stress and phytohormone responses in rice. Plant Mol Biol. 2009;71(3):291–305.
Zhang Z, Huang R. Analysis of Malondialdehyde, Chlorophyll, Proline, Soluble sugar, Glutathione content in seedling. 2013;3(14):e817. Avaliable: http://www.bio-protocol.org/e817. Jan 2015
Zhang L, Xiao S, Li W, Feng W, Li J, Wu Z, et al. Overexpression of a Harpin-encoding gene hrf1 in rice enhances drought tolerance. Journal of experimental botany. 2011;62(12):4229–38.
Baryla A, Laborde C, Montillet JL, Triantaphylides C, Chagvardieff P. Evaluation of lipid peroxidation as a toxicity bioassay for plants exposed to copper. Environ Pollut. 2000;109:131–5.
Tiwari V, Chaturvedi AK, Mishra A, Jha B. Introgression of the SbASR-1 Gene Cloned from a Halophyte Salicornia brachiata Enhances Salinity and Drought Endurance in Transgenic Groundnut (Arachis hypogaea) and Acts as a Transcription Factor. PLoS One. 2015;10(7):e0131567.
Xu K, Huang X, Wu M, Wang Y, Chang Y, Liu K, et al. A Rapid, Highly Efficient and Economical Method of Agrobacterium-Mediated In planta Transient Transformation in Living Onion Epidermis. PLoS One. 2014;9(1):e83556. doi:10.1371/journal.pone.0083556.
Zou M, Guan Y, Ren H, Zhang F, Chen F. A bZIP transcription factor, OsABI5, is involved in rice fertility and stress tolerance. Plant Mol Biol. 2008;66(6):675–83.
Kumar GR, Sakthivel K, Sundaram RM, Neeraja CN, Balachandran SM, Rani NS, et al. Allele mining in crops: prospects and potentials. Biotechnol Adv. 2010;28(4):451–61.
Parida SK, Mukerji M, Singh AK, Singh NK, Mohapatra T. SNPs in stress-responsive rice genes: validation, genotyping, functional relevance and population structure. BMC Genomics. 2012;13(1):426.
Fujii H, Verslues PE, Zhu JK. Identification of two protein kinases required for abscisic acid regulation of seed germination, root growth, and gene expression in Arabidopsis. Plant Cell. 2007;19(2):485–94.
Ng LM, Soon FF, Zhou XE, West GM, Kovach A, Suino-Powell KM, et al. Structural basis for basal activity and autoactivation of abscisic acid (ABA) signaling SnRK2 kinases. Proc Natl Acad Sci. 2011;108(52):21259–64.
Belin C, De Franco PO, Bourbousse C, Chaignepain S, Schmitter JM, Vavasseur A, et al. Identification of features regulating OST1 kinase activity and OST1 function in guard cells. Plant Physiol. 2006;141(4):1316–27.
Vilela B, Moreno A, Capellades M, Pagès M, Lumbreras V. ZmSnRK2.8 responds to ABA through the SnRK2-PP2C complex. Maydica. 2012;57(1):11–8.
Uno Y, Furihata T, Abe H, Yoshida R, Shinozaki K, Yamaguchi-Shinozaki K. Arabidopsis basic leucine zipper transcription factors involved in an abscisic acid-dependent signal transduction pathway under drought and high-salinity conditions. Proc Natl Acad Sci USA. 2000;97:11632–7.
Kobayashi Y, Murata M, Minami H, Yamamoto S, Kagaya Y, Hobo T, et al. Abscisic acid‐activated SNRK2 protein kinases function in the gene‐regulation pathway of ABA signal transduction by phosphorylating ABA response element‐binding factors. Plant J. 2005;44(6):939–49.
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(4):1755–68.
Ouyang SQ, Liu YF, Liu P, Lei G, He SJ, Ma B, et al. Receptor‐like kinase OsSIK1 improves drought and salt stress tolerance in rice (Oryza sativa) plants. Plant J. 2010;62(2):316–29.
Gao Y, Jiang W, Dai Y, Xiao N, Zhang C, Li H, et al. A maize phytochrome-interacting factor 3 improves drought and salt stress tolerance in rice. Plant Mol Biol. 2015;87(4–5):413–28.
Cui M, Zhang W, Zhang Q, Xu Z, Zhu Z, Duan F, Wu R. Induced over-expression of the transcription factor OsDREB2A improves drought tolerance in rice. Plant Physiol Biochem. 2011;49(12):1384–91.
Xiang Y, Huang Y, Xiong L. Characterization of stress-responsive CIPK genes in rice for stress tolerance improvement. Plant Physiol. 2007;144(3):1416–28.
Lawlor DW. Genetic engineering to improve plant performance under drought: physiological evaluation of achievements, limitations, and possibilities. J Exp Bot. 2013;64(1):83–108.
Sun SJ, Qi GN, Gao QF, Wang HQ, Yao FY, Hussain J, et al. 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.
Fan X, Wu J, Chen T, Tie W, Chen H, Zhou F, et al. 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.
Mittler R. Oxidative stress, antioxidants and stress tolerance. Trends Plant Sci. 2002;7(9):405–10.
Hobo T, Kowyama Y, Hattori T. A bZIP factor, TRAB1, interacts with VP1 and mediates abscisic acid-induced transcription. Proc Natl Acad Sci. 1999;96(26):15348–53.
Xiang Y, Tang N, Du H, Ye H, Xiong L. Characterization of OsbZIP23 as a key player of the basic leucine zipper transcription factor family for conferring abscisic acid sensitivity and salinity and drought tolerance in rice. Plant Physiol. 2008;148:1938–52.
Yamaguchi-Shinozaki K, Mundy J, Chua NH. Four tightly linked rab genes are differentially expressed in rice. Plant Mol Biol. 1990;14(1):29–39.
Mundy J, Chua NH. Abscisic acid and water-stress induce the expression of a novel rice gene. EMBO J. 1988;7(8):2279.
Liu C, Wu Y, Wang X. bZIP transcription factor OsbZIP52/RISBZ5: a potential negative regulator of cold and drought stress response in rice. Planta. 2012;235(6):1157–69.
Xiao B, Huang Y, Tang N, Xiong L. Over-expression of a LEA gene in rice improves drought resistance under the field conditions. Theor Appl Genet. 2007;115(1):35–46.
Finkelstein RR, Gampala SS, Rock CD. Abscisic acid signaling in seeds and seedlings. Plant Cell. 2002;14(Suppl):S15–45.
Yang A, Dai X, Zhang, WH. A R2R3-type MYB gene, OsMYB2, is involved in salt, cold, and dehydration tolerance in rice. Journal of experimental botany. 2012;63(7):2541–56.
Zhang X, Li J, Liu A, Zou J, Zhou X, Xiang J, et al. Expression Profile in Rice Panicle: Insights into Heat Response Mechanism at Reproductive Stage. PLoS One. 2012;7(11):e49652.
Hirano K, Aya K, Hobo T, Sakakibara H, Kojima M, Shim RA, et al. Comprehensive transcriptome analysis of phytohormone biosynthesis and signaling genes in microspore/pollen and tapetum of rice. Plant Cell Physiol. 2008;49(10):1429–50.
This work is dedicated to the memory of Late Professor Soumitra K. Sen who was involved with this work from the time it was conceived and designed. We express our gratitude to Prof. Saumen Hajra, Centre of Biomedical Research, Lucknow for his valuable guidance and encouragement throughout this study. We are thankful to Ms. Anuttama Dutta for providing help and support during a few experimental works. The technical help from Mr. Meghnath Prasad, Mr. Partha Das, Mr. Sona Dogra, Mr. Uttam Dogra, Mrs. Gayatri Aditya and Mr. Manoj Aditya are also acknowledged. Department of Science and Technology, Government of India is sincerely acknowledged for providing INSPIRE Fellowship to AD during his research tenure.
This work was supported by Indian Council of Agricultural Research, Government of India (NAIP-PAR-C4/C/30033) to Late Professor Soumitra K. Sen.
Availability of data and materials
All relevant data are within this article and its supporting information files. Gene sequences have been deposited in the GenBank database.
AD designed and performed the major experiments, analysed the results and wrote the manuscript. MKS performed few experiments and interpreted the data. SG designed few experiments and analysed the data. MKM conceived and designed the study, interpreted the results and corrected the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
List of primers used in this study. (DOCX 21 kb)
Schematic representation of the genetic constructs based on pCAMBIA1301 plasmid used for Agrobacterium-mediated transformation of drought-sensitive indica rice cultivar IR20. (A) The gene overexpression (OE) construct of SAPK9 carrying 1086 bp CDS from O. rufipogon. The developed transgenic rice lines were designated as SAOE#1, 2, 3 etc. (B) The RNAi-mediated gene silencing (RNAi) construct of endogenous SAPK9 gene. The 605 bp 5′-part of SAPK9 CDS from O. rufipogon was cloned in sense and antisense orientation flanking an arbitrary 200 bp DNA linker. The developed transgenic rice lines were designated as RNAi#1, 2, 3 etc. (TIF 453 kb)
PCR-amplified SAPK9 CDS from the chosen rice genotypes. Lanes 1-11 represents O. rufipogon, O. nivara, Nagina22, Manipuri, Vandana, Swarna, IR20, IR36, IR64, IR72, HRC300, respectively. Lane M- standard molecular weight marker. (TIF 1083 kb)
Multiple sequence alignment for analysis of amino acid polymorphism in the SAPK9 CDS from the selected rice genotypes using Jalview software. (TIF 3333 kb)
Amino acid sequence alignment of the isolated SAPK9 CDS from wild rice Oryza rufipogon (accession no. KT387673) with the reported sequence of japonica rice cultivar (accession no.AB125310). (TIF 487 kb)
Multiple alignments showing amino acid sequence similarity of SAPK9 with Arabidopsis OST1/SnRK2.6 (PDB ID: 3UC4) and SnRK2.3 (PDB ID: 3UC3). (TIF 644 kb)
The relative expression level of SAPK9 gene in reproductive stage (panicle initiation) was analysed through real-time PCR from leaf tissues of three sets of plants, i.e. OE, RNAi and NT plants. For internal reference, rice polyubiquitin1 (OsUbi1) gene was used. Error bars represent the mean ± SD of triplicate measurements. Student’s t-test was performed to find out statistically significant differences (**P < 0.01). (TIF 240 kb)
Relative expression level of SAPK8 and SAPK10 gene in SAPK9 overexpressed (OE) and gene silenced (RNAi) transgenic rice lines. Analysis of real-time PCR depicting transcript level of (A) SAPK8 and (B) SAPK10 in leaf tissues of three sets of plants, i.e. OE, RNAi and non-transgenic (NT) plants under drought stress. For internal reference, rice polyubiquitin1 (OsUbi1) gene was used. Error bars represent the mean ± SD of triplicate measurements. Student’s t-test was performed to find out statistically significant differences (*P < 0.01). (TIF 339 kb)
Possible molecular mechanism of SAPK9-regulated drought stress tolerance in rice. Upon drought stress, the SAPK9 is activated by ABA signaling pathway, which in turn phosphorylates the bZIP transcription factors activating the expression of downstream genes involved in maintaining the physiological conditions for cellular homeostasis. Similarly, the SAPK9 also transcriptionally activates (might require an unknown factor) for upregulation of several hierarchically downstream genes including the anion channel genes leading to stomatal closure. (TIF 899 kb)
About this article
Cite this article
Dey, A., Samanta, M.K., Gayen, S. et al. 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 16, 158 (2016). https://doi.org/10.1186/s12870-016-0845-x