Caffeine Content and Related Gene Expression: Novel Insight into

Mar 4, 2019 - Caffeine Content and Related Gene Expression: Novel Insight into Caffeine Metabolism in Camellia Plants Containing Low, Normal, and High...
0 downloads 0 Views 2MB Size
Subscriber access provided by WEBSTER UNIV

Chemistry and Biology of Aroma and Taste

Caffeine content and related gene expression: novel insight into caffeine metabolism in Camellia plants containing low, normal and high caffeine concentrations Biying Zhu, Lin-Bo Chen, Mengqian Lu, Jing Zhang, Jieyun Han, Wei-Wei Deng, and Zheng-Zhu Zhang J. Agric. Food Chem., Just Accepted Manuscript • DOI: 10.1021/acs.jafc.9b00240 • Publication Date (Web): 04 Mar 2019 Downloaded from http://pubs.acs.org on March 5, 2019

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 42

Journal of Agricultural and Food Chemistry

1

Caffeine content and related gene expression: novel insight into caffeine

2

metabolism in Camellia plants containing low, normal and high caffeine

3

concentrations

4 5

Biying Zhu#1, Lin-Bo Chen#2, Mengqian Lu1, Jing Zhang1, Jieyun Han1, Wei-Wei

6

Deng1*, Zheng-Zhu Zhang1*

7 8

1State

9

Science & Technology, Anhui Agricultural University, Hefei, Anhui 230036, China

Key Laboratory of Tea Plant Biology and Utilization, School of Tea and Food

10

2Tea

11

Yunnan 666201, China

Research Institute, Yunnan Academy of Agricultural Sciences, Menghai,

12 13

# These

14

*Corresponding

15

Zhang ([email protected]) ; Tel/fax: +86 551 65785471

authors (Biying Zhu and Lin-Bo Chen) contribute equally. author: Wei-Wei Deng ([email protected]), Zheng-Zhu

16 17 18 19 20 21 22 23

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

24

Abstract

25

Caffeine is a crucial secondary metabolic product in tea plants. Although the

26

presence of caffeine in tea plants has been identified, the molecular mechanisms

27

regulating relevant caffeine metabolism remain unclear. To elucidate the caffeine

28

biosynthesis and catabolism in Camellia plants, fresh, germinated leaves from four

29

Camellia plants with low (2), normal (1) and high (1) caffeine concentrations, namely

30

low-caffeine tea 1 (LCT1, Camellia crassicolumna), low-caffeine tea 2 (LCT2, C.

31

crassicolumna), Shuchazao (SCZ, C. sinensis) and Yunkang 43 (YK43, C. sinensis)

32

were used in this research. Transcriptome and purine alkaloids analyses of these

33

Camellia leaves were performed using RNA-Seq and liquid chromatography–mass

34

spectrometry (LC-MS). Moreover,

35

the metabolic fate of caffeine in leaves of these plants. Caffeine content was

36

correlated with related genes expression levels, and a quantitative real-time (qRT)

37

PCR analysis of specific genes showed a consistent tendency with the obtained

38

transcriptomic analysis. Based on the results of stable isotope-labelled tracer

39

experiments, we discovered a degradation pathway of caffeine to theobromine. These

40

findings could assist researchers in understanding the caffeine-related mechanisms in

41

Camellia plants containing low, normal, and high caffeine content and be applied to

42

caffeine regulation and breeding improvement in future research.

43

Keywords: Camellia plants, caffeine metabolism, catechins, differentially expressed

44

genes, RNA-Seq, theanine

15N-caffeine

tracing was performed to determine

45

ACS Paragon Plus Environment

Page 2 of 42

Page 3 of 42

Journal of Agricultural and Food Chemistry

46

Introduction

47

The tea plant (Camellia sinensis) is one of evergreen perennial plants with a

48

lifespan of more than 100 years. Tea, as the second most popular (after water) natural

49

nonalcoholic beverage, is processed from the leaves of tea plants.1 Numerous studies

50

have suggested that tea can prevent cancer and other neurodegenerative or

51

cardiovascular diseases.2-4 However, due to its high caffeine content, high tea intake

52

can affect those sensitive to caffeine by increasing their anxiety, blood pressure,

53

insomnia, etc.5-6

54

Caffeine, 1,3,7-trimethylxanthine, is found in tea, mate, cocoa, coffee and some

55

other plant species.7 These purine alkaloids are accumulated in young leaves,

56

cotyledons, seeds and fruits. In tea leaves, these alkaloids typically contain 2%–5%

57

(w/w) caffeine. The typical biosynthetic pathways of caffeine are: xanthosine (XR) →

58

7-methyxanthosine (7mXR) → 7-methylxanthine (7mX) → theobromine (Tb) →

59

caffeine (Cf) as the major pathway; 7-methylxanthine (7mX) → paraxanthine (Px)

60

→ caffeine (Cf) as a minor route in leaves of tea plants.8 N-methyltransferases (NMTs)

61

were reported to catalyze the methylation steps (the 1st, 3rd and 4th steps) in the

62

major

63

S-adenosyl-L-methionine (SAM).9 Among them, caffeine synthase (CS) has ability to

64

catalyze the final two steps (7-methylxanthine → theobromine → caffeine) and shows

65

a bifunction of two NMTs.10 The paraxanthine NMT (catalyzing the steps of

66

paraxanthine → caffeine) was reported to be existed in tea chloroplasts.

67

activity of caffeine biosynthesis has been investigated in the petals and stamens

caffeine

biosynthetic

pathway

with

ACS Paragon Plus Environment

the

methyl

donor

11

of

A high

Journal of Agricultural and Food Chemistry

68

(before flowering) of tea plants.12 A tea CS gene that encodes CS (TCS1; GenBank

69

accession no. AB031280), cloned from young leaves of tea plant (Camellia sinensis),

70

was reported by Kato et al.13 An analysis of the expression patterns of TCS1 (encoded

71

TCS) in the different organs of tea seedlings revealed that the highest expression was

72

found in young leaves. The biosynthesis of caffeine might be closely related with the

73

expression levels of TCS1.14 The association analysis on the concentrations of

74

caffeine and the expressions of TCS1 was elucidated in tea and other related species.15

75

It was reported that there were six types in the alleles of TCS1, TCS1a~f. Among

76

them, TCS1a showed as the predominant position of allele; TCS1b-f appeared as the

77

rare alleles which only existed in some wild species.16 Hongyacha, as a kind of wild

78

tea species, was reported to show a big difference in the characteristics of the

79

morphology with a previous reported caffeine-free Cocoa tea (Camellia ptilophylla,

80

CCT) in China.17

81

In contrast to caffeine biosynthesis, purine alkaloid catabolism and caffeine

82

degradation in Camellia plants is relatively unknown. A degradation pathway of

83

caffeine has been previously reported in coffee. Caffeine can be subsequently

84

catabolized to theophylline, 3-methylxanthine and xanthine. And xanthine enters the

85

purine catabolism, and is finally degraded to NH3 and CO2 [xanthine (X) → uric acid

86



87

pathway has been noted in some microorganisms.19-21 And the activity of demethylase

88

has been reported in some microorganisms, but not been investigated in plants,

89

especially in coffee or tea thus far.

allantoin → allantoic acid → urea → NH3 and CO2].18 The caffeine degradation

ACS Paragon Plus Environment

Page 4 of 42

Page 5 of 42

Journal of Agricultural and Food Chemistry

90

To determine the caffeine metabolism in Camellia plants, four plants LCT1 (C.

91

crassicolumna), LCT2 (C. crassicolumna), SCZ (C. sinensis) and YK43 (C. sinensis),

92

with low, normal, and high caffeine content in leaves were performed in this research

93

(Figure 1A). RNA-Seq was used to analyze the transcriptomes and identify the related

94

genes in the pathway of caffeine metabolism. To elucidate the caffeine metabolic fate

95

in these Camellia plants, a stable isotope tracer experiment involving [15N2]-caffeine

96

was conducted. The obtained results can assist researchers in understanding the

97

caffeine mechanisms in Camellia plants and be applied to caffeine regulation and

98

breeding improvement.

99 100

Materials and methods

101

Plant materials

102

Fresh young leaves (one apical bud with two terminal leaves) from these Camellia

103

plants (LCT1, LCT2, SCZ, and YK43) were picked at the germplasm resource garden

104

in the Tea Research Institute of Yunnan Academy of Agricultural Sciences, Menghai,

105

Yunnan, China, in June 2017. The freshly plucked leaf samples were stored at −80 °C

106

(immediately frozen in liquid nitrogen first) until they were used for further analyses.

107 108

The extraction of RNA, construction of library, and sequencing

109

Total RNAs were extracted from the leaf samples by using an RNA-prep Pure Plant

110

Kit (TianGen, Beijing, China). Total RNA quantity was evaluated by a bioanalyzer of

111

Agilent 2100 (Agilent, Santa Clara, CA, USA). A NanoDropTM ultraviolet

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

112

spectrophotometer (Thermo, Waltham, MA, USA) was used to estimate the RNA

113

integrity. cDNA libraries were established by a kit of NEBNext Ultra RNA Library

114

Perp (Gene, Beijing, China) and sequenced by the Illumina HiSeq 4000 platform

115

(Illumina, San Diego, CA) in BGI, Shenzhen, China.

116 117

Data filtering and read mapping

118

High-quality clean reads were acquired by removing adaptor sequences, duplicated

119

sequences, ambiguous (with the ratio of “N”>5%) and low quality reads (the reads in

120

which quality values of ≤15 comprised more than 20%). All the analyses were

121

performed with the clean high quality data.22 The clean reads with high quality from

122

samples were mapped to a tea plant genome

123

Spliced Alignment of Transcripts).24

23

by HISAT (Hierarchical Indexing for

124 125

Functional annotation and analysis of differentially expressed genes

126

In order to get functional annotations of the proteins, the unigene sequences of

127

samples were searched by using BLASTX against nonredundant protein (NR)

128

database, nonredundant nucleotide (NT) database, clusters of orthologous groups of

129

proteins (COG) database, Kyoto Encyclopedia of Genes and Genomes (KEGG), and

130

Swiss-Prot annotated protein sequence database.25 The program of Blast2GO (version

131

2.5.0) was employed to receive the gene ontology (GO) annotations based on NR.26

132

We used InterProScan (version 5.11-51.0) to annotate the unigenes.27

133

After annotations were obtained, expression levels of the unigenes were calculated

ACS Paragon Plus Environment

Page 6 of 42

Page 7 of 42

Journal of Agricultural and Food Chemistry

134

with RSEM software28 and the values of FPKM (fragments per kilobase of transcript

135

sequence per millions of base pairs sequenced).29 Transcriptome data were also

136

analyzed to identify the DEGs (differentially expressed genes) using DEseq2.30 The

137

FDR (false discovery rate), the correction parameter, was used to determine DEGs.31

138

Genes with a |log2ratio| of ≥1 and FDR of 0.90 were considered to indicate

184

good correlation and used to construct TF-gene-metabolite network. Moreover, the

185

network was displayed using Cytoscape (version 3.6.0).37 The expression of the target

186

genes in caffeine metabolic pathway was included for further analysis.

187 188 189

Stable isotope-labeling tracer experiment for [15N2]-caffeine [15N2]-Caffeine (1, 3-15N2, 99%) was obtained from Cambridge Isotope 15N-tracer

190

Laboratories, MA, USA. The

was conducted according to the publication

191

of Ito and Ashihara,38 with a slight modification. The segments of fresh leaves

192

(approximately 100 mg) were placed in a flasks with 2 mL of K-Pi buffer (30 mM,

193

pH 5.6, comprising 20 mM sucrose, 1% sodium ascorbate), and 15N2-labelled caffeine

194

(10 mM). For comparison, caffeine samples (10 mM [without 15N labeling] and 0 mM)

195

were also used as controls. The flasks were treated at 26 °C for 12, 24, and 48 h. After

196

incubation, samples were collected and frozen in liquid N2 for further analysis. For

197

the analysis of 15N-metabolites, samples were homogenized with 80% aqueous MeOH.

198

The fraction of MeOH soluble was concentrated and examined using an LC–MS

199

system (Palo Alto, CA, USA). The detection method was the same as the method that

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

200

previously described.

201 202

Determination of other main metabolites

203

Catechins were analyzed with the reported method by Jiang39 with minor

204

modifications. A 0.06 g freeze-dried sample was ground with 1 mL of 80% methanol

205

containing 1% acetic acid, and then added to 2 mL of the extraction solution. After

206

the extraction was centrifuged at 16060 g for 10 min at 4 °C, the residue was

207

re-extracted as aforementioned. Subsequently, supernatant was set to 4 mL and

208

filtered by a membrane of 0.22 μm. Standard compounds and extracted samples were

209

determined by a Waters e2695 HPLC system (Milford, MA, USA). A Phenomenex

210

C18 column (5 μm, 25 cm × 4.6 mm, Torrance, California, USA), was performed at a

211

flow rate of 1.0 mL/min, and at 35 °C. The injection volume of sample was set to 10

212

μL. And the mobile phase (100% methanol and 0.2% acetic acid in water) was

213

employed with a gradient program of the acetic acid in water, which was increased

214

linearly from 95% to 80% in 2 min, to 75% in 12 min, to 58% in 6 min, to 58% in 2

215

min, to 0% in 6 min, held at 0% until 31 min, to 95% in 4 min, and then, held at 95%

216

until 38 min. Eluate was analyzed by absorbance at 278 nm.

217

Free amino acids were extracted from samples according to the published method

218

from Wang et al.40 Freeze-dried sample (0.20 g) was powdered and extracted with 2

219

mL of 4% sulfonyl salicylic acid by ultrasonic extraction for 30 min. The extraction

220

was centrifuged at 13 680 g for 30 min, and the supernatant was metered and filtered

221

by a membrane with 0.45 μm. Concentrations of amino acids were determined

ACS Paragon Plus Environment

Page 10 of 42

Page 11 of 42

Journal of Agricultural and Food Chemistry

222

according to the method of Liao 41 by a Hitachi amino acid analyzer (L-8900, Tokyo,

223

Japan) with a 4.60 mm × 60 mm column (Hitachi ion exchange resin).

224 225

Statistical analysis

226

Three independent biological replicates and technical replicates were used for

227

calculating the mean value and the standard deviation (SD) of the metabolite

228

concentrations. SPSS 17.0 software was used to determine the significant differences

229

by Duncan's multiple-range test at the 5% level. Double coordinate figure was

230

produced with Prism 7.0 software. Heat maps were performed with R package. The

231

correlations were analyzed via Pearson correlation and the network was displayed

232

using Cytoscape (version 3.6.0).

233 234

Results

235

RNA sequencing and reference genome alignment

236

For RNA-Seq, high-quality total RNAs of samples were reverse transcribed into

237

cDNAs, and amplified to construct 12 cDNA libraries. Totally, 80.36 Gb clean reads

238

were constructed with an average value of 6.72 Gb each sample. Q20 (base-calling

239

error probability of 99%) values were more than 96.48% and the Q30 values were

240

more than 89.61% (base-calling error probability of 99.90%; Table 1). The average

241

ratio of the sample map to the genome was 76.82%. These results showed the

242

obtained high-quality transcriptomic data could be used for the further analysis.

243

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

244

DEGs in these Camellia plants

245

The DEGs between the control and comparable sample were defined by the

246

fold-change value of the FPKM. As SCZ was the plant material used for genome

247

sequencing23 and SCZ leaves contain normal caffeine concentrations, we chose SCZ

248

as the control for this experiment. In our transcriptome data, 14 117 DEGs were

249

detected among these four Camellia plants. In SCZ-VS-LCT1, 5701 DEGs were

250

detected, of which 2725 were upregulated and 2976 were downregulated. In

251

SCZ-VS-LCT2, 8240 DEGs were identified, of which 3580 were upregulated and

252

4660 were downregulated. In SCZ-VS-YK43, 9257 DEGs were noted, of which 3923

253

were upregulated and 5334 were downregulated (Figures 1B and S1).

254 255

Functional enrichment analysis of DEGs using GO and KEGG

256

All unigenes of these Camellia plants were annotated using BLASTX. And we

257

conducted KEGG and GO enrichment analyses by using the reference genes as the

258

background to know the function of DEGs. GO is an international standard system for

259

gene functional classification26 that can fully describe the biological category of genes

260

and reflect information about genes involved in an organism’s metabolism. GO

261

categories include molecular function (MF), cellular component (CC) and biological

262

process (BP).42 In the BP category, the metabolic process, the cellular process, and the

263

single-organism process were the top three enriched GO terms. In the CC category,

264

cell, cell part, and membrane were the top terms. And in the category of MF, catalytic

265

activity and binding were the top terms.

ACS Paragon Plus Environment

Page 12 of 42

Page 13 of 42

Journal of Agricultural and Food Chemistry

266

In SCZ-VS-LCT1, 4937, 3153 and 5122 DEGs were enriched to categories of CC,

267

MF and BP, respectively. In SCZ-VS-LCT2, 7034, 4626 and 7740 DEGs were

268

enriched to categories of CC, MF and BP, respectively. In SCZ-VS-YK43, 7682,

269

5276 and 8681 DEGs were enriched to categories of CC, MF and BP, respectively.

270

GO classification and the upregulated and downregulated DEGs in the three

271

comparisons are shown in Figure S2.

272

We also constructed the KEGG analysis with the obtained unigenes against the

273

reference canonical pathways. The results revealed 131, 132, and 132 metabolic

274

pathways in the SCZ-VS-LCT1, SCZ-VS-LCT2, and SCZ-VS-YK43 samples,

275

respectively. These pathways included carbon metabolism, flavone and flavonol

276

biosynthesis, purine metabolic pathway, and caffeine metabolic pathway. Among

277

them, the target pathways related to the concerned ones were selected (Figure 2). A

278

substantial number of DEGs were noted in the metabolic pathways and secondary

279

metabolite biosynthesis. Alanine, aspartate, and glutamate metabolism; anthocyanin

280

biosynthesis; secondary metabolite biosynthesis; caffeine metabolism; flavonoid

281

biosynthesis; flavone and flavonol biosynthesis; phenylalanine metabolism; and

282

purine metabolism resulted in different degrees of enrichment in SCZ-VS-LCT1,

283

SCZ-VS-LCT2, and SCZ-VS-YK43. Of them, 22 genes were involved in caffeine

284

metabolism, of which 4, 8, and 4 DEGs were identified from SCZ-VS-LCT1,

285

SCZ-VS-LCT2, and SCZ-VS-YK43, respectively.

286 287

Purine alkaloid content in these Camellia plants

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 14 of 42

288

LC–MS was used to determine purine alkaloids with a high sensitivity and

289

accuracy. The contents of purine alkaloids in these Camellia leaves were shown in

290

Table 2. Among these four Camellia plants, YK43 showed the highest total purine

291

alkaloids content (46.68 mg/g), followed by SCZ (35.49 mg/g), LCT1 (15.02 mg/g),

292

and LCT2 (8.57 mg/g). Among these purine alkaloids, the caffeine concentrations in

293

YK43, SCZ, LCT1, and LCT2 were 42.82, 34.17, 8.33 and 0.59 mg/g DW,

294

respectively. Moreover, the caffeine concentrations in LCT1 and LCT2 were

295

considerably lower than that in YK43 and SCZ. However, the concentration of

296

another purine alkaloid, theobromine, in LCT1 (5.43 mg/g) and LCT2 (7.70 mg/g)

297

was considerably higher than that in YK43 (3.50 mg/g) and SCZ (1.27 mg/g). In

298

addition, the concentration of theophylline was substantially higher in the leaf

299

samples of SCZ (4.77 μg/g) and LCT2 (2.46 μg/g) than in those of LCT1 (0.29 μg/g)

300

and YK43 (0.53 μg/g). In the leaves of LCT2, 13.02 μg/g of xanthine was detected,

301

which was considerably higher than the xanthine concentration in the other three

302

Camellia plants.

303 304

Expressions of the genes in the caffeine metabolism

305

In addition to the pathway of “caffeine (Cf) → theophylline (Tp) →

306

3-methylxanthine (3mX) → xanthine (X) → uric acid → allantoin → allantoic acid

307

→ urea → CO2 + NH3”, 43 we obtained the other five degradation pathways annotated

308

in

309

(https://www.kegg.jp/kegg-bin/highlight_pathway?scale=1.0&map=map00232&keyw

KEGG

ACS Paragon Plus Environment

Page 15 of 42

Journal of Agricultural and Food Chemistry

310

ord=caffeine): “caffeine (Cf) → paraxanthine (Px) → 1,7-dimethyluric acid, ”

311

“caffeine (Cf) → paraxanthine (Px) → 7-methylxanthine (7-mX) → 7-methyluric

312

acid, ”

313

1-methyluric acid,” “caffeine (Cf) → theobromine (Tb) → 3,7-dimethyuric acid,”

314

and “caffeine (Cf) → 1,3,7-trimethyluric acid → 3,6,8-trimethylallantoin” (Figure 3).

315

The identification of these pathways is of considerable importance in research related

316

to the metabolism of caffeine in Camellia plants.

“caffeine (Cf) → paraxanthine (Px) → 1-methylxanthine (1-mX) →

317

The expression levels of genes in caffeine’s metabolic pathway were determined to

318

ascertain the association between the related gene expression in caffeine metabolism

319

and purine alkaloid content in these Camellia plants. Genes involved in caffeine

320

metabolism could be grouped into five categories, namely de novo, purine salvage,

321

caffeine degradation, caffeine synthesis and methyl donor synthesis44-45. And the

322

expression patterns of these unigenes in caffeine metabolism were shown in Figure 4.

323

TEA022617.1 (encoded adenylosuccinate synthase, ASS) in the de novo pathway,

324

TEA006735.1 (encoded S-adenosylmethionine synthetase, SAMS) in the methyl donor

325

synthesis, and TEA019288.1 (encoded adenosine kinase, ADK) were highly expressed

326

in these four Camellia plants. In the case of purine alkaloids synthesis pathway, the

327

expression

328

TEA012581.1 (encoded tea caffeine synthase, TCS) were highly expressed in

329

Camellia plants with high caffeine concentration. The general expression trends of

330

unigenes involved in caffeine degradation pathways were higher in the Camellia

331

plants with low caffeine content (Figure 4). From the KEGG enrichment analysis, 4, 8,

patterns

of

TEA015791.1,

TEA028050.1,

ACS Paragon Plus Environment

TEA028051.1,

and

Journal of Agricultural and Food Chemistry

332

and 4 DEGs were obtained in SCZ-VS-LCT1, SCZ-VS-LCT2, and SCZ-VS-YK43,

333

respectively. These DEGs could be annotated as cytochrome P450 family 1 subfamily

334

A polypeptide 2 (CYP1A2), xanthine oxidase (XO), and urate oxidase (UOX), all of

335

which are in the degradation pathway of caffeine (Table S2).

336

To further validate the gene expression levels, 20 genes related to purine salvage,

337

caffeine degradation, caffeine synthesis, and the methyl donor synthesis pathways in

338

the caffeine metabolism were selected for verification by qRT-PCR. These obtained

339

results indicated that the gene expressions showed a consistent tendency with the

340

transcriptome data (Figure 5).

341 342

Metabolic fate of [15N2]-caffeine in the leaves of these Camellia plants

343

Stable isotope-labeling tracer experiment was an effective approach for

344

investigating the metabolic pathways that has been comprehensive reported.46-48 To

345

investigate the mechanism of caffeine catabolism in these Camellia leaves, a stable

346

isotope-labeling tracer was conducted to identify the metabolites involved in caffeine

347

degradation. 15N-labelled metabolites were analyzed by LC–MS, and identified by the

348

standards or the reported tR, ([M+H]+/[M-H]-), major fragment ions, and λmax.49

349

[M+H]+ at m/z 195, and fragment ions at m/z 110 and 138 were determined to be

350

caffeine. According to the second-order MS results, the [M+H]+ of [15N2]-caffeine

351

was at m/z 197, and fragment ions were at m/z 139 and 111. Furthermore, [M+H]+ at

352

m/z 181, and fragment ions at m/z 110 and 138 were considered to be theobromine.

353

Therefore, we deduced that [M+H]+ at 183 and fragment ions at m/z 139, 111 were

ACS Paragon Plus Environment

Page 16 of 42

Page 17 of 42

Journal of Agricultural and Food Chemistry

354

[15N2]-theobromine. In the [15N2]-caffeine stable isotope-labelled tracing experiments,

355

[15N2]-theobromine, not theophylline, was detected in the metabolites when

356

[15N2]-caffeine was added (Figure 6).

357 358

TF regulation network in caffeine metabolism

359

Transcriptional control plays an important role in the regulation of secondary

360

metabolites in plants.37 In this study, a correlation analysis between the accumulation

361

of metabolites (caffeine, theobromine, and theophylline) and the expression of related

362

genes and TFs was performed (Figure 7). In total, 53 positively (Figure 7A) and 55

363

negatively (Figure 7B) regulated TF families were associated with unigenes involved

364

in these purine alkaloids metabolites. Most of the TFs can be categorized to bHLH,

365

WRKY, GRAS, MYB/MYB-related, and NAC families. Numerous key genes

366

involved in caffeine metabolism were associated with many TFs, demonstrating that

367

transcriptional regulation of caffeine metabolism was complicated. Pearson

368

correlation coefficients of >0.90 (between metabolites and gene expression, gene

369

expression and TF) were selected to construct a network diagram for positive

370

regulation. A significant positive correlation was revealed in the expression levels of

371

TEA022940.1 (encoded adenine phosphoribosyltransferase, APRT), TEA012581

372

(TCS), TEA028050.1 (TCS), and TEA002100.1 (encoded 5'-nucleotidase, 5'-NT) and

373

caffeine content. Furthermore, theophylline content has a positive correlation with the

374

expression of TEA029871.1 (encoded AMP deaminase, AMPD) and TEA010278.1

375

(AMPD). Theobromine content was shown a positive correlation with the expressions

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

376

of TEA012514.1 (SAMS) and TEA006735.1 (SAMS), TEA033814.1 (CYP1A2),

377

TEA028881.1 (5'-NT), TEA022785.1 (encoded ribokinase, RK), TEA015523.1

378

(APRT), and BGI_novel_G008251 (encoded allantoinase, ALN).

379 380

Contents of Amino acids and catechins, and gene expression profiles in theanine

381

and catechins metabolisms

382

The concentrations of amino acids in these Camellia leaves are listed in Table 3.

383

Among these plants, the content of total free amino acid of YK43 (19.25 mg/g) was

384

significantly higher than those of the other three Camellia plants, with LCT2 (6.45

385

mg/g) having the lowest content. The concentrations of glutamic acid (Glu),

386

glutamine (GluNH2) and theanine in YK43 and SCZ were higher than those in the

387

other two Camellia plants LCT1 and LCT2. The theanine content in YK43 (10.43

388

mg/g) was considerably higher than that in LCT1 (2.71 mg/g), LCT2 (1.43 mg/g). To

389

better understand theanine accumulation, the expression of key genes encoding

390

alanine aminotransferase (ALT), theanine synthetase (TS), glutamate dehydrogenase

391

(GDH), glutamate synthase (GOGAT), glutamine synthetase (GS) and alanine

392

decarboxylase (ADC) in the theanine biosynthetic pathway was investigated (Figure

393

8). The expression of TS (TEA015198.1) in YK43 was not altered significantly,

394

whereas the expression of GS (TEA015580.1) and ADC (TEA032991.1) was higher;

395

this probably enhanced the theanine biosynthesis process.

396

Catechins can be classified into epicatechin (EC), EC gallate (ECG), catechin (C),

397

gallocatechin (GC), epigalocatechin (EGC), epigallocatechin gallate (EGCG), and

ACS Paragon Plus Environment

Page 18 of 42

Page 19 of 42

Journal of Agricultural and Food Chemistry

398

gallocatechin gallate (GCG) in fresh tea leaves.50 The concentrations of catechins and

399

gallic acid in the leaves of different Camellia plants are listed in Table 4. Among

400

them, total catechins and gallic acid concentrations in YK43 (189.90 mg/g) and SCZ

401

(176.10 mg/g) were higher, whereas LCT1 (55.29 mg/g) had the lowest catechins and

402

gallic acid concentrations. The concentrations of EGCG, ECG, and EGC were higher

403

in YK43 and SCZ than in the other two Camellia plants. The key genes in catechins

404

biosynthetic pathway were also investigated (Figure 9).

405 406

Discussion

407

Caffeine is an important flavor compound in some soft drinks, particularly in coffee

408

and tea.51 According to the recent studies, caffeine content was less affected by

409

environments, and influenced by genotypic factors;52 caffeine content is different

410

between tea varieties.53 The caffeine distribution was investigated in 23 Camellia

411

species. And caffeine concentrations in young leaves were different.54 The catabolism

412

and biosynthesis of purine alkaloids in Camellia were also reported.55-56 Among these

413

Camellia plants, caffeine biosynthesis and caffeine synthase genes in tea plant and its

414

related species were extensively investigated.15-16 These Camellia resources provide

415

us a novel approach for understanding caffeine metabolism, not only biosynthesis but

416

also catabolism.

417

RNA-Seq is an efficient approach that is widely used for the transcripts

418

quantification and discovery.57-58 From the results of the enriched KEGG pathway,

419

DEGs involved in the caffeine metabolic pathways were found in the comparison

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

420

groups SCZ-VS-LCT1, SCZ-VS-LCT2, and SCZ-VS-YK43. Concentrations of

421

caffeine and theobromine among these Camellia plants were different. LCT1 and

422

LCT2 had the lower caffeine concentrations and relatively higher theobromine

423

concentrations than the other two plants, while SCZ and YK43 had the higher caffeine

424

concentration with a relatively lower theobromine concentration. The variation of

425

caffeine was closely related to the expression of related genes in metabolic pathways

426

(Table S3). The decreased caffeine content in the low-caffeine plants was probably

427

correlated with the high expressions of related genes in caffeine degradation pathway

428

(Figure 4, Table S4). In addition, many TFs could efficiently regulate the gene

429

expression. And the mechanism through which target genes are regulated by TFs,

430

further affecting the metabolite content, should be verified. The major degradation

431

pathways of caffeine have been clarified. From the results of the stable isotope tracer

432

experiment, a catabolic pathway from caffeine to theobromine in the leaves of

433

Camellia plants was confirmed and detected using LC-MS. It is different from the

434

previous publication that caffeine was first degraded to theophylline by Ashihara’s

435

laboratory using the radioisotope method.59

436

Theanine and catechins were important secondary metabolites in tea, and also

437

greatly contributed to the quality of tea. 60,61 The concentration of theanine was much

438

higher in YK43 than other three Camellia plants. The expression of GS

439

(TEA015580.1) was markedly upregulated in YK43. The sequences of GS and TS are

440

highly homologous,62 and in some plant species, theanine can be synthesized if

441

substrate exists. In addition, the high expression level of ADC in YK43 could enable

ACS Paragon Plus Environment

Page 20 of 42

Page 21 of 42

442

Journal of Agricultural and Food Chemistry

abundant theanine accumulation.

443

In this study, leaves form Camellia plants with low, normal, and high caffeine

444

concentrations were carried out RNA-Seq, characteristic compound determination,

445

qRT-PCR validation, stable isotope-labelled tracing, and network analysis of

446

TF-gene-metabolite for understanding the caffeine metabolism. Caffeine content

447

variation in these Camellia plants was greatly correlated with the gene expression in

448

caffeine metabolic pathway. For the results of stable isotope-labelled tracer, caffeine

449

was firstly degraded to theobromine, not theophylline, in the leaves of these Camellia

450

plants. The obtained results could provide new insights to understanding caffeine

451

metabolism in Camellia plants, and could also be applied to caffeine regulation and

452

breeding improvement in subsequent research.

453 454

Acknowledgments

455

This work was supported by the National Natural Science Foundation of China

456

(31570692, 31870679), the Changjiang Scholars and Innovative Research Team in

457

University (IRT_15R01), and the Natural Science Foundation of Anhui Province

458

(1608085QC60).

459 460

Authors’ contributions

461

WD and ZZ guided this research. LC prepared the samples. BZ, ML, JZ, and JH

462

analyzed the transcriptomic data and performed all experiments. BZ, WD and ZZ

463

drafted and revised the manuscript. All authors have read and approved the final

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

464

manuscript.

465 466 467

Ethical standards The authors declare compliance with the ethical standards.

468 469 470

Competing interests The authors declare no competing interests.

471 472 473 474

Data achieving statement The clean data will be available in NCBI SRA under the BioProject (PRJNA522339), if the manuscript is accepted for publication.

475 476

Legends for figures

477

Figure 1 A, Young leaves (one apical bud with two terminal leaves) of these Camellia

478

plants: low-caffeine tea 1 (LCT1), low-caffeine tea 2 (LCT2), Shuchazao (SCZ), and

479

Yunkang 43 (YK43); B, Venn diagram of DEGs in the three comparison groups: total

480

DEGs, upregulated DEGs, and downregulated DEGs

481 482

Figure 2 Enriched KEGG pathways of DEGs. The RichFactor in a pathway is the

483

ratio of the number of DEGs and all annotated genes. The larger enrichment factor is,

484

the higher enrichment degree is; the smaller Q value is, the more significant

485

enrichment is

ACS Paragon Plus Environment

Page 22 of 42

Page 23 of 42

Journal of Agricultural and Food Chemistry

486 487

Figure 3 Caffeine metabolic pathway. Red represents the pathway of caffeine

488

synthesis, blue represents the main pathway of caffeine degradation reported

489

previously,

490

KEGG. GMPS, GMP synthase; IMPDH, IMP dehydrogenase; GDA, guanine

491

deaminase; URE, urease; 3-NMT, 1-NMT, caffeine synthase (TCS); N-MeNase,

492

N-methyl nucleosidase; ALLC, allantoicase; Anase, adenosine nucleosidase; GK,

493

guanosine

494

phosphoribosyltransferase; NDM, N-demethylase

42

and green represents the pathway of caffeine degradation annotated in

kinase;

Gnase,

guanosine

nucleosidase;

GPRT,

guanine

495 496

Figure 4 Heat maps of key genes expressions in the caffeine metabolic pathway in

497

these Camellia plants. Each column represents the expression in leaves of different

498

Camellia plants and each row represents one gene. FPKM is shown on the logarithmic

499

scale. The red means higher expression and blue denotes lower expression. HPRT,

500

hypoxanthine phosphoribosyltransferase; ADSL, adenylosuccinate lyase

501 502

Figure 5 qRT-PCR validation for the key genes involved in caffeine metabolic

503

pathway. Relative expression levels from qRT-PCR and FPKM from the

504

transcriptome data are shown

505 506

Figure 6 Mass chromatograms of labelled theobromine (Tb) in leaves of Camellia

507

plants LCT1, LCT2, SCZ, and YK43, with [15N2]-caffeine. m/z 138, m/z 110 are

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 24 of 42

508

characteristic ions of non-labelled Tb in LC-MS analysis; m/z 139, m/z 111 are

509

characteristic ions of labelled Tb. The pattern of non-labelled Tb was used to assign

510

the fragments of labeled Tb

511 512

Figure 7 Transcription factor (TF) regulation network of caffeine’s metabolic pathway.

513

The red nodes represent metabolites, the yellow nodes represent genes, and green

514

nodes represent the enzyme annotated by each gene. TF is represented by pink and

515

blue nodes, and a line between TF and gene indicates the expressions are correlated.

516

A, positive regulation between metabolites and genes; B, negative regulation between

517

metabolites and genes

518 519

Figure 8 Expressions of the key genes in the pathway of theanine biosynthesis. A,

520

theanine biosynthetic pathway; B, heat map of key genes

521 522

Figure 9 Expressions of key genes in pathway of flavonoid biosynthesis. A, flavonoid

523

biosynthetic pathway; B, heat map of key genes. C4H, cinnamate 4-hydroxylase; PAL,

524

phenylalanine ammonia-lyase; 4CL, 4-coumarate-CoA ligase; CHI, chalcone

525

isomerase; ANR, anthocyanidin reductase; DFR, dihydroflavonol 4-reductase; CHS,

526

chalcone synthase; LAR, leucocyanidin reductase; F3H, flavanone 3-hydroxylase;

527

ANS,

528

acyltransferases

leucoanthocyanidin

oxidase;

SCPL,

serine

529 530

Supplementary material

ACS Paragon Plus Environment

carboxypeptidase-like

Page 25 of 42

Journal of Agricultural and Food Chemistry

531

Figure S1 DEG statistics. The X axis represents the comparison method for each

532

group; the Y axis represents the DEGs number. Red denotes upregulated DEGs and

533

blue represents downregulated DEGs

534

Figure S2 GO classification of DEGs in the three comparison groups: A,

535

SCZ-VS-LCT1; B, SCZ-VS-LCT2; C, SCZ-VS-YK43. The X axis represents the GO

536

terms; the Y axis represents the number of upregulated and downregulated DEGs

537

Table S1 The primers used for qRT-PCR verification

538

Table S2 DEGs in caffeine metabolism in the comparison groups of SCZ-VS-LCT1,

539

SCZ-VS-LCT2 and SCZ-VS-YK43

540

Table S3 The correlation coefficient of caffeine and gene expression levels

541

Table S4 The FPKM value of unigene in the caffeine degradation pathway, and the

542

correlation coefficient of caffeine and gene expression levels

543 544

Reference

545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560

1.

Fang, W. P.; Meinhardt, L. W.; Tan, H. W.; Zhou, L.; Mischke, S.; Zhang, D., Varietal identification

of tea (Camellia sinensis) using nanofluidic array of single nucleotide polymorphism (SNP) markers. Horticulture Research 2014, 1, 14035. 2.

Bordoni, A.; Hrelia, S.; Angeloni, C.; Giordano, E.; Guarnieri, C.; Caldarera, C. M.; Biagi, P. L.,

Green tea protection of hypoxia/reoxygenation injury in cultured cardiac cells. Journal of Nutritional Biochemistry 2002, 13 (2), 103-111. 3.

Nie, G.; Jin, C.; Cao, Y.; Shen, S.; Zhao, B., Distinct effects of tea catechins on

6-hydroxydopamine-induced apoptosis in PC12 cells. Archives of Biochemistry & Biophysics 2002, 397 (1), 84-90. 4.

Lambert, J. D.; Hong, J.; Yang, G. Y.; Liao, J.; Yang, C. S., Inhibition of carcinogenesis by

polyphenols: evidence from laboratory investigations. American Journal of Clinical Nutrition 2005, 81 (1), 284s-291s. 5.

Ogita, S.; Uefuji, H.; Yamaguchi, Y.; Koizumi, N.; Sano, H., Producing decaffeinated coffee plants.

Nature 2003, 423 (6942), 823. 6.

Denis, S.; Augur, C.; Marin, B.; Roussos, S., A new HPLC analytical method to study fungal

caffeine metabolism. Biotechnology Techniques 1998, 12 (5), 359-362.

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

561 562 563 564 565 566 567 568 569 570 571 572 573 574 575 576 577 578 579 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598 599 600 601 602 603 604

7.

Ashihara, H.; Crozier, A., Caffeine: a well known but little mentioned compound in plant science.

Trends in Plant Science 2001, 6 (9), 407-413. 8.

Kato, M.; Kanehara, T.; Shimizu, H.; Suzuki, T.; Gillies, F. M.; Crozier, A.; Ashihara, H., Caffeine

biosynthesis in young leaves of Camellia sinensis: In vitro studies on N-methyltransferase activity involved in the conversion of xanthosine to caffeine. Physiologia Plantarum 1996, 98 (3), 629–636. 9.

Ashihara, H.; Crozier, A., Biosynthesis and metabolism of caffeine and related purine alkaloids in

plants. Advances in Botanical Research 1999, 30 (08), 117-205. 10. Kato, M.; Mizuno, K.; Fujimura, T.; Iwama, M.; Irie, M.; Crozier, A.; Ashihara, H., Purification and characterization of caffeine synthase from tea leaves. Plant Physiology 1999, 120 (2), 579-86. 11. Kato, A.; Crozier, A.; Ashihara, H., Subcellular localization of the N-3 methyltransferase involved in caffeine biosynthes in tea. Phytochemistry 1998, 48 (5), 777-779. 12. Fujimori, N.; Ashihara, H., Adenine metabolism and the synthesis of purine alkaloids in flowers of Camellia. Phytochemistry 1990, 29 (11), 3513-3516. 13. Kato, M.; Mizuno, K.; Crozier, A.; Fujimura, T.; Ashihara, H., Caffeine synthase gene from tea leaves. Nature 2000, 406 (6799), 956-957. 14. Li, Y.; Ogita, S.; Keya, C. A.; Ashihara, H., Expression of caffeine biosynthesis genes in tea (Camellia sinensis). Zeitschrift Für Naturforschung C 2008, 63 (3-4), 267. 15. Jin, J. Q.; Yao, M. Z.; Ma, C. L.; Ma, J. Q.; Chen, L., Association mapping of caffeine content with TCS1 in tea plant and its related species. Plant Physiology & Biochemistry Ppb 2016, 105, 251-259. 16. Jin, J. Q.; Yao, M. Z.; Ma, C. L.; Ma, J. Q.; Chen, L., Natural allelic variations of TCS1 play a crucial role in caffeine biosynthesis of tea plant and its related species. Plant Physiology & Biochemistry 2016, 100, 18-26. 17. Jin, J. Q.; Chai, Y. F.; Liu, Y. F.; Zhang, J.; Yao, M. Z.; Chen, L., Hongyacha, a naturally caffeine-free tea plant from Fujian, China. Journal of Agricultural and Food Chemistry 2018, 66 (43), 11311-11319. 18. Ashihara, H.; Monteiro, A. M.; Moritz, T.; Gillies, F. M.; Crozier, A., Catabolism of caffeine and related purine alkaloids in leaves of Coffea arabica L. Planta 1996, 198 (3), 334-339. 19. Buerstmayr, H.; Lemmens, M., Breeding healthy cereals: genetic improvement of Fusarium resistance and consequences for mycotoxins. World Mycotoxin J 2015, 8 (5), 591-602. 20. Asano, Y.; Komeda, T.; Yamada, H., Microbial production of theobromine from caffeine. Bioscience, Biotechnology, and Biochemistry 2014, 57 (8), 1286-1289. 21. Brand, D.; Pandey, A.; Roussos, S.; Soccol, C. R., Biological detoxification of coffee husk by filamentous fungi using a solid state fermentation system. Enzyme and Microbial Technology 2000, 27 (1), 127-133. 22. Götz, S.; Garcíagómez, J. M.; Terol, J.; Williams, T. D.; Nagaraj, S. H.; Nueda, M. J.; Robles, M.; Talón, M.; Dopazo, J.; Conesa, A., High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Research 2008, 36 (10), 3420-3435. 23. Wei, C.; Yang, H.; Wang, S.; Zhao, J.; Liu, C.; Gao, L.; Xia, E.; Lu, Y.; Tai, Y.; She, G., Draft genome sequence of Camellia sinensis var. sinensis provides insights into the evolution of the tea genome and tea quality. Proceedings of the National Academy of Sciences of the United States of America 2018, 115 (18), 201719622. 24. Kim, D.; Langmead, B.; Salzberg, S. L., HISAT: a fast spliced aligner with low memory requirements. Nature Methods 2015, 12 (4), 357-360. 25. Kanehisa, M.; Goto, S.; Sato, Y.; Furumichi, M.; Tanabe, M., KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Research 2012, 40, 109-114.

ACS Paragon Plus Environment

Page 26 of 42

Page 27 of 42

605 606 607 608 609 610 611 612 613 614 615 616 617 618 619 620 621 622 623 624 625 626 627 628 629 630 631 632 633 634 635 636 637 638 639 640 641 642 643 644 645 646 647 648

Journal of Agricultural and Food Chemistry

26. Conesa, A.; Gotz, S.; Garciagomez, J. M.; Terol, J.; Talon, M.; Robles, M., Gene ontology database Blast2GO:a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 2005, 21(18), 3674–3676. 27. Quevillon, E.; Silventoinen, V.; Pillai, S.; Harte, N.; Mulder, N.; Apweiler, R.; Lopez, R., InterProScan: protein domains identifier. Nucleic Acids Research 2005, 33, 116-20. 28. Li, B.; Dewey, C. N., RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics 2011, 12 (1), 323-323. 29. Trapnell, C.; Williams, B. A.; Pertea, G.; Mortazavi, A.; Kwan, G.; Van Baren, M. J.; Salzberg, S. L.; Wold, B. J.; Pachter, L., Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nature Biotechnology 2010, 28 (5), 511-515. 30. Love, M. I.; Huber, W.; Anders, S., Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biology 2014, 15 (12), 550. 31. Audic, S.; Claverie, J. M., The significance of digital gene expression profiles. Genome Research 1997, 7 (10), 986-995. 32. Anders, S.; Huber, W., Differential expression analysis for sequence count data. Genome Biology 2010, 11(10), 106 33. Trapnell, C.; Hendrickson, D. G.; Sauvageau, M.; Goff, L.; Rinn, J. L.; Pachter, L., Differential analysis of gene regulation at transcript resolution with RNA-seq. Nature Biotechnology 2013, 31 (1), 46-53. 34. Tai, Y.; Wei, C.; Yang, H.; Zhang, L.; Chen, Q.; Deng, W.; Wei, S.; Zhang, J.; Fang, C.; Ho, C., Transcriptomic and phytochemical analysis of the biosynthesis of characteristic constituents in tea (Camellia sinensis) compared with oil tea (Camellia oleifera). BMC Plant Biology 2015, 15 (1), 190. 35. Wu, Z. J.; Tian, C.; Jiang, Q.; Li, X. H.; Zhuang, J., Selection of suitable reference genes for qRT-PCR normalization during leaf development and hormonal stimuli in tea plant (Camellia sinensis). Scientific Reports 2016, 6, 19748. 36. Zheng, X. Q.; Ashihara, H., Distribution, biosynthesis and function of purine and pyridine alkaloids in Coffea arabica seedlings. Plant Science 2004, 166 (3), 807-813. 37. Vom, E. D.; Kijne, J. W.; Memelink, J., Transcription factors controlling plant secondary metabolism: what regulates the regulators? Phytochemistry 2002, 61 (2), 107-114. 38. Ito, E.; Ashihara, H., Contribution Purine nucleotide biosynthesis de novo to the formation of caffeine in Young Tea (Camellia sinensis ) Leaves. Journal of Plant Physiology 1999, 154 (2), 145-151. 39. Jiang, X.; Liu, Y.; Li, W.; Zhao, L.; Meng, F.; Wang, Y.; Tan, H.; Yang, H.; Wei, C.; Wan, X.; Gao, L.; Xia, T., Tissue-specific, development-dependent phenolic compounds accumulation profile and gene expression pattern in tea plant [Camellia sinensis]. PLoS One 2013, 8 (4), e62315. 40. Wang, H. F.; Tsai, Y. S.; Lin, M. L.; Ou, A., Comparison of bioactive components in GABA tea and green tea produced in Taiwan. Food Chemistry 2006, 96 (4), 648-653. 41. Liao, Y.; Zeng, L.; Li, P.; Sun, T.; Wang, C.; Li, F.; Chen, Y.; Du, B.; Yang, Z., Influence of plant growth retardants on quality of codonopsis radix. Molecules 2017, 22 (10), 1655. 42. Sherlock, G., Gene Ontology: tool for the unification of biology. Canadian Institute of Food Science & Technology Journal 2009, 22 (4), 415. 43. P, M., Catabolism of caffeine in plants and microorganisms. Frontiers in Bioscience A Journal & Virtual Library 2004, 9 (1-3), 1348. 44. Ashihara, H.; Yokota, T.; Crozier, A., Biosynthesis and catabolism of purine alkaloids. Advances in

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

649 650 651 652 653 654 655 656 657 658 659 660 661 662 663 664 665 666 667 668 669 670 671 672 673 674 675 676 677 678 679 680 681 682 683 684 685 686 687 688 689 690 691

Botanical Research 2013, 68 (68), 111-138. 45. Deng, W. W.; Ashihara, H., Profiles of purine metabolism in leaves and roots of Camellia sinensis seedlings. Plant and Cell Physiology 2010, 51 (12), 2105-2118. 46. Boatright, J.; Negre, F.; Chen, X.; Kish, C. M.; Wood, B.; Peel, G.; Orlova, I.; Gang, D.; Dudareva, R. N., Understanding in vivo benzenoid metabolism in petunia petal tissue. Plant Physiology 2004, 135(4), 1993-2011. 47. Eisenreich, W.; Bacher, A.; Arigoni, D.; Rohdich, F., Biosynthesis of isoprenoids via the non-mevalonate pathway. Cellular & Molecular Life Sciences 2004, 61 (12), 1401. 48. Hiroshi, H.; Toshiyuki, O.; Haruka, I.; Kensuke, T.; Miwa, S.; Masakazu, H.; Naoharu, W., Functional characterization of aromatic amino acid aminotransferase involved in 2-phenylethanol biosynthesis in isolated rose petal protoplasts. Journal of Plant Physiology 2012, 169 (5), 444-451. 49. Long-Ze, L.; Pei, C.; Harnly, J. M., New phenolic components and chromatographic profiles of green and fermented teas. Journal of Agricultural & Food Chemistry 2008, 56 (17), 8130. 50. Liang, Y.; Ma, W.; Lu, J.; Ying, W., Comparison of chemical compositions of Ilex latifolia Thumb and Camellia sinensis L. Food Chem. 2001, 75 (3), 339-343. 51. Mohanpuria, P.; Kumar, V.; Yadav, S. K., Tea caffeine: Metabolism, functions, and reduction strategies. Food Science & Biotechnology 2010, 19 (2), 275-287. 52. Wei, K.; Wang, L.-Y.; Zhou, J.; He, W.; Zeng, J.-M.; Jiang, Y.-W.; Cheng, H., Comparison of catechins and purine alkaloids in albino and normal green tea cultivars (Camellia sinensis L.) by HPLC. Food Chem. 2012, 130 (3), 720-724. 53. Wang, L. Y.; Wei, K.; Jiang, Y. W.; Cheng, H.; Zhou, J.; He, W.; Zhang, C. C., Seasonal climate effects on flavanols and purine alkaloids of tea (Camellia sinensis L.). European Food Research & Technology 2011, 233 (6), 1049-1055. 54. Nagata, T. S., S., Difference in caffeine, flavanols and amino acids contents in leaves of cultivated species and hybrids in the genus Camellia. Japan Agricultural Research Quarterly 1986, 19, 276-280. 55. Ashihara, H.; Kubota, H., Biosynthesis of purine alkaloids in Camellia plants. Plant and cell physiology 1987, 28 (3), 535-539. 56. Kato, M.; Ashihara, H., Biosynthesis and catabolism of purine alkaloids in Camellia plants. Natural Product Communications 2008, 3 (9), 1429-1435. 57. Cole, T.; Williams, B. A.; Geo, P.; Ali, M.; Gordon, K.; Baren, M. J., Van; Salzberg, S. L.; Wold, B. J.; Lior, P., Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nature Biotechnology 2010, 28 (5), 511-515. 58. Mortazavi, A.; Williams, B. A.; Mccue, K.; Schaeffer, L.; Wold, B., Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nature Methods 2008, 5(7),621-628. 59. Ashihara, H.; Gillies, F. M.; Crozier, A., Metabolism of caffeine and related purine alkaloids in leaves of tea (Camellia sinensis L.). Plant & Cell Physiology 1997, 38 (4), 413-419. 60. Yamaguchi, S.; Ninomiya, K., Umami and food palatability. The Journal of Nutrition 2000, 130(4), 921S-926S. 61. Fei, G.; Guo, Y.; Pu, W.; Yu, W.; Ni, D., Transcriptional profiling of catechins biosynthesis genes during tea plant leaf development. Planta 2017, 246 (1), 1-14.

62. Li, C. F.; Zhu, Y.; Yu, Y.; Zhao, Q. Y.; Wang, S. J.; Wang, X. C.; Yao, M. Z.; Luo, D.; Li, X.; Chen, L., Global transcriptome and gene regulation network for secondary metabolite biosynthesis of tea plant (Camellia sinensis). BMC Genomics 2015, 16 (1), 560.

ACS Paragon Plus Environment

Page 28 of 42

Page 29 of 42

Journal of Agricultural and Food Chemistry

Table 1. Summary of sequencing reads Sample

Total Raw Reads (Mb)

Total Clean Reads (Mb)

Total Clean Bases (Gb)

Clean Reads Q20 (%)

Clean Reads Q30 (%)

Clean Reads Ratio (%)

LCT1_1 LCT1_2 LCT1_3 LCT2_1 LCT2_2 LCT2_3 SCZ_1 SCZ_2 SCZ_3 YK43_1 YK43_2 YK43_3

60.43 57.16 55.53 60.43 60.43 55.53 60.43 58.80 55.53 57.16 55.53 55.53

44.58 44.65 44.60 44.56 45.04 44.57 45.16 45.11 44.55 44.55 45.13 44.97

6.69 6.70 6.69 6.68 6.76 6.69 6.77 6.77 6.68 6.68 6.77 6.75

97.42 97.81 97.85 96.48 96.77 97.61 96.49 97.47 97.80 97.56 97.81 97.70

92.66 93.66 93.75 89.61 90.91 93.10 90.16 92.70 93.48 92.95 93.53 93.27

73.77 78.10 80.31 73.73 74.53 80.27 74.73 76.72 80.22 77.93 81.27 80.98

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 30 of 42

Table 2. Determination of purine alkaloids in the leaves of Camellia plants Compounds Purine alkaloids (μg/g) Xanthine (X) 1-Methylxanthine (1-mX) 3-Methylxantine (3- mX) Paraxanthine (Px) Theophylline (Tp) Purine alkaloids (mg/g) 7-Methylxantine (7- mX) Theobromine(Tb) Caffeine (Cf) Total purine alkaloids (mg/g)

C. sinensis Shuchazao (SCZ)

C. crassicolumna low caffeine tea 1 (LCT1)

C. crassicolumna low caffeine tea 2 (LCT2)

C. sinensis var. assamica YK43 (YK43)

0.82 ± 0.14 0.33 ± 0.05 n.d. 1.17 ± 0.11 4.77 ± 0.19

1.49 0.04

± ± n.d. 0.72 ± 0.29 ±

0.41** 0.00**

0.33 0.07

0.04** 0.06**

13.02 ± 2.26** 0.12 ± 0.03** tr. 0.73 ± 0.08** 2.46 ± 0.26**

0.05 ± 0.01 1.27 ± 0.21 34.17 ± 3.7

0.14 5.43 8.33

± ± ±

0.00** 0.39** 0.49**

0.26 ± 0.02** 7.7 ± 0.74** 0.59 ± 0.09**

0.22 3.50 42.82

± ± ±

0.04** 0.13** 0.85**

35.49 ± 3.79

15.02

±

2.33**

8.57 ± 0.81**

46.68

±

0.71**

± ± n.d. 1.11 ± 0.53 ±

0.05** 0.01** 0.19 0.10**

Values are expressed as μg/g or mg/g dry weight. n.d., not detected; tr., trace. The significance of the difference of contents between comparison groups (SCZ-VS-LCT1, SCZ-VS-LCT2, and SCZ-VS-YK43) is shown by *P < 0.05, **P < 0.01.

ACS Paragon Plus Environment

Page 31 of 42

Journal of Agricultural and Food Chemistry

Table 3. Determination of amino acids in the leaves of Camellia plants Compounds

C. sinensis Shuchazao (SCZ)

Main amino acids (μg/g) Asp 606.96 Thr 228 .00 Ser 681.06 Glu 2614.20 GluNH2 1125.21 Theanine 9482.10 Ala 180.14 Leu 5.34 Trp 1363.32 Arg 25.92 Pro 33.45 Total free amino 16.64 acids (mg/g)

C. crassicolumna low caffeine tea 1 (LCT1)

± ± ± ± ± ± ± ± ± ± ±

16.32 1.40 8.56 10.38 53.89 256.08 17.61 0.31 14.61 0.85 5.35

577.64 162.40 371.36 1455.40 994.31 2711.90 121.09 221.55 820.74 178.82 49.13

± ± ± ± ± ± ± ± ± ± ±

19.87 7.79* 20.64* 64.86 40.67 156.76** 4.47** 7.21** 23.79** 11.11** 4.48**

±

0.33**

8.61

±

0.33**

C. crassicolumna low caffeine tea 2 (LCT2) 680.96 192.92 282.57 1310.70 697.52 1428.50 127.17 1.39 718.06 5.36 55.15

± ± ± ± ± ± ± ± ± ± ±

13.26 7.33 11.74** 9.49** 29.05** 14.78** 4.01** 0.03 22.26** 0.49* 3.44**

6.45

±

0.08**

C.sinensis var. assamica YK43 (YK43) 864.21 211.76 587.91 2613.90 1841.20 10429.00 183.45 2.47 1059.10 90.21 48.58

± ± ± ± ± ± ± ± ± ± ±

120.15** 25.18 64.38* 192.28** 128.09** 589.95** 15.92 0.28 106.01** 10.21** 2.54**

19.25

±

1.23**

Values are expressed as μg/g or mg/g dry weight. The significance of the difference of contents between comparison groups (SCZ-VS-LCT1, SCZ-VS-LCT2, and SCZ-VS-YK43) is shown by *P < 0.05, **P < 0.01.

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 32 of 42

Table 4. Determination of catechins and gallic acid in the leaves of Camellia plants Compounds (A) Catechins (mg/g) Catechin (C) Gallocatechin (GC) Gallocatechin gallate (GCG) Epi-catechin (EC) Epi-gallocatechin (EGC) Epi-catechin-3-gallate (ECG) Epi-gallocatechin-3-gallate (EGCG) (B) Gallic acid (GA) (mg/g) Total catechins and gallic acid (mg/g)

C. sinensis C. crassicolumna C. crassicolumna Shuchazao low caffeine tea 1 low caffeine tea 2 (SCZ) (LCT1) (LCT2)                 3.74 ± 0.24 3.22 ± 0.37 4.59 ± 0.70 3.78 ± 0.08 n.d. n.d. 1.99 ± 0.00 2.67 ± 0.18 2.95 ± 0.18 16.99 ± 0.46 9.69 ± 0.81** 8.94 ± 1.49** 34.25 ± 0.70 10.32 ± 0.20** 13.64 ± 1.11** 33.53 ± 2.39 9.62 ± 1.16** 14.28 ± 0.63**

C.sinensis var. assamica YK43 (YK43)    6.54 ± 0.64** 4.14 ± 0.13** 6.14 ± 0.66** 29.5 ± 1.08** 18.35 ± 1.35** 42.31 ± 0.97**

80.91

±

4.29

19.59

±

1.64**

40.21

±

1.80**

82.48

±

3.06

0.34

±

0.03

0.19

±

0.01**

0.21

±

0.02**

0.43

±

0.03**

176.10

±

8.45

55.29

±

2.61**

85.75

±

5.43**

189.90

±

3.25*

Values are expressed as mg/g dry weight. n.d., not detected. The significance of the difference of contents between comparison groups (SCZ-VS-LCT1, SCZ-VS-LCT2, and SCZ-VS-YK43) is shown by *P < 0.05, **P < 0.01.

ACS Paragon Plus Environment

Page 33 of 42

Journal of Agricultural and Food Chemistry

Figure 1

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 34 of 42

Figure 2

ACS Paragon Plus Environment

Page 35 of 42

Journal of Agricultural and Food Chemistry

Figure 3

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 36 of 42

Figure 4

ACS Paragon Plus Environment

Page 37 of 42

Journal of Agricultural and Food Chemistry

Figure 5

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 38 of 42

Figure 6

ACS Paragon Plus Environment

Page 39 of 42

Journal of Agricultural and Food Chemistry

Figure 7

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 40 of 42

Figure 8

ACS Paragon Plus Environment

Page 41 of 42

Journal of Agricultural and Food Chemistry

Figure 9

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Graphic abstract (8.44cm wide and 4.73 cm tall)

ACS Paragon Plus Environment

Page 42 of 42