Systematic Identification and Analysis of Lysine Succinylation in

Characterization of Polysaccharides from Feijoa Fruits (Acca sellowiana Berg.) and Their Utilization as Growth Substrates by Gut Commensal Bacteroides...
0 downloads 0 Views 1MB Size
Subscriber access provided by University of Sunderland

Omics Technologies Applied to Agriculture and Food

Systematic identification and analysis of lysine succinylation in strawberry stigmata Xianping Fang, Ya Xin, Zheliang Sheng, Hui Liu, Aili Jiang, Fang Wang, Jian Yang, Xiaojun Xi, Qian Zha, Liqing Zhang, Liangying Dai, Chengqi Yan, and jianping chen J. Agric. Food Chem., Just Accepted Manuscript • DOI: 10.1021/acs.jafc.8b02708 • Publication Date (Web): 27 Aug 2018 Downloaded from http://pubs.acs.org on August 28, 2018

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 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 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.

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 36

Journal of Agricultural and Food Chemistry

1

Systematic identification and analysis of lysine succinylation in

2

strawberry stigmata

3

Xianping Fang1, 2, Ya Xin3, Zheliang Sheng4, Hui Liu3, Aili Jiang2, Fang Wang5, Jian Yang1,

4

Xiaojun Xi2, Qian Zha2, Liqing Zhang2, Liangying Dai6, Chengqi Yan5* and Jianping

5

Chen1*

6

1

Institute of Plant Virology, Ningbo University, Ningbo 315211, China.

7

2

Institute of Forestry and Pomology, Shanghai Academy of Agricultural Sciences, Shanghai

8

201403, China.

9

3

Hangzhou Academy of Agricultural Sciences, Hangzhou 310024, China

10

4

Institute of Mountain Hazards and Environment, Chinese Academy of Sciences, Chengdu

11 12

610041, China. 5

13 14

Institute of Biotechnology, Ningbo Academy of Agricultural Sciences, Ningbo 315040, China.

6

College of Plant Protection, Hunan Agricultural University, Changsha 410128, China.

15 16

* Correspondence:

17

Jianping Chen, 86-574- 87609779, [email protected];

18

Chengqi Yan, 86-574-89184038, [email protected];

19 20 21 22 23 24 25 26 27 1

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 2 of 36

28

Abstract: The various post-translational modifications of plant proteins have important

29

regulatory roles in development. We therefore examined various modified proteins from

30

strawberry stigmata and found that succinylation of lysine residues was the most abundant

31

type of modification. We then subjected proteins from strawberry stigmata to an efficient

32

enrichment method for succinylated peptides and identified 200 uniquely succinylated

33

lysines in 116 proteins. A bioinformatics analysis revealed that these proteins are involved

34

in important biological processes, including stress responses, vesicular transport, and

35

energy metabolism. Proteomics, combined with immunoprecipitation and immunoblotting,

36

revealed an obvious increase in succinylation of the assembly polypeptide 2 (AP2) and

37

clathrin from 0.5 to 2 h after pollination, suggesting that succinylation is involved in the

38

recognition of pollen-stigma signaling substances and vesicular transport. These results

39

suggest that AP2/clathrin-mediated vesicular transport processes are regulated by lysine

40

succinylation during pollen recognition.

41

Keywords: strawberry stigmata, succinylation, proteomics, AP2, clathrin

42 43 44 45 46 47 48 49 50 51 52 2

ACS Paragon Plus Environment

Page 3 of 36

Journal of Agricultural and Food Chemistry

53

Introduction

54

The successful completion of sexual reproduction in flowering plants relies on

55

pollen-stigma interactions. The stigma of angiosperms screens pollens that falls onto it, and

56

will usually recognize compatible pollen. Subsequently, the pollen absorbs water from the

57

stigma and this hydration increases its metabolic activity. Although sometimes the pollen

58

tube growth starts even in case of the incompatible or foreign pollen, the growth can be

59

stopped later. Thus the screening process is the first step of recognition between pollen and

60

pistil, and consists of a mutual exchange of signal substances between the two.1 The stigma

61

plays the most important role in this process of exchange and recognition,2 which is mostly

62

mediated by interactions between proteins.3 Some proteases in the stigma may be involved

63

in processes of signal transduction, signal recognition and stress resistance,3 but their exact

64

biological function is still not fully understood.4

65

Flower stigmata can effectively promote pollen germination by allowing pollen tubes to

66

enter their tissues, so the research work on stigmata is of great significance.5

67

Comprehensive analysis on spatiotemporal expression profiles in stigmata on

68

transcriptomic levels have been performed in Arabidopsis thaliana,6 Oryza sativa,7

69

Nicotiana tabacum.8 Besides, Goldman et al. used a specific promoter of the secretory zone

70

cells in tobacco stigmata to drive cytotoxin gene expression in the tobacco stigma. As a

71

result, they obtained tobacco plants without secretory zone cells. This transgenic tobacco

72

plant did not produce any secretion on the stigma surface, resulting in female sterility. 9

73

Research by Kim et al. on the lily stigmata indicated that the chemocyanin protein in lily

74

stigma extracts could guide the polarized growth of the pollen tube, and the existence of

75

stigma cysteine-rich adhesin could enhance the guidance of chemocyanin to the polarized

76

growth of the pollen tube.10

3

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 4 of 36

77

Many cellular proteins need to be modified to some extent after expression to perform

78

their respective functions.11 Protein post-translational modification (PTM) can play

79

important roles in regulating protein activity, structure and function, and their significance

80

is widely recognized.12-16 The study of PTM has become increasingly interesting but

81

complex with the identification of many PTM types, including phosphorylation,

82

ubiquitination, methylation, succinylation, sumoylation, etc.11, 17 Such studies are useful for

83

unraveling functionality and in-depth studies of the regulatory processes of proteins with

84

PTMs (as opposed to “naked” proteins) will become more and more critical.18 PTM study

85

has the potential to reveal the diversity and active states of proteins as well as the mode and

86

mechanisms of their functionality. Finally, these methods provide a basis for the study of

87

the molecular function involved in plant development process.11 There are suggestions that

88

several proteins are involved in the pollen-stigma recognition process, or assist pollen tubes

89

to grow and stretch into stigma tissues.10 Studies on protein PTMs have important

90

significance for revealing the biological function of the stigma and to identify the signal

91

pathway(s) related to pollination recognition in the stigma. Some phosphorylated proteins

92

have been identified in Oryza sativa pistil.19 Furthermore, research on Brassica napus

93

stigmata also suggested that the armadillo repeat–containing 1 protein can ubiquitinate the

94

exocyst complex Exo70A1, thus leading to the rejection of compatible pollen as the

95

self-incompatibility response.20

96

Succinylation is another important type of PTM15 and occurs frequently in bacteria (e.g.

97

Escherichia coli and Mycobacterium tuberculosis),21,

22

98

gondii)23 and mammalian cells (e.g. human and mouse).24 Succinylation of the bacterial

99

polysaccharide can also play a vital role in nodule invasion and possibly nodule

100

development.25 Studies of succinylation in plants began recently when 325 succinylation

101

modification siteswere identified in Taxusmedia,26 and 347 succinylation modification sites

protozoans (e.g. Toxoplasma

4

ACS Paragon Plus Environment

Page 5 of 36

Journal of Agricultural and Food Chemistry

102

were identified in tomato plants by liquid chromatography-tandem mass spectrometry

103

(LC-MS/MS).27 Some computational methods were also developed to identify the

104

succinylated proteins and sites.28 Furthermore, it has been reported that the vesicular

105

transport pathway has a positive effect on the directional pollen tube growth.29 The

106

assembly polypeptide 2 (AP2) complex is a key player in clathrin-mediated vesicular

107

transport process.30 Amazingly, we found that AP2 and clathrin are both highly

108

succinylated after pollination in strawberry stigmata.

109

To reveal which protein PTM is the most widespread in strawberry stigmata, we used

110

high-quality modifiable pan antibodies (broad-spectrum antibodies against each PTM) with

111

different protein modification types in western blotting (WB) studies to comprehensively

112

screen for six common protein modification types: phosphorylation (tyrosine), acetylation,

113

succinylation, ubiquitination, methylation and sumoylation. We then conducted functional

114

analyses of the most highly expressed of these modified proteins, which have provided a

115

model for the signal regulatory and molecular mechanisms underpinning the stigma-pollen

116

interaction process.

117

Materials and methods

118

Strawberry stigmata collection and protein extraction

119

The Fragaria ananassa cultivar Hongyu was obtained from the Hangzhou Academy of

120

Agricultural Sciences, Zhejiang, China. Plants were cultivated as described.31 Briefly, plants

121

were kept in a growth chamber under a 10-h light/14-h dark cycle at 30 °C under 150 mol

122

m–2 s–1 light or at 26 °C in the dark with a relative humidity of 60%. On the day before

123

flowering, we emasculated and bagged the flowers before dehiscence. For pollination, we

124

tried to brush around to bring all pollen from the stamens into the pistils in the center of the

125

flower. Fresh stigmata (300 in total) from 100 strawberry plants served as a biological 5

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 6 of 36

126

replicate with three replicates used. Stigmata (2 g in total per replicate) were pulverized in

127

liquid nitrogen using a mortar and pestle. The powder was then transferred to a 50-mL

128

centrifuge tube and solubilized in cold 10% (w/v) trichloroacetic acid/acetone containing 50

129

mM dithiothreitol (DTT), 0.1% (w/v) Protease Inhibitor Cocktail Set VI (Merck Millipore,

130

Billerica, USA), and polyvinylpolypyrrolidone powder for 2 h at –20 °C. After

131

centrifugation at 20,000  g at 4 °C for 10 min, the supernatant was discarded. The pellet

132

was washed twice with cold acetone supplemented with 50 mM DTT and 1 mM

133

phenylmethanesulfonyl fluoride (PMSF), then air dried and resuspended in 8 M urea

134

containing 2 mM EDTA, 10 mM DTT and 0.1% (w/v) Protease Inhibitor Cocktail Set VI

135

(Merck Millipore, Billerica, USA). The sample was sonicated three times on ice using a

136

Scientz high-intensity sonicator. Debris was removed by centrifugation at 20,000  g at

137

4 °C for 10 min. The supernatant was transferred into a new tube, and the protein content

138

was determined using reagents from a 2-D Quant kit (GE Healthcare). Stigmata were

139

collected at three time points (2 h, 0.5 h, and 0 h) before pollination and two time points (0.5

140

h and 2 h) after pollination. Only stigma protein at 0 h was used for the WB screening of the

141

various PTMs and the immunoprecipitation (IP) and subsequent succinylome analysis.

142

Strawberry stigmata samples collected at all five different time points were used for IP and

143

WB experiments of clathrin, AP2 and glyceraldehyde 3-phosphate dehydrogenase (GAPDH)

144

(Figure 1).

145

Protein tryptic digestion

146

Protein disulfides were chemically reduced with 10 mM DTT for 1 h at 56 °C and then

147

alkylated with 55 mM iodoacetamide for 45 min at room temperature in the dark. Next,

148

proteins were precipitated with three volumes of pre-chilled acetone for 30 min at –20 °C.

149

After centrifugation at 13,000  g at 4 °C for 20 min, the pellet was dissolved in 0.5 M

150

triethylamonium bicarbonat (TEAB) and sonicated for 5 min. Following a second 6

ACS Paragon Plus Environment

Page 7 of 36

Journal of Agricultural and Food Chemistry

151

centrifugation step as described above, the supernatant was collected. Protein in the

152

supernatant (~3 mg) was digested with trypsin (Promega) overnight at 37 °C (1:50

153

trypsin-to-protein mass ratio).

154

Affinity enrichment of succinylated peptides

155

Affinity enrichment of succinylated peptides was performed as described.32 Briefly, to

156

enrich for succinylated peptides, tryptic peptides, prepared as described above, were

157

dissolved in NETN buffer (100 mM NaCl, 1 mM EDTA, 50 mM Tris-HCl, pH 8.0, 0.5%

158

(w/v) NP-40) and incubated with a pre-washed pan-anti-succinylated lysines (SuKs)

159

antibody conjugated to agarose beads (PTM-402, PTM Biolabs, Hangzhou, China) overnight

160

at 4 °C with gentle shaking. The beads were carefully washed four times with NETN buffer

161

and twice with distilled H2O. Bound peptides were eluted from the beads in 0.1% (v/v)

162

trifluoroacetic acid. Eluted fractions were combined and vacuum-dried in a SpeedVac

163

(Thermo Fisher Scientific). The peptides were passed through C18 ZipTip resin (Millipore)

164

and then subjected to LC-MS/MS.

165

LC-MS/MS

166

Peptides (1 μg) were dissolved in solvent A (0.1% (v/v) formic acid (FA) in ddH2O) and

167

loaded directly onto a reversed-phase pre-column (Acclaim PepMap 100, C18, 3-m particle

168

size, 100 Å pore size, 75 m ×2 cm column diameter and length, Thermo Fisher Scientific).

169

Peptides were separated using a reversed-phase analytical chromatography column (Acclaim

170

PepMap RSLC, C18, 3-m particle size, 100 Å pore size, 75 m × 15 cm column diameter

171

and length, Thermo Fisher Scientific) connected to an EASY-nLC 1000 UPLC system. The

172

eluent system was a linear gradient of 5–20% solvent B (0.1% (v/v) FA in acetonitrile) over a

173

30-min period and then 20–35% solvent B for a 10-min period at a flow rate of 300 nL/min.

7

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 8 of 36

174

The resulting peptides were subjected to MS/MS using a Q Exactive hybrid

175

quadrupole-Orbitrap mass spectrometer (Thermo Fisher Scientific) as described below.

176

The peptides were initially fragmented in a nanospray ionization source coupled to a Q

177

Exactive Hybrid Quadrupole-Orbitrap mass spectrometer for MS/MS. Parent peptide ions

178

were detected in the Orbitrap at a resolution of 70,000 m/z. Peptides were selected for

179

MS/MS using a 25% normalized collisional energy with a 12% stepped-normalized

180

collisional energy, and the fragments produced from the low, medium, and high collisional

181

energies were combined and detected simultaneously.33 Ion fragments were detected in the

182

Orbitrap at a resolution of 17,500. A data-dependent procedure that alternated between one

183

mass spectral scan followed by 20 tandem mass spectral scans was applied for the 20 most

184

intense precursor ions that had a threshold ion count of 3  104 in the MS survey scan and a

185

15.0-s dynamic exclusion. The electrospray ionization voltage was 2.0 kV. An automatic

186

gain control was used to prevent overfilling the ion trap; 105 ions were accumulated for

187

generation of the MS/MS spectra. For the MS scans, the m/z range was 350 to 1600. For the

188

MS/MS scans, the smallest m/z value included was 100.

189

Database search

190

The MS/MS data were searched using MaxQuant as part of the integrated Andromeda

191

search engine (v.1.4.0.5). Tandem mass spectra were searched against the NCBI Fragaria

192

vesca

193

https://www.ncbi.nlm.nih.gov/protein/?term=txid101020%5BOrganism%5D) concatenated

194

with a reverse-decoy database and sequences of common protein contaminants. The MS raw

195

data files and corresponding search parameter files can be found on the ProteomeXchange

196

Consortium34

197

(http://proteomecentral.proteomexchange.org/cgi/GetDataset?ID=PXD008475). Trypsin/P

198

was specified as the cleavage enzyme; at most, two missed cleavages, four modifications per

database

website

with

(23,319

the

dataset

sequences,

identifier

PXD008475

8

ACS Paragon Plus Environment

Page 9 of 36

Journal of Agricultural and Food Chemistry

199

peptide, and five charges were allowed. Mass error was set to 5 ppm for precursor ions and

200

0.02 Da for the fragment ions. Carbamidomethylated cysteines were specified as a fixed

201

modification. Methionine oxidation and SuK residues were specified as variable

202

modifications. The false-discovery rate thresholds for proteins, peptides, and modification

203

sites were taken as 0.01. The minimum peptide length was set to seven. We removed any of

204

the SuKs that were identified with a localization probability (modification (K) probability)

205

of < 0.75 and contaminant protein sequences.

206

Protein annotation, classification and subcellular location prediction

207

The Gene Ontology (GO) annotated proteome for the succinylated proteins was derived

208

from the UniProt-GOA database (http://www.ebi.ac.uk/GOA/). First, the IDs of the

209

succinylated proteins were converted into the corresponding UniProt IDs and then mapped

210

onto the GO IDs according to their protein IDs. When a succinylated protein could not be

211

annotated according to a GO term in the UniProt-GOA database, InterProScan software was

212

then used to attempt to annotate the GO function of the protein based on protein sequence

213

alignment with proteins of known function.35 Then the succinylated proteins were classified

214

by GO annotation according to the ‘biological process’ and ‘molecular function’ categories.

215

The Kyoto Encyclopedia of Genes and Genomes (KEGG) database was used to annotate the

216

pathways associated with the succinylated proteins. First, the KEGG online service

217

‘Automatic Annotation Server’ was used to annotate the KEGG database description of each

218

protein.36 The annotation was then mapped onto the KEGG pathway database using the

219

KEGG online service tool KEGG mapper. WoLF PSORT was used for predicting subcellular

220

localization.37

221

Functional enrichment analysis

9

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 10 of 36

222

As noted above, GO function enrichment analyses of the three ontologies, namely

223

‘biological process’, ‘cellular component’, and ‘molecular function’, and the KEGG pathway

224

enrichment analyses were performed to gain further insight into the function and pathways of

225

the

226

bioinformation-analysis-tools

227

(https://github.com/JonesMacke/Bioinformation-anlaysis-tools.git)

228

operating systems and Perl language. All of protein annotation information from NCBI

229

Fragaria vesca database was used as background for performing GO/KEGG enrichment

230

analysis. Fisher’s exact test was then used to check for enrichment or depletion (right-tailed

231

test) of specific annotation terms among the members of the resulting protein clusters.

232

Derived p-values were further adjusted to address multiple hypothesis testing by the method

233

proposed by Benjamini and Hochberg.38 Terms with adjusted p-value < 0.05 in each cluster

234

were treated as significant.

succinylated

proteins.

GO/KEGG

analysis

were

performed

using

by

the

Unix-based

235

Analysis of sequences surrounding succinylated lysines

236

Motif-x software was used to characterize the amino acid sequences surrounding each SuK

237

residue (10 residues upstream and downstream of the SuK were included) in all protein

238

sequences to develop a consensus sequence for the SuK-containing sites. Protein sequences

239

downloaded from the NCBI Fragaria vesca database were used as background database

240

parameters in Motif-x, and all other Motif-x parameters were taken as their default settings.

241

Analysis of protein-protein interaction networks

242

The protein-protein interaction (PPI) network containing the identified succinylated

243

proteins was characterized using Cytoscape software (http://www.cytoscape.org/).39 The

244

proteins involved in the interaction network were obtained from the STRING database

245

(https://string-db.org/).40 The species Arabidopsis thaliana of STRING database was used to 10

ACS Paragon Plus Environment

Page 11 of 36

Journal of Agricultural and Food Chemistry

246

perform PPI analysis by protein sequences alignment. STRING defines a metric called the

247

“confidence score” to define the confidence that an interaction is real. We included all

248

interactions with a confidence score ≥ 0.7 (high confidence) to develop the network, which

249

was then was visualized using Cytoscape. A recently developed graphical, theoretical

250

clustering algorithm, denoted "Molecular Complex Detection" (MCODE),41 can be used to

251

identify densely connected regions in large PPI networks that may represent molecular

252

complexes, with MCODE being part of the plug-in tool kit for the network analysis and

253

visualization software Cytoscape. The nodal density is an important parameter that can be

254

used to evaluate the importance of a node in a particular network, and we therefore

255

calculated this parameter for each SuK-containing protein.

256

Immunoprecipitation, western blotting screening and validation

257

WB screening analysis was performed to reveal which protein PTM is the most

258

widespread in strawberry stigmata. The levels of AP2 and clathrin lysine succinylation were

259

also detected by combination of IP with WB at three time points (2 h, 0.5 h, and 0 h) before

260

pollination and two time points (0.5 h and 2 h) after pollination.

261

Antibodies used are listed as follows: anti-phosphotyrosine antibody (PTM-702, PTM

262

Biolabs, Hangzhou, China), anti-acetyllysine antibody (PTM-105, PTM Biolabs, Hangzhou,

263

China), anti-succinyllysine antibody (PTM-401, PTM Biolabs, Hangzhou, China),

264

anti-diglycine lysine antibody (PTM-1101, PTM Biolabs, Hangzhou, China), anti-mono-,

265

dimethyllysine antibody (PTM-602, PTM Biolabs, Hangzhou, China), anti-SUMO-1

266

(ab5316, Abcam, Cambridge, England) and anti-SUMO-3 (ab5317, Abcam, Cambridge,

267

England). Antibodies for detecting AP2, clathrin and GAPDH (reference) were raised in

268

rabbits

269

126-CKRGKRGLMEDRYSA-140,

270

Biotechnology, Hangzhou, China). Briefly, 0.5 mg of peptide fragments were emulsified

using

the

respective

peptide and

fragments

72-CMQTENLELKKLVYL-86,

18-CHRSQASCVGLQHSS-31

(HuaAn

11

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 12 of 36

271

with an equal volume of Freund's complete adjuvant (#F5881, Sigma-Aldrich) with a

272

sonicator and administered subcutaneously to the back of the healthy New Zealand white

273

rabbits at 4 sites (200 μL/site). Prior to the first immunization, about 5 mL of blood was

274

collected from the ear vein of each rabbit to obtain pre-immune serum, which was used for

275

corresponding control experiments. Three booster injections were administered 3, 5 and 7

276

weeks later using 0.25 mg of peptide fragments but emulsified with an equal volume of

277

Freund's incomplete adjuvant (#F5506, Sigma-Aldrich). Seven days after the final booster

278

injection, the blood was collected from the marginal vein in the ear and the serum was

279

separated for subsequent experiments. The pre-immune and post-immune serums were

280

assayed via enzyme-linked immunosorbent assays (ELISA), and the high titer of serum were

281

then purified. The immunoreactivity and immunogenicity of purified serum was further

282

validated by ELISA and WB. The entire process of antibody development mainly includes

283

steps of immunogen preparation, immunization, serum collection, ELISA screening,

284

purification, further evaluation and validation (Supplementary Table S1).

285

IP analysis was performed using an IP kit (#11719394001, Roche) following the

286

manufacturer’s instructions. Briefly, 400 µg of protein from total cell lysate was incubated

287

with 2 µg of antibody and subsequent to this incubation 20 µg of protein-G agarose was

288

added. The immune complex was incubated at 4 °C for 1 h,

289

centrifugation at 12,500 g for 1 min and washed four times with the buffer (20 mM Hepes /

290

NaOH (pH 7.4), 150 mM NaCl, 1 mM Na3VO4, 10 mM EDTA, 0.02 % NP-40, 10 µg/mL

291

aprotinin10 µg/mL leupeptin and 1 mM PMSF). Protein solutions were subjected by

292

SDS-PAGE, and then probed in a PVDF membrane using AP2, clathrin, GAPDH or pan

293

PTMs antibodies in a 1: 1,000 ~ 1,500 dilution. The HRP-labelled goat anti-rabbit secondary

294

antibody was diluted to 1:2500, andthe bands were probed by an ECL kit (Multisciences

295

Biotech). The immunoblot bands were quantified by band density using an ImageQuant LAS

then recovered by

12

ACS Paragon Plus Environment

Page 13 of 36

Journal of Agricultural and Food Chemistry

296

4000 with ImageQuant TL 7.0 software (GE Healthcare Life Sciences, Piscataway, NJ).

297

Each experiment was repeated two (soly WB screening) or three times (IP/WB).

298

Results

299

Identification of lysine succinylation in strawberry stigmata

300

Although WB analysis using various PTM-specific pan antibodies showed that

301

succinylation is the most evident type of modification among various PTMs

302

(phosphorylation (Tyr), acetylation, succinylation, ubiquitination, methylation and

303

sumoylation) in strawberry stigmata (Figure 2A), phosphorylation (Ser/Thr) level is not

304

determined in this study and it will be studied in our further research. Thus we cannot

305

ignore the possibility that phosphorylation (Ser/Thr) might be more evident than

306

succinylation. However, high abundance of succinylated level in strawberry stigma proteins

307

is an interesting phenomenon. Two hundred succinylation sites in 116 proteins were

308

successfully identified with high confidence (Figure 2B, Supplementary Table S2). The MS

309

proteomics dataset has been deposited into the ProteomeXchange Consortium34 via the

310

iProX partner repository under the dataset identifier PXD008475/IPX0001091000

311

(http://proteomecentral.proteomexchange.org/cgi/GetDataset?ID=PXD008475). The MS

312

data were validated by determining the mass error of all identified peptides. The

313

distribution of mass errors was near zero, and most values deviated < 0.02 Da, indicating

314

that the mass accuracy of the MS data met the requirement of proteomics analysis (Figure

315

2C). The identified succinylated proteins involved in many aspects of plant tissues and

316

organs (Figure 2D),26, 27, 42, 43 suggesting that succinylation might have important regulatory

317

roles in different plant species.

318

Comparison between lysine succinylation and acetylation in strawberry plants

13

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 14 of 36

319

Of various PTMs, acetylation has been extensively mapped in diverse organisms. We

320

compared the identified succinylation sites in strawberry stigmata of this study to our

321

previously reported acetylation sites in strawberry leaves.44 Interestingly, we found that a

322

number of 53 strawberry stigmata proteins (Figure 3A) harboring 66 succinylated sites

323

(Figure 3B) could be acetylated at the same site(s) in strawberry leaves. Two hundred

324

succinylated sites were mapped on 200 peptides and 116 proteins (Figure 3C). Among the

325

116 succinylated proteins, 83 (72%) had one succinylation site and 13 (11%) had two

326

succinylation sites (Figure 3D). All succinylated peptides had lengths of 7–26 amino acids

327

(Figure 3E). Similarly, 473 (69%) AcK proteins had one acetylated site (Figure 3D). The

328

identified acetylated peptides ranged from 7 to 29 amino acids in length (Figure 3E).

329

Functional characterization of lysine succinylated proteins

330

The transport related terms including vesicular transport, hydrogen transport and anion

331

transport, were obviously enriched (Figure 4A). In the cellular component category, the

332

cytoplasm and protein complex were obviously enriched (Figure 4A). KEGG pathway

333

enrichment analysis identified many enriched pathways related to energy metabolism,

334

including carbon metabolism, pyruvate metabolism, glycolysis/gluconeogenesis and the

335

citrate cycle (Supplementary Table S3).

336

Through GO function classification analysis, we found that biological process category,

337

stress response and transport related proteins were the major succinylated proteins

338

accounting respectively for 21%, 19%, 15% and 12% of all the lysine succinylated proteins

339

(Supplementary Figure S1A). In the molecular function analysis, proteins associated with

340

catalytic activity, binding, oxidoreductase and transporter activity accounted for nearly 60%

341

(21%, 15%, 13% and 10%, respectively) of all the identified succinylated proteins

342

(Supplementary Figure S1B). In the cellular component category, the proportions of

343

membrane, organelle and macromolecular complex related proteins were 34%, 31% and 14

ACS Paragon Plus Environment

Page 15 of 36

Journal of Agricultural and Food Chemistry

344

24%, respectively (Supplementary Figure S1C). The subcellular localization analysis of the

345

identified proteins suggested that most were predicted to be localized in the cytosol (43%)

346

(Supplementary Figure S1D), suggesting important roles for lysine succinylation in cytosol.

347

Moreover, many proteins were localized in the nuclei (20%).

348

We then compared the function category of the identified proteins in this study to the

349

previously reported succinylated proteins in other plants.26,

27, 42, 43

350

strawberry stigma succinylated proteins have a much higher percentage of transport process

351

(Figure 4B), oxidoreductase activity, transporter activity (Figure 4C), nuclear, plasma and

352

extracellular (Figure 4D) than the succinylated proteins of other four plants, indicating that

353

succinylation might be closely involved in the transport and nuclear signal transduction

354

process in strawberry stigmata.

We found that the

355

Sequence-motif analysis for the succinylated peptides

356

To identify possible specific sequence motifs proximal to the SuK residues, we examined

357

the relative frequencies of amino acids in specific positions of peptides containing 21

358

residues and centered around each SuK residue and then compared those positions with

359

those of 21-mers that did not contain an SuK. We identified four motifs that were

360

significantly enriched in SuK residues from 200 unique SuK-containing sites, which

361

accounted for 95% of the identified sites. The four consensus sequence motifs were FSuK,

362

ISuK, LSuK, and VSuK (Table 1 and Supplementary Figure S2).

363

Analysis of interaction networks of succinylated proteins in strawberry stigmata

364

To investigate how the succinylated proteins are related and how they are involved in

365

different interacting pathways, we constructed a PPI network for all identified succinylated

366

proteins using the STRING database and Cytoscape software (Figure 5). Using the

367

MCODE plug-in toolkit, we extracted several highly enriched interaction clusters from the 15

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 16 of 36

368

complete interaction network. Most of the proteins were clustered into one of five groups,

369

including the TCA cycle (succinyl-CoA ligase, malate dehydrogenase, etc.), oxidative

370

phosphorylation (ATP synthase CF1 alpha subunit, ATP synthase subunit beta, etc.) and

371

glycolysis (dihydrolipoyl dehydrogenase 1, fructose-bisphosphate aldolase, etc.) (Figure 5).

372

Detailed information concerning the proteins involved in these networks is presented in

373

Supplementary Table S4.

374

Pollination leads to increased hypersuccinylation of AP2 and clathrin

375

Among the different transport processes, the vesicular transport process has the highest

376

degree of enrichment (Figure 4A). AP2 and clathrin are two major types of vesicular

377

transport related proteins and the four peptides from these two proteins are obviously

378

succinylated in strawberry stigmata (Supplementary Figure S3), so we selected these two

379

proteins for subsequent validation experiments. To examine the levels of AP2 and clathrin

380

succinylation in vivo, we immunoprecipitated endogenous AP2 and clathrin by IP from

381

strawberry stigmata at five time points before and after pollination and examined protein

382

succinylation by WB (Figure 6). Quantitation of three independent experiments confirmed

383

that succinylation of AP2 and clathrin is significantly higher in pollinated stigmata when

384

compared with non-pollinated stigmata (Figure 6). The succinylation level of AP2

385

markedly increased (2.3-fold change at 0.5 h) when the stigmata were pollinated manually.

386

In contrast to AP2, the level of clathrin succinylation slightly increased at 0.5 h, yet with a

387

subsequent significant increase at 2 h after pollination. However, these two proteins levels

388

didn’t exhibit an obvious change after pollination, indicating that hypersuccinylation of

389

AP2 and clathrin is closely associated with strawberry stigma-pollen recognition process.

390

Discussion

16

ACS Paragon Plus Environment

Page 17 of 36

Journal of Agricultural and Food Chemistry

391

Protein succinylation modification refers to the process of covalent binding of a

392

succinyl-group to lysine residues in an enzymatic or non-enzymatic manner,45 so protein

393

succinylation can lead to a substantial change in the chemical properties of many proteins.

394

Importantly, the succinylation at the lysine residue induces charge mutations that, in

395

physiological pH, range from +1 to -1, thereby promoting the structural and functional

396

adjustment of the matrix proteins.15 At present, many studies have confirmed that protein

397

modification is ubiquitous in prokaryotes and eukaryotes, suggesting that lysine

398

succinylation regulates many important cellular metabolic processes in organisms.23, 32, 46

399

In this present study, we found highly significant succinylation amongst the modified

400

strawberry stigmata proteins and identified 200 succinylation modification sites in 116

401

proteins. These abundant succinylated proteins work as stress tolerance response proteins,

402

vesicular transport proteins, energy and metabolism related proteins, or are related to

403

synthesis and modification (Figure 4A). Upon flowering of the strawberry plant, the stigma

404

is directly exposed to the outside, and the pollen tube can enter the stigma tissue. Many

405

succinylated proteins (peroxidase, ascorbate peroxidase, catalase isozyme 1, etc.) are

406

related to adverse resistance and stress response process.31 This suggests that succinylation

407

is likely an important part of the signal pathway of the self-defense response that would

408

enhance resistance of the stigma to pathogen invasion, and actively create a good

409

pollination environment. After pollen germination, the stigma tissue exchanges signals with

410

the pollen tube, and may provide the raw materials needed by cells for the elongation of

411

pollen tubes. It appears that the abundant expression of succinylation-modified proteins

412

related to energy and metabolism in stigmata ensures a smooth process.

413

Besides, our study also showed that lysine succinylation is involved in the Calvin cycle,

414

oxidative phosphorylation and glycolysis process. Similarly, many proteins associated with

415

these metabolic pathways have also been found succinylated in other plants (succinyl-CoA 17

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 18 of 36

416

ligase, malate dehydrogenase, dihydrolipoyl dehydrogenase, etc.).26, 27, 42, 43 These findings

417

suggest that lysine succinylation has a positive effect on photosynthesis and the Calvin

418

cycle. Our PTM proteome data supports the conclusion that mature stigmata of strawberry

419

plants possess a vigorous metabolism. The stigma cells provide material needed for future

420

pollen-stigma recognition, pollen germination and pollen tube elongation via genetic

421

transcription as well as protein synthesis and modification.47

422

More importantly, compared with the distribution of the succinylated proteins in Taxus

423

media,26 Solanum lycopersicum,27 Oryza sativa42 and Brachypodium distachyon,43 the

424

higher distribution of biological process, molecular function and subcellular location in

425

strawberry stigmata is transport process, transport activity and nuclear/extracellular,

426

respectively (Figure 4B-D). Thus the succinylated proteins involving in transport response

427

and locating in nuclear or extracellular might have important biological functions in the

428

stigma. We identified 14 succinylated proteins involving in different transport processes,

429

and 10 of these proteins have transport activity. These transporter related proteins include

430

F-type H+-transporting ATPase, coatomer subunit beta-1-like, coatomer subunit beta-2-like,

431

clathrin, AP2 complex, ran-binding protein 1, ATP synthase, small ubiquitin-related

432

modifier 2-like, metal-nicotianamine transporter, etc. (Supplementary Table S2). The

433

interaction of pollens and stigmata requires the transport and transmission of signaling

434

agents. The fact that a large number of transporter related proteins are succinylated

435

indicates that succinylation might regulate the pollen-stigma recognition process.

436

Among the four highly enriched transport processes (vesicular transport, hydrogen

437

transport, regulation of anion transport and ion transport), the vesicular transport has the

438

highest enrichment degree (Figure 4A). Vesicular transport is not only a basic process of

439

life, but also an extremely complicated dynamic biological process, related to various

440

proteins and regulatory molecules, and employing processes that are highly conserved in 18

ACS Paragon Plus Environment

Page 19 of 36

Journal of Agricultural and Food Chemistry

441

higher eukaryotes. In addition, vesicular transport pathways are important for directional

442

pollen tube growth process.29 Studies on self-incompatible plants also indicated that

443

endocytosis,

444

pollination-recognition process.20,

445

clathrin-mediated vesicular transport process.30

446

exocytosis

and

vesicular 48-50

transport

are

all

important

for

the

The AP2 complex is a key player in

Thus, based on the previous reports, our proteomics study and the up-regulated

447

succinylation

448

AP2/clathrin-mediated vesicular transport process plays a vital role in the pollination

449

recognition process of strawberries. Molecules promoting the germination and polarized

450

growth of compatible pollen as well as those that suppress the germination and growth of

451

incompatible pollen are transferred to the cell surface via exocytosis. Simultaneously,

452

signal recognition molecules from the pollen surface are introduced into the endomembrane

453

system via endocytosis, completing the recognition process in the stigma outer layers cells.

454

We identified numerous proteins that are involved in the vesicular transport in strawberry

455

stigmata, including two kinds of the most essential vesicular coat proteins: coatomer and

456

clathrin (Figure 7 and Supplementary Table S2). Vesicular transport mediated by coatomer

457

is considered to mediate non-selective material transport. This includes transport from

458

endoplasmic reticulum to the Golgi apparatus, between Golgi vesicles, and from the Golgi

459

apparatus to plasma membranes. However, clathrin can not only mediate transport from the

460

Golgi apparatus to the lysosome, plant vacuole or to the outside of the plasma membrane,

461

but also transport foreign substances into the cytoplasm, or from the endosome to the

462

lysosome. The vesicular coat of clathrin is the grid structure constructed by clathrin fibers,

463

and many different adapter proteins fill the frame between the clathrin structure and

464

envelope. These proteins mediate the connection between clathrin and the envelope

465

transmembrane protein receptor (the transport carrier), and form the clathrin

level

of

AP2/clathrin

after

pollination,

we

speculate

that

19

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 20 of 36

466

vesicle-integrated structure together with dynamin. It is well known that four adapter

467

proteins, AP1, AP2, AP3 and AP4, selectively combine with different receptor-transporter

468

complexes, enveloping the foreign substances and forming specific transport vesicles.

469

In our study, we identified high-level succinylation in clathrin and in the adapter protein

470

AP2 (Figure 6 and Supplementary Table S2). We speculate that the succinylation of

471

clathrin/AP2 may assist in driving the cytoplasmic membrane to form a capsule and capture

472

various transmembrane specific receptor proteins (these enable the vesicles to envelop

473

specific molecules) in the vesicular transport system. A combination of plasma membrane

474

receptor and specific ligands (such as pollen proteins or other informational molecules

475

secreted from the pollen cells), enables the stigma outer layers cells memberanes to be

476

selective in AP2/clathrin via succinylation. AP2/clathrin in un-succinylated state will not

477

combine with the receptor, and succinylated AP2/clathrin can be further combined with the

478

receptor-ligand complex to form the vesicles. Then the vesicles can be budded, released and

479

fused into the organelles membranes. Thus the signal substances in pollen cells are

480

successfully transported into the stigma cell to further complete the pollination recognition

481

process (Figure 7). So, through a mechanism of interaction and recognition between pollen,

482

informational molecules and strawberry stigma, the stigma is able to respond smoothly to

483

the up-stream signal and promote successful pollen germination and pollen tube elongation.

484

Succinylation as well as the process of directional vesicular transport likely start this

485

process.

486

In this work, by combining comprehensive IP/WB analysis with high-affinity enrichment

487

of succinylated peptides as well as very sensitive mass spectrometry and advanced

488

bioinformatic tools, we systematically defined the lysine succinylome in strawberry

489

stigmata. The identification of 200 SuK sites in 116 proteins expands the lysine

490

succinylome catalog, especially for non-model plants. Functional characterization of 20

ACS Paragon Plus Environment

Page 21 of 36

Journal of Agricultural and Food Chemistry

491

succinylated proteins indicated that lysine succinylation is involved in diverse biological

492

processes and are present in various cellular components having diverse functions. Clathrin

493

and AP2 are highly succinylated, suggesting that succinylation of vesicular transport

494

proteins is important for the stigma-pollen recognition process in strawberry. In addition,

495

motif analysis extracted four consensus sequence motifs. Our research reveals that

496

succinylation modification is likely to play an important regulatory role in the

497

pollen-stigma recognition process and thus provides an important direction for the study of

498

strawberry reproductive development. The major vesicular transporters AP2 and clathrin

499

are highly succinylated after pollination and are likely to regulate the transport of signaling

500

substances via succinylation modification. It would be more meaningful if we carry out the

501

quantitative succinylome analysis at all the time points, though we have made a

502

global succinylome analyses at 0 h and found out the possible candidate succinylation sites.

503

Our future research will perform succinyl-null or succinyl-mimic mutation (i.e., a

504

lysine-to-arginine change) of the corresponding sites by PCR reaction to observe the

505

functional phenotype of plants, and further determine the role of succinylation in the

506

pollen-stigma signal transduction process.

507

Abbreviations

508

PTM, post-translational modification; LC-MS/MS, liquid chromatography-tandem mass

509

spectrometry; WB, western blotting; AP2, assembly polypeptide 2; PMSF,

510

phenylmethanesulfonyl fluoride; IP, immunoprecipitation; GAPDH, glyceraldehyde

511

3-phosphate dehydrogenase; NCE, normalized collisional energy; GO, gene ontology;

512

KEGG, kyoto encyclopedia of genes and genomes; PPI, protein-protein interaction; ELISA,

513

enzyme-linked immune-sorbent assays.

514

Funding sources 21

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 22 of 36

515

This research was supported by grants from the Hi-Tech Program (‘863’ Program) of

516

China, Ministry of Science and Technology (Grant 2014A0A603-15); the Scientific and

517

Technological Project of Hangzhou (20170432B14); the Key Program of Zhejiang

518

Provincial Foundation for Natural Science (LZ14C140001, LZ16C130002, NB2016C11017);

519

and the National Key Research and Development Program of China (2016YFD0200804).

520

Conflict of interest statement

521

522

The authors declare no conflict of interest. Acknowledgments

523

We thank M. J. Adams, Minehead, UK for his help in correcting the English of the

524

manuscript.

525

Supporting information

526

The supplementary materials including four supplementary tables and three

527

supplementary figures for this article can be found online at https://pubs.acs.org/. The mass

528

spectrometry proteomics dataset has been deposited to the ProteomeXchange Consortium

529

via the iProX partner repository under the dataset identifier PXD008475/IPX0001091000

530

(http://proteomecentral.proteomexchange.org/cgi/GetDataset?ID=PXD008475).

531 532

References

533 534 535 536 537 538 539 540 541

[1]

T. Gaude, S. Mccormick, Signaling in pollen-pistil interactions, Cell Dev. Biol. 10 (1999) 139-147.

[2]

V.E. Franklin-Tong, The difficult question of sex: the mating game, Curr. Opin. Plant Biol. 5 (2002) 14-18.

[3]

M.J. Wheeler, V.E. Franklin-Tong, F.C.H. Franklin, The molecular and genetic basis of pollen-pistil interactions, New Phytol. 151 (2001) 565-584.

[4]

S.M. Mcinnis, D.C. Emery, R. Porter, R. Desikan, J.T. Hancock, S.J. Hiscock, The role of stigma peroxidases in flowering plants: insights from further characterization of a stigma-specific peroxidase SSP from Senecio squalidus Asteraceae, J. Exp. Bot. 57 (2006) 1835-1846.

[5]

M. Wolters-Arts, L.V.D. Weerd, A.C.V. Aelst, J.V.D. Weerd, H. Vanas, C. Marian, Water-conducting properties of lipids during pollenhydration, Plant Cell Environ. 25 (2002) 513-519. 22

ACS Paragon Plus Environment

Page 23 of 36

542 543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561 562 563 564 565 566 567 568 569 570 571 572 573 574 575 576 577 578 579 580 581 582 583 584

Journal of Agricultural and Food Chemistry

[6]

R. Swanson, T. Clark, D. Preuss, Expression profiling of Arabidopsis stigma tissue identifies stigma-specific genes, Sex. Plant Reprod. 18 (2005) 163-171.

[7]

M.N. Li, W.Y. Xu, W.Q. Yang, Z.S. Kong, Y.B. Xue, Genome·wide gene expression profiling reveals conserved and novel molecular functions of the stigma in rice, Plant Physiol. 144 (2007) 1797-1812.

[8]

A.C. Quiapim, M.S. Brito, L.A.S. Bernardes, I. DaSilva, I. Malavazi, H. DePaoli, Analysis of the Nicotiana tabacum stigma/style transcriptome reveals gene expression differences between wet and dry stigma species, Plant Physiol. 149 (2009) 1211-1230.

[9]

M.H.S. Goldman, R.B. Goldberg, C. Mariani, Female sterile tobacco plants are produced by stigma-specific cell ablation, Embo J. 13 (1994) 2976-2984.

[10] S. Kim, J.C. Mollet, J. Dong, K. Zhang, S.Y. Park, E.M. Lord, Chemocyanin, a small basic protein from the lily stigma, induces pollen tube chemotropism, Proc. Natl. Acad. Sci. USA. 100 (2003) 16125-16130. [11] V.O. Jesper, M. Matthias, Status of large-scale analysis of posttranslational modifications by mass spectrometry, Mol. Cell Proteomics 12 (2013) 3444-3452. [12] J.V. Olsen, B. Blagoev, F. Gnad, B. Macek, C. Kumar, P. Mortensen, Global, in vivo, and site-specific phosphorylation dynamics in signaling networks, Cell 127 (2006) 635– 648. [13] M. Schwickart, X. Huang, J.R. Lill, J. Liu, R. Ferrando, D.M. French, Deubiquitinase USP9X stabilizes MCL1 and promotes tumour cell survival, Nature 463 (2010) 103–107. [14] Z. Zha, X. Han, M.D. Smith, Y. Liu, P.M. Giguere, D. Kopanja, A non-canonical function of Gβ as a subunit of E3 ligase in targeting GRK2 ubiquitylation. Mol. Cell 58 (2015) 794-803. [15] Z. Zhang, M. Tan, Z. Xie, L. Dai, Y. Chen, Y. Zhao, Identification of lysine succinylation as a new post-translational modification, Nat. Chem. Biol. 7 (2011) 58–63. [16] S. Zhao, W. Xu, W. Jiang, W. Yu, Y. Lin, T. Zhang, Regulation of cellular metabolism by protein lysine acetylation, Science 327 (2010) 1000-1004. [17] V. Schwämmle, S. Sidoli, C. Ruminowicz, X. Wu, C.F. Lee, K. Helin, Systems level analysis of histone H3 post-translational modifications (PTMs) reveals features of PTM crosstalk in chromatin regulation, Mol. Cell Proteomics 15 (2016) 2715-2729. [18] M. Matthias, N.J. Ole, Proteomic analysis of post-translational modifications, Nat. Biotechnol. 21 (2013) 255-261. [19] K. Wang, Y. Zhao, M. Li, F. Gao, M.K. Yang, P.F. Yang, Analysis of phosphoproteome in rice pistil, Proteomics 14 (2014) 2319-2334. [20] M.A. Samuel, Y.T. Chong, K.E. Haasen, M.G. Aldea-Brydges, D.L. Stone, D.R. Goring, Cellular pathways regulating responses to compatible andself-incompatible pollen in Brassica and Arabidopsis stigmas intersect at Exo70A1, a putative component of the exocyst complex, Plant Cell 21 (2009) 2655–2671. [21] G. Colak, Z. Xie, A.Y. Zhu, L. Dai, Z. Lu, Y. Zhang, Identification of lysine succinylation substrates and the succinylation regulatory enzyme CobB in Escherichia coli, Mol. Cell Proteomics 12 (2013) 3509–3520. [22] L. Xie, W. Liu, Q. Li, S. Chen, M. Xu, Q. Huang, The first succinyl-proteome profiling of extensively drug-resistant Mycobacterium tuberculosis revealed involvement of succinylation in cellular physiology. J. Proteome Res. 14 (2015) 107–119. [23] X. Li, X. Hu, Y. Wan, G. Xie, X. Li, D. Chen, Systematic identification of the lysine succinylation in the protozoan parasite Toxoplasma gondii, J. Proteome Res. 13 (2014) 6087–6095. [24] J. Du, Y. Zhou, X. Su, J.J. Yu, S. Khan, H. Jiang, Sirt5 is a NAD-dependent protein lysine demalonylase and desuccinylase, Science 334 (2011) 806–809. [25] J.A. Leigh, J.W. Reed, J.F. Hanks, A.M. Hirsch, G.C. Walker, Rhizobium meliloti mutants that fail to succinylate their calcofluor-binding exopolysaccharide are defective in nodule invasion, Cell 51 (1987) 579-587.

23

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

585 586 587 588 589 590 591 592 593 594 595 596 597 598 599 600 601 602 603 604 605 606 607 608 609 610 611 612 613 614 615 616 617 618 619 620 621 622 623 624 625 626 627 628

Page 24 of 36

[26] C. Shen, J. Xue, T. Sun, H. Guo, L. Zhang, Y. Meng, Succinyl-proteome profiling of a high taxol containing hybrid Taxus species (Taxus × media) revealed involvement of succinylation in multiple metabolic pathways, Sci. Rep. 6 (2016) 21764. [27] W. Jin, F. Wu, Proteome-wide Identification of lysine succinylation in the proteins of tomato (Solanumly copersicum), PLoS ONE 11(2016) e0147586. [28] Y. Xu, Y.X. Ding, J. Ding, Y.H. Lei, L.Y. Wu, N.Y. Deng, iSuc-PseAAC: predicting lysine succinylation in proteins by incorporating peptide position-specific propensity, Sci. Rep. 5 (2015) 10184. [29] A. Krichevsky, S.V. Kozlovsky, G.W. Tian, M.H. Chen, A. Zaltsman, V. Citovsky, How pollen tubes grow, Dev Biol. 303 (2007) 405-420. [30] B.T. Kelly, S.C. Graham, N. Liska, P.N. Dannhauser, S. Höning, E.J. Ungewickell, AP2 controls clathrin polymerization with a membrane-activated switch, Science 345 (2014) 459–463. [31] X.P. Fang, W.Y. Chen, Y. Xin, H.M. Zhang, C.Q. Yan, H. Liu, Proteomic analysis of strawberry leaves infected with Colletotrichum fragariae, J. Proteomics 75 (2012) 4074-4090. [32] M. Yang, Y. Wang, Y. Chen, Z. Cheng, J. Gu, J. Deng, Succinylome analysis reveals the involvement of lysine succinylation in metabolism in pathogenic Mycobacterium tuberculosis, Mol. Cell Proteomics 14 (2015) 796-811. [33] J.K. Diedrich, A.F. Pinto, J.R. Yates, Energy dependence of HCD on peptide fragmentation: stepped collisional energy finds the sweet spot, J. Am. Soc. Mass Spectr. 24 (2013) 1690-1699. [34] J.A. Vizcaíno, E.W. Deutsch, R. Wang, A. Csordas, F. Reisinger, D. Ríos, ProteomeXchange provides globally co-ordinated proteomics data submission and dissemination, Nat. Biotechnol. 30 (2014) 223–226. [35] E.C. Dimmer, R.P. Huntley, Y. Alam-Faruque, T. Sawford, C. Donovan, M. Martin, The UniProt-GO annotation database in 2011, Nucleic Acids Res. 40 (2012) D565-D570. [36] Y. Moriya, M. Itoh, S. Okuda, A.C. Yoshizawa, M.Kanehisa, KAAS: an automatic genome annotation and pathway reconstruction server, Nucleic Acids Res. 35 (2007) W182-185. [37] P. Horton, K. Park, T. Obayashi, N. Fujita, H. Harada, C.J. Adama-Collier, WoLF PSORT: protein localization predictor, Nucleic Acids Res. 35 (2007) W585-W587. [38] D.W. Huang, B.T. Sherman, R.A. Lempicki, Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists, Nucleic Acids Res. 37 (2009) 1-13. [39] P. Shannon, A. Markiel, O. Ozier, N. Baliga, J. Wang, D. Ramage, Cytoscape: a software environment for integrated models of biomolecular interaction networks, Genome Res. 13 (2003) 2498-2504. [40] D. Szklarczyk, A. Franceschini, M. Kuhn, M. Simonovic, A. Roth, P. Minguez, The STRING database in 2011: functional interaction networks of proteins, globally integrated and scored, Nucleic Acids Res. 39 (2011) 561-568. [41] D.B. Gary, W.H. Christopher, An automated method for finding molecular complexes in large protein interaction networks, BMC Bioinformatics 4 (2003) 2. [42] D.L. He, Q. Wang, M. Li, R.N. Damaris, X.L. Yi, Z.Y. Cheng, Global proteome analyses of lysine acetylation and succinylation reveal the widespread involvement of both modification in metabolism in the embryo of germinating rice seed, J. Proteome Res. 15 (2016) 879-890. [43] S.M. Zhen, X. Deng, J. Wang, G.R. Zhu, H. Cao, L.L. Yuan, First Comprehensive proteome analyses of lysine acetylation and succinylation in seedling leaves of Brachypodium distachyon L., Sci. Rep. 6 (2016) 31576. [44] X.P. Fang, W.Y. Chen, Y. Zhao, H.M. Zhang, C.Q. Yan, Global analysis of lysine acetylation in strawberry leaves, Front. Plant Sci. 6 (2015) 739. [45] K.N. Papanicolaou, B. O'Rourke, D.B. Foster, Metabolism leaves its mark on the powerhouse: recent progress in post-translational modifications of lysine in mitochondria, Front. Physiol. 5 (2014) 301. [46] J. Pan, R. Chen, C. Li, W. Li, Z. Ye, Global analysis of protein lysine succinylation profiles and their overlap with lysine acetylation in the marine bacterium Vibrio parahemolyticus, J. Proteome Res. 14 (2015) 4309-4318. 24

ACS Paragon Plus Environment

Page 25 of 36

629 630 631 632 633 634 635

Journal of Agricultural and Food Chemistry

[47] N. Nazemof, P. Couroux, T. Xing, L. S. Robert, Proteomic analysis of the mature Brassica stigma reveals proteins with diverse roles in vegetative and reproductive development, Plant Sci. 250 (2016)51-58. [48] R. Ivanov, T. Gaude, Endocytosis and endosomal regulation of the S-receptor kinase during the self-incompatibility response in Brassica oleracea, Plant Cell 21(2009) 2107-2117. [49] A. Kumar, B. Mcclure, Pollen-pistil interactions and the endomembrane system, J. Exp. Bot. 61(2010) 2001-2013. [50] M. Iwano, H. Shiba, K. Matoba, T. Miwa, M. Funato, T. Entani, Actin dynamics in papilla cells of Brassica rapa during self and cross-pollination, Plant Physiol. 144 (2007) 72-81.

636 637 638 639 640 641 642 643 644 645 646 647 648 649 650 651 652 653 654 655 656 657 658 25

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 26 of 36

659

Figure captions

660

Figure 1. The work flow of screening, enrichment and identification of strawberry

661

stigmata lysine succinylome by LC-MS/MS. An integrated strategy for global mapping of

662

lysine succinylation in strawberry stigmata.

663

Figure 2. Western blotting screening and identification of strawberry stigmata lysine

664

succinylome. (A) Equal amounts of total protein lysates (20 μg) were extracted from the

665

strawberry stigmata. Coomassie blue staining was used for the loading control. WB

666

experiments using the pan anti-phosphotyrosine, anti-acetyllysine, anti-succinyllysine,

667

anti-diglycine lysine, anti-mono-, dimethyllysine, anti-SUMO-1 or anti-SUMO-3 antibodies

668

were performed to assess phosphorylation (PTyr), acetylation (AcK), succinylation (SuK),

669

ubiquitination (GG), methylation (Me) and sumoylation (SUMO-1 / SUMO-3) levels in the

670

stigma protein, respectively. Equal amounts of protein were loaded into each lane. Each

671

experiment was performed in duplicate. (B) Venn diagram analysis of all the lysine

672

succinylated peptides and proteins from three biological replicates. (C) The distribution of

673

succinylated peptides based on their mass error. (D) The number of succinylated proteins

674

and sites identified in this study compared with previous studies.

675

Figure 3. Profile of identified succinylated and acetylated sites, peptides and proteins in

676

strawberry plants. (A) Venn diagrams illustrating the overlap between the SuK and AcK

677

proteins. (B) Venn diagrams illustrating the overlap between the SuK and AcK sites. (C)

678

Numbers of SuK and AcK sites, SuK and AcK peptides and SuK and AcK proteins. (D)

679

Distribution of SuK and AcK peptides in one protein. (E) Distribution of SuK and AcK

680

peptides based on their length.

681

Figure 4. Characterization of strawberry stigmata lysine succinylome. (A) Enrichment

682

analysis of all the identified lysine succinylated proteins based on the classification of GO 26

ACS Paragon Plus Environment

Page 27 of 36

Journal of Agricultural and Food Chemistry

683

annotation in terms of biological process, molecular function and cellular component. The

684

distribution and function of succinylated proteins identified in strawberry stigmata compared

685

with previous reported plant tissues and organs based on biological processes (B), molecular

686

functions (C) and subcellular locations (D).

687

Figure 5. PPI network of all the lysine succinylated proteins. Different geometric shapes

688

indicate various significantly enriched sub-clusters. Most proteins were clustered into one of

689

five groups, including the TCA cycle, oxidative phosphorylationand glycolysis.

690

Figure 6. Pollination Leads to increased hypersuccinylation of AP2 and clathrin. WB

691

showing succinylation of AP2 and clathrin immunoprecipitated from strawberry stigmata

692

protein at three time points (2 h, 0.5 h, and 0 h) before pollination and two time points (0.5 h

693

and 2 h) after pollination. Integrated density values werecalculated and are shown relative to

694

strawberry stigmata at 0 h. Each experiment was repeated three times IP:

695

Immunoprecipitation; WB: Western blotting; AP2: Assembly polypeptide 2; GAPDH:

696

Glyceraldehyde 3-phosphate dehydrogenase; SuK: Succinylated lysine.

697

Figure

698

clathrin/AP2–mediated vesicular transport of signal molecules. AP2/clathrin in

699

un-succinylated state will not combine with the receptor, and succinylated AP2/clathrin can

700

be further combined with the receptor-ligand complex to form the vesicles. Then the vesicles

701

can be budded, released and fused into the organelles membranes of strawberry stigma outer

702

layers cells. Thus the signal substances in pollen cells are successfully transported into the

703

stigma cell to further complete the pollination recognition process.

7.

A

proposed

model

depicting

how

succinylation

promotes

the

704 705

27

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 28 of 36

Table 1. Motif analysis of the identified succinylated peptides in strawberry stigmata. Motif

a

Foreground

Motif a

Score

Background

Matches

Size

Matches

Size

Fold Increaseb

.........FK..........

4.7

21

232

111

3495

2.85

.........IK..........

3.71

26

211

193

3384

2.16

.........LK..........

3.58

35

185

326

3191

1.85

.........VK.........

2.87

26

150

265

2865

1.87

Motif score, the probalility of the succinylated sites in some motifs and the value was set as -Log10(P

value). The confident cut-off score is 2.00. b

Fold increase is set as [matches /size in foregound]  [matches/size in background].

28

ACS Paragon Plus Environment

Page 29 of 36

Journal of Agricultural and Food Chemistry

Figure 1

1

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 30 of 36

Figure 2

1

ACS Paragon Plus Environment

Page 31 of 36

Journal of Agricultural and Food Chemistry

Figure 3

2

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 32 of 36

Figure 4

1

ACS Paragon Plus Environment

Page 33 of 36

Journal of Agricultural and Food Chemistry

Figure 5

1

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 34 of 36

Figure 6

2

ACS Paragon Plus Environment

Page 35 of 36

Journal of Agricultural and Food Chemistry

Figure 7

1

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 36 of 36

Graphic for table of contents

2

ACS Paragon Plus Environment