Proteomic Analysis of Phosphoproteins in the Rice Nucleus During the

Jun 2, 2015 - The early stage of seed germination is the first step in the plant life cycle without visible morphological change. To investigate the m...
0 downloads 9 Views 2MB Size
Subscriber access provided by NEW YORK UNIV

Article

Proteomic Analysis of Phosphoproteins in the Rice Nucleus during Early Stage of Seed Germination Ming Li, Xiaojian Yin, Katsumi Sakata, Pingfang Yang, and Setsuko Komatsu J. Proteome Res., Just Accepted Manuscript • Publication Date (Web): 02 Jun 2015 Downloaded from http://pubs.acs.org on June 3, 2015

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Journal of Proteome Research is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 50

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

Proteomic Analysis of Phosphoproteins in the Rice Nucleus during Early Stage of Seed Germination

Ming Li1, 2,#, Xiaojian Yin2,#, Katsumi Sakata3, Pingfang Yang1, and Setsuko Komatsu2,*

1 Key Laboratory of Plant Germplasm Enhancement and Specialty Agriculture, Wuhan Botanical Garden, Chinese Academy of Sciences, Wuhan 430074, China 2 National Institute of Crop Science, National Agriculture and Food Research Organization, Tsukuba 305-8518, Japan 3 Department of Life Science and Informatics, Maebashi Institute of Technology, Maebashi 371-0816, Japan

Corresponding Author *Setsuko Komatsu, National Institute of Crop Science, National Agriculture and Food Research Organization, Tsukuba 305-8518, Japan. Tel: +81-29-838-8693. Fax: +81-29-838-8694. E-mail: [email protected].

Running title: Nuclear phosphoproteomics during rice germination

ABBREVIATIONS ABA, abscisic acid; GAs, gibberellins; LC, liquid chromatography; MS, mass spectrometry; PolyMAC, Polymer-based Metal-ion Affinity Capture; qRT-PCR, quantitative reverse transcription polymerase chain reaction

1 ACS Paragon Plus Environment

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ABSTRACT: The early stage of seed germination is the first step in the plant life cycle without visible morphological change. To investigate the mechanism controlling the early stage of rice seed germination, a gel-free/label-free nuclear phosphoproteomics was performed. A total of 3,467 phosphopeptides belonging to 102 nuclear phosphoproteins of rice embryos were identified. Protein synthesis related proteins were mainly phosphorylated. During the first 24 h following imbibition, 115 nuclear phosphoproteins were identified and significant changes in phosphorylation level over time were observed in 29 phosphoproteins. Cluster analysis indicated that nucleotide binding proteins and zinc finger CCCH/BED type proteins increased in abundance during the first 12 h of imbibition and then decreased. The in silico protein-protein interactions for 29 nuclear phosphoproteins indicated that Sas10/Utp3 protein, which functions in snoRNA bind and gene silence, was the center of the phosphoproteins network in nuclei. The germination rate of seeds was significantly slow down with phosphatase-inhibitor treatment. The mRNA expression of zinc finger CCCH type protein did not change and zinc finger BED type was up-regulated in rice embryos during early stage of germination with phosphatase-inhibitor treatment. These results suggest that phosphorylation and dephosphorylation of nuclear protein involve in rice seed germination. Furthermore, transcription factors such as zinc finger CCCH/BED type protein might play a key role through nuclear phosphoproteins, and Sas10/Utp3 protein might be interacted with nuclear phosphoproteins in rice embryos to mediate the early stage of seed germination.

KEYWORDS: rice, nucleus, phosphoproteomics, germination, phosphonetwork

2 ACS Paragon Plus Environment

Page 2 of 50

Page 3 of 50

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

INTRODUCTION For many species of plants, seeds are necessary to reproduce and spread offspring, and seed production is also a good strategy to overcome severe environmental conditions 1. Seed germination, which usually starts when water is imbibed by the quiescent dry seed and ends with the elongation of the radicle, is a key stage in the plant life cycle 2,3. Based on the increase in seed weight, rice seed germination is divided into three phases: first, rapid water uptake without visible morphological change (Phase I), followed by a plateau stage (Phase II) without any apparent change in weight, and finally, rapid water uptake (Phase III) accompanied by growth of the radical and coleoptile 3. Transcriptomic analysis of rice indicated that the switch involved in germination may happen in Phase I, with the onset of mRNA biosynthesis 4. Among three phases, Phase I is the important stage during seed germination because mRNA biosynthesis starts at this stage. To understand the mechanism of rice seed germination, a proteomic technique was used. Exploring gibberellins (GAs) and abscisic acid (ABA) responsive functional proteins such as DELLA in the GA and ABA signaling pathways suggested many proteins are involved in GA and ABA modulation during seed germination 5,6. Profiling of rice seed protein abundance during germination revealed comprehensive metabolic and regulatory pathways during seed germination 7. Study of embryos revealed that starch biosynthetic enzymes play major roles in rice seed germination 8. Investigation of coleoptiles indicated that they respond to anoxia dramatically, with an increase in abundance of proteins related to stress responses and redox metabolism 9-11. Analysis of the aleuronic layer has demonstrated that thioredoxin activates embryo-specific protein 2 with a concurrent unfolding of its substrate during rice seed germination 12-14. Because the function of a protein is largely correlated with its cellular location and posttranslational modifications, research on subcellular proteomics and posttranslational 3 ACS Paragon Plus Environment

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

modifications might be helpful in understanding the mechanism controlling rice seed germination. In order to clarify this mechanism, mitochondrial proteomics and phosphoproteomics studies have been performed during germination of rice seeds. In rice seed mitochondrial proteomics, it was found that mitochondria matured during germination in anaerobic conditions, and respiratory functionality was repressed by reduced protein abundances of tricarboxylic acid cycle components and cytochrome-containing complexes of the respiratory chain 15,16. In rice seed phosphoproteomics, it was found that brassinosteroid signaling might participate in the triggering of rice seed germination 17,18. However, proteomic studies of other organelles or posttranslational modifications are still lacking. The nucleus is responsible for storing most genetic information and controlling the plant's activities 19. All nuclear proteins are synthesized in the cytoplasm and imported into the nucleus to execute their functions 20,21. Furthermore, after the plasma membrane perceives changes in the extracellular environment, signals must be communicated to the nucleus to induce or suppress gene transcription 21. Phosphorylation plays a key role in the signalling cascade through kinases and phosphatases to activate or inactivate enzyme activity in biological pathways 22. As a reversible protein modification, phosphorylation enhances or inhibits proteins transport between nucleus and cytoplasm 23. Some studies on regulation and activation of nuclear localized proteins during seed germination have focused on a few individual proteins using genetic techniques 24-26. However, it is difficult to carry out large-scale identification of phosphorylation sites in nuclear proteins during rice seed germination. Phosphorylation plays a key role in protein transport between nucleus and cytoplasm 23. To understand the role of nuclear proteins and their phosphorylation during germination of rice seeds, a gel-free/label-free proteomic technique was used. After extraction and 4 ACS Paragon Plus Environment

Page 4 of 50

Page 5 of 50

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

digestion of nuclear proteins, phosphopeptides were enriched and analyzed using a phosphopeptide enrichment reagent. These nuclear phosphopeptides were analyzed using nano liquid chromatography (LC) coupled tandem mass spectrometry (MS). In addition, to identify the key nuclear phosphoproteins involved in rice seed germination, in silico protein-protein interactions were analyzed using hierarchical clustering and the S-system differential equation. Furthermore, gene expression levels of key nuclear phosphoproteins were analyzed using quantitative reverse-transcription polymerase chain reaction (qRT-PCR) and were compared with the abundance of phosphoproteins.

EXPERIMENTAL PROCEDURES Plant Material, Seed Germination, and Treatment Seeds of rice (Oryza sativa cultivar Nipponbare) were dehusked and placed in culture plates (10 cm x 10 cm) for germination. Rice seeds were germinated on moist mud bed in agriculture and wet filter paper in this experiment 7,27. Seeds were imbibed for 12 h or 24 h in deionized water in a dark growth chamber at 25ºC. Embryos were separated from seeds after 0, 12, or 24 h of imbibition using a knife blade and collected on ice. In total, 50 seeds were prepared for proteomic experiments and 1,400 for nuclear phosphoproteomics experiments. For RNA experiments, 50 seeds were prepared. Rice seeds were treated with 0.01% dimethyl sulfoxide (control), 1% phosphatase inhibitor (Sigma, St. Louis, MO, USA), 10 µM GAs (Wako, Osaka, Japan), and 10 µM ABA (Tokyo Chemical Industry, Tokyo, Japan) solutions, and embryos were collected after 0, 12 and 24 h of treatments. Three independent experiments were performed as biological replicates for all experiments (Supplemental Figure 1).

Isolation of Nuclei and Extraction of Nuclear Proteins 5 ACS Paragon Plus Environment

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nuclei were isolated using a Plant Nuclei Isolation/Extraction Kit (Sigma) following the manufacturer’s protocol with some modifications. Rice embryos (2 g) were ground with 5 mL 1x Nuclei Isolation buffer (Sigma) using a mortar and pestle. The homogenate was filtered through a double layer of Filter Mesh 100 (Sigma) and transferred to a 15 mL tube, then centrifuged at 1,300 x g for 10 min at 4°C. The resulting pellet was resuspended in Nuclei Isolation buffer containing protease inhibitor mixture (Roche, Werk Penzberg, Germany) and phosphatase-inhibitor mixture, then layered on top of a 60% Percoll/2.3 M sucrose cushion. Percoll was diluted to 60% using Nuclei Isolation buffer, and 2.3 M sucrose was supplied in kit. After centrifugation at 3,200 x g for 30 min at 4°C, the middle layer between 60% Percoll and 2.3 M sucrose was carefully collected with a Pasteur pipette. The purified nuclei were washed 2 times using Nuclei Isolation buffer. Nuclei were vortexed at 4°C for 20 min with extraction buffer (Sigma) containing protease inhibitor mixture and phosphatase-inhibitor mixture, and then sonicated in ice water for 20 min. Vortexing and sonication were repeated once. After sonication, the homogenate was centrifuged at 12,000 x g for 30 min at 4°C and the supernatant was collected as nuclear protein (Supplemental Figure 2). Protein concentrations were measured using the Bradford assay with bovine serum albumin as the standard 28.

Extraction of Total Proteins Rice embryos (0.1 g) were ground to powder in liquid nitrogen using a mortar and pestle. The powder was transferred to an acetone solution containing 10% trichloroacetic acid and 0.07% 2-mercaptoethanol, and the resulting mixture was vortexed and then sonicated for 10 min. The suspension was incubated for 1 h at -20°C with vortexing every 15 min and then centrifuged at 9,000 × g at 4°C for 20 min. The supernatant was discarded and the pellet was washed twice with 0.07% 2-mercaptoethanol in acetone. The pellet was dried using a 6 ACS Paragon Plus Environment

Page 6 of 50

Page 7 of 50

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

Speed-Vac concentrator (Savant Instruments, Hickville, NY, USA) and then resuspended in lysis buffer consisting of 7 M urea, 2 M thiourea, 5% CHAPS, and 2 mM tributylphosphine by vortexing for 1 h at 25°C. The suspension was centrifuged at 20,000 × g for 20 min at 25°C and the supernatant was collected as total protein. Protein concentrations were measured using the Bradford assay as above.

Western Blot Analysis Proteins were extracted using SDS sample buffer containing 60 mM Tris-HCl (pH 6.8), 2 % SDS, 10% glycerol, and 5% 2-mercaptoethanol 29. The protein concentration was measured with a Pierce 660 nm Protein Assay Kit with Ionic Detergent Compatibility Reagent (Thermo Fisher Scientific, Rockford, IL, USA). Proteins were separated by electrophoresis on a 12% SDS-polyacrylamide gel and then transferred onto a polyvinylidene difluoride membrane using a semidry transfer blotter. The blotted membrane was incubated overnight at 4°C in blocking buffer consisting of 20 mM Tris-HCl (pH 7.5), 500 mM NaCl, and 5% skim milk (Difco, Sparks, MD, USA). After blocking, the membrane was incubated with a 1:8,000 dilution of anti-histone H3 antibody (Santa Cruz Biotechnology, Santa Cruz, CA, USA) or with a 1:3,000 dilution of anti-ascorbate peroxidase antibody (provided by Proteomic Laboratory in National Institute of Crop Science, Tsukuba, Japan) for 1 h at room temperature. Anti-rabbit IgG conjugated with horseradish peroxidase (Bio-Rad, Hercules, CA, USA) was used as the secondary antibody. After a 1 h incubation with the secondary antibody, signals were detected using Chemi-Lumi One Super Kit (Nacalai Tesque, Kyoto, Japan) following the manufacturer’s protocol, and the signals were visualized using a LAS-3000 luminescent image analyzer (Fujifilm, Tokyo, Japan). Coomassie brilliant blue staining was used as a loading control. The relative intensities of bands were calculated using Quantity One software (version 4.5, Bio-Rad). 7 ACS Paragon Plus Environment

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Purification and Digestion of Proteins Proteins (100 µg) were purified by phase separation in the organic layer. After adjusted the volume to 150 µL (100 µg), 600 µL methanol was added and the resulting solution was thoroughly mixed. Subsequently, 150 µL chloroform was added to the solution, which was again mixed. To induce phase separation, 450 µL water was added to the solution, and the resulting mixture was vortexed and then centrifuged at 20,000 x g for 5 min. The upper aqueous phase was discarded, and 450 µL methanol was added to the organic phase. The samples were centrifuged at 20,000 x g for 5 min, the resulting supernatants were discarded, and the pellets were dried and resuspended in 50 mM NH4HCO3. Proteins were reduced with 50 mM dithiothreitol for 30 min at 56°C, alkylated with 50 mM iodoacetamide for 30 min at 37°C in darkness, and digested with trypsin at a 1:100 enzyme/protein concentration at 37°C for 16 h. The treated peptides were desalted using a MonoSpin C18 column (GL Sciences, Tokyo, Japan).

Phosphopeptide Enrichment Polymer-based Metal-ion Affinity Capture (PolyMAC) Phosphopeptide Enrichment Reagent (Tymora, West Lafayette, IN, USA) was used to enrich phosphopeptides. Desalted peptides (200 µg) were resuspended in 100 µL loading buffer consisting of 100 mM glycolic acid, 1% trifluoroacetic acid, and 50% acetonitrile, to which 10 µL PolyMAC-Ti reagent was added. After shaking the resulting mixture for 30 min, 200 µL 300 mM HEPES-NaOH (pH 7.7) was added to increase the pH of the sample above 6.3. After brief vortexing, 50 µL of magnetic capture beads were added to the sample, which was shaken vigorously for 10 min. The beads were separated from the solution using a magnetic separator rack (Invitrogen Dynal AS, Oslo, Norway). Captured beads were washed once in loading buffer and twice in washing 8 ACS Paragon Plus Environment

Page 8 of 50

Page 9 of 50

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

buffer, consisting of 100 mM acetic acid, 1% trifluoroacetic acid, and 80% acetonitrile. Phosphopeptides were then eluted by washing the beads twice with 100 µL 400 mM NH4HCO3. After drying the eluent, the pellet was resuspended well in 10 µL 0.25% formic acid. The phosphopeptide solution was collected after centrifugation.

Mass Spectrometric Analysis Peptides in 0.1% formic acid were loaded onto an Ultimate 3,000 nano-LC system (Dionex, Germering, Germany) equipped with a C18 PepMap trap column (300 µm ID × 5 mm, Dionex). The peptides were eluted from the trap column and separated using 0.1% formic acid in acetonitrile at a flow rate of 200 nL/min on a C18 Tip column (75 µm 1D × 120 mm, Nikkyo Technos, Tokyo, Japan) with a spray voltage of 1.8 kV. The peptide ions in the spray were analyzed on a nanospray LTQ XL Orbitrap mass spectrometer (Thermo Fisher Scientific, San Jose, CA, USA) operated in data-dependent acquisition mode with the installed Xcalibur software (version 2.0.7, Thermo Fisher Scientific). Full-scan mass spectra were acquired in the Orbitrap mass spectrometer over 400-1,500 m/z with a resolution of 30,000. A lock mass function was used to obtain high mass accuracy 30. As the lock mass, the ions C24H39O4+ (m/z 391.28429), C14H46NO7Si7+ (m/z 536.16536), and C16H52NO8Si8+ (m/z 610.18416) were used. Values for ion isolation window were set as follows: activation type was collision-induced dissociation, default charge state was 2, isolation width was 2.0 m/z, normalized collision energy was 35.0, and activation time was 30.000 ms. Values for dynamic exclusion were determined as follows: repeat count was 2, repeat duration was 30 sec, exclusion list size was 500, exclusion duration was 90 sec, and exclusion mass width was ± 1.5 Da. The three most intense precursor ions above a threshold value of 500 were selected for collision-induced fragmentation. The acquired mass spectra were used for protein identification. 9 ACS Paragon Plus Environment

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Protein Identification Using Acquired Mass Spectrometry Data Protein identification was performed using the Mascot search engine (version 2.5.1, Matrix Science, London, UK) and a rice protein database (50,253 sequences and 15,266,515 residues) obtained from The Rice Annotation Project Database (RAP-DB, http://rapdb.dna.affrc.go.jp), including protein sequences supported by FL-cDNA and EST data (IRGSP-1.0_protein_2013-4-24) and protein sequence predicted computationally (IRGSP-1.0_predicted protein_2013-3-9) 31,32. DTA files were generated from acquired raw data files and converted to Mascot generic files using Proteome Discoverer software (version 1.4.0.288, Thermo Fisher Scientific). The parameters used in the Mascot searches were: carbamidomethylation of cysteine was set as a fixed modification, and oxidation of methionine was set as a variable modification. For peptides after phosphopeptide enrichment, carbamidomethylation of cysteine was set as a fixed modification, and phosphoST, phosphoY, and oxidation of methionine were set as variable modifications. For both normal peptides and peptides after phosphopeptide enrichment, trypsin was specified as the proteolytic enzyme, and one missed cleavage was allowed. Peptide mass tolerance was set at 10 ppm, fragment mass tolerance was set at 0.8 Da, and peptide charges were set at +2, +3, and +4. An automatic decoy database search was performed as part of the search. Mascot results were filtered with the Mascot percolator to improve the accuracy and sensitivity of peptide identification 33. False discovery rates for peptide identification of all searches were less than 1.0%. Peptides with a percolator ion score of more than 13 (p0.9 were considered candidate interactions. In the model protein interaction diagram, a red arrow indicates an inductive interaction and corresponds to gij > 0 in the S-system differential equation, and a blue T-bar indicates a suppressive interaction and corresponds to gij < 0 in the S-system differential equation 39.

RNA Extraction and Quantitative Reverse Transcription Polymerase Chain Reaction Analysis Rice embryos (0.1 g) were ground into fine powder in liquid nitrogen with a sterilized mortar and pestle. Total RNA was extracted from the powder using an RNeasy Plant Mini kit (Qiagen, Valencia, CA, USA). RNA was reverse-transcribed using an iScript cDNA synthesis kit (Bio-Rad) according to the manufacturer’s instructions. qRT-PCR was performed in a 10 µL reaction using SsoAdvanced SYBR Green Supermix (Bio-Rad) and a MyiQ single-color real-time PCR detection system (Bio-Rad). The PCR conditions were as follows: 95°C for 30 sec, followed by 45 cycles of 95°C for 10 sec and 60°C for 30 sec. Gene expression was normalized using 18S rRNA as an internal control. The primers were designed using the Primer3 web interface (http://frodo.wi.mit.edu) (Supplemental Table 1). Specificity of primers 12 ACS Paragon Plus Environment

Page 12 of 50

Page 13 of 50

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

was checked by BLASTN search against the RAP-DB with the designed primers as queries, by melt curve analysis and by agarose gel electrophoresis of amplified fragments.

Statistical Analysis The statistical significance of results was evaluated using Student’s t-test when only two groups were compared. The statistical significance of results from multiple groups was evaluated using a one-way ANOVA test. Both calculations were performed using GraphPad software (version 5.0, GraphPad Software, San Diego, CA, USA). A p value < 0.05 was considered to be statistically significant.

RESULTS Purity of Nuclear Proteins Isolated from Rice Embryos To understand the role of nuclear proteins and their phosphorylation during rice seed germination, nuclear proteins were isolated. Nuclei were isolated from rice embryos collected after 24 h of imbibition (Supplemental Figure 2). In total, 100 µg of nuclear proteins was extracted from 1 g of rice embryos. Western blotting with an anti-histone H3 and ascorbate peroxidase antibodies was performed to check the purity of nuclear proteins. The antiascorbate peroxidase antibody cross-reacted with 31-kDa proteins and the anti-histone H3 antibody cross-reacted with 17-kDa nuclear proteins (Figure 1). The relative intensity of cytosolic ascorbate peroxidase was markedly stronger in the total protein fraction than in the nuclear protein fraction. The strong intensity of bands corresponding to histone H3 indicated that the nuclei were significantly enriched. Total protein and nuclear proteins extracted from rice embryos were analyzed using nanoLC-MS/MS. A total of 349 total proteins (Supplemental Table 2) and 362 nuclear proteins (Supplemental Table 3) were identified. Comparing total proteins with nuclear 13 ACS Paragon Plus Environment

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

proteins, 149 proteins were in common, 200 were specific to total proteins, and 213 were apparently specific to nuclear proteins (Figure 2A). In order to evaluate the purity of nuclear protein, the localization of each protein from the total protein and nuclear protein fractions was predicted using intracellular targeting prediction programs (Figure 2B); 71 (20%) of the proteins in the total protein and 196 (54%) of the nuclear proteins were predicted to be located in nuclei (Figure 2B). Proteins contaminating the nuclear fraction mainly belonged to the cytoplasm (16%) and chloroplast (18%) (Figure 2B). These results indicated that nuclear proteins in the nucleus fraction were enriched more than 2.5 times compared to total protein by the procedure.

Function of Nuclear Proteins in Rice Embryos after 24 h of Imbibition In order to understand functions of 362 nuclear proteins, they were categorized using MapMan bin codes. In total, 362 proteins were classified into 21 functional categories (Figure 3). Categories related to protein (39%), RNA (13%), development (6%), and DNA (6%) represented the largest number of proteins. The functional category of protein contained 141 proteins and this category was further classified into 4 groups, which were protein synthesis (83%), protein degradation (9%), protein targeting (6%), and others (2%) (Figure 3). Moreover, 47 proteins and 22 proteins respectively belonged to the RNA and DNA categories, with 16 nucleotide-, RNA-, and DNA-binding proteins belonging to the RNA category, and 8 histones belonging to the DNA category (Supplemental Table 3).

Phosphoproteomic Analysis of Nuclear Proteins in Rice Embryos after 24 h of Imbibition To understand the phosphorylation status of nuclear proteins, a phosphoproteomics study was performed. Nuclear proteins were extracted from rice embryos after 24 h of imbibition. Nuclear phosphopeptides were enriched and analyzed by nanoLC-MS/MS. 14 ACS Paragon Plus Environment

Page 14 of 50

Page 15 of 50

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

Based on this technique, 3,467 phosphopeptides (Supplemental Table 4) belonging to 102 phosphorylated proteins (Supplemental Table 5) were identified. Overall, 65% and 30% of phosphorylated peptides respectively contained a single and double phosphoryl group (Figure 4A). Serine, threonine, and tyrosine phosphorylated peptides accounted for 64, 31, and 5%, respectively (Figure 4B). Nuclear phosphoproteins were classified into 14 functional categories (Figure 4C). Among these categories, protein, RNA, DNA, and cell contained the largest number of proteins, comprising 25, 21, 10, and 5%, respectively. The functional category of protein contained 25 proteins, and this category was further classified into 4 groups, which were protein synthesis (64%), protein degradation (24%), protein targeting (8%), and protein folding (4%) (Figure 4C).

Changes in Abundance of Nuclear Phosphoproteins in Rice Embryos during Early Stage of Seed Germination Phosphorylation is a reversible posttranslational modification involved mainly in signal transduction and protein localization. To understand the role of nuclear phosphoproteins during rice seed germination after 0, 12, and 24 h, gel-free/label-free phosphoproteomic analysis was performed. In total, 115 nuclear phosphoproteins were identified during rice seed germination (Supplemental Table 6). After quantitative normalization of nuclear phosphoprotein abundance and one-way ANOVA, 29 nuclear phosphoproteins were identified as significantly changed in abundance during rice seed germination (Table 1, Supplemental Table 7). Based on the changes in abundance of nuclear phosphoproteins at different time points during germination, 29 nuclear phosphoproteins were divided into 5 clusters (Figure 5). Cluster 1 comprised 4 nuclear phosphoproteins that exhibited increased phosphorylation 15 ACS Paragon Plus Environment

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

during germination. Cluster 2 contained 11 nuclear phosphoproteins that reached a peak phosphorylation level at 12 h of imbibition and then decreased. Clusters 3, 4, and 5 consisted of 3, 10, and 1 nuclear phosphoproteins, respectively, that had a sustained decrease in phosphorylation level. DNA and nucleotide binding proteins, zinc finger proteins, and WRC domain containing proteins, which belonged to Cluster 2, reached a peak phosphorylation level at 12 h of imbibition and then decreased (Figure 5, Table 1).

In Silico Protein-Protein Interaction Analysis of Nuclear Phosphoproteins in Rice Embryos during Early Stage of Seed Germination To determine which of the 29 nuclear phosphoproteins that significantly changed in abundance may interact during rice seed germination, time-dependent changes in the abundance of each phosphoprotein were used in conjunction with an S-system based modeling approach. After protein-protein interaction analysis, the network of these nuclear phosphoproteins was revealed (Figure 6). A total of 29 nuclear phosphoproteins were predicted to interact with other proteins. Nuclear phosphoproteins in Cluster 2 mainly induced expression of peptidase S1C_HrtA/DegP2/Q/S protein (protein No. 1) and AT-hook protein 1 (protein No. 2). Sas10/Utp3 protein (protein No. 23) was positioned in the center of network of nuclear phosphoproteins. The abundance of Sas10/Utp3 protein was suppressed by peptidase S1C HrtA/DegP2/Q/S protein and AT-hook protein 1 during rice seed germination.

RNA Expression Level of Nuclear Phosphoproteins in Rice Embryos during Early Stage of Seed Germination Among 29 nuclear phosphoproteins which exhibited significant changes in abundance during rice seed germination, 8 of them were selected for qRT-PCR to check their RNA 16 ACS Paragon Plus Environment

Page 16 of 50

Page 17 of 50

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

expression levels in rice embryos during early stage of seed germination. These 8 nuclear phosphoproteins were both predicted to locate in nuclei using WoLF PSORT and NucPred softwares, and the ratio of protein abundance during the germination process was above or below 10 times. Eight qRT-PCR products were analyzed by agarose gel electrophoresis and the sizes of products were consistent with design (Supplemental Figure 3, numbers 2, 3, 4, 5, 6, 7, 8, and 9). Among 8 examined genes, the mRNA expression levels of 5 genes including nucleotide binding protein, WRC domain containing protein, multidomain cyclophilin type peptidyl prolyl cis-trans isomerase, GTP binding signal recognition particle SRP54 G-domain containing protein, and Sas10/Utp3 protein did not change in rice embryos during early stage of seed germination. However, the mRNA expression level of AT hook DNA binding protein was up-regulated in rice embryos at 12 h and down-regulated at 24 h of imbibition. The mRNA expression level of zinc finger CCCH type protein was down-regulated in rice embryos during early stage of seed germination. In contrast, the mRNA expression level of ribosomal RNA processing protein 7 was up-regulated in rice embryos during early stage of seed germination (Figure 7).

Analysis of the Role of Zinc Finger Proteins and Sas10/Utp3 Protein in Rice Embryos during Early Stage of Seed Germination In 29 nuclear phosphorylated proteins predicted interaction network, two zinc finger proteins which are zinc finger CCCH type and zinc finger BED type located on the upstream of network to interact with other proteins (Figure 6). In Arabidopsis, some zinc finger proteins were confirmed that they involved in seed germination through the interaction between zinc finger proteins and phytohormones such as GA and ABA 40,41. To investigate the role of phosphorylation in rice embryos during early stage of seed germination, the effect of phosphatase inhibitor on seed germination was examined. 17 ACS Paragon Plus Environment

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Furthermore, the mRNA expression levels of two zinc finger proteins were analyzed in rice embryos during early stage of seed germination with phytohormones and phosphatase-inhibitor treatments. Rice seeds were treated with phosphatase inhibitor, GAs, and ABA for 24, 48, and 72 h of imbibition. The phenotypes of seed germination were similar between phosphatase inhibitor and ABA treatments (Figure 8A). The germination rate of seeds treated with phosphatase inhibitor was significantly slow down (Figure 8A). The water uptake of seeds treated with phosphatase inhibitor or ABA was significantly lower than that of control (Figure 8B). The mRNA expression level of zinc finger CCCH type protein was down-regulated in rice embryos during early stage of seed germination with phosphatase-inhibitor treatment (Figure 9A). In contrast, the mRNA expression level of zinc finger BED type protein was up-regulated in rice embryos during early stage of seed germination with phosphatase inhibitor and ABA treatments (Figure 9B). Additionally, the mRNA expression of Sas10/Utp3 protein was analyzed (Figure 10). It was no change in rice embryos during early stage of seed germination with GAs and ABA treatments. However, it was up-regulated with phosphatase inhibitor in rice embryos during early stage of seed germination (Figure 10).

DISCUSSION Seed germination, which usually starts with the uptake of water by dry seed and terminates with elongation of radicle 2,3. Rice seeds were germinated on moist mud bed in agriculture and wet filter paper in experiment 7,27. In this study, some hypoxic inducible proteins were detected by proteomic analysis for total proteins (Supplemental Table 2). In both germination conditions, some part of seeds but not the whole seeds were submerged in water. It indicated that rice can germinate in low oxygen condition 42. Howell et al. 15,16 18 ACS Paragon Plus Environment

Page 18 of 50

Page 19 of 50

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

characterized changes in water content and metabolic activity in embryo during rice seed germination under aerobic and anaerobic conditions, and furthermore Narsai et al. 43 defined 10 metabolites and 1136 transcripts as aerobic responders, and 13 metabolites and 730 transcripts as anaerobic responders in rice embryos and young seedlings. In order to assess whether seeds germinate on wet filter paper are under anaerobic condition, the total proteins identified in this nuclear phosphoproteomics were compared with transcriptomic data in research by Narsai et al. 43. In 349 total proteins, there are 8 and 4 genes which encoding proteins were identified in anaerobic and aerobic responders respectively (Supplemental Table 8). It might indicate that there is hypoxic stress in rice embryos during early stage of seed germination to some extent. However, it might because at the early stage of germination, anaerobic respiration is the major source of ATPs 44,45. Seed germination is a well-regulated process including a series of signal transduction and gene expression steps 46. The first 24 h of rice seed imbibition is usually regarded as the most important stage because germination-related metabolic reactions are reactivated during this period 3. Transcriptomic analysis of germinating rice seeds indicated that the onset of mRNA biosynthesis starts at an early stage of germination 4. To better understand the phosphorylation of nuclear proteins in rice embryos at this stage, the phosphorylation status of nuclear proteins was analyzed. In this study, 22 transcription factors were identified as nuclear phosphoproteins expressed in rice embryos after 24 h of imbibition (Supplemental Table 5). Phosphoproteomics of total protein indicated that a number of transcription factors change phosphorylation status in rice embryos during the early stage of seed germination 17. However, most of the transcription factors identified as phosphoproteins in nucleus (supplemental table 6) in this study didn’t exist in previous study 17. It indicated that subcellular proteomic analysis of the nucleus during rice seed germination is useful. Of the transcription factors, two AT-hook DNA binding proteins, which were AT-hook 19 ACS Paragon Plus Environment

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

DNA binding protein 1 (Os08t0512400-02) and AT-hook DNA binding conserved site domain containing protein (Os04t0501600-01), were identified in the nucleus and demonstrated significantly increased phosphorylation levels after 12 h of imbibition, which then dramatically decreased (Table 1, Figure 5). AT-hook DNA binding proteins are AT-rich DNA binding proteins belonging to the high mobility group non-histone chromosomal protein HMG-I(Y) family 47. In plants, an AT-hook motif was found to bind to special regions of chromosomal DNA called nuclear matrix attachment regions 48. In Arabidopsis, an AT-hook protein called AGF1 responds to the GA-controlled negative feedback of AtGA3ox1, encoding GA3-oxidase, at the posttranscriptional or posttranslational level 49. In this study, two AT-hook DNA binding proteins had phosphorylation sites (Supplemental Table 4) and the abundance of these proteins were increased after 12 h of imbibition and then dramatically decreased. These results suggest that they might be involved in the homeostasis of GAs through changing protein phosphorylation status during rice seed germination. Based on analysis of in silico protein-protein interactions of 29 nuclear phosphoproteins, Sas10/Utp3 was interacted with many proteins and it might be positioned in the center of the network of nuclear phosphoproteins. Furthermore, its abundance was suppressed by peptidase S1C HrtA/DegP2/Q/S (Os06t0234100) protein and AT-hook protein 1 during germination (Figure 6). Peptidase S1C HrtA/DegP2/Q/S, which contains a recognizable PDZ domain, functions in mRNA export from the nucleus and protein import into the nucleus 50,51. In Arabidopsis, rice, and poplar, the HrtA/DegP2 family contains 16, 15, and 17 homolog proteins, respectively 52-54. However, only one DegP protein (Deg9, At5g40200) has been identified as a nuclear protein experimentally and by a software algorithm in Arabidopsis 55-57, and it has two orthologs (Os02t0742500 and Os06t0234100) in rice 51. It has been hypothesized that this protein is involved in ribosomal DNA 20 ACS Paragon Plus Environment

Page 20 of 50

Page 21 of 50

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

transcription and ribosomal RNA processing/modification 50. In this study, peptidase S1C HrtA/DegP2/Q/S protein demonstrated an increased phosphorylation level during germination, which indicates that DegP protein in the nucleus of rice embryos might be involved in seed germination. Sas10/Utp3 protein contains a Sas10/Utp3 domain and is a component of the U3 ribonucleoprotein complex. Furthermore, it is a regulator of chromatin silencing in yeast 58 and is involved in 18S ribosomal RNA biogenesis 59. CANu1 protein, which contains a Sas10/Utp3 domain, was observed to be localized in nuclei in humans 60. This is consistent with our results using two types of software that predicted localization of Sas10/Utp3 protein in nuclei in rice embryos during germination (Table 1). However, there are not many studies of Sas10/Utp3 protein in plants. This study identified Sas10/Utp3 protein as positioned in the center of the phosphorylation network in nuclei of rice embryos (Figure 5), and the mRNA expression level of Sas10/Utp3 didn’t change in embryos during early stage of seed germination. Furthermore, the mRNA expression level of Sas10/Utp3 did not change with phytohormones, but it was up-regulated with phosphatase inhibitor in rice embryos during early stage of seed germination (Figure 10). It is suggested that Sas10/Utp3 protein might play through its phosphorylation and dephosphorylation switching during the early stage of seed germination without the relation to phytohormone. In this study, two zinc finger proteins exhibited high phosphorylation levels after 12 h of imbibition and then decreased (Table 1, Figure 5). Zinc finger proteins containing conserved cysteine and histidine ligands were first recognized in Xenopus transcription factor IIIA 61. These proteins have various functions ranging from DNA binding and RNA binding to protein-protein interactions and membrane association 62. OsLOL1, a C2C2 type zinc finger protein, promotes seed germination through increasing the GA level in rice. OsKO2, an important biosynthetic enzyme in GA biosynthesis, is upregulated through 21 ACS Paragon Plus Environment

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

OsLOL1, which interacts with OsbZIP58 during germination 41. AtTZF4/5/6, CCCH zinc finger proteins, are reportedly involved in seed germination in Arabidopsis. AtTZF4/5/6 are up-regulated by ABA and downregulated by GA, and their proteins also act as a positive regulator for ABA- and negative regulator for GA-mediated seed germination responses 40. In this study, the mRNA expression level of zinc finger CCCH type protein (Os02t0161200-01) did not change in rice embryos during early stage of seed germination with phosphatase inhibitor and ABA treatments (Figure 9A). Zinc finger BED domain containing protein were thought to function as either transcription activators or repressors by modifying local chromatin structure on binding to GC-rich sequences 63. It also was reported that zinc BED type proteins involve in potato and barley disease resistance 64,65. However, there is no report about zinc finger BED type protein in seed during germination. In this study, the mRNA expression level of zinc finger BED type protein (Os03t0733400-01) was up-regulated in rice embryos during early stage of seed germination with phosphatase inhibitor and ABA treatments (Figure 9B). It might be because seed needs more dephosphorylating status of zinc finger BED type proteins to promote germination. These results indicate that zinc finger BED type and CCCH type proteins might mediate rice seed germination through phosphorylation and posttranslational modifications.

CONCLUSIONS The early stage of seed germination is the first step in rice life cycle without visible morphological change 3. To investigate the mechanism controlling the physiological status shifting in the early stage of rice seed germination, nuclear phosphoproteomics was performed. The main findings of this nuclear phosphoproteomics are as follows: (i) Out of 115 identified nuclear phosphoproteins, 29 nuclear phosphoproteins significantly changed 22 ACS Paragon Plus Environment

Page 22 of 50

Page 23 of 50

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

the phosphorylation level during early stage of seed germination; (ii) Study of in silico protein-protein interactions for 29 nuclear phosphoproteins indicated that Sas10/Utp3 protein was interacted with most phosphoproteins in nuclei, and its mRNA expression level didn’t change in rice embryos during early stage of seed germination; (iii) Zinc finger proteins which are zinc finger CCCH type protein and zinc finger BED type protein play important roles through their phosphorylation/dephosphorylation during early stage of seed germination. To investigate the role of phosphorylated status of zinc finger proteins and Sas10/Utp3 proteins in rice embryos during early stage of seed germination, the effect of phosphatase inhibitor on seed germination was examined. Physiological experiments indicated that protein dephosphorylation might promote rice seed germination. Taken together, it is suggested that phosphorylation and dephosphorylation of nuclear proteins involve in rice seed germination, and transcription factors such as zinc finger CCCH type protein and zinc finger BED type protein might have key role through phosphorylation/dephosphorylation. Furthermore, Sas10/Utp3 protein might be interacted with most nuclear phosphoproteins in rice embryos to mediate the early stage of seed germination.

ASSOCIATED CONTENT Supporting Information Supplemental Table 1. Primer sequences of genes selected for qRT-PCR. Supplemental Table 2. Total proteins identified in rice embryos after 24 h of imbibition. Supplemental Table 3. Nuclear proteins identified in rice embryos after 24 h of imbibition. Supplemental Table 4. The probabilities of phosphorylation site assignment for nuclear phosphopeptides in rice embryos after 24 h of imbibition. Supplemental Table 5. Nuclear phosphoproteins identified in rice embryos after 24 h of 23 ACS Paragon Plus Environment

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

imbibition. Supplemental Table 6. Nuclear phosphoproteins identified in rice embryos during early stage of germination. Supplemental Table 7. The probabilities of phosphorylation site assignment for nuclear phosphopeptides in rice embryos during early stage of germination. Supplemental Table 8. Comparison of the proteome data in this study with transcriptome data in Narsai et al. 43 in rice embryo. Supplemental Figure 1. Experimental design of nuclear phosphoproteomic study. Supplemental Figure 2. Procedure for isolation of nuclei. Supplemental Figure 3. Agarose gel electrophoresis of qRT-PCR products. This material is available free of charge via the Internet at http://pubs.acs.org.

Accession Codes The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium (http://proteomecentral.proteomexchange.org) via the PRIDE partner repository 66

with the dataset identifier PXD001774.

AUTHOR INFORMATION Corresponding Author * (S.K.) Tel: +81-29-838-8693. Fax: +81-29-838-8694. E-mail: [email protected]. Author Contribution # M.L. and X.Y. contributed equally to this work.

Notes The authors declare no competing financial interest. 24 ACS Paragon Plus Environment

Page 24 of 50

Page 25 of 50

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

ACKNOWLEDGMENTS The authors thank Dr. K. Sugimoto at the National Institute of Agrobiological Sciences for providing rice seeds. The authors thank Dr. Y. Nanjo and Dr. K. Nishizawa at the National Institute of Crop Science for useful suggestion and comments. The authors also thank Ms. C. Kawashima and Ms. A. Sato at the National Institute of Crop Science for experimental support during this research. This work was supported by the National Natural Science Foundation of China (31271805) and by the National Agriculture and Food Research Organization.

REFERENCES (1) Bazzaz, F. A.; Ackerly, D. D.; Reekie, E. G. Reproductive Allocation in Plants. In Seeds: The Ecology of Regeneration in Plant Communities, 2nd Edition; Fenner, M., Eds.; CABI Publishing: New York, NY, USA, 2000; pp 1-22. (2) Bewley, J.D.; Black, M. Seeds: Physiology of Development and Germination; Plenum Press: New York, NY, USA, 1994; pp 377-420. (3) Bewley, J. D. Seed germination and dormancy. Plant Cell 1997, 9, 1055-1066. (4) Howell, K. A.; Narsai, R.; Carroll, A.; Ivanova, A.; Lohse, M.; Usadel, B.; Harvey Millar A. Whelan, J. Mapping metabolic and transcript temporal switches during germination in rice highlights specific transcription factors and the role of RNA instability in the germination process. Plant Physiol. 2009, 149, 961-980. (5) Kim, S. T.; Kang, S. Y.; Wang, Y.; Kim, S. G.; Hwang, D. H.; Kang, K. Y. Analysis of embryonic proteome modulation by GA and ABA from germinating rice seeds. Proteomics 2008, 8, 3577-3587. (6) Liu, S. J.; Xu, H. H.; Wang, W. Q.; Li, N.; Wang, W. P.; Møller, I. M.; Song, S. Q. A 25 ACS Paragon Plus Environment

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

proteomic analysis of rice seed germination as affected by high temperature and ABA treatment. Physiol. Plant. 2014, doi:10.1111/ppl.12292 (7) He, D. L.; Han, C.; Yao, J. L.; Shen, S. H.; Yang, P. F. Constructing the metabolic and regulatory pathways in germinating rice seeds through proteomic approach. Proteomics 2011, 11, 2693-2713. (8) Han, C.; He, D. L.; Li, M.; Yang, P. F. In-depth proteomic analysis of rice embryo reveals its important roles in seed germination. Plant Cell Physiol. 2014a, 55, 1826-1847. (9) Perata, P.; Pozueta-Romero, J.; Akazawa, T.; Yamaguchi, J. Effect of anoxia on starch breakdown in rice and wheat seeds. Planta 1992, 188, 611-618. (10) Sadiq, I.; Fanucchi, F.; Paparelli, E.; Alpi, E.; Bachi, A.; Alpi, A.; Perata, P. Proteomic identification of differentially expressed proteins in the anoxic rice coleoptile. J. Plant Physiol. 2011, 168, 2234-2243. (11) Shingaki-Wells, R. N.; Huang, S. B.; Taylor, N. L.; Carroll, A. J.; Zhou, W. X.; Millar, A. H. Differential molecular responses of rice and wheat coleoptiles to anoxia reveal novel metabolic adaptations in amino acid metabolism for tissue tolerance. Plant Physiol. 2011, 156, 1706-1724. (12) Yano, H.; Wong, J. H.; Lee, Y. M.; Cho, M. J.; Buchanan, B. B. A strategy for the identification of proteins targeted by thioredoxin. Proc. Natl. Acad. Sci. USA 2001, 98, 4794-4799. (13) Yano, H.; Kuroda, M. Disulfide proteome yields a detailed understanding of redox regulations: A model study of thioredoxin‐linked reactions in seed germination. Proteomics 2006, 6, 294-300. (14) Ferrari, F.; Fumagalli, M.; Profumo, A.; Viglio, S.; Sala, A.; Dolcini, L.; Temporini, C. Nicolis, S.; Merli, D.; Corana, F.; Casado, B.; Iadarola, P. Deciphering the proteomic 26 ACS Paragon Plus Environment

Page 26 of 50

Page 27 of 50

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

profile of rice (Oryza sativa) bran: a pilot study. Electrophoresis 2009, 30, 4083-4094. (15) Howell, K. A.; Millar, A. H.; Whelan, J. Ordered assembly of mitochondria during rice germination begins with promitochondrial structures rich in components of the protein import apparatus. Plant Mol. Biol. 2006, 60, 201-223. (16) Howell, K. A.; Cheng, K.; Murcha, M. W.; Jenkin, L. E.; Millar, A. H.; Whelan, J. Oxygen initiation of respiration and mitochondrial biogenesis in rice. J. Biol. Chem. 2007, 282, 15619-15631. (17) Han, C.; Yang, P. F.; Sakata, K.; Komatsu, S. Quantitative proteomics reveals the role of protein phosphorylation in rice embryos during early stages of germination. J. Proteome Res. 2014b, 13, 1766-1782. (18) Han, C.; Wang, K.; Yang, P. F. Gel-based comparative phosphoproteomic analysis on rice embryo during germination. Plant Cell Physiol. 2014c, 55, 1376-1394. (19) Henfrey, R. D.; Slater, R. J. Isolation of Plant Nuclei. In New Nucleic Acid Techniques; Walker, J. M., Eds.; Humana Press: New York, NY, USA, 1988; pp 447-452. (20) Steidl, S.; Tüncher, A.; Goda, H.; Guder, C.; Papadopoulou, N.; Kobayashi, T.; Tsukagoshi, N.; Kato, M. Brakhage, A. A. A single subunit of a heterotrimeric CCAAT-binding complex carries a nuclear localization signal: piggy back transport of the pre-assembled complex to the nucleus. J. Mol. Biol. 2004, 342, 515-524. (21) Jones, A. M.; MacLean, D.; Studholme, D. J.; Serna-Sanz, A.; Andreasson, E.; Rathjen, J. P.; Peck, S. C. Phosphoproteomic analysis of nuclei-enriched fractions from Arabidopsis thaliana. J. Proteomics 2009, 72, 439-451. (22) Deribe, Y. L.; Pawson, T.; Dikic, I. Post-translational modifications in signal integration. Nat. Struct. Biol. 2010, 17, 666-672. (23) Kaffman, A.; O'Shea, E. K. Regulation of nuclear localization: a key to a door. Annu. Rev. Cell Dev. Biol. 1999, 15, 291-339. 27 ACS Paragon Plus Environment

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(24) Chiatante, D.; Onell, E.; Patrignani, G.; Scippa, G. S. Localization of a nuclear protein (QP47) in embryonic meristems during seed maturation and germination andits distribution among crop plants. J. Exp. Bot. 1995, 46, 815-821. (25) Huq, E.; Tepperman, J. M.; Quail, P. H. GIGANTEA is a nuclear protein involved in phytochrome signaling in Arabidopsis. Proc. Natl. Acad. Sci. USA 2000, 97, 9789-9794. (26) Cho, J. N.; Ryu, J. Y.; Jeong, Y. M.; Park, J.; Song, J. J.; Amasino, R. M.; Noh, B.; Noh, Y. S. Control of seed germination by light-induced histone arginine demethylation activity. Dev. Cell 2012, 22, 736-748. (27) Yamauchi, M.; Aragones, D. V.; Casayuran, P. R.; Cruz, P. C.; Asis, C. A.; Cruz, R. T. Seedling establishment and grain yield of tropical rice sown in puddled soil. Agron. J. 2000, 92, 275-282. (28) Bradford, M. M. A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal. Biochem. 1976, 72, 248-254. (29) Laemmli, U. K. Cleavage of structural proteins during the assembly of the head of bacteriophage T4. Nature 1970, 227, 680-685. (30) Olsen, J.V.; de Godoy, L. M.; Li, G.; Macek, B.; Mortensen, P.; Pesch, R.; Makarov, A.; Lange, O.; Horning, S.; Mann, M. Parts per million mass accuracy on an orbitrap mass spectrometer via lock mass injection into a C-trap. Mol. Cell. Proteomics 2005, 4, 2010-2021. (31) Kawahara, Y.; Bastide, M.; Hamilton, J. P.; Kanamori, H.; McCombie, R.; Ouyang, S.; Schwartz, D. C.; Tanaka, T.; Wu, J. Z.; Zhou, S. G.; Childs, K. L.; Davidson, R. M.; Lin, H. N.; Quesada-Ocampo, L.; Vaillancourt, B.; Sakai, H.; Lee, S. S.; Kim, J.; Numa, H.; Itoh, T.; Buell, C. R.; Matsumoto, T. Improvement of the Oryza sativa Nipponbare reference genome using next generation sequence and optical map data. Rice 2013, 6, 4. 28 ACS Paragon Plus Environment

Page 28 of 50

Page 29 of 50

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

(32) Sakai, H.; Lee, S. S.; Tanaka, T.; Numa, H.; Kim, J.; Kawahara, Y.; Wakimoto, H.; Yang, C. C.; Iwamoto, M.; Abe, T.; Yamada, Y.; Muto, A.; Inokuchi, H.; Ikemura, T.; Matsumoto, T.; Sasaki, T.; Itoh, T. Rice Annotation Project Database (RAP-DB): an integrative and interactive database for rice genomics. Plant Cell Physiol. 2013, 54, e6. (33) Brosch, M.; Yu, L.; Hubbard, T.; Choudhary, J. Accurate and sensitive peptide identification with Mascot percolator. J. Proteome Res. 2009, 8, 3176-3181. (34)Usadel, B.; Nagel, A.; Thimm, O.; Redestig, H.; Blaesing, O.E.; Palacios-Rofas, N.; Selbig, J.; Hannemann, J.; Piques, M. C.; Steinhauser, D.; Scheible, W. R.; Gibon, Y.; Morcuende, R.; Weicht, D.; Meyer, S.; Stitt, M. Extension of the visualization tool MapMan to allow statistical analysis of arrays, display of corresponding genes, and comparison with known responses. Plant Physiol. 2005, 138, 1195-1204. (35) Horton, P.; Park, K. J.; Obayashi, T.; Fujita, N.; Harada, H.; Adams-Collier, C. J.; Nakai, K. WoLF PSORT: protein localization predictor. Nucleic Acids Res. 2007, 35, W585-W587. (36) Brameier, M.; Krings, A.; MacCallum, R. M. NucPred-predicting nuclear localization of proteins. Bioinformatics 2007, 23, 1159-1160. (37) Sturn, A.; Quackenbush, J.; Trajanoski, Z. Genesis: cluster analysis of microarray data. Bioinformatics 2002, 18, 207-208. (38) Voit, E. O. Models of biochemical systems. Computational analysis of biochemical systems: a practical guide for biochemists and molecular biologists. Cambridge University Press, Cambridge, UK, 2000; pp 37-75. (39) Tanaka, N.; Mitsui, S.; Nobori, H.; Yanagi, K.; Komatsu, S. Expression and function of proteins during development of the basal region in rice seedlings. Mol. Cell. Proteomics 2005, 4, 796-808. (40) Bogamuwa, S.; Jang, J. C. The Arabidopsis tandem CCCH zinc finger proteins AtTZF4, 29 ACS Paragon Plus Environment

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

5 and 6 are involved in light, abscisic acid and gibberellic acid mediated regulation of seed germination. Plant, cell & environ. 2013, 36, 1507-1519. (41) Wu, J. H.; Zhu, C. F.; Pang, J. H.; Zhang, X. R.; Yang, C. L.; Xia, G. X.; Tian, Y. C.; He, C. Z. OsLOL1, a C2C2-type zinc finger protein, interacts with OsbZIP58 to promote seed germination through the modulation of gibberellin biosynthesis in Oryza sativa. Plant J. 2014, 80, 1118-1130. (42) Magneschi, L.; Perata, P. Rice germination and seedling growth in the absence of oxygen. Ann. Bot., 2009, 103, 181-196. (43) Narsai, R.; Howell, K. A.; Carroll, A.; Ivanova, A.; Millar, A. H.; Whelan, J. Defining core metabolic and transcriptomic responses to oxygen availability in rice embryos and young seedlings. Plant Physiol, 2009, 151, 306-322. (44) Raymond, P.; Al-Ani, A.; Pradet, A. ATP production by respiration and fermentation, and energy charge during aerobiosis and anaerobiosis in twelve fatty and starchy germinating seeds. Plant Physiol, 1985, 79, 879-884. (45) Yang, P. F.; Nouri, M. Z.; Komatsu, S. Proteomics and applications to food science in rice. In Proteomics in foods. Toldrá, F., Eds.; Springer: New York, USA, 2013, pp383-384 (46) Rajjou, L.; Duval, M.; Gallardo, K.; Catusse, J.; Bally, J.; Job, C.; Job, D. Seed germination and vigor. Annu. Rev. Plant Biol. 2012, 63, 507-533. (47) Reeves, R.; Nissen, M.S. The AT-DNA-binding domain of mam-malian high mobility group I chromosomal proteins. J. Biol. Chem. 1990, 265, 8573–8582 (48) Fujimoto, S.; Matsunaga, S.; Yonemura, M.; Uchiyama, S.; Azuma, T.; Fukui; K. Identification of a novel plant MAR DNA binding protein localized on chromosomal surfaces. Plant Mol. Biol. 2004, 56: 225–239. (49) Matsushita, A.; Furumoto, T.; Ishida, S.; Takahashi, Y. AGF1, an AT-Hook Protein, Is 30 ACS Paragon Plus Environment

Page 30 of 50

Page 31 of 50

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

Necessary for the Negative Feedback of AtGA3ox1 Encoding GA 3-Oxidase. Plant Physiol. 2007, 143, 1152-1162. (50) Schuhmann, H.; Adamska, I. Deg proteases and their role in protein quality control and processing in different subcellular compartments of the plant cell. Physiol. Plant. 2012, 145, 224-234. (51) Schuhmann, H.; Huesgen, P.; Adamska, I. The family of Deg/HtrA proteases in plants. BMC Plant Biology 2012, 12, 52. (52) Adam, Z.; Adamska, I.; Nakabayashi, K.; Ostersetzer, O.; Haussuhl, K.; Manuell,A.; Zheng, B.; Vallon,O.; Roderme, S. R.; Shinozaki, K.; Clarke, A. K. Chloroplast and Mitochondrial Proteases in Arabidopsis. A proposed nomenclature. Plant Physiol. 2001, 125, 1912-1918. (53) García-Lorenzo, M.; Sjödin, A.; Jansson, S.; Funk, C. Protease gene families in Populus and Arabidopsis. BMC Plant Biology 2006, 6, 30. (54) Tripathi, L. P.; Sowdhamini, R. Cross genome comparisons of serine proteases in Arabidopsis and rice. BMC genomics 2006, 7, 200. (55) Pendle, A. F.; Clark, G. P.; Boon, R.; Lewandowska, D.; Lam, Y. W.; Andersen, J.; Mann, M.; Lamond, A. I.; Brown, J. W. S.; Shaw, P. J. Proteomic analysis of the Arabidopsis nucleolus suggests novel nucleolar functions. Mol. Boil. Cell 2005, 16, 260-269. (56) Huesgen, P. F. Functional genomics of Deg and GCP proteases in photosynthetic organisms. Ph.D. Thesis. 2007, University of Konstanz, Germany. (57) Schuhmann, H. Deg proteasea in Arabidopsis thaliana. Ph.D. Thesis. 2008, University of Konstanz, Germany. (58) Kamakaka, R. T.; Rine, J. Sir-and silencer-independent disruption of silencing in Saccharomyces by Sas10p. Genetics 1998,149, 903-914. 31 ACS Paragon Plus Environment

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(59) Dragon, F.; Gallagher, J. E.; Compagnone-Post, P. A.; Mitchell, B. M.; Porwancher, K. A.; Wehner, K. A.; Wormsley, S.; Settlage, R. E.; Shabanowitz, J.; Osheim, Y.; Beyer, A. L.; Hunt, D. F.; Baserga, S. J. A large nucleolar U3 ribonucleoprotein required for 18S ribosomal RNA biogenesis. Nature 2002, 417, 967-970. (60) Sihn, C. R.; Lee, Y. S.; Jeong, J. S.; Park, K.; Kim, S. H. CANu1, a novel nucleolar protein, accumulated on centromere in response to DNA damage. Genes to Cells 2008, 13, 787-796. (61) Miller, J.; McLachlan, A. D.; Klug, A. Repetitive zinc-binding domains in the protein transcription factor IIIA from Xenopus oocytes. EMBO J. 1985, 4, 1609-1614. (62) Laity, J. H.; Lee, B. M.; Wright, P. E. Zinc finger proteins: new insights into structural and functional diversity. Curr. Opin. Struct. Biol. 2001, 11, 39-46. (63) Aravind, L. The BED finger, a novel DNA-binding domain in chromatin boundary element binding proteins and transposases. Trends Biochem. Sci. 2000, 25, 421-423. (64) Bendahmane, A.; Querci, M.; Kanyuka, K.; Baulcombe, D. C. Agrobacterium transient expression system as a tool for the isolation of disease resistance genes: application to the Rx2 locus in potato. Plant J. 2000, 21, 73-81. (65) Shen, Q. H.; Zhou, F.; Bieri, S.; Haizel, T.; Shirasu, K.; Schulze-Lefert, P. Recognition specificity and RAR1/SGT1 dependence in barley Mla disease resistance genes to the powdery mildew fungus. Plant Cell 2003, 15, 732-744. (66) Vizcaíno, J. A.; Côté, R. G.; Csordas, A.; Dianes, J. A.; Fabregat, A.; Foster, J. M.; Griss, J.; Alpi, E.; Birim, M.; Contell, J.; O’Kelly, G.; Schoenegger, A.; Ovelleiro, D.; Pérez-Riverol, Y.; Reisinger, F.; Ríos, D.; Wang, R.; Hermjakob, H. The PRoteomics IDEntifications (PRIDE) database and associated tools: status in 2013. Nucleic Acids Res. 2013, 41, D1063-D1069.

32 ACS Paragon Plus Environment

Page 32 of 50

Page 33 of 50

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

Figure legends Figure 1. Western blot analysis of nuclear proteins purified from rice embryos. Total (TP) and nuclear (N) proteins were purified from rice embryos after 24 h of imbibition. Proteins were separated by SDS-polyacrylamide gel electrophoresis and transferred onto a polyvinylidene difluoride membrane. The membrane was incubated with anti-histone or ascorbate peroxidase antibodies. Anti-rabbit IgG conjugated with horseradish peroxidase was used as the secondary antibody. Signals were detected using an ECL system. “M” denotes the marker lane. Coomassie brilliant blue (CBB) staining of the gel was used as a loading control and the relative intensities of bands were calculated. Data are means ± SD from 3 independent biological replicates (R1, R2, and R3). Asterisks indicate significant changes according to Student’s t-test (*P