- Research article
- Open Access
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
© The Author(s). 2016
- Received: 16 March 2016
- Accepted: 5 July 2016
- Published: 13 July 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.
- Abscisic acid (ABA)
- Drought tolerance
- Gene silencing
- Grain yield
- Osmotic potential
- Rice crop
- Stomatal closure
- Stress-responsive gene
- Sucrose non-fermenting 1-related kinase 2 (SnRK2)
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).
The single copy SAPK9 gene displays differential transcript expression profile and allelic polymorphism in drought-tolerant and drought-sensitive 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 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
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
Overexpression of SAPK9 improves drought tolerance in rice by increasing the water retention capacity in transgenic plants through osmotic adjustment and stomatal closure
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.
Transgenic rice lines overexpressing SAPK9 show reduced lipid peroxidation and increased antioxidant activity
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
Overexpression of SAPK9 increases grain yield in transgenic rice under drought stress by improving pollen maturation, spikelet fertility, and panicle weight
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.
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
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.
- Cutler SR, Rodriguez PL, Finkelstein RR, Abrams SR. Abscisic acid: emergence of a core signaling network. Annu Rev Plant Biol. 2010;61:651–79.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Boudsocq M, Laurière C. Osmotic signaling in plants. Multiple pathways mediated by emerging kinase families. Plant Physiol. 2005;138(3):1185–94.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.Google Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Sambrook J, Fritsch EF, Maniatis T. Molecular Cloning: A Laboratory Manual. Cold Spring Harbor: Cold Spring Harbor Laboratory; 1989.Google Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.Google Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Mittler R. Oxidative stress, antioxidants and stress tolerance. Trends Plant Sci. 2002;7(9):405–10.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Mundy J, Chua NH. Abscisic acid and water-stress induce the expression of a novel rice gene. EMBO J. 1988;7(8):2279.PubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Finkelstein RR, Gampala SS, Rock CD. Abscisic acid signaling in seeds and seedlings. Plant Cell. 2002;14(Suppl):S15–45.PubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar