A Comparative Proteome Approach Reveals Metabolic Changes

4 days ago - In some industrial processes, cold and light stresses are recognized as two important environmental triggers for the transformation of my...
4 downloads 6 Views 1MB Size
Subscriber access provided by - Access paid by the | UCSB Libraries

Omics Technologies Applied to Agriculture and Food

Comparative proteome approach reveals metabolic changes associated with Flammulina velutipes mycelia in response to cold and light stress Jing-yu Liu, Ming-chang Chang, Jun-long Meng, Cui-ping Feng, and Yu Wang J. Agric. Food Chem., Just Accepted Manuscript • DOI: 10.1021/acs.jafc.8b00383 • Publication Date (Web): 27 Mar 2018 Downloaded from http://pubs.acs.org on March 27, 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 28

Journal of Agricultural and Food Chemistry

1

Comparative proteome approach reveals metabolic changes

2

associated with Flammulina velutipes mycelia in response to

3

cold and light stress

4

Jing-yu Liu,*,†,‡ Ming-chang Chang,*,†,‡ Jun-long Meng,†,‡ Cui-ping Feng,†,‡ Yu Wang†

5



6



College of Food Science and Engineer, Shanxi Agricultural University, Taigu 030801, China

Shanxi Engineering Research Center of Edible Fungi, Taigu 030801, China

7 8

ABSTRACT: In some industrial processes, cold and light stresses are recognized as two important

9

environmental triggers for the transformation of mycelia into fruit-bodies via intermediate primordia in

10

Flammulina velutipes cultivation. To gain insights into the regulation mechanism of F. velutipes mycelia

11

in response to cold and light stress, proteins expressed abundantly and characteristically at particular

12

stress states were investigated by using the isobaric tags for relative and absolute quantitation (iTRAQ)

13

labeling technique. Among the 1046 non-redundant proteins identified with high confidence, 264 proteins,

14

which were detected as differentially expressed proteins, were associated with 176 specific KEGG

15

pathways. In-depth data analysis revealed that the regulatory network underlying the cold and light

16

response mechanisms of F. velutipes mycelia were complex and multifaceted, as it included varied

17

functions such as rapid energy supply, the biosynthesis of lysine, phenylalanine, tyrosine and γ-

18

aminobutyric acid, calcium signal transduction process, dynein-dependent actin and microtubule

19

cytoskeleton formation, autolysis, oxidative stress adaption, pigment secretion, tissue and organ

20

morphogenesis, and other interesting stress-related processes. Insights about the proteins might shed light

21

on an intuitive understanding of the cold and light stress response mechanism underlying the fruiting

22

processes of F. velutipes. Furthermore, the data might also provide further insights into the stress

1

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 2 of 28

23

response mechanism of macro-fungi and valuable information for scientific improvement of some

24

mushroom cultivation techniques in practice.

25 26

KEYWORDS: Flammulina velutipes, environmental stress, proteome, iTRAQ, LC-MS/MS

27 28 29

■ INTRODUCTION

30

Due to high nutritional and medicinal value, mushrooms comprise an important macro-fungi class of

31

significant commercial relevance.1 However, basic stress physiological studies of mushrooms are given

32

less attention than their lower relatives, such as yeasts and filamentous fungi. 2-4 Flammulina velutipes

33

(Curt. ex Fr.) Sing, commonly known as golden-needle mushroom, is one of the most popular

34

industrialized cultivation mushrooms in the global market. 5 Meanwhile, its culture conditions can be

35

conveniently established in a biology laboratory by using some common refrigeration and light apparatus.

36

F. velutipes has been recognized as a potentially excellent model fungal species for basic and applied

37

studies on macro-fungi.6-9 Preliminary study has indicated that temperature reduction with weak light

38

irradiation is an effective regulation method for controlling fruiting process of theses fungi. 9-11 In some

39

industrial processes, cold and light stress are environmental triggers that has been routinely used to

40

induce transformation of mycelia to primordia and control both yield and quality of F. velutipes fruit-

41

bodies. Our previous study also showed that cold stress is one of the pivotal environmental factors that

42

influencing mycelium growth and fruit-body formation in F. velutipes, but few primordia came into being

43

during F. velutipes mycelium growth under cold stress in complete darkness. 12 This indicates that F.

44

velutipes responds to the stimuli of cold and light stress by orchestrating a range of physiological and

45

metabolic processes. An understanding of these processes at the molecular level will not only be useful 2

ACS Paragon Plus Environment

Page 3 of 28

Journal of Agricultural and Food Chemistry

46

in the cultivation practice of F. velutipes, but also be valuable in providing basic knowledge on the

47

fruiting mechanism of F. velutipes. Unfortunately, the molecular biology of F. velutipes has been studied

48

poorly until present. For example, only a handful of genes involved the fruiting process of F. velutipes

49

have been isolated and preliminarily characterized.7,8,13 Current knowledge of the cold and light stress

50

response mechanism underlying the fruiting processes of F. velutipes remains fragmented. We are still

51

far from understanding the mechanism of mycelia in response to cold and light stress in these fungi.

52

Proteomics techniques have emerged as powerful tools for the global assessment and measurement

53

of systematic flux changes of total cell protein expression at particular biological states. The standard

54

gel-based method, one-/two-dimensional (1D/2D) gel electrophoresis followed by mass spectrometry

55

(MS), has been chosen and widely practiced for decades.14 Recently, with the advent of chromatography

56

and MS techniques, high-throughput shotgun proteomics, liquid chromatograph (LC) followed by MS,

57

are being increasingly practiced.15 The isobaric tags for relative and absolute quantitation (iTRAQ)

58

labeling coupled to 2D LC-MS/MS technique, one of the LC-based proteome techniques, has been

59

developed and proven to be a more suitable tool for analyzing chronological proteomic changes in

60

complex developmental processes of fungi.4,9,12,16 In this study, we used the iTRAQ labeling coupled to

61

2D LC-MS/MS technique to assess the global chronological changes in protein expression patterns of F.

62

velutipes mycelia response to cold and light stress. Furthermore, these protein expression patterns were

63

systematically investigated, analyzed and compared to better understand the stress response mechanism

64

in F. velutipes mycelia.

65

■ MATERIALS AND METHODS

66

Strains and Culture Conditions

67

Brown fruit-body dikaryon strain Fv13 of F. velutipes were cultured on culture plates (potato 20%,

3

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 4 of 28

68

glucose 2% and agar powder 2%) at 23-25 °C in darkness. A total of 285 plates were divided randomly

69

into three groups: control group, cold stress (CS) group, and cold and light stress (CLS) group. 20 days

70

after incubation, the whole plate mycelia of the control group were collected as control samples.

71

Subsequently, 95 plates of CS group were transferred and cultured at 12-15 ℃ for 22 days in darkness,

72

and the whole plate mycelia were collected as CS samples. 95 plates of CLS group were transferred and

73

cultured at 12-15 ℃ in low light (200-300 lx) for 22 days, and the whole plate mycelia were collected as

74

CLS samples, not including primordia. All samples were stored at -80 °C until protein extraction.

75

Protein Preparation and iTRAQ Labeling

76

Protein preparation was performed following published protocols with minor modification. 9 Briefly,

77

the control, CS and CLS samples were extracted by STD buffer (SDS 4%, Tris-HCl 150 mM, DTT 1

78

mM, pH 8.0) and UT buffer (Urea 8 M, Tris-HCl 150 mM, pH 8.0). Total protein concentration was

79

determined by agarose gel and Bradford Assay Kit (Applied Biosystem). Protein digestion was

80

performed by trypsin solution (5% trypsin, Promega) overnight at 37 °C. The resulting tryptic peptides

81

were collected and incubated with 8-plex iTRAQ reagents application kit (AB Sciex) according to the

82

manufacturer’s instructions. CLS samples were labeled with reagent 115; control samples and CS

83

samples were labeled with reagent 113 and 114, respectively.

84

2D LC-MS/MS Analysis and Data Analysis

85

2D LC-MS/MS analysis was performed following published protocols with minor modification. 12

86

Briefly, the peptides were fractionated using a C18 column (waters bed 2.1 × 50 mm, 1.7 μm, Applied

87

Biosystem). Peptides were eluted with solvent A (20 mM ammonium formate) adjusted to pH 10 and a

88

linear gradient solvent B (ACN, 5-35 %) at a flow rate of 600 μL/min. The absorbance was monitored at

89

214 nm. The peptide fractions were collected together to obtain 20 final fractions and subsequently

4

ACS Paragon Plus Environment

Page 5 of 28

Journal of Agricultural and Food Chemistry

90

separated by nano-HPLC (Eksigent Technologies) on the secondary RP analytical column (Eksigent, 3

91

μm, 75 μm × 150 mm) equipped with Triple TOF 5600 MS (Thermo Finnigan). The gradient conditions

92

were a linear phase B (98% ACN, 0.1% formic acid) gradient of 5-45 % from 5-100 min at a total

93

controlled flow rate of 300 nL/min. Electrospray voltage of the mass spectrometer was 2.5 kV versus the

94

inlet. The scanning range of the tandem mass spectrometry was 100 to 1250 m/z. A dynamic exclusion

95

time of 25 s per cycle was used to select 25 most intense precursors for fragmentation.

96

Tandem mass spectrometry extraction was processed by Proteome Wizard version 3.0.10185 64-bit.

97

The MS data of three above mentioned samples were processed using MASCOT (Matrix Science,

98

version 2.3.0) software. Agaricales. Fasta in UniProt (http://www.uniprot.org/) was used as a reference

99

database. Data were normalized by Scaffold Q+ (software version Scaffold-4.4.6). Protein identifications

100

and protein probabilities were determined by the Scaffold Local FDR algorithm, a false discovery rate

101

(FDR) < 0.05, and the protein prophet algorithm, respectively. 17 In this study, FDR of 2 were chosen as the criteria for protein identification. In all samples, the matrix

103

[0.000,0.0188,0.931,0.0490,0.001000];

[0.000,0.0282,0.932,0.0390,0.000700];

104

[0.000,0.000,0.929,0.0689,0.00220];

[0.000,0.00940,0.930,0.0590,0.00160];

105

[0.000600,0.0377,0.933,0.0288,0.000];

[0.00140,0.0566,0.933,0.00870,0.000];

106

[0.000900,0.0471,0.933,0.0188,0.000];

[0.00270,0.0744,0.921,0.00180,0.000]

107

[0.000,0.000,0.000,0.000,0.000] were adopted to correct channels according to the published algorithm

108

in i-Tracker.18

109

Functional Categorization Analysis

110

For functional characterization of the differentially expressed proteins (DEPs), the abundance of

111

DEPs were detected by p < 0.05 (one-way ANOVA test) and 1.5 fold cut-off. Gene Ontology (GO)

5

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 6 of 28

112

enrichment of down- or up-regulated DEPs were subsequently conducted according to the information

113

from Swiss-Prot and EBI database (Agaricus). To gain GO mapping and annotation information, data of

114

the DEPs was sequentially uploaded to agriGO (http://bioinfo.cau.edu.cn/agriGO/) by using Blast2GO

115

(software version 2.6.6). Enzyme codes of DEPs were sequentially uploaded by our in-house compiled

116

PERL program KEGG.pl to KEGG (the Kyoto Encyclopedia of Genes and Genomes) website

117

(http://www.genome.jp/kegg/), and mapped to annotated sequences and metabolic pathways in KEGG.

118

RNA Extraction, Reverse Transcription and qRT-PCR analysis

119

Total RNA of the three samples, control, CS and CLS was extracted according to the manufacturer’s

120

instructions of TRIZOL Reagent (Invitrogen). Genomic DNA was removed using DNase I treatment

121

(Promega). Quality of purified RNA was accepted if they had a ratio of 1.8-2.0 at A260/A280. Total RNA

122

was reverse transcribed according to the manufacturer’s instructions of Prime-Script® RT reagent Kit

123

(Takara). Quantitative real-time Polymerase Chain Reaction (qRT-PCR) was sequentially conducted

124

using SYBR Premix EX TaqTM II (Takara) and a Step One Real-time PCR System (Applied Biosystems).

125

To validate our proteome data, the expression patterns of 18 mRNAs, each corresponding to 18

126

representative DEPs including nine up-regulated and nine down-regulated proteins, respectively. The

127

data were then normalized to the expression levels of the endogenous control gene, GAPDH. Primer pairs

128

were designed for qRT-PCR analysis using the Primer 5.0 software. Primer data of the 18 representative

129

DEPs are summarized in Table S1. Our data indicated that the amplifications were specific.

130

■ RESULTS

131

Morphological Changes of Flammulina velutipes Mycelia in Response to Cold and Light Stress

132

After being incubated at 23-25°C in darkness for 20 days, the plates of the control group (Fig. 1A)

133

were fully colonized with vigorous aerial mycelia (Fig. 1a). After CS treatment, aerial mycelia of CS

6

ACS Paragon Plus Environment

Page 7 of 28

Journal of Agricultural and Food Chemistry

134

group plates (Fig. 1B) became flattened and secreted some yellow physiological water (Fig. 1b). After

135

CLS treatment, a few primordia came into being, but the aerial mycelia of the plates secreted less yellow

136

physiological water (Fig. 1C, 1c).

137

B

A

C

138 139 140 141 142 143

a

144

c

b

145 146 147

Figure 1. Morphology changes of Flammulina velutipes mycelia in response to cold and light stress. A: control

148

sample; B: cold stress sample; C: cold and light stress sample; a: aerial hyphae; b: yellow physiological water; c:

149

primordia. All strains were inoculated at the same time.

150

Protein Identification and Quantification

151

According to the above mentioned criteria for protein identification, 1046 non-redundant proteins

152

were obtained from three biological replicates for further comparative analysis. According to a criterion

153

of 95% significance and 1.5 fold cutoff, following which a total of 550 non-redundant proteins were

154

identified as DEPs. The Venn diagram analysis was subjected to illustrate the 160 up-regulated and 409

155

down-regulated DEPs distributions and their overlap between CS and CLS samples (Figure 2).

156 157 158 159 160 161 162 163

A

B

12

CS

1

147

63

CLS

CS

37

309

CLS

Figure 2. Venn diagrams and expression of 551 DEPs identified under cold and light stress in mycelia of Flammulina

7

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 8 of 28

164

velutipes. The numbers of DEPs with up-regulation (A) and down-regulation (B) under given development stages

165

are shown in different segments. CS: cold stress; CLS: cold and light stress.

166

Among the 550 DEPs, 286 DEPs, which were characterized as predicted, uncharacterized or putative

167

proteins under the reference database, were discarded. In general, two clearly different expression

168

patterns in CS and CLS samples of F. velutipes were generalized among the remainder 264 DEPs. A total

169

of 70 DEPs were up-regulated in at least one of two samples (CS or CLS), including one common DEPs

170

(acetyltransferase component of pyruvate dehydrogenase complex) in two samples (CS and CLS), while

171

two and 66 DEPs were significantly up-regulated in CS and CLS respectively. Among 204 down-

172

regulated DEPs in at least one of the two samples, 21 DEPs were common to two samples (CS and CLS),

173

while 33 and 149 DEPs were particularly down-regulated in CS and CLS, respectively. The data of up-

174

regulated DEPs are summarized in Table 1 and Table 2. The data of down-regulated DEPs are

175

summarized in Supplementary Table S2, Table S3 and Table S4.

176

Table 1. List of up-regulated proteins under only cold stress in mycelia of Flammulina velutipes. Uniprot ID a

No. 1

D6RQF5_COPC7

Protein name Aminoadipate-semialdehyde

Gene Name

Fold change b CS1/C1

CS2/C2

CS3/C3

CC1G_15694

1.56

1.91

1.91

dehydrogenase 2

V2X4P6_MONRO

Acetyl-hydrolase

Moror_11650

1.71

1.91

1.75

3

D8PKQ0_SCHCM

Acetyltransferase component

SCHCODRAFT_81014

1.67

1.86

1.82

of pyruvate dehydrogenase complex

177 178 179 180 181 No. 1

a

Uniprot: http://www.uniprot.org/.

b

CS: cold stress; C: Control.

Table 2. List of up-regulated proteins under only cold and light stress in mycelia of Flammulina velutipes. Uniprot ID a Q5EGJ1_FLAVE

Protein name Translation elongation factor

Gene Name

Fold change b CLS1/C1

CLS2/C2

CLS3/C3

tef1

1.67

1.57

1.71

EF1-alpha (Fragment) 2

V2WVN8_MONRO

Hsp70 chaperone

Moror_9965

2.50

2.33

2.33

3

D8QLU7_SCHCM

Calmodulin

SCHCODRAFT_80005

1.53

1.67

1.65

8

ACS Paragon Plus Environment

Page 9 of 28

Journal of Agricultural and Food Chemistry

4

D8PR32_SCHCM

40S ribosomal protein S1

RPS1

2.00

1.67

1.84

5

B0DI83_LACBS

Pyruvate carboxylase

PYC

1.60

1.57

1.60

6

B0DHG8_LACBS

Beta-tubulin 1 tubb1

LbTUBb1

2.63

2.30

1.85

7

V2XAN6_MONRO

Calnexin

Moror_7233

1.63

1.59

1.82

8

A0A067P3W4_PLEOS

Histone H4

PLEOSDRAFT_1088653

1.63

1.67

1.60

9

A0A067NVC3_PLEOS

Histone H2B

PLEOSDRAFT_1088190

1.82

1.57

1.56

10

A0A067SYC7_9AGAR

60S ribosomal protein L20

GALMADRAFT_142989

2.33

2.71

2.17

11

V2WVN7_MONRO

Inorganic diphosphatase

Moror_16434

3.00

1.86

1.50

12

A0A067SEN5_9AGAR

Serine hydroxymethyltransferase

GALMADRAFT_257154

1.61

1.67

1.60

13

B0CPE4_LACBS

Cell division control/GTP binding protein

LACBIDRAFT_291428

2.09

2.53

2.31

14

B0DZY3_LACBS

Catalase

LbCAT

4.25

3.75

4.50

15

A8NSX2_COPC7

Ribosomal protein rpl23a

CC1G_12270

2.00

3.33

2.60

16

V2YIE8_MONRO

Alpha-1,4 glucan phosphorylase

Moror_2611

1.53

1.71

1.50

17

D8PXF2_SCHCM

Succinate dehydrogenase

SCHCODRAFT_66862

1.71

1.57

1.80

[ubiquinone] iron-sulfur subunit, mitochondrial 18

R4HKU2_FLAVE

Myosin regulatory light chain

cdc4

1.88

2.00

1.53

19

A0A067NJV8_PLEOS

Autophagy-related protein

PLEOSDRAFT_1097703

1.75

1.20

2.00

20

D8PW86_SCHCM

Phosphoglycerate kinase

SCHCODRAFT_74370

1.63

2.40

2.33

21

V2YDU9_MONRO

Scf complex subunit skp1

Moror_859

2.25

2.11

2.13

22

A8NZV3_COPC7

Superoxide dismutase

CC1G_06963

3.00

1.67

3.00

23

D2JY86_FLAVE

Initiation factor 5a

--

1.60

1.58

1.55

24

V2XI69_MONRO

Proteasome subunit alpha type

Moror_1703

1.50

1.50

1.50

25

V2XUP9_MONRO

2-methylcitrate dehydratase

Moror_6806

1.63

2.17

3.00

26

A0A067N8U4_PLEOS

6-phosphogluconate

PLEOSDRAFT_41725

3.00

3.20

2.71

dehydrogenase, decarboxylating 27

V2XD70_MONRO

60s ribosomal protein l17

Moror_13629

2.33

2.75

2.67

28

A0A067TK21_9AGAR

Proteasome subunit alpha type

GALMADRAFT_235810

1.75

2.00

1.72

29

B0CQH3_LACBS

Acetyltransferase component of

LACBIDRAFT_378984

1.89

1.95

--

pyruvate dehydrogenase complex 30

V2XYL9_MONRO

Cyanide hydratase

Moror_894

1.75

2.00

1.81

31

A8NBN6_COPC7

60S ribosomal protein L27a

CC1G_02496

2.71

2.71

1.78

32

A8N0P0_COPC7

Peptidyl-prolyl cis-trans

CC1G_04343

1.50

2.50

2.17

gpd

1.63

1.57

1.64

isomerase 33

D2SRR0_9AGAR

Glyceraldehyde-3-phosphate

34

V2X183_MONRO

ATP synthase d subunit

Moror_9086

3.50

2.33

2.67

35

B0CU44_LACBS

26S proteasome regulatory

PATPA16201

1.63

1.82

1.71

dehydrogenase (Fragment)

complex, ATPase RPT4 36

A0A067NIX1_PLEOS

Malate dehydrogenase

PLEOSDRAFT_1090652

2.25

2.47

2.32

37

V2WVA6_MONRO

Aryl-alcohol dehydrogenase

Moror_4055

2.83

2.39

2.41

38

A8P7B9_COPC7

Calcium-dependent protein

CC1G_08194

1.75

1.86

1.82

serine/threonine phosphatase 39

V2WZD6_MONRO

Transketolase

Moror_2319

1.62

1.58

1.57

40

V2WRP8_MONRO

60s ribosomal protein l2

Moror_1359

1.63

1.72

1.77

9

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

41

D8QCI6_SCHCM

Peptidyl-prolyl cis-trans

Page 10 of 28

SCHCODRAFT_16915

2.25

3.33

2.50

isomerase 42

D8Q077_SCHCM

Isocitrate lyase

SCHCODRAFT_66738

1.55

1.67

1.73

43

V2XH67_MONRO

Elongation factor 2

Moror_4797

2.64

1.75

1.86

44

V2XTA5_MONRO

Proteasome 26s subunit ATPase 3

Moror_12489

4.20

2.43

2.75

45

B0CUA0_LACBS

Actin filament-coating protein

LACBIDRAFT_176790

3.08

3.11

2.78

tropomyosin 46

V2XPR5_MONRO

Dynein light chain cytoplasmic

Moror_14124

1.60

1.50

1.50

47

A0A067SXJ7_9AGAR

Proliferating cell nuclear antigen

GALMADRAFT_102448

1.78

1.50

2.00

48

K5WMW0_AGABU

Transcription elongation factor

AGABI1DRAFT_115721

2.71

3.38

2.60

TEF EF1B 49

V2XT09_MONRO

Fe-containing alcohol

Moror_7341

2.42

2.75

1.90

50

V2X429_MONRO

Glyceraldehyde-3-phosphate

Moror_13000

2.67

3.00

2.00

dehydrogenase 51

A8N1B9_COPC7

Ribosomal protein L4/L1

CC1G_10540

3.67

2.66

--

52

V2XDI0_MONRO

Het-c2 protein

Moror_16482

2.00

1.91

1.67

53

V2X892_MONRO

Sec14 cytosolic factor

Moror_3846

1.81

1.69

1.67

54

B0D3H6_LACBS

Glutamate decarboxylase

LACBIDRAFT_315921

1.95

1.55

1.86

55

A8NHX8_COPC7

Alpha-mannosidase

CC1G_01530

1.65

1.57

1.59

56

A8NN10_COPC7

Class V chitinase ChiB1

CC1G_09921

3.00

3.60

3.22

57

A8NCY6_COPC7

60S ribosomal protein L10

CC1G_08610

2.50

2.00

1.71

58

V2XI54_MONRO

Arg-6 protein

Moror_4407

--

1.83

1.64

59

V2YBK4_MONRO

Microtubule binding protein

Moror_5349

1.80

2.30

2.00

60

B0CVX0_LACBS

Isocitrate lyase

LACBIDRAFT_171643

--

2.67

3.00

61

V2XG31_MONRO

Proteasome subunit beta type

Moror_369

1.60

1.74

1.71

62

V2XWY1_MONRO

Heat shock protein

--

2.07

--

2.13

63

G8A534_FLAVE

Beta-galactosidase

Moror_17812

1.67

1.50

1.56

64

A0A067T146_9AGAR

Eukaryotic translation initiation

HCR1

2.43

2.22

2.06

65

A0A067NJ33_PLEOS

INT6

1.71

1.64

1.58

factor 3 subunit J Eukaryotic translation initiation factor 3 subunit E 66 67

V2XAU1_MONRO

Serine hydroxymethyltransferase

Moror_7089

1.72

1.53

1.52

D8PRY7_SCHCM

Tubulin alpha-1A chain

SCHCODRAFT_9815

5.75

6.67

5.88

182 183

a

Uniprot: http://www.uniprot.org/.

b

CLS: cold and light stress; C: Control.

184

Functional Categorization Analysis

185

Among the 550 DEPs, functional information for the 264 DEPs was characterized by searching for

186

homologous proteins against the genomics information category of Agaricales. Fasta (UniProt) database.

187

The results of enrichment analysis showed that under the category of the cellular component, DEPs were

188

mainly found in cytoplasm, ribosome, membrane, nucleus, and proteasome complex. In the biological

10

ACS Paragon Plus Environment

Page 11 of 28

Journal of Agricultural and Food Chemistry

189

process, DEPs were mainly enriched in translation, protein folding, carbohydrate metabolic process, fatty

190

acid biosynthetic process, response to stress and signal transduction. In molecular function, DEPs were

191

mostly related to binding, enzyme activity and structural constituent of ribosome. Unfortunately, 13 of

192

the DEPs had unknown locations. The GO mapping and annotation data of up-regulated and down-

193

regulated DEPs are summarized and shown in Figure 3A and 3B, respectively.

A

Cellular component

B

Cellular component

194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222

Biological process

Molecular function

Biological process

Molecular function

223

Figure 3. Bioinformatics analysis of the above mentioned differentially expressed proteins through Gene Ontology

224

(GO) in three domains: cellular component, biological process and molecular function. The statistics at GO level 2

11

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

225

is shown in this figure. A: Up-regulated proteins; B: Down-regulated proteins.

226

Metabolic Pathway Analysis

Page 12 of 28

227

The data of GO analysis demonstrated that the DEPs during the transformation of mycelia into

228

primordia of F. velutipes under cold and light stress were influenced by a crucial and complicated cellular

229

process, particularly the metabolic processes. To gain further insights into pathway information of the

230

metabolic processes, data of the 264 DEPs were sequentially uploaded to KEGG website to annotate

231

sequences and metabolic pathways. The result showed that 139 (52.65 %) of the 264 DEPs were

232

associated with 176 specific KEGG pathways. Further analysis of 70 up-regulated DEPs showed that 42

233

DEPs (60 %) were mainly distributed in 99 specific KEGG pathways. In general, one (aminoadipate-

234

semialdehyde dehydrogenase) of the two up-regulated DEPs at CS samples were associated with

235

biosynthesis of amino acids (path: ko01230), lysine biosynthesis (path: ko00300) and lysine degradation

236

(path: ko00310). 40 of the 67 DEPs that were significantly up-regulated in CLS sample were distributed

237

in 97 specific KEGG pathways and mainly involved in signaling pathway (path: ko04010, path: ko04014,

238

path: ko04020, path: ko04022, path: ko04068, path: ko04152, path: ko04310, path: ko04350,

239

path:ko04910, path: ko04921, path: ko04922), amino acid biosynthesis and metabolism (path: ko00250,

240

path: ko00330, path:ko00460, path: ko00520, path: ko01230), energy metabolism (path: ko00010, path:

241

ko00020, path: 00030, path: ko00190, path: ko00630), glutathione metabolism (path: ko00480) (Figure

242

4), proteasome (path: ko03050), ribosome (path: ko03010) and cell cycle (path: ko04111). Further

243

functions of the DEPs are discussed in the following discussion sections.

244 245

12

ACS Paragon Plus Environment

Page 13 of 28

Journal of Agricultural and Food Chemistry

246 247 248 249 250 251 252 253 254 255 256 257 258 259

Figure 4. Representative metabolic pathway maps of differentially expressed proteins involved in glutathione

260

metabolism in KEGG. The colored enzyme codes were noted as follows: 6-phosphogluconate dehydrogenase (EC:

261

1.1.1.44 and 1.1.1.343), ribonucleoside-diphosphate reductase subunit (EC: 1.17.4.1), spermidine synthase (EC:

262

2.5.1.16), gamma-glutamyltranspeptidase / glutathione hydrolase (EC: 2.3.2.2 and 3.4.19.13). Red color: Up-

263

regulated proteins; Green color: Down-regulated proteins; Blue color: other identified enzymes.

264

Transcriptional Expression Pattern Analysis of Representative Differentially Abundant Proteins

265

Due to lack of specific antibodies for most proteins in basidiomycete studies, the western blotting

266

technique was limited to confirm only the most edible mushroom proteome data. In the study, we

267

determined and compared transcriptional expression patterns of 18 mRNAs, each corresponding to 18

268

representative DEPs, by qRT-PCR technique. Our data showed that 13 genes, encoding 60S ribosomal

269

protein L20, Serine hydroxymethyltransferase, malate dehydrogenase, elongation factor 2, proteasome

270

26s subunit ATPase 3, heat shock protein sks2, Tryptophan synthase, 2-cysteine peroxiredoxin, glycoside

271

hydrolase family 13 protein, V-type proton ATPase subunit a, alcohol oxidase-like protein, Glutathione-

13

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 14 of 28

272

disulfide reductase and S-adenosylmethionine synthase exhibited similar expression patterns with

273

representative protein expression patterns in the control, CS and CLS samples (Figure 5A). However, the

274

mRNA expression patterns of three genes, encoding inorganic diphosphatase, transketolase and 2-

275

methylcitrate dehydratase, were not consistent with their corresponding up-regulated proteins. Two of

276

the 18 genes, encoding superoxide dismutase and Hsp70 chaperone, showed opposite expression patterns

277

between their mRNAs and proteins (Figure 5B). This result is in accordance with previous reports due

278

to translational or post-translational regulation.9,12,19

279 280 281 282 283 284 285 286 287 288 289

Figure 5. Transcriptional expression analysis of representative proteins as revealed by qRT-PCR. A: Relative mRNA

290

expression patterns of representative up-regulated proteins; B: Relative mRNA expression patterns of representative

291

down-regulated proteins. CS: cold stress; CLS: cold and light stress. The data were then normalized to the expression

292

levels of the endogenous control gene, GAPDH (n = 3, mean ±standard error of the mean, *p-value < 0.05).

293

■ DISCUSSION

294

In all three life super-kingdoms (archaea, bacteria and eukaryotes), the cellular stress response have

295

been linked to essential processing and stability of protein and DNA. 20 Development in proteomic

296

analysis during the last decade provided insights into a systematic understanding of the complex

297

developmental events of filamentous fungi and yeast at the molecular level in response to cell stress

298

exposure.20,21 However, unlike other lower model fungi, the stress proteome of micro-fungi, especially

14

ACS Paragon Plus Environment

Page 15 of 28

Journal of Agricultural and Food Chemistry

299

mushrooms, remains largely uncharacterized. In micro-fungi, F. velutipes, a tetrapolar basidiomycete

300

mushroom, has emerged as a potential model system for studying basic mechanisms of genetics,

301

development, and cellular stress response.6,9,11,12 In this study, the stress proteome associated with F.

302

velutipes mycelia in response to cold and light stress were investigated by using iTRAQ labeling coupled

303

to 2D LC-MS/MS technique. Our data showed that some interesting up-regulated DEPs might play a key

304

role in the response mechanism of F. velutipes mycelia. The biological relevance of these up-regulated

305

DEPs is analyzed and discussed under some of the common central metabolism tasks, including carbon

306

and energy metabolism, biosynthesis and metabolism of amino acids, signal transduction pathways, and

307

other interesting stress-related processes.

308

Increased Abundance of Stress-related Proteins of Carbon and Energy Metabolism

309

Our results indicated complex flux changes of cell protein expression systems associated with F.

310

velutipes mycelia in response to cold and light stress in carbon and energy metabolism. The data showed

311

that acetyltransferase component of pyruvate dehydrogenase complex (ACPDC) had higher expression

312

levels under both CS and CLS samples. ACPDC participates in pyruvate metabolism and converts

313

pyruvate to acetyl CoA. Acetyl-CoA is an important molecule that participates in carbohydrate, protein

314

and lipid metabolism. Acetyl-hydrolase, which catalyzes the conversion of acetyl-CoA and H2O into

315

Coenzyme A and acetate, had higher expression levels in mycelia in response to CS, but their expression

316

levels returned to normal under CLS. CoA is the main input for the citric acid cycle (TCA cycle) and

317

fatty acid biosynthesis which are required for fungal mycelia growth and fruit-body formation.22 In fungi,

318

accumulation of acetate also reduced growth by repressing mating and filamentation, and promoted the

319

building block for fatty acid biosynthesis.23 Aldehyde dehydrogenase also serves as the last enzyme in

320

the formation reaction of Neurospora carotenoid.24 This finding revealed that acetyl-hydrolase might play

15

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 16 of 28

321

an important role in improving the mycelium adaptability, secreting pigment and reducing damage under

322

CS environment in F. velutipes.

323

We found that some crucial proteins participated in processes like central energy metabolisms for

324

sugar, glycolysis, citric acid cycle, pentose phosphate pathway and glucose release from starch, also had

325

higher expression levels in CLS sample, such as α-1,4 glucan phosphorylase, 6-phosphogluconate

326

dehydrogenase,

327

dehydrogenase, isocitrate lyase, 2-methylcitrate dehydratase, malate dehydrogenase, succinate

328

dehydrogenase [ubiquinone] iron-sulfur subunit and transketolase. The catalysates of these enzymes

329

serve as important intermediates for gluconeogenesis, ATP and fatty acid biosynthesis that play key roles

330

in several processes of fungi, such as the synthesis of biological membranes, cellular structure, protein

331

localization, energy production and storage for fungal growth. 22,25,26 On the basis that low activity of the

332

enzymes in CS sample, high activities of ATP synthase d subunit and ribosomal protein, and

333

morphological characters of F. velutipes mycelia in CLS samples. One plausible explanation for the

334

result is light irradiation stimulation plays a crucial role in the rapid energy supply for the transformation

335

of mycelia into primordia in F. velutipes.

336

Increased Abundance of Stress-related Proteins of Amino Acid Biosynthesis and Metabolism

phosphoglycerate

kinase,

pyruvate

carboxylase,

glyceraldehyde-3-phosphate

337

Amino acids function as building blocks and energy in fungus growth and development. Stress-

338

related metabolites of amino acids also have the potential to integrate environmental stress conditions

339

into physiological outputs and intracellular signals. In this study, our data showed that six proteins

340

associated with amino acid metabolism, including aminoadipate-semialdehyde dehydrogenase

341

(AASDH), aryl-alcohol dehydrogenase (AAD), serine hydroxymethyltransferase (SHMT), glutamate

342

decarboxylase (GD), cyanide hydratase (CHT) and peptidyl-prolyl cis-trans isomerase, increased in

16

ACS Paragon Plus Environment

Page 17 of 28

Journal of Agricultural and Food Chemistry

343

abundance during mycelia response to cold and light stress in F. velutipes. Among the six proteins,

344

AASDH had higher expression levels in CS samples. AASDH is an important enzyme for the

345

biosynthesis of lysine in yeasts and higher fungi.27 Accumulation of lysine is beneficial in the convert of

346

pyridoxal phosphate (PLP) from vitamin B6, where PLP plays a key role as a coenzyme in transamination

347

reactions of amino acids and amino sugar.28

348

Our data showed that the activities of AAD, SHMT, GD, and CHT were remarkably increased in

349

CLS samples. In wood-saprophyte fungi and some yeast, AAD serves as an important enzyme for

350

detoxifying aromatic aldehydes released during lignin degradation.29 The catalysate of AAD are involved

351

in the synthesis of phenylalanine and tyrosine that might be advantageous to stress responses in

352

microorganism.30 SHMT, one of PLP-enzymes, catalyzes the inter-conversion of L-Ser and

353

tetrahydrofolate (THF) to Gly and 5,10-CH2-THF. This reaction plays an important role in providing the

354

largest one-carbon fragments available to the cell.31 GD is an enzyme that catalyzes the conversion of

355

glutamate to γ-aminobutyric acid (GABA). GABA levels could serve as the key signal connecting

356

cellular carbon metabolism to developmental or environmental responses of various biotic and abiotic

357

stresses in plants and fungi.32,33 In Stagonospora nodorum, GABA has positive effects on asexual

358

development of mycelia in response to environmental stresses.34 Recent work showed that fruit-bodies

359

of F. velutipes contained the highest amount of GABA among 20 different mushrooms. 35 CHTs that are

360

the key enzymes catalyzing the degradation of cyanide are widely found in filamentous fungi and few

361

other members of Basidiomycota36 and were postulated to serve as nitrogenase.37 In this study, we first

362

found the CHT increased in abundance under light induction in F. velutipes mycelia. This finding

363

revealed that these enzymes mediated regulation of the biosynthesis of lysine, phenylalanine, tyrosine

364

and GABA, the degradation of cyanide, and one-carbon fragments available, which might play important

17

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 18 of 28

365

roles in developmental and environmental responses F. velutipes mycelia under CLS environment.

366

Increased Abundance of Stress-related Proteins of Signaling Pathways and Other Biological

367

Processes

368

In fungi, some conserved signaling pathways that control the adaptation to the environmental changes

369

is crucial for cell cycle regulation, reproduction and stress response. Here the result showed a complex

370

picture of signaling pathways associated with mycelia response to light stress in F. velutipes. Our data

371

showed that some interesting stress-related DEPs, such as α-mannosidase, calnexin, calmodulin,

372

calcium-dependent protein serine/threonine phosphatase, cell division control/GTP binding protein, beta-

373

tubulin 1 tubb1, tubulin alpha-1A chain, actin filament-coating protein tropomyosin, microtubule binding

374

protein, dynein light chain cytoplasmic and myosin regulatory light chain, remarkably increased

375

abundance in CLS samples. In the filamentous fungus Aspergillus oryzae, α-mannosidases act on N-

376

linked mannooligosaccharides that are involved in the formation of the membrane. 38 Calnexin is a

377

molecular chaperone that modulates the Ca2+ homeostasis and calcium signal transduction process in

378

endoplasm reticulum. Calnexin is required for mycelium growth and conidiation in the filamentous

379

fungus Aspergillus nidulans under some endoplasmic reticulum stress.39 Calmodulin and calcium-

380

dependent protein serine/threonine phosphatase play critical roles in the regulation of stress responses

381

via the calcium-calcineurin pathway that controls fundamental aspects of survival, morphogenesis and

382

signaling cascades in fungi.40 In the fungus Ustilago, calcium signaling was proved to participate in

383

dynein-dependent microtubule organization that drives tissue morphogenesis.41 GTP-binding proteins

384

coordinate cell division and signaling and participate in forming filaments in eukaryotic cells. In fungi,

385

the proteins function via scaffolding of cytosolic proteins and associate with actin and microtubule

386

cytoskeleton networks to form higher-order structures.42 In eukaryotes, alpha and beta -tubulin are

18

ACS Paragon Plus Environment

Page 19 of 28

Journal of Agricultural and Food Chemistry

387

necessary for the assembly of dynamic microtubules.43 Actin filament-coating protein tropomyosin and

388

microtubule binding protein contribute to the binding of the filaments and regulate the diverse range of

389

actin and microtubule cytoskeleton structures, respectively. The interdependent collaboration of actin

390

and microtubule cytoskeletons is crucial for cytokinesis, polarity maintenance, membrane remodeling

391

and tissue morphogenesis during fungal growth.44 Myosin regulatory light chain functions via the action

392

of myosin kinases and reversible phosphorylation to modulate basal actomyosin oscillations that have

393

emerged as key regulators of tissue and organ morphogenesis.45 Cytoplasmic dynein light chain functions

394

in peroxisome biogenesis that attributes to assemble the dynein motor in yeast.46 Cytoplasmic dynein is

395

required for the microtubule-based transport and spatial organization in filamentous fungi.47 Considered

396

together, these findings suggested the coordinated organization of calcium signal transduction process

397

and dynein-dependent actin and microtubule cytoskeletons formation might play a crucial role in the

398

regulation of survival, tissue morphogenesis and signaling cascades during F. velutipes mycelia under

399

CLS environment.

400

In CLS samples, we also found a dramatic increase in the activities of class V chitinase ChiB1, β-

401

galactosidase, autophagy-related protein and Hsp70 chaperone. In synergy with the β-glucanases, ChiB1,

402

a chitinase, may play an important role in the developing fruit-body of Coprinopsis cinerea, especially

403

pileus autolysis.48 In fungi, autophagy-related proteins are essential for triggering autophagy-mediated

404

processes to regulate nuclear dynamics and hyphal differentiation during vegetative growth, hyphal

405

fusion, fruit-body development and sexual reproduction.49,50 Heat-shock protein 70 (Hsp70) chaperones

406

play important roles in stabilizing partially folded protein structures and rescuing the cells from apoptosis

407

during stress that affects a wide range of normal cellular physiology. 51 Our data showed that the activity

408

of superoxide dismutase and catalase were remarkably increased in CLS samples suggesting F. velutipes

19

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 20 of 28

409

mycelium adaption to increased oxidative stress induced by light stress. Superoxide dismutase and

410

catalase, located in the peroxisome, are two crucial enzymes for protecting cells and tissues from

411

oxidative damage by reactive oxygen species (ROS) in many fungi. 52 In addition, we observed that

412

expression levels of some initiation and elongation factors associated with the initiation of primordia

413

from mycelia and sexual development in fungi, such as translation elongation factor EF1-alpha53,

414

elongation factors8,54 and transcription factors55, were significantly up-regulated in CLS sample. More

415

broadly, understanding of expression pattern of the DEPs associated with survival, oxidative stress

416

adaption, autolysis and sexual development could be helpful to provide further insights into the response

417

mechanism of F. velutipes mycelia under CLS environment.

418

In conclusion, here, we try to systematically analyze and investigate the stress proteome associated

419

with mycelia in response to cold and light stress in F. velutipes by using iTRAQ labeling coupled to 2D

420

LC-MS/MS technique. Our data showed that three interesting up-regulated DEPs, including acetyl-

421

hydrolase, AASDH and ACPDC, may play an important role in improving the mycelium adaptability

422

and the ability to secrete pigment under CS environment in F. velutipes. The regulatory network

423

underlying the response mechanism of F. velutipes mycelia under CLS environment was complex and

424

multifaceted, as it included varied functions such as the rapid energy supply, the biosynthesis of lysine,

425

phenylalanine, tyrosine and GABA, calcium signal transduction process, dynein-dependent actin and

426

microtubule cytoskeletons formation, autolysis, oxidative stress adaption, tissue morphogenesis, and

427

other interesting stress-related processes. The data might provide further insights into the regulation

428

mechanism of mycelia in response to cold and light stress underlying mycelium growth, development

429

and fruit-body formation in F. velutipes. Furthermore, the data might also shed new light on an intuitive

430

understanding of the stress response mechanism in macro-fungi and provide valuable evidence for

20

ACS Paragon Plus Environment

Page 21 of 28

Journal of Agricultural and Food Chemistry

431

scientific improvement of some mushroom cultivation techniques in practice.

432

■ ASSOCIATED CONTENT

433

Table S1: the summarized primer data of 18 representative differentially expressed proteins.

434

Table S2: list of down-regulated proteins under both cold stress and cold and light stress in mycelia of

435

Flammulina velutipes

436

Table S3: list of down-regulated proteins under only cold stress in mycelia of Flammulina velutipes

437

Table S4: list of down-regulated proteins under only cold and light stress in mycelia of Flammulina

438

velutipes

439

■ AUTHOR INFORMATION

440

Corresponding Author

441

*Tel: +86 (0354) 6288325; fax: +86 (0354) 6288325; e-mail: [email protected], [email protected].

442

Author Contributions

443

All authors have given approval to the final version of the manuscript. Notes: The authors declare no

444

competing financial interest.

445

■ ACKNOWLEDGMENTS

446

This study was financially supported by grants from the National Natural Science Foundation of China

447

(no. 31301826), the Key Scientific and Technological Project from Shanxi Province (no. FT2014-03-01),

448

and the National Key Research and Development Program of China (no. 2017YFD0400200). We thank

449

MogoEdit (http://wwq.mogoedit.com) for its linguistic assistance during the preparation of this

450

manuscript.

451

■ REFERENCES

452

(1) Murat, C.; Mello, A.; Abbà, S.; Vizzini, A.; Bonfante, P. Edible mycorrhizal fungi: identification, life cycle and morphogenesis.

21

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

453 454 455 456 457 458 459 460 461

Page 22 of 28

Mycorrhiza 2008, 33, 707-732.

(2) Gasch, A. P.; Werner-Washburne M. The genomics of yeast responses to environmental stress and starvation. Funct. Integr.

Genomic. 2002, 2, 181–192.

(3) Idnurm, A.; Heitman, J. Light controls growth and development via conserved pathway in the fungal kingdom. PLoS Biol.

2005, 3, e95.

(4) Oliveira, J. M. P. F. D.; Graaff, L. H. D. Proteomics of industrial fungi: trends and insights for biotechnology. Appl. Microbiol.

Biotechnol. 2010, 89, 225-237

(5) Wu, Z.; Peng, W.; He X.; Wang, B.; Gan, B.; Zhang, X. Mushroom tumor: a new disease on Flammulina velutipes caused by

Ochrobactrum pseudogrignonense. FEMS Microbiol. Lett. 2016, 363: fnv226.

462

(6) Park, Y. J.; Baek, J. H.; Lee, S.; Kim, C.; Rhee, H.; Kim, H.; Seo, J. S.; Park, H. R.; Yoon, D. E.; Nam, J. Y.; Kim, H. I.; Kim,

463

J. G.; Yoon, H.; Kang, H. W.; Cho, J. Y.; Song, E. S.; Sung, G. H.; Yoo, Y. B.; Lee, C. S.; Lee, B. M.; Kong, W. S. Whole genome

464

and global gene expression analyses of the model mushroom Flammulina velutipes reveal a high capacity for lignocellulose

465

degradation, PLoS One 2014, 9, e93560.

466

(7) Yamada, M.; Kurano, M.; Inatomi, S.; Taguchi, G.; Okazaki, M.; Shimosaka, M. Isolation and characterization of a gene

467

coding for chitin deacetylase specifically expressed during fruiting body development in the basidiomycete Flammulina velutipes

468

and its expression in the yeast Pichia pastoris. FEMS Microbiol. Lett. 2008, 289, 130-137.

469

(8) Sekiya, S.; Yamada, M.; Shibata, K.; Okuhara, T.; Yoshida, M.; Inatomi, S.; Taguchi, G.; Shimosaka, M. Characterization of

470

a gene coding for a putative adenosine deaminase-related growth factor by RNA interference in the basidiomycete Flammulina

471

velutipes. J Biosci. Bioeng. 2013, 115, 360-365.

472 473 474

(9) Liu, J. Y.; Chang, M. C.; Meng, J. L.; Feng, C. P.; Zhao H., Zhang M. L. Comparative proteome reveals metabolic changes

during the fruiting process in Flammulina velutipes. J Agr. Food Chem. 2017, 65, 5091-5100.

(10) Sakamoto, Y.; Ando, A.; Tamai, Y.; Yajima, T. Influences of temperature and light on morphological changes during fruit

22

ACS Paragon Plus Environment

Page 23 of 28

475 476 477 478 479 480 481

Journal of Agricultural and Food Chemistry

body formation in Flammulina velutipes, Mycoscience 2004, 45, 333-339.

(11) Sakamoto, Y.; Ando, A.; Tamai, Y.; Yajima, T. Pileus differentiation and pileus-specific protein expression in Flammulina

velutipes, Fungal Genet. Biol. 2007, 44, 14-24.

(12) Liu, J. Y.; Chang, M. C.; Meng, J. L.; Feng, C. P.; Yuan, L. G. iTRAQ-based quantitative proteome revealed metabolic

changes of Flammulina velutipes mycelia in response to cold stress. J Proteomics 2017, 156, 75-84.

(13) Wang, W.; Liu, F.; Jiang, Y.; Wu, G.; Guo, L.; Chen, R.; Chen, B.; Lu, Y.; Dai, Y.; Xie, B. The multigene family of fungal

laccases and their expression in the white rot basidiomycete Flammulina velutipes. Gene 2015, 563, 142-149.

482

(14) Patterson, S. D.; Aebersold, R. H. Proteomics: the first decade and beyond. Nat. Genet. 2003, 33, 311-323.

483

(15) Thelen, J. J.; Miernyk, J. A. The proteomic future: where mass spectrometry should be taking us. Biochem. J 2012, 444, 169-

484 485 486 487 488 489 490 491 492 493 494

181.

(16) Zieske, L R. A perspective on the use of iTRAQTM reagent technology for protein complex and profiling studies. J Exp. Bot.

2006, 57, 1501-1508

(17) Nesvizhskii, A. I.; Keller, A.; Kolker, E.; Aebersold, R. A statistical model for identifying proteins by tandem mass

spectrometry. Anal. Chem. 2003, 75, 4646-4658.

(18) Shadforth, I. P.; Dunkley, T. P.; Lilley, K. S.; Bessant, C. i-Tracker: For quantitative proteomics using iTRAQ™. BMC

Genomics 2005, 6, 145.

(19) Ma, C.; Zhou, J.; Chen, G.; Bian, Y.; Lv, D.; Li, X; Wang, Z.; Yan, Y. iTRAQ-based quantitative proteome and phosphoprotein

characterization reveals the central metabolism changes involved in wheat grain development. BMC Genomics 2014, 15, 1029.

(20) Kültz, D. Evolution of the cellular stress proteome: from monophyletic origin to ubiquitous function. J Exp. Biol. 2003, 206,

3119-3124.

495

(21) Kim, Y.; Nandakumar, M. P.; Marten, M. R. Proteomics of filamentous fungi. Trends Biotechnol. 2007, 25, 395-400.

496

(22) Strijbis, K.; Distel, B. Intracellular acetyl unit transport in fungal carbon metabolism. Eukaryot Cell 2010, 9, 1809-1815.

23

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

497 498 499 500 501 502 503 504

Page 24 of 28

(23) Lambie, S. C. Lipid and acetate metabolism influence host colonization by the fungal plant pathogen Ustilago maydis.

Thesis/Dissertation; University of British Columbia 2014; doi: 10.14288/1.0135647.

(24) Estrada, A. F.; Youssar, L.; Scherzinger, D.; Al-Babili, S.; Avalos, J. The ylo-1 gene encodes an aldehyde dehydrogenase

responsible for the last reaction in the Neurospora carotenoid pathway. Mol. Microbiol. 2008, 69, 1207-1220.

(25) Jitrapakdee, S.; Maurice, M. S.; Rayment, I.; Cleland, W. W.; Wallace, J. C.; Attwood, P. V. Structure, mechanism and

regulation of pyruvate carboxylase. Biochem. J 2008, 413, 369-387.

(26) Bräsen, C.; Esser, D.; Rauch, B.; Siebers, B. Carbohydrate metabolism in archaea: current insights into unusual enzymes and

pathways and their regulation. Microbiol. Mol. Biol. Rev. 2014, 78, 89-175.

505

(27) Zabriskie, T. M.; Jackson, M. D. Lysine biosynthesis and metabolism in fungi. Nat. Prod. Rep. 2000, 17, 85-97.

506

(28) Cerqueira, N. M. F. S. A.; Fernandes, P. A.; Ramos, M. J. Computational mechanistic studies addressed to the transimination

507

reaction present in all pyridoxal 5-phosphate-requiring enzymes. J Chem. Theory Comput. 2011, 7, 1356-1368.

508

(29) Yang, D. D.; de Billerbeck, G. M.; Zhang, J. J.; Rosenzweig, F.; Francois, J. M. Deciphering the origin, evolution and

509

physiological function of the subtelomeric aryl-alcohol dehydrogenase gene family in the yeast Saccharomyces cerevisiae. Appl.

510

Environ. Microbiol. 2018, 84, e01553-17.

511 512 513 514

(30) Hosseini, N. M.; Hussain, M. A.; Britz, M. L. Stress responses in probiotic Lactobacillus casei. Crit. Rev. Food Sci. 2015,

55, 740-749.

(31) Appaji Rao N, Ambili M, Jala VR, Subramanya HS, Savithri HS. Structure-function relationship in serine

hydroxymethyltransferase. Biochim. Biophys. Acta. 2003, 1647, 24-29.

515

(32) Ramesh, S. A.; Tyerman, S. D.; Xu, B.; Bose, J.; Kaur, S.; Conn, V.; Domingos, P.; Ullah, S.; Wege, S.; Shabala, S.; Feijo J.

516

A.; Ryan, P. R.; Gilliha, M. GABA signalling modulates plant growth by directly regulating the activity of plant-specific anion

517

transporters. Nat. Commun. 2015, 6, 7879.

518

(33) Baldy, P. Metabolisme du gamma-aminobutyrate chez Agaricus bisporus. I. la L-glutamate-1-carboxy-lyase. Physiol.

24

ACS Paragon Plus Environment

Page 25 of 28

519 520 521 522 523

Journal of Agricultural and Food Chemistry

Plantarum 2010, 34, 365-372.

(34) Mead, O; Thynne, E.; Winterberg, B; Solomon, P. S. Characterising the role of GABA and its metabolism in the wheat

pathogen Stagonospora nodorum. PLoS One 2013, 8, e78368.

(35) Chen, S. Y.; Ho, K. J.; Hsieh, Y. J.; Wang, L. T.; Mau, J. L. Contents of lovastatin, γ-aminobutyric acid and ergothioneine in

mushroom fruiting bodies and mycelia. LWT-Food Sci. Technol. 2012, 47, 274-278.

524

(36) Veselá, A. B.; Rucká, L.; Kaplan, O.; Pelantová, H.; Nešvera, J.; Pátek, M.; Martínková, L. Bringing nitrilase sequences from

525

databases to life: the search for novel substrate specificities with a focus on dinitriles. Appl. Microbiol. Biotechnol. 2016, 100,

526

2193-2202.

527 528 529 530

(37) Kao, C. M.; Liu, J. K.; Lou, H. R.; Lin, C. S.; Chen, S. C. Biotransformation of cyanide to methane and ammonia by Klebsiella

oxytoca. Chemosphere 2003, 50, 1055-1061.

(38) Filamentous fungus Aspergillus oryzae has two types of α-1, 2-mannosidases, one of which is a microsomal enzyme that

removes a single mannose residue from Man9GlcNAc2. Glycoconjugate J 2000, 17, 745-748.

531

(39) Zhang, S.; Zheng, H.; Chen, Q.; Chen, Y.; Wang, S.; Lu, L.; Zhang S. The lectin chaperone calnexin is involved in

532

endoplasmic reticulum stress response by regulating Ca2+ homeostasis in Aspergillus nidulans. Appl. Environ. Microbiol. 2017, 83,

533

e00673-17.

534

(40) Rispail, N.; Soanes, D. M.; Ant, C.; Czajkowski, R.; Grünler, A.; Huguet, R.; Perez-Nadales, E.; Poli, A.; Sartorel, E.;

535

Valiante, V.; Yang, M.; Beffa, R.; Brakhage, A. A.; Gow, N. A.; Kahmann, R.; Lebrun, M. H.; Lenasi, H.; Perez-Martin, J.; Talbot,

536

N. J.; Wendland, J.; Di Pietro A. Comparative genomics of MAP kinase and calcium–calcineurin signaling components in plant

537

and human pathogenic fungi. Fungal Genet. Biol. 2009, 46, 287-298.

538 539 540

(41) Adamí ková, L.; Straube, A.; Schulz, I.; Steinberg, G. Calcium signaling is involved in dynein-dependent microtubule

organization. Mol. Biol. Cell 2004, 15, 1969-1980.

(42) Bridges, A. A.; Gladfelter, A. S. Septin form and function at the cell cortex. J. Biol. Chem. 2015, 290, 17173-17180.

25

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

541 542 543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561 562

Page 26 of 28

(43) Gunning P. W.; Ghoshdastider, U.; Whitaker, S.; Popp, D.; Robinson, R. C. The evolution of compositionally and

functionally distinct actin filaments. J Cell Sci. 2015, 128: 2009-2019.

(44) Takeshita, N.; Manck, R.; Grun, N.; de Vega, S. H.; Fischer, R. Interdependence of the actin and the microtubule

cytoskeleton during fungal growth. Curr. Opin. Microbiol. 2014, 20, 34-41.

(45) Andrea, V. E.; Inna, G; Míguez, D. G.; Acaimo, G. R.; Martín-Bermudo, M. D. Myosin light-chain phosphatase regulates

basal actomyosin oscillations during morphogenesis. Nat. Commun. 2016, 7, 10746.

(46) Chang, J.; Tower, R. J.; Lancaster, D. L.; Rachubinski, R. A. Dynein light chain interaction with the peroxisomal import

docking complex modulates peroxisome biogenesis in yeast. J Cell Sci. 2013, 126: 4698-4706.

(47) Egan, M. J.; McClintock, M. A.; Hollyer, I. H. L.; Elliott, H. L.; Reck-Peterson, S. L. Cytoplasmic dynein is required for

the spatial organization of protein aggregates in filamentous fungi. Cell Rep. 2015, 11, 201-209.

(48) Zhou, Y.; Kang, L.; Niu, X.; Wang, J.; Liu, Z.; Yuan, S. Purification, characterization and physiological significance of a

chitinase from the pilei of Coprinopsis cinerea fruiting bodies, FEMS Microbiol. Lett. 2016, 363, fnw120.

(49) Corral-Ramos, C.; Roca M. G.; Pietro, A. D.; Roncero, M. I. G.; Ruiz-Roldán, C. Autophagy contributes to regulation of

nuclear dynamics during vegetative growth and hyphal fusion in Fusarium oxysporum. Autophagy 2015, 11, 131-144.

(50) Voigt, O.; Poggeler, S. Autophagy genes Smatg8 and Smatg4 are required for fruiting-body development, vegetative growth

and ascospore germination in the filamentous ascomycete Sordaria macrospora. Autophagy 2013, 9, 33-49

(51) Mashaghi, A.; Bezrukavnikov, S.; Minde, D. P.; Wentink, A. S.; Kityk, R.; Zachmann-Brand, B.; Mayer, M. P.; Kramer, G.;

Bukau, B.; Tans, S. J. Alternative modes of client binding enable functional plasticity of Hsp70. Nature 2016, 539, 448-453.

(52) Isobe, K.; Inoue, N.; Takamatsu, Y.; Kamada, K.; Wakao, N. Production of catalase by fungi growing at low pH and high

temperature, J Biosci. Bioeng. 2006, 101, 73-76.

(53) Chen, L.; Zhang, B. B.; Cheung, P. C. K. Comparative proteomic analysis of mushroom cell wall proteins among the

different developmental stages of Pleurotus tuber-regium. J Agric. Food Chem. 2012, 60, 6173-6182.

26

ACS Paragon Plus Environment

Page 27 of 28

563 564 565 566

Journal of Agricultural and Food Chemistry

(54) Morrissey, C.; Schwefel, D.; Ennis-Adeniran, V.; Taylor, I. A.; Crow, Y. J.; Webb, M. The eukaryotic elongation factor

eEF1A1 interacts with SAMHD1. Biochem. J. 2015, 466, 69-76.

(55) Ohm, R. A.; de Jong, J. F.; de Bekker, C.; Wosten, H. A.; Lugones, L. G. Transcription factor genes of Schizophyllum

commune involved in regulation of mushroom formation. Mol. Microbiol. 2011, 81, 1433-1445.

27

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Table of Contents Graphic:

ACS Paragon Plus Environment

Page 28 of 28