Subscriber access provided by University of Winnipeg Library
Omics Technologies Applied to Agriculture and Food
Quantitative Trait Loci Mapping for Theobromine and Caffeine Contents in Tea Plant (Camellia sinensis) Jian-Qiang Ma, Ji-Qiang Jin, Ming-Zhe Yao, Chun-Lei Ma, Yan-Xia Xu, Wan-Jun Hao, and Liang Chen J. Agric. Food Chem., Just Accepted Manuscript • DOI: 10.1021/acs.jafc.8b05355 • Publication Date (Web): 28 Nov 2018 Downloaded from http://pubs.acs.org on November 30, 2018
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 28
Journal of Agricultural and Food Chemistry
Quantitative Trait Loci Mapping for Theobromine and Caffeine Contents in Tea Plant (Camellia sinensis) Jian-Qiang Ma#, Ji-Qiang Jin#, Ming-Zhe Yao, Chun-Lei Ma, Yan-Xia Xu, WanJun Hao, Liang Chen* Key Laboratory of Tea Plant Biology and Resources Utilization, Ministry of Agriculture, Tea Research Institute of the Chinese Academy of Agricultural Sciences (TRICAAS), 9 South Meiling Road, Hangzhou 310008, China #Authors *E-mail:
contributed equally to this work.
[email protected]. Phone: +86 571 86652835. Fax: +86 571
86653866.
1
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
1
Abstract
2
Understanding the genetic basis of theobromine and caffeine accumulation in
3
the tea plant is important due to their contribution to tea flavor. In this study,
4
quantitative trait loci (QTL) analyses were carried out to identify genetic variants
5
associated with theobromine and caffeine contents and ratio, using a pseudo-
6
testcross population derived from an intervarietal cross between two varieties of
7
Camellia sinensis. A total of ten QTL controlling caffeine content (CAF),
8
theobromine content (TBR), sum of caffeine and theobromine (SCT), and caffeine-
9
to-theobromine ratio (CTR) were identified over four measurement years. The
10
major QTL controlling CAF, qCAF1, was mapped onto LG01 and validated across
11
years, explaining an average of 20.1% of the phenotypic variance. The other QTL
12
were detected in one or two years, and of them there were four, two and three for
13
TBR, SCT and CTR, respectively. The present results provide valuable information
14
for further fine mapping and cloning functional genes, and for genetic improvement
15
in tea plant.
16
Keywords: tea plant, caffeine, theobromine, quantitative trait loci
2
ACS Paragon Plus Environment
Page 2 of 28
Page 3 of 28
Journal of Agricultural and Food Chemistry
17
INTRODUCTION
18
Tea is one of the most popular non-alcoholic beverages, made by infusing of
19
processed tender shoots of tea plant (Camellia sinensis) and related species, which
20
all belong to the section Thea of the genus Camellia in the family Theaceace. These
21
plants have been utilized more than two thousand years in China where they were
22
originally used for traditional medicine and subsequently for food and drinking.
23
Recent research on phytochemical analysis of tea extracts have identified many
24
bioactive compounds, such as flavonoids, methylxanthines, L-theanine and
25
polysaccharides. 1-6 These versatile natural metabolites have not only contributed to
26
the unique taste and flavor of tea, but also proved to have powerful health-promoting
27
effects,7-10 which raises a growing interest and consumption of tea and relative
28
products containing tea constituents.
29
Methylxanthines are a group of purine-like alkaloids.11-13 Caffeine (1,3,7-
30
trimethylxanthine) and theobromine (3,7-dimethylxanthine) are the major
31
methylxanthines accumulated in the tea plant, contributing up to about 6% of the
32
dry weight of ‘two and a bud’ (i.e. the apical bud and two adjacent leaves harvested
33
for tea producing).14-16 Caffeine had a wide range of beneficial effects, such as
34
central nervous system stimulation, diuresis and cardiovascular function, while
35
theobromine has similar but less biological activity.2,17 Caffeine and theobromine
36
are bitter tastants which elicit bitterness through a family of TAS2R receptors,18-20
37
playing an important role on tea flavor.21-24 For example, the previous studies
38
showed caffeine was the main contributor to bitterness in green tea.25-30 Scharbert
39
and Hofmann31 demonstrated caffeine was as one of the key bitter compounds in
40
black tea. Pérez-Burillo et al.32 found a high correlation between caffeine and
3
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
41
bitterness in white tea. Furthermore, the interaction of caffeine and flavonoids
42
increased the intensity of bitterness and astringency in tea infusion.22,26,31 Therefore,
43
understanding the genetic basis of caffeine and theobromine accumulation in the tea
44
plant is important for breeding efforts to improve tea quality.33-36
45
Previous investigations have proved theobromine and caffeine are products of
46
purine metabolic pathway in tea plant.11-13 The main biosynthetic route is
47
xanthosine→ 7-methylxanthosine → 7-methylxanthine → theobromine → caffeine.
48
The three methylation steps are catalyzed by the N-methyltransferases (NMTs),
49
using S-adenosyl-L-methionine (SAM) as the methyl donor. The tea caffeine
50
synthase 1 (TCS1), a bifunctional NMT with both 1-N- and 3-N-methylation
51
activity, plays a crucial role in the biosynthesis of theobromine and caffeine.37 The
52
studies we conducted recently revealed that the tea germplasms containing specific
53
TCS1 alleles with low transcription level or enzymatic activity of encoded proteins
54
showed less caffeine accumulation than those with normal alleles.38,39 However, the
55
contents of theobromine and caffeine are a complicated quantitative trait affected
56
by both genetic and environmental factors, and to date, there is still limited
57
information defining the underlying genetic mechanism.
58
In this study, we used an F1 full-sib family to investigate the genetic basis of
59
variation for theobromine and caffeine accumulation in the tea plant through
60
quantitative trait loci (QTL) mapping of four associated traits, i.e. caffeine content
61
(CAF), theobromine content (TBR), sum of caffeine and theobromine (SCT), and
62
caffeine-to-theobromine ratio (CTR). The findings will help to elucidate the
63
molecular mechanism controlling theobromine and caffeine accumulation and to
64
accelerate the genetic improvement of tea plant.
4
ACS Paragon Plus Environment
Page 4 of 28
Page 5 of 28
Journal of Agricultural and Food Chemistry
65
MATERIALS AND METHODS
66
Mapping population. An F1 pseudo-testcross population was generated by
67
crossing C. sinensis var. sinensis ‘Yingshuang’ (‘YS’) and C. sinensis var.
68
pubilimba ‘Beiyue Danzhu’ (‘BD’). The population composed of 148 individuals
69
was genotyped using 406 simple sequence repeats (SSR) and 6042 SNP markers for
70
genetic map construction. Details were fully described in the previous studies.40
71
Phenotypic data collection. Field experiments were conducted at the
72
Experimental Station of TRICAAS as described previously.40,41 In briefly, each
73
parent and progeny were clonally propagated and grown in a randomized complete
74
block design with three replications. The ‘two and a bud’ of first flush in the spring
75
were harvested and dried for phytochemical analysis using high performance liquid
76
chromatography (HPLC) in 2010, 2011, 2014 and 2015.
77
The HPLC method was carried out following the international standard ISO
78
14502-2:2005 with minor modifications. Firstly, a 0.2 g of dry ‘two and a bud’
79
powder was extracted twice in the 70 °C water bath for 10 minutes using 5 mL 70%
80
methanol. The supernatant was collected and made up the volume to 10 mL with 70%
81
methanol. Then 1 mL of the supernatant was diluted to 5 mL with stabilizing
82
solution (10% v/v acetonitrile, 500 mg/mL EDTA and ascorbic acid), and
83
subsequently filtered through a hydrophilic nylon membrane filter with a 0.45 µm
84
pore size. Afterwards, the solution was detected using HPLC-UVD (Waters 2695-
85
2489) with a SunFire C18 column (5 µm, 4.6 mm × 250 mm, Milford, MA, USA ).
86
The mobile flow rate was 1.0 mL/min, and column oven temperature maintained at
87
35 °C. Solvents were (A) 1% v/v formic acid and (B) acetonitrile. A binary gradient
88
elution system was performed as follow: 93.5 to 85% A (0-16 min), 85 to 75% A
5
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
89
(16-20 min), 75.0 to 93.5% A (20-26 min), isocratic 93.5% A (26-30 min). Spectral
90
data were collected at 280 nm wavelength. The chromatographic peaks of caffeine
91
and theobromine were identified by comparing the retention times of reference
92
substances. Calibration curves for quantification were generated by plotting the
93
peak area against the known concentration of reference substance.
94
Phenotypic data analysis.
Statistical analyses of phenotypic data were
95
conducted using PASW Statistics version 18.0 (IBM SPSS Inc., Chicago, IL, USA).
96
The significant differences between the mean values of parental traits were analyzed
97
using t-test. The mean, standard deviation (SD), coefficient of variation (CV),
98
kurtosis and skewness in the F1 population were calculated for each trait. Pearson's
99
correlations were estimated using the combined datasets over years. Variance
100
components were calculated using multi-factor analysis of variance. The years of
101
phenotypic data collection were considered as different environments, and
102
heritability was estimated as variance among genotypes divided by the total of
103
genotypes, genotype-by-environment interaction and residual variances.
104
QTL analysis. Linkage map information was described in the previous
105
studies.40 QTL mapping was performed using MapQTL 6 with the combination of
106
interval mapping (IM) and restricted multiple QTL model (rMQM) methods.42
107
Logarithm of odds (LOD) threshold for QTL declaration was estimated for each
108
trait and year using a permutation test (1,000 replications). The location of
109
maximum LOD score exceeded the genome-wide significance threshold of 95%
110
somewhere on the linkage group was considered as the position of the QTL, and the
111
region with a LOD score greater than the threshold was considered as the confidence
112
interval. The QTL detected across years for each trait was assigned to be the same
6
ACS Paragon Plus Environment
Page 6 of 28
Page 7 of 28
Journal of Agricultural and Food Chemistry
113
stable QTL if confidence intervals overlapped, and QTL with over 10% of
114
phenotypic variance explained (PVE) was defined as the major QTL.
115
RESULTS AND DISCUSSION
116
Phenotypic variation. The phenotypic values of the traits measured of the
117
parents and F1 population over four years are shown in Table 1. The two parents
118
significantly differed in all traits investigated. The paternal parent ‘BD’ had higher
119
CAF, TBR and SCT, but lower CTR than the maternal parent ‘YS’. The overall
120
averages of CAF, TBR, SCT and CTR across years were 33.9 mg▪g-1, 6.7 mg▪g-1,
121
40.6 mg▪g-1 and 6.0 for ‘YS’, while 41.0 mg▪g-1, 10.4 mg▪g-1, 51.4 mg▪g-1 and 4.6
122
for ‘BD’, respectively. The previous investigation of the representative Chinese tea
123
germplasms revealed that there was a high level of variation in CAF and TBR, with
124
a range from 2.3 to 53.4 mg▪g-1 and 0.5 to 46.7 mg▪g-1, respectively.14 Meanwhile,
125
the three varieties of C. sinensis was significantly different in CAF, and generally
126
the varieties of C. sinensis var. pubilimba and var. assamica had a higher CAF than
127
var. sinensis.14 The maternal parent ‘YS’ was developed from the progeny by
128
crossing the sinensis and assamica varieties, and hence had an average level of CAF
129
and TBR. In comparison, the paternal parent ‘BD’, a germplasm of the pubilimba
130
variety, showed a slightly higher CAF and TBR.
131
A continuous distribution with transgressive bidirectional segregation was
132
observed among the F1 population for all traits, and which indicated quantitative
133
genetic architecture (Figure 1). The phenotypic values within the progeny across
134
years varied from 21.3 to 49.0 mg▪g-1 for CAF, 4.2 to 16.1 mg▪g-1 for TBR, 29.1 to
135
55.5 mg▪g-1 for SCT, and 1.7 to 25.4 for CTR, respectively (Table 1). The mean
136
values were slightly smaller than the corresponding mid-parent values for all traits,
7
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
137
with the exception of CTR. The population coefficient of variation (CV) for TBR
138
and CTR was considerably higher than that for CAF and SCT. The high level of
139
genetic variation among the progeny for the aforesaid traits suggested that
140
significant improvements can be achieved through hybridization and selection.
141
The Pearson correlations between each trait over years are presented in Table
142
2. CAF and TBR showed a significant negative correlation (r = -0.466, p < 0.001).
143
SCT was more highly correlated with CAF (r = 0.681, p < 0.001) than with TBR (r
144
= 0.330, p < 0.001). CTR showed a significant positive correlation with CAF (r =
145
0.550, p < 0.001), and a significant negative correlation with TBR (r = -0.860, p