Phosphoproteomic Analysis of Paper Mulberry Reveals

Mar 30, 2017 - On the basis of the protein–protein interaction network analysis, a protein kinase and phosphatases hub protein (P1959) were found to...
2 downloads 16 Views 3MB Size
Subscriber access provided by UB + Fachbibliothek Chemie | (FU-Bibliothekssystem)

Article

Phosphoproteomic Analysis of Paper Mulberry Reveals Phosphorylation Functions in Chilling Tolerance Zhi Pi, Meiling Zhao, Xian-Jun Peng, and Shi-Hua Shen J. Proteome Res., Just Accepted Manuscript • DOI: 10.1021/acs.jproteome.6b01016 • Publication Date (Web): 30 Mar 2017 Downloaded from http://pubs.acs.org on March 30, 2017

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 67

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

Phosphoproteomic Analysis of Paper Mulberry Reveals Phosphorylation Functions in Chilling Tolerance Zhi Pi,†,‡ Mei-Ling Zhao, †,‡ Xian-Jun Peng† and Shi-Hua Shen†,* †

Key Laboratory of Plant Resources, Institute of Botany, the Chinese Academy of

Sciences, Beijing 100093, China ‡

University of Chinese Academy of Sciences, Beijing 100049, China

KEYWORDS: chilling tolerance, phosphorylation, paper mulberry, hub protein, CBF/DREB-responsive pathway

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: Paper mulberry is a valuable woody species with a good chilling tolerance. In this study, phosphoproteomic analysis, physiological measurement and mRNA quantification were employed to explore the molecular mechanism of chilling (4 °C) tolerance in paper mulberry. After chilling for 6 hours, 427 significantly changed phosphoproteins were detected in paper mulberry seedlings without obvious physiological injury. When obvious physiological injury occurred after chilling for 48 hours, a total of 611 phosphoproteins were found to be significantly changed at the phosphorylation level. Several protein kinases, especially CKII, were possibly responsible for these changes, according to conserved sequence analysis. The results of GO analysis showed that phosphoproteins were mainly responsible for signal transduction, protein modification and translation during chilling. Additionally, transport and cellular component organization were enriched after chilling for 6 and 48 hours, respectively. Based on the protein-protein interaction network analysis, a protein kinases and phosphatases hub protein (P1959) was found to be involved in cross-talk between Ca2+, BR, ABA and ethylene-mediated signaling pathways. We also highlighted the phosphorylation of BpSIZ1 and BpICE1 possibly impacted on the CBF/DREB-responsive pathway. From these results, we developed a schematic for the chilling tolerance mechanism at phosphorylation level. INTRODUCTION Low temperature, as an extreme environment, is responsible for 30-40 % yield reduction in temperate growing areas.1 Along with temperature drop, the plants suffer 2

ACS Paragon Plus Environment

Page 2 of 67

Page 3 of 67

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

both water imbalance and oxidative damage, which are respectively induced by decreasing root hydraulic conductance and reactive oxygen species (ROS) production.2, 3

To adapt to chilling, complex mechanisms such as osmolyte biosynthesis and

transportation,4 membrane lipid unsaturation5 and ROS detoxification6 have evolved. These responses are due to the reprogramming of gene expression, which is mediated by Ca2+ and phytohormones such as ABA,7 auxin,8 BR9 and ethylene.10 Recently, the novel gene COLD1, which interacts with G-proteins leading to the activation of the Ca2+ channel for temperature sensing, has been identified in rice. A SNP of COLD1 caused the adaptation to a cold environment in japonica rice.11 In Arabidopsis thaliana, gene microarray studies reveal that the large majority (approximately 30 %) of the most highly chilling-induced genes belong to the CBF/DREB regulons, suggesting that they play a major role in configuring the low-temperature transcriptome.12, 13 The activity and stabilization of transcription factors such as ICE1 and ABI5 implicated in CBF/DREB

and

phosphorylation14,

ABA-responsive

pathways

are

strictly

dependent

on

15

, ubiquitination16 and SUMOylation.16 Thus, post-translational

modification is an essential strategy for regulation of cellular processes in response to chilling. Phosphorylation is one of most frequent post-translational modifications due to over thousand kinase genes occupying 3-4 % of functional genes in plant.17 Benefiting from the advancement of functional genomic and mass spectrum tools, identification and quantification of phosphorylation is feasible at large scale level. It has been widely 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

Page 4 of 67

employed for understanding phosphorylation networks in subcellular organelles, such as chloroplast,18 mitochondria19 and plasma membrane;20 downstream substrates of key protein kinases including STN8,21 SnRK22 and OXI1;23 and molecular mechanism of biotic/abiotic stress24,

25

and hormone response.26,

27

To understand how

phosphorylation functions in chilling response and defense, a total of 42 phosphoproteins in leaf sheath of rice have been explored by 2-DE with 32P-labeling visualization and mass spectrometric identification, of which 4 phosphoproteins involved in glycolytic metabolism processes and Ca2+-mediated signaling were demonstrated to be significantly changed after chilling.28 In the following years, 2-DE coupled with Pro-Q Diamond stain was employed to explore chilling-related phosphoproteins in rice roots29 and cell suspensions of A. thaliana.30 Although this information suggested the importance and ubiquity of phosphorylation in response to chilling, it merely focuses on model herbaceous plants and is still far from sufficient to understand the complex mechanism of chilling response owing to the limitation of 2-DE separation ability and staining sensitivity. The development of specific phosphopeptide enrichment and high resolution mass spectrum allows us to explore thousands of phosphoproteins simultaneously.18 Therefore, we applied the TiO2 enrichment and LTQ-Orbitrap mass spectrometry for exploring the phosphoproteome characterization of paper mulberry under chilling. Paper mulberry (Broussonetia papyrifera L.), which belongs Moraceae family, is a dioecious, perennial and woody species.31 The whole plant has high value in paper 4

ACS Paragon Plus Environment

Page 5 of 67

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

making, stock farming and medicine. Thus, most studies are mainly focus on analyzing fiber structure,32 nutrients33 and medicinal ingredients.34 However, its environmental adaptation characterization has been relatively ignored. It natively possesses wide temperature adaption owing to natural distribution in China from Hainan (19°12′N 109°42′E) to Beijing (39°55′N 116°23′E).35 Individuals are tolerant of frost and can overwinter in the place with -25 °C. Unfortunately, there are few studies focus on molecular mechanism of chilling response and defense in paper mulberry, including several previous studies on chilling-related transcriptomics and proteomics.36, 37 Low temperature is still the major limiting factor that significantly restrict the extension of the paper mulberry planting, which presents a challenge to the further production of the paper mulberry. In this study, quantitative phosphoproteomics has been applied to explore the chilling response and defense mechanism in paper mulberry at phosphorylation level. It not only broadens the understanding of how mulberry paper resists to low temperature stress, but also provide the opportunity to discover the candidate phosphorylated proteins and the potential pathway for chilling in paper mulberry.

MATERIALS AND METHODS Plant material and chilling treatments Paper mulberry (Broussonetia kazinoki × Broussonetia papyifera) were cultured on Murashige and Skoog medium in an artificial climatic chamber kept at 26 °C with a 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

14/10 hours photoperiod (day/night), quantum flux density of 80 µmol/(m 2 s). After 25 days, plantlets were transferred to 4 °C for chilling treatments. Only the last three fully expanded leaves of paper mulberry with the similar growth vigor were harvested in this study, and a mixed sampling strategy was adopted to eliminate differences between individuals. Both short-term (0, 6, 12, 24, 36 and 48 hours) and long-term (0, 2, 4, 6, 8 and 10 days) chilling treatments were set for detecting physiological changes. Samples were rapidly frozen in liquid nitrogen and stored in -80 °C for further protein and RNA extraction. Three treatments (0, 6 and 48 hours treatment) were also prepared for protein extraction. Additionally, the leaves with 0, 2, 6, 12 and 24 hours chilling were sampled for RNA-seq analysis. For quantitative RT-PCR analysis, chilled paper mulberry leaves were respectively sampled at different time points (0, 1, 3, 6, 12, 24, 24 and 48 hours). Three independent biological replicates of control and treated samples were analyzed for physiological and phosphoproteomic analysis. RNA-seq analysis was performed without replicate, while the expression of interesting genes was checked by Q-PCR with three biological replicates. Physiological characteristics measurement Chlorophyll fluorescence was analyzed using a MAXI-IMAGING-PAM chlorophyll fluorometer (Heinz Walz, Effeltrich, Germany). Briefly, the whole plantlets were dark adapted for 30 minutes after treatment, and the maximum quantum yield of PSII (Fv/Fm) was determined for each sample by analyzing fully expanded intact leaves. The method described by Yan et al.38 was employed to measure relative electrolytic 6

ACS Paragon Plus Environment

Page 6 of 67

Page 7 of 67

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

leakage (REL). Four leaves from each treatment were cut off and immersed in tubes with 10 mL deionized water. The tubes were placed on a shaker at 28 °C. Two hours later, the electrical conductivity of bathing solution (Lt) was measured by conductivity meter (INESA instrument, Shanghai, China). Then the tubes were incubated in boiling water for 20 minutes and subsequently cooled down to 28 °C. The electrical conductivity (Lo) was measured again. The relative electrolyte leakage was calculated by the formula: REL= Lt/Lo × 100 Malondialdehyde (MDA) content was determined by the thiobarbituric acid reaction. 39

First, 0.3 g leaves was ground with an ice-cooled mortar and pestle and then extracted

with 3 mL 10 % (w/v) trichloroacetic acid. The homogenate was then centrifuged at 8,000 g for 10 minutes at 4 °C. Then, 3 mL supernatant and 3 mL 0.05 % (w/v) TBA were mixed and heated for 30 min in a boiling water bath and then cooled on ice to stop the reaction. The concentration of MDA was calculated as follows: MDA=6.45(A532 − A600) − 0.56A450 Protein extraction and digestion Protein was extracted using a modified method according to Shen et al.40 Briefly, 0.5-1 g leaves were homogenized in 2 mL of the homogenization buffer containing 20 mM Tris-HCl (pH 7.5), 250 mM sucrose, 10 mM ethylene glycol tetraacetic acid, 5 % 2-mercaptoethanol, 1 % Triton X-100, 1 mM phenylmethylsulfonyl fluoride and PhosSTOP Phosphatase Inhibitor Mixture (Roche, Mannnheim, Germany). Then, the homogenate was centrifuged at 10,000 g for 10 minutes at 4 °C. The supernatant was 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

mixed with a quarter of the volume of 50 % (w/v) trichloroacetic acid to attain 10 % (w/v) trichloroacetic acid in final mixture. After incubating in an ice bath for 30 minutes, the mixture was centrifuged at 15,000 g for 10 minutes at 4 °C, and the supernatant was discarded. After washing three times with acetone, the pellet was collected by centrifugation and air dried. Then, 3 mg total protein from each sample was resuspended in 1500 µL of 4 % SDS and 0.1 M DTT in 0.1 M Tris-HCl, pH 7.6, at room temperature, and following the filter-aided sample preparation (FASP) protocol described elsewhere.41 Briefly, for every 500 µg protein, 500 µL of 8 M urea in 0.1 M Tris/HCl, pH 8.5 (UA) was added, and samples were centrifuged at 12,000 g at 20 °C for 30 minutes. This step was performed three times. Then, 100 µL of 0.05 M iodoacetamide in 8 M urea was added to the filters, and the samples were incubated in darkness for 20 minutes. Filters were washed twice with 400 µL of 8 M UA, followed by two washes with 400 µL of 50 mM NH4HCO3. Finally, 1 µg/µL trypsin (Promega, Madison, USA) was added in 150 µL of 50 mM NH4HCO3 to each filter. The protein to enzyme ratio was 100:1 (w/w). Samples were incubated overnight at 37 °C, and released peptides were collected by centrifugation. Phosphopeptide enrichment The phosphopeptide enrichment procedure was described elsewhere with some modification.42, 43 In detail, the TiO2 beads (GL Sciences, Tokyo, Japan) were incubated in 800 µL loading buffer containing 65 % ACN and 2 % TFA saturated by glutamic acid. Tryptic peptides (3mg) were resuspended in 800 µL loading buffer. After shaking the 8

ACS Paragon Plus Environment

Page 8 of 67

Page 9 of 67

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

mixture of tryptic peptide and TiO2 beads (1:1, w/w) for 40 minutes and collecting TiO2 beads by centrifugation, 800 µL wash buffer (65 % ACN and 0.1 % TFA) was added and shaken for 40 minutes twice. Subsequently, TiO2 beads were eluted twice with 800 µL elution buffer (500 mM NH4OH and 60 % ACN). The eluates were dried and reconstituted in 0.1 % formic acid for MS analysis. LC-MS/MS Analysis LC was performed on an Easy-nLC System (Thermo Fisher Scientific, Odense, Denmark). Peptides were separated on a 15-cm fused silica emitter packed in-house with the reverse phase material ReproSil-Pur C18AQ, 3 µm resin (Dr. Maisch, Ammerbuch, Germany) with a 55 minutes gradient at 300 nL/min from 5 % to 22 %, and 5 minutes from 22 % to 90 % of 100 % ACN/0.1 % acetic acid (v/v). The quadrupole Orbitrap mass spectrometer (Thermo Fisher Scientific, Bremen, Germany) was operated in the positive ion mode using a data dependent “top 12” method. Survey scans (m/z 300-1500) and MS/MS scans were acquired at a resolution of 70,000 and 17,500 at m/z 200. Up to the top 12 most abundant isotope patterns with charge ≥ 2 from the survey scan were selected with an isolation window of 2 Th and fragmented by HCD with normalized collision energies of 25. For accurate mass measurements, the lock-mass option was employed (445.120025).44 The maximum ion injection times for the survey scan and the MS/MS scans were 50 ms and 100 ms, respectively, and the ion target values were set to 1E6 and 1E5, respectively. Selected sequenced ions were dynamically excluded for 30 seconds. 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 and phosphopeptide quantification Raw mass spectrometric data were analyzed with the MaxQuant software (version 1.5.2.8).45 A time-dependent mass recalibration algorithm was used for recalibration to improve the mass accuracy of precursor ions.46 MS/MS spectra were searched by the Andromeda search engine, which is incorporated into the MaxQuant software suite, against the paper mulberry protein database (28,908 entries), whose protein sequences were predicted based on genome and transcriptome sequencing. For the search, trypsin allowing for cleavage N-terminal to proline (trypsin/P) was chosen as enzyme specificity. Cysteine carbamidomethylation was selected as a fixed modification, while protein N-terminal acetylation, methionine oxidation and phosphorylation on serine/threonine/tyrosine were selected as variable modifications. Maximally, two missed cleavages were allowed. For MS and MS/MS, the tolerances of the main search for peptides were set at 7 ppm and 20 ppm, respectively. A false discovery rate (FDR) of 0.01 for proteins and peptides and a minimum peptide length of 7 amino acids were required. Identified peptides that are all shared between two proteins were combined and reported as one protein group. Proteins matching to the reverse database were filtered out. For high-confident phosphorylation site localization, the following two criteria were defined: their localization probability for the assignment was at least 0.75, and the PTM score difference from the second possible localization assignment was 5 or higher. The median localization probability of such class I sites is almost always higher than 0.99.47 Phosphopeptides were quantified by their extracted ion 10

ACS Paragon Plus Environment

Page 10 of 67

Page 11 of 67

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

chromatograms (‘Intensity’ in MaxQuant),47,

48

enabling the ‘match between runs’

option (time window 2 minutes). Then, normalization followed with imputation as Karpievitch et al.49 described were performed to remove systematic biases and obtain a complete matrix of intensities. Briefly, phosphopeptide intensity distributions were log-transformed to base 2, normalized across biological replicates by linear regression, and standardized by median absolute deviation and mean centering across samples with InfernoRDN.50 Subsequently, missing intensity value was imputed by lowest observed value of each treatment. The number of imputations for each phosphopeptide was restricted to three. Only the phosphopeptides that met the following restrictions were regarded as significantly changed phosphoproteins: (1) phosphorylation site localization probability greater than 0.75 and phosphorylation site score difference not less than 5; (2) the p value of one-way ANOVA analysis < 0.05; and (3) the fold change between each treatments > 2. The mass spectrometry proteomics data were deposited to the ProteomeXchange Consortium via the PRIDE partner repository with the dataset identifier PXD005516.51 Transcriptomics analysis The transcriptomics analyses were prepared according to our former study.52 Briefly, total RNAs were isolated with TRIzol Reagent (Life Technologies, Shanghai, China) from each sample according to the manufacturer’s instructions. RNAs were treated with RNase-free DNase I (Takara, Dalian, China) to remove the DNA residues. The quality and purity of RNAs were assessed with OD260/230 ratio and RNA integrity 11

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

Page 12 of 67

number (RIN) using the NanoDrop 2000 (Thermo Fisher, Waltham, USA) and the Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, USA), respectively. mRNAs were purified from the mixed high-quality total RNA using Oligo (dT) RNA purification. First strand cDNAs were synthesized using Oligo (dT) primer, then second strand cDNAs were synthesized using RNase H and DNA polymerase I. Two different adapters were ligated to the cDNA fragments, which were then purified using second strand AMPure XP beads (Beckman Coulter, California, USA). cDNA templates were enriched by multiplexing PCR with primers 1.0 and 2.0. Barcodes were added to each library during amplification. Raw sequence data were generated by the Illumina pipeline and are available in NCBI’ s Short Read Archive (SRA) database under accession number SRP029966. All of the clean reads were pooled together and assembled to form the global transcriptome of the paper mulberry. To obtain the transcriptome of each sample, raw data were also assembled for each sample. The expression level of every transcript in each sample was calculated by quantifying the number of Illumina reads that mapped to the paper mulberry transcriptome. The raw gene expression counts were normalized using the RPKM (reads per kb per million reads) value. For screening of differentially expressed genes, p values that correspond to differentially expressed genes (DEGs) were obtained via a general

Chi

squared

test

that

was

performed

using

IDEG6

(http://telethon.bio.unipd.it/bioinfo/IDEG6/). The threshold of p values in multiple tests was checked through manipulating the false discovery rate (FDR) value. Among 12

ACS Paragon Plus Environment

Page 13 of 67

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 five samples, the transcripts with the highest RPKM value more than 5 and the ratio of RPKM between samples of more than 3 (fold change > 3) and an FDR < 0.01 were considered to have significant changes in expression in response to chilling. Quantitative RT-PCR analysis Total RNA was isolated from leaf tissue using TransZol Plant Mini kit (TransGen, Beijing, China), and the DNA was digested using DNase I (Takara, Dalian, China) according to the manufacturer's instructions. RNA integrity was checked on an agarose gel and quantified by using the NanoDrop 2000 (Thermo Fisher, Waltham, USA ). The first-strand cDNA for each sample was prepared from 1 µg of total RNA using SuperScript II reverse transcriptase (Takara, Dalian, China) following the manufacturer’s instructions, and template cDNAs were diluted 5-fold for PCR. Gene-specific primers were designed using the Primer 5 program and are listed in Table S-1. Quantitative RT-PCR analysis was performed by our former study.36 Samples and standards were run in triplicate on each plate using the SYBR® Premix Ex Taq™ II kit (Takara, Dalian, China) on a StepOneTM real-time PCR system (Applied Biosystems, Foster City, USA) following the manufacturer’s recommendations. Real-Time PCR was performed in a 20 µL reaction volume containing 6.8 µL of ddH2O, 10 µL of SYBR® Premix Ex Taq II, 0.4 µL of ROX Reference Dye II, 0.4 µL of forward primer (10 µmol/L), 0.4 µL of reverse primer (10 µmol/L), and 2 µL of template cDNA. The PCR programs were run as follows: 30 seconds of pre-denaturation at 95 °C, 40 cycles of 5 seconds at 95 °C and 30 seconds at 60 °C, followed by steps for dissociation curve 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

Page 14 of 67

generation. Dissociation curves for each amplicon were carefully examined to confirm the lack of multiple amplicons at different melting temperatures. Relative transcript levels for each sample were obtained using the comparative cycle threshold method using the cycle threshold value of glyceraldehyde-3-phosphate dehydrogenase gene as a control. Bioinformatics Protein

function

was

annotated

through

Blast2GO

software

(https://www.blast2go.com/home). Enrichment analysis was conducted using the Singular

Enrichment

Analysis

(SEA)

tool

in

the

AgriGO

toolkit

(http://bioinfo.cau.edu.cn/agriGO/).53 Customized annotated reference was selected and uploaded based on our genomic sequencing. Under advanced options, the gene ontology type chosen was Plant GO slim. Significantly changed phosphorylation sites were subjected to k-means clustering using MultiExperiment Viewer (MeV). The Euclidean distance of log-ratios normalized to their respective maxima per phosphorylation site was used to measure the distance between regulated phosphorylation sites. The iceLogo program (version 1.2) was used for the analysis of sequence features adjacent to the significantly changed phosphorylation sites.54 The 15 amino acids surrounding each phosphorylation site were extracted and aligned as an experimental set. The reference dataset consisted of non-phosphorylated peptides, which were generated by random sampling from the full sequences of significantly changed 14

ACS Paragon Plus Environment

Page 15 of 67

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

phosphoproteins as Schmidt et al.55 described. In addition, another reference dataset containing all of the non-changed phosphopeptides that were identified in our mass spectrometric analysis was used for calculation. The analysis was visualized using iceLogo figures with percentage difference as scoring system (p < 0.05). The protein-protein interaction of the phosphoproteins was analyzed by the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) (http://string-db.org/) database of physical and functional interactions. The confidence score was set at the high level (> 0.700). The protein-protein interaction network was arranged and displayed using Cytoscape (version3.3.0) software (http://www.cytoscape.org/). Protein domains were searched by Pfam (http://pfam.xfam.org/) and PlantsP (http://plantsp.genomics.purdue.edu/index.html).

RESULTS Physiological changes of paper mulberry after chilling After 25 days growth, paper mulberry plantlets were transferred to 4 °C for 10 days chilling treatment. The growth of paper mulberry was hardly inhibited by chilling (Figure 1A). A slight wilting of leaves occurred after chilling for 2 days and then was gradually aggravated. However, necrotic areas and external discoloration were not caused. According to chlorophyll fluorescence (Figure 1B), leaf age and position were correlated with chilling sensitivity. The top three fully expanded leaves of paper mulberry were highly sensitive to chilling and were selected for further measurement. 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

In contrast, the older leaves located at the bottom did not show obvious change until 4 days after chilling. The change of Fv/Fm value in the top three fully expanded leaves was calculated and shown in Figure 1C. It was markedly suppressed with prolonged chilling time, and the Fv/Fm nearly reached the lowest point (approximately 0.2) after chilling for 6 days. To understand when PSII was first damaged by chilling, a short-term series of chilling treatments was performed. The results show that no significant change occurred within 1 day, and values of Fv/Fm were decreased from 0.827 to 0.631 after chilling for 2 days. The value of relative electrolytic leakage soared from 8.29 % to 86.99 % during chilling for 10 days (Figure 1D). The significant increase was first observed at 1 day after chilling. It is worth noting that the degree of relative electrolytic leakage increase within the last 4 days after chilling was nearly double times that of the first 6 days. These findings implied that the plasma membrane of paper mulberry was relatively stable within 1 day after chilling, and drastic collapse of plasma membrane occurred after chilling for 6 days. MDA content of leaves presented a gradual rising tendency as chilling time increased (Figure 1E). The significant increase of MDA content was detected within 1 day after chilling, which was earlier than plasma membrane leakage and photoinhibition. At last, MDA content had been increased approximately 2.4-fold comparing to control. Comparative analysis of paper mulberry phosphoproteomics after chilling. 16

ACS Paragon Plus Environment

Page 16 of 67

Page 17 of 67

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

To investigate the phosphorylation involved in chilling tolerance, we evaluated the changes in the phosphoproteome of paper mulberry leaves after chilling for 6 and 48 hours, when non-significant and significant physiological injury was detected. In total, 2,478 phosphopeptides with 3105 phosphorylation sites corresponding to 1,509 phosphoproteins were identified (Table S-2), and the proportions of pS, pT, and pY sites were estimated as 84.7 %, 15.0 %, and 0.3 %, respectively. To evaluate the quality of the complete experiment, correlation analyses for replicates and PCA analysis were performed. Three samples were individually clustered with their replicates with Pearson correlation coefficients between 0.74 and 0.93 (Figure S-1). A total of 933 significantly changed phosphorylation sites corresponding to 719 phosphoproteins were screened out (Figure 2A). Compared to untreated leaves, 509 and 768 phosphorylation sites corresponding to 427 and 611 phosphoproteins were differentially changed after chilling for 6 hours and 48 hours, respectively (Table S-3 and 4). There were 253 and 256 phosphorylation sites up-regulated and down-regulated, respectively, at 6 hours under chilling (Figure 2B). After chilling for 48 hours, a total of 628 and 140 phosphorylation sites were respectively up-regulated and down-regulated at the phosphorylation level (Figure 2B). Between 6 hours treatment and 48 hours treatment, there were 683 phosphorylation sites corresponding to 555 phosphoproteins screened out as significantly changed at the phosphorylation level (Table S-5), of which 568 phosphorylation sites were up-regulated and 115 phosphorylation sites were down-regulated at the phosphorylation level. A total of 16, 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

22 and 168 phosphorylated sites were uniquely detected after chilling for 0, 6 and 48 hours, respectively (Figure 2B). Then, all of the significantly changed phosphorylation sites were classified into eight groups with four up-regulated and down-regulated patterns (Figure 2C). IceLogo analysis of phosphorylation sites in each cluster. To eliminate any amino acid frequency bias, non-phosphorylated peptides from full-length phosphoproteins were used as a reference dataset (Figure 3A). Two clear consensus motifs ([SP] and [RxxS]) were enriched in the whole clusters and occupied a relatively high percentage. An acidic motif ([SxE]) was also observed in Cluster 4. However, only [SxE] was overrepresented compared to non-changed phosphopeptides (Figure 3B). Additionally, histidine at position -2 presented a preference of Cluster 1, 2, 4 and 7 against both different reference datasets. Glycine residue at position -1 was enriched in Cluster 5 and 8, containing down-regulated phosphopeptides after chilling. Of note, the percentage difference of these residues compared with two reference datasets was always below 10 %, and some did not belong to clear consensus motifs. GO annotation and enrichment analysis of significant changed phosphoproteins Blast2GO and AgriGO were applied to annotate and analyze significantly changed phosphoproteins. From the biological process perspective, protein modification, signal transduction and translation were significantly overrepresented after chilling for 6 and 48 hours (Figure 4A). Amongst protein modification, 77 protein kinases and 15 phosphatases possessed over 70 % of phosphoproteins (Table 1 and Figure 5A). For 18

ACS Paragon Plus Environment

Page 18 of 67

Page 19 of 67

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

signal transduction, the ABA, BR, light, Ca2+ and phosphatidylinositol-mediated signaling accounted for over half of signal transduction terms under 48 hours treatment (Figure 5B). In terms of translation, phosphorylation mainly happened in ribosomal proteins, RNA binding proteins and translational initiation factors after chilling (Figure 5C). Transport was merely enriched after chilling for 6 hours. The majority was protein and ion transport, followed by hormone and water transport (Figure 5D). With chilling duration, cellular component organization was also overrepresented. Cellular component organization was found to contain up to 22 % cytoskeleton organization, when chilling injury was aggravated after 48 hours (Figure 5E). The other phosphoproteins were also involved in the term of cell morphogenesis, chromatin organization, protein complex organization and chloroplast organization, which had 16 %, 15 %, 13 % and 12 %, respectively, of cellular component organization. In addition, protein modification processes were enriched in both Cluster 3 and Cluster 4. Cellular component organization and signal transduction were also overrepresented in Cluster 3. From the molecular function perspective, kinase activity as well as binding for RNA, protein, nucleotide and lipid was overrepresented (Figure 4B). Among these, kinase activity and protein binding were enriched after chilling for 6 hours and 48 hours. Additionally, lipid binding was overrepresented after chilling for 6 hours, and binding for nucleotide and RNA were only enriched after chilling for 48 hours. Regarding cellular component perspective for 6 hour and 48 hour chilling treatments, significantly 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

changed phosphoproteins were mainly located in the nucleus, plasma membrane, cytosol, thylakoid and ribosome. Transcriptional expression analysis of paper mulberry after chilling After stringent quality checks and data cleaning, we obtained 12, 992, 767, 550 raw reads, totaling 9.164 Gb. The sequencing data output and quality assessment were listed in Table S-6. By applying screening thresholds of 3-fold changes and FDR < 0.01, a total of 5,800 DEGs were detected as being involved in chilling tolerance (Figure 6A). Then, a total of 274 differentially expressed protein kinases and phosphatases were screened out at the transcriptional level in response to chilling, accounting for approximately 5 % of total DEGs (Figure 6B). Relative to the control, 24 protein kinases and phosphatases with 13 up-regulated genes and 11 down-regulated genes were detected after chilling for 2 hours. There were 203 protein kinases and phosphatases detected between control and 6 hours treatment. Among these, 158 genes were up-regulated and 45 genes were down-regulated. In addition, 25 and 40 protein kinases and phosphatases were up-regulated after chilling for 12 and 24 hours, respectively, while 12 and 9 protein kinases and phosphatases were down-regulated. To understand the transcriptional expression pattern corresponding to significantly changed phosphoproteins, RNA-seq of the paper mulberry transcriptome under chilling was performed. It is notable that the expression patterns of 706 corresponding mRNAs contained 554 non-differentially and 152 differentially expressed mRNAs was 20

ACS Paragon Plus Environment

Page 20 of 67

Page 21 of 67

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

determined (Figure 6C and Table S-7). Subsequently, q-PCR was performed to reassess the expression of interesting genes (Figure 6D). The expression of 11 genes, including KEG, SIZ1 and PLC2, was not induced by chilling, while corresponding proteins were significantly changed at the phosphorylation level. The up-regulation of the BSL3 gene occurred after chilling for 24 hours, and changes at the phosphorylation level were detected after chilling for 6 hours. The transcriptional expression of CAT was down-regulated after chilling, while the corresponding protein was significantly phosphorylated. Additionally, the expression pattern of BpDREB1 was analyzed, which increased to the maximum within 6 hours after chilling, and then drastically decreased (Figure S-2). Protein-protein interaction analysis of protein kinases and phosphatases Considering the essential role of protein kinases and phosphatases, 363 protein kinases and phosphatases significantly changed in response to chilling at the transcriptional or phosphorylation levels were used for establishing a protein-protein interaction network (Table 1 and S-8). Finally, a protein-protein interaction network containing 93 nodes and 303 edges was demonstrated through Cytoscape (Figure 7). These proteins showed relationships with MAPK cascades as well as light, ABA, cyclin, ethylene, BR and Ca2+-mediated signaling pathways. A potential hub protein (P1959) was noted in the center of network with 52 first neighbors; these neighbors were spread in a wide distribution of signal transduction, including MAPK cascades as well as ABA, ethylene, BR and Ca2+-mediated signaling pathway. Among these, STN8 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

with 25 neighbors was associated with cyclin-mediated signaling pathway. For further understanding the features of P1959, functional domain analysis was performed. P1959 contained three domains, namely a protein phosphatase 2C domain, a cyclic nucleotide-binding domain and a protein kinase domain located in order from the N-terminal to C-terminal (Figure S-3). To estimate conservation of phosphorylation sites, sequence alignment was computed with homologous sequence of Morus alba, P. trichocarpa, A. thaliana and Glycine max. The result indicated that three domains and phosphopeptides demonstrate high conservation in all five species.

DISCUSSION Phosphorylation was essential for chilling tolerance in paper mulberry Unlike chilling sensitive crops such as rice, cotton and maize, there was no obvious physiological injury of cellular membrane and chloroplast occurring within 24 hours (Figure 1). The Fv/Fm of chilling sensitive crops such as rice, maize and cotton were generally decreased over 50 % at 7 °C within 20 hours.56 The significant increase of REL and MDA were also reported within several hours after chilling.38, 57 These results suggested that paper mulberry had good chilling resistance, although it failed to resist long-term chilling. Through a combination of phosphoproteomic and transcriptomic approaches, we performed a comprehensive analysis of chilling tolerance in paper mulberry leaves for the first time. Phosphorylation changes widely occurred in chilling tolerance and the 22

ACS Paragon Plus Environment

Page 22 of 67

Page 23 of 67

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

numbers drastically increased with chilling duration, especially after physiological injury occurrence (Figure 2). It is notable that a high proportion of significantly changed phosphoproteins were not induced at the transcriptional level (Figure 6C), indicating the essential role of phosphorylation in chilling tolerance. Subsequently, we employed iceLogo for sequence analysis against non-phosphorylated peptides and non-changed phosphopeptides (Figure 3). The most obvious difference between the two methods was that [SP] and [RxxS] were highly enriched against with non-phosphorylated peptides, but not with non-changed phosphopeptides. A possible reason was that the MAPKs,58 CDKs,58 or CDPKs59 recognizing [SP] or [RxxS] might catalyze not only chilling-responsive proteins but also other proteins. In contrast, [SxE] was constantly overrepresented by different references. Thus, CKIIs might act as an important protein kinase causing increasing protein phosphorylation under chilling.60 Previous phosphoproteomics studies in abiotic stress have showed that phosphorylation changes mainly happened in regulatory proteins involving signal transduction, transcription, translation and transport.61, 62 Similarly, phosphorylation events were implicated in these functional categories against chilling. In addition, chilling-induced phosphorylation contributed to cellular component organization and protein modification (Figure 4). A recent phosphoproteomics study in a chill-tolerant fly showed that the term of cytoskeleton organization was highly enriched, and phosphorylation may play a role in cytoskeleton depolymerization.63 Together with the above phosphoproteomic results, a complex cellular network regulated by 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

phosphorylation

was concluded,

as

illustrated

in

Page 24 of 67

Figure

8.

Briefly,

the

chilling-responsive signaling was transduced through protein phosphorylation, which was invoked by phytohormones, Ca2+ light and phosphoinositide. The expression of chilling-related proteins was controlled by phosphorylation, mainly involving ribosome and translation initiation factors. During early chilling, phosphorylation also played an important role in the regulation of transport for homeostasis. With physiological injury caused by chilling, the proteins related to cytoskeleton and chloroplast organization were phosphorylated against chilling. A hub protein (P1959) was associated with a number of chilling-responsive signaling pathways. The term of protein modification was highly enriched after chilling for both 6 and 48 hours, containing over 70 % proteins that participated in protein phosphorylation or dephosphorylation (Figure 4A and 5A). Interestingly, a protein kinase and phosphatase hub protein (P1959) interacted with dozens of protein kinases and phosphatases including CIPK23, SRK2E, KEG, BSL1 and CTR1. Its A. thaliana orthologue (AT2G20050) has also been noticed as hub protein among protein kinases in response to oligogalacturonides, which is a class of plant damage-associated molecular patterns.20 The N-terminal and C-terminal part of P1959 exhibited homology with protein phosphatase 2C and cGMP/cAMP-dependent protein kinase, respectively (Figure S-3). In plants, they are considered essential components of Ca2+, ABA, cAMP and cGMP-mediated signaling pathways.64-66 This could account for the interaction of 24

ACS Paragon Plus Environment

Page 25 of 67

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

P1959 with such a large number of protein kinases and phosphatases. However, there were few studies of its molecular functions hitherto, and the impact of P1959 phosphorylation on the protein kinase and phosphatase network is still worthy of future study. The neighbors of P1959 were associated with a number of chilling-responsive signaling and regulated multiple biological processes. When plants are exposed to chilling temperature, Ca2+ is generally accumulated in cytosol67 and regulates several metabolism processes, including ROS scavenging68 and K+ transport.69 The phosphoproteomic analysis in rice has been investigated, and phosphorylation of OsCDPK13 is an important signaling component in the response to cold stress.28 Consistent with this, two BpCDPKs (P315 and P847) were found to be up-regulated at the phosphorylation level after chilling. For the ROS scavenging system, two antioxidant enzymes, namely CAT (P1091) and APX (P1886), were phosphorylated in response to chilling, whose phosphorylation sites belong to [K/RxxS]. This indicated that CDPKs as possible upstream protein kinase play a role in controlling the activity of antioxidant enzymes against chilling. BpCIPK23 (P1406) was also found to exhibit phosphorylation of Thr 190 located within the activation loop. It has been reported to function in phosphorylation of the K+ channel and contributes to K+ uptake.70 It might be responsible for the phosphorylation of two potassium transporters (P1512 and P1669) detected after chilling for 48 hours. In addition to Ca2+-mediated signaling pathway, phosphorylation changes were implicated in phytohormone-mediated 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

signaling. The BpSRK2E (P2045), BpSRK2I (P1113) and BpKEG (P2676) involved in ABA-mediated signaling pathway were observed to be up-regulated after chilling for 48 hours. The phosphorylation sites of BpSRK2E and BpSRK2I were positioned in the activation loop, where they are possibly responsible for activation of downstream transcription factors.71 BpKEG, as a negative regulator, would be self-degraded after phosphorylation.72 In BR and ethylene-mediated signaling pathway, several protein kinases and phosphatases were identified and shown to be phosphorylated. Although whether these phosphorylation sites were regulatory was unclear, the downstream substrates BpBES1/BZR1 (P675)73 and BpEIN2 (P428)74 exhibited phosphorylation changes at regulatory sites. This finding implied that phosphorylation was implicated and activated these signaling pathways in response to chilling. Phosphorylation was involved in chloroplast movement against chilling In the protein kinase and phosphatase network, a small group contained two phototropins associated with light-mediated signaling (Figure 7). In A. thaliana, phototropins are considered to mediate not only phototropism and stomatal opening, but also chloroplast movement.75, 76 Chloroplast movement is an efficient strategy to alleviate PSII damage and maintain maximal photosynthetic output.77 In plants, many studies have revealed that low temperature can facilitate light-induced chloroplast movement, which is one of the vital mechanisms for chilling tolerance plants.78, 79 In current study, mass spectrometry showed that multiple phosphorylation sites of 26

ACS Paragon Plus Environment

Page 26 of 67

Page 27 of 67

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

BpPHOT1 (P285) and BpPHOT2 (P2400) were up-regulated at the phosphorylation level after chilling for 6 and 48 hours. Autophosphorylation is frequently observed accompanied with activation of phototropins in A. thaliana.80 For AtPHOT1, residues in the hinge region between LOV1 and LOV2 are significantly more phosphorylated under strong light than in darkness.81 A total of five chilling-induced phosphorylation sites in BpPHOT1 were located in the hinge region between the LOV1 and LOV2 domain (Figure S-4). Among these, Ser 386 and Ser 426, which were respectively found to correspond with Ser 350 for AtPHOT1 and Ser 358 for VfPHOT1, underwent autophosphorylation by blue light. The phosphorylation at these sites was also essential for 14-3-3 protein binding.82 The 14-3-3 proteins was thought to bind to the phosphorylated target leading to a change in structure that regulates activity.83 In addition to signaling response, chloroplast movement depends on different numbers of chloroplast-actin filaments on the front and rear regions of the moving chloroplast.84 After chilling, we observed chilling-induced change of CHUP1 (P1546), PMI1 (P319), JAC1 (P2615) and KCA2 (P2751) at the phosphorylation level. They are implicated in multiple steps of chloroplast movement, including generation, maintenance and reorganization of chloroplast-actin filaments.85-87 Although there is relatively little physiological and molecular evidence for explaining chloroplast movement through phosphorylation, chloroplast movement-related proteins were frequently phosphorylated under strong light, suggesting phosphorylation plays an important role.81 AtKCA1 can be phosphorylated at Ser 842 and Ser 846 by CDKA;1, 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

Page 28 of 67

causing a conformational change in the structure of KCA with implications in folding and dimerization.88 Two significant changed phosphorylation sites of BpKCA2 (Ser 680 and Ser 685) were located in the coiled-coil stalk domain, which contributes to KCA homo- or hetero- dimerization. The phosphorylation site of BpPMI1 was close to the N-terminal side of the cleavage signal, which was only phosphorylated after chilling for 48 hours. It suggested that phosphorylation possibly played a role in PMI1 translocation. Together with these, chloroplast movement could be modulated by phosphorylation and partly account for the gentle decrease of Fv/Fm in the paper mulberry after chilling (Figure 1C). Phosphorylation participated in the regulation of phosphoinositide metabolism. Phosphoinositides, although occupying a small fraction of the plasma membrane, are believed to play key roles in signaling transduction and cytoskeleton organization.89 In the mulberry tree, an increase of phosphoinositides was observed from autumn to winter, suggesting it was closely associated with cold acclimation.90 A total of eight participants in phosphoinositide metabolism were differentially phosphorylated. They were responsible for reactions leading to the generation of different phosphoinositide species. Amongst them, PLC is capable of catalyzing the hydrolysis of phosphatidylinositol

4,5-bisphosphate

into

diacylglycerol

and

inositol

1,4,5-trisphosphate regarded as two important second messenger molecules. The production can not only result in release of Ca2+ into the cytoplasm but also cytoskeleton rearrangement and growth regulation.91, 92 Recent studies have shown that 28

ACS Paragon Plus Environment

Page 29 of 67

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

AtPLC2, as a primary phosphoinositide-specific phospholipase C, contributes to phosphoinositide metabolism, which is not induced by environmental stresses such as dehydration, salinity and low temperature.92, 93 In paper mulberry, no obvious change in the expression of BpPLC2 (P482) was detected during chilling (Figure 6D). Two residues (Ser 270 and Ser 357) were significantly phosphorylated in response to chilling, of which Ser 357 was was located in the catalytic Y domain (Figure S-4). Similarly, corresponding phosphorylation sites have been detected in the catalytic Y domain

of

AtPLC2,

suggesting

regulation

of

enzymatic

activity

though

phosphorylation.93 The present study also showed that phosphoinositide phosphatase SAC2 (P1500) was significantly phosphorylated after chilling for 6 hours. Along with continuous chilling, significant phosphorylation of BpSAC1 (P723) and BpSAC3 (P1569) were detected. The mutant of AtSAC1 is found to cause aberrant organization of the actin cytoskeleton.94 The phosphorylation site of three phosphoinositide phosphatase SACs were located in SAC domain and nearby phosphatase catalytic core, suggesting effect on its activity. Coinciding with phosphoinositide metabolism, the phosphorylation of two phosphatidylinositol transfers (P1963 and P2583) were respectively up-regulated after chilling for 6 and 48 hours, which was required to deliver phosphatidylinositol to plasma membrane.95 It has been reported that phosphorylation of soybean phosphatidylinositol

transfer

alter

interaction

with

membrane

leading

to

phosphatidylinositol transport.96 Thus, phosphorylation was seen as a determinant of 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

Page 30 of 67

phosphoinositide metabolism and responsible for chilling perception and cytoskeleton organization. Phosphorylation of ICE1 and SIZ1 impacted on BpDREB1 expression. The

phosphorylation

events

in

the

regulatory

components

of

the

CBF/DREB-responsive pathway also caught our attention due to its essential role in chilling tolerance. SIZ1 interacts with ICE1, enhances its stabilization and then facilitates expression of CBF3/DREB1A.97 However, the regulation of SIZ1 was still enigmatic. AtSIZ1 expression was constitutive and unaffected by chilling.97 Consistent with this, the expression of BpSIZ1 in paper mulberry was also not induced by chilling (Figure 6D). Recent study in overexpression and antisense MdSIZ1 mutants has been reported that the SUMO conjugation level does not depend on transcriptional and translational expression.99 These studies implied that the activity of E3 SUMO ligase SIZ1 might be regulated at post-translational level. Phosphoproteomic analysis showed that Ser 495 and Ser 497 of BpSIZ1 (P2340) were significantly phosphorylated in response to chilling. These sites were adjacent to SP-RING domain, where are identified as SUMO attachment and regarded as the essential of multiple self-SUMOyation sites (Figure 9).98 Phosphorylation in the extension of the SUMOylation motif was known to interplay with SUMOylation in both animals and plants.100,

101

Thus, up-regulation of phosphorylation seems to impact on

self-SUMOylation of BpSIZ1 and probably caused the expression changes of BpDREB1 (Figure S-2). 30

ACS Paragon Plus Environment

Page 31 of 67

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

Additionally, the Ser 446 of BpICE1 (P1160) corresponding to Ser 403 of AtICE1 was significantly dephosphorylated after chilling for 6 and 48 hours. The single point mutant of Ser 403 revealed that this residue is necessary for the regulation of transactivation and stabilization.102 It is also speculated to be an O-linked glycosylation site and linked with ubiquitination. Our mass spectrometry analysis provided evidence for the phosphorylation of Ser 446 in ICE1 of paper mulberry (Figure S-4). It was possible that competition between phosphorylation and glycosylation was happened in Ser 446. The dephosphorylation of Ser 446 was enable O-linked glycosylation causing for degradation of ICE1 and further decreasing of BpDREB1 expression (Figure S-2).

CONCLUSION It is important to keep in mind that the changes observed after chilling for 48 hours could also be caused by protein abundance changes. Some details might contribute to distinguish phosphorylation or abundance changes, indirectly. If a protein was significantly expressed or suppressed, this might cause up- or down-regulation of all phosphopeptides. Some phosphoproteins such as P1959, BpEIN2, BpKCA2 with multiple phosphorylation sites contained both unchanged and changed sites, suggesting phosphorylation changes occurred during chilling tolerance. In addition, interesting phosphoproteins such as PLC292 and SIZ197 were reported to be constitutively expressed and not found to be significantly changed at the protein level 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

in the other species after chilling.103, 104 Despite this limitation, our results showed phosphorylation, an essential component, were responsible for chilling perception and defense in paper mulberry. The phosphorylation changes after chilling for 6 and 48 hours probably depended on CKIIs, MAPKs, CDKs and CDPKs and implicated in Ca2+, phytohormone, phosphoinositide and light-mediated signaling pathway for chilling perception. A hub protein (P1959) was observed to be involved in cross-talk with these signaling pathways. Additionally, transport, cytoskeleton organization and chloroplast movement were regulated against chilling at the phosphorylation level. We also highlighted the phosphorylation of BpSIZ1 and BpICE1, which might interplay with SUMOylation or glycosylation, possibly impacted on the CBF/DREB-responsive pathway and deserved for further study. Finally, we developed a schematic of chilling-related pathways in paper mulberry, which would broaden the insight into chilling response and defense mechanisms at the phosphorylation level. ASSOCIATED CONTENT Supporting Information Figure S-1: The quality control of phosphproteomics analysis. Figure S-2: The expression pattern of BpDREB1 after chilling. Figure S-3: The homology, function domain and phosphorylation sites of P1959. Figure S-4: The phosphorylation sites and function domains of ICE, PLC2 and PHOT1. Table S-1: qRT-PCR primers for the genes corresponding to significantly changed phosphoproteins. Table S-2: All of the 32

ACS Paragon Plus Environment

Page 32 of 67

Page 33 of 67

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

phosphorylation sites identified in paper mulberry leaves. Table S-3: The significantly changed phosphorylation sites between 0 h and 6 h treatment. Table S-4: The significantly changed phosphorylation sites between 0 h and 48 h treatment. Table S-5: The significantly changed phosphorylation sites between 6 h and 48 h treatment. Table S-6: Summary for the sequencing outcomes and the results of transcriptome assembly. Table S-7: The transcriptional expression of genes corresponding to significantly changed phosphoproteins. Table S-8: Differentially expressed protein kinases and phosphatases determined by RNA-seq. Table S-9: The regulatory phosphorylation sites in response to chilling. AUTHOR INFORMATION Corresponding Author * Tel and Fax: +86-10-62836545. E-mail: [email protected]. Author Contributions Z. P. and M. L. Z. performed the experimental work and wrote the manuscript. S. H. S. and X. J. P. jointly directed this work. Notes The authors declare no competing financial interest.

ACKNOWLEDGMENTS This work was supported by Poverty Relief Project of the Chinese Academy of Sciences (KFJ-FP-24), Huimin Technology Demonstration Project of the National 33

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

Modern

Agricultural

Science

and

Technology

Page 34 of 67

Achievements

City

(Z151100001015008), Knowledge Innovation Program through the Chinese Academy of Sciences (KZCX2-YW-359-2) and the National Natural Science Foundation of China (31270653).

ABBREVIATIONS MAPK, Mitogen-activated protein kinase; CDPK, Calcium-dependent protein kinase; CDK, Cyclin-dependent protein kinase; CKII, Casein kinase II; SIZ1, E3 SUMO-protein ligase SIZ1; ICE, Transcription factor ICE1; PLC, Phosphoinositide phospholipase C; PHOT, Phototropin; CAT, Catalase; APX, L-ascorbate peroxidase; CIPK, CBL-interacting serine/threonine-protein kinase; SRK2, SNF1-related protein kinase 2; KEG, E3 ubiquitin-protein ligase KEG; BSL, Serine/threonine-protein phosphatase BSL; CHUP1, Chloroplast unusual positioning 1; KCA, Kinesin-like protein for action based chloroplast movement; PMI, Plastid movement impaired.

REFERENCE (1) Thakur, P.; Kumar, S.; Malik, J. A.; Berger, J. D.; Nayyar, H. Cold stress effects on reproductive development in grain crops: an overview. Environ. Exp. Bot. 2010, 67, 429-443. (2) Aroca, R.; Tognoni, F.; Irigoyen, J. J.; Sánchez-Díaz, M.; Pardossi, A. Different root low temperature response of two maize genotypes differing in chilling 34

ACS Paragon Plus Environment

Page 35 of 67

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

sensitivity. Plant Physiol. Biochem. 2001, 39, 1067-1073. (3) Mahajan, S.; Tuteja, N. Cold, salinity and drought stresses: an overview. Arch. Biochem. Biophys. 2005, 444, 139-158. (4) Ashraf, M.; Foolad, M. Roles of glycine betaine and proline in improving plant abiotic stress resistance. Environ. Exp. Bot. 2007, 59, 206-216. (5) Shahandashti, S. S. K.; Amiri, R. M.; Zeinali, H.; Ramezanpour, S. S. Change in membrane fatty acid compositions and cold-induced responses in chickpea. Mol. Biol. Rep. 2013, 40, 893-903. (6) Lyons, J. M., Chilling injury in plants. Annu. Rev. Plant Physiol. 1973, 24, 445-466. (7) Tähtiharju, S.; Palva, T., Antisense inhibition of protein phosphatase 2C accelerates cold acclimation in Arabidopsis thaliana. Plant J. 2001, 26, 461-470. (8) Shibasaki, K.; Uemura, M.; Tsurumi, S.; Rahman, A., Auxin response in Arabidopsis under cold stress: underlying molecular mechanisms. Plant Cell 2009, 21, 3823-3838. (9) Divi, U. K.; Krishna, P., Overexpression of the brassinosteroid biosynthetic gene AtDWF4 in Arabidopsis seeds overcomes abscisic acid-induced inhibition of germination and increases cold tolerance in transgenic seedlings. J. Plant Growth Regul. 2010, 29, 385-393. (10) Zhao, D.; Shen, L.; Fan, B.; Yu, M.; Zheng, Y.; Lv, S.; Sheng, J. Ethylene and cold participate in the regulation of LeCBF1 gene expression in postharvest tomato 35

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

fruits. FEBS Lett. 2009, 583, 3329-3334. (11) Ma, Y.; Dai, X.; Xu, Y.; Luo, W.; Zheng, X.; Zeng, D.; Pan, Y.; Lin, X.; Liu, H.; Zhang, D. COLD1 confers chilling tolerance in rice. Cell 2015, 160, 1209-1221. (12) Thomashow, M. F. Molecular basis of plant cold acclimation: insights gained from studying the CBF cold response pathway. Plant Physiol. 2010, 154, 571-577. (13) Vogel, J. T.; Zarka, D. G.; Van Buskirk, H. A.; Fowler, S. G.; Thomashow, M. F. Roles of the CBF2 and ZAT12 transcription factors in configuring the low temperature transcriptome of Arabidopsis. Plant J. 2005, 41, 195-211. (14) Barrero-Gil, J.; Salinas, J. Post-translational regulation of cold acclimation response. Plant Sci. 2013, 205, 48-54. (15) Shi, Y.; Yang, S. ABA regulation of the cold stress response in plants. In Abscisic Acid: Metabolism, Transport and Signaling, Springer: 2014; pp 337-363. (16) Qin, F.; Shinozaki, K.; Yamaguchi-Shinozaki, K. Achievements and challenges in understanding plant abiotic stress responses and tolerance. Plant Cell Physiol. 2011, 52, 1569-1582. (17) Lehti-Shiu, M. D.; Shiu, S.-H. Diversity, classification and function of the plant protein kinase superfamily. Philos. Trans. R. Soc. Lond B. Biol. Sci. 2012, 367 2619-2639. (18) Reiland, S.; Messerli, G.; Baerenfaller, K.; Gerrits, B.; Endler, A.; Grossmann, J.; Gruissem, W.; Baginsky, S. Large-scale Arabidopsis phosphoproteome profiling reveals novel chloroplast kinase substrates and phosphorylation networks. Physiol. 36

ACS Paragon Plus Environment

Page 36 of 67

Page 37 of 67

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

2009, 150, 889-903. (19) Havelund, J. F.; Thelen, J. J.; Møller, I. M., Biochemistry, proteomics, and phosphoproteomics of plant mitochondria from non-photosynthetic cells. Front Plant Sci. 2013, 4, 51. (20) Mattei, B.; Spinelli, F.; Pontiggia, D.; De Lorenzo, G. Comprehensive analysis of the membrane phosphoproteome regulated by oligogalacturonides in Arabidopsis thaliana. Front Plant Sci. 2016, 7, 1107. (21) Reiland, S.; Finazzi, G.; Endler, A.; Willig, A.; Baerenfaller, K.; Grossmann, J.; Gerrits, B.; Rutishauser, D.; Gruissem, W.; Rochaix, J.-D. Comparative phosphoproteome profiling reveals a function of the STN8 kinase in fine-tuning of cyclic electron flow (CEF). P. Natl. Acad. Sci. USA 2011, 108, 12955-12960. (22) Umezawa, T.; Sugiyama, N.; Takahashi, F.; Anderson, J. C.; Ishihama, Y.; Peck, S. C.; Shinozaki, K., Genetics and phosphoproteomics reveal a protein phosphorylation network in the abscisic acid signaling pathway in Arabidopsis thaliana. Sci. Signal. 2013, 6, rs8. (23) Howden, A. J.; Salek, M.; Miguet, L.; Pullen, M.; Thomas, B.; Knight, M. R.; Sweetlove, L. J. The phosphoproteome of Arabidopsis plants lacking the oxidative signal-inducible1 (OXI1) protein kinase. New Phytol. 2011, 190, 49-56. (24) Hu, X.; Wu, L.; Zhao, F.; Zhang, D.; Li, N.; Zhu, G.; Li, C.; Wang, W. Phosphoproteomic analysis of the response of maize leaves to drought, heat and their combination stress. Front Plant Sci. 2015, 6, 298. 37

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

(25) Yu, B.; Li, J.; Koh, J.; Dufresne, C.; Yang, N.; Qi, S.; Zhang, Y.; Ma, C.; Duong, B. V.; Chen, S. Quantitative proteomics and phosphoproteomics of sugar beet monosomic addition line M14 in response to salt stress. J. Proteomics 2016, 143, 286-297. (26) Minkoff, B. B.; Stecker, K. E.; Sussman, M. R., Rapid phosphoproteomic effects of abscisic acid (ABA) on wild-type and ABA receptor-deficient A. thaliana mutants. Mol. Cell. Proteomics 2015, 14, 1169-1182. (27) Zhang, H.; Zhou, H.; Berke, L.; Heck, A. J.; Mohammed, S.; Scheres, B.; Menke, F. L. Quantitative phosphoproteomics after auxin-stimulated lateral root induction identifies an SNX1 protein phosphorylation site required for growth. Mol. Cell. Proteomics 2013, 12, 1158-1169. (28) Khan, M.; Takasaki, H.; Komatsu, S. Comprehensive phosphoproteome analysis in rice and identification of phosphoproteins responsive to different hormones/stresses. J. Proteome Res. 2005, 4, 1592-1599. (29) Chen, J.; Lin, T.; Xu, H.; Tian, D.; Luo, Y.; Ren, C.; Yang, L.; Shi, J. Cold-induced changes of protein and phosphoprotein expression patterns from rice roots as revealed by multiplex proteomic analysis. Plant Omics 2012, 5, 194. (30) El-Khatib, R. T.; Good, A. G.; Muench, D. G., Analysis of the Arabidopsis cell suspension phosphoproteome in response to short-term low temperature and abscisic acid treatment. Physiol. Plant. 2007, 129, 687-697. (31) Chang, C. S.; Liu, H. L.; Moncada, X.; Seelenfreund, A.; Seelenfreund, D.; Chung, 38

ACS Paragon Plus Environment

Page 38 of 67

Page 39 of 67

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

K. F. A holistic picture of Austronesian migrations revealed by phylogeography of Pacific paper mulberry. P. Natl. Acad. Sci. USA 2015, 112, 13537-13542. (32) Zhang, X.; Qu, L. J.; Guo, X. Q. In Study on Structure and Property of Broussonetia papyrifera (BP) Fiber, Advanced Materials Research, 2011; Trans Tech Publ: 2011; pp 1593-1596. (33) Yan, T.; Qi Y. D.; Rong, Z.; Gui L. Y.; Wei, X. Analysis on feed nutritive value of hybrid Broussonetia papyrifera leaf. Pratacult. Sci. 2009, 26, 136-139. (34) Zhao, M.; Ma, Y.; Pan, Y. H.; Zhang, C. H.; Yuan, W. X. A hevein-like protein and a class I chitinase with antifungal activity from leaves of the paper mulberry. Biomed. Chromatogr. 2011, 25 908-912. (35) Whistler, W. A.; Elevitch, C. R. Broussonetia papyrifera (paper mulberry). In Species Profiles for Pacific Island Agroforestry, Permanent Agriculture Resources: Hawaii, 2006; pp 1-13. (36) Peng, X. J.; Wu, Q. Q.; Teng, L. H.; Tang, F.; Pi, Z.; Shen, S. H., Transcriptional regulation of the paper mulberry under cold stress as revealed by a comprehensive analysis of transcription factors. BMC Plant Biol. 2015, 15, 108. (37) Peng, X. J.; Teng, L. H.; Yan, X.; Zhao, M. L.; Shen, S. H. The cold responsive mechanism of the paper mulberry: decreased photosynthesis capacity and increased starch accumulation. BMC genomics 2015, 16, 898. (38) Yan, S. P.; Zhang, Q. Y.; Tang, Z. C.; Su, W. A.; Sun, W. N. Comparative proteomic analysis provides new insights into chilling stress responses in rice. Mol. 39

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

Cell. Proteomics 2006, 5, 484-496. (39) Pi, Z.; Stevanato, P.; Yv, L. H.; Geng, G.; Guo, X. L.; Yang, Y.; Peng, C. X.; Kong, X. S. Effects of potassium deficiency and replacement of potassium by sodium on sugar beet plants. Russ. J. Plant Physiol. 2014, 61, 224-230. (40) Shen, S.; Jing, Y.; Kuang, T. Proteomics approach to identify wound-response related proteins from rice leaf sheath. Proteomics 2003, 3, 527-535. (41) Wisniewski, J. R.; Zougman, A.; Nagaraj, N.; Mann, M. Universal sample preparation method for proteome analysis. Nat. Methods 2009, 6, 359-362. (42) Wu, J.; Shakey, Q.; Liu, W.; Schuller, A.; Follettie, M. T., Global Profiling of Phosphopeptides by Titania Affinity Enrichment. J. Proteome Res. 2007, 6, 4684-4689. (43) Zanivan, S.; Gnad, F.; Wickstrom, S. A.; Geiger, T.; Macek, B.; Cox, J.; Fassler, R.; Mann, M. Solid Tumor Proteome and Phosphoproteome Analysis by High Resolution Mass Spectrometry. J. Proteome Res. 2008, 7, 5314-5326. (44) 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. (45) Cox, J.; Mann, M. MaxQuant enables high peptide identification rates, individualized p.p.b.-range mass accuracies and proteome-wide protein quantification. Nat. Biotechnol. 2008, 26, 1367-1372. 40

ACS Paragon Plus Environment

Page 40 of 67

Page 41 of 67

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

(46) Cox, J.; Mann, M., Computational Principles of Determining and Improving Mass Precision and Accuracy for Proteome Measurements in an Orbitrap. J. Am. Soc. Mass Spectrom 2009, 20, 1477-1485. (47) Gnad, F.; Forner, F.; Zielinska, D. F.; Birney, E.; Gunawardena, J.; Mann, M. Evolutionary constraints of phosphorylation in eukaryotes, prokaryotes and mitochondria. Mol. Cell. Proteomics 2010, 9, 2642-2653. (48) Lundby, A.; Secher, A.; Lage, K.; Nordsborg, N. B.; Dmytriyev, A.; Lundby, C.; Olsen, J. V. Quantitative maps of protein phosphorylation sites across 14 different rat organs and tissues. Nat. Commun. 2012, 3, 876. (49) Karpievitch, Y. V.; Dabney, A. R.; Smith, R. D. Normalization and missing value imputation for label-free LC-MS analysis. BMC Bioinf. 2012, 13, S5. (50) Taverner, T.; Karpievitch, Y. V.; Polpitiya, A. D.; Brown, J. N.; Dabney, A. R.; Anderson, G. A.; Smith, R. D. DanteR: an extensible R-based tool for quantitative analysis of -omics data. Bioinformatics 2012, 28, 2404-2406. (51) Vizcaíno, J. A.; Deutsch, E. W.; Wang, R.; Csordas, A.; Reisinger, F.; Ríos, D.; Dianes, J. A.; Sun, Z.; Farrah, T.; Bandeira, N. ProteomeXchange provides globally coordinated proteomics data submission and dissemination. Nat. Biotechnol. 2014, 32, 223-226. (52) Peng X. J.; Teng L. H.; Wang X. M.; Wang Y. C.; Shen S. H. De Novo assembly of expressed transcripts and global transcriptomic analysis from seedlings of the Paper Mulberry (Broussonetia kazinoki × Broussonetia papyifera). PLoS One 41

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

2014, 9, e97487. (53) Du, Z.; Zhou, X.; Ling, Y.; Zhang, Z.; Su, Z. agriGO: a GO analysis toolkit for the agricultural community. Nucleic Acids Res. 2010, 38, W64-W70. (54) Colaert, N.; Helsens, K.; Martens, L.; Vandekerckhove, J.; Gevaert, K. Improved visualization of protein consensus sequences by iceLogo. Nat. Methods 2009, 6, 786. (55) Schmidt, A.; Trentini, D. B.; Spiess, S.; Fuhrmann, J.; Ammerer, G.; Mechtler, K.; Clausen, T. Quantitative phosphoproteomics reveals the role of protein arginine phosphorylation in the bacterial stress response. Mol. Cell. Proteomics 2014, 13, 537. (56) Hetherington, S. E.; He, J.; Smillie, R. M. Photoinhibition at low temperature in chilling-sensitive and-resistant plants. Plant Physiol. 1989, 90, 1609-1615. (57) Duan, M.; Feng, H. L.; Wang, L. Y.; Li, D.; Meng, Q. W. Overexpression of thylakoidal ascorbate peroxidase shows enhanced resistance to chilling stress in tomato. J. Plant physiol. 2012, 169, 867-877. (58) Prasad, T. K.; Goel, R.; Kandasamy, K.; Keerthikumar, S.; Kumar, S.; Mathivanan, S.; Telikicherla, D.; Raju, R.; Shafreen, B.; Venugopal, A. Human protein reference database-2009 update. Nucleic acids Res. 2009, 37, D767-D772. (59) Wang, X.; Bian, Y.; Cheng, K.; Gu, L. F.; Ye, M.; Zou, H.; Sun, S. S. M.; He, J. X. A large-scale protein phosphorylation analysis reveals novel phosphorylation motifs and phosphoregulatory networks in Arabidopsis. J. Proteomics 2013, 78, 42

ACS Paragon Plus Environment

Page 42 of 67

Page 43 of 67

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

486-498. (60) Roach, P., Multisite and hierarchal protein phosphorylation. J. Biol. Chem. 1991, 266, 14139-14142. (61) Lv, D. W.; Subburaj, S.; Cao, M.; Yan, X.; Li, X.; Appels, R.; Sun, D. F.; Ma, W.; Yan, Y. M. Proteome and phosphoproteome characterization reveals new response and defense mechanisms of Brachypodium distachyon leaves under salt stress. Mol. Cell. Proteomics 2014, 13, 632-652. (62) Zhang, M.; Ma, C. Y.; Lv, D. W.; Zhen, S. M.; Li, X. H.; Yan, Y. M. Comparative phosphoproteome analysis of the developing grains in bread wheat (Triticum aestivum L.) under well-watered and water-deficit conditions. J. Proteome Res. 2014, 13, 4281-4297. (63) Teets, N. M.; Denlinger, D. L. Quantitative phosphoproteomics reveals signaling mechanisms associated with rapid cold hardening in a chill-tolerant fly. J. Proteome Res. 2016, 15, 2855-2862. (64) Zhang, Y.; McCormick, S. AGCVIII kinases: at the crossroads of cellular signaling. Trends Plant Sci. 2009, 14, 689-695. (65) Ohta, M.; Guo, Y.; Halfter, U.; Zhu, J. K. A novel domain in the protein kinase SOS2 mediates interaction with the protein phosphatase 2C ABI2. P. Natl. Acad. Sci. 2003, 100, 11771-11776. (66) Yoshida, R.; Umezawa, T.; Mizoguchi, T.; Takahashi, S.; Takahashi, F.; Shinozaki, K. The regulatory domain of SRK2E/OST1/SnRK2.6 interacts with ABI1 and 43

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

integrates abscisic acid (ABA) and osmotic stress signals controlling stomatal closure in Arabidopsis. J. Biol. Chem. 2006, 281, 5310-5318. (67) Wang, H.; Jian, L. C. Elevation difference of Ca2+ levels in young leaf cells of bromegrass and induced cold-tolerant enhancement under different controlled chilling temperatures. Forestry Studies in China 2008, 10, 229-234. (68) Zou, J. J.; Li, X. D.; Ratnasekera, D.; Wang, C.; Liu, W. X.; Song, L. F.; Zhang, W. Z.; Wu, W. H. Arabidopsis CALCIUM-DEPENDENT PROTEIN KINASE8 and CATALASE3 function in abscisic acid-Mediated signaling and H2O2 homeostasis in stomatal guard cells under drought Stress. Plant Cell 2015, 27, 1445. (69) Ronzier, E.; Corratgé-Faillie, C.; Sanchez, F.; Prado, K.; Briere, C.; Leonhardt, N.; Thibaud, J. B.; Xiong, T. C. CPK13, a non-canonical CPK, specifically inhibits KAT2 and KAT1 Shaker channels and reduces stomatal opening. Plant Physiol. 2014, 166, 314-326. (70) Hashimoto, K.; Eckert, C.; Anschütz, U.; Scholz, M.; Held, K.; Waadt, R.; Reyer, A.; Hippler, M.; Becker, D.; Kudla, J. Phosphorylation of calcineurin B-like (CBL) calcium sensor proteins by their CBL-interacting protein kinases (CIPKs) is required for full activity of CBL-CIPK complexes toward their target proteins. J. Biol. Chem. 2012, 287, 7956-7968. (71) Belin, C.; De Franco, P. O.; Bourbousse, C.; Chaignepain, S.; Schmitter, J. M.; Vavasseur, A.; Giraudat, J.; Barbier-Brygoo, H.; Thomine, S. Identification of features regulating OST1 kinase activity and OST1 function in guard cells. Plant 44

ACS Paragon Plus Environment

Page 44 of 67

Page 45 of 67

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

Physiol. 2006, 141, 1316-1327. (72) Liu, H.; Stone, S. L. Abscisic acid increases Arabidopsis ABI5 transcription factor levels by promoting KEG E3 ligase self-ubiquitination and proteasomal degradation. Plant Cell 2010, 22, 2630-2641. (73) Mora-García, S.; Vert, G.; Yin, Y.; Caño-Delgado, A.; Cheong, H.; Chory, J. Nuclear protein phosphatases with Kelch-repeat domains modulate the response to brassinosteroids in Arabidopsis. Genes Dev. 2004, 18, 448-460. (74) Ju, C.; Yoon, G. M.; Shemansky, J. M.; Lin, D. Y.; Ying, Z. I.; Chang, J.; Garrett, W. M.; Kessenbrock, M.; Groth, G.; Tucker, M. L. CTR1 phosphorylates the central regulator EIN2 to control ethylene hormone signaling from the ER membrane to the nucleus in Arabidopsis. P. Natl. Acad. Sci. 2012, 109, 19486-19491. (75) Sakai T.; Kagawa T.; Kasahara K.; Swartz T. E.; Christie J. M.; Briggs W. R.; Wada M, W.; Okada K. Arabidopsis nph1 and npl1: Blue light receptors that mediate both phototropism and chloroplast relocation. P. Natl. Acad. Sci. 2001, 98, 6969-6974. (76) Jarillo, J. A.; Gabrys, H.; Capel, J.; Alonso, J. M.; Ecker, J. R.; Cashmore, A. R. Phototropin-related NPL1 controls chloroplast relocation induced by blue light. Nature 2001, 410, 952-954. (77) Davis, P. A.; Hangarter, R. P. Chloroplast movement provides photoprotection to plants by redistributing PSII damage within leaves. Photosynth. Res. 2012, 112, 45

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

153-161. (78) Kodama, Y.; Tsuboi, H.; Kagawa, T.; Wada, M. Low temperature-induced chloroplast relocation mediated by a blue light receptor, phototropin 2, in fern gametophytes. J. Plant Res. 2008, 121, 441-448. (79) Łabuz, J.; Hermanowicz, P.; Gabryś, H. The impact of temperature on blue light induced chloroplast movements in Arabidopsis thaliana. Plant Sci. 2015, 239, 238-249. (80) Salomon, M.; Eisenreich, W.; Dürr, H.; Schleicher, E.; Knieb, E.; Massey, V.; Rüdiger, W.; Müller, F.; Bacher, A.; Richter, G. An optomechanical transducer in the blue light receptor phototropin from Avena sativa. P. Natl. Acad. Sci. USA 2001, 98, 12357-12361. (81) Boex-Fontvieille, E.; Jossier, M.; Davanture, M.; Zivy, M.; Hodges, M.; Tcherkez, G. Differential protein phosphorylation regulates chloroplast movement in response to strong light and darkness in Arabidopsis thaliana. Plant Mol. Biol. Rep. 2014, 32, 987-1001. (82) Sullivan, S.; Thomson, C. E.; Kaiserli, E.; Christie, J. M. Interaction specificity of Arabidopsis 14-3-3 proteins with phototropin receptor kinases. FEBS lett. 2009, 583, 2187-2193. (83) DeLille, J. M.; Sehnke, P. C.; Ferl, R. J. The Arabidopsis 14-3-3 family of signaling regulators. Plant Physiol. 2001, 126, 35-38. (84) Kadota, A.; Yamada, N.; Suetsugu, N.; Hirose, M.; Saito, C.; Shoda, K.; Ichikawa, 46

ACS Paragon Plus Environment

Page 46 of 67

Page 47 of 67

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

S.; Kagawa, T.; Nakano, A.; Wada, M. Short actin-based mechanism for light-directed chloroplast movement in Arabidopsis. P. Natl. Acad. Sci. USA 2009, 106, 13106-13111. (85) Kasahara, M.; Kagawa, T.; Oikawa, K.; Suetsugu, N.; Miyao, M.; Wada, M. Chloroplast avoidance movement reduces photodamage in plants. Nature 2002, 420, 829-832. (86) Suetsugu, N.; Yamada, N.; Kagawa, T.; Yonekura, H.; Uyeda, T. Q.; Kadota, A.; Wada, M., Two kinesin-like proteins mediate actin-based chloroplast movement in Arabidopsis thaliana. P. Natl. Acad. Sci. USA 2010, 107, 8860-8865. (87) Kodama, Y.; Suetsugu, N.; Kong, S. G.; Wada, M. Two interacting coiled-coil proteins, WEB1 and PMI2, maintain the chloroplast photorelocation movement velocity in Arabidopsis. P. Natl. Acad. Sci. USA 2010, 107, 19591-19596. (88) Vanstraelen, M.; Acosta, J. A. T.; De Veylder, L.; Inzé, D.; Geelen, D. A plant-specific subclass of C-terminal kinesins contains a conserved a-type cyclin-dependent kinase site implicated in folding and dimerization. Plant Physiol. 2004, 135, 1417-1429. (89) Paolo, G.; Camilli, P. Phosphoinositides in cell regulation and membrane dynamics. Nature 2006, 443, 651-657. (90) Yoshida, S. Chemical and Biophysical Changes in the Plasma Membrane during Cold Acclimation of Mulberry Bark Cells (Morus bombycis Koidz. cv Goroji). Plant Physiol. 1984, 76, 257-265. 47

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

(91) Hunt, L.; Otterhag, L.; Lee, J. C.; Lasheen, T.; Hunt, J.; Seki, M.; Shinozaki, K.; Sommarin, M.; Gilmour, D. J.; Pical, C., Gene-specific expression and calcium activation of Arabidopsis thaliana phospholipase C isoforms. New Phytol. 2004, 162, 643-654. (92) Kanehara, K.; Yu, C. Y.; Cho, Y.; Cheong, W. F.; Torta, F.; Shui, G.; Wenk, M. R.; Nakamura, Y. Arabidopsis AtPLC2 is a primary phosphoinositide-specific phospholipase c in phosphoinositide metabolism and the endoplasmic reticulum stress response. PLoS Genet. 2015, 11, e1005511. (93) Singh, A.; Bhatnagar, N.; Pandey, A.; Pandey, G. K. Plant phospholipase C family: Regulation and functional role in lipid signaling. Cell calcium 2015, 58, 139-146. (94) Zhong, R.; Burk, D. H.; Nairn, C. J.; Woodjones, A.; Rd, M. W.; Ye, Z. H. Mutation of SAC1, an Arabidopsis SAC domain phosphoinositide phosphatase, causes alterations in cell morphogenesis, cell wall synthesis, and actin organization. Plant Cell 2005, 17, 1449-1466. (95) Di, P. G.; De, C. P., Phosphoinositides in cell regulation and membrane dynamics. Nature 2006, 443, 651-657. (96) Monks, D. E.; Aghoram, K.; Courtney, P. D.; DeWald, D. B.; Dewey, R. E. Hyperosmotic stress induces the rapid phosphorylation of a soybean phosphatidylinositol transfer protein homolog through activation of the protein kinases SPK1 and SPK2. Plant Cell 2001, 13, 1205-1219. (97) Miura, K.; Jin, J. B.; Lee, J.; Yoo, C. Y.; Stirm, V.; Miura, T.; Ashworth, E. N.; 48

ACS Paragon Plus Environment

Page 48 of 67

Page 49 of 67

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

Bressan, R. A.; Yun, D. J.; Hasegawa, P. M. SIZ1-mediated sumoylation of ICE1 controls CBF3/DREB1A expression and freezing tolerance in Arabidopsis. Plant Cell 2007, 19, 1403-1414. (98) Garcia-Dominguez, M.; March-Diaz, R.; Reyes, J. C. The PHD domain of plant PIAS proteins mediates sumoylation of bromodomain GTE proteins. J. Biol. Chem. 2008, 283, 21469-21477. (99) Zhang, R. F.; Guo, Y.; Li, Y. Y.; Zhou, L. J.; Hao, Y. J.; You, C. X. Functional identification of MdSIZ1 as a SUMO E3 ligase in apple. J. Plant Physiol. 2016, 198, 69-80. (100) Gareau, J. R.; Lima, C. D. The SUMO pathway: emerging mechanisms that shape specificity, conjugation and recognition. Nat. Rev. Mol. Cell Biol. 2010, 11, 861-871. (101) Khan, M.; Rozhon, W.; Unterholzner, S. J.; Chen, T.; Eremina, M.; Wurzinger, B.; Bachmair, A.; Teige, M.; Sieberer, T.; Isono, E. Interplay between phosphorylation and SUMOylation events determines CESTA protein fate in brassinosteroid signalling. Nat. Commun. 2014, 5. (102) Miura, K.; Ohta, M.; Nakazawa, M.; Ono, M.; Hasegawa, P. M. ICE1 Ser403 is necessary for protein stabilization and regulation of cold signaling and tolerance. Plant J. 2011, 67, 269-279. (103) Cheng, T.; Chen, J.; Ef, A. A.; Wang, P.; Wang, G.; Hu, X.; Shi, J. Quantitative proteomics analysis reveals that S-nitrosoglutathione reductase (GSNOR) and 49

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

nitric oxide signaling enhance poplar defense against chilling stress. Planta 2015, 242, 1361-1390. (104) Neilson K. A.; Mariani M.; Haynes P. A. Quantitative proteomic analysis of cold-responsive proteins in rice. Proteomics 2011, 11, 1696-706. FIGURES LEGENDS Figure 1. The morphological and physiological changes of paper mulberry after chilling. (A) The morphology changes were photographed under a fluorescent lamp. (B) The chlorophyll fluorescence changes were captured by chlorophyll fluorometer. (C) The changes of Fv/Fm after chilling. (D) The changes of MDA during chilling. (E) The changes of REL during chilling. Error bars represent SD of three replicates at least. Different letters denote significant differences (p < 0.05). Photograph courtesy of Zhi Pi. Copyright 2016. Figure 2. Statistics analysis of phosphoproteomics of paper mulberry after chilling. (A) Venn diagram of significantly changed phosphorylation sites and phosphoproteins distributed in each sample. (B) The up-regulated and down-regulated significantly changed phosphorylation sites in every pairwise comparison. Black columns represent the number of significantly changed phosphorylation sites. Gray columns represent 16, 22 and 168 phosphorylation sites that were uniquely detected after chilling for 0, 6 and 48 hours, respectively corresponding the number shown in brackets. (C) Clustering analysis of significantly changed phosphorylation sites among three treatments. Eight groups were clustered based on k-mean clustering using MultiExperiment Viewer 50

ACS Paragon Plus Environment

Page 50 of 67

Page 51 of 67

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

(MeV). The fold change compared to unchilled treatment was log-transformed to base 2. Figure 3. IceLogo analysis of significantly changed peptides. (A) IceLogo generated after analyzing aligned sequences of each cluster using non-phosphorylated peptides for background correction. (B) IceLogo generated after analyzing aligned sequences of each cluster as compared with non-changed phosphorylation peptides. Figure 4. AgriGO analysis of the significantly changed phosphoproteins after chilling for 6 and 48 hours. (A) Biological process. (B) Molecular function. (C) Cellular component. Significantly overrepresented terms (FDR adjusted p < 0.05) are marked with an asterisk. Figure 5. The pie chart of overrepresented terms. (A) Protein modification. (B) Translation. (C) Transport. (D) Signal transduction. (E) Cellular component organization. Protein modification and translation are significantly enriched after chilling for 6 and 48 hours. Transport, signal transduction and cellular component organization were respectively overrepresented after chilling for 6 and 48 hours. Figure 6. The transcriptional expression of chilling responsive genes in paper mulberry. (A) The up-regulated and down-regulated DEGs in every comparison between control and chilling treatments. (B) The up-regulated and down-regulated protein kinases and phosphatases in every comparison between control and chilling treatments. “Percentage” referred to the percentage of the differentially expressed protein kinases and phosphatases accounting for the total differentially expressed genes of each 51

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

Page 52 of 67

pairwise comparison. (C) The number of non-differentially and differentially expressed genes corresponding to significantly changed phosphoproteins. (D) The heat map of selected

phosphoproteins

responding

to

chilling

at

the

transcription

and

phosphorylation level. Black letters mean expression of genes did not significantly change after chilling. Up and down-regulated genes are shown in red and green letters, respectively. The number in brackets represented the code number of phosphorylation sites corresponding to Table S-3 and 4. Figure 7. Protein-protein interaction analysis of protein kinases and phosphatases in response to chilling. The diameters of nodes depend on the number of first neighbors. Red and blue lines present edges of P1959 and STN8, respectively. Blue and yellow nodes mean proteins that change in response to chilling at the transcription or phosphorylation levels corresponding to Table S-8 and Table 1, respectively. Green nodes mean kinases and phosphatases responding at both the transcription and phosphorylation levels. Figure 8. Schematic presentation of systematic chilling response and defense mechanisms in paper mulberry leaves at phosphorylation level. Red and green balls represent up- and down-regulation at the phosphorylation level after chilling. Among these, the phosphorylation regulatory sites were showed in Table S-9. Solid lines represent catalytic reaction, while dotted lines represent material transportation. Red and green arrows indicate phosphorylation and dephosphorylation reactions, respectively. Yellow arrows represent ubiquitination or SUMOylation. Black arrows 52

ACS Paragon Plus Environment

Page 53 of 67

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

mean the other catalytic reaction. The area outlined with red dotted lines represents an unclear phosphorylation cascade containing hub protein P1959. Figure 9. Sequence alignment, structural composition, three-dimensional structure, and tandem mass spectrometry (MS/MS) spectra map for SIZ1 in paper mulberry. (A) BpSIZ1 sequence alignment with Morus alba, Populus trichocarpa, Arabidopsis thaliana and Glycine max. Phosphopeptide detected by MS/MS was underlined. The phosphorylation and SUMOylation sites were labeled by asterisk. (B) Function domain distribution of BpSIZ1. (C) Phosphopeptide MS/MS spectra map showed the phosphorylation sites at Ser 495 and Ser 497.

Table 1. The protein kinases and phosphatases in response to chilling at the phosphorylation level Protein Phosphosite

Description

Peptide

Cluster No.a

ID

ID

P114

77

Probable receptor kinase

_GM(ox)IPLTQIPNVS(ph)R_

3

P206

135

Probable serine/threonine protein kinase

_VVVYS(ph)DR_

3

P206

136

Probable serine/threonine protein kinase

_VVVYSDRLS(ph)SGESR_

3

P285

202

Phototropin-1

_RNS(ph)HAGTR_

1

P285

203

Phototropin-1

_RS(ph)FMGLIR_

1

P285

198

Phototropin-1

_GFPRVS(ph)EDLK_

4

P285

200

Phototropin-1

_RNSES(ph)VAPPR_

4

P285

194

Phototropin-1

_ALS(ph)ESTNRPFIR_

2

P295

211

_VESSSS(ph)AR_

3

Probable LRR receptor-like serine/threonine protein kinase

P315

229

Calcium-dependent protein kinase 13

_FNSLS(ph)VK_

1

P358

280

Probable inactive receptor protein kinase

_AATESES(ph)LSGDYSK_

4

P361

285

Serine/threonine protein kinase dst4 isoform X1

_FSS(ph)LELIGR_

2

P385

303

Phosphatase 2A 55 kDa regulatory subunit B

_GAES(ph)PGVDANGNSFD FTTK_

53

ACS Paragon Plus Environment

4

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

Page 54 of 67

P385

304

Phosphatase 2A 55 kDa regulatory subunit B

_SLS(ph)SSITR_

5

P453

377

Probable protein phosphatase 2C 6

_VNS(ph)LLSLPR_

4

P501

429

Probable protein phosphatase 2C 39

_AFGDKS(ph)LK_

5

P515

435

Serine/threonine protein kinase PBS1 isoform X1

_LGPVGDNTHVST(ph)R_

2

P515

437

Serine/threonine protein kinase PBS1 isoform X1

_(ac)SCFSCFS(ph)PGRK_

3

P530

450

CDPK-related protein kinase 4

_WPLPPPS(ph)PAK_

1

P530

449

CDPK-related protein kinase 4

_TES(ph)GIFR_

5

P544

459

Probable inactive receptor protein kinase

_APLQS(ph)PTR_

1

P544

461

Probable inactive receptor protein kinase

P544

460

Probable inactive receptor protein kinase

P546

468

Serine/threonine protein kinase BLUS1 isoform X1

P546

465

Serine/threonine protein kinase BLUS1 isoform X1

_NGS(ph)ITEATTSTSEK_

4

P546

469

Serine/threonine protein kinase BLUS1 isoform X1

_RTPS(ph)FSGPLMLPNR_

4

P570

486

Serine/threonine protein kinase dst1

_MSSS(ph)SIPESVIR_

5

Dual specificity tyrosine-phosphorylation-regulated

_VIDLGSSCFETDHLCS(ph)

kinase 4

YVQSR_

_MIEEIRQS(ph)DSENRPSSE ENK_ _AS(ph)AEVLGK_ _QTYSGPLMPGAVLSHSVS (ph)ER_

4 8 3

3

P579

493

P591

2751

Probable receptor protein kinase

_LLDSGESHITT(ph)R_

4

P598

513

Serine/threonine protein kinase EDR1

_AFS(ph)DLNCDR_

3

P628

538

Cyclin-dependent protein kinase G-2 isoform X1

_VVENGTKS(ph)PAER_

5

P664

575

Serine/threonine protein kinase ATM

_EHS(ph)PDEMLSK_

3

P667

580

Probable inactive receptor protein kinase

P682

603

Serine/threonine protein kinase fray2

_GIS(ph)AWNFNLEDLK_

2

P682

602

Serine/threonine protein kinase fray2

_FLSGT(ph)LLPDNALSAK_

3

P682

2774

Serine/threonine protein kinase fray2

_FLSGT(ph)LLPDNALSAK_

5

P682

601

Serine/threonine protein kinase fray2

_FKVTSAELS(ph)PK_

6

P710

642

_VDSGGLS(ph)PDSWR_

3

P710

641

P766

686

P847

748

P893

G-type lectin S-receptor-like serine/threonine protein kinase

_IEELCRS(ph)SLREDPQSDL VK_

G-type lectin S-receptor-like serine/threonine protein

_LLQALVSGDS(ph)FHGVK

kinase

_

Probable LRR receptor-like serine/threonine protein

3

4

_LNEEENTHIS(ph)TR_

5

Calcium-dependent kinase 32

_DGS(ph)LKLDNESR_

3

787

Leucine-rich repeat receptor protein kinase SRF6

_HKS(ph)FDDDDFSK_

3

P893

2808

Leucine-rich repeat receptor protein kinase SRF6

_T(ph)IGVDQGTSR_

5

P920

3129

Mitogen-activated protein kinase 12

P1085

986

Serine/threonine protein kinase GRIK2

_S(ph)EEILHFR_

4

P1113

1008

Serine/threonine protein kinase SRK2I

_SSVLHS(ph)QPK_

3

kinase

_VSFNDAPSAIFWTDY(ph)V ATR_

54

ACS Paragon Plus Environment

6

Page 55 of 67

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

_ISYPETGGGGGS(ph)SHGS

P1186

1075

Probable serine/threonine protein kinase

P1186

1073

Probable serine/threonine protein kinase

_ASHNNPTIPNVS(ph)K_

5

P1192

1081

Probable inactive receptor protein kinase

_DFGS(ph)GVHDAEK_

5

P1192

2862

Probable inactive receptor protein kinase

_MIEDLRPQLT(ph)K_

5

P1204

1095

Mitogen-activated protein kinase kinase kinase 13-A

_S(ph)DDSFGSQFLK_

3

P1204

1096

Mitogen-activated protein kinase kinase kinase 13-A

_NEGLGS(ph)TNQR_

4

P1234

1120

LRR receptor-like serine/threonine protein kinase ERL1 _ITENLS(ph)EK_

8

P1250

1132

CDPK-related protein kinase 5-like

_NS(ph)PALAPASADR_

2

GEAR_

3

P1257

2870

Serine/threonine protein kinase HT1-like

_VWT(ph)MDKNK_

5

P1267

1143

Ethylene response sensor 1

_FTQPLSGS(ph)GR_

6

P1300

1186

SNF1-related kinase regulatory subunit beta-1

_S(ph)NGELGAR_

5

P1326

1205

Probable protein phosphatase 2C 12

_S(ph)EHHTVPLSVLLK_

3

P1406

2894

CBL-interacting serine/threonine protein kinase 23

_EDGLLHT(ph)TCGTPNYV APEVINNK_

3

P1413

1278

Protein-tyrosine-phosphatase MKP1

_S(ph)LDEWPK_

3

P1413

1279

Protein-tyrosine-phosphatase MKP1

_S(ph)LKLPVLAEK_

3

P1429

2901

Probable serine/threonine protein kinase

P1446

1307

Calmodulin-binding receptor-like cytoplasmic kinase 1

P1455

1315

Receptor protein kinase FERONIA

P1455

1316

Receptor protein kinase FERONIA

P1558

2928

Serine/threonine protein kinase STN8

P1568

1422

P1586

1443

P1617

1491

Serine/threonine protein kinase tricorner

_IADFGLASFFDPHQNQPLT (ph)SR_ _RAAEQDIGATHIS(ph)T(ph) QVK_ _GKDPDAS(ph)PGFDGTVT DSR_ _HFS(ph)FSEIK_ _AAPETCAEFLGSFVADKT NT(ph)QFTK_

Serine/threonine protein kinase minibrain-like isoform _VDHHPGGGHWFAAGLS(p X1

h)PNIPGR_

Leucine-rich repeat receptor-like serine/threonine protein kinase

_AS(ph)AETLGR_ _DFSVGNNLS(ph)GALQSD GR_

3

3

3 5 4

6

1

5

P1650

1513

Probable protein phosphatase 2C 52

_VNS(ph)LVQLPR_

3

P1674

1532

Serine/threonine protein kinase, isoform 1

_ASCFAGRPS(ph)ER_

3

P1674

1537

Serine/threonine protein kinase, isoform 1

P1674

1533

Serine/threonine protein kinase, isoform 1

_LSFLPELIAS(ph)AK_

4

P1689

1546

Casein kinase I isoform delta-like

_PVVNPGPS(ph)AER_

7

_VAENTDAENIGS(ph)MTK _

3

P1694

1550

Mitogen-activated protein kinase kinase kinase YODA _AENPTS(ph)PGSR_

3

P1730

1598

Probable receptor protein kinase

_AES(ph)APELR_

3

P1780

1633

Probable LRR receptor-like serine/threonine protein

_LGQSLS(ph)VQK_

3

55

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

Page 56 of 67

kinase P1824

1682

Serine/threonine-protein kinase EDR1

P1826

1683

Protein kinase 2B

P1855

1721

Probable inactive receptor protein kinase

_LHKGDMS(ph)PEK_

3

P1855

1720

Probable inactive receptor protein kinase

_DGLGIGFS(ph)GK_

5

P1858

2978

Probable protein phosphatase 2C 33

_VNT(ph)LLTLPR_

3

P1858

1723

Probable protein phosphatase 2C 33

P1879

1734

Probable inactive leucine-rich repeat receptor kinase

P1921

2983

Serine/threonine protein phosphatase 7

P1922

2984

Casein kinase I isoform alpha

P1926

1782

P1926

1783

P1944

1808

P1959

1822

P1991

Receptor-like serine/threonine protein kinase ALE2 isoform X1 Receptor-like serine/threonine protein kinase ALE2 isoform X1 Cyclin-dependent protein kinase D-3-like

_HSS(ph)PLQR_ _VDTAQSSLVTSAS(ph)GVS K_

_S(ph)PHASSSPLSPNLGVK _ _LANIAGGS(ph)PTLQSSR_ _ET(ph)SVSNENKESVPNEK _ _LAAT(ph)EDEGSTSPLPER _

3 3

4 5 4

3

_GLSGGS(ph)VSEHGFTLR_

2

_GLSGGSVS(ph)EHGFTLR_

7

_NSEFNPHEGPTVLS(ph)PP R_

1

Phosphatase 2C and cyclic nucleotide-binding kinase

_VPVPQVLEVTGSES(ph)PS

domain-containing protein

TFSWISK_

1853

Protein phosphatase 2C 66 isoform 1

_ALLECAS(ph)SR_

5

P2002

1863

Probable serine/threonine protein kinase

_DIS(ph)VVSATR_

7

P2015

3000

Serine/threonine protein kinase STY46 ame

P2045

1905

Serine/threonine-protein kinase SRK2E

_SALLHS(ph)QPK_

4

P2065

1915

Probable receptor protein kinase

_FHSTS(ph)LVVK_

1

P2101

3011

Probable serine/threonine protein kinase

P2101

3012

Probable serine/threonine protein kinase

_HPWPT(ph)FQTSPVYETD QTR_

_APMQT(ph)WSGPLVDPAA VGAPR_ _IADFGLASFFDPNNKQPM T(ph)SR_

3

3

3

3

P2110

1972

Phytochrome A

_VSNS(ph)VSADQQPR_

4

P2128

1987

Receptor-like serine/threonine protein kinase ALE2

_EATEGS(ph)QHISTR_

4

P2167

2029

Serine/threonine protein kinase EDR1 isoform X1

_MLSDDLKRPS(ph)NVEK_

3

P2167

2027

Serine/threonine protein kinase EDR1 isoform X1

_HST(ph)FLSSR_

4

P2178

2049

Serine/threonine protein phosphatase BSL1

_QLS(ph)LDQFENESR_

4

P2221

2082

Serine/threonine protein kinase STY17

_PLMS(ph)PNPMHVR_

7

P2296

2149

Protein-tyrosine-phosphatase MKP1

P2296

2150

Protein-tyrosine-phosphatase MKP1

_KEEAGS(ph)PGRDEVATNI GER_ _VHAVPAS(ph)PNSMLR_

56

ACS Paragon Plus Environment

4 5

Page 57 of 67

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

_MIEEIRPSDSESRPS(ph)PE

P2309

2165

Probable inactive receptor protein kinase

P2400

2258

Phototropin-2

_SQS(ph)HDILSK_

4

P2400

2257

Phototropin-2

_NRLS(ph)ENTEQQSAK_

4

P2410

2270

Proline-rich receptor protein kinase PERK1

_FSFDTDTHVST(ph)R_

4

P2424

2284

Calmodulin-binding receptor-like cytoplasmic kinase 3 _YSMS(ph)PQLSR_

3

P2424

2282

Calmodulin-binding receptor-like cytoplasmic kinase 3 _LGS(ph)VHLNMSQVAK_

5

P2424

2283

Calmodulin-binding receptor-like cytoplasmic kinase 3 _VTAS(ph)PFR_

7

P2433

2296

Probable leucine-rich repeat receptor protein kinase

_AAS(ph)RGVGIPQLK_

3

P2443

2300

Serine/threonine protein kinase PBS1

_ILKNEEGGGS(ph)GR_

3

P2443

2301

Serine/threonine protein kinase PBS1

_LGPTGDKS(ph)HVSTR_

3

P2507

2373

Probable serine/threonine protein kinase

P2536

3082

Cyclin-dependent protein kinase C-1

_SFSNDHDANLT(ph)NR_

4

P2536

2394

Cyclin-dependent protein kinase C-1

_EIVTS(ph)PGPEK_

5

P2565

2412

Receptor-like protein kinase TMK3

_LAPEGKAS(ph)FETR_

3

_FTSGLGFRLS(ph)PK_

3

Phosphatidylinositol 3,4,5-trisphosphate 3-phosphatase

DNK_

_EHPPS(ph)PNTVTENWPSS R_

3

3

P2605

2448

P2676

2518

E3 ubiquitin-protein ligase KEG

_S(ph)PPASPDNDLTK_

3

P2681

2526

Serine/threonine protein kinase ATG1c

_KVEGFPLSACS(ph)PAR_

5

P2694

2547

Serine/threonine protein phosphatase BSL3

_SVVTS(ph)PR_

3

and -tyrosine-phosphatase PTEN2A

P2694

2544

Serine/threonine protein phosphatase BSL3

P2694

3105

Serine/threonine protein phosphatase BSL3

P2699

2554

Serine/threonine protein kinase CTR1

_LILFGGATALEGNSAASG TPS(ph)SAGSAGIR_ _ITPFGEPPT(ph)PR_ _S(ph)AAGTPEWMAPEVLR _

4 6 4

P2723

3110

Probable receptor protein kinase

_EHLNDPLRT(ph)PLNWR_

7

P2742

2599

Probable phosphatase 2C 4

_SFS(ph)QGGFAFR_

1

P2742

2598

Probable phosphatase 2C 4

_GFMS(ph)GPLDR_

3

a The tendencies of phosphorylation change corresponding to Figure 2C

57

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

For TOC only

Photograph courtesy of Zhi Pi. Copyright 2016.

58

ACS Paragon Plus Environment

Page 58 of 67

Page 59 of 67

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 1. The morphological and physiological changes of paper mulberry after chilling. (A) The morphology changes were photographed under a fluorescent lamp. (B) The chlorophyll fluorescence changes were captured by chlorophyll fluorometer. (C) The changes of Fv/Fm after chilling. (D) The changes of MDA during chilling. (E) The changes of REL during chilling. Error bars represent SD of three replicates at least. Different letters denote significant differences (p < 0.05). Photograph courtesy of Zhi Pi. Copyright 2016. Figure 1 110x115mm (300 x 300 DPI)

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

Figure 2. Statistics analysis of phosphoproteomics of paper mulberry after chilling. (A) Venn diagram of significantly changed phosphorylation sites and phosphoproteins distributed in each sample. (B) The upregulated and down-regulated significantly changed phosphorylation sites in every pairwise comparison. Black columns represent the number of significantly changed phosphorylation sites. Gray columns represent 16, 22 and 168 phosphorylation sites that were uniquely detected after chilling for 0, 6 and 48 hours, respectively corresponding the number shown in brackets. (C) Clustering analysis of significantly changed phosphorylation sites among three treatments. Eight groups were clustered based on k-mean clustering using MultiExperiment Viewer (MeV). The fold change compared to unchilled treatment was log-transformed to base 2. Figure 2 168x160mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 60 of 67

Page 61 of 67

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 3. IceLogo analysis of significantly changed peptides. (A) IceLogo generated after analyzing aligned sequences of each cluster using non-phosphorylated peptides for background correction. (B) IceLogo generated after analyzing aligned sequences of each cluster as compared with non-changed phosphorylation peptides. Figure 3 80x35mm (300 x 300 DPI)

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

Figure 4. AgriGO analysis of the significantly changed phosphoproteins after chilling for 6 and 48 hours. (A) Biological process. (B) Molecular function. (C) Cellular component. Significantly overrepresented terms (FDR adjusted p < 0.05) are marked with an asterisk. Figure 4 91x80mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 62 of 67

Page 63 of 67

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 5. The pie chart of overrepresented terms. (A) Protein modification. (B) Translation. (C) Transport. (D) Signal transduction. (E) Cellular component organization. Protein modification and translation are significantly enriched after chilling for 6 and 48 hours. Transport, signal transduction and cellular component organization were respectively overrepresented after chilling for 6 and 48 hours. Figure 5 102x61mm (300 x 300 DPI)

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

Figure 6. The transcriptional expression of chilling responsive genes in paper mulberry. (A) The up-regulated and down-regulated DEGs in every comparison between control and chilling treatments. (B) The upregulated and down-regulated protein kinases and phosphatases in every comparison between control and chilling treatments. “Percentage” referred to the percentage of the differentially expressed protein kinases and phosphatases accounting for the total differentially expressed genes of each pairwise comparison. (C) The number of non-differentially and differentially expressed genes corresponding to significantly changed phosphoproteins. (D) The heat map of selected phosphoproteins responding to chilling at the transcription and phosphorylation level. Black letters mean expression of genes did not significantly change after chilling. Up and down-regulated genes are shown in red and green letters, respectively. The number in brackets represented the code number of phosphorylation sites corresponding to Table S-3 and 4. Figure 6 115x74mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 64 of 67

Page 65 of 67

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 7. Protein-protein interaction analysis of protein kinases and phosphatases in response to chilling. The diameters of nodes depend on the number of first neighbors. Red and blue lines present edges of P1959 and STN8, respectively. Blue and yellow nodes mean proteins that change in response to chilling at the transcription or phosphorylation levels corresponding to Table S-8 and Table 1, respectively. Green nodes mean kinases and phosphatases responding at both the transcription and phosphorylation levels. Figure 7 150x127mm (300 x 300 DPI)

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

Figure 8. Schematic presentation of systematic chilling response and defense mechanisms in paper mulberry leaves at phosphorylation level. Red and green balls represent up- and down-regulation at the phosphorylation level after chilling. Among these, the phosphorylation regulatory sites were showed in Table S-9. Solid lines represent catalytic reaction, while dotted lines represent material transportation. Red and green arrows indicate phosphorylation and dephosphorylation reactions, respectively. Yellow arrows represent ubiquitination or SUMOylation. Black arrows mean the other catalytic reaction. The area outlined with red dotted lines represents an unclear phosphorylation cascade containing hub protein P1959. Figure 8 140x110mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 66 of 67

Page 67 of 67

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 9. Sequence alignment, structural composition, three-dimensional structure, and tandem mass spectrometry (MS/MS) spectra map for SIZ1 in paper mulberry. (A) BpSIZ1 sequence alignment with Morus alba, Populus trichocarpa, Arabidopsis thaliana and Glycine max. Phosphopeptide detected by MS/MS was underlined. The phosphorylation and SUMOylation sites were labeled by asterisk. (B) Function domain distribution of BpSIZ1. (C) Phosphopeptide MS/MS spectra map showed the phosphorylation sites at Ser 495 and Ser 497. Figure 9 81x37mm (300 x 300 DPI)

ACS Paragon Plus Environment