Mass spectrometry imaging of glucosinolates in arabidopsis flowers and siliques

Glucosinolates are multi-functional plant secondary metabolites which play a vital role in plant defence and are, as dietary compounds, important to human health and livestock well-being. Knowledge of the tissue-specific regulation of their biosynthesis and accumulation is essential for plant breeding programs. Here, we report that in Arabidopsis thaliana, glucosinolates are accumulated differentially in specific cells of reproductive organs. Using matrix-assisted laser desorption/ionization (MALDI) mass spectrometry imaging (MSI), distribution patterns of three selected compounds, 4-methylsulfinylbutyl (glucoraphanin), indol-3-ylmethyl (glucobrassicin), and 4-benzoyloxybutyl glucosinolates, were mapped in the tissues of whole flower buds, sepals and siliques. The results show that tissue localization patterns of aliphatic glucosinolate glucoraphanin and 4-benzoyloxybutyl glucosinolate were similar, but indole glucosinolate glucobrassicin had different localisation, indicating a possible difference in function. The high resolution images obtained by a complementary approach, cryo-SEM Energy Dispersive X-ray analysis (cryo-SEM-EDX), confirmed increased concentration of sulphur in areas with elevated amounts of glucosinolates, and allowed identifying the cell types implicated in accumulation of glucosinolates. High concentration of sulphur was found in S-cells adjacent to the phloem in pedicels and siliques, indicating the presence of glucosinolates. Moreover, both MALDI MSI and cryo-SEM-EDX analyses indicated accumulation of glucosinolates in cells on the outer surface of the sepals, suggesting that a layer of glucosinolate-accumulating epidermal cells protects the whole of the developing flower, in addition to the S-cells, which protect the phloem. This research demonstrates the high potential of MALDI MSI for understanding the cell-specific compartmentation of plant metabolites and its regulation.

Cesium toxicity in Arabidopsis

Cesium (Cs) is chemically similar to potassium (K). However, although K is an essential element, Cs is toxic to plants. Two contrasting hypotheses to explain Cs toxicity have been proposed: (1) extracellular Cs+ prevents K+ uptake and, thereby, induces K starvation; and (2) intracellular Cs+ interacts with vital K+-binding sites in proteins, either competitively or noncompetitively, impairing their activities. We tested these hypotheses with Arabidopsis (Arabidopsis thaliana). Increasing the Cs concentration in the agar (Cs(agar)) on which Arabidopsis were grown reduced shoot growth. Increasing the K concentration in the agar (K(agad)) increased the Cs(agar) at which Cs toxicity was observed. However, although increasing Cs(agar) reduced shoot K concentration (K(shoot)), the decrease in shoot growth appeared unrelated to K(shoot) per se. Furthermore, the changes in gene expression in Cs-intoxicated plants differed from those of K-starved plants, suggesting that Cs intoxication was not perceived genetically solely as K starvation. In addition to reducing K(shoot) increasing Cs(agar) also increased shoot Cs concentration (Cs(shoot)), but shoot growth appeared unrelated to Cs(shoot) per se. The relationship between shoot growth and Cs(shoot)/Kt(shoot) suggested that, at a nontoxic Cs(shoot) growth was determined by K(shoot) but that the growth of Cs-intoxicated plants was related to the Cs(shoot)/K(shoot) quotient. This is consistent with Cs intoxication resulting from competition between K+ and Cs+ for K+-binding sites on essential proteins.

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.

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.

Zea mays annexins modulate cytosolic free Ca2+ and generate a Ca2+-permeable conductance

Regulation of reactive oxygen species and cytosolic free calcium ([Ca2+](cyt)) is central to plant function. Annexins are small proteins capable of Ca2+-dependent membrane binding or membrane insertion. They possess structural motifs that could support both peroxidase activity and calcium transport. Here, a Zea mays annexin preparation caused increases in [Ca2+] cyt when added to protoplasts of Arabidopsis thaliana roots expressing aequorin. The pharmacological profile was consistent with annexin activation (at the extracellular plasma membrane face) of Arabidopsis Ca2+-permeable nonselective cation channels. Secreted annexins could therefore modulate Ca2+ influx. As maize annexins occur in the cytosol and plasma membrane, they were incorporated at the intracellular face of lipid bilayers designed to mimic the plasma membrane. Here, they generated an instantaneously activating Ca2+-permeable conductance at mildly acidic pH that was sensitive to verapamil and Gd3+ and had a Ca2+-to-K+ permeability ratio of 0.36. These results suggest that cytosolic annexins create a Ca2+ influx pathway directly, particularly during stress responses involving acidosis. A maize annexin preparation also demonstrated in vitro peroxidase activity that appeared independent of heme association. In conclusion, this study has demonstrated that plant annexins create Ca2+-permeable transport pathways, regulate [Ca2+] cyt, and may function as peroxidases in vitro.

Quantitative proteomics using uniform 15N-labeling, mascot and the trans-proteomic pipeline

Stable isotope labeling combined with MS is a powerful method for measuring relative protein abundances, for instance, by differential metabolic labeling of some or all amino acids with 14N and 15N in cell culture or hydroponic media. These and most other types of quantitative proteomics experiments using high-throughput technologies, such as LC-MS/MS, generate large amounts of raw MS data. This data needs to be processed efficiently and automatically, from the mass spectrometer to statistically evaluated protein identifications and abundance ratios. This paper describes in detail an approach to the automated analysis of uniformly 14N/15N-labeled proteins using MASCOT peptide identification in conjunction with the trans-proteomic pipeline (TPP) and a few scripts to integrate the analysis workflow. Two large proteomic datasets from uniformly labeled Arabidopsis thaliana were used to illustrate the analysis pipeline. The pipeline can be fully automated and uses only common or freely available software.

Mechanisms and function of flower and inflorescence reversion

Flower and inflorescence reversion involve a switch from floral development back to vegetative development, thus rendering flowering a phase in an ongoing growth pattern rather than a terminal act of the meristem. Although it can be considered an unusual event, reversion raises questions about the nature and function of flowering. It is linked to environmental conditions and is most often a response to conditions opposite to those that induce flowering. Research on molecular genetic mechanisms underlying plant development over the last 15 years has pinpointed some of the key genes involved in the transition to flowering and flower development. Such investigations have also uncovered mutations which reduce floral maintenance or alter the balance between vegetative and floral features of the plant. How this information contributes to an understanding of floral reversion is assessed here. One issue that arises is whether floral commitment (defined as the ability to continue flowering when inductive conditions no longer exist) is a developmental switch affecting the whole plant or is a mechanism which assigns autonomy to individual meristems. A related question is whether floral or vegetative development is the underlying default pathway of the plant. This review begins by considering how studies of flowering in Arabidopsis thaliana have aided understanding of mechanisms of floral maintenance. Arabidopsis has not been found to revert to leaf production in any of the conditions or genetic backgrounds analysed to date. A clear-cut reversion to leaf production has, however, been described in Impatiens balsamina. It is proposed that a single gene controls whether Impatiens reverts or can maintain flowering when inductive conditions are removed, and it is inferred that this gene functions to control the synthesis or transport of a leaf-generated signal. But it is also argued that the susceptibility of Impatiens to reversion is a consequence of the meristem-based mechanisms controlling development of the flower in this species. Thus, in Impatiens, a leaf-derived signal is critical for completion of flowering and can be considered to be the basis of a plant-wide floral commitment that is achieved without accompanying meristem autonomy. The evidence, derived from in vitro and other studies, that similar mechanisms operate in other species is assessed. It is concluded that most species (including Arabidopsis) are less prone to reversion because signals from the leaf are less ephemeral, and the pathways driving flower development have a high level of redundancy that generates meristem autonomy even when leaf-derived signals are weak. This gives stability to the flowering process, even where its initiation is dependent on environmental cues. On this interpretation, Impatiens reversion appears as an anomaly resulting from an unusual combination of leaf signalling and meristem regulation. Nevertheless, it is shown that the ability to revert can serve a function in the life history strategy (perenniality) or reproductive habit (pseudovivipary) of many plants. In these instances reversion has been assimilated into regular plant development and plays a crucial role there.

The assessment of enriched apoplastic extracts using proteomic approaches

In plant tissues the extracellular environment or apoplast, incorporating the cell wall, is a highly dynamic compartment with a role in many important plant processes including defence, development, signalling and assimilate partitioning. Soluble apoplast proteins from Arabidopsis thaliana, Triticum aestivum and Oryza sativa were separated by two-dimensional electrophoresis. The molecular weights and isoelectric points for the dominant proteins were established prior to excision, sequencing and identification by matrix-assisted laser-desorption ionisation time of flight mass spectrometry (MALDI - TOF MS). From the selected spots, 23 proteins from O. sativa and 25 proteins from A. thaliana were sequenced, of which nine identifications were made in O. sativa (39%) and 14 in A. thaliana (56%). This analysis revealed that: (i) patterns of proteins revealed by two-dimensional electrophoresis were different for each species indicating that speciation could occur at the level of the apoplast, (ii) of the proteins characterised many belonged to diverse families reflecting the multiple functions of the apoplast and (iii), a large number of the apoplast proteins could not be identified indicating that the majority of extracellular proteins are yet to be assigned. The principal proteins identified in the aqueous matrix of the apoplast were involved in defence, i.e. germin-like proteins or glucanases, and cell expansion, i.e. β-D-glucan glucohydrolases. This study has demonstrated that proteomic analysis can be used to resolve the apoplastic protein complement and to identify adaptive changes induced by environmental effectors.

Quantitative proteomics using uniform N-15-labeling, MASCOT, and the trans-proteomic pipeline

Stable isotope labeling combined with MS is a powerful method for measuring relative protein abundances, for instance, by differential metabolic labeling of some or all amino acids with N-14 and N-15 in cell culture or hydroponic media. These and most other types of quantitative proteomics experiments using high-throughput technologies, such as LC-MS/MS, generate large amounts of raw MS data. This data needs to be processed efficiently and automatically, from the mass spectrometer to statistically evaluated protein identifications and abundance ratios. This paper describes in detail an approach to the automated analysis of Uniformly N-14/N-15-labeled proteins using MASCOT peptide identification in conjunction with the trans-proteomic pipeline (TPP) and a few scripts to integrate the analysis workflow. Two large proteomic datasets from uniformly labeled Arabidopsis thaliana were used to illustrate the analysis pipeline. The pipeline can be fully automated and uses only common or freely available software.

Low relative humidity triggers RNA-directed de novo DNA methylation and suppression of genes controlling stomatal development

Environmental cues influence the development of stomata on the leaf epidermis, and allow plants to exert plasticity in leaf stomatal abundance in response to the prevailing growing conditions. It is reported that Arabidopsis thaliana ‘Landsberg erecta’ plants grown under low relative humidity have a reduced stomatal index and that two genes in the stomatal development pathway, SPEECHLESS and FAMA, become de novo cytosine methylated and transcriptionally repressed. These environmentally-induced epigenetic responses were abolished in mutants lacking the capacity for de novo DNA methylation, for the maintenance of CG methylation, and in mutants for the production of short-interfering non-coding RNAs (siRNAs) in the RNA-directed DNA methylation pathway. Induction of methylation was quantitatively related to the induction of local siRNAs under low relative humidity. Our results indicate the involvement of both transcriptional and post-transcriptional gene suppression at these loci in response to environmental stress. Thus, in a physiologically important pathway, a targeted epigenetic response to a specific environmental stress is reported and several of its molecular, mechanistic components are described, providing a tractable platform for future epigenetics experiments. Our findings suggest epigenetic regulation of stomatal development that allows for anatomical and phenotypic plasticity, and may help to explain at least some of the plant’s resilience to fluctuating relative humidity.

Cell walls of developing wheat starchy endosperm: comparison of composition and RNA-Seq transcriptome

The transcriptome of the developing starchy endosperm of hexaploid wheat (Triticum aestivum) was determined using RNA-Seq isolated at five stages during grain fill. This resource represents an excellent way to identify candidate genes responsible for the starchy endosperm cell wall, which is dominated by arabinoxylan (AX), accounting for 70% of the cell wall polysaccharides, with 20% (1,3; 1,4)-beta-D-glucan, 7% glucomannan, and 4% cellulose. A complete inventory of transcripts of 124 glycosyltransferase (GT) and 72 glycosylhydrolase (GH) genes associated with cell walls is presented. The most highly expressed GT transcript (excluding those known to be involved in starch synthesis) was a GT47 family transcript similar to Arabidopsis (Arabidopsis thaliana) IRX10 involved in xylan extension, and the second most abundant was a GT61. Profiles for GT43 IRX9 and IRX14 putative orthologs were consistent with roles in AX synthesis. Low abundances were found for transcripts from genes in the acyl-coA transferase BAHD family, for which a role in AX feruloylation has been postulated. The relative expression of these was much greater in whole grain compared with starchy endosperm, correlating with the levels of bound ferulate. Transcripts associated with callose (GSL), cellulose (CESA), pectin (GAUT), and glucomannan (CSLA) synthesis were also abundant in starchy endosperm, while the corresponding cell wall polysaccharides were confirmed as low abundance (glucomannan and callose) or undetectable (pectin) in these samples. Abundant transcripts from GH families associated with the hydrolysis of these polysaccharides were also present, suggesting that they may be rapidly turned over. Abundant transcripts in the GT31 family may be responsible for the addition of Gal residues to arabinogalactan peptide.

Boosting the globalization of plant proteomics through INPPO: current developments 3 and future prospects

The International Plant Proteomics Organization (INPPO) is a non-profit-organization consisting of people who are involved or interested in plant proteomics. INPPO is constantly growing in volume and activity, which is mostly due to the realization among plant proteomics researchers worldwide for the need of such a global platform. Their active participation resulted in the rapid growth within the first year of INPPO’s official launch in 2011 via its website (www.inppo.com) and publication of the ‘viewpoint paper’ in a special issue of PROTEOMICS (May 2011). Here, we will be highlighting the progress achieved in the year 2011 and the future targets for the year 2012 and onwards. INPPO has achieved a successful administrative structure, the Core Committee (CC; composed of President, Vice-President, and General Secretaries), Executive Council (EC), and General Body (GB) toward achieving the INPPO objectives by its proposed initiatives. Various committees and subcommittees are in the process of being functionalized via discussion amongst scientists around the globe. INPPO’s primary aim to popularize the plant proteomics research in biological sciences has also been recognized by PROTEOMICS where a new section has been introduced to plant proteomics starting January 2012, following the very first issue of this journal devoted to plant proteomics in May 2011. To disseminate organizational activities to the scientific community, INPPO has launched a biannual (in January & July) newsletter entitled “INPPO Express: News & Views” with the first issue published in January 2012. INPPO is also planning to have several activities in 2012, including programs within the Education Outreach committee in different countries, and the development of research ideas and proposals with priority on crop and horticultural plants, while keeping tight interactions with proteomics programs on model plants such as Arabidopsis thaliana, rice, or Medicago truncatula. Altogether, the INPPO progress and upcoming activities are because of immense support, dedication, and hard work of all members of the INPPO family, and also due to the wide encouragement and support from the communities (scientific and non-scientific).

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.