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.

Optimising the analysis of transcript data using high density oligonucleotide arrays and genomic DNA-based probe selection

Background: Affymetrix GeneChip arrays are widely used for transcriptomic studies in a diverse range of species. Each gene is represented on a GeneChip array by a probe- set, consisting of up to 16 probe-pairs. Signal intensities across probe- pairs within a probe-set vary in part due to different physical hybridisation characteristics of individual probes with their target labelled transcripts. We have previously developed a technique to study the transcriptomes of heterologous species based on hybridising genomic DNA (gDNA) to a GeneChip array designed for a different species, and subsequently using only those probes with good homology. Results: Here we have investigated the effects of hybridising homologous species gDNA to study the transcriptomes of species for which the arrays have been designed. Genomic DNA from Arabidopsis thaliana and rice (Oryza sativa) were hybridised to the Affymetrix Arabidopsis ATH1 and Rice Genome GeneChip arrays respectively. Probe selection based on gDNA hybridisation intensity increased the number of genes identified as significantly differentially expressed in two published studies of Arabidopsis development, and optimised the analysis of technical replicates obtained from pooled samples of RNA from rice. Conclusion: This mixed physical and bioinformatics approach can be used to optimise estimates of gene expression when using GeneChip arrays.

A comparison of the Thlaspi caerulescens and Thlaspi arvense shoot transcriptomes

Whole-genome transcriptome profiling is revealing how biological systems are regulated at the transcriptional level. This study reports the development of a robust method to profile and compare the transcriptomes of two nonmodel plant species, Thlaspi caerulescens, a zinc (Zn) hyperaccumulator, and Thlaspi arvense, a nonhyperaccumulator, using Affymetrix Arabidopsis thaliana ATH1-121501 GeneChip (R) arrays (Affymetrix, Santa Clara, CA, USA). Transcript abundance was quantified in the shoots of agar- and compost-grown plants of both species. Analyses were optimized using a genomic DNA (gDNA)-based probe-selection strategy based on the hybridization efficiency of Thlaspi gDNA with corresponding A. thaliana probes. In silico alignments of GeneChip (R) probes with Thlaspi gene sequences, and quantitative real-time PCR, confirmed the validity of this approach. Approximately 5000 genes were differentially expressed in the shoots of T. caerulescens compared with T. arvense, including genes involved in Zn transport and compartmentalization. Future functional analyses of genes identified as differentially expressed in the shoots of these closely related species will improve our understanding of the molecular mechanisms of Zn hyperaccumulation.

Tandem quadruplication of HMA4 in the zinc (Zn) and cadmium (Cd) hyperaccumulator Noccaea caerulescens

Zinc (Zn) and cadmium (Cd) hyperaccumulation may have evolved twice in the Brassicaceae, in Arabidopsis halleri and in the Noccaea genus. Tandem gene duplication and deregulated expression of the Zn transporter, HMA4, has previously been linked to Zn/Cd hyperaccumulation in A. halleri. Here, we tested the hypothesis that tandem duplication and deregulation of HMA4 expression also occurs in Noccaea. A Noccaea caerulescens genomic library was generated, containing 36,864 fosmid pCC1FOS (TM) clones with insert sizes similar to 20-40 kbp, and screened with a PCR-generated HMA4 genomic probe. Gene copy number within the genome was estimated through DNA fingerprinting and pooled fosmid pyrosequencing. Gene copy numbers within individual clones was determined by PCR analyses with novel locus specific primers. Entire fosmids were then sequenced individually and reads equivalent to 20-fold coverage were assembled to generate complete whole contigs. Four tandem HMA4 repeats were identified in a contiguous sequence of 101,480 bp based on sequence overlap identities. These were flanked by regions syntenous with up and downstream regions of AtHMA4 in Arabidopsis thaliana. Promoter-reporter beta-glucuronidase (GUS) fusion analysis of a NcHMA4 in A. thaliana revealed deregulated expression in roots and shoots, analogous to AhHMA4 promoters, but distinct from AtHMA4 expression which localised to the root vascular tissue. This remarkable consistency in tandem duplication and deregulated expression of metal transport genes between N. caerulescens and A. halleri, which last shared a common ancestor > 40 mya, provides intriguing evidence that parallel evolutionary pathways may underlie Zn/Cd hyperaccumulation in Brassicaceae.

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.

Transgenerational, dynamic methylation of stomata genes in response to low relative humidity

Transgenerational inheritance of abiotic stress-induced epigenetic modifications in plants has potential adaptive significance and might condition the offspring to improve the response to the same stress, but this is at least partly dependent on the potency, penetrance and persistence of the transmitted epigenetic marks. We examined transgenerational inheritance of low Relative Humidity-induced DNA methylation for two gene loci in the stomatal developmental pathway in Arabidopsis thaliana and the abundance of associated short-interfering RNAs (siRNAs). Heritability of low humidity-induced methylation was more predictable and penetrative at one locus (SPEECHLESS, entropy ≤ 0.02; χ2 &lt; 0.001) than the other (FAMA, entropy ≤ 0.17; χ2 ns). Methylation at SPEECHLESS correlated positively with the continued presence of local siRNAs (r2 = 0.87; p = 0.013) which, however, could be disrupted globally in the progeny under repeated stress. Transgenerational methylation and a parental low humidity-induced stomatal phenotype were heritable, but this was reversed in the progeny under repeated treatment in a previously unsuspected manner.

Investigations into plant biochemical wound-response pathways involved in the production of aphid-induced plant volatiles

Feeding damage to plants by insect herbivores induces the production of plant volatiles, which are attractive to the herbivores natural enemies. Little is understood about the plant biochemical pathways involved in aphid-induced plant volatile production. The aphid parasitoid Diaeretiella rapae can detect and respond to aphid-induced volatiles produced by Arabidopsis thaliana. When given experience of those volatiles, it can learn those cues and can therefore be used as a novel biosensor to detect them. The pathways involved in aphid-induced volatile production were investigated by comparing the responses of D. rapae to volatiles from a number of different transgenic mutants of A. thaliana, mutated in their octadecanoid, ethylene or salicylic acid wound-response pathways and also from wild-type plants. Plants were either undamaged or infested by the peach-potato aphid, Myzus persicae. It is demonstrated that the octadecanoid pathway and specifically the COI1 gene are required for aphid-induced volatile production. The presence of salicylic acid is also involved in volatile production. Using this model system, in combination with A. thaliana plants with single point gene mutations, has potential for the precise dissection of biochemical pathways involved in the production of aphid-induced volatiles

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.

‘Symptomless’ infection by Botrytis cinerea

The study was carried out to clarify the nature of symptomless infection by Botrytis cinerea and to what extent it differs from aggressive necrotic infection in Lactuca sativa (lettuce) and Arabidopsis thaliana. Symptomless plants were produced by dry spore inoculation in plants growing in controlled environmental conditions or in glasshouses. Plating out of surface-disinfected and non-surface-disinfected samples of inoculated, apparently healthy, plants on selective medium revealed that the fungus was spreading from the initial inoculation site to newly developing plant organs both internally and externally. Similar findings were obtained in microscope experiments in which host plants were inoculated with GFP labelled B. cinerea and symptomless spreading was monitored under confocal laser scanning microscope. Spore germination on leaf surface was followed by development of sub-cuticular vesicles and plant cell damage in the infected epidermal cell and a few nearby cells. Sparsely branched long hyphae arose from the vesicles and spread on the leaf surface; spread was mostly on the outer surface of the epidermal layer but occasionally below the cuticle or epidermal cells. In the late symptomless phase, mycelium arising from single vesicles formed several mycelial networks on leaves. Experiments were carried out to compare the extent of gene expression in symptomless and necrotic infections, using RT-qPCR. Expression of selected genes was quantified in tissue samples based on the amount of mRNA of the respective genes found. In both host species, the mRNA concentration of signalling genes bcg1, bmp1 and calcineurin, and the pathogenicity genes bcsod1 and bcpg1 were similar to or slightly greater in symptomless samples than in necrotic samples. The mRNA of the signalling gene bac and pathogenicity genes bcbot1 and bcnep1, were not detected or detected in lower abundance than in necrosis. In lettuce, the leaves developing distant from the site of inoculation showed similar results to A. thaliana, but in healthy leaves close to the site of inoculation mRNA concentrations of bac and bcnep1 were similar to necrotic samples. Thus, in both host species, the fungus grew along with the plant and moved to newly growing plant parts without producing symptoms; during this growth some pathogenicity genes were less expressed than in necrotic infection.

Detection of Oxidation Products of 5-Methyl-2′-Deoxycytidine in Arabidopsis DNA

Epigenetic regulations play important roles in plant development and adaptation to environmental stress. Recent studies from mammalian systems have demonstrated the involvement of ten-eleven translocation (Tet) family of dioxygenases in the generation of a series of oxidized derivatives of 5-methylcytosine (5-mC) in mammalian DNA. In addition, these oxidized 5-mC nucleobases have important roles in epigenetic remodeling and aberrant levels of 5-hydroxymethyl-29-deoxycytidine (5-HmdC) were found to be associated with different types of human cancers. However, there is a lack of evidence supporting the presence of these modified bases in plant DNA. Here we reported the use of a reversed-phase HPLC coupled with tandem mass spectrometry method and stable isotope-labeled standards for assessing the levels of the oxidized 5-mC nucleosides along with two other oxidatively induced DNA modifications in genomic DNA of Arabidopsis. These included 5- HmdC, 5-formyl-29-deoxycytidine (5-FodC), 5-carboxyl-29-deoxycytidine (5-CadC), 5-hydroxymethyl-29-deoxyuridine (5- HmdU), and the (59S) diastereomer of 8,59-cyclo-29-deoxyguanosine (S-cdG). We found that, in Arabidopsis DNA, the levels of 5-HmdC, 5-FodC, and 5-CadC are approximately 0.8 modifications per 106 nucleosides, with the frequency of 5-HmdC (per 5-mdC) being comparable to that of 5-HmdU (per thymidine). The relatively low levels of the 5-mdC oxidation products suggest that they arise likely from reactive oxygen species present in cells, which is in line with the lack of homologous Tetfamily dioxygenase enzymes in Arabidopsis.