Genome-wide Association Mapping of Polyphenol ... - ACS Publications

May 10, 2016 - A genome-wide association study (GWAS) was conducted for free (FP) and bound (BP) phenolic compounds and their antioxidant capacities...
0 downloads 0 Views 377KB Size
Subscriber access provided by ORTA DOGU TEKNIK UNIVERSITESI KUTUPHANESI

Article

Genome-wide association mapping of polyphenol contents and antioxidant capacity in whole grain rice Feifei Xu, Jinsong Bao, Tae-sung Kim, and Yong-Jin Park J. Agric. Food Chem., Just Accepted Manuscript • DOI: 10.1021/acs.jafc.6b01289 • Publication Date (Web): 10 May 2016 Downloaded from http://pubs.acs.org on May 14, 2016

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 free 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 accessible to all readers and 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.

Journal of Agricultural and Food Chemistry 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 31

Journal of Agricultural and Food Chemistry

1

Genome-wide association mapping of polyphenol contents and

2

antioxidant capacity in whole grain rice

3

Feifei Xu†, ‡, ||, Jinsong Bao†, ‡, ||, Tae-Sung Kim† and Yong-Jin Park*, †,

§

4 5



6

University, Yesan, Republic of Korea

7



8

Zhejiang University, Huajiachi Campus, Hangzhou, 310029, China

9

§

10

Department of Plant Resources, College of Industrial Science, Kongju National

Institute of Nuclear Agricultural Science, College of Agriculture and Biotechnology,

Legume Bio-Resource Center of Green Manure (LBRCGM), Kongju National

University, Yesan, Republic of Korea

11 12

|| These authors equally contribute to this work

13

* To whom correspondence should be addressed. E-mail:[email protected];

14

Tel: +82 41-330-1201; Fax: +82- 41-330-1209

15

1

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

16

ABSTRACT

17

Whole grains contain various bioactive phytochemicals including phenolic acids, and

18

consumption of whole grain may provide desirable health benefits and reduce the

19

risks of chronic diseases due to their antioxidant activities. In this study, we qualified

20

and quantified 9 bound phenolic compounds in 32 red and 88 white pericarp

21

accessions of rice. Genome-wide association study (GWAS) was conducted for free

22

(FP) and bound (BP) phenolic compounds and their antioxidant capacities with

23

high-quality single nucleotide polymorphisms (SNPs) in two colored grain panels and

24

the whole panel. Rc was detected for all FP and antioxidant capacities in the whole

25

panel. Three loci (chr.1:30970095, chr.6: 24392269 and chr.9: 6670223) for more than

26

five phenolic-related traits, two loci (chr.4: 34120529 and chr.11: 28947480) for more

27

than two FP-related traits, and one locus (chr.11: 23220681) for ferulic acid detected

28

in two panels were potentially new genes that are valuable for further gene cloning.

29

Overall, this study increases our understanding on the genetics of phenolic acid

30

biosynthesis in the phenylpropanoid pathway.

31

KEYWORDS: Polyphenol content, antioxidant capacity, phenolic acid, GWAS, rice

32

2

ACS Paragon Plus Environment

Page 2 of 31

Page 3 of 31

Journal of Agricultural and Food Chemistry

33

INTRODUCTION

34

Whole grain cereals are known as functional foods or nutraceuticals, which can

35

reduce the risk of coronary heart disease,1 insulin resistance, and type-2 diabetes2 due

36

to the antioxidant capacities of phytochemicals present in the grains. Therefore,

37

consumers are interested in healthy diets of whole grain foods, and breeders have

38

attempted to breed new varieties with high antioxidant activities.

39

Pigmented whole grain rice has several colors, ranging from white to red to black

40

(purple),3 which is determined based on the composition of anthocyanin pigments that

41

are ubiquitous throughout the plant kingdom.4 The majority of antioxidant compounds

42

are mainly distributed in bran (pericarp, seed coat, nucellus, and aleurone), followed

43

by the embryo, and a small amount is detected in the endosperm.1, 5-7 In general,

44

antioxidant capacity is strongly correlated with pericarp color; black and red rice

45

varieties exhibit higher antioxidant activities than white rice varieties due to the

46

presence of anthocyanins in black rice and proanthocyanidins in red rice.8-11 However,

47

some study showed that in red and black rice cultivars, anthocyanin pigments

48

contributed little to the total antioxidant capacity, and other phenolic compounds

49

contributed a lot to their antioxidant capacity.4 Besides anthocyanins and

50

proanthocyanidins, phenolic acids play a role in determining antioxidant activities in

51

white rice varieties. It seems that phenolic acids are mainly in the form of insoluble

52

forms, while flavonoids and anthocyanins are mainly in soluble forms.8, 12 Ferulic acid

53

(56~77% of total phenolic acids) is the most abundant phenolic acid that mainly

54

existed in the endosperm, bran, and whole grain, followed by p-coumaric (8~24%),

55

sinapic (2~12%), gallic (1~6%), protocatechuic (1~4%), p-hydroxybenzoic (1~2%),

56

vanillic (1%), and syringic acids (1%). While caffeic acid, chlorogenic acid, cinnamic

57

acid, and ellagic acid are minor constitutes that account for less than 1% of total 3

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 4 of 31

58

phenolic acids.8 Phenolic compounds and antioxidant properties in whole grain rice

59

vary with the seed development stages,13 varieties,14 and grain shape.15 Despite

60

detailed information on the structure, component, distribution, and dynamic changes

61

in polyphenols, only a few studies have been reported regarding the genetics of FP

62

and BP compounds, as well as the antioxidant capacities of whole grain.

63

Two genes affecting proanthocyanidin synthesis in red- and brown-colored rice

64

were elucidated.16, 17 The Rc gene encodes a basic helix-loop-helix (bHLH) protein,

65

and the recessive white rice allele rc has a 14-bp deletion within exon 6 that knocks

66

out the bHLH domain of the protein, resulting in the white color. Rd, encoding

67

dihydroflavonol-4-reductase, is a gene related to anthocyanin biosynthesis. Rc, in the

68

absence of Rd, produces brown seeds, whereas Rd alone has no phenotype. One gene

69

regulating anthocyanin synthesis in black-colored rice was characterized;18 the black

70

grain trait originated from ectopic expression of the Kala4 bHLH gene due to a

71

rearrangement in the promoter region. Both Rc and Kala4 genes activate upstream or

72

downstream flavonol biosynthesis genes to produce specific pigments. Kim et al.

73

investigated the expression of five anthocyanin biosynthetic genes and showed that

74

their expression was enhanced during seed maturation and correlated with ambient

75

temperature during seedling growth. Based on association mapping, Shao et al.

76

identified 41 marker loci for grain color, phenolic content, flavonoid content, and

77

antioxidant capacity, and confirmed that Ra and Rc genes are main-effect loci for rice

78

grain color and polyphenol traits.20 The RM316 locus, located on chromosome 9, has

79

a significant association with phenolic content. Shao et al. evaluated the genotype ×

80

environment interactions for polyphenols and antioxidant capacity of rice, 3 finding

81

that genotype dominantly controls polyphenol and antioxidant capacities with > 94%

82

phenotype variation. They also detected 23 quantitative trait loci (QTLs); two QTLs 4

ACS Paragon Plus Environment

19

Page 5 of 31

Journal of Agricultural and Food Chemistry

83

located close to Ra and Rc were detected in both years, and others that were sensitive

84

to the environment were detected only in one year. However, excluding Rc and Rd, all

85

QTLs for polyphenols were localized at a large interval, and high-resolution

86

resequencing data is required for the detection of new QTLs for polyphenols.

87

Genome-wide association studies (GWAS) is a comprehensive approach that can

88

systemically search for causal genetic variation inrice genome,21 and succeed in

89

detecting QTLs for grain weight,22 agronomic traits, 23 polyphenol contents,3 among

90

others. However, identification of QTLs for phenolic acids, such as ferulic acid,

91

p-coumaric acid, and sinapic acid, etc. has not been reported.

92

The objective of the present investigation was to (1) identify and quantify bound

93

phenolic acids in whole grain rice and (2) detect QTLs for TPC, FPC, TFC, and BPC

94

compounds, as well as their antioxidant capacities in three panels, including the red

95

grain panel, white grain panel, and whole panel (red grain panel plus white grain

96

panel).

97 98

MATERIALS AND METHODS

99

Plant Materials. A core collection of 137 accessions was established previously

100

based on SSR markers using a heuristic approach with PowerCore software.24,

101

Among them, 120 accessions, including 88 white pericarp (white grain panel) and 32

102

red pericarp accessions (red grain panel), were used in this study. The whole panel

103

was planted in 2014 during the growing season from late April to late September at

104

Kongju National University farm, Yeasan Gun, South Korea. The grains were

105

air-dried to a moisture content of about 12%, stored at room temperature for 3 months.

106

The samples were de-hulled into brown rice, and then ground and passed through a

107

100-mesh sieve to obtain whole grain flour. All the flour samples were stored at 4 °C 5

ACS Paragon Plus Environment

25

Journal of Agricultural and Food Chemistry

108

prior to extraction.

109 110

Extraction of Soluble-free Phenolic Compounds. The extraction of soluble-free

111

phenolic compounds in whole grain materials was performed according to the method

112

described by Yu et al. 26, with minor modifications. Briefly, 0.5 g of whole grain flour

113

was extracted twice with 5 mL of 80% methanol. After shaking at 250 rpm for 1 h

114

using an oscillator (HZ- 9210K Desktop frozen oscillator, Jiangsu, China), the

115

mixture was centrifuged at 3000rpmfor 10 min. The supernatants of each duplicate

116

extraction were collected and combined together, and then stored at −20°C until use.

117

Duplicate extractions were performed for each sample.

118 119

Extraction of Insoluble-bound Phenolic Compounds. Bound phenolics extraction

120

was performed as described by Yu et al.

121

Residue obtained from the crude extraction was washed twice with 20 mL of ddH2O

122

to eliminate organic solvent and then subjected to alkaline hydrolysis with 20mL of

123

4M NaOH on a shaker (HZ-9210K Desktop frozen oscillator, Jiangsu, China) (250

124

rpm) for 2 h. The resultant mixture was adjusted to a pH between 1.5 and 2.0 with 6 N

125

HCl, and then extracted with 40mLof ethyl acetate three times (20, 10, and 10mL).

126

The supernatant was collected and combined into an Erlenmeyer flask and then

127

evaporated to dryness at 35°C under vacuum using a rotary evaporator and water bath

128

(RE-2000A, Ya Rong Biochemistry Instrument Factory, Shanghai, China). The dried

129

residue was re-dissolved in 2 mL of 50% methanol and stored at -20°C. Prior to

130

HPLC analysis, they were filtered through a 0.45-µm syringe filter.

26

and Qiu et al.10, with minor modifications.

131 132

Determination of Total Phenolic Content (TPC). Folin-Ciocalteu assay was 6

ACS Paragon Plus Environment

Page 6 of 31

Page 7 of 31

Journal of Agricultural and Food Chemistry

133

performed to determine TPC in the soluble-free (FPC) and insoluble-bound (BPC)

134

fractions according to Shao et al.3 First, 10-fold-diluted Folin-Ciocalteu reagent and

135

appropriately diluted samples (red rice free phenolics: 6-fold-diluted; all bound

136

phenolics: 5-fold-diluted, and white rice free phenolics: not diluted) were prepared

137

prior to measurements. Briefly, 200 µL of sample extracts or standards solution was

138

mixed completely with 1.5 mL of Folin-Ciocalteu reagent and reacted for 5 min. Next,

139

1.5 mL of saturated sodium carbonate (75 g/L) was added to the mixture. After 2 h in

140

the dark at room temperature, the absorbance of the resulting blue color was recorded

141

using a spectrophotometer at 725 nm. Duplicate determinations were conducted for

142

each extract. The calibration curve was prepared using gallic acid solutions. Phenolic

143

content is expressed as milligrams of gallic acid equivalent per 100 g (mg GAE/100g)

144

of dry weight whole grain rice flour.

145 146

Determination of Total Flavonoid Content (TFC). Soluble-free phenolic extracts

147

were used to determine the TFC. Solution A (5% NaNO2), solution B

148

(10%AlCl3·6H2O), and 1M NaOH was prepared before measurements.3 First, 150 µL

149

of solution A was mixed completely with 1mL of appropriately diluted samples or

150

standards solution and reacted for 5 min. Second, 150 µL of solution B was mixed

151

completely with the above mixture and reacted for 5 min. Third, 1mL of 1M NaOH

152

solution was added to the mixture, after which 3 mL of ddH2O was immediately

153

added to the mixture. Finally, after 45 min in the dark at room temperature, the

154

absorbance of the resulting yellow color was recorded using a spectrophotometer at

155

510 nm. TFC is expressed as milligrams of catechin acid equivalent per 100 g of rice

156

flour (mg CAE/100 g).

157

7

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 8 of 31

158

Determination of ABTS·+(2, 2-azino-bis-(3-ehylbenzothiazoline-6-sulphonic acid)

159

diammonium salt) Radical Scavenging Activity. The antioxidant capacities of rice

160

extracts

161

spectrophotometer based on the improved ABTS radical cation color decolorization

162

assay27 with minor modifications.3 The working solution of ABTS was prepared by

163

mixing 7 mM ABTS and 2.45 mM potassium persulfate (v/v, 2:1) and reacted in the

164

dark at room temperature for 12~16 h before use. The ABTS•+ working solution was

165

diluted in methanol to adjust the absorbance to 0.700 ± 0.020 at 734 nm. Then, 3.9 mL

166

of ABTS•+ working solution was completely mixed with 100 µL of sample or

167

standards and reacted in the dark at room temperature for 45 min. Subsequently, the

168

absorbance was recorded using a spectrophotometer at 734 nm. Results are expressed

169

as µmol of Trolox equivalent antioxidant capacity per 100 g of dry weight rice flour

170

(µmol TE/100 g).

(free

phenolic

and

bound

phenolic)

were

determined

using

a

171 172

Determination of DPPH (2, 2-diphenyl-1-picrylhydrazyl) Radical Scavenging

173

Activity. Free phenolic and bound phenolic antioxidant capacities were also

174

determined using the DPPH assay.3 Briefly, 3 mL of 100 µM DPPH solution was

175

completely mixed with 200 µL of sample or standard solution and reacted in the dark

176

at room temperature for 45 min. The absorbance was recorded using a

177

spectrophotometer at 517 nm. The DPPH scavenging activity (%) of both samples and

178

standard (Trolox) was calculated, as follows:

179 180 181

DPPH (%)= (1-Asample/Acontrol) ×100 The results are expressed as Trolox equivalent DPPH radical scavenging activity per 100 g of sample (µmol TE/100 g).

182

8

ACS Paragon Plus Environment

Page 9 of 31

Journal of Agricultural and Food Chemistry

183

Identification and Quantification of Phenolic Acids with High Performance

184

Liquid Chromatography (HPLC). Identification and quantification of bound

185

phenolic compounds were performed on an HPLC (Waters 2695, Milford, MA)

186

equipped with a photodiode array (PDA) detector (Waters 2998) and an autosampler

187

(Waters 717 plus) according to the method described by Yu et al. 26 and Shen et al. 28.

188

A column with the following specifications, 150 mm × 4.6 mm, 5 µm C18 110 A

189

column (Phenomenex, Torrance, CA), was used for HPLC separation. The

190

temperature was set at 35°C and the injection volume was 10 µL using an auto

191

sampler. The mobile phase consisted of A (0.1% acetic acid in water) and B (0.1%

192

acetic acid in methanol), and the flow rate was 0.9 mL/min. A 70 min linear gradient

193

was set as follows: 0~11 min, 9~14% B, 11–14 min, 14~15% B, 14~17 min, 15% B,

194

17~24 min, 15~16.5% B, 24~28 min, 16.5~19% B, 28~30 min, 19~25% B, 30~36

195

min, 25~ 26% B, 36~38 min, 26~28% B, 38~41 min, 28~35% B, 41~46 min, 35~40%

196

B, 46~48 min, 40~48% B, 48~53 min, 48~53% B, 53~65 min, 53~70% B, 65~66 min,

197

70~9% B, 66~70 min, and 9% B. Phenolic acids were detected at a wavelength of 280

198

nm. A total of 11 phenolic acid standards were used in this study (Figure S1),

199

corresponding to gallic acid, protocatechuic, p-hydroxybenzoic acid, vanillic acid,

200

caffeic acid, syringic acid, p-coumaric acid, ferulic acid, sinapic acid, iso-ferulic acid,

201

and o-coumaric, respectively. The contents of phenolic acids were quantified based on

202

the standard curves of the corresponding phenolic acids at the wavelength of 280 nm.

203 204

Statistical Analysis and Genome-wide Association Mapping. All statistical

205

analyses were performed with SAS program version 9.1 (SAS Institute Inc, Cary, NC,

206

USA). Proc Corr was used to examine Pearson correlation between different traits,

207

and the soluble-free phenolic-related traits and insoluble-bound phenolic-related traits 9

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

208

were performed separately. Analysis of variance (ANOVA) was performed using the

209

general linear model procedure (Proc glm), and Duncan’s new multiple range test was

210

used to determine significant differences between red and white grain rice panels.

211

Genome resequencing data of the 120 accessions (whole panel) were collected

212

from the 295 rice whole genome resequencing data, which were obtained previously

213

using HiSeq 2000 and HiSeq 2500.29 To control for false-positives, the P+K model

214

performed based on the mixed linear model (MLM) was used for genome-wide

215

association mapping of phenolic acid-related traits in three panels by GAPIT.30 SNP

216

sites with the lowest P value in the peak region (P