Glucosinolates in brassica vegetables: the influence of the food supply chain on intake, bioavailability and human health

Glucosinolates (GLSs) are found in Brassica vegetables. Examples of these sources include cabbage, Brussels sprouts, broccoli, cauliflower and various root vegetables (e.g. radish and turnip). A number of epidemiological studies have identified an inverse association between consumption of these vegetables and the risk of colon and rectal cancer. Animal studies have shown changes in enzyme activities and DNA damage resulting from consumption of Brassica vegetables or isothiocyanates, the breakdown products (BDP) of GLSs in the body. Mechanistic studies have begun to identify the ways in which the compounds may exert their protective action but the relevance of these studies to protective effects in the human alimentary tract is as yet unproven. In vitro studies with a number of specific isothiocyanates have suggested mechanisms that might be the basis of their chemoprotective effects. The concentration and composition of the GLSs in different plants, but also within a plant (e.g. in the seeds, roots or leaves), can vary greatly and also changes during plant development. Furthermore, the effects of various factors in the supply chain of Brassica vegetables including breeding, cultivation, storage and processing on intake and bioavailability of GLSs are extensively discussed in this paper.

The potential to intensify sulforaphane formation in cooked broccoli (Brassica oleracea var. italica) using mustard seeds (Sinapis alba)

Sulforaphane, a naturally occurring cancer chemopreventive, is the hydrolysis product of glucoraphanin, the main glucosinolate in broccoli. The hydrolysis requires myrosinase isoenzyme to be present in sufficient activity; however processing leads to its denaturation and hence reduced hydrolysis. In this study, the effect of adding mustard seeds, which has a more resilient isoform of myrosinase, to processed broccoli was investigated with a view to intensify the formation of sulforaphane. Thermal inactivation of myrosinase from both broccoli and mustard seeds was studied. Thermal degradation of broccoli glucoraphanin was investigated in addition to the effects of thermal processing on the formation of sulforaphane and sulforaphane nitrile. Limited thermal degradation of glucoraphanin (less than 12 %) was observed when broccoli was placed in vacuum sealed bag (sous vide) and cooked in a water bath at 100 ºC for 8 and 12 min. Boiling broccoli in water prevented the formation of any significant levels of sulforaphane due to inactivated myrosinase. However, addition of powdered mustard seeds to the heat processed broccoli significantly increased the formation of sulforaphane.

Shoot calcium and magnesium concentrations differ between subtaxa, are highly heritable, and associate with potentially pleiotropic loci in Brassica oleracea

Calcium (Ca) and magnesium (Mg) are the most abundant group II elements in both plants and animals. Genetic variation in shoot Ca and shoot Mg concentration (shoot Ca and Mg) in plants can be exploited to biofortify food crops and thereby increase dietary Ca and Mg intake for humans and livestock. We present a comprehensive analysis of within-species genetic variation for shoot Ca and Mg, demonstrating that shoot mineral concentration differs significantly between subtaxa (varietas). We established a structured diversity foundation set of 376 accessions to capture a high proportion of species-wide allelic diversity within domesticated Brassica oleracea, including representation of wild relatives (C genome, 1n = 9) from natural populations. These accessions and 74 modern F-1 hybrid cultivars were grown in glasshouse and field environments. Shoot Ca and Mg varied 2- and 2.3-fold, respectively, and was typically not inversely correlated with shoot biomass, within most subtaxa. The closely related capitata (cabbage) and sabauda (Savoy cabbage) subtaxa consistently had the highest mean shoot Ca and Mg. Shoot Ca and Mg in glasshouse-grown plants was highly correlated with data from the field. To understand and dissect the genetic basis of variation in shoot Ca and Mg, we studied homozygous lines from a segregating B. oleracea mapping population. Shoot Ca and Mg was highly heritable (up to 40). Quantitative trait loci (QTL) for shoot Ca and Mg were detected on chromosomes C2, C6, C7, C8, and, in particular, C9, where QTL accounted for 14 to 55 of the total genetic variance. The presence of QTL on C9 was substantiated by scoring recurrent backcross substitution lines, derived from the same parents. This also greatly increased the map resolution, with strong evidence that a 4-cM region on C9 influences shoot Ca. This region corresponds to a 0.41-Mb region on Arabidopsis (Arabidopsis thaliana) chromosome 5 that includes 106 genes. There is also evidence that pleiotropic loci on C8 and C9 affect shoot Ca and Mg. Map-based cloning of these loci will reveal how shoot-level phenotypes relate to Ca 21 and Mg 21 uptake and homeostasis at the molecular level.

Shoot zinc (Zn) concentration varies widely within Brassica oleracea L. and is affected by soil Zn and phosphorus (P) levels

The low availability of zinc (Zn) in soils and crops affects dietary Zn intake worldwide. This study sought to determine if the natural genetic variation in shoot Zn concentrations (Zn(shoot)) is sufficient to pursue a crop improvement breeding strategy in a leafy vegetable crop. The gene-pool of Brassica oleracea L. was sampled using a large (n = 376) diversity foundation set (DFS), representing almost all species-wide common allelic variation, and 74 commercial varieties (mostly F(1)). The DFS genotypes were grown at low and high soil phosphorus (P) levels under glasshouse and field conditions, and also in a Zn-deficient soil, with or without Zn-fertilisation, in a glasshouse. Despite the large variation in Zn(shoot) among genotypes, environment had a profound effect on Zn(shoot) The heritability of Zn(shoot) was significant, but relatively low, among 90 doubled-haploid (DH) lines from a mapping population. While several quantitative trait loci (QTL) associated with Zn(shoot) occurred on chromosomes C2, C3, C5, C7, and C9, these were generally weak and conditional upon growth conditions. Breeding for Zn(shoot) in B. oleracea is therefore likely to be challenging. Shoot P concentrations increased substantially in all genotypes under low soil Zn conditions. Conversely, only some genotypes had increased Zn(shoot) at low soil P levels. Sufficient natural genetic variation may therefore exist to study some of the interactions between Zn and P nutrition.

Shoot yield drives phosphorus use efficiency in Brassica oleracea and correlates with root architecture traits

The environmental and financial costs of using inorganic phosphate fertilizers to maintain crop yield and quality are high. Breeding crops that acquire and use phosphorus (P) more efficiently could reduce these costs. The variation in shoot P concentration (shoot-P) and various measures of P use efficiency (PUE) were quantified among 355 Brassica oleracea L. accessions, 74 current commercial cultivars, and 90 doubled haploid (DH) mapping lines from a reference genetic mapping population. Accessions were grown at two or more external P concentrations in glasshouse experiments; commercial and DH accessions were also grown in replicated field experiments. Within the substantial species-wide diversity observed for shoot-P and various measures of PUE in B. oleracea, current commercial cultivars have greater PUE than would be expected by chance. This may be a consequence of breeding for increased yield, which is a significant component of most measures of PUE, or early establishment. Root development and architecture correlate with PUE; in particular, lateral root number, length, and growth rate. Significant quantitative trait loci associated with shoot-P and PUE occur on chromosomes C3 and C7. These data provide information to initiate breeding programmes to improve PUE in B. oleracea.

Regulatory hotspots are associated with plant gene expression under varying soil phosphorus supply in Brassica rapa

Gene expression is a quantitative trait that can be mapped genetically in structured populations to identify expression quantitative trait loci (eQTL). Genes and regulatory networks underlying complex traits can subsequently be inferred. Using a recently released genome sequence, we have defined cis- and trans-eQTL and their environmental response to low phosphorus (P) availability within a complex plant genome and found hotspots of trans-eQTL within the genome. Interval mapping, using P supply as a covariate, revealed 18,876 eQTL. trans-eQTL hotspots occurred on chromosomes A06 and A01 within Brassica rapa; these were enriched with P metabolism-related Gene Ontology terms (A06) as well as chloroplast-and photosynthesis-related terms (A01). We have also attributed heritability components to measures of gene expression across environments, allowing the identification of novel gene expression markers and gene expression changes associated with low P availability. Informative gene expression markers were used to map eQTL and P use efficiency-related QTL. Genes responsive to P supply had large environmental and heritable variance components. Regulatory loci and genes associated with P use efficiency identified through eQTL analysis are potential targets for further characterization and may have potential for crop improvement.

A Brassica exon array for whole-transcript gene expression profiling

Affymetrix GeneChip (R) arrays are used widely to study transcriptional changes in response to developmental and environmental stimuli. GeneChip (R) arrays comprise multiple 25-mer oligonucleotide probes per gene and retain certain advantages over direct sequencing. For plants, there are several public GeneChip (R) arrays whose probes are localised primarily in 39 exons. Plant whole-transcript (WT) GeneChip (R) arrays are not yet publicly available, although WT resolution is needed to study complex crop genomes such as Brassica, which are typified by segmental duplications containing paralogous genes and/or allopolyploidy. Available sequence data were sampled from the Brassica A and C genomes, and 142,997 gene models identified. The assembled gene models were then used to establish a comprehensive public WT exon array for transcriptomics studies. The Affymetrix GeneChip (R) Brassica Exon 1.0 ST Array is a 5 mu M feature size array, containing 2.4 million 25-base oligonucleotide probes representing 135,201 gene models, with 15 probes per gene distributed among exons. Discrimination of the gene models was based on an E-value cut-off of 1E(-5), with <= 98 sequence identity. The 135 k Brassica Exon Array was validated by quantifying transcriptome differences between leaf and root tissue from a reference Brassica rapa line (R-o-18), and categorisation by Gene Ontologies (GO) based on gene orthology with Arabidopsis thaliana. Technical validation involved comparison of the exon array with a 60-mer array platform using the same starting RNA samples. The 135 k Brassica Exon Array is a robust platform. All data relating to the array design and probe identities are available in the public domain and are curated within the BrassEnsembl genome viewer at http://www.brassica.info/BrassEnsembl/index.html.

Distribution of calcium (Ca) and magnesium (Mg) in the leaves of Brassica rapa under varying exogenous Ca and Mg supply

Background and Aims Leafy vegetable Brassica crops are an important source of dietary calcium (Ca) and magnesium (Mg) and represent potential targets for increasing leaf Ca and Mg concentrations through agronomy or breeding. Although the internal distribution of Ca and Mg within leaves affects the accumulation of these elements, such data are not available for Brassica. The aim of this study was to characterize the internal distribution of Ca and Mg in the leaves of a vegetable Brassica and to determine the effects of altered exogenous Ca and Mg supply on this distribution. Methods Brassica rapa ssp. trilocularis ‘R-o-18’ was grown at four different Ca:Mg treatments for 21 d in a controlled environment. Concentrations of Ca and Mg were determined in fully expanded leaves using inductively coupled plasma-mass spectrometry (ICP-MS). Internal distributions of Ca and Mg were determined in transverse leaf sections at the base and apex of leaves using energy-dispersive X-ray spectroscopy (EDS) with cryo-scanning electron microscopy (cryo-SEM). Key Results Leaf Ca and Mg concentrations were greatest in palisade and spongy mesophyll cells, respectively, although this was dependent on exogenous supply. Calcium accumulation in palisade mesophyll cells was enhanced slightly under high Mg supply; in contrast, Mg accumulation in spongy mesophyll cells was not affected by Ca supply. Conclusions The results are consistent with Arabidopsis thaliana and other Brassicaceae, providing phenotypic evidence that conserved mechanisms regulate leaf Ca and Mg distribution at a cellular scale. The future study of Arabidopsis gene orthologues in mutants of this reference B. rapa genotype will improve our understanding of Ca and Mg homeostasis in plants and may provide a model-to-crop translation pathway for targeted breeding.

Genetic analysis of potassium use efficiency in Brassica oleracea

Potassium (K) fertilizers are used in intensive and extensive agricultural systems to maximize production. However, there are both financial and environmental costs to K-fertilization. It is therefore important to optimize the efficiency with which K-fertilizers are used. Cultivating crops that acquire and/or utilize K more effectively can reduce the use of K-fertilizers. The aim of the present study was to determine the genetic factors affecting K utilization efficiency (KUtE), defined as the reciprocal of shoot K concentration (1/K(shoot)), and K acquisition efficiency (KUpE), defined as shoot K content, in Brassica oleracea. Genetic variation in K(shoot) was estimated using a structured diversity foundation set (DFS) of 376 accessions and in 74 commercial genotypes grown in glasshouse and field experiments that included phosphorus (P) supply as a treatment factor. Chromosomal quantitative trait loci (QTL) associated with K(shoot) and KUpE were identified using a genetic mapping population grown in the glasshouse and field. Putative QTL were tested using recurrent backcross substitution lines in the glasshouse. More than two-fold variation in K(shoot) was observed among DFS accessions grown in the glasshouse, a significant proportion of which could be attributed to genetic factors. Several QTL associated with K(shoot) were identified, which, despite a significant correlation in K(shoot) among genotypes grown in the glasshouse and field, differed between these two environments. A QTL associated with K(shoot) in glasshouse-grown plants (chromosome C7 at 62 center dot 2 cM) was confirmed using substitution lines. This QTL corresponds to a segment of arabidopsis chromosome 4 containing genes encoding the K(+) transporters AtKUP9, AtAKT2, AtKAT2 and AtTPK3. There is sufficient genetic variation in B. oleracea to breed for both KUtE and KUpE. However, as QTL associated with these traits differ between glasshouse and field environments, marker-assisted breeding programmes must consider carefully the conditions under which the crop will be grown.

Effects of planting density and nitrogen application on seed yield and other morphological traits of the leafy vegetable kale (Brassica oleracea).

The relationship between food security and sustainable land use is considered to be of the uttermost importance to increase yields without having to increase the agricultural land area over which crops are grown. In the present study nitrogen concentration (25 and 85 kg ha-1) and planting density (6.7, 10 and 25 plants m-2) were investigated for their effect on whole plant physiology and pod seed yield in kale (Brassica oleracea), to determine if the fruit (pod) yield could be manipulated agronomically. Nitrogen concentration did not significantly affect seed yield and it is therefore recommended that the lower concentration be used commercially. Conversely planting density did have a significant effect with increases in seed yield observed at the highest planting density of 25 plants m-2, therefore this high planting density would be recommended commercially to maximise area efficiency, highlighting that simple agronomic changes are capable of increasing crop yields over a set area.

Genetical and comparative genomics of Brassica under altered Ca supply identifies Arabidopsis Ca-transporter orthologs

Although Ca transport in plants is highly complex, the overexpression of vacuolar Ca2+ transporters in crops is a promising new technology to improve dietary Ca supplies through biofortification. Here, we sought to identify novel targets for increasing plant Ca accumulation using genetical and comparative genomics. Expression quantitative trait locus (eQTL) mapping to 1895 cis- and 8015 trans-loci were identified in shoots of an inbred mapping population of Brassica rapa (IMB211 × R500); 23 cis- and 948 trans-eQTLs responded specifically to altered Ca supply. eQTLs were screened for functional significance using a large database of shoot Ca concentration phenotypes of Arabidopsis thaliana. From 31 Arabidopsis gene identifiers tagged to robust shoot Ca concentration phenotypes, 21 mapped to 27 B. rapa eQTLs, including orthologs of the Ca2+ transporters At-CAX1 and At-ACA8. Two of three independent missense mutants of BraA.cax1a, isolated previously by targeting induced local lesions in genomes, have allele-specific shoot Ca concentration phenotypes compared with their segregating wild types. BraA.CAX1a is a promising target for altering the Ca composition of Brassica, consistent with prior knowledge from Arabidopsis. We conclude that multiple-environment eQTL analysis of complex crop genomes combined with comparative genomics is a powerful technique for novel gene identification/prioritization.

Foraging strategies of natural enemies of the peach-potato aphid Myzus persicae

The foraging strategies of two natural enemies of the peach-potato aphid Myzus persicae: the seven-spot ladybird Coccinella septempunctata and the parasitoid wasp Diaeretiella rapae, were investigated. Specifically the roles of plant semiochemicals in the location of plants infested with M. persicae by these natural enemies were examined. I investigated the olfactory responses of female C. septempunctata to volatiles collected from M. persicae and four Brassica cultivars; Brassica rapa, B. juncea, B. napus cultivar ‘Apex’ and B. napus cultivar ‘Courage’ and wild-type Arabidopsis thaliana that were: undamaged, previously infested by M. persicae and infested with M. persicae. C. septempunctata showed no attraction to volatiles from M. persicae alone. C. septempunctata significantly changed its searching behaviour in response to plant volatiles from B. rapa, B. napus cv. ‘Apex’ and Arabidopsis infested with M. persicae. C. septempunctata was also found to display a significant turning bias when foraging on a branching horizontal wire stem. A model was developed to investigate how turning biases affect the foraging efficiency of C. septempunctata in dichotomous branched environments. Simulations using this model indicated that turning biases could potentially increase searching efficiency. D. rapae showed a significant preference for volatiles from M. persicae infested wild-type Arabidopsis but no preference to volatiles from M. persicae alone or M. persicae honeydew. Volatile emissions by Arabidopsis were shown to be localised to the area of aphid-infestation rather than systemic. Using gas chromatography plants infested with M. persicae were shown to emit a quantitatively different volatile blend than undamaged plants. In experiments with jasmonate mutants of Arabidopsis the jasmonate (octadecanoid) wound response pathway was implicated as being important for the production of M. persicae induced volatiles, attractive to D. rapae. Other wound response pathways were also found to be involved in the production of the full blend of M. persicae induced volatiles.

Premature germination of resting spores as a means of Protecting brassica crops from Plasmodiophora brassicae Wor., (clubroot)

Premature germination of resting spores as a means of protecting brassica crops from Plasmodiphora brassicae Wor., (Clubroot). Crop Protection. Clubroot disease causes substantial yield and quality losses in broadacre oil seed and intensive vegetable brassica crops worldwide. The causal microbe Plasmodiophora brassicae Wor., perennates as soil-borne dormant resting spores. Their germination is triggered by exudates from host roots. A valuable addition to sustainable integrated control strategies could be developed by identifying and synthesising the molecules responsible for stimulating resting spore germination. This paper reports experiments in which stimulatory exudates were collected from brassica roots following exposure to infective stages of P. brassicae. Analyses identified a germination signalling molecule of circa 1 kDa formed of glucose sub-units. Mass spectral analyses showed this to be a complex hexasaccharide carbohydrate with structural similarities to the components of plant cell walls. This is the first report of a host generated hexasaccharide which is capable of stimulating the germination of resting spores of P. brassicae. The implications for environmentally benign control of clubroot are discussed briefly.

Thermal and pressure stability of myrosinase enzymes from black mustard (Brassica nigra L. W.D.J Koch. var. nigra), brown mustard (Brassica juncea L. Czern. var. juncea) and yellow mustard (Sinapsis alba L. Subsp Maire) seeds

This study investigates the effects of temperature and pressure on inactivation of myrosinase extracted from black, brown and yellow mustard seeds. Brown mustard had higher myrosinase activity (2.75 un/mL) than black (1.50 un/mL) and yellow mustard (0.63 un/mL). The extent of enzyme inactivation increased with pressure (600-800 MPa) and temperature (30-70 °C) for all the mustard seeds. However, at combinations of lower pressures (200-400 MPa) and high temperatures (60-80 °C), there was less inactivation. For example, application of 300 MPa and 70 °C for 10 minutes retained 20%, 80% and 65% activity in yellow, black and brown mustard, respectively, whereas the corresponding activity retentions when applying only heat (70 °C, 10min) were 0%, 59% and 35%. Thus, application of moderate pressures (200-400 MPa) can potentially be used to retain myrosinase activity needed for subsequent glucosinolate hydrolysis.

Exploring the effects of genotypical and phenotypical variations in bitter taste sensitivity on perception, liking and intake of brassica vegetables in the UK

Brassicaceous vegetables (BV) have chemoprotective effects and yet consumption of BV in the UK is low. Previous studies suggest perception, liking and intake of BV are influenced by bitter taste sensitivity which this study further explores. Phenotypical taste sensitivity of 136 subjects was classified using propythiouracil (PROP) and sodium chloride and fungiform papillae density (FPD) was measured from tongue images. Polymorphisms of TAS2R38 and gustin (CA6) genes were analysed. Liking and bitterness of four raw vegetables (two BV (broccoli and white cabbage) and two non-BV (spinach and courgette)), as well as habitual consumption, were evaluated. There was a significant association between TAS2R38 genotype and PROP taster status (p<0.0001) and between FPD and PROP taster status (p=0.029). Individuals with greater sensitivity for PROP predominantly had TAS2R38 PAV/PAV genotype and greater FPD. BV were perceived as more bitter than non-BV (p<0.0001) with PAV/PAV subjects perceiving significantly stronger bitter intensity. There was a significant difference in liking for the four vegetables (p=0.002), and between consumers of different TAS2R38 genotype (p=0.0024). Individuals with TAS2R38 AVI/AVI genotype liked BV more. Regarding intake, both PAV/PAV and AVI/AVI individuals consumed more total vegetables and BV than PAV/AVI. Although PROP nontasters tended to consume more vegetables and BV than the other two phenotype groups, liking and vegetable intake were not significantly affected by taste phenotype. Although there was not a significant effect of CA6 genotype on bitterness ratings, there was a significant interaction between CA6 and TAS2R38, and in addition CA6 genotype was significantly associated with BV intake. However, these effects require validation as the proportions of the population with the CA6 G/G genotype was extremely small (7%). Our results confirmed that bitter taste perception in vegetables was influenced by both genotype and phenotype of bitter taste sensitivity. Moreover, our findings demonstrated that neither genotype nor phenotype of taste sensitivity alone accurately predict vegetable liking and intake as demographic factors were found to have a substantial influence.