Diversity of Kale (Brassica oleracea var. sabellica ... - ACS Publications

Mar 30, 2016 - Diversity of Kale (Brassica oleracea var. sabellica): Glucosinolate. Content and Phylogenetic Relationships. Christoph Hahn,*,†. Anja...
1 downloads 0 Views 2MB Size
Subscriber access provided by NORTHERN KENTUCKY UNIV STEELY

Article

The diversity of kale (Brassica oleracea var. sabellica): Glucosinolate content and phylogenetic relationships Christoph Hahn, Anja Müller, Nikolai Kuhnert, and Dirk Albach J. Agric. Food Chem., Just Accepted Manuscript • DOI: 10.1021/acs.jafc.6b01000 • Publication Date (Web): 30 Mar 2016 Downloaded from http://pubs.acs.org on April 5, 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 37

Journal of Agricultural and Food Chemistry

The diversity of kale (Brassica oleracea var. sabellica): Glucosinolate content and phylogenetic relationships

Christoph Hahn * 1 , Anja Müller 2 , Nikolai Kuhnert 2 , Dirk Albach 1

1

Carl von Ossietzky University Oldenburg, Institute for Biology and Environmental

Sciences, Carl von Ossietzky Str. 9-11, 26111 Oldenburg, Germany

2

Jacobs University Bremen, School of Engineering and Science, Campus Ring 8,

Research III, 28759 Bremen, Germany

* corresponding author: E-mail: [email protected] Phone: 0049-441-798-3343

1 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

1

Abstract

2

Recently, kale has become popular due to nutritive components beneficial for human

3

health. It is an important source of phytochemicals such as glucosinolates that trigger

4

associated cancer-preventive activity. However, nutritional value varies among

5

glucosinolates and among cultivars. Here, we start a systematic determination of the

6

content of five glucosinolates in 25 kale varieties and 11 non-kale Brassica oleracea

7

cultivars by HPLC-DAD-ESI-MSn and compare the profiles with results from the

8

analysis of SNPs derived from a KASP genotyping assay. Our results demonstrate that

9

the glucosinolate levels differ markedly among varieties of different origin. The

10

comparison of the phytochemical data with phylogenetic relationships revealed that the

11

common name kale refers to at least three different groups. German, American and

12

Italian kales differ morphologically and phytochemically. Landraces do not show

13

outstanding glucosinolate levels. Our results demonstrate the diversity of kale and the

14

importance of preserving a broad genepool for future breeding purposes.

Page 2 of 37

15 16 17

Keywords: kale, cabbage, glucosinolates, HPLC-ESI-MS, phylogenetic relationships,

18

KASP assay, SNPs

2 ACS Paragon Plus Environment

Page 3 of 37

Journal of Agricultural and Food Chemistry

19

Introduction

20

Brassica oleracea L. is economically one of the most important species of the

21

Brassicaceae family because cabbage is a widely consumed vegetable all over the

22

world. Worldwide production is 70 million tons.1 The major reason for this success is

23

large variation in morphology and use by humans of this vegetable. The morphological

24

variation ranges from white and red cabbage over broccoli, cauliflower, kale, kohlrabi,

25

Brussels sprout to lesser known varieties. All domesticated cabbage varieties are

26

believed to have originated from wild cabbage (Brassica oleracea var. oleracea) native

27

to different coastal habitats at the Mediterranean and Atlantic.2 The presently known

28

cabbage cultivars are likely to have originated independently in different regions of

29

Europe from various wild forms. Romans already cultivated several different kinds of

30

cabbage.3 Evidence for this hypothesis is given by the study of glucosinolates, for

31

example.4 Song et al. suggested that ancient forms of Brassica oleracea may be

32

thousand head kale and Chinese kale / kai-lan.5 The most recent study on the phylogeny

33

of Brassica reported B. oleracea to be phylogenetically close to B. juncea (L.) Czern.

34

(Indian mustard), B. incana Ten. (Mediterranean mustard), B. rapa L. (turnip rape) and

35

B. napus L. (rapeseed, hybrid of wild cabbage and turnip rape), as well as B. cretica

36

Lam. and B. montana DC. (other wild cabbage types).6 Surprisingly little is, however,

37

known about relationships of varieties within B. oleracea. Louarn et al. demonstrated

38

polyphyly of white cabbage using AFLP (Amplified Fragment-length Polymorphism)

39

but did not resolve relationships among varieties.7 In contrast, Mei et al. resolved

40

broccoli as early branching in B. oleracea and sister-group relationships of kohlrabi and

41

savoy cabbage, as well as white cabbage and Chinese kale but sampling with B.

42

oleracea was too small to allow further inferences.8 Including all varieties is a huge task

43

given the enormous diversity in the species.

3 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 4 of 37

44

Here, we focus on kale (B. oleracea var. sabellica L.), which is a valuable fresh

45

vegetable during wintertime given its tolerance against frost. It is mostly assumed to

46

have originated from wild cabbage occurring on cliffs and rocky coasts along the

47

Atlantic but this is still under dispute.2 Typical German kale differs from other cabbages

48

in their curly leaves. However, the term kale is rather widely used for different

49

morphological types of large-leafed B. oleracea varying in leaf morphology and

50

sometimes described as different varieties, for example. In general, kale is characterized

51

by high levels of various nutrients and further constitutional metabolites like

52

glucosinolates, flavonoids and carotenoids.9 In this study we investigate different kale

53

varieties from various origins and countries (like Northern Germany, Italy, USA) that

54

differ in traits like growth height, shape, leaf color and curling. However, what is not

55

known is whether these geographical origins also form cohesive phytochemical and

56

phylogenetic groups. The main question is, therefore, whether kale constitutes a

57

homogenous group or whether regional groups of kale can be distinguished

58

phytochemically and genetically.

59 60

The phytochemical analysis focusses on five frequently occurring glucosinolates in the

61

leaves in late autumn. Glucosinolates are a large group of anionic, hydrophilic, sulphur-

62

containing secondary plant metabolites. They make a significant contribution to the

63

specific spicy cabbage flavor and bitter taste in Brassica crops.10 More than 100

64

different glucosinolates are known11 and about 20 are predominant in Brassica sp.12

65

They are derived from amino acids and can be grouped into aliphatic, aromatic and

66

indole glucosinolates. The endogenous enzyme myrosinase (thioglucoside

67

glucohydrolase) catalyzes the hydrolysis of glucosinolates to yield glucose, sulfate and

68

aglycones. The aglycones then decompose to breakdown products like isothiocyanates,

4 ACS Paragon Plus Environment

Page 5 of 37

Journal of Agricultural and Food Chemistry

69

thiocyanates, indoles and organic cyanides (nitriles) – depending on the substrate, pH

70

and availability of certain ions.13-15 Several isothiocyanates are found to exhibit

71

anticancerogenic effects against various types of cancer, such as breast, lung and colon

72

cancer.16-17 Compounds are bioavailable and show interference in early stages of cancer

73

development by inducing detoxification enzymes (phase II enzymes), inhibiting

74

activation enzymes (phase I enzymes) and triggering of cell cycle arrest or apoptosis.18-

75

21

76

cultivars is therefore an important aspect of Brassica breeding. We, here, investigate

77

quantitatively the occurrence of five glucosinolates, i.e., the aliphatics gluconapin,

78

progoitrin and glucoraphanin, the aromatic gluconasturtiin and the indole glucobrassicin

79

(Table 1, Fig. 1). Most of these glucosinolates lead to metabolites that are associated

80

with beneficial anticancerogenic properties (especially glucoraphanin22). Additionally,

81

progoitrin and the aromatic gluconasturtiin contribute to the bitterness of cabbage.13, 23

82

Progoitrin in particular is associated with the detrimental effect of potentially causing

83

goiter (through the hydrolytic product oxazolidine-2-thione)24. It is generally assumed

84

that bitter compounds are eliminated by breeding from commercially cultured plants.

85

Therefore, we hypothesized that landraces have increased levels of these compounds

86

compared with commercial cultivars. Glucosinolate profiles are known to differ among

87

Brassica oleracea cultivars.25 Information on a broad sample of kale cultivars and

88

determining specific glucosinolate profiles for certain groups of cultivars may therefore

89

be beneficial for future breeding and choice of kale types for improved nutrition.

90

Whereas glucosinolate levels are known to be influenced by various factors such as life

91

cycle, temperature and fertilization, common origin is still considered a crucial aspect of

92

phytochemical profiles.12

The occurrence and differences in specific glucosinolates among various Brassica

93

5 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 6 of 37

94

For investigating relationships between different kales and their relation to other

95

cabbage cultivars we made use of the single-plex SNP genotyping platform KASPTM

96

(Kompetitive Allele Specific Polymerase Chain Reaction (PCR)). Genetic variation is

97

measured on the basis of Single Nucleotide Polymorphisms (SNPs). They are useful for

98

investigating genetic diversity and discriminating between cultivars and have been

99

commonly applied in research on Brassica crops such as B. napus (e.g. Dalton-Morgan

100

et al 26). SNPs are one of the most common type of genetic variation in eukaryotic

101

genomes and most commonly found in the DNA between genes. The advantages of this

102

assay are low average genotyping errors, cheap genotyping costs, as well as scalable

103

flexibility in applications.27

104

Here, we use SNP and phytochemical data with the aim to identify glucosinolate

105

profiles specific to certain kale varieties. We hypothesize that different groups of kale

106

have evolved in different parts of the world differing in their ability to synthesize

107

different glucosinolates. Further, similarities between phytochemical profiles and

108

genetic data allow to infer the relationships of kale varieties, especially the origin of the

109

American kale varieties.

110 111

Material and Methods

112

Plant Material

113

In March 2014, seeds of 25 varieties of kale (Table 2) as well as 11 non-kale Brassica

114

oleracea cultivars were seeded in a greenhouse under natural daylight. When plants

115

reached the 5-leaf stage, plants were transplanted in a randomized scheme into a field at

116

the University of Oldenburg botanical garden. Water and fertilizer were applied

117

according to standard cultural practices. Leaves were harvested in August/September

118

2014.

6 ACS Paragon Plus Environment

Page 7 of 37

Journal of Agricultural and Food Chemistry

119 120

Quantification of the Glucosinolate Content

121

Extraction from Plant Material. For each kale variety, 10 g fresh material of the

122

youngest fully developed leaves were weighed (excluding the thick midnerves), treated

123

at 120 °C for 2h15 (to inactivate the enzyme myrosinase), ground to a fine powder using

124

a mixer mill, transferred to a screwtop vial and mixed with 8 mL of 70% v/v

125

methanol/water (of HPLC grade). The sealed vials were sonicated and then stored at

126

-20 °C until further analysis. Two individual plants per variety were sampled.

127 128

Analysis. Kale leaf extracts were obtained after thermal deactivation of myrosinase as

129

aqueous methanolic extracts. Extracts were profiled using an adapted liquid

130

chromatography–mass spectrometry method (LC-MS) from the literature.28, 29 Five

131

glucosinolates (i.e., glucobrassicin, glucoraphanin, gluconasturtiin, gluconapin and

132

progoitrin; Fig. 1) were detected and quantified with high sensitivity by HPLC-qTOF-

133

MS with negative ion electrospray ionization (ESI). A Bruker impact HD qTOF mass

134

spectrometer (Bruker Daltonik GmbH, Bremen, Germany) coupled to an Agilent 1260

135

Series HPLC system was used (Agilent Technologies, Santa Clara CA, USA),

136

consisting of: Online Degaser G1322A, Bin Pump G1312B, High Performance

137

Autosampler HiP-ALS G1367E, Thermostatted column compartment G1316A, Diode

138

Array Detector DAD-VL G1315D. Before analysis, 2 mL of the extracted samples were

139

filtered (0.45 μm pore size) into 2 mL screwtop vials. The HPLC-column was a 150 mm

140

x 2 mm (2.8 μm particle size) Pursuit XRs Ultra 2.8 Diphenyl column (Agilent

141

Technologies, Santa Clara CA, USA). The mobile phase was 0.005% formic acid

142

(HCOOH) with an acetonitrile/water gradient (CH3CN) over 75 min at a flow rate of 0.5

143

mL/min. The gradient profile increased from 5 to 75% linearly in 45 min followed by a 7 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 8 of 37

144

return to 5% and 30 min isocratic to re-equilibrate. The eluate absorbance was

145

monitored by an UV-detector at 229 nm. MS-scans were performed in the negative ion

146

mode on the glucosinolate [M-H]- ions and spectra were recorded between m/z (mass-

147

to-charge ratio) 100 and 700 (capillary voltage: 4500 V, nebulizer: 1.8 bar, dry gas

148

temperature: 200 °C, dry gas flow: 9 L/min). For quantification eight point calibration

149

curves (see Supplementary Figure 1) were obtained for each glucosinolate from

150

authentic reference substances (analytical standard glucosinolates from PhytoLab

151

GmbH & Co. KG, Vestenbergsgreuth, Germany). Quantification was carried out in the

152

negative ion mode on extracted ion chromatograms of each respective glucosinolate m/z

153

value with isolation width of 0.002 Da. All Pearson correlation coefficients R2 were

154

found to be above 0.99 and repeat injections yielding a relative standard deviation

155

(RSD) below 5%. The lowest quantifiable amount is indicated for each glucosinolate by

156

the signal-to-noise ratio (S/N) (see Table 3). Mass calibration was carried out using a

157

0.01 mM sodium formate solution injected prior to each chromatographic run using an

158

enhanced quadratic calibration algorithm. Spiking of the kale extract with authentic

159

reference compounds was carried out with a RSD below 5% indicating that no

160

appreciable matrix effect was operating.

161 162

Data Evaluation. For data evaluation, ESI Compass 1.3 Data Analysis software was

163

used (Version 4.0.234.0, Bruker Daltonik GmbH, Bremen, Germany). Statistical

164

analyses were conducted using the software R (v3.0.3, R Core Team30), specifically

165

Analysis of Variance (ANOVA), Tukey’s HSD test and Pearson’s product-moment

166

correlation analysis.

167 168

Phylogenetic Analyses Using KASP-Assay

8 ACS Paragon Plus Environment

Page 9 of 37

Journal of Agricultural and Food Chemistry

169

DNA-Extraction and KASP Assay Analysis. For each variety, one intact, fully

170

developed leaf was cut, placed in a plastic bag with silica beads inside and stored until it

171

was completely dried. DNA was extracted from leaves using the innuPREP Plant DNA

172

Kit (Analytik Jena AG, Jena, Germany) and afterwards stored at -20 °C until further

173

analysis. The subsequent KASPTM assay, developed by KBioscience/LGC Genomics

174

(Hoddesdon, Hertfordshire, UK), was carried out by TraitGenetics (TraitGenetics

175

GmbH, Gatersleben, Stadt Seeland, Germany). Each sample was analyzed with 100

176

KASP markers well distributed over the genome of Brassica oleracea. KASP assays

177

rely on allele-specific oligo extension and Fluorescence Resonance Energy Transfer

178

(FRET) for signal generation.31 In the KASP cycle, an allele-specific primer matches a

179

target SNP and, with a common reverse primer, amplifies the target region. In a second

180

PCR round, the reverse primer binds to the previously generated product and thus

181

makes a complementary copy. In a third PCR round, a fluorescence-labeled oligo (FAM

182

or HEX from FRET cassette) binds to the amplicon from the previous round.32, 33

183 184

Data Analysis. Phylogenetic relationships were inferred using Bayesian Inference (BI)

185

and Neighbor-Net methods. The BI analysis was conducted using MrBayes (v3.2.4,

186

Ronquist et al.34) sampling across the entire GTR (general time reversible) model space,

187

including gamma. Two independent runs of 10,000,000 generations were completed

188

with four chains each (three heated, one cold). Trees were sampled every 100

189

generations, the first 25% of generations were discarded as burn-in. A majority-rule

190

consensus of the remaining trees from the two runs was produced with posterior

191

probabilities. Siberian kale (B. napus var. pabularia DC.) was used as outgroup taxon.

192

Compared to a phylogenetic tree, phylogenetic networks can often better illustrate the

193

evolutionary history, e.g. in the case of hybridization. They display the extent of

9 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 10 of 37

194

conflicting signals in the data.35 Neighbor-Net is a distance based method for

195

constructing phylogenetic networks, based on the Neighbor Joining (NJ) algorithm.36 A

196

Neighbor-Net was constructed using SplitsTree (v4.13.1, Huson & Bryant37) with

197

uncorrected P-distances and ambiguous states were matched. A detailed description of

198

the underlying mechanism of the method is given by Bryant & Moulton.38

199 200

Results

201

Quantification of the Glucosinolate Content

202

The five glucosinolates previously described in kale gluconapin, progoitrin,

203

glucoraphanin, gluconasturtiin and glucobrassicin were identified in the chromatogram

204

based on their high resolution m/z value of their [M-H] pseudomolecular ion, which all

205

showed errors below 5 ppm. Additionally the structures were confirmed using a tandem

206

MS measurement with all glucosinolates showing characteristic fragment ions at m/z

207

259.01 (C6H12O9S). Based on this characteristic fragmentation mechanism a search for

208

previously unreported glucosinolates was carried out by construction of extracted ion

209

chromatograms in All MSn mode (see Supplementary Figure 2). Isothiocyanate products

210

(e.g., erucin, sulforaphane) are absent from the LC-MS chromatograms (see

211

Supplementary Figure 4).

212 213

The content of five glucosinolates were analyzed in different kale and cabbage samples,

214

with two individual plants being measured for each variety with analytical duplicate

215

measurements. The resulting levels, received with the help of reference standard

216

measurements, are summarized in the Supplementary Table 1. The gluconapin content

217

varied between 0.005 mg/100g fresh weight (FW) (Redbor) and 19.6 mg/100g FW

218

(Vates) (Fig. 2). The American varieties (especially Georgia Southern, Champion and

10 ACS Paragon Plus Environment

Page 11 of 37

Journal of Agricultural and Food Chemistry

219

Vates) showed significantly higher gluconapin contents compared to all others (p
0.05). Additionally, the general lack of coherence

408

among German landraces and German commercial cultivars indicates multiple origins

409

of German commercial cultivars from various Northern German landraces. Holterfehn

410

and Jellen as well as Buss Bunde, Rote Palme and Diepholzer each have very similar

411

glucosinolate contents. However, looking at the phylogenetic results, none of them are

412

placed in a sister-group relationship. A lot of these Northern varieties have their native

413

region around the Northern German city “Leer” (Holterfehn, Jellen, Buss Bunde,

414

Lammertsfehn, Schatteburg, Neuefehn). In the BI phylogeny (Supplementary Figure 3)

415

they cluster nearly perfectly together. In the Neighbor-Net (Fig. 4) we find them also

416

quite near to each other.

417

18 ACS Paragon Plus Environment

Page 19 of 37

Journal of Agricultural and Food Chemistry

418

Flower-Sprout, a hybrid between kale and Brussels sprouts, is a newly developed hybrid

419

cultivar originating in England (Gärtner Pötschke GmbH, Kaarst, Germany). It is not

420

known which kale variety was used as parent for hybridization, but with reference to

421

Flower-Sprout’s red colored leaves it is imaginable that the parental variety also had red

422

leaves (like Redbor or Rote Palme).

423

Siberian kale was used as outgroup taxon in phylogenetic analyses because it belongs to

424

Brassica napus ssp. pabularia. The leaves of Siberian resemble those of kale a bit,

425

because it is believed that Brassica napus originated from hybridization between turnip

426

(Brassica rapa) and kale. Carlson et al. reported a high progoitrin and glucobrassicin

427

content in this variety.48 We instead found a high gluconasturtiin content, which is

428

consistent with studies on Brassica rapa.49, 50

429 430

Implications for health-oriented breeding

431

While we have assembled data on different aspects of the glucosinolate content and

432

relationships of kale and demonstrated differences between regional origins, a number

433

of mechanisms may explain why glucosinolate levels differ between varieties. Several

434

sets of genes regulate the content of specific glucosinolates in different cultivars. Also

435

additional factors may influence the glucosinolate composition (soil properties, drought,

436

temperature, light etc.12). Since some of the glucosinolate metabolites exhibit cancer

437

chemopreventive or even anticancerogenic properties and others have rather detrimental

438

effects and can decrease food quality through bitter taste (due to high progoitrin levels,

439

for example23), there has to be a trade-off between beneficial and unfavorable

440

properties. A combination of different metabolites in plants should be aimed at. Most

441

potent breakdown products in terms of anticancerogenic effects are the ones of

442

glucoraphanin, like the isothiocyanate sulforaphane.13, 51 Therefore kale varieties with

19 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 20 of 37

443

elevated glucoraphanin levels are desirable as functional food with added health benefit.

444

Frostara and Black Tuscany as well as red cabbage and Romanesco produce much

445

glucoraphanin. In comparison typical broccoli varieties contain around 10 mg

446

glucoraphanin per 100 g FW.28 As our quantitative data show glucoraphanin levels in

447

some kale varieties can be as high as ten times that of an average broccoli, which is

448

currently the gold standard with respect to provision of glucoraphanin dietary exposure.

449

In contrast, the American varieties, Palmizio and Lerchenzungen have rather low

450

amounts. Frostara at the same time produces comparatively much progoitrin, and so do

451

Siberian kale, the American kale, red and wild cabbage. These ones may therefore taste

452

more bitter. The aromatic gluconasturtiin is mostly present in low amounts in Brassica

453

oleracea, Siberian kale (B. napus) instead contains much higher volumes. The wide

454

range of variability in glucosinolates content among kale can be used as an information

455

for breeding and developing new varieties with enhanced benefits.

456 457

Abbreviations Used

458

AFLP, amplified fragment-length polymorphism; ANOVA, analysis of variance; BI,

459

bayesian inference; Da, Dalton; ESI, electrospray ionization; FRET, fluorescence

460

resonance energy transfer; FW, fresh weight; GTR, general time reversible; HD, high

461

definition; HPLC, high-pressure liquid chromatography; KASP, kompetitive allele

462

specific polymerase chain reaction; MS, mass spectrometry; m/z, mass-to-charge ratio;

463

NJ, neighbor joining; PCR, polymerase chain reaction; ppm, parts per million; qTOF,

464

quadrupole time of flight; RSD, relative standard deviation; S/N, signal-to-noise ratio;

465

SNP, single nucleotide polymorphism

466 467

20 ACS Paragon Plus Environment

Page 21 of 37

Journal of Agricultural and Food Chemistry

468

Acknowledgements

469

We thank Vera Mageney for help with DNA extraction, Reinhard Lühring for seeds of

470

local Northern German kale and Holger Ihler and Gertrud Genz for assistance in the

471

cultivation of the plants.

472

Funding / Notes

473

The authors declare no competing financial interest.

474 475

Supporting Information Available

476

Supplementary Figure 1: calibration curves of glucoraphanin and glucobrassicin

477

Supplementary Figure 2: base peak chromatogram (BPC) and extracted ion

478

chromatograms (EIC) of kale variety Frostara

479

Supplementary Figure 3: Bayesian Inference phylogeny

480

Supplementary Figure 4: extracted ion chromatogram (EIC) of sulforaphane and

481 482

glucoraphanin in kale variety Buss Bunde Supplementary Table 1: glucosinolate values of analyzed kale and cabbage samples

21 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 22 of 37

483

References

484

(1) FAOSTAT. Food and Agriculture Organization of the United Nations – Statistics

485

Division; http://faostat3.fao.org/download/Q/QC/E (accessed December 15, 2015).

486

(2) Körber-Grohne, U. Nutzpflanzen in Deutschland: Kulturgeschichte und Biologie, 2nd

487

ed.; Konrad Theiss: Stuttgart, Germany, 1988.

488

(3) Wittstein, G. C. Die Naturgeschichte des Cajus Plinius Secundus; Gressner &

489

Schramm: Leipzig, Germany, 1881.

490

(4) Mithen, R. F.; Lewis, B. G.; Heaney, R. K.; Fenwick, G. R. Glucosinolates of wild

491

and cultivated Brassica species. Phytochemistry 1987, 26, 1969-1973.

492

(5) Song, K. M.; Osborn, T. C.; Williams, P. H. Brassica taxonomy based on nuclear

493

restriction fragment length polymorphisms (RFLPs). 2. Preliminary analysis of

494

subspecies within B. rapa (syn. campestris) and B. oleracea. Theor. Appl. Genet. 1988,

495

76, 593-600.

496

(6) Arias, T.; Pires, C. A fully resolved chloroplast phylogeny of the brassica crops and

497

wild relatives (Brassicaceae: Brassiceae): Novel clades and potential taxonomic

498

implications. Taxon 2012, 61 (5), 980-988.

499

(7) Louarn S.; Torp, A. M.; Holme, I. B.; Andersen, S. B.; Jensen, B. D. Database derived

500

microsatellite markers (SSRs) for cultivar differentiation in Brassica oleracea. Genet.

501

Resour. Crop Evol. 2007, 54, 1717-1725.

502

(8) Mei, J.; Li, Q.; Yang, X.; Qian, L.; Liu, L.; Yin, J.; Frauen, M.; Li, J.; Qian, W.

503

Genomic relationships between wild and cultivated Brassica oleracea L. with emphasis

504

on the origination of cultivated crops. Genet. Resour. Crop Evol. 2010, 57, 687-692.

22 ACS Paragon Plus Environment

Page 23 of 37

Journal of Agricultural and Food Chemistry

505

(9) Ferioli, F.; Giambanelli, E.; D’Antuono, F.; Costa, H. S.; Albuquerque, T. G.; Silva,

506

A. S.; Hayran, O.; Kocaoglu, B. Comparison of leafy kale populations from Italy,

507

Portugal, and Turkey for their bioactive compound content: phenolics,

508

glucosinolates, carotenoids, and chlorophylls. J. Sci. Food Agric. 2013, 93, 3478-3489.

509

(10) MacLeod, A. J. Volatile flavour compounds of the Cruciferae. In The biology and

510

chemistry of the Cruciferae; Vaughan, J. G.; MacLeod, A. J.; Jones, B. M. G., Eds.;

511

Academic Press: London, UK, 1976; pp 307-330.

512

(11) Kjaer, A.; Skrydstrup, T. Selenoglucosinolates: Synthesis and Enzymatic

513

Hydrolysis. Acta Chem. Scand., Ser. B 1987, 41, 29-33.

514

(12) Fahey, J. W.; Zalcmann, A. T.; Talalay, P. The chemical diversity and distribution

515

of glucosinolates and isothiocyanates among plants. Phytochemistry 2001, 56, 5-51.

516

(13) Kushad, M. M.; Brown, A. F.; Kurilich, A. C.; Juvik, J. A.; Klein, B. P.; Wallig,

517

M. A.; Jeffery, E. H. Variation of Glucosinolates in Vegetable Crops of Brassica

518

oleracea. J. Agric. Food Chem. 1999, 47, 1541-1548.

519

(14) Verkerk, R.; Dekker, M.; Jongen, W. M. F. Glucosinolates: natural toxicants in

520

food. In Natural Toxicants in Food; Watson, D., Ed.; Academic Press: New York, USA,

521

1998; pp 29-53.

522

(15) Song, L.; Morrison, J. J.; Botting, N. P.; Thornalley, P. J. Analysis of

523

glucosinolates, isothiocyanates, and amine degradation products in vegetable extracts

524

and blood plasma by LC-MS/MS. Anal. Biochem. 2005, 347, 234-243.

525

(16) Devi, J. R.; Thangam, E. B. Extraction and separation of glucosinolates from

526

Brassica oleracea var. rubra. Adv. Biol. Res. 2010, 4 (6), 309-313.

527

(17) Beecher, C. W. W. Cancer preventive properties of varieties of Brassica oleracea:

528

a review. Am. J. Clin. Nutr. 1994, 59 (5), 1166 -1170. 23 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 24 of 37

529

(18) Mithen, R.; Faulkner, K.; Magrath, R.; Rose, P.; Williamson, G.; Marquez, J.

530

Development of isothiocyanate-enriched broccoli, and its enhanced ability to induce

531

phase 2 detoxification enzymes in mammalian cells. Theor. Appl. Genet. 2003, 106,

532

727-734.

533

(19) Clarke, J. D.; Dashwood, R. H.; Ho, E. Multi-targeted prevention of cancer by

534

sulforaphane. Cancer Lett. 2008, 269, 291-304.

535

(20) Yoxall, V.; Kentish, P.; Coldham, N.; Kuhnert, N.; Sauer, M. J.; Ioannides, C.

536

Modulation of hepatic cytochromes P450 and phase II enzymes by dietary doses of

537

sulforaphane in rats: Implications for its chemopreventive activity. Int. J. Cancer 2005,

538

117, 356-362.

539

(21) Hanlon, N.; Coldham, N.; Gielbert, A.; Kuhnert, N.; Sauer, M. J.; King, L. J.;

540

Ioannides, C. The absolute bioavailability and dose dependent pharmacokinetic

541

behaviour of dietary doses of the chemopreventive isothiocyanate sulforaphane in rat.

542

Br. J. Nutr. 2008, 99, 559-564.

543

(22) Iori, R.; Bernardi, R.; Gueyrard, D.; Rollin, P.; Palmieri, S. Formation of

544

Glucoraphanin by Chemoselective Oxidation of Natural Glucoerucin: A

545

Chemoenzymatic Route to Sulforaphane. Bioorg. Med. Chem. Lett. 1999, 9, 1047-1048.

546

(23) van Doorn, H. E.; van der Kruk, G. C.; van Holst, G.-J.; Raaijmakers-Ruijs, N. C.

547

M. E.; Postma, E.; Groeneweg, B.; Jongen, W. H. F. The glucosinolates sinigrin and

548

progoitrin are important determinants for taste preferences and bitterness of Brussels

549

sprouts. J. Sci. Food Agric. 1998, 78, 30-38.

550

(24) Liu, Z.; Hirani, A. H.; McVetty, P. B. E.; Daayf, F.; Quiros, C. F.; Li, G. Reducing

551

progoitrin and enriching glucoraphanin in Brassica napus seeds through silencing of the

552

GSL-ALK gene family. Plant. Mol. Biol. 2012, 79, 179-189. 24 ACS Paragon Plus Environment

Page 25 of 37

Journal of Agricultural and Food Chemistry

553

(25) Nilsson, J.; Olsson, K.; Engqvist, G.; Ekvall, J.; Olsson, M.; Nyman, M.; Akesson,

554

B. Variation in the content of glucosinolates, hydroxycinnamic acids, carotenoids, total

555

antioxidant capacity and low-molecular-weight carbohydrates in Brassica vegetables. J.

556

Sci. Food Agric. 2006, 86, 528-538.

557

(26) Dalton-Morgan, J.; Hayward, A.; Alamery, S.; Tollenaere, R.; Mason, A. S.;

558

Campbell, E.; Patel, D.; Lorenc, M. T.; Yi, B.; Long, Y.; Meng, J.; Raman, R.; Raman,

559

H.; Lawley, C.; Edwards, D.; Batley, J. A high-throughput SNP array in the

560

amphidiploid species Brassica napus shows diversity in resistance genes. Funct. Integr.

561

Genomics 2014, 14 (4), 643-655.

562

(27) Semagn, K.; Babu, R.; Heame, S.; Olsen, M. Single nucleotide polymorphism

563

genotyping using Kompetitive Allele Specific PCR (KASP): overview of the

564

technology and its application in crop improvement. Mol. Breed. 2014, 33, 1-14.

565

(28) Song, L.; Thornalley, P. J. Effect of storage, processing and cooking on

566

glucosinolate content of Brassica vegetables. Food Chem. Toxicol. 2007, 45 (2), 216-

567

224.

568

(29) Mellon, F. A.; Bennett, R. N.; Holst, B.; Williamson, G. Intact Glucosinolate

569

Analysis in Plant Extracts by Programmed Cone Voltage Electrospray LC/MS:

570

Performance and Comparison with LC/MS/MS Methods. Anal. Biochem. 2002, 306, 83-

571

91.

572

(30) R Core Team. R: A Language and Environment for Statistical Computing. R

573

Foundation for Statistical Computing, Vienna, Austria; http://www.R-project.org

574

(accessed December 15, 2015).

575

(31) Kumpatla, S. P.; Buyyarapu, R.; Abdurakhmonov, I. Y.; Mammadov, J. A.

576

Genomics-Assisted Plant Breeding in the 21st Century: Technological Advances and 25 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 26 of 37

577

Progress. In Plant Breeding; Abdurakhmonov, I., Ed.; InTech: Rijeka, Croatia, 2012; pp

578

131-184.

579

(32) LGC Genomics. KASP genotyping technology. LGC Limited, Teddington, UK;

580

http://www.lgcgroup.com/kasp/#.VOtBuCzKmkg (accessed December 15, 2015).

581

(33) He, C.; Holme, J.; Anthony, J. SNP Genotyping: The KASP Assay. Methods Mol.

582

Biol. 2014, 1145, 75-86.

583

(34) Ronquist, F.; Teslenko, M.; van der Mark, P.; Ayres, D. L.; Darling, A.; Höhna, S.;

584

Larget, B.; Liu, L.; Suchard, M. A.; Huelsenbeck, J. P. MRBAYES 3.2: Efficient

585

Bayesian phylogenetic inference and model selection across a large model space. Syst.

586

Biol. 2012, 61 (3), 539-542.

587

(35) Penny, D.; Hendy, M. D.; Holland, B. R. Phylogenetic Networks. In Handbook of

588

Statistical Genetics, Vol. 1; Balding, D. J.; Bishop, M.; Cannings, C., Ed.; Wiley & Sons:

589

Chichester, UK, 2008; 3rd ed., pp 520-523.

590

(36) Saitou, N.; Nei, M. The Neighbor-Joining method: a new method for

591

reconstruction of phylogenetic trees. Mol. Biol. Evol. 1987, 4, 406-425.

592

(37) Huson, D. H.; Bryant, D. Application of Phylogenetic Networks in Evolutionary

593

Studies. Mol. Biol. Evol. 2006, 23 (2), 254-267.

594

(38) Bryant, D.; Moulton, V. Neighbor-Net: An Agglomerative Method for the

595

Construction of Phylogenetic Networks. Mol. Biol. Evol. 2004, 21 (2), 255-265.

596

(39) Ciska, E.; Martyniak-Przybyszewska, B.; Kozlowska, H. Content of Glucosinolates

597

in Cruciferous Vegetables Grown at the Same Site for Two Years under Different

598

Climatic Conditions. J. Agric. Food Chem. 2000, 48, 2862-2867.

26 ACS Paragon Plus Environment

Page 27 of 37

Journal of Agricultural and Food Chemistry

599

(40) Kliebenstein, D. J.; Kroymann, J.; Brown, P.; Figuth, A.; Pedersen, D.;

600

Gershenzon, J.; Mitchell-Olds, T. Genetic Control of Natural Variation in Arabidopsis

601

Glucosinolate Accumulation. Plant Physiol. 2001, 126, 811-825.

602

(41) Mewis, I.; Tokuhisa, J. G.; Schultz, J. C.; Appel, H. M.; Ulrichs, C.; Gershenzon, J.

603

Gene expression and glucosinolate accumulation in Arabidopsis thaliana in response to

604

generalist and specialist herbivores of different feeding guilds and the role of defense

605

signaling pathways. Phytochemistry 2006, 67, 2450-2462.

606

(42) Falk, K. L.; Kästner, J.; Bodenhausen, N.; Schramm, K.M.; Paetz, C.; Vassao, D.

607

G.; Reichelt, M.; von Knorre, D.; Bergelson, J.; Erb, M.; Gershenzon, J.; Meldau, S.

608

The role of glucosinolates and the jasmonic acid pathway in resistance of Arabidopsis

609

thaliana against molluscan herbivores. Mol. Ecol. 2014, 23, 1188-1203.

610

(43) Cartea, M. E.; Velasco, P.; Obregon, S.; Padilla, G.; de Haro, A. Seasonal

611

Variation in Glucosinolate Content in Brassica oleracea Crops Grown in Northwestern

612

Spain. Phytochemistry 2008, 69 (2), 403-410.

613

(44) Mikkelsen, M. D.; Petersen, B. L.; Glawischnig, E.; Jensen, A. B.; Andreasson, E.;

614

Halkier, B. A. Modulation of CYP79 genes and glucosinolate profiles in Arabidopsis by

615

defense signaling pathways. Plant Physiol. 2003, 131, 298-308.

616

(45) Gladis, T.; Hammer, K. Nomenclatural notes on the Brassica oleracea-group.

617

Genet. Resour. Crop Evol. 2001, 48, 7-11.

618

(46) Huxley, A. J. The new Royal Horticultural Society dictionary of gardening;

619

Stockton Press: New York, USA, 1992.

620

(47) Velasco, P.; Francisco, M.; Moreno, D. A.; Ferreres, F.; Garcia-Viguera, C.;

621

Cartea, M. E. Phytochemical Fingerprinting of Vegetable Brassica oleracea and

27 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 28 of 37

622

Brassica napus by Simultaneous Identification of Glucosinolates and Phenolics.

623

Phytochem. Anal. 2011, 22, 144-152.

624

(48) Carlson, D. G.; Daxenbichler, M. E.; VanEtten, C. H.; Kwolek, W. F.; Williams, P.

625

H. Glucosinolates in Crucifer Vegetables: Broccoli, Brussels Sprouts, Cauliflower,

626

Collards, Kale, Mustard Greens, and Kohlrabi. J. Am. Soc. Hortic. Sci. 1987, 112 (1),

627

173-178.

628

(49) Yang, B.; Quiros, C. F. Survey of glucosinolate variation in leaves of Brassica

629

rapa crops. Genet. Resour. Crop Evol. 2010, 57, 1079-1089.

630

(50) Kang, J. Y.; Ibrahim, K. E.; Juvik, J. A.; Kim, D. H.; Kang, W. J. Genetic and

631

Environmental Variation of Glucosinolate Content in Chinese Cabbage. Hortic. Sci.

632

2006, 41 (6), 1382-1385.

633

(51) Fahey, J. W.; Stephenson, K. K.; Talalay, P. Glucosinolates, myrosinase, and

634

isothiocyanates: three reasons for eating Brassica vegetables. ACS Symp. Ser. 1998,

635

701, 16-22.

636

(52) The Royal Society of Chemistry. ChemSpider. Chemical structure database.

637

Burlington House, London; http://www.chemspider.com/ (accessed December 15,

638

2015).

28 ACS Paragon Plus Environment

Page 29 of 37

Journal of Agricultural and Food Chemistry

Figure captions

Figure 1. Chemical structures of glucosinolates quantified. The structures were obtained by using ChemDraw software (v15.0.0.106, PerkinElmer Informatics, Massachusetts, USA).

Figure 2. Glucosinolate content (mg/100g FW) of different kale and cabbage varieties. Two individual plants per variety were measured (black and grey bars). Note the different scales.

Figure 3. Differences in glucosinolate content (mg/100g FW) between the groups of kale and cabbage varieties: (A) gluconapin content, (B) progoitrin content, (C) glucoraphanin content, (D) gluconasturtiin content, (E) glucobrassicin content. Note the different scales. Different letters above each boxplot indicate significant differences in a post-hoc Tukey’s HSD test (p < 0.05).

Figure 4. Neighbor-Net of different kale and cabbage varieties based on the KASPanalysis. Sample names are color coded according to their origin. Numbers beside branches are Bayesian posterior probabilities (support >70) for the respective clade. Leaf margin curliness of the varieties is indicated.

29 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 30 of 37

Tables

Table 1. The Glucosinolates Analyzed In This Study, Listed With Their Relevant Characteristics. Molecular Formulas Were Obtained From The Royal Society Of Chemistry52. group

gluconapin

m/z experimental [M-H] 372.0399

m/z theoretical [M-H] 372.0418

error [ppm] 4.9

molecular formula C11 H19 N1 O9 S2

retention time [min] 7.0

progoitrin

388.0368

388.0366

2.5

C11 H19 N1 O10 S2

4.4

glucoraphanin

436.0396

436.0400

3.5

C12 H23 N1 O10 S3

4.3

aromatic

gluconasturtiin

422.0578

422.0574

1.7

C15 H21 N1 O9 S2

17.0

indole

glucobrassicin

447.0496

447.0537

2.6

C16 H20 N2 O9 S2

15.8

aliphatic

trivial name

30 ACS Paragon Plus Environment

Page 31 of 37

Journal of Agricultural and Food Chemistry

Table 2. Summary Of The Kale And Cabbage Varieties Used In This Study. The Assignment To The Analyses Is Displayed As Well As The Origin Of The Seeds (e.g., Company The Varieties Have Been Obtained From). included in glucosinolate analysis

variety

non-kale cabbage cultivars

local Northern German kale

American varieties

Italian varieties

German kale varieties (B. oleracea convar. acephala var. sabellica)

Lerchenzungen

x

Vitessa

included in phylogenetic analyses x x

Vitessa

x

Lage Fijngekrulde

x

Niedriger Grüner Krauser

x

Reflex F1

x

x

Redbor F1

x

x

Frostara

x

x

Winnetou F1

x

x

Flower-Sprout "Petit Posy Mix" (Brussels sprouts × kale) Black Tuscany

x

x

x

x

Palmizio Senza Testa

x

x

Negro Romano

x

x

Georgia Southern

x

x

Morris Heading

x

x

Champion

x

x

Vates

x

x

Buss Bunde

x

x

Niedriger von Rosenweide

x

x

Rote Palme

x

x

Rote Palme Holterfehn

x

x

Schatteburg

x

x

Jellen

x

x

Neuefehn Hainwatjes

x

x

Lammertsfehn

x

x

Diepholzer Dickstrunk

x

x

Red cabbage "Roodkop" B. oleracea convar. capitata var. rubra Savoy "Vertas 2" B. oleracea convar. capitata var. sabauda Butter cabbage B. oleracea var. costata Kohlrabi "Delikateß blauer" (blue) B. oleracea convar. acephala var. gongylodes Romanesco broccoli B. oleracea convar. botrytis var. botrytis Brussels sprouts "Hilds Ideal" B. oleracea convar. gemmifera var. gemmifera Marrow-stem kale "Westfälischer Furchenkohl" B. oleracea convar. acephala var. medullosa Marrow-stem kale "Walking Stick" B. oleracea convar. acephala var. medullosa Galician cabbage B. oleracea convar. acephala Wild cabbage „Helgoländer“ B. oleracea L. Siberian kale B. napus ssp. pabularia

x x x x x x x

x

x

x

x

x

x

x

x

x

seed origin Kiepenkerl, Bruno Nebelung GmbH, Everswinkel, Germany N.L.Chrestensen Erfurter Samen- und Pflanzenzucht GmbH, Erfurt, Germany Centrala Nasienna, Saatbau Polska Sp. z o.o., Środa Śląska, Poland N.L.Chrestensen Erfurter Samen- und Pflanzenzucht GmbH, Erfurt, Germany PNOS (Polskie Nasiennictwo Ogrodnictwo Szkółkarstwo Sp. z o.o.), Ożarów Mazowiecki, Poland Dürr Samen S. Schwenk e.K., Reutlingen, Germany Kiepenkerl, Bruno Nebelung GmbH, Everswinkel, Germany Kiepenkerl, Bruno Nebelung GmbH, Everswinkel, Germany Kiepenkerl, Bruno Nebelung GmbH, Everswinkel, Germany Gärtner Pötschke GmbH, Kaarst, Germany Thompson & Morgan (UK) Ltd., Ipswich, UK Thompson & Morgan (UK) Ltd., Ipswich, UK Reinhard Lühring, private farmer, Rhauderfehn, Germany

Sustainable Seed Company, Covelo CA, USA

Reinhard Lühring, private farmer, Rhauderfehn, Germany

Gartenland GmbH Aschersleben, Essen, Germany toom, Centor-Warenhandels-GmbH, Köln, Germany Magic Garden Seeds, Andreas Fái-Pozsár, Regensburg, Germany N.L.Chrestensen Erfurter Samen- und Pflanzenzucht GmbH, Erfurt, Germany Magic Garden Seeds, Andreas Fái-Pozsár, Regensburg, Germany toom, Centor-Warenhandels-GmbH, Köln, Germany Kiepenkerl, Bruno Nebelung GmbH, Everswinkel, Germany Thompson & Morgan (UK) Ltd., Ipswich, UK Reinhard Lühring, private farmer, Rhauderfehn, Germany Botanical Garden of the University of Oldenburg, Oldenburg, Germany Exotic-Samen Wolfgang Meier, Fürstenwalde, Germany

31 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 32 of 37

Table 3. Parameters For Calibration Curves: Slope, Axis Intercept, Pearson Correlation Coefficient R2, As Well As Relative Standard Deviation Of Analytical Repeat Injections (RSD) And Signal-To-Noise Ratio (S/N) Of The Five Glucosinolate Standards. glucosinolate

slope

axis intercept

gluconapin progoitrin glucoraphanin gluconasturtiin glucobrassicin

0.0012 0.0006 0.0039 0.0005 0.0012

-0.1089 -0.3175 0.0663 -0.0386 0.1321

coefficient of determination R2 0.9998 0.9980 0.9998 0.9999 1.0

RSD (%)

S/N

3.97 3.01 4.16 4.37 3.23

45 230 45 195 55

32 ACS Paragon Plus Environment

Page 33 of 37

Journal of Agricultural and Food Chemistry

Figure graphics

Figure 1:

33 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 34 of 37

Figure 2:

34 ACS Paragon Plus Environment

Page 35 of 37

Journal of Agricultural and Food Chemistry

Figure 3:

35 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 36 of 37

Figure 4:

36 ACS Paragon Plus Environment

Page 37 of 37

Journal of Agricultural and Food Chemistry

Graphic for table of contents

37 ACS Paragon Plus Environment