Using genomic DNA-based probe-selection to improve the sensitivity of high-density oligonucleotide arrays when applied to heterologous species

High-density oligonucleotide (oligo) arrays are a powerful tool for transcript profiling. Arrays based on GeneChip® technology are amongst the most widely used, although GeneChip® arrays are currently available for only a small number of plant and animal species. Thus, we have developed a method to improve the sensitivity of high-density oligonucleotide arrays when applied to heterologous species and tested the method by analysing the transcriptome of Brassica oleracea L., a species for which no GeneChip® array is available, using a GeneChip® array designed for Arabidopsis thaliana (L.) Heynh. Genomic DNA from B. oleracea was labelled and hybridised to the ATH1-121501 GeneChip® array. Arabidopsis thaliana probe-pairs that hybridised to the B. oleracea genomic DNA on the basis of the perfect-match (PM) probe signal were then selected for subsequent B. oleracea transcriptome analysis using a .cel file parser script to generate probe mask files. The transcriptional response of B. oleracea to a mineral nutrient (phosphorus; P) stress was quantified using probe mask files generated for a wide range of gDNA hybridisation intensity thresholds. An example probe mask file generated with a gDNA hybridisation intensity threshold of 400 removed > 68 % of the available PM probes from the analysis but retained >96 % of available A. thaliana probe-sets. Ninety-nine of these genes were then identified as significantly regulated under P stress in B. oleracea, including the homologues of P stress responsive genes in A. thaliana. Increasing the gDNA hybridisation intensity thresholds up to 500 for probe-selection increased the sensitivity of the GeneChip® array to detect regulation of gene expression in B. oleracea under P stress by up to 13-fold. Our open-source software to create probe mask files is freely available http://affymetrix.arabidopsis.info/xspecies/ webcite and may be used to facilitate transcriptomic analyses of a wide range of plant and animal species in the absence of custom arrays.

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.

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.

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.

AtTRB1, a telomeric DNA-binding protein from Arabidopsis, is concentrated in the nucleolus and shows highly dynamic association with chromatin

AtTRB1, 2 and 3 are members of the SMH (single Myb histone) protein family, which comprises double-stranded DNA-binding proteins that are specific to higher plants. They are structurally conserved, containing a Myb domain at the N-terminus, a central H1/H5-like domain and a C-terminally located coiled-coil domain. AtTRB1, 2 and 3 interact through their Myb domain specifically with telomeric double-stranded DNA in vitro, while the central H1/H5-like domain interacts non-specifically with DNA sequences and mediates protein–protein interactions. Here we show that AtTRB1, 2 and 3 preferentially localize to the nucleus and nucleolus during interphase. Both the central H1/H5-like domain and the Myb domain from AtTRB1 can direct a GFP fusion protein to the nucleus and nucleolus. AtTRB1–GFP localization is cell cycle-regulated, as the level of nuclear-associated GFP diminishes during mitotic entry and GFP progressively re-associates with chromatin during anaphase/telophase. Using fluorescence recovery after photobleaching and fluorescence loss in photobleaching, we determined the dynamics of AtTRB1 interactions in vivo. The results reveal that AtTRB1 interaction with chromatin is regulated at two levels at least, one of which is coupled with cell-cycle progression, with the other involving rapid exchange.

Plant U13 orthologues and orphan snoRNAs identified by RNomics of RNA from Arabidopsis nucleoli

Small nucleolar RNAs (snoRNAs) and small Cajal body-specific RNAs (scaRNAs) are non-coding RNAs whose main function in eukaryotes is to guide the modification of nucleotides in ribosomal and spliceosomal small nuclear RNAs, respectively. Full-length sequences of Arabidopsis snoRNAs and scaRNAs have been obtained from cDNA libraries of capped and uncapped small RNAs using RNA from isolated nucleoli from Arabidopsis cell cultures. We have identified 31 novel snoRNA genes (9 box C/D and 22 box H/ACA) and 15 new variants of previously described snoRNAs. Three related capped snoRNAs with a distinct gene organization and structure were identified as orthologues of animal U13snoRNAs. In addition, eight of the novel genes had no complementarity to rRNAs or snRNAs and are therefore putative orphan snoRNAs potentially reflecting wider functions for these RNAs. The nucleolar localization of a number of the snoRNAs and the localization to nuclear bodies of two putative scaRNAs was confirmed by in situ hybridization. The majority of the novel snoRNA genes were found in new gene clusters or as part of previously described clusters. These results expand the repertoire of Arabidopsis snoRNAs to 188 snoRNA genes with 294 gene variants.

Glucosinolate-accumulating S-cells in Arabidopsis leaves and flower stalks undergo programmed cell death at early stages of differentiation

Plant secondary metabolites glucosinolates (GSL) have important functions in plant resistance to herbivores and pathogens. We identified all major GSL that are accumulated in S-cells in Arabidopsis by MALDI-TOF MS, and estimated by LC-MS that the total GSL concentration in these cells is above 130 mM. The precise locations of the S-cells outside phloem bundles in rosette and cauline leaves and in flower stalks were visualised using sulphur mapping by cryo-SEM/EDX. S-cells contain up to 40% of total sulphur in flower stalk tissues. S-cells in emerging flower stalks and developing leaf tissues show typical signs of Programmed Cell Death (PCD) or apoptosis, such as chromatin condensation in the nucleus and blebbing of the membranes. TUNEL staining for DNA double strand breaks confirmed PCD in S-cells in postmeristematic tissues in the flower stalk as well as in the leaf. Our results show that S-cells in postmeristematic tissues proceed to an extreme degree of metabolic specialisation besides PCD. Accumulation and maintenance of a high concentration of GSL in these cells are accompanied by degradation of a number of cell organelles. The substantial changes in the cell composition during S-cell differentiation indicate the importance of this particular GSL-based phloem defence system. The specific anatomy of the S-cells and ability to accumulate specialised secondary metabolites is similar to that of the non-articulated laticifer cells in latex plants and thus indicates a common evolutionary origin.

Changes in gene expression in Arabidopsis shoots during phosphate starvation and the potential for developing smart plants

Our aim was to generate and prove the concept of "smart" plants to monitor plant phosphorus (P) status in Arabidopsis. Smart plants can be genetically engineered by transformation with a construct containing the promoter of a gene up-regulated specifically by P starvation in an accessible tissue upstream of a marker gene such as beta-glucuronidase (GUS). First, using microarrays, we identified genes whose expression changed more than 2.5-fold in shoots of plants growing hydroponically when P, but not N or K, was withheld from the nutrient solution. The transient changes in gene expression occurring immediately (4 h) after P withdrawal were highly variable, and many nonspecific, shock-induced genes were up-regulated during this period. However, two common putative cis-regulatory elements (a PHO-like element and a TATA box-like element) were present significantly more often in the promoters of genes whose expression increased 4 h after the withdrawal of P compared with their general occurrence in the promoters of all genes represented on the microarray. Surprisingly, the expression of only four genes differed between shoots of P-starved and -replete plants 28 h after P was withdrawn. This lull in differential gene expression preceded the differential expression of a new group of 61 genes 100 h after withdrawing P. A literature survey indicated that the expression of many of these "late" genes responded specifically to P starvation. Shoots had reduced P after 100 h, but growth was unaffected. The expression of SQD1, a gene involved in the synthesis of sulfolipids, responded specifically to P starvation and was increased 100 h after withdrawing P. Leaves of Arabidopsis bearing a SQD1::GUS construct showed increased GUS activity after P withdrawal, which was detectable before P starvation limited growth. Hence, smart plants can monitor plant P status. Transferring this technology to crops would allow precision management of P fertilization, thereby maintaining yields while reducing costs, conserving natural resources, and preventing pollution.

Overexpression of coconut AINTEGUMENTA-like gene, CnANT, promotes in vitro regeneration in transgenic Arabidopsis

Knowledge of the molecular biological changes underlying the process of embryogenesis is important for the improvement of somatic embryogenesis of coconut. Among the transcription factors that control the transition from vegetative to embryogenic growth, members of APETALA2/Ethylene-responsive element binding protein domain family play an important role in promoting embryo development. Significant insights into the role of AP2 genes have been obtained by the ectopic expression of AP2 sub family genes in transgenic Arabidopsis. A homolog of the AINTEGUMENTA-like gene that encodes the two AP2 domains and the linker region was identified in the coconut genome. Phylogenetic analysis showed that this gene, CnANT, encodes a protein that branched with BABY BOOM/PLETHORA clade in the AINTEGUMENTA-like major clade and was similar to the oil palm EgAP2-1 protein. According to real time RT-PCR results, higher expression of CnANT was observed in more mature zygotic embryos. Also, high CnANT expression was recorded in embryogenic callus compared to other stages of somatic embryogenesis. We examined the effect of ectopic CnANT expression on the development and regenerative capacity of transgenic Arabidopsis. Overexpression of CnANT in Arabidopsis induced hormone free regeneration of explants. Furthermore, ectopic expression of CnANT enhanced regeneration in vitro and suggested a role for this gene in cell proliferation during in vitro culture.