Skip to main content

Main menu

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
    • Front Matter Portal
    • Journal Club
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses
  • Submit
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Accessibility Statement
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian

User menu

  • Log in
  • My Cart

Search

  • Advanced search
Home
Home
  • Log in
  • My Cart

Advanced Search

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
    • Front Matter Portal
    • Journal Club
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses
  • Submit
Research Article

Plant Raf-like kinase integrates abscisic acid and hyperosmotic stress signaling upstream of SNF1-related protein kinase2

Masashi Saruhashi, Totan Kumar Ghosh, Kenta Arai, Yumiko Ishizaki, Kazuya Hagiwara, Kenji Komatsu, Yuh Shiwa, Keiichi Izumikawa, Harunori Yoshikawa, Taishi Umezawa, Yoichi Sakata, and Daisuke Takezawa
  1. aGraduate School of Science and Engineering, Saitama University, Saitama 338-8570, Japan;
  2. bDepartment of Bioproduction Technology, Junior College of Tokyo University of Agriculture, Tokyo 156-8502, Japan;
  3. cGenome Research Center, NODAI Research Institute, Tokyo University of Agriculture, Tokyo 156-8502, Japan;
  4. dFaculty of Agriculture, Tokyo University of Agriculture and Technology, Tokyo 183-8538, Japan;
  5. eCentre for Gene Regulation and Expression, School of Life Sciences, University of Dundee, Dundee DD1 5EH, Scotland, United Kingdom;
  6. fGraduate School of Bio-Applications and Systems Engineering, Tokyo University of Agriculture and Technology, Tokyo 184-8588, Japan;
  7. gPrecursory Research for Embryonic Science and Technology (PRESTO), Japan Science and Technology Agency, Tokyo 102-0007 Japan;
  8. hDepartment of BioScience, Tokyo University of Agriculture, Tokyo 156-8502, Japan;
  9. iInstitute for Environmental Science and Technology, Saitama University, Saitama 338-8570, Japan

See allHide authors and affiliations

PNAS November 17, 2015 112 (46) E6388-E6396; first published November 4, 2015; https://doi.org/10.1073/pnas.1511238112
Masashi Saruhashi
aGraduate School of Science and Engineering, Saitama University, Saitama 338-8570, Japan;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Totan Kumar Ghosh
aGraduate School of Science and Engineering, Saitama University, Saitama 338-8570, Japan;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Kenta Arai
aGraduate School of Science and Engineering, Saitama University, Saitama 338-8570, Japan;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Yumiko Ishizaki
aGraduate School of Science and Engineering, Saitama University, Saitama 338-8570, Japan;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Kazuya Hagiwara
aGraduate School of Science and Engineering, Saitama University, Saitama 338-8570, Japan;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Kenji Komatsu
bDepartment of Bioproduction Technology, Junior College of Tokyo University of Agriculture, Tokyo 156-8502, Japan;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Yuh Shiwa
cGenome Research Center, NODAI Research Institute, Tokyo University of Agriculture, Tokyo 156-8502, Japan;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Keiichi Izumikawa
dFaculty of Agriculture, Tokyo University of Agriculture and Technology, Tokyo 183-8538, Japan;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Harunori Yoshikawa
dFaculty of Agriculture, Tokyo University of Agriculture and Technology, Tokyo 183-8538, Japan;
eCentre for Gene Regulation and Expression, School of Life Sciences, University of Dundee, Dundee DD1 5EH, Scotland, United Kingdom;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Taishi Umezawa
fGraduate School of Bio-Applications and Systems Engineering, Tokyo University of Agriculture and Technology, Tokyo 184-8588, Japan;
gPrecursory Research for Embryonic Science and Technology (PRESTO), Japan Science and Technology Agency, Tokyo 102-0007 Japan;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Yoichi Sakata
hDepartment of BioScience, Tokyo University of Agriculture, Tokyo 156-8502, Japan;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Daisuke Takezawa
aGraduate School of Science and Engineering, Saitama University, Saitama 338-8570, Japan;
iInstitute for Environmental Science and Technology, Saitama University, Saitama 338-8570, Japan
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: takezawa@mail.saitama-u.ac.jp
  1. Edited by Julian I. Schroeder, University of California, San Diego, La Jolla, CA, and approved October 8, 2015 (received for review June 10, 2015)

  • Article
  • Figures & SI
  • Info & Metrics
  • PDF
Loading

Significance

Plants can sense loss of water caused by drought and stimulate internal mechanisms for protecting cells from damage with the aid of the stress hormone abscisic acid (ABA). Analysis of a mutant of the basal land plant, the moss Physcomitrella patens, revealed that an impairment of a protooncogene Raf-like protein kinase, designated “ARK” (for “ABA and abiotic stress-responsive Raf-like kinase”), causes a loss of both ABA sensitivity and osmotic stress tolerance. We show evidence that ARK has a role in integrating ABA and osmotic signals upstream of the sucrose nonfermenting 1-related protein kinase2, known to be a central regulator of stress signaling in plants.

Abstract

Plant response to drought and hyperosmosis is mediated by the phytohormone abscisic acid (ABA), a sesquiterpene compound widely distributed in various embryophyte groups. Exogenous ABA as well as hyperosmosis activates the sucrose nonfermenting 1 (SNF1)-related protein kinase2 (SnRK2), which plays a central role in cellular responses against drought and dehydration, although the details of the activation mechanism are not understood. Analysis of a mutant of the moss Physcomitrella patens with reduced ABA sensitivity and reduced hyperosmosis tolerance revealed that a protein kinase designated “ARK” (for “ABA and abiotic stress-responsive Raf-like kinase”) plays an essential role in the activation of SnRK2. ARK encoded by a single gene in P. patens belongs to the family of group B3 Raf-like MAP kinase kinase kinases (B3-MAPKKKs) mediating ethylene, disease resistance, and salt and sugar responses in angiosperms. Our findings indicate that ARK, as a novel regulatory component integrating ABA and hyperosmosis signals, represents the ancestral B3-MAPKKKs, which multiplied, diversified, and came to have specific functions in angiosperms.

  • abscisic acid
  • Physcomitrella patens
  • Raf-like kinase
  • SnRK2

Plant responses to water deficit stress are among the most important adaptations required for life on land. Ever since the colonization of land by the algal ancestor in the Middle Ordovician, embryophytes have developed physiological mechanisms for sensing and transducing environmental signals to protect themselves from a forthcoming severe water deficit. Abscisic acid (ABA), which is found in various embryophyte groups, is one of key molecules that play a role in water stress responses in the terrestrial environment. In angiosperms, ABA induces closure of stomata in leaves for prevention of water loss, and it maintains dormancy and desiccation tolerance in maturing seeds. In the moss Physcomitrella patens ABA induces tolerance to freezing, hyperosmosis, desiccation of protonema cells, and the accumulation of a number of stress-related transcripts such as those encoding late embryogenesis abundant (LEA)-like hydrophilic proteins, which are thought to be important for protecting cells from dehydration-induced damage to membranes and other cellular macromolecules (1⇓–3). These facts suggest conserved roles of ABA in dehydration stress tolerance, although whether common signaling mechanisms operate in various land plant groups has not been clarified.

Studies of the signal-transduction process in angiosperms have shed light on a critical role of sucrose nonfermenting 1-related protein kinase2 (SnRK2) in the ABA and dehydration responses. SnRK2 was first identified as an ABA-activated protein kinase in fava beans that mediates ABA-induced closure of leaf stomata (4) and also was identified as a factor encoded by OPEN STOMATA1 (OST1) in Arabidopsis thaliana, mutation of which causes an ABA-hyposensitive stomata closure response (5). Further analysis revealed that SnRK2 plays a pivotal role in the signaling of ABA by phosphorylating key transcription factors and ion channels (6, 7). SnRK2s also are activated by hyperosmotic stress, indicating a role in drought tolerance associated with the induction of stress-responsive genes (8, 9). SnRK2s of angiosperms are categorized into three subclasses: I, II, and III (8, 10). In-gel kinase assays using histone IIIS as a substrate revealed that subclass I SnRK2s are activated by osmotic stress but not by ABA, whereas the subclass II and subclass III SnRK2s are activated by both ABA and osmotic stress, with subclass III being activated more strongly than subclass II by ABA. Analysis of an Arabidopsis mutant null for three subclass III SnRK2 isoforms (srk2dsrk2esrk2i or snrk2.2snrk2.3snrk2.6) has shown that subclass III members play an essential role in the signal transduction of the ABA response (11, 12). Phosphopeptide mapping of subclass III SnRK2 revealed that ABA induced phosphorylation of Ser and Thr residues in the activation loop between kinase domains VII and VIII (12, 13). Furthermore, discovery of soluble ABA receptors (PYR/PYL/RCAR) has led to the establishment of the core model for the activation of SnRK2. In this model, activity of subclass III SnRK2s is suppressed by group A protein phosphatase 2C (PP2Cs) in the absence of ABA, and binding of ABA to the PYR/PYL/RCAR inhibits PP2Cs, thus allowing activation of SnRK2 (14⇓–16). Phosphorylation in the activation loop of SnRK2 is thought to be completed by autophosphorylation (11), although the details of the mechanism remain to be clarified.

Entire genome sequencing of P. patens (17) has revealed that four PYR/PYL/RCARs, two group A PP2Cs, and four SnRK2s are encoded in the moss genome, indicating that the core members needed for the ABA response are conserved from bryophytes to angiosperms (18). All SnRK2s encoded in the P. patens genome were assigned to the subclass III category, suggesting that this subclass represents the ancestral members of SnRK2. The P. patens PpOST1-1 gene, encoding SnRK2, was shown to complement the ABA-insensitive stomatal response of ost1 in transgenic Arabidopsis, and disruption of PpOST1-1 in P. patens resulted in an ABA-insensitive stomatal response in the sporophyte (19). That study suggested possible conservation of SnRK2 function in land plants. A recent study also indicated a role of group A PP2Cs as negative regulators of ABA signaling in bryophytes. A transgenic line lacking both group A PP2C genes (PpABI1A and PpABI1B) exhibited constitutive desiccation tolerance with the accumulation of a number of stress-associated transcripts regardless of ABA treatment (20). In this knockout line, however, SnRK2 was activated by exogenous ABA in a manner similar to that in the WT line (20), indicating group A PP2Cs have a minor role in the regulation of SnRK2 and the presence of a previously unidentified mechanism for its activation.

To understand better the molecular mechanisms underlying the ABA response in bryophytes, we mutagenized protonema cells of P. patens by UV irradiation and obtained mutant lines that could grow better than the WT line on medium containing ABA. One of the mutants, designated “AR7,” showed an ABA-insensitive growth response and reduced freezing-stress tolerance (21). AR7 also showed reduced responses to dehydration and cold (22), indicating that the molecule integrating ABA and stress-signaling processes is impaired in AR7. This report characterizes AR7 further and identifies the gene altered in AR7. Comparative genomic sequence analysis between AR7 and the WT line revealed that AR7 has a mutation in the gene encoding a group B3 Raf-like MAP kinase kinase kinase (B3-MAPKKK), subsequently designated “ARK” (for “ABA and abiotic stress responsive Raf-like kinase”). Further analyses indicated that ARK plays a key role in SnRK2-mediated ABA and hyperosmosis responses in the moss. On the basis of these collective findings, we propose that, in addition to PYR/PYL/RCAR and PP2C, ARK is an essential element of signaling components for the regulation of SnRK2 in basal land plants.

Results

Analysis of Transcriptome and Protein Kinase Activity in AR7.

The AR7 line has been characterized as a mutant with reduced ABA sensitivity in various aspects such as growth inhibition, morphological changes leading to brood cell formation, accumulation of soluble sugars and LEA-like boiling-soluble proteins, and tolerance of freezing and dehydration stress (21, 22). Microarray analysis of transcript profiles revealed that the expression of a large portion of the ABA-responsive genes was affected in AR7. Expression of 518 of 579 ABA up-regulated genes in the WT line was reduced significantly in AR7 (Fig. 1A, Left). These genes included 27 LEA-like genes encoding highly hydrophilic polypeptides, which are considered to be associated with dehydration stress tolerance (Table S1). RNA gel-blot analysis of four representative LEA-like transcripts, 17B9 (Pp1s294_52V6.1), 6A5 (Pp1s267_21V6.1), LEAIII (Pp1s118_232V6.1), and 19C6 (Pp1s52_261V6.1), confirmed impaired ABA-induced gene expression in AR7 (Fig. 1B). Furthermore, levels of 150 of 165 ABA down-regulated transcripts were higher in AR7 than in the WT line (Fig. 1A, Right). These results indicated that the gene impaired in AR7 encodes an important positive regulator of ABA-responsive gene expression in P. patens.

Fig. 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 1.

ABA-induced gene expression and SnRK2 activity in the WT line and AR7. (A) Results of comparative microarray analysis showing ABA up-regulated and down-regulated transcripts. (B) RNA gel blot analysis of WT and AR7 protonemata. Total RNA isolated from protonemata treated with 10 µM ABA for the indicated times was probed with radiolabeled cDNA fragments of the representative ABA-inducible transcripts, 17B9 (Pp1s294_52V6.1), 6A5 (Pp1s267_21V6.1), LEAIII (Pp1s118_232V6.1), and 19C6 (Pp1s52_261V6.1) (43). (C and D) In-gel kinase assays of WT and AR7 protonema for detection of SnRK2 activity. P. patens protonemata were treated with 10 µM ABA for the indicated times (C) or were subjected to various treatments (D). ABA, 10 µM ABA for 30 min; Cold, 0 °C for 1 h; Desic, desiccation in a laminar flow for 30 min; Man, 0.5 M mannitol for 30 min.

View this table:
  • View inline
  • View popup
Table S1.

Expression of 42 LEA genes detected in the microarray analysis

In angiosperms, most of the physiological processes required for the ABA response, including gene expression for various stress-associated proteins, are controlled by subclass III SnRK2s (10). ABA-induced activation of the protein kinase corresponding to the expected size of SnRK2 (39 kDa) has been demonstrated recently in P. patens protonemata by in-gel kinase assays using histone IIIS as a substrate (20). We found that protonemata of AR7 showed little activity of the ABA-induced kinase compared with WT protonemata (Fig. 1C). We further compared the kinase activity in the WT line and AR7 in response to other abiotic stresses and found that induction of the kinase activity by cold, desiccation, and hyperosmosis generated by 0.5 M mannitol treatment was also defective in AR7 (Fig. 1D). From these results, we speculated that the gene mutated in AR7 is involved in ABA-induced activation of SnRK2, which controls levels of the majority of transcripts during the response to exogenous ABA.

Identification of the Gene Mutated in AR7.

To identify the mutation locus in the genome of AR7, comparative genomic sequence analysis of AR7 and the WT line was carried out using the Illumina HiSeq 2000 sequencer. More than 286 million and 101 million reads of 100-bp nucleotide fragments of the WT line and AR7, covering roughly 58 and 21 times the genome size, respectively, were analyzed against the reference genome sequence of P. patens subsp. patens V1.6. The analysis revealed a total of 2,068 mutations in the AR7 genome, causing 47 nonsynonymous substitutions of amino acids in the encoded proteins (Dataset S1). To examine whether mutations of any of these 47 candidate genes caused ABA-insensitive phenotypes in AR7, the cDNA clones for these genes were isolated from the WT line and transiently expressed in the AR7 protonemata by particle bombardment. The candidate cDNA was cobombarded with the ABA-induced Em promoter fused to the β-glucuronidase reporter gene (Em-GUS) (23, 24) to determine whether any of the clones restore ABA-induced GUS expression in AR7, which otherwise shows very little GUS expression. Introduction of one cDNA clone corresponding to the gene model Pp1s462-10V6.1 was found to restore the ABA insensitivity of AR7. The Pp1s462-10V6.1 gene consisted of 17 exons, and analysis of the longest cDNA revealed that it encodes a 1,148-amino acid protein having a C-terminal protein kinase domain similar to plant Raf-like MAPKKKs belonging to B3-MAPKKKs (Fig. 2A). The stretch of amino acids from positions 86–185 in the N-terminal nonkinase domain contained a putative PAS domain possibly involved in molecular interactions (25), which typically is conserved in group B2 MAPKKKs (B2-MAPKKKs) categorized in angiosperms (26). The Pp1s462-10V6.1 gene was designated as “ARK.” Sequence analysis revealed that Ser-532 (TCC) located in a position N-terminal to the kinase domain was changed to Phe-532 (TTC) by a single C-to-T transition in the fifth exon of ARK in AR7 (Fig. 2A).

Fig. 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 2.

Characterization of Pp1s462-10V6.1 of P. patens. (A) Structures of the Pp1s462-10V6.1 gene mutated in AR7 and the encoded protein ARK deduced from the sequence of the isolated cDNA (GenBank KT581394). The C-to-T transition in the fifth exon caused an amino acid change from Ser to Phe at position 532 N-terminal to the protein kinase domain (KD). PD indicates a putative PAS domain. (B) Growth of protonemata with or without ABA of WT, AR7, the AR7 knock-in line with the ARK sequence changed to that of WT (AR7/ARKWT), and the WT knock-in line with the ARK sequence changed to that of AR7 type (WT/ARKS532F). (C) Transient assays for detection of ABA-induced gene expression. The protonemata of the lines used in B were bombarded with the Em-GUS and Ubi-LUC constructs and incubated with or without 10 µM ABA for 1 d, followed by the GUS and LUC assays. Levels of gene expression are represented by GUS/LUC ratio. Error bars indicate the SE (n = 3). ***P < 0.01 compared with ABA-treated WT (t-test). (D) ABA-induced freezing tolerance of the lines used in B. Protonemata were treated for 1 d with 10 µM ABA and either were kept at 4 °C (UF) or were subjected to equilibrium freezing to −3 °C or −4 °C. Survival of the protonemata was estimated by measuring electrolyte leakage from the injured cells. Error bars indicate the SE (n = 3). ***P < 0.001 and **P < 0.01 compared with the ABA-negative (−) tissue of the same line (t-test). (E) Osmotic stress tolerance of the protonemata. Pieces of protonemata were grown on BCDAT medium containing 0.4 M mannitol. (F) Cold-acclimation capacities of the protonemata. The cells were incubated on BCDAT medium at 0 °C for the indicated number of days and then were subjected to freezing at −3 °C to determine electrolyte leakage. Error bars indicate the SE (n = 3). *P < 0.05 and **P < 0.01 compared with the nonacclimated tissue of the same line (t-test). (G) Accumulation of stress-associated transcripts. RNA extracted from protonemata after treatment with ABA, cold, or hyperosmotic 0.5 M mannitol was electrophoresed, blotted onto a nylon membrane, and reacted with radiolabeled probes of the indicated transcripts. Ethidium bromide-stained rRNA is shown for verification of equal loading.

To determine why the Ser-to-Phe mutation at position 532 found in AR7 caused ABA insensitivity, we analyzed the accumulation of ARK proteins in the WT line and AR7 by immunoblot analysis using an antibody raised against the C-terminal amino acids of ARK. We detected immunoreactive bands of ∼120 kDa consistent with the length of 1,148 amino acids of ARK. However, the accumulation of ARK proteins appeared to be similar in the WT line and AR7 and did not change significantly during ABA treatment (Fig. S1). To confirm that the Ser-to-Phe mutation at position 532 caused the ABA-insensitive phenotype, a gene knock-in experiment was conducted. A genomic fragment of ARK from the WT line was fused with the neomycin phosphotransferase gene cassette so the marker would be inserted into a position downstream of the polyadenylation site. The resultant construct was introduced into AR7 by polyethylene glycol-mediated protoplast transformation to obtain transgenic plants resistant to the antibiotic G418. After confirmation of gene targeting to the ARK locus (Fig. S2), the generated knock-in line (AR7/ARKWT) was used for various physiological analyses of ABA responses. Exogenous ABA inhibited protonemal growth of AR7/ARKWT and facilitated subsequent formation of brood cells similar to that in the WT line (Fig. 2B). Transient expression analysis showed ABA-induced Em-GUS expression in the AR7/ARKWT protonemata similar to that in the WT line (Fig. 2C). Furthermore, impaired ABA-induced freezing tolerance of AR7 was recovered in AR7/ARKWT to a level similar to that in the WT line (Fig. 2D) concomitant with accumulation of LEA-like boiling-soluble proteins (Fig. S3).

Fig. S1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. S1.

(Upper) Immunoblot analysis of the WT line and AR7 for detection of ARK. Protonemata were incubated with 10 µM ABA for the indicated times and were homogenized in a buffer containing 50 mM Tris⋅Cl (pH 8.0), 150 mM NaCl, and 1% Triton X-100 with a 1/100 volume of proteinase inhibitor mixture (P9599; Sigma,) on ice. After centrifugation at 14,000 × g for 10 min at 4 °C, the supernatants were recovered and used for electrophoresis. The proteins (30 µg) were electrophoresed on a 6% (wt/vol) SDS-polyacrylamide gel, blotted onto a PVDF membrane, and reacted with affinity-purified anti-ARK antibody. (Lower) Coomassie Brilliant Blue staining of the large subunit of the ribulose 1,5-bisphosphate carboxylase/oxygenase.

Fig. S2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. S2.

Generation of the ARK knock-in line of P. patens. (A) Schematic representation of the gene knock-in construct of ARK for replacement of Phe (F) at the amino acid position 532 (aa532) with Ser (S). The gene cassette 35Sp-NPT-35Ster was inserted into the position 3′ to the polyadenylation site to select transformants on G418-containing medium. Arrows indicate positions of primers used for PCR analysis. (B) PCR analysis to confirm targeting using the ARK genome-specific primers (primers 1 and 2) and NPT primers (primers 3 and 4).

Fig. S3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. S3.

Accumulation of boiling-soluble proteins associated with ABA-induced freezing tolerance in protonemata of WT, AR7, and two knock-in lines, AR7/ARKWT and WT/ARKS532F. (A) Profiles of boiling-soluble proteins of protonemata cultured with or without 10 µM ABA for 1 d. Protonemata were homogenized in a solution containing 50 mM Tris⋅Cl (pH 7.5), 100 mM NaCl, 1 mM EDTA, and 1 mM DTT with a 1/100 volume of proteinase inhibitor mixture (P9599; Sigma) on ice. After centrifugation at 14,000 × g for 10 min at 4 °C, the supernatant was recovered. Total soluble proteins adjusted to equal concentrations were boiled for 1 min, and insoluble materials were removed by centrifugation at 14,000 × g for 10 min at 4 °C. These proteins in the supernatant were designated as boiling-soluble proteins. The boiling-soluble proteins corresponding to 20 μg of total soluble proteins were used for gel electrophoresis. The proteins were electrophoresed on 12% (wt/vol) SDS-polyacrylamide gels, and the gels were stained with silver. Sizes of molecular weight markers are shown in kilodaltons. (B) Immunoblot analysis of boiling-soluble proteins used in A. The proteins were reacted with an antibody against dehydrin, a group 2 LEA protein (PLA-100D; Enzo Life Sciences).

Hyperosmotic- and cold-stress responses also were examined in the AR7/ARKWT line. Although growth of AR7 protonemata was susceptible to 0.4 M mannitol contained in the medium, AR7/ARKWT protonemata grew on the medium in a manner similar to that of the WT line (Fig. 2E). Furthermore, AR7 showed reduced cold acclimation capacity (i.e., cold-induced freezing tolerance) compared with that of the WT line (22), but AR7/ARKWT acclimated to cold in a manner similar to that of the WT.

To confirm the phenotypic reversion by the change in a single amino acid of ARK, a similar gene knock-in experiment was carried out to introduce an AR7-type mutation (Ser-532 to Phe) into the ARK locus of the WT line. The resultant WT/ARKS532F line was found to mimic AR7 phenotypes in protonemata growth. In addition, that line was insensitive to ABA, cold, and hyperosmotic treatment as observed in AR7 (Fig. 2 A–F). When expression profiles of representative transcripts of WT, AR7, AR7/ARKWT, and WT/ARKS532F were compared in terms of responsiveness to ABA, hyperosmosis, and cold, the results were consistent with the growth and stress-tolerance phenotypes in these lines (Fig. 2G). These results indicated that ARK plays a crucial role in cellular responses to ABA and abiotic stresses in P. patens and that the single Ser-to-Phe change at the amino acid position 532 of ARK was sufficient to alter the responses.

ARK Kinase Activity Toward SnRK2.

To examine changes in ABA-induced protein kinase activity in the ARK knock-in lines, in-gel kinase assays were carried out using histone IIIS as a substrate. Protonemata of AR7/ARKWT had ABA-induced kinase activity similar to that of the WT line, whereas protonemata of the WT/ARKS532F line showed very little kinase activity, comparable to that of AR7 (Fig. 3A). From these results, we speculated that ARK is either directly or indirectly involved in the activation of SnRK2s. We found that PpSnRK2B (Pp1s240_91V6.1), one of the SnRK2s of P. patens, is activated by ABA, as shown by in-gel kinase assays of transgenic P. patens expressing PpSnRK2B fused to GFP (PpSnRK2B-GFP). PpSnRK2B also was activated by hyperosmosis (Fig. S4). Cellular localization of ARK-GFP and of PpSnRK2B-GFP was analyzed in P. patens protonema cells by confocal microscopy. ARK-GFP was localized mainly in the cytosol, and PpSnRK2B-GFP was localized in the cytosol and nucleus in the protonema cells, suggesting possible colocalization of PpSnRK2B with ARK in the cytosol (Fig. 3B). Similar experiments using the onion epidermal cells also indicated possible localization of PpSnRK2B and ARK in the cytosol (Fig. S5).

Fig. 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 3.

ARK as a possible regulator of SnRK2. (A, Upper) In-gel kinase assays of protonemata of WT, AR7, and ARK knock-in lines for the detection of ABA-activated protein kinases. Proteins were electrophoresed on an SDS-polyacrylamide gel containing histone IIIS. After denaturation and renaturation procedures, the proteins were reacted with γ-[32P]ATP. (Lower) Coomassie Brilliant Blue (CBB) staining of ribulose 1,5-bisphosphate carboxylase/oxygenase. (B) The ARK-GFP construct fused to the CaMV35 promoter was introduced into protonema cells of P. patens. Localization of GFP and PpSnRK2-GFP is shown for comparison. (C) Kinase activity of GST-ARKKD on MBP-PpSnRK2B. Recombinant proteins of MBP-PpSnRK2B (1 µg) (lanes 1 and 3) and GST-ARKKD (0.05 µg) (lanes 2 and 3) were reacted with γ-[32P]ATP and electrophoresed on an SDS-polyacrylamide gel, followed by Coomassie Brilliant Blue staining and autoradiography (32P). Closed and open arrowheads indicate the positions of MBP-PpSnRK2B and GST-ARKKD, respectively. (D) In-gel kinase assays of the recombinant MBP-PpSnRK2B (1 µg) with or without autophosphorylation (Auto-P) or phosphorylation by GST-ARKKD (0.05 µg). Positions of molecular-weight markers are shown in kilodaltons. Results of staining with Coomassie Brilliant Blue and immunoblot analysis using the anti-ARK antibody (Anti-ARK) are shown for comparison. Closed and open arrowheads indicate the positions of MBP-PpSnRK2B and GST-ARKKD, respectively. (E) Phosphorylation sites in the activation loop of MBP-PpSnRK2B identified by phosphopeptide mapping. The amino acid sequence and probable phosphorylation sites (asterisks) of the identified phosphopeptide (m/z = 924.46) and its ion-current chromatograms are shown. The MBP-PpSnRK2B protein was either phosphorylated with GST-ARKKD or autophosphorylated (Auto-P) for 15 min. (Details are shown in Fig. S6.)

Fig. S4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. S4.

Analysis of transgenic P. patens expressing PpSnRK2B-GFP. (A) Immunoblot analysis of control and transgenic P. patens expressing PpSnRK2B-GFP using the anti-GFP antibody. (B) In-gel kinase assay of protonemata of transgenic P. patens treated with ABA (1 or 10 µM, 30 min), mannitol (0.25 or 0.5 M, 30 min), desiccation (30 min), or cold (1 h). Electrophoresis was carried out using a polyacrylamide gel containing histone IIIS as a substrate. After denaturation and renaturation treatment, the proteins were reacted with γ-32P ATP. After washing, the gel was dried and exposed to X-ray film. Note that the film is overexposed for native kinase signals for better visualization of PpSnRK2B-GFP signals.

Fig. S5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. S5.

Cellular localization of ARK-GFP in onion epidermal cells. The ARK-GFP construct was fused to the cauliflower mosaic virus 35S promoter and bombarded into the cells of epidermal peel of the onion bulb. The results were compared with the localization with GFP alone and with PpSnRK2B-GFP.

To determine the role of ARK in the possible regulation of PpSnRK2B, phosphorylation assays were carried out using recombinant proteins. The protein kinase domain of ARK fused to GST (GST-ARKKD) and PpSnRK2B fused to maltose-binding protein (MBP-PpSnRK2B) were expressed in Escherichia coli, and affinity-purified recombinant proteins were used for in vitro phosphorylation assays. Results of assays indicated that MBP-PpSnRK2B showed some autophosphorylation activity, but GST-ARKKD phosphorylated MBP-PpSnRK2B more strongly (Fig. 3C). To examine the effect of this phosphorylation of MBP-PpSnRK2B on its kinase activity, in-gel kinase assays of nonphosphorylated, autophosphorylated, and ARK-phosphorylated MBP-PpSnRK2B were carried out using histone IIIS as a substrate. Results of assays indicated that the activity of MBP-PpSnRK2B was increased by phosphorylation catalyzed by GST-ARKKD but was not affected by autophosphorylation (Fig. 3D). Phosphopeptide mapping revealed that GST-ARKKD phosphorylated specific serine residues (Ser-165 and Ser-169) of the MBP-PpSnRK2B protein (Fig. 3E and Fig. S6). The phosphorylation site corresponding to Ser-169, which lies between subdomains VII and VIII in the activation loop, has been identified as the critical Ser residue for the activation of Arabidopsis SnRK2s (27).

Fig. S6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. S6.

Phosphorylation of PpSnRK2B analyzed with a nanoLC-MS/MS system. (A) MS/MS spectrum of the phosphopeptide (m/z = 924.46) of MBP-PpSnRK2B generated by GST-ARKKD phosphorylation. (B) Ion current chromatograms of the phosphopeptide from MBP-PpSnRK2B phosphorylated by GST-ARKKD (Upper) and total tryptic peptides (Lower). (C) Chromatograms of the peptide derived from autophosphorylated MBP-PpSnRK2B (Upper) and total tryptic peptides (Upper). (D–F) Results of a reference peptide (m/z = 691.87) by the analyses shown in A–C.

To explore possible regulation of ARK, we also examined the phosphorylation status of ARK. Proteins extracted from transgenic P. patens expressing ARK-GFP (Fig. 3B) were subjected to purification by magnetic beads with the anti-GFP antibody attached, and the purified proteins were used for phosphopeptide mapping. We identified phosphorylation of a peptide corresponding to amino acids 1029–1043 of ARK, with phosphorylation on either Ser-1029 or Ser-1030 (Fig. 4A). These residues were located within the activation loop in the kinase domain, comparable with phosphorylation sites of other eukaryotic MAPKKKs (Fig. 4B). To determine the role of phosphorylation of ARK, mutation constructs were generated by replacing both Ser-1029 and Ser-1030, or either, with Ala, and the constructs were used for transient assays in AR7 cells with Em-GUS to determine their capacity to restore ABA responsiveness (Fig. 4C). We found that ARK with Ala substitution of Ser-1029 (S1029A) does not restore ABA-induced Em-GUS expression in AR7. When both Ser-1029 and Ser-1030 were replaced with Ala (S1029A/S1030A), the GUS expression level was as low as in the Ser-1029 mutation, whereas replacing only Ser-1030 with Ala (S1030A) restored Em-GUS expression to a level similar to that in nonmutated ARK. These results indicate that phosphorylation of ARK at Ser-1029 in the activation loop is essential for its activation.

Fig. 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 4.

Phosphorylation in the activation loop of ARK is essential for ABA response. (A) MS/MS spectrum of the ARK-derived phosphopeptide found to be increased by ABA treatment in protonemata. Proteins were purified from ABA-treated protonemata of P. patens expressing ARK-GFP by using magnetic beads with the anti-GFP antibody attached. The purified proteins were subjected to digestion with Lys-C and trypsin and analysis by a nanoLC-MS-MS system. Either Ser-1029 or Ser-1030 was predicted to be the site of phosphorylation. (B) Sequence alignment of the amino acids in the activation loop of ARK and mammalian MAPKKKs. Predicted phosphorylation sites for ARK (S1029 and S1030), TAK1 (46), MLK3 (47), MEKK1 (48), and cRAF (49) are indicated by bold letters. (C) Transient assays of AR7 protonemata using constructs of ARK-GFP with or without mutations in the predicted phosphorylation sites. The cDNA was fused to the rice actin promoter and introduced into the AR7 cells with Em-GUS and Ubi-LUC, and the cells were incubated for 1 d with or without 10 µM ABA before protein extraction for the GUS and LUC assays. Levels of gene expression are represented by GUS/LUC ratio. Error bars indicate the SE (n = 3). *P < 0.05 compared with the ARK-GFP–introduced cells (t-test).

Functional Complementation of the Ark Mutation by B3-MAPKKKs from Other Land Plant Sources.

To determine the diversity of ARK-related genes in land plants, phylogenetic analysis of group B1–B4 Raf-like MAPKKKs from P. patens, A. thaliana, and Selaginella moellendorffii (Fig. 5A) was conducted. Three and six members of B3-MAPKKKs are encoded in the genomes of S. moellendorffii and A. thaliana, respectively, whereas ARK is the only B3-MAPKKK in P. patens. Furthermore, the S. moellendorffii and A. thaliana genomes encode two and six B2-MAPKKKs, respectively, but they were not found in the P. patens genome. The genes for group B1 MAPKKK (B1-MAPKKK) appeared to be the closest relatives of ARK in P. patens. To determine the functional relevance of the B1-MAPKKK genes to ARK, we examined whether the introduction of these B1-MAPKKK genes can restore ABA-induced GUS expression of AR7 by a transient assay. Introduction of two B1-MAPKKKs (Pp1s63_65V6.1 and Pp1s75_1V6.1) did not restore ABA-induced Em-GUS expression in AR7, indicating that the response is specific to ARK belonging to group B3 (Fig. 5B). We further tested whether B3-MAPKKK genes from A. thaliana and S. moellendorffii can complement ABA insensitivity of AR7 by transient assays. Fig. 5C shows that A. thaliana cDNA clones for At1g73660 and At4g24480 and the S. moellendorffii cDNA clone for Sm269874 restored ABA-induced Em-GUS expression in AR7. Results of experiments using other B3-MAPKKK clones indicated that At1g18160 also restored the GUS expression in AR7, but At5g11850, At5g03730, and At1g08720 had little effect (Fig. S7). These results indicated that specific genes of vascular plants encoding B3-MAPKKK can functionally complement the impairment of ARK in the moss.

Fig. 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 5.

Phylogenetic and functional analyses of group B3 Raf-like MAPKKKs in land plants. (A) Phylogenetic analysis of plant Raf-like MAPKKKs belonging to groups B1–B4 in P. patens, S. moellendorffii, and A. thaliana. (B and C) Transient assays of P. patens protonemata using various MAPKKK DNA clones. The cDNA clones were fused to the rice actin promoter and were cointroduced with Em-GUS and Ubi-LUC into protonemata. The protonemata were incubated with or without ABA for 1 d and were used for GUS and LUC assays. Error bars indicate the SE (n = 3). ***P < 0.001 and **P < 0.01 compared with the vector control (t-test). (B) Results for ARK and two group B1 MAPKKKs, Pp1s63_65V6.1 and Pp1s75_1V6.1. (C) Results for ARK, two Arabidopsis genes, At1g73660 and At4g24480, and one Selaginella gene, Sm269874.

Fig. S7.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. S7.

Transient assays of AR7 protonemata using B3-MAPKKKs of Arabidopsis. The cDNAs fused to the rice actin promoter were cointroduced with Em-GUS and Ubi-LUC into AR7protonemata. The protonemata were incubated with or without ABA for 1 d and were used for GUS and LUC assays. Error bars indicate the SE (n = 3). ***P < 0.001 and *P < 0.05 compared with the pFL vector control (t-test).

Discussion

Although subclass III SnRK2 is recognized as a key regulator of ABA signaling, the mode of kinase activity regulation has not been fully elucidated. A generally accepted mechanism of SnRK2 activation is that ABA suppresses PP2C activity via binding to PYR/PYL/RCAR receptors, thus facilitating enhancement of autophosphorylation of SnRK2 (11). In Arabidopsis, activity of SnRK2.2, SnRK2.3, and SnRK2.6 was increased in the triple mutant of group A PP2C (abi1hab1pp2ca), indicating that inhibition of the PP2C by the ABA-receptor complex is the key mechanism for activation of SnRK2 by autophosphorylation (28). The present study of ARK, a group B3 Raf-like kinase, provides evidence of positive regulation of ABA signaling in addition to negative regulation by group A PP2Cs in bryophytes. Plant Raf-like kinases have been shown to function in various types of biotic and abiotic signaling, but they are not thought to function in a typical MAP kinase pathway involving MAP kinase kinases and MAP kinases, and there is little biochemical evidence regarding their target substrates (29). This result is consistent with our observation of little ABA-induced kinase activity when myelin basic protein, a substrate commonly used for detection of MAP kinases, was used for in-gel kinase assays. Our study showing that a specific Raf-like kinase functions for activation of SnRK2 provides new insights into a signaling pathway for a possible connection between these kinases through unidentified mechanisms operating in plant ABA response.

Our study also indicated the possible mechanism for hyperosmosis-induced activation of plant SnRK2s mediated by an upstream protein kinase. SnRK2s representing all three subclasses in Arabidopsis and rice are activated by hyperosmosis (8, 10), but the mechanism of activation has not been elucidated. Results of a previous report indicated that hyperosmosis-induced SnRK2 activation might involve unidentified upstream protein kinases that phosphorylate SnRK2s (27). The reduced hyperosmosis-induced kinase activity in AR7 indicated that ARK might be involved in the activation of SnRK2 (Fig. 1C) toward osmostress tolerance (Fig. 2E). Hyperosmosis response by both ABA-dependent and independent mechanisms in the moss has been demonstrated recently using an ABA-deficient mutant of P. patens (30). Given that both mechanisms are mediated by ARK, this kinase is possibly a regulatory target of ABA and hyperosmotic signals upstream of SnRK2s.

Phosphopeptide mapping and mutational analysis of transient assays indicated that ARK itself might be activated by phosphorylation in the activation loop (Fig. 4 A and C). Activation by phosphorylation in the activation loop is well defined in mammalian and yeast MAPKKKs (Fig. 4B). Phosphorylation and activation events in these kinases typically are modulated through their nonkinase domain by dimerization, binding of other regulatory proteins, or phosphorylation by upstream MAPKKK kinases (31). The loss of ARK function with the Ser-to-Phe substitution at position 532 in the nonkinase domain (Fig. 2 A–G) suggests that the region around Ser-532 functions as a regulatory domain, and the Ser-to-Phe alteration might have resulted in a drastic change in microenvironment of the region because of the introduction of a bulky hydrophobic residue, leading to the abolishment of the kinase activity by an unidentified mechanism. The role of the N-terminal nonkinase domain in the regulation of kinase activity is well defined in Ste11p MAPKKK in budding yeast (32), wherein an intramolecular autoinhibitor sequence in the N-terminal domain interacts with the C-terminal kinase domain to inhibit its catalytic activity. The inhibition is enhanced by the binding of Ste50p to the N-terminal domain and is ameliorated upon signal perception by phosphorylation at the specific Ser and Thr residues by an upstream kinase Ste20p. The role of the N-terminal domain in stimulating ARK phosphorylation and how ABA facilitates this process have not been clarified.

ARK is the sole B3-MAPKKK in P. patens, but angiosperms have multiple B3-MAPKKKs with various functions. CONSTITUTIVE TRIPLE RESPONSE1 (CTR1, At5g03730) is well known as a negative regulator of ethylene signaling (33). ENHANCED DISEASE RESISTANCE1 (EDR1, At1g08720) and At1g73660 have been demonstrated to be involved in salicylic acid-induced defense responses (34) and in salt tolerance (35), respectively. Evidence of signal crosstalk between these apparently unrelated stimuli has been provided by analysis of Arabidopsis mutants. For instance, the ctr1 mutant allele has been isolated as a mutant that enhances ABA insensitivity of abi1-1 (36). The ctr1 allele also has been identified as a mutant defective in sugar response (37, 38). From studies using various other mutants, signal crosstalk among ABA, sugar, and ethylene in Arabidopsis has been proposed (39), although details of the mechanism have not been clarified. At1g73660, identified as the at6 locus that improves salt-stress tolerance (35), also has been identified as the gene responsible for sugar insensitive8 (sis8) showing sugar-resistant seedling growth (40, 41). The loss-of-function mutant of At1g73660, however, exhibited a WT response to ABA, the gibberellic acid biosynthesis inhibitor paclobutrazol, and the ethylene precursor aminocyclopropane carboxylic acid in germination experiments (35, 40), indicating that At1g73660 has a distinct function specific to salt and sugar. We hypothesize that the prototype B3-MAPKKK (ARK), mainly engaged in ABA and hyperosmosis responses in ancestral land plants, had diversified in vascular plants by acquiring more specialized functions toward sugar, salt, ethylene, and pathogenic signals. It should be noted that restoration of ABA-induced gene expression in AR7 by some of Arabidopsis B3-MAPKKKs (Fig. 5C) does not necessarily reflect physiological events during the ABA response in angiosperms. Whether any of the B3-MAPKKKs in angiosperms are involved in the positive regulation of SnRK2 would be clarified by detailed analysis of the mutants of B3-MAPKKKs with respect to ABA response and SnRK2 activation.

Experimental Procedures

Plant Materials, Growth Conditions, and Acclimation Treatment.

Protonema tissues of P. patens (Hedw.) Bruch & Schimp. were grown on cellophane-overlaid 0.8% agar plates of BCDAT medium supplemented with 0.5% glucose (42). The protonemata were cultured in a controlled environment growth chamber at 25 °C under continuous illumination (35 µmol photons⋅m−2⋅s−1). ABA, cold, and hyperosmotic treatments of the protonemata were carried out as described previously (22, 30).

Transcriptome Analysis.

Microarray analysis was carried out using a P. patens custom oligonucleotide microarray in which 33,942 gene-specific probes were mounted in an Agilent format microarray as described previously (20). Protonemata of the WT line and AR7 cultured on normal BCDAT medium or medium containing 10 µM ABA for 6 h were used for RNA extraction. One microgram of total RNA for each sample was used for RNA synthesis using the Amino Allyl MessageAmp II aRNA Amplification Kit (Life Technologies). The aRNA were labeled with Cy5 and Cy3 using the same kit. Fluorescent probes were hybridized on the array at 65 °C for 16 h by a dye-swap experimental procedure. The microarray was scanned with a GenePix 4000B system (Molecular Devices). Digitalization of signal intensities and calculation of the Cy5/Cy3 ratio were performed with Genepix Pro-6.0 software (Molecular Devices). In this process, quality control was carried out with signal values; probes with saturated intensity values (65,535) or a low signal-to-noise ratio (<3) were removed. Normalization was done by the Lowess (Locally Weighted Scatter plot Smoother) method. After filtering for intensity value, a total of 14,034 probes were subjected to further analysis. Fold-change analysis was performed on double replicate samples to determine the genes decreased in AR7. A twofold expression cutoff was considered to determine the genes decreased in AR7. ABA-regulated genes used for the analysis were described previously (20). The microarray data were deposited in Gene Expression Omnibus (accession GSE68914).

Transformation of P. patens.

A gene knock-in construct was made by amplification of the DNA fragment of genomic regions of ARK and fusion with the resistant cassette CaMV35S::NptII::NOS of pTN82 (a generous gift from Mitsuyasu Hasebe, National Institute for Basic Biology, Okazaki, Japan). Transformation of P. patens was carried out by the polyethylene glycol-mediated method as described previously (42).

Total RNA Preparation and RNA Gel Blot Analysis.

RNA extraction, electrophoresis, and blotting onto a nylon membrane were carried out according to the method described previously (22). The RNAs blotted on the membrane were hybridized with radiolabeled probes of 17B9 (Pp1s294_52V6.1), LEAIII (Pp1s118_232V6.1), Blt (Pp1s2_779V6.1), LEA172 (Pp1s370_29V6.1), and 6A5 (Pp1s267_21V6.1) (22, 30, 43). Hybridization was carried out using Ultrahyb hybridization buffer (Ambion) at 42 °C overnight. Posthybridization washings were carried out at 65 °C for 1 h in a buffer containing 0.2× SSC (1× SSC containing 150 mM NaCl and 15 mM Na citrate) and 0.2% SDS. Hybridizing signals were detected by exposure to X-ray film.

Comparative Genomic Analysis.

Sequencing libraries were prepared using the NEBNext DNA Sample Prep Reagent Set 1 (New England Biolabs) following the manufacturer’s protocols. Briefly, 3 µg of genomic DNA was fragmented to an average length of 300 bp by the Covaris S2 system (Covaris). After repair of the fragmented DNA, a single A nucleotide was ligated to the 3′ end. Illumina Index PE adaptors (Illumina) were ligated to the fragments, and fragments between 350 and 450 bp were selected with a 1.5% Pippin prep gel (Sage Science). The size-selected product was amplified by PCR for 12 cycles with primers InPE1.0, InPE2.0, and Index primer containing a unique index tag for each individual sample. The final product was validated by Agilent Bioanalyzer 2100 (Agilent). Pooled libraries were sequenced on HiSeq 2000 (Illumina) following the manufacturer’s protocol, generating 100-bp paired-end reads and 6-bp index tags. All FASTQ files were imported to CLC Genomics Workbench software ver. 5.5.1 (CLC bio). After the quality of each FASTQ file had been evaluated, reads were mapped to the reference genome of P. patens (downloaded from genome.jgi.doe.gov/PhytozomeV9/download/_JAMO/52b9c81b166e730e43a34f83/Ppatens_152.fa.gz) using default parameters. Variants were detected using the probabilistic variant detection function with following parameters: minimum coverage = 5, maximum expected variants = 1, ignore quality scores = no, ignore nonspecific matches = yes, ignore broken pairs = no, filter 454/Ion homopolymer indels = no, variant probability = 90.0, require presence in both forward and reverse reads = yes. Nonsynonymous substitutions were detected using the “Amino Acid Changes” tool of CLC Genomics Workbench with reference to the P. patens ver.1.6 coding sequences.

Transient Assays of P. patens Protonemata.

Protonemata of P. patens cultured on the BCDAT agar medium were bombarded with Em-GUS and Ubi-LUC by using the PDS-1000He particle delivery system (Bio-Rad) with 1-µm gold particles and 1,100-psi rupture discs under a vacuum of 26 mmHg. The bombarded protonemata were cultured at 25 °C under continuous light for 24 h on BCDAT agar medium with or without 10 µM ABA. Proteins extracted from the protonemata were used for GUS and LUC assays to determine the GUS/LUC ratio as previously described (44).

Tests for Freezing Tolerance and Hyperosmotic-Stress Tolerance.

The freezing tolerance of P. patens protonemata was determined as described previously (30). In brief, protonema tissues placed in a glass test tube were seeded with ice at −1 °C, kept at −1 °C for 50 min, and then cooled at a rate of −2.4 °C/h to desired temperatures. Electrolyte leakage from the damaged tissues was measured after thawing by using a conductivity meter and was represented as the percentage of total ions released after subsequent boiling (1). For osmotic stress treatment, the tissues were placed on a medium containing 0.4 M mannitol and were grown for 7 d under continuous light.

Analysis Using GFP-Fusion Constructs.

GFP-fusion constructs of ARK and PpSnRK2B were generated by fusing the cDNA to CaMV35S::sGFP(S65T) (45) using the NcoI site at start codon of GFP. Five-day-old P. patens protonemata were bombarded with the GFP-fusion constructs, and cells were observed under a confocal laser-scanning microscope.

Preparation of Recombinant Proteins.

For preparation of MBP-PpSnRKB, the entire coding sequence of PpSnRKB cDNA (Pp1s240_91V6.1) was amplified from reverse transcripts of P. patens protonemata and fused in-frame to the NdeI-digested pMAL-c5X vector (New England Biolabs). The MBP-fusion protein was purified using amylose resin according to the manufacturer’s instructions. For GST-ARKKD, the region of the protein kinase domain (amino acids 867–1148) was fused in-frame to the BamHI/NotI-digested pGEX-5 × 3 vector (GE Healthcare). The GST-fusion protein was purified using glutathione-Sepharose resin. Both recombinant proteins were subjected to ultrafiltration using a Nanosep 30-kDa size-exclusion column (Pall) for concentration and removal of low-molecular-weight materials.

Immunoblot Analysis.

Proteins electrophoresed on 10% (wt/vol) SDS-polyacrylamide gels were transferred onto PVDF membranes, reacted with primary and secondary antibodies, and detected by either nitroblue tetrazolium/5-bromo-4-chloro-3-indolyl phosphate reagent or chemiluminescent reagents. For preparation of ARK antibodies, antiserum raised against the C-terminal 15 amino acids (LGGTPKSGLSDRDL) of ARK was subjected to affinity purification using the recombinant GST-ARKKD protein fixed on a PVDF membrane.

Protein Kinase Assays.

In vitro phosphorylation reactions were carried out using affinity-purified recombinant proteins at 30 °C for 15 min in 50 mM Tris⋅Cl (pH 7.5), 10 mM MgCl2, 1 mM DTT, and 50 µM γ-[32P]ATP. The reaction mixture was electrophoresed on an SDS-polyacrylamide gel and stained with Coomassie Brilliant Blue, and phosphoproteins were detected by autoradiography. In-gel kinase assays were carried out essentially as described previously (20). In brief, 50 µg of total soluble proteins extracted from ABA-treated and nontreated protonemata were electrophoresed using SDS-polyacrylamide gel containing histone IIIS (H5055; Sigma). After denaturation and renaturation procedures, the gel was incubated with 50 µM γ-32P ATP. The gel was washed, dried, and exposed to X-ray film to detect kinase activity. For recombinant proteins, the kinase reaction was carried out at 30 °C for 15 min in 50 mM Tris⋅Cl (pH 7.5), 10 mM MgCl2, and 1 mM DTT with 50 µM nonradioactive ATP, and the reaction mixture was loaded onto the histone IIIS-polymerized polyacrylamide gel for in-gel kinase assays.

Phosphopeptide Mapping.

Protonemata of P. patens expressing ARK-GFP were used for protein extraction. The ARK-GFP proteins were purified by using anti-GFP magnetic microbeads (Miltenyi Biotech) according to the manufacturer's instructions. Determination of phosphorylation sites of proteins was carried out essentially as described previously (12). Phosphorylated proteins digested with Lys-C and trypsin were analyzed by a nanoLC-MS-MS system for peptide identification, quantification, and prediction of phosphorylation sites. For enrichment of phosphopeptides, the digested peptides were subjected to hydroxyl acid-associated metal oxide affinity chromatography.

Acknowledgments

We thank Ralph Quatrano for Em-GUS, Tuan-hua David Ho for Ubi-LUC, Yasuo Niwa for 35S-sGFP(S65T), Mitsuyasu Hasebe for pTN82, and the NODAI Genome Research Center of Tokyo University of Agriculture for the use of Illumina HiSeq 2000. This study was supported in part by Grants-in-Aid for Scientific Research 26291054, 23119504, and 25119705 from the Ministry of Education, Culture, Sports, Science, and Technology of Japan (to D.T.), 24570058 and 16780235) (to Y. Sakata), and 25840103 and 15H04383 (to T.U.); by Grant for Basic Science Research Projects 110493 from The Sumitomo Foundation (to Y. Sakata); and by Precursory Research for Embryonic Science and Technology (PRESTO) from Japan Science and Technology Agency (T.U.).

Footnotes

  • ↵1To whom correspondence should be addressed. Email: takezawa{at}mail.saitama-u.ac.jp.
  • Author contributions: Y. Sakata and D.T. designed research; M.S., T.K.G., K.A., Y.I., K.H., K.K., Y. Shiwa, K.I., H.Y., T.U., Y. Sakata, and D.T. performed research; M.S., T.K.G., Y.I., K.K., Y. Shiwa, K.I., H.Y., T.U., Y. Sakata, and D.T. analyzed data; and M.S. and D.T. wrote the paper.

  • The authors declare no conflict of interest.

  • This article is a PNAS Direct Submission.

  • Data deposition: The sequence reported in this paper has been deposited in the GenBank database (accession no. KT581394), and the microarray data have been deposited in the Gene Expression Omnibus database (accession no. GSE68914).

  • This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1511238112/-/DCSupplemental.

Freely available online through the PNAS open access option.

References

  1. ↵
    1. Nagao M,
    2. Minami A,
    3. Arakawa K,
    4. Fujikawa S,
    5. Takezawa D
    (2005) Rapid degradation of starch in chloroplasts and concomitant accumulation of soluble sugars associated with ABA-induced freezing tolerance in the moss Physcomitrella patens. J Plant Physiol 162(2):169–180
    .
    OpenUrlCrossRefPubMed
  2. ↵
    1. Cuming AC,
    2. Cho SH,
    3. Kamisugi Y,
    4. Graham H,
    5. Quatrano RS
    (2007) Microarray analysis of transcriptional responses to abscisic acid and osmotic, salt, and drought stress in the moss, Physcomitrella patens. New Phytol 176(2):275–287
    .
    OpenUrlCrossRefPubMed
  3. ↵
    1. Khandelwal A, et al.
    (2010) Role of ABA and ABI3 in desiccation tolerance. Science 327(5965):546
    .
    OpenUrlAbstract/FREE Full Text
  4. ↵
    1. Li J,
    2. Wang XQ,
    3. Watson MB,
    4. Assmann SM
    (2000) Regulation of abscisic acid-induced stomatal closure and anion channels by guard cell AAPK kinase. Science 287(5451):300–303
    .
    OpenUrlAbstract/FREE Full Text
  5. ↵
    1. Assmann SM
    (2003) OPEN STOMATA1 opens the door to ABA signaling in Arabidopsis guard cells. Trends Plant Sci 8(4):151–153
    .
    OpenUrlCrossRefPubMed
  6. ↵
    1. Umezawa T, et al.
    (2010) Molecular basis of the core regulatory network in ABA responses: Sensing, signaling and transport. Plant Cell Physiol 51(11):1821–1839
    .
    OpenUrlAbstract/FREE Full Text
  7. ↵
    1. Fujita Y,
    2. Yoshida T,
    3. Yamaguchi-Shinozaki K
    (2013) Pivotal role of the AREB/ABF-SnRK2 pathway in ABRE-mediated transcription in response to osmotic stress in plants. Physiol Plant 147(1):15–27
    .
    OpenUrlCrossRef
  8. ↵
    1. Kobayashi Y,
    2. Yamamoto S,
    3. Minami H,
    4. Kagaya Y,
    5. Hattori T
    (2004) Differential activation of the rice sucrose nonfermenting1-related protein kinase2 family by hyperosmotic stress and abscisic acid. Plant Cell 16(5):1163–1177
    .
    OpenUrlAbstract/FREE Full Text
  9. ↵
    1. Umezawa T,
    2. Yoshida R,
    3. Maruyama K,
    4. Yamaguchi-Shinozaki K,
    5. Shinozaki K
    (2004) SRK2C, a SNF1-related protein kinase 2, improves drought tolerance by controlling stress-responsive gene expression in Arabidopsis thaliana. Proc Natl Acad Sci USA 101(49):17306–17311
    .
    OpenUrlAbstract/FREE Full Text
  10. ↵
    1. Boudsocq M,
    2. Barbier-Brygoo H,
    3. Laurière C
    (2004) Identification of nine sucrose nonfermenting 1-related protein kinases 2 activated by hyperosmotic and saline stresses in Arabidopsis thaliana. J Biol Chem 279(40):41758–41766
    .
    OpenUrlAbstract/FREE Full Text
  11. ↵
    1. Fujii H, et al.
    (2009) In vitro reconstitution of an abscisic acid signalling pathway. Nature 462(7273):660–664
    .
    OpenUrlCrossRefPubMed
  12. ↵
    1. Umezawa T, et al.
    (2009) Type 2C protein phosphatases directly regulate abscisic acid-activated protein kinases in Arabidopsis. Proc Natl Acad Sci USA 106(41):17588–17593
    .
    OpenUrlAbstract/FREE Full Text
  13. ↵
    1. Yoshida R, et al.
    (2006) The regulatory domain of SRK2E/OST1/SnRK2.6 interacts with ABI1 and integrates abscisic acid (ABA) and osmotic stress signals controlling stomatal closure in Arabidopsis. J Biol Chem 281(8):5310–5318
    .
    OpenUrlAbstract/FREE Full Text
  14. ↵
    1. Ma Y, et al.
    (2009) Regulators of PP2C phosphatase activity function as abscisic acid sensors. Science 324(5930):1064–1068
    .
    OpenUrlAbstract/FREE Full Text
  15. ↵
    1. Park S-Y, et al.
    (2009) Abscisic acid inhibits type 2C protein phosphatases via the PYR/PYL family of START proteins. Science 324(5930):1068–1071
    .
    OpenUrlAbstract/FREE Full Text
  16. ↵
    1. Soon F-F, et al.
    (2012) Molecular mimicry regulates ABA signaling by SnRK2 kinases and PP2C phosphatases. Science 335(6064):85–88
    .
    OpenUrlAbstract/FREE Full Text
  17. ↵
    1. Rensing SA, et al.
    (2008) The Physcomitrella genome reveals evolutionary insights into the conquest of land by plants. Science 319(5859):64–69
    .
    OpenUrlAbstract/FREE Full Text
  18. ↵
    1. Sakata Y,
    2. Komatsu K,
    3. Takezawa D
    (2014) ABA as a universal plant hormone. Prog Bot 75:57–96
    .
    OpenUrlCrossRef
  19. ↵
    1. Chater C, et al.
    (2011) Regulatory mechanism controlling stomatal behavior conserved across 400 million years of land plant evolution. Curr Biol 21(12):1025–1029
    .
    OpenUrlCrossRefPubMed
  20. ↵
    1. Komatsu K, et al.
    (2013) Group A PP2Cs evolved in land plants as key regulators of intrinsic desiccation tolerance. Nat Commun 4:2219
    .
    OpenUrlPubMed
  21. ↵
    1. Minami A,
    2. Togawa S,
    3. Nagao M,
    4. Takezawa D
    (2006) Altered freezing tolerance in the Physcomitrella patens mutant with reduced sensitivity to abscisic acid. Cryobiol Cryotechnol 52(2):135–139
    .
    OpenUrl
  22. ↵
    1. Bhyan SB, et al.
    (2012) Cold acclimation in the moss Physcomitrella patens involves abscisic acid-dependent signaling. J Plant Physiol 169(2):137–145
    .
    OpenUrlCrossRefPubMed
  23. ↵
    1. Marcotte WR Jr,
    2. Russell SH,
    3. Quatrano RS
    (1989) Abscisic acid-responsive sequences from the em gene of wheat. Plant Cell 1(10):969–976
    .
    OpenUrlAbstract/FREE Full Text
  24. ↵
    1. Knight CD, et al.
    (1995) Molecular responses to abscisic acid and stress are conserved between moss and cereals. Plant Cell 7(5):499–506
    .
    OpenUrlAbstract/FREE Full Text
  25. ↵
    1. Hefti MH,
    2. Françoijs K-J,
    3. de Vries SC,
    4. Dixon R,
    5. Vervoort J
    (2004) The PAS fold. A redefinition of the PAS domain based upon structural prediction. Eur J Biochem 271(6):1198–1208
    .
    OpenUrlPubMed
  26. ↵
    1. Ichimura K, et al., MAPK Group
    (2002) Mitogen-activated protein kinase cascades in plants: A new nomenclature. Trends Plant Sci 7(7):301–308
    .
    OpenUrlCrossRefPubMed
  27. ↵
    1. Boudsocq M,
    2. Droillard MJ,
    3. Barbier-Brygoo H,
    4. Laurière C
    (2007) Different phosphorylation mechanisms are involved in the activation of sucrose non-fermenting 1 related protein kinases 2 by osmotic stresses and abscisic acid. Plant Mol Biol 63(4):491–503
    .
    OpenUrlCrossRefPubMed
  28. ↵
    1. Ng LM, et al.
    (2011) Structural basis for basal activity and autoactivation of abscisic acid (ABA) signaling SnRK2 kinases. Proc Natl Acad Sci USA 108(52):21259–21264
    .
    OpenUrlAbstract/FREE Full Text
  29. ↵
    1. Menges M, et al.
    (2008) Comprehensive gene expression atlas for the Arabidopsis MAP kinase signalling pathways. New Phytol 179(3):643–662
    .
    OpenUrlCrossRefPubMed
  30. ↵
    1. Takezawa D, et al.
    (2015) Epoxycarotenoid-mediated synthesis of abscisic acid in Physcomitrella patens implicating conserved mechanisms for acclimation to hyperosmosis in embryophytes. New Phytol 206(1):209–219
    .
    OpenUrlCrossRefPubMed
  31. ↵
    1. Kolch W
    (2000) Meaningful relationships: The regulation of the Ras/Raf/MEK/ERK pathway by protein interactions. Biochem J 351(Pt 2):289–305
    .
    OpenUrlAbstract/FREE Full Text
  32. ↵
    1. Drogen F, et al.
    (2000) Phosphorylation of the MEKK Ste11p by the PAK-like kinase Ste20p is required for MAP kinase signaling in vivo. Curr Biol 10(11):630–639
    .
    OpenUrlCrossRefPubMed
  33. ↵
    1. Kieber JJ,
    2. Rothenberg M,
    3. Roman G,
    4. Feldmann KA,
    5. Ecker JR
    (1993) CTR1, a negative regulator of the ethylene response pathway in Arabidopsis, encodes a member of the raf family of protein kinases. Cell 72(3):427–441
    .
    OpenUrlCrossRefPubMed
  34. ↵
    1. Frye CA,
    2. Innes RW
    (1998) An Arabidopsis mutant with enhanced resistance to powdery mildew. Plant Cell 10(6):947–956
    .
    OpenUrlAbstract/FREE Full Text
  35. ↵
    1. Gao L,
    2. Xiang CB
    (2008) The genetic locus At1g73660 encodes a putative MAPKKK and negatively regulates salt tolerance in Arabidopsis. Plant Mol Biol 67(1-2):125–134
    .
    OpenUrlCrossRefPubMed
  36. ↵
    1. Beaudoin N,
    2. Serizet C,
    3. Gosti F,
    4. Giraudat J
    (2000) Interactions between abscisic acid and ethylene signaling cascades. Plant Cell 12(7):1103–1115
    .
    OpenUrlAbstract/FREE Full Text
  37. ↵
    1. Zhou L,
    2. Jang JC,
    3. Jones TL,
    4. Sheen J
    (1998) Glucose and ethylene signal transduction crosstalk revealed by an Arabidopsis glucose-insensitive mutant. Proc Natl Acad Sci USA 95(17):10294–10299
    .
    OpenUrlAbstract/FREE Full Text
  38. ↵
    1. Gibson SI,
    2. Laby RJ,
    3. Kim D
    (2001) The sugar-insensitive1 (sis1) mutant of Arabidopsis is allelic to ctr1. Biochem Biophys Res Commun 280(1):196–203
    .
    OpenUrlCrossRefPubMed
  39. ↵
    1. Rolland F,
    2. Baena-Gonzalez E,
    3. Sheen J
    (2006) Sugar sensing and signaling in plants: Conserved and novel mechanisms. Annu Rev Plant Biol 57:675–709
    .
    OpenUrlCrossRefPubMed
  40. ↵
    1. Laby RJ,
    2. Kincaid MS,
    3. Kim D,
    4. Gibson SI
    (2000) The Arabidopsis sugar-insensitive mutants sis4 and sis5 are defective in abscisic acid synthesis and response. Plant J 23(5):587–596
    .
    OpenUrlCrossRefPubMed
  41. ↵
    1. Huang Y,
    2. Li CY,
    3. Qi Y,
    4. Park S,
    5. Gibson SI
    (2014) SIS8, a putative mitogen-activated protein kinase kinase kinase, regulates sugar-resistant seedling development in Arabidopsis. Plant J 77(4):577–588
    .
    OpenUrlCrossRefPubMed
  42. ↵
    1. Nishiyama T,
    2. Hiwatashi Y,
    3. Sakakibara I,
    4. Kato M,
    5. Hasebe M
    (2000) Tagged mutagenesis and gene-trap in the moss, Physcomitrella patens by shuttle mutagenesis. DNA Res 7(1):9–17
    .
    OpenUrlAbstract
  43. ↵
    1. Chen T,
    2. Uemura M
    1. Minami A,
    2. Nagao M,
    3. Arakawa K,
    4. Fujikawa S,
    5. Takezawa D
    (2006) Physiological and morphological alteration associated with freezing tolerance in the moss Physcomitrella patens. Cold Hardiness in Plants: Molecular Genetics, Cell Biology and Physiology, eds Chen T, Uemura M (CABI, Wallingford), pp 138–152
    .
  44. ↵
    1. Komatsu K, et al.
    (2009) Functional analyses of the ABI1-related protein phosphatase type 2C reveal evolutionarily conserved regulation of abscisic acid signaling between Arabidopsis and the moss Physcomitrella patens. Plant Mol Biol 70(3):327–340
    .
    OpenUrlCrossRefPubMed
  45. ↵
    1. Niwa Y
    (2003) A synthetic green fluorescent protein gene for plant biotechnology. Plant Biotechnol 20(1):1–11
    .
    OpenUrlCrossRef
  46. ↵
    1. Kishimoto K,
    2. Matsumoto K,
    3. Ninomiya-Tsuji J
    (2000) TAK1 mitogen-activated protein kinase kinase kinase is activated by autophosphorylation within its activation loop. J Biol Chem 275(10):7359–7364
    .
    OpenUrlAbstract/FREE Full Text
  47. ↵
    1. Leung IW,
    2. Lassam N
    (2001) The kinase activation loop is the key to mixed lineage kinase-3 activation via both autophosphorylation and hematopoietic progenitor kinase 1 phosphorylation. J Biol Chem 276(3):1961–1967
    .
    OpenUrlAbstract/FREE Full Text
  48. ↵
    1. Deak JC,
    2. Templeton DJ
    (1997) Regulation of the activity of MEK kinase 1 (MEKK1) by autophosphorylation within the kinase activation domain. Biochem J 322(Pt 1):185–192
    .
    OpenUrlAbstract/FREE Full Text
  49. ↵
    1. Kolch W, et al.
    (1993) Protein kinase C α activates RAF-1 by direct phosphorylation. Nature 364(6434):249–252
    .
    OpenUrlCrossRefPubMed
PreviousNext
Back to top
Article Alerts
Email Article

Thank you for your interest in spreading the word on PNAS.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Plant Raf-like kinase integrates abscisic acid and hyperosmotic stress signaling upstream of SNF1-related protein kinase2
(Your Name) has sent you a message from PNAS
(Your Name) thought you would like to see the PNAS web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Plant Raf-like kinase for ABA and stress signaling
Masashi Saruhashi, Totan Kumar Ghosh, Kenta Arai, Yumiko Ishizaki, Kazuya Hagiwara, Kenji Komatsu, Yuh Shiwa, Keiichi Izumikawa, Harunori Yoshikawa, Taishi Umezawa, Yoichi Sakata, Daisuke Takezawa
Proceedings of the National Academy of Sciences Nov 2015, 112 (46) E6388-E6396; DOI: 10.1073/pnas.1511238112

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Request Permissions
Share
Plant Raf-like kinase for ABA and stress signaling
Masashi Saruhashi, Totan Kumar Ghosh, Kenta Arai, Yumiko Ishizaki, Kazuya Hagiwara, Kenji Komatsu, Yuh Shiwa, Keiichi Izumikawa, Harunori Yoshikawa, Taishi Umezawa, Yoichi Sakata, Daisuke Takezawa
Proceedings of the National Academy of Sciences Nov 2015, 112 (46) E6388-E6396; DOI: 10.1073/pnas.1511238112
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Mendeley logo Mendeley

Article Classifications

  • Biological Sciences
  • Plant Biology
Proceedings of the National Academy of Sciences: 112 (46)
Table of Contents

Submit

Sign up for Article Alerts

Jump to section

  • Article
    • Abstract
    • Results
    • Discussion
    • Experimental Procedures
    • Acknowledgments
    • Footnotes
    • References
  • Figures & SI
  • Info & Metrics
  • PDF

You May Also be Interested in

Reflection of clouds in the still waters of Mono Lake in California.
Inner Workings: Making headway with the mysteries of life’s origins
Recent experiments and simulations are starting to answer some fundamental questions about how life came to be.
Image credit: Shutterstock/Radoslaw Lecyk.
Depiction of the sun's heliosphere with Voyager spacecraft at its edge.
News Feature: Voyager still breaking barriers decades after launch
Launched in 1977, Voyagers 1 and 2 are still helping to resolve past controversies even as they help spark a new one: the true shape of the heliosphere.
Image credit: NASA/JPL-Caltech.
Drop of water creates splash in a puddle.
Journal Club: Heavy water tastes sweeter
Heavy hydrogen makes heavy water more dense and raises its boiling point. It also appears to affect another characteristic long rumored: taste.
Image credit: Shutterstock/sl_photo.
Mouse fibroblast cells. Electron bifurcation reactions keep mammalian cells alive.
Exploring electron bifurcation
Jonathon Yuly, David Beratan, and Peng Zhang investigate how electron bifurcation reactions work.
Listen
Past PodcastsSubscribe
Panda bear hanging in a tree
How horse manure helps giant pandas tolerate cold
A study finds that giant pandas roll in horse manure to increase their cold tolerance.
Image credit: Fuwen Wei.

Similar Articles

Site Logo
Powered by HighWire
  • Submit Manuscript
  • Twitter
  • Facebook
  • RSS Feeds
  • Email Alerts

Articles

  • Current Issue
  • Special Feature Articles – Most Recent
  • List of Issues

PNAS Portals

  • Anthropology
  • Chemistry
  • Classics
  • Front Matter
  • Physics
  • Sustainability Science
  • Teaching Resources

Information

  • Authors
  • Editorial Board
  • Reviewers
  • Subscribers
  • Librarians
  • Press
  • Cozzarelli Prize
  • Site Map
  • PNAS Updates
  • FAQs
  • Accessibility Statement
  • Rights & Permissions
  • About
  • Contact

Feedback    Privacy/Legal

Copyright © 2021 National Academy of Sciences. Online ISSN 1091-6490