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.

Tackling drought stress: RECEPTOR-LIKE KINASES present new approaches

Global climate change and a growing population require tackling the reduction in arable land and improving biomass production and seed yield per area under varying conditions. One of these conditions is suboptimal water availability. Here, we review some of the classical approaches to dealing with plant response to drought stress and we evaluate how research on RECEPTOR-LIKE KINASES (RLKs) can contribute to improving plant performance under drought stress. RLKs are considered as key regulators of plant architecture and growth behavior, but they also function in defense and stress responses. The available literature and analyses of available transcript profiling data indeed suggest that RLKs can play an important role in optimizing plant responses to drought stress. In addition, RLK pathways are ideal targets for nontransgenic approaches, such as synthetic molecules, providing a novel strategy to manipulate their activity and supporting translational studies from model species, such as Arabidopsis thaliana, to economically useful crops.

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.

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.

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 < 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.

Global scale transcriptome analysis of Arabidopsis embryogenesis in vitro

Background Somatic embryogenesis (SE) in plants is a process by which embryos are generated directly from somatic cells, rather than from the fused products of male and female gametes. Despite the detailed expression analysis of several somatic-to-embryonic marker genes, a comprehensive understanding of SE at a molecular level is still lacking. The present study was designed to generate high resolution transcriptome datasets for early SE providing the way for future research to understand the underlying molecular mechanisms that regulate this process. We sequenced Arabidopsis thaliana somatic embryos collected from three distinct developmental time-points (5, 10 and 15 d after in vitro culture) using the Illumina HiSeq 2000 platform. Results This study yielded a total of 426,001,826 sequence reads mapped to 26,520 genes in the A. thaliana reference genome. Analysis of embryonic cultures after 5 and 10 d showed differential expression of 1,195 genes; these included 778 genes that were more highly expressed after 5 d as compared to 10 d. Moreover, 1,718 genes were differentially expressed in embryonic cultures between 10 and 15 d. Our data also showed at least eight different expression patterns during early SE; the majority of genes are transcriptionally more active in embryos after 5 d. Comparison of transcriptomes derived from somatic embryos and leaf tissues revealed that at least 4,951 genes are transcriptionally more active in embryos than in the leaf; increased expression of genes involved in DNA cytosine methylation and histone deacetylation were noted in embryogenic tissues. In silico expression analysis based on microarray data found that approximately 5% of these genes are transcriptionally more active in somatic embryos than in actively dividing callus and non-dividing leaf tissues. Moreover, this identified 49 genes expressed at a higher level in somatic embryos than in other tissues. This included several genes with unknown function, as well as others related to oxidative and osmotic stress, and auxin signalling. Conclusions The transcriptome information provided here will form the foundation for future research on genetic and epigenetic control of plant embryogenesis at a molecular level. In follow-up studies, these data could be used to construct a regulatory network for SE; the genes more highly expressed in somatic embryos than in vegetative tissues can be considered as potential candidates to validate these networks.

Improving power and accuracy of genome-wide association studies via a multi-locus mixed linear model methodology

Genome-wide association studies (GWAS) have been widely used in genetic dissection of complex traits. However, common methods are all based on a fixed-SNP-effect mixed linear model (MLM) and single marker analysis, such as efficient mixed model analysis (EMMA). These methods require Bonferroni correction for multiple tests, which often is too conservative when the number of markers is extremely large. To address this concern, we proposed a random-SNP-effect MLM (RMLM) and a multi-locus RMLM (MRMLM) for GWAS. The RMLM simply treats the SNP-effect as random, but it allows a modified Bonferroni correction to be used to calculate the threshold p value for significance tests. The MRMLM is a multi-locus model including markers selected from the RMLM method with a less stringent selection criterion. Due to the multi-locus nature, no multiple test correction is needed. Simulation studies show that the MRMLM is more powerful in QTN detection and more accurate in QTN effect estimation than the RMLM, which in turn is more powerful and accurate than the EMMA. To demonstrate the new methods, we analyzed six flowering time related traits in Arabidopsis thaliana and detected more genes than previous reported using the EMMA. Therefore, the MRMLM provides an alternative for multi-locus GWAS.

‘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.