Quantitative Trait Loci Mapping for Theobromine and Caffeine

Nov 28, 2018 - Quantitative trait loci (QTL) analyses were carried out to identify genetic variants associated with theobromine and caffeine contents ...
0 downloads 0 Views 596KB Size
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