Molecular Characterization of WRKY Transcription Factors that Act as

Oct 16, 2018 - Molecular Characterization of WRKY Transcription Factors that Act as Negative Regulators of O-methylated Catechins Biosynthesis in Tea ...
0 downloads 0 Views 3MB Size
Subscriber access provided by UNIV OF LOUISIANA

Agricultural and Environmental Chemistry

Molecular Characterization of WRKY Transcription Factors that Act as Negative Regulators of O-methylated Catechins Biosynthesis in Tea Plant (Camellia Sinensis L.) Yong Luo, Shuangshuang Yu, Juan Li, Qin Li, Kunbo Wang, Jianan Huang, and Zhonghua Liu J. Agric. Food Chem., Just Accepted Manuscript • Publication Date (Web): 16 Oct 2018 Downloaded from http://pubs.acs.org on October 16, 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 26

Journal of Agricultural and Food Chemistry

1

Molecular Characterization of WRKY Transcription Factors that Act as

2

Negative Regulators of O-methylated Catechins Biosynthesis in Tea Plant

3

(Camellia Sinensis L.)

4

Authors: Yong Luo1, Shuangshuang Yu1, Juan Li1,2, Qin Li1,2, Kunbo Wang1,2*, Jianan Huang1,2,*, Zhonghua

5

Liu1,2,*

6 7

1

8

China

9

2

Key Laboratory of Tea Science of Ministry of Education, Hunan Agricultural University, Changsha, 410128, P. R.

National Research Center of Engineering and Technology for Utilization of Botanical Functional Ingredients &

10

Hunan Co-innovation Center for Utilization of Botanical Functional Ingredients, Hunan Agricultural University,

11

Changsha, 410128, P. R. China

12 13 14 15 16 17 18 19 20 21 22 23 1

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

24 25

ABSTRACT: Tea O-methylated catechins especially (-)-epigallocatechin-3-O-(3-O-methyl)-gallate

26

(EGCG3"Me) have been attracted much attention due to their positive health effects. The

27

transcription regulators of O-methylated catechins biosynthesis remains elusive. In this study,

28

expression pattern of genes related to O-methylated catechins biosynthesis including CsLAR, CsANS,

29

CsDFR, CsANR and CCoAOMT in three tea cultivars with different content of EGCG3"Me was

30

investigated. Two WRKY transcription factors (TFs), designated as CsWRKY31 and CsWRKY48,

31

belonging to group IIc and IIb of the WRKY family respectively, were further identified.

32

CsWRKY31 and CsWRKY48 were nuclear-localized proteins, and possessed transcriptional

33

repression ability. Furthermore, expression of CsWRKY31 and CsWRKY48 showed negative

34

correlation with CsLAR, CsDFR, and CCoAOMT, during EGCG3"Me accumulation in tea leaves.

35

More importantly, W-box (C/T)TGAC(T/C) elements were located in the promoter of CsLAR,

36

CsDFR, and CCoAOMT, and further assays revealed that CsWRKY31 and CsWRKY48 were

37

capable of repressing the transcription of CsLAR, CsDFR, and CCoAOMT, via the attachment of

38

their promoters to the W-box elements. Collectively, our findings identify two novel negative

39

regulators of O-methylated catechins biosynthesis in tea plant, which might provide a potential

40

strategy to breed high quality tea cultivar.

41

KEYWORDS: Camellia sinensis, Catechins biosynthesis, EGCG3"Me, WRKY, Transcriptional regulator

42 43 44 2

ACS Paragon Plus Environment

Page 2 of 26

Page 3 of 26

Journal of Agricultural and Food Chemistry

45 46

INTRODUCTION

47

The tea plant (Camellia Sinensis L.) is one of the most important commercial plants. It has been cultivated

48

worldwide for at least 2000 years due to it stimulating and soothing effects as well as positive effects on human

49

health. Numerous studies showed that tea can prevent diseases including anti-cancer, anti-cardiovascular,

50

anti-neurodegenerative, and oxidative stress1-3. In general, the bioactive properties of tea are catechins. A new

51

catechin, (-)-epigallocatechin-3-O-(3-O-methyl)-gallate (EGCG3"Me) that has been reported recently to have

52

stronger anti-allergic and anti-hypertensive effects than epigallocatechin-3-O-gallate (EGCG)4-6. Another

53

biologically important characteristic of EGCG3"Me is absorbed more easily than EGCG4-6. The quality and taste of

54

tea is predominantly determined by flavonoids and catechins. Previously, we have identified two tea cultivars

55

‘Jinmudan’ and ‘Jinguanyin’ with high EGCG3"Me content7, but the molecular mechanism(s) concerning

56

EGCG3"Me biosynthesis in these two cultivars remain unclear. Therefore, exploring the regulatory mechanisms of

57

catechins biosynthesis by diverse factors may contribute to greater understanding of mechanisms for increasing

58

catechins accumulation in tea plants by effective strategies such as environmental regulation and genetic breeding.

59

In general, catechins and EGCG3"Me are thought to be biosynthesized through the flavonoid pathways5, 8-9. As

60

shown in Fig. 1a, the biosynthesis of catechins and EGCG3"Me is regulated by the dynamic balance of the

61

expression of biosynthesis-related genes leucocyanidin reductase (LAR), anthocyanidin synthase (ANS),

62

dihydroflavonol-4-reductase(DFR)

63

3-O-methyltransferase (CCoAOMT) has been suggested to be directly associated with EGCG3"Me accumulation5.

and

anthocyanidin

reductase

(ANR).

Moreover,

caffeoyl-CoA

64

Recently, many transcription regulators including MYB, bHLH, MADS and WD40 were found to be associated

65

with catechins metabolic regulation by controlling the expression of important enzyme genes participate in 3

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

66

catechins biosynthetic pathway10-11. As a major transcription factor(TF) familys in plants, WRKYs are key

67

regulators of plant development, and stress defenses12-14. The WRKY family proteins have one or two WRKY

68

domains consisting of a 60-amino-acid region that possesses a highly conserved WRKYGQK in their N-termini,

69

and C2HC zinc finger in C-termini12. WRKY transcription factors can regulate downstream transcription by

70

specifically binding to the W-box (C/T)TGAC(T/C)15. Upon the completion of genome sequencing, WRKY TFs

71

have been isolated and identified in many plants, including Arabidopsis, Glycine max, Oryza sativa and Hordeum

72

vulgare16-18. WRKY TFs are also reported to play regulatory roles in plant metabolite biosynthesis, including

73

phenolic, phenylpropanoids, terpenes and alkaloids by regulating the downstream genes19. For example, a opium

74

poppy (Papaver somniferum) WRKY was found to affect the biosynthesis of alkaloids20. Similarly, OsWRKY76

75

was reported to be regulating the production of terpenes and phenylpropanoids21. WRKYs were reported to

76

up-regulate the biosynthesis of phenylpropanoids, such as flavonols, but few down-regulated regulators was

77

reported from various plants22-23. Knowledge regarding the involvement of WRKY TFs in EGCG3"Me

78

biosynthesis is far from complete.

79

In this study, two tea cultivars Jinmudan and Jinguanyin with high levels of EGCG3"Me were used. The

80

expression of EGCG3"Me biosynthesis-related genes including CsLAR, CsDFR and CCoAOMT were performed.

81

Moreover, two potential WRKY TFs termed CsWRKY31 and CsWRKY48 that negatively associated with

82

EGCG3"Me biosynthesis were identified and characterized. Further experiments revealed that CsWRKY31 and

83

CsWRKY48 directly bound to the promoter of CsLAR, CsDFR and CCoAOMT, and repressed their expressions.

84

Our findings demonstrate that CsWRKY31 and CsWRKY48 might negatively participate in catechins biosynthesis

85

by down-regulating CsLAR, CsDFR and CCoAOMT expressions. Our study disseminate novel knowledge about the

86

transcriptional regulatory of plant catechins biosynthesis and clues for the breeding of tea cultivar resources with 4

ACS Paragon Plus Environment

Page 4 of 26

Page 5 of 26

Journal of Agricultural and Food Chemistry

87

high EGCG3"Me levels.

88

MATERIALS AND METHODS

89

Plant materials and treatments.

90

Three different tea cultivars (Camellia sinensis (L.) O. Kuntze.) including Fudingdabai, Jinmudan and

91

Jinguanyin, were cultivated in the Gao qiao tea farm of the Tea Research Institute of Hunan Agricultural Sciences,

92

Changsha, China. For each cultivar, leaves at standard stage with “one bud and two leaves” from ten different tea

93

plants were collected (Figure 1b), then frozen in liquid nitrogen and stored at −80 °C until later use for

94

physiological and molecular analysis.

95

Quantification of EGCG and EGCG3"Me.

96

The EGCG3"Me contents in tea leaves in the different three cultivars were detected by LC-10ATVP HPLC

97

system (Shimadzu, Tokyo, Japan). The tea samples were separated using a reversed-phase column (Welchorm C18

98

200×4.6, 5mm). Standards of EGCG and EGCG3"Me(cas:83104-874) were purchased from Sigma-Aldrich (St.

99

Louis, MO, USA).

100

Gene Isolation, Sequence and Expression Analysis.

101

Each tea leaf sample was used for the extraction of total RNA using the RNeasy Mini kit (Qiagen, Hilden,

102

Germany) according to the manufacturer’s instructions. RNA quality and quantity were monitored by gel

103

electrophoresis and spectrophotometry. The DNA-free total RNA was reversely transcribed into the first strand of

104

cDNAs using a Prime ScriptRT Reagent Kit (Takara, Dalian, China). Three genes including CsLAR, CsDFR and

105

CCoAOMT related to catechins and EGCG3"Me biosynthesis and two WRKY genes CsWRKY31 and CsWRKY48

106

were identified and selected from the Camellia sinensis genome (http://www.plantkingdomgdb.com/tea_tree/) and

107

our RNA-seq database by using leaves from different cultivars (unpublished data). The full-length sequences of 5

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

108

CsWRKY31 and CsWRKY48 were amplified with specific primers (Primers are shown in Table S1), and blasted in

109

the National Center for Biotechnology Information (NCBI). The sequence alignment and phylogenetic tree were

110

constructed using ClustalW and MEGA5, respectively.

111

Gene sequences we used in this article were downloaded from NCBI, including DFR (GenBank Accession

112

No. AY648027.1), ANS (GenBank Accession No. AY830416.1), ANR (GenBank Accession No. GU944768.1),

113

LAR(GenBank Accession No. EF205148.1), CCoAOMT(GenBank Accession No. KF268598.1). Gene expression

114

level was performed using Go Taq qPCR Master Mix Kit (Promega, Madison, WI, USA) on a Bio-Rad CFX96

115

real-time PCR system. The PCR conditions started with an pre-denaturation step at 94°C for 5 min, then followed

116

by 40 cycles of 94°C for 10s, 60°C for 30s, and 72°C for 30s, in a 20-mL reaction. Primers are shown in Table S1.

117

Subcellular Localization Analysis.

118

The coding region of CsWRKY31 or CsWRKY48 was cloned into the pEAQ-GFP vector (Primers used are shown

119

in Table S1)and was verified by further sequencing. About 4 to 6 weeks-old N. benthamiana leaves were transiently

120

transformed with the well constructed vectors through Agrobacterium tumefaciens strain GV3101 as described

121

previously24. GFP fluorescence was recorded by a fluorescence microscope after 2-3 days of infiltration.

122

Transcriptional Activation Assay.

123

CsWRKY31 and CsWRKY48 were inserted into the pGBKT7 vector (Clontech, USA). Then the fusion construct

124

pGBKT7-CsWRKY31 or -CsWRKY48, positive control (p53+T-antigen) and negative control (pGBKT7 vector)

125

were transformed separately into yeast cells. The transcriptional activation ability of CsWRKY31 and CsWRKY48

126

was evaluated basing on the growth status and α-galactosidase activity of yeast cells that grow on SD medium

127

(SD/-Trp), or (SD/-Trp-His-Ade).

128

Promoter Isolation and Analysis. 6

ACS Paragon Plus Environment

Page 6 of 26

Page 7 of 26

Journal of Agricultural and Food Chemistry

129

The promoter sequences of CsLAR, CsDFR and CCoAOMT were amplified by PCR (primers are shown in Table

130

S1) using genomic DNA of tea leaves as the template according to the DNeasy Plant Mini Kit (Qiagen). Conserved

131

W-box motifs presented in the promoters were predicted and identified through the Plant-CARE online software

132

(http://bioinformatics.psb.ugent.be/webtools/plantcare/html/).

133

Protein Expression and Electrophoretic Mobility Shift Assay (EMSA).

134

The sequence of CsWRKY31 N-terminus (from 475 to 753 bp including the WRKY domain) and the coding

135

region of CsWRKY48 were cloned into the pGEX-4T-1 (Amersham Biosciences) vector and transformed into

136

Escherichia coli strain BM Rosetta (DE3). The induced GST-CsWRKY31-N and GST-CsWRKY48 protein was

137

induced and purified using a Glutathione-Superflow Resin (Clontech) and stored at −80 °C. The 5’ ends of

138

synthesized oligonucleotide probes containing W-box in the promoters of CsLAR, CsDFR and CCoAOMT were

139

labeled with biotin. The Light Shift Chemiluminescent EMSA Kit (Thermo Scientific) was used for EMSA

140

experiment as previously described25. Briefly, GST-CsWRKY31-N and GST-CsWRKY48 protein and

141

biotin-labeled probes were incubated together, then the assay mixtures were analyzed by 6% native polyacrylamide

142

gel electrophoresis, after transferred onto nylon membrane, the protein-DNA complexes was detected using a

143

ChemiDoc™ MP Imaging System (Bio-Rad, USA). GST protein alone as well as the unlabeled and mutated probes

144

were used as negative control.

145

Dual Luciferase Reporter Assay.

146

To investigate the transcriptional activity of CsWRKY31 and CsWRKY48, their coding sequences were

147

cloned into pBD vector. The double-reporter vector included a GAL4-LUC and an internal control REN as

148

previously described25. To assess the binding activity of CsWRKY31 and CsWRKY48 to the promoters of CsLAR,

149

CsDFR and CCoAOMT. CsWRKY31 and CsWRKY48 were cloned into pEAQ vector as effectors. The promoters 7

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 8 of 26

150

of CsLAR, CsDFR and CCoAOMT were recombined to pGreenII 0800-LUC vector as reporters. The constructed

151

effectors and reporters plasmids with different combinations were co-transformed into the Nicotiana benthamiana

152

leaves .

153

Firefly luciferase and Renilla luciferase were quantified after 2 d of infiltration with the dual luciferase assay

154

kit (Promega, USA). The transactivation ability of CsWRKY31 and CsWRKY48 was assessed by the LUC to REN

155

ratio. At least six independent experiments were included for each combination.

156

Statistical Analyses.

157

All data are presented as means ± standard errors (S.E.) of three or six independent biological replicates. To

158

compare the statistical difference and determine the significance of experimental means, the data were evaluated by

159

student's t-test at P value < 0.05 or 0.01.

160

RESULTS AND DISCUSSION

161

The

162

Biosynthesis-Related Genes among Cultivars

Variation

of

EGCG

and

EGCG3"Me

Contents,

and

the

Expression

of

EGCG3"Me

163

As a bioactive constituent, EGCG3"Me was recently reported to occur naturally in tea leaves, but the quantities

164

of EGCG3"Me was extremely limited and detected in few tea cultivars5. In this study, EGCG and EGCG3"Me

165

levels in tea leaves of three cultivars including ‘Fudingdabai’, ‘Jinmudan’ and ‘Jinguanyin’ at the standard stage

166

with “one bud and two young leaves” were measured by HPLC. As shown in Figure 2a, EGCG3"Me content in the

167

leaves of ‘Jinguanyin’ and ‘Jinmudan’ was higher than that in ‘Fudingdabai’. EGCG3"Me level in ‘Jinguanyin’ and

168

‘Jinmudan’ was 6.12 ± 0.35 mg/g and 7.81 ± 0.30 mg/g respectively, whereas was almost undetectable in

169

‘Fudingdabai’. The contents of EGCG in leaves of ‘Jinguanyin’ and ‘Jinmudan’ was also higher than in

170

‘Fudingdabai’ (Figure 2b). Lv et al. (2014)26 reported that different tea cultivars contained different EGCG3"Me 8

ACS Paragon Plus Environment

Page 9 of 26

171

Journal of Agricultural and Food Chemistry

content and only 4 tea cultivars were found to be rich in EGCG3"Me.

172

LAR, ANS, DFR, ANR and CCoAOMT are important enzymes in catechins and EGCG3"Me biosynthesis

173

(Figure 1a). To further investigate the association of LAR, ANS, DFR, ANR and CCoAOMT with EGCG3"Me

174

biosynthesis, their gene expression levels in leaves of ‘Fudingdabai’, ‘Jinmudan’ and ‘Jinguanyin’ were compared.

175

Some significant correlations were noted. As shown in Figure 2c, transcript levels of CsLAR, CsDFR and

176

CCoAOMT were obviously higher in ‘Jinguanyin’ and ‘Jinmudan’ than in ‘Fudingdabai’, while no obvious

177

difference among these three cultivars was observed in CsANR and CsANS expression.

178

Identification and Bioinformatics Analysis of CsWRKY31 and CsWRKY48.

179

To date, ~72, 109, 197 and 45 WRKY members have been identified in Arabidopsis, Oryza sativa, Glycine max

180

and Hordeum vulgare respectively16, 18-19, 27. About 50 WRKY genes have been identified in tea plants28, and two

181

WRKY genes attracted our attention as their expression was obviously down-regulated in ‘Jinguanyin’ and

182

‘Jinmudan’ according to our RNA-seq database. Therefore, the full-length of these two genes were cloned and

183

designated as CsWRKY31 and CsWRKY48 since the sequence of these two WRKY genes had high degree of

184

homology with AtWRKY31 (NP_567644.1) (58%) and AtWRKY48 (NP_199763.1) (54%) respectively.

185

CsWRKY31 and CsWRKY48 cDNA contains an Open Reading Frame (ORF) of 1260 and 486 bp in length,

186

encoding 419 and 161 amino acids that resulted in calculated molecular mass of 45.13 and 17.87 kDa, respectively.

187

WRKY proteins can be categorized into three major groups (I-III) according to the number of WRKY DNA

188

binding domain and Zinc-finger motif. And group II can be further divided into IIa, IIb, IIc, IId, and IIe12, 29.

189

Multiple sequence alignment showed that CsWRKY31 and CsWRKY48 consisted of a highly conserved WRKY

190

DNA binding domain (DBD) and a zinc finger (C2H2 motif) (Figure 3a), which shared a high degree of homology

191

with the group II WRKY family. Phylogenetic analysis further indicated that CsWRKY31 and CsWRKY48 was 9

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 10 of 26

192

clustered with AtWRKY31 and AtWRKY48 into the group IIb and group IIc WRKY family, respectively (Figure

193

3b). The group IIb Solanum lycopersicum SlWRKY73 and group IIc AtWRKY23 have been reported to regulate

194

plant-specialized metabolism, such as the biosynthesis of important pharmaceutical, aromatherapy, biofuel, and

195

industrial components30-31. Similarly, pitaya HpWRKY44 could be clustered with the Arabidopsis thaliana Group I

196

WRKYs including AtWRKY44 and AtWRKY33. Notably, the Arabidopsis thaliana WRKY44 and WRKY33 were

197

reported to regulate the biosynthesis of indole alkaloid and phenylpropanoid. Accordingly, HpWRKY44 was

198

reported to be involved in the biosynthesis of secondary metabolites in pitaya24. Thus, it could be speculated that

199

CsWRKY31 and CsWRKY48 might be associated with the biosynthesis of EGCG3"Me in tea plants.

200

Expression and Molecular Characterization of CsWRKY31 and CsWRKY48.

201

Similar to the RNA-seq data, RT-qPCR analysis exhibited that CsWRKY31 and CsWRKY48 transcript levels

202

were remarkably lower in ‘Jinguanyin’ and ‘Jinmudan’ than in ‘Fudingdabai’ (Figure 4a), showing the opposite

203

pattern with CsLAR, CsDFR and CCoAOMT expression, and EGCG3"Me accumulation, implied that CsWRKY31

204

and CsWRKY48 might be negatively involved in EGCG3"Me biosynthesis in the tea plant. WRKYs are usually

205

nuclear proteins27,

206

GFP-empty proteins with these two WRKY proteins transiently expressed into tobacco leaves. As shown in Figure

207

4b, the 35S:GFP-empty signal was found in the nucleus and cytoplasm, whereas the 35S:CsWRKY31-GFP and

208

35S:CsWRKY48-GFP fusion proteins were predominately observed in the nucleus, indicating that CsWRKY31

209

and CsWRKY48 are nucleus-localized proteins.

32-33.

To find out the subcellular location of CsWRKY31 and CsWRKY48, we fused the

210

Transcriptional ability of TFs has a key role in regulating downstream genes34-35. The transcriptional activities

211

of CsWRKY31 and CsWRKY48 were further investigated using a GAL4-responsive reporter system in yeast cells.

212

The positive control (pGBKT7-53 + pGADT7-T) had normal growth in the SD plates without tryptophan, histidine, 10

ACS Paragon Plus Environment

Page 11 of 26

Journal of Agricultural and Food Chemistry

213

and adenine and showed α-galactosidase activity, whereas the yeasts transformed with the pGBKT7-empty vector,

214

pGBKT7-CsWRKY31 and pGBKT7-CsWRKY48 did not grow and lacked the α-galactosidase activity (Figure 5a),

215

indicating that CsWRKY31 and CsWRKY48 did not possess any transcriptional activation activity in the yeast

216

cells, and they might be transcriptional repressors. The transcriptional repression activities of CsWRKY31 and

217

CsWRKY48 were further verified by dual-luciferase assays (Figure 5b and 5c). Compared with the pBD-empty

218

control, both CsWRKY31 and CsWRKY48 significantly repressed the LUC reporter activities (Figure 5c),

219

suggesting that CsWRKY31 and CsWRKY48 are likely act as transcriptional repressors.

220

CsWRKY31 and CsWRKY48 Target CsLAR, CsDFR and CCoAOMT Promoters via the W-box Element.

221

It is well-known that WRKYs regulate their target genes expressions by binding to consensus sequence

222

(C/T)TGAC(T/C), known as W-box presented in the promoters15,

223

CsLAR, CsDFR and CCoAOMT promoters (Text S1). EMSA using the purified recombinant GST-CsWRKY31-N

224

and GST-CsWRKY48 protein (Figure 6a and 6c) was performed to confirm whether CsWRKY31 and CsWRKY48

225

could target CsLAR, CsDFR and CCoAOMT promoters via the W-box. As expected, both GST-CsWRKY31-N and

226

CsWRKY48 fusion proteins could directly bind to labeled CsLAR, CsDFR and CCoAOMT fragment containing

227

W-box and caused mobility shifts. The mobility shift was effectively abolished due to unlabeled fragment being

228

used as a cold probe in a dose-dependent manner, but not by the mutated probes (Figure 6b and 6d). In addition,

229

when the probes was incubated with GST alone, the mobility shift was not observed. These data reveal that

230

CsWRKY31 and CsWRKY48 target the W-box motif in the CsLAR, CsDFR and CCoAOMT promoters.

231

Trans-repression of CsWRKY31 and CsWRKY48 on CsLAR, CsDFR and CCoAOMT Promoters

19.

Indeed, the W-box motifs were found in

232

To further investigate whether the promoter activity of CsLAR, CsDFR and CCoAOMT could be repressed by

233

CsWRKY31 or CsWRKY48, the well-established transient dual luciferase assay was employed, using a double 11

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 12 of 26

234

reporter plasmid containing the LUC driven by CsLAR, CsDFR and CCoAOMT promoters and the REN driven by

235

the CaMV35S promoter, together with an effector plasmid expressing CsWRKY31 or CsWRKY48 (Figure 7a).

236

The results showed that LUC/REN ratio was remarkably decreased when the CsLAR, CsDFR and CCoAOMT

237

pro-LUC reporter construct was co-transfected with CsWRKY31 or CsWRKY48. Compared with the empty

238

control (Figure 7b), suggesting that CsWRKY31 and CsWRKY48 repressed the promoter activity of CsLAR,

239

CsDFR and CCoAOMT. Taken together, our data support the notion that CsWRKY31 and CsWRKY48 act as

240

transcriptional repressors of CsLAR, CsDFR and CCoAOMT through direct binding to their promoters. In tea plant,

241

the R2R3-MYB, bHLH and WD40 were reported to be involved in flavonoid biosynthesis10,36. For instance,

242

MYB4a, a R2R3-MYB TF, was reported negatively regulate the phenylpropanoid and shikimate pathways37.

243

WRKY TFs were mostly reported in plant responses to stress. For example, Arabidopsis WRKY52 was able to

244

confer resistance against bacterial pathogen Ralstonia solanacearum38. In tea plants, WRKY2 is shown to function

245

in cold and drought stress responses39. Nevertheless, several studies suggest that WRKY TFs regulate the

246

production of several secondary metabolites by regulating the genes within the metabolite pathway. For instance,

247

Arabidopsis WRKY12, WRKY23 and WRKY44 played an important role in regulating the production of lignin,

248

flavonol and tannin, respectively19, 30, 40-41. While Rice (Oryza sativa) OsWRKY45 up-regulated the accumulation

249

of oryzalexin, phytocassane and momilactone, by priming the target gene expression22. Overexpression of

250

OsWRKY13 up-regulates phenylpropanoid pathway related-genes23. Moreover, it has been well documented that

251

regulatory proteins such as TFs, seldom act alone. Numerous studies reveal that WRKY TFs physically interact

252

with a wide range of proteins including themselves, as well as VQ and MAPKs42-45. Therefore, it would be

253

meaningful to focus on determining whether CsWRKY31 and CsWRKY48 can interact with each other or other

254

proteins to co-regulate EGCG3"Me biosynthesis-related genes in the near future studies. 12

ACS Paragon Plus Environment

Page 13 of 26

Journal of Agricultural and Food Chemistry

255

In summary, two tea cultivars Jinmudan and Jinguanyin with high contents EGCG3"Me were used to explore the

256

underlying mechanism of EGCG3"Me biosynthesis. The expression levels of CsLAR, CsDFR and CCoAOMT were

257

obviously higher in Jinmudan and Jinguanyin than Fudingdabai. Moreover, we characterized two transcriptional

258

repressors CsWRKY31 and CsWRKY48, and revealed that CsWRKY31 and CsWRKY48 repressed CsLAR,

259

CsDFR and CCoAOMT transcription via directly targeting their promoters. We propose that CsWRKY31 and

260

CsWRKY48 might be acting as negatively regulators of EGCG3"Me biosynthesis in tea plants through directly

261

repressing catechins biosynthetic genes. These findings help us better understanding the transcriptional regulation

262

of EGCG3"Me biosynthesis involving WRKY TFs.

263

ABBREVIATIONS USED

264

EGCG, (-)-epigallocatechin-3-gallate; EGCG3"Me, Epigallocatechin-3-O-(3-O-methyl) gallate ; HPLC, high

265

performance liquid chromatography; DFR, dihydroflavonol; 4-reductase; LAR, leucoanthocyanidinase; reductase;

266

ANS, anthocyanidin synthase; ANR, anthocyanidin reductase; CCoAOMT, Caffeoyl-CoA 3-O-methyltransferase;

267

bHLH, basic Helix-Loop-Helix; TF, transcription factor; ORF, Open Reading Frame; NCBI, National Center for

268

Biotechnology Information; qRT-PCR, Quantitative real-time PCR; GFP,

269

RT-PCR, Reverse transcriptionpolymerase chain reaction.

270

ASSOCIATED CONTENT

271

Supporting Information

272

Primers used in this study (Table S1), and nucleotide sequences of CCoAOMT, CsLAR and CsDFR promoters

273

(Text S1)

274

Corresponding Author

275

* Telephone: +86-73184635304. Fax: +86-73184635304, email: [email protected] 13

ACS Paragon Plus Environment

Green fluorescence protein;

Journal of Agricultural and Food Chemistry

276

* Telephone: +86-73184635306. Fax: +86-73184635304, email: [email protected]

277

* Telephone: +86-73184635306. Fax: +86-73184635304, email: [email protected]

Page 14 of 26

278 279

Author Contributions

280

Yong Luo, Shuangshuang Yu, Kunbo Wang, Jianan Huang and Zhonghua Liu designed and performed the

281

experiments, analyzed the data, and co-wrote the manuscript. Juan Li and Qin Li helped with collecting materials,

282

extracting total RNA, and performing qRT-PCR.

283

Funding

284

The work was financially supported by the National Natural Science Foundation of China (31470692, 31670691

285

and 31500567).

286

Notes

287

The authors declare no competing financial interest.

288 289

REFERENCES

290

1.

291

stress-related brain aging in C57BL/6J mice. Brain Res. 2010, 1353 (2), 28-35.

292

2.

293

neurodegenerative diseases]. Journal of Hygiene Research 2010, 39 (1), 123-126.

294

3.

295

cancer and cardiovascular disease: Epidemiological evidence for multiple targeting prevention. Biofactors 2000, 13 (1-4),

296

49-54.

297

4.

298

O-methylated EGCG. J. Agric. Food Chem.2010, 58 (3), 1903-1908.

299

5.

300

O-methyltransferase from Camellia sinensis and analysis of its catalytic activity. J. Zhejiang Univ Sci. B. 2015, 16 (2),

301

103-112.

302

6.

303

containing O-methylated catechin and ginger extract enhancement. Cytotechnology 2007, 55 (2-3), 135-42.

304

7.

305

Yeast Two-Hybrid Library. Hunan Agricultural University(Master dissertation) 2017.

306

8.

Li, Q.; Zhao, H.; Ming, Z.; Zhang, Z.; Yong, L., Chronic green tea catechins administration prevents oxidative Li, Q.; Li, Y., [Review on the neuroprotective effects of green tea polyphenols for the treatment of Nakachi, K.; Matsuyama, S.; Miyake, S.; Suganuma, M.; Imai, K., Preventive effects of drinking green tea on

Kurita, I.; Maeda-Yamamoto, M.; Tachibana, H.; Kamei, M., Antihypertensive effect of Benifuuki tea containing Zhang, Y.; Lv, H. P.; Ma, C. Y.; Guo, L.; Tan, J. F.; Peng, Q. H.; Lin, Z., Cloning of a caffeoyl-coenzyme A

Maeda-Yamamoto, M.; Ema, K.; Shibuichi, I., In vitro and in vivo anti-allergic effects of 'benifuuki' green tea Yu, S., Comparative Transcriptome Analysis of Genes Involved in Tea Catechins Biosynthesis and Construction of Wang, W.; Zhou, Y.; Wu, Y.; Dai, X.; Liu, Y.; Qian, Y.; Li, M.; Jiang, X.; Wang, Y.; Gao, L.; Xia, T., Insight into 14

ACS Paragon Plus Environment

Page 15 of 26

Journal of Agricultural and Food Chemistry

307

Catechins Metabolic Pathways of Camellia sinensis Based on Genome and Transcriptome Analysis. J. Agric. Food

308

Chem. 2018, 66 (16), 4281–4293

309

9.

310

catechin synthetic genes in Camellia sinensis at different developmental stages. Bot. Stud. 2016, 57 (1), 31.

311

10. Zhao, L.; Gao, L.; Wang, H.; Chen, X.; Wang, Y.; Yang, H.; Wei, C.; Wan, X.; Xia, T., The R2R3-MYB, bHLH,

312

WD40, and related transcription factors in flavonoid biosynthesis. Funct. Integr. Genomics. 2013, 13 (1), 75-98.

313

11. Guo, F.; Guo, Y.; Wang, P.; Wang, Y.; Ni, D., Transcriptional profiling of catechins biosynthesis genes during tea

314

plant leaf development. Planta 2017, 246 (6), 1139-1152.

315

12. Agarwal, P.; Chikara, J., WRKY: its structure, evolutionary relationship, DNA-binding selectivity, role in stress

316

tolerance and development of plants. Mol. Biol. Rep. 2011, 38 (6), 3883-3896.

317

13. Tripathi, P.; Rabara, R. C.; Rushton, P. J., A systems biology perspective on the role of WRKY transcription factors

318

in drought responses in plants. Planta 2014, 239 (2), 255-266.

319

14. Eulgem, T.; Rushton, P. J.; Robatzek, S.; Somssich, I. E., The WRKY superfamily of plant transcription factors.

320

Trends Plant Sci. 2000, 5 (5), 199-206.

321

15. Luo, D. L.; Ba, L. J.; Shan, W.; Kuang, J. F.; Lu, W. J.; Chen, J. Y., Involvement of WRKY Transcription Factors

322

in Abscisic-Acid-Induced Cold Tolerance of Banana Fruit. J. Agric. Food Chem. 2017, 65 (18), 3627-3635.

323

16. Mangelsen, E.; Kilian, J.; Berendzen, K. W.; Kolukisaoglu, Ü. H.; Harter, K.; Jansson, C.; Wanke, D., Phylogenetic

324

and comparative gene expression analysis of barley ( Hordeum vulgare ) WRKY transcription factor family reveals

325

putatively retained functions between monocots and dicots. BMC Genomics 2008, 9 (1), 194.

326

17. Rushton, P. J.; Somssich, I. E.; Ringler, P.; Shen, Q. J., WRKY transcription factors. Trend Plant Sci. 2010, 15 (5),

327

247-258.

328

18. Schmutz, J.; Cannon, S. B.; Schlueter, J.; Ma, J.; Mitros, T.; Nelson, W.; Hyten, D.; Song, Q.; Thelen, J.; Cheng, J.,

329

Schmutz J, Cannon S, Schlueter J, Ma J, Mitros T, et al. (2010 ) Genome sequence of the palaeopolyploid soybean.

330

Nature 2010, 463(7278):178-183

331

19. Schluttenhofer, C.; Yuan, L., Regulation of specialized metabolism by WRKY transcription factors. Plant Physiol.

332

2015, 167 (2), 295-306.

333

20. Mishra, S.; Triptahi, V.; Singh, S.; Phukan, U. J.; Gupta, M. M.; Shanker, K.; Shukla, R. K., Wound induced

334

tanscriptional regulation of benzylisoquinoline pathway and characterization of wound inducible PsWRKY transcription

335

factor from Papaver somniferum. PloS One 2013, 8 (1), e52784.

336

21. Yokotani, N.; Sato, Y.; Tanabe, S.; Chujo, T.; Shimizu, T.; Okada, K.; Yamane, H.; Shimono, M.; Sugano, S.;

337

Takatsuji, H., WRKY76 is a rice transcriptional repressor playing opposite roles in blast disease resistance and cold

338

stress tolerance. J. Exp. Bot.2013, 64 (16), 5085-5097.

339

22. Shen, H.; Liu, C.; Zhang, Y.; Meng, X.; Zhou, X.; Chu, C.; Wang, X., OsWRKY30 is activated by MAP kinases to

340

confer drought tolerance in rice. Plant Mol. Biol. 2012, 80 (3), 241-253.

341

23. Qiu, D.; Xiao, J.; Xie, W.; Liu, H.; Li, X.; Xiong, L.; Wang, S., Rice gene network inferred from expression

342

profiling of plants overexpressing OsWRKY13, a positive regulator of disease resistance. Mol. plant 2008, 1 (3),

343

538-551.

344

24. Cheng, M. N.; Huang, Z. J.; Hua, Q. Z.; Shan, W.; Kuang, J. F.; Lu, W. J.; Qin, Y. H.; Chen, J. Y., The WRKY

345

transcription factor HpWRKY44 regulates CytP450-like1 expression in red pitaya fruit (Hylocereus polyrhizus). Hortic.

346

Res.2017, 4, 17039.

347

25. Fan, Z. Q.; Tan, X. L.; Shan, W.; Kuang, J. F.; Lu, W. J.; Chen, J. Y., Characterization of a Transcriptional

348

Regulator BrWRKY6 that Associates with Gibberellin-Suppressed Leaf Senescence of Chinese Flowering Cabbage. J.

Zhang, L. Q.; Wei, K.; Cheng, H.; Wang, L. Y.; Zhang, C. C., Accumulation of catechins and expression of

15

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 16 of 26

349

Agric. Food Chem. 2018, 66 (8), 1791-1799

350

26. Lv, H. P.; Yang, T.; Ma, C. Y.; Wang, C. P.; Shi, J.; Zhang, Y.; Peng, Q. H.; Tan, J. F.; Guo, L.; Lin, Z., Analysis of

351

naturally occurring 3″-Methyl-epigallocatechin gallate in 71 major tea cultivars grown in China and its processing

352

characteristics. J. Funct. Foods 2014, 7 (1), 727-736.

353

27. Dai, X.; Wang, Y.; Zhang, W. H., OsWRKY74, a WRKY transcription factor, modulates tolerance to phosphate

354

starvation in rice. J. Exp. Bot.2016, 67 (3), 947-960.

355

28. Wu, Z. J.; Li, X. H.; Liu, Z. W.; Li, H.; Wang, Y. X.; Zhuang, J., Transcriptome-wide identification of Camellia

356

sinensis WRKY transcription factors in response to temperature stress. Mol. Genet Genomics 2016, 291 (1), 255-269.

357

29. Phukan, U. J.; Jeena, G. S.; Shukla, R. K., WRKY Transcription Factors: Molecular Regulation and Stress

358

Responses in Plants. Front. Plant Sci. 2016, 7 (807560), 760.

359

30. W, G.; I, D. S.; DR, L.; C, L.; L, J.; G, G.; R, V. B.; M, K.; B, D. R.; B, V., Transcription factor WRKY23 assists

360

auxin distribution patterns during Arabidopsis root development through local control on flavonol biosynthesis. Proc.

361

Natl. Acad. Sci. U.S.A.2012, 109 (5), 1554-1559.

362

31. Spyropoulou, E. A.; Haring, M. A.; Schuurink, R. C., RNA sequencing on Solanum lycopersicum trichomes

363

identifies transcription factors that activate terpene synthase promoters. BMC Genomics 2014, 15 (1), 402.

364

32. Li, P.; Song, A.; Gao, C.; Wang, L.; Wang, Y.; Sun, J.; Jiang, J.; Chen, F.; Chen, S., Chrysanthemum WRKY gene

365

CmWRKY17 negatively regulates salt stress tolerance in transgenic chrysanthemum and Arabidopsis plants. Plant Cell

366

Rep. 2015, 34 (8), 1365-1378.

367

33. Tan, X.; Fan, Z.; Li, L.; Wu, Y.; Kuang, J.; Lu, W.; Chen, J., Molecular Characterization of a Leaf

368

Senescence-Related Transcription Factor BrWRKY75 of Chinese Flowering Cabbage. Horticultural Plant Journal 2016,

369

2 (5), 272-278.

370

34. Birkenbihl, R. P.; Kracher, B.; Somssich, I. E., Induced Genome-Wide Binding of Three Arabidopsis WRKY

371

Transcription Factors during Early MAMP-Triggered Immunity. Plant Cell 2017, 29 (1), tpc.00681.2016.

372

35. Liu, S.; Kracher, B.; Ziegler, J.; Birkenbihl, R. P.; Somssich, I. E., Negative regulation of ABA signaling by

373

WRKY33 is critical for Arabidopsis immunity towards Botrytis cinerea 2100. Elife 2015, 4 (e07295).

374

36. Sun, B.; Zhu, Z.; Cao, P.; Chen, H.; Chen, C.; Zhou, X.; Mao, Y.; Lei, J.; Jiang, Y.; Meng, W., Purple foliage

375

coloration in tea (Camellia sinensisL.) arises from activation of the R2R3-MYB transcription factor CsAN1. Sci Rep

376

2016, 6, 32534.

377

37. Li, M.; Li, Y.; Guo, L.; Gong, N.; Pang, Y.; Jiang, W.; Liu, Y.; Jiang, X.; Zhao, L.; Wang, Y.; Xie, D. Y.; Gao, L.;

378

Xia, T., Functional Characterization of Tea (Camellia sinensis) MYB4a Transcription Factor Using an Integrative

379

Approach. Front. Plant Sci. 2017, 8, 943.

380

38. Phukan, U. J.; Jeena, G. S.; Shukla, R. K., WRKY Transcription Factors: Molecular Regulation and Stress

381

Responses in Plants. Front. Plant Sci. 2016, 7, 760.

382

39. Wang, Y.; Shu, Z.; Wang, W.; Jiang, X.; Li, D.; Pan, J.; Li, X., CsWRKY2, a novel WRKY gene from Camellia

383

sinensis, is involved in cold and drought stress responses. Biologia Plantarum 2016, 60 (3), 443-451.

384

40. Johnson, C. S.; Kolevski, B.; Smyth, D. R., TRANSPARENT TESTA GLABRA2, a trichome and seed coat

385

development gene of Arabidopsis, encodes a WRKY transcription factor. Plant Cell 2002, 14 (6), 1359-1375.

386

41. Wang, H.; Dixon, R. A., Mutation of WRKY transcription factors initiates pith secondary wall formation and

387

increases stem biomass in dicotyledonous plants. Proc. Natl. Acad. Sci. U.S.A.2010, 107 (51), 22338-22343.

388

42. Akagi, A.; Fukushima, S.; Okada, K.; Jiang, C. J.; Yoshida, R.; Nakayama, A.; Shimono, M.; Sugano, S.; Yamane,

389

H.; Takatsuji, H., WRKY45-dependent priming of diterpenoid phytoalexin biosynthesis in rice and the role of cytokinin

390

in triggering the reaction. Plant Mol. Biol. 2014, 86 (1-2), 171-183. 16

ACS Paragon Plus Environment

Page 17 of 26

Journal of Agricultural and Food Chemistry

391

43. D, R.; Y, L.; KY, Y.; L, H.; G, M.; J, G.; S, Z., A fungal-responsive MAPK cascade regulates phytoalexin

392

biosynthesis in Arabidopsis. Proc. Natl. Acad. Sci. U.S.A. 2008, 105 (14), 5638-5643.

393

44. Ye, Y. J.; Xiao, Y. Y.; Han, Y. C.; Shan, W.; Fan, Z. Q.; Xu, Q. G.; Kuang, J. F.; Lu, W. J.; Prakash, L.; Chen, J.

394

Y., Banana fruit VQ motif-containing protein5 represses cold-responsive transcription factor MaWRKY26 involved in

395

the regulation of JA biosynthetic genes. Sci. Rep.2013, 6, 23632.

396

45. Sébastien, B.; Li, J.; Tapio, P. E., WRKY54 and WRKY70 co-operate as negative regulators of leaf senescence

397

inArabidopsis thaliana. J. Exp. Bot. 2012, 63 (7), 2667-2679.

398 399 400 401 402 403 404

Figure legends

405 406 407

Figure 1. The biosynthetic pathways of EGCG3"Me, and the plant materials of tea leaves.(a)Possible biosynthetic pathways of EGCG3"Me. (b) Plant materials of tea leaves (a terminal bud and two young leaves).

408 409 410 411 412 413 414 415 416 417 418

Figure 2. The contents of EGCG3"Me and EGCG, and the expression pattern of EGCG3"Me biosynthesis related-genes. (a)The content of EGCG3"Me in different tea cultivars by HPLC in dry weight of leaf sample, “n.d.”indicates not detected. (b)The content of EGCG in different tea cultivars. (c)Relative expression analyzed by real time quantitative reverse transcription-PCR (qRT-PCR) of EGCG3"Me biosynthesis related-genes, β-actin was used as an internal control. Each value represents the mean ± SE of three replicates. *P ≤ 0.05, **P ≤ 0.01, compared with the Fudingdabai, respectively. “NS” indicate no differences by Student’s t test. (Note: PAL, phenylalanine ammonia-lyase; C4H, cinnamate 4-hydroxylase;4CL, 4-coumarate:CoA ligase; CHS, chalcone synthase; CHI, chalcone isomerase; F3H, flavanone 3b-hydroxylase; F3’H, flavonoid 30-hydroxylase; F3’5H, flavonoid 3’5’-hy-droxylase; DFR, dihydroflavonol reductase; LAR, leuacoanthocyanidin reductase; ANS, anthocyanin synthase; ANR, anthocyanidin reductase; CCoAOMT, caffeoyl-CoA 3-O-methyltransferase.)

419 420 421 422 423 424 425 426 427

Figure 3. Identification and amino acid analysis of CsWRKY31 and CsWRKY48 involved in EGCG3"Me biosynthesis in tea tree (Camellia sinensis). (a)Comparative sequence analysis and dendrogram of CsWRKY31 and CsWRKY48 protein with the other WRKY proteins. The conserved WRKY domain and the C2H2 zinc-finger motifs are indicated by red letters. (b) Phylogenetic tree analysis of CsWRKY31 and CsWRKY48 protein comparison with amino acid sequences of the WRKY family isolated from Arabidopsis. The full-length amino acid sequences were downloaded from the institute for Genomic Research (http://www.tigr.org) and the National Center for Biotechnology information (http://www.ncbi.nlm.nih.gov). The amino acid sequences were aligned with ClustalW and the phylogenetic tree was constructed by MEGA 5.0 using the neighbor-joining method .

428 429 430 431

Figure 4. Expression pattern of CsWRKY31 and CsWRKY48 in different tea cultivars and their subcellular localizations. (a) Relative expressions analyzed by real time quantitative reverse transcription-PCR (qRT-PCR). **Significant differences in values (P < 0.01) by Student's t-test. (b)Subcellular localizations of CsWRKY31 and CsWRKY48 in 17

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

432 433 434

Page 18 of 26

tobacco leaves. CsWRKY31 and CsWRKY48 fused with the GFP or GFP control were infiltrated into tobacco leaves via Agrobacterium tumefaciens strain GV3101. After 48 h of the infiltration, GFP fluorescence was visualized using a fluorescence microscope, Scale bars=50µm.

435 436 437 438 439 440 441 442 443 444 445 446 447 448

Figure 5. Molecular Characterization of CsWRKY31 and CsWRKY48. (a)The transcriptional activity analysis of CsWRKY31 and CsWRKY48 in yeast cell assays. Both the full length of the two WRKYs are fused with pGBKT7, and transformed yeasts were selected on SD/-Trp or SD/-Trp/-His/-Ade/ X-α-gal media for 3-5 d at 30°C. Transcription activation was monitored by detecting yeast growth and the expression of a-galactosidase (a-Gal) activity. (b & c)Transcriptional activation ability of CsWRKY31 and CsWRKY48 in tobacco leaves. The double-reporter plasmids contain 5×GAL4 and mini CaMV35S fused to firefly luciferase (LUC) and renilla luciferase (REN) driven by CaMV35S. The effector plasmids contain the CsWRKY31 and CsWRKY48 genes fused to GAL4BD driven by the CaMV35S. The dual REN/LUC reporter and effectors were co-transformed into tobacco leaves by Agrobacterium tumefaciensstrain GV3101. After 48-72h of the infiltration, LUC and REN luciferase activities were assayed, and the transcription activation ability of CsWRKY31 and CsWRKY48 is indicated by the ratio of LUC to REN. Each value represents the means of six biological replicates, and vertical bars represent the S.E. **Significant differences in values (P < 0.01) by Student's t-test, compared with pBD-empty.

449 450 451 452 453 454 455 456 457

Figure 6. Electrophoretic mobility shiſt assay (EMSA) showing both CsWRKY31 and CsWRKY48 can specifically bind to the promoter region of EGCG3"Me biosynthetic genes including CCoAOMT, CsLAR, and CsDFR. (a&c) SDS–PAGE gel stained with Coomassie blue demonstrating affinity purification of the recombinant GST-CsWRKY31-N and GST-CsWRKY48 protein used for the EMSA. (b&d) Recombinant GST-CsWRKY31-N and GST-CsWRKY48 protein bind directly to the promoters of CCoAOMT, CsLAR, and CsDFR containing W-box(T/CTGACT/C) element, respectively. Biotin-labeled DNA probe from the promoters or mutant probe was incubated with the two recombinant protein, and the DNA-protein complexes were separated on 6% native polyacrylamide gels. +and++indicate increasing amounts unlabeled probes for competition.

458 459 460 461 462 463 464

Figure 7. CsWRKY31 and CsWRKY48 repressed the promoter activity of genes involved in the EGCG3"Me biosynthesis pathways.(a)Constructs used in the transient transactivation assays. (b)The LUC/REN ratio was calculated as the final transcriptional activity. Empty vector was used as the effector in the control assay. The promoters of CCoAOMT, CsLAR, and CsDFR genes were used in dual luciferase assays, and each value represents the means of six biological replicates. Vertical bars represent the S.E.. Significant differences in values (*P < 0.05, **P < 0.01) by Student's t-test.

465 466 467 468 469 18

ACS Paragon Plus Environment

Page 19 of 26

Journal of Agricultural and Food Chemistry

Figure 1. The biosynthetic pathways of EGCG3"Me, and the plant materials of tea leaves.(a)Possible biosynthetic pathways of EGCG3"Me. (b) Plant materials of tea leaves (a terminal bud and two young leaves). 419x217mm (300 x 300 DPI)

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Figure 2. The contents of EGCG3"Me and EGCG, and the expression pattern of EGCG3"Me biosynthesis related-genes. (a)The content of EGCG3"Me in different tea cultivars by HPLC in dry weight of leaf sample, “n.d.”indicates not detected. (b)The content of EGCG in different tea cultivars. (c)Relative expression analyzed by real time quantitative reverse transcription-PCR (qRT-PCR) of EGCG3"Me biosynthesis relatedgenes, β-actin was used as an internal control. Each value represents the mean ± SE of three replicates. *P ≤ 0.05, **P ≤ 0.01, compared with the Fudingdabai, respectively. “NS” indicate no differences by Student’s t test. (Note: PAL, phenylalanine ammonia-lyase; C4H, cinnamate 4-hydroxylase;4CL, 4-coumarate:CoA ligase; CHS, chalcone synthase; CHI, chalcone isomerase; F3H, flavanone 3b-hydroxylase; F3’H, flavonoid 30-hydroxylase; F3’5H, flavonoid 3’5’-hy-droxylase; DFR, dihydroflavonol reductase; LAR, leuacoanthocyanidin reductase; ANS, anthocyanin synthase; ANR, anthocyanidin reductase; CCoAOMT, caffeoyl-CoA 3-O-methyltransferase.) 422x130mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 20 of 26

Page 21 of 26

Journal of Agricultural and Food Chemistry

Figure 3. Identification and amino acid analysis of CsWRKY31 and CsWRKY48 involved in EGCG3"Me biosynthesis in tea tree (Camellia sinensis). (a)Comparative sequence analysis and dendrogram of CsWRKY31 and CsWRKY48 protein with the other WRKY proteins. The conserved WRKY domain and the C2H2 zinc-finger motifs are indicated by red letters. (b) Phylogenetic tree analysis of CsWRKY31 and CsWRKY48 protein comparison with amino acid sequences of the WRKY family isolated from Arabidopsis. The full-length amino acid sequences were downloaded from the institute for Genomic Research (http://www.tigr.org) and the National Center for Biotechnology information (http://www.ncbi.nlm.nih.gov). The amino acid sequences were aligned with ClustalW and the phylogenetic tree was constructed by MEGA 5.0 using the neighbor-joining method . 170x203mm (300 x 300 DPI)

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Figure 4. Expression pattern of CsWRKY31 and CsWRKY48 in different tea cultivars and their subcellular localizations. (a) Relative expressions analyzed by real time quantitative reverse transcription-PCR (qRTPCR). **Significant differences in values (P < 0.01) by Student's t-test. (b)Subcellular localizations of CsWRKY31 and CsWRKY48 in tobacco leaves. CsWRKY31 and CsWRKY48 fused with the GFP or GFP control were infiltrated into tobacco leaves via Agrobacterium tumefaciens strain GV3101. After 48 h of the infiltration, GFP fluorescence was visualized using a fluorescence microscope, Scale bars=50µm. 170x59mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 22 of 26

Page 23 of 26

Journal of Agricultural and Food Chemistry

Figure 5. Molecular Characterization of CsWRKY31 and CsWRKY48. (a)The transcriptional activity analysis of CsWRKY31 and CsWRKY48 in yeast cell assays. Both the full length of the two WRKYs are fused with pGBKT7, and transformed yeasts were selected on SD/-Trp or SD/-Trp/-His/-Ade/ X-α-gal media for 3-5 d at 30°C. Transcription activation was monitored by detecting yeast growth and the expression of agalactosidase (a-Gal) activity. (b & c)Transcriptional activation ability of CsWRKY31 and CsWRKY48 in tobacco leaves. The double-reporter plasmids contain 5×GAL4 and mini CaMV35S fused to firefly luciferase (LUC) and renilla luciferase (REN) driven by CaMV35S. The effector plasmids contain the CsWRKY31 and CsWRKY48 genes fused to GAL4BD driven by the CaMV35S. The dual REN/LUC reporter and effectors were co-transformed into tobacco leaves by Agrobacterium tumefaciensstrain GV3101. After 48-72h of the infiltration, LUC and REN luciferase activities were assayed, and the transcription activation ability of CsWRKY31 and CsWRKY48 is indicated by the ratio of LUC to REN. Each value represents the means of six biological replicates, and vertical bars represent the S.E. **Significant differences in values (P < 0.01) by Student's t-test, compared with pBD-empty. 170x109mm (300 x 300 DPI)

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Figure 6. Electrophoretic mobility shiſt assay (EMSA) showing both CsWRKY31 and CsWRKY48 can specifically bind to the promoter region of EGCG3"Me biosynthetic genes including CCoAOMT, CsLAR, and CsDFR. (a&c) SDS–PAGE gel stained with Coomassie blue demonstrating affinity purification of the recombinant GST-CsWRKY31-N and GST-CsWRKY48 protein used for the EMSA. (b&d) Recombinant GSTCsWRKY31-N and GST-CsWRKY48 protein bind directly to the promoters of CCoAOMT, CsLAR, and CsDFR containing W-box(T/CTGACT/C) element, respectively. Biotin-labeled DNA probe from the promoters or mutant probe was incubated with the two recombinant protein, and the DNA-protein complexes were separated on 6% native polyacrylamide gels. +and++indicate increasing amounts unlabeled probes for competition. 135x181mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 24 of 26

Page 25 of 26

Journal of Agricultural and Food Chemistry

Figure 7. CsWRKY31 and CsWRKY48 repressed the promoter activity of genes involved in the EGCG3"Me biosynthesis pathways.(a)Constructs used in the transient transactivation assays. (b)The LUC/REN ratio was calculated as the final transcriptional activity. Empty vector was used as the effector in the control assay. The promoters of CCoAOMT, CsLAR, and CsDFR genes were used in dual luciferase assays, and each value represents the means of six biological replicates. Vertical bars represent the S.E.. Significant differences in values (*P < 0.05, **P < 0.01) by Student's t-test. 256x323mm (300 x 300 DPI)

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

TOC 84x47mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 26 of 26