A proposed framework for the description of plant metabolomics experiments and their results

The study of the metabolite complement of biological samples, known as metabolomics, is creating large amounts of data, and support for handling these data sets is required to facilitate meaningful analyses that will answer biological questions. We present a data model for plant metabolomics known as ArMet (architecture for metabolomics). It encompasses the entire experimental time line from experiment definition and description of biological source material, through sample growth and preparation to the results of chemical analysis. Such formal data descriptions, which specify the full experimental context, enable principled comparison of data sets, allow proper interpretation of experimental results, permit the repetition of experiments and provide a basis for the design of systems for data storage and transmission. The current design and example implementations are freely available (http://www.armet.org/). We seek to advance discussion and community adoption of a standard for metabolomics, which would promote principled collection, storage and transmission of experiment data.

Use of mass spectrometry for imaging metabolites in plants

We discuss and illustrate recent advances that have been made to image the distribution of metabolites among cells and tissues of plants using different mass spectrometry technologies. These technologies include matrix-assisted laser desorption ionization, desorption electrospray ionization, and secondary ion mass spectrometry. These are relatively new technological applications of mass spectrometry and they are providing highly spatially resolved data concerning the cellular distribution of metabolites. We discuss the advantages and limitations of each of these mass spectrometric methods, and provide a description of the technical barriers that are currently limiting the technology to the level of single-cell resolution. However, we anticipate that advances in the next few years will increase the resolving power of the technology to provide unprecedented data on the distribution of metabolites at the subcellular level, which will increase our ability to decipher new knowledge concerning the spatial organization of metabolic processes in plants.

Modern plant metabolomics: advanced natural product gene discoveries, improved technologies, and future prospects

Plant metabolomics has matured and modern plant metabolomics has accelerated gene discoveries and the elucidation of a variety of plant natural product biosynthetic pathways. This review covers the approximate period of 2000 to 2014, and highlights specific examples of the discovery and characterization of novel genes and enzymes associated with the biosynthesis of natural products such as flavonoids, glucosinolates, terpenoids, and alkaloids. Additional examples of the integration of metabolomics with genome-based functional characterizations of plant natural products that are important to modern pharmaceutical technology are also reviewed. This article also provides a substantial review of recent technical advances in mass spectrometry imaging, nuclear magnetic resonance imaging, integrated LC-MS-SPE-NMR for metabolite identifications, and X-ray crystallography of microgram quantities for structural determinations. The review closes with a discussion on the future prospects of metabolomics related to crop species and herbal medicine.

Reverse genetic characterization of cytosolic acetyl-CoA generation by ATP-citrate lyase in Arabidopsis

Acetyl-CoA provides organisms with the chemical flexibility to biosynthesize a plethora of natural products that constitute much of the structural and functional diversity in nature. Recent studies have characterized a novel ATP-citrate lyase (ACL) in the cytosol of Arabidopsis thaliana. In this study, we report the use of antisense RNA technology to generate a series of Arabidopsis lines with a range of ACL activity. Plants with even moderately reduced ACL activity have a complex, bonsai phenotype, with miniaturized organs, smaller cells, aberrant plastid morphology, reduced cuticular wax deposition, and hyperaccumulation of starch, anthocyanin, and stress-related mRNAs in vegetative tissue. The degree of this phenotype correlates with the level of reduction in ACL activity. These data indicate that ACL is required for normal growth and development and that no other source of acetyl-CoA can compensate for ACL-derived acetyl-CoA. Exogenous malonate, which feeds into the carboxylation pathway of acetyl-CoA metabolism, chemically complements the morphological and chemical alterations associated with reduced ACL expression, indicating that the observed metabolic alterations are related to the carboxylation pathway of cytosolic acetyl-CoA metabolism. The observations that limiting the expression of the cytosolic enzyme ACL reduces the accumulation of cytosolic acetyl-CoA-derived metabolites and that these deficiencies can be alleviated by exogenous malonate indicate that ACL is a nonredundant source of cytosolic acetyl-CoA.

Molecular characterization of a heteromeric ATP-citrate lyase that generates cytosolic acetyl-coenzyme A in Arabidopsis

Acetyl-coenzyme A (CoA) is used in the cytosol of plant cells for the synthesis of a diverse set of phytochemicals including waxes, isoprenoids, stilbenes, and flavonoids. The source of cytosolic acetyl-CoA is unclear. We identified two Arabidopsis cDNAs that encode proteins similar to the amino and carboxy portions of human ATP-citrate lyase (ACL). Coexpression of these cDNAs in yeast (Saccharomyces cerevisiae) confers ACL activity, indicating that both the Arabidopsis genes are required for ACL activity. Arabidopsis ACL is a heteromeric enzyme composed of two distinct subunits, ACLA (45 kD) and ACLB (65 kD). The holoprotein has a molecular mass of 500 kD, which corresponds to a heterooctomer with an A(4)B(4) configuration. ACL activity and the ACLA and ACLB polypeptides are located in the cytosol, consistent with the lack of targeting peptides in the ACLA and ACLB sequences. In the Arabidopsis genome, three genes encode for the ACLA subunit (ACLA-1, At1g10670; ACLA-2, At1g60810; and ACLA-3, At1g09430), and two genes encode the ACLB subunit (ACLB-1, At3g06650 and ACLB-2, At5g49460). The ACLA and ACLB mRNAs accumulate in coordinated spatial and temporal patterns during plant development. This complex accumulation pattern is consistent with the predicted physiological needs for cytosolic acetyl-CoA, and is closely coordinated with the accumulation pattern of cytosolic acetyl-CoA carboxylase, an enzyme using cytosolic acetyl-CoA as a substrate. Taken together, these results indicate that ACL, encoded by the ACLA and ACLB genes of Arabidopsis, generates cytosolic acetyl-CoA. The heteromeric organization of this enzyme is common to green plants (including Chlorophyceae, Marchantimorpha, Bryopsida, Pinaceae, monocotyledons, and eudicots), species of fungi, Glaucophytes, Chlamydomonas, and prokaryotes. In contrast, all known animal ACL enzymes have a homomeric structure, indicating that a evolutionary fusion of the ACLA and ACLB genes probably occurred early in the evolutionary history of this kingdom.

Phylogenetic and experimental characterization of an acyl-ACP thioesterase family reveals significant diversity in enzymatic specificity and activity

Background

Acyl-acyl carrier protein thioesterases (acyl-ACP TEs) catalyze the hydrolysis of the thioester bond that links the acyl chain to the sulfhydryl group of the phosphopantetheine prosthetic group of ACP. This reaction terminates acyl chain elongation of fatty acid biosynthesis, and in plant seeds it is the biochemical determinant of the fatty acid compositions of storage lipids.

Results

To explore acyl-ACP TE diversity and to identify novel acyl ACP-TEs, 31 acyl-ACP TEs from wide-ranging phylogenetic sources were characterized to ascertain their in vivo activities and substrate specificities. These acyl-ACP TEs were chosen by two different approaches: 1) 24 TEs were selected from public databases on the basis of phylogenetic analysis and fatty acid profile knowledge of their source organisms; and 2) seven TEs were molecularly cloned from oil palm (Elaeis guineensis), coconut (Cocos nucifera) and Cuphea viscosissima, organisms that produce medium-chain and short-chain fatty acids in their seeds. The in vivo substrate specificities of the acyl-ACP TEs were determined in E. coli. Based on their specificities, these enzymes were clustered into three classes: 1) Class I acyl-ACP TEs act primarily on 14- and 16-carbon acyl-ACP substrates; 2) Class II acyl-ACP TEs have broad substrate specificities, with major activities toward 8- and 14-carbon acyl-ACP substrates; and 3) Class III acyl-ACP TEs act predominantly on 8-carbon acyl-ACPs. Several novel acyl-ACP TEs act on short-chain and unsaturated acyl-ACP or 3-ketoacyl-ACP substrates, indicating the diversity of enzymatic specificity in this enzyme family.

Conclusion

These acyl-ACP TEs can potentially be used to diversify the fatty acid biosynthesis pathway to produce novel fatty acids.

The role of pyruvate dehydrogenase and acetyl-coenzyme A synthetase in fatty acid synthesis in developing Arabidopsis seeds

Acetyl-coenzyme A (acetyl-CoA) formed within the plastid is the precursor for the biosynthesis of fatty acids and, through them, a range of important biomolecules. The source of acetyl-CoA in the plastid is not known, but two enzymes are thought to be involved: acetyl-CoA synthetase and plastidic pyruvate dehydrogenase. To determine the importance of these two enzymes in synthesizing acetyl-CoA during lipid accumulation in developing Arabidopsis seeds, we isolated cDNA clones for acetyl-CoA synthetase and for the ptE1alpha- and ptE1beta-subunits of plastidic pyruvate dehydrogenase. To our knowledge, this is the first reported acetyl-CoA synthetase sequence from a plant source. The Arabidopsis acetyl-CoA synthetase preprotein has a calculated mass of 76,678 D, an apparent plastid targeting sequence, and the mature protein is a monomer of 70 to 72 kD. During silique development, the spatial and temporal patterns of the ptE1beta mRNA level are very similar to those of the mRNAs for the plastidic heteromeric acetyl-CoA carboxylase subunits. The pattern of ptE1beta mRNA accumulation strongly correlates with the formation of lipid within the developing embryo. In contrast, the level of mRNA for acetyl-CoA synthetase does not correlate in time and space with lipid accumulation. The highest level of accumulation of the mRNA for acetyl-CoA synthetase during silique development is within the funiculus. These mRNA data suggest a predominant role for plastidic pyruvate dehydrogenase in acetyl-CoA formation during lipid synthesis in seeds.

High-spatial and high-mass resolution imaging of surface metabolites of Arabidopsis thaliana by laser desorption-ionization mass spectrometry using colloidal silver

High-spatial resolution and high-mass resolution techniques are developed and adopted for the mass spectrometric imaging of epicuticular lipids on the surface of Arabidopsis thaliana. Single cell level spatial resolution of approximately 12 mum was achieved by reducing the laser beam size by using an optical fiber with 25 mum core diameter in a vacuum matrix-assisted laser desorption ionization-linear ion trap (vMALDI-LTQ) mass spectrometer and improved matrix application using an oscillating capillary nebulizer. Fine chemical images of a whole flower were visualized in this high spatial resolution showing substructure of an anther and single pollen grains at the stigma and anthers. The LTQ-Orbitrap with a MALDI ion source was adopted to achieve MS imaging in high mass resolution. Specifically, isobaric silver ion adducts of C29 alkane (m/z 515.3741) and C28 aldehyde (m/z 515.3377), indistinguishable in low-resolution LTQ, can now be clearly distinguished and their chemical images could be separately constructed. In the application to roots, the high spatial resolution allowed molecular MS imaging of secondary roots and the high mass resolution allowed direct identification of lipid metabolites on root surfaces.

Platform biochemicals for a biorenewable chemical industry

The chemical industry is currently reliant on a historically inexpensive, petroleum-based carbon feedstock that generates a small collection of platform chemicals from which highly efficient chemical conversions lead to the manufacture of a large variety of chemical products. Recently, a number of factors have coalesced to provide the impetus to explore alternative renewable sources of carbon. Here we discuss the potential impact on the chemical industry of shifting from non-renewable carbon sources to renewable carbon sources. This change to the manufacture of chemicals from biological carbon sources will provide an opportunity for the biological research community to contribute fundamental knowledge concerning carbon metabolism and its regulation. We discuss whether fundamental biological research into metabolic processes at a holistic level, made possible by completed genome sequences and integrated with detailed structural understanding of biocatalysts, can change the chemical industry from being dependent on fossil-carbon feedstocks to using biorenewable feedstocks. We illustrate this potential by discussing the prospect of building a platform technology based upon a concept of combinatorial biosynthesis, which would explore the enzymological flexibilities of polyketide biosynthesis.

Molecular characterization of meiotic recombination across the 140-kb multigenic a1-sh2 interval of maize

The 140-kb a1-sh2 interval of the maize genome contains at least four genes (a1, yz1, x1, and sh2). Partial sequence analysis of two haplotypes has revealed many single nucleotide polymorphisms and InDel polymorphisms, including several large structural polymorphisms. The physical positions of 101 meiotic recombination breakpoints are not distributed uniformly across the interval and are instead concentrated within three recombination hot spots. Two of these recombination hot spots are genic (a1 and yz1) and one is apparently nongenic. The x1 gene is not a recombination hot spot. Thus, these results suggest that not all hot spots are genes and indicate that not all genes are hot spots. Two of the 101 recombination events arose by means of either noncrossover events involving conversion tract lengths of at least 17 kb or double-crossover events. Only one recombination breakpoint mapped to the approximately 80-kb distal portion of the a1-sh2 interval that contains large amounts of repetitive DNA including retrotransposons; in this region the ratio of genetic to physical distance is less than 0.5% of the genome's average. These results establish that the retrotransposon faction of the maize genome is relatively inert recombinationally.

The relationship between genetic and physical distances in the cloned a1-sh2 interval of the Zea mays L. genome

A 470-kb segment from the long arm of chromosome 3 of Zea mays (inbred LH82), encompassing the a1-sh2 interval, was cloned as a yeast artificial chromosome. Comparison of the sizes of the restriction fragments generated from the cloned DNA fragment and from the DNA isolated from the maize inbred line LH82 established the colinearity of the a1-sh2 interval in these DNAs. By utilizing a chromosome fragmentation technique, a yeast artificial chromosome encompassing the a1-sh2 interval was separately fragmented at the a1 and sh2 loci. Comparison of the sizes of these fragmentation products established the physical distance between the a1 and sh2 loci to be 140 kb. Furthermore, these fragmentation experiments established the physical orientation of the a1 and sh2 genes relative to the maize centromere. The molecular cloning of the contiguous region between the a1 and sh2 loci made it possible to define the relationship between physical and genetic distances over a relatively large segment of the maize genome. In this interval, the relationship between physical and genetic distances is 1560 kb/centimorgan, which compares with 1460 kb/centimorgan for the entire maize genome, and 217 kb/centimorgan for a 1-kb segment within the a1 locus. Therefore, these findings are consistent with the hypothesis that genes per se are preferred sites for meiotic recombination rather than the hypothesis that genes reside in large recombinationally active segments of the genome.

A global approach to analysis and interpretation of metabolic data for plant natural product discovery

Discovering molecular components and their functionality is key to the development of hypotheses concerning the organization and regulation of metabolic networks. The iterative experimental testing of such hypotheses is the trajectory that can ultimately enable accurate computational modelling and prediction of metabolic outcomes. This information can be particularly important for understanding the biology of natural products, whose metabolism itself is often only poorly defined. Here, we describe factors that must be in place to optimize the use of metabolomics in predictive biology. A key to achieving this vision is a collection of accurate time-resolved and spatially defined metabolite abundance data and associated metadata. One formidable challenge associated with metabolite profiling is the complexity and analytical limits associated with comprehensively determining the metabolome of an organism. Further, for metabolomics data to be efficiently used by the research community, it must be curated in publicly available metabolomics databases. Such databases require clear, consistent formats, easy access to data and metadata, data download, and accessible computational tools to integrate genome system-scale datasets. Although transcriptomics and proteomics integrate the linear predictive power of the genome, the metabolome represents the nonlinear, final biochemical products of the genome, which results from the intricate system(s) that regulate genome expression. For example, the relationship of metabolomics data to the metabolic network is confounded by redundant connections between metabolites and gene-products. However, connections among metabolites are predictable through the rules of chemistry. Therefore, enhancing the ability to integrate the metabolome with anchor-points in the transcriptome and proteome will enhance the predictive power of genomics data. We detail a public database repository for metabolomics, tools and approaches for statistical analysis of metabolomics data, and methods for integrating these datasets with transcriptomic data to create hypotheses concerning specialized metabolisms that generate the diversity in natural product chemistry. We discuss the importance of close collaborations among biologists, chemists, computer scientists and statisticians throughout the development of such integrated metabolism-centric databases and software.

Maize Mu transposons are targeted to the 5' untranslated region of the gl8 gene and sequences flanking Mu target-site duplications exhibit nonrandom nucleotide composition throughout the genome

The widespread use of the maize Mutator (Mu) system to generate mutants exploits the preference of Mu transposons to insert into genic regions. However, little is known about the specificity of Mu insertions within genes. Analysis of 79 independently isolated Mu-induced alleles at the gl8 locus established that at least 75 contain Mu insertions. Analysis of the terminal inverted repeats (TIRs) of the inserted transposons defined three new Mu transposons: Mu10, Mu 11, and Mu12. A large percentage (>80%) of the insertions are located in the 5' untranslated region (UTR) of the gl8 gene. Ten positions within the 5' UTR experienced multiple independent Mu insertions. Analyses of the nucleotide composition of the 9-bp TSD and the sequences directly flanking the TSD reveals that the nucleotide composition of Mu insertion sites differs dramatically from that of random DNA. In particular, the frequencies at which C's and G's are observed at positions -2 and +2 (relative to the TSD) are substantially higher than expected. Insertion sites of 315 RescueMu insertions displayed the same nonrandom nucleotide composition observed for the gl8-Mu alleles. Hence, this study provides strong evidence for the involvement of sequences flanking the TSD in Mu insertion-site selection.

Plants contain multiple biotin enzymes: discovery of 3-methylcrotonyl-CoA carboxylase, propionyl-CoA carboxylase and pyruvate carboxylase in the plant kingdom

Acetyl-CoA carboxylase is the sole biotin enzyme previously reported in plants. Western analysis with 125I-streptavidin of proteins extracted from carrot somatic embryos visualized six biotin-containing polypeptides, the relative molecular masses of which are 210,000, 140,000, 73,000, 50,000, 39,000, and 34,000. This multiplicity of the biotin-containing polypeptides can be partly explained by the discovery of 3-methylcrotonyl-CoA carboxylase, propionyl-CoA carboxylase, and pyruvate carboxylase in extracts of somatic carrot embryos, biotin enzymes previously unknown in the plant kingdom. These biotin enzymes seem to be widely distributed in the plant kingdom.

Rosmarinic acid in Prunella vulgaris ethanol extract inhibits lipopolysaccharide-induced prostaglandin E2 and nitric oxide in RAW 264.7 mouse macrophages

Prunella vulgaris has been used therapeutically for inflammation-related conditions for centuries, but systematic studies of its anti-inflammatory activity are lacking and no specific active components have been identified. In this study, water and ethanol extracts of four P. vulgaris accessions were applied to RAW 264.7 mouse macrophages, and the ethanol extracts significantly inhibited lipopolysaccharide (LPS)-stimulated prostaglandin E2 (PGE2) and nitric oxide (NO) production at 30 microg/mL without affecting cell viability. Extracts from different accessions of P. vulgaris were screened for anti-inflammatory activity to identify accessions with the greatest activity. The inhibition of PGE2 and NO production by selected extracts was dose-dependent, with significant effects seen at concentrations as low as 10 microg/mL. Fractionation of ethanol extracts from the active accession, Ames 27664, suggested fractions 3 and 5 as possible major contributors to the overall activity. Rosmarinic acid (RA) content in P. vulgaris was found to independently inhibit inflammatory response, but it only partially explained the extracts' activity. LPS-induced cyclooxygenase-2 (COX-2) and nitric oxide synthase (iNOS) protein expression were both attenuated by P. vulgaris ethanol extracts, whereas RA inhibited only COX-2 expression.

Sequence analysis of the cloned glossy8 gene of maize suggests that it may code for a beta-ketoacyl reductase required for the biosynthesis of cuticular waxes

The gl8 locus of maize (Zea mays L.) was previously defined by a mutation that reduces the amount and alters the composition of seedling cuticular waxes. Sixty independently derived gl8 mutant alleles were isolated from stocks that carried the Mutator transposon system. A DNA fragment that contains a Mu8 transposon and that co-segregates with one of these alleles, gl8-Mu3142, was identified and cloned. DNA flanking the Mu8 transposon was shown via allelic cross-referencing experiments to represent the gl8 locus. The gl8 probe revealed a 1.4-kb transcript present in wild-type seedling leaves and, in lesser amounts, in other organs and at other developmental stages. The amino acid sequence deduced from an apparently full-length gl8 cDNA exhibits highly significant sequence similarity to a group of enzymes from plants, eubacteria, and mammals that catalyzes the reduction of ketones. This finding suggests that the GL8 protein probably functions as a reductase during fatty acid elongation in the cuticular wax biosynthetic pathway.

Characterization of two GL8 paralogs reveals that the 3-ketoacyl reductase component of fatty acid elongase is essential for maize (Zea mays L.) development

Prior analyses established that the maize (Zea mays L.) gl8a gene encodes 3-ketoacyl reductase, a component of the fatty acid elongase required for the biosynthesis of very long chain fatty acids (VLCFAs). A paralogous gene, gl8b, has been identified that is 96% identical to gl8a. The gl8a and gl8b genes map to syntenic chromosomal regions, have similar, but not identical, expression patterns, and encode proteins that are 97% identical. Both of these genes are required for the normal accumulation of cuticular waxes on seedling leaves. The chemical composition of the cuticular waxes from gl8a and gl8b mutants indicates that these genes have at least overlapping, if not redundant, functions in cuticular wax biosynthesis. Although gl8a and gl8b double mutant kernels have endosperms that cannot be distinguished from wild-type siblings, these kernels are non-viable because their embryos fail to undergo normal development. Double mutant kernels accumulate substantially reduced levels of VLCFAs. VLCFAs are components of a variety of compounds, for example, cuticular waxes, suberin, and sphingolipids. Consistent with their essential nature in yeast, the accumulation of the ceramide moiety of sphingolipids is substantially reduced and their fatty acid composition altered in gl8a and gl8b double mutant kernels relative to wild-type kernels. Hence, we hypothesize that sphingolipids or other VLCFA-containing compounds are essential for normal embryo development.

Articulation of three core metabolic processes in Arabidopsis: fatty acid biosynthesis, leucine catabolism and starch metabolism

BACKGROUND

Elucidating metabolic network structures and functions in multicellular organisms is an emerging goal of functional genomics. We describe the co-expression network of three core metabolic processes in the genetic model plant Arabidopsis thaliana: fatty acid biosynthesis, starch metabolism and amino acid (leucine) catabolism.

RESULTS

These co-expression networks form modules populated by genes coding for enzymes that represent the reactions generally considered to define each pathway. However, the modules also incorporate a wider set of genes that encode transporters, cofactor biosynthetic enzymes, precursor-producing enzymes, and regulatory molecules. We tested experimentally the hypothesis that one of the genes tightly co-expressed with starch metabolism module, a putative kinase AtPERK10, will have a role in this process. Indeed, knockout lines of AtPERK10 have an altered starch accumulation. In addition, the co-expression data define a novel hierarchical transcript-level structure associated with catabolism, in which genes performing smaller, more specific tasks appear to be recruited into higher-order modules with a broader catabolic function.

CONCLUSION

Each of these core metabolic pathways is structured as a module of co-expressed transcripts that co-accumulate over a wide range of environmental and genetic perturbations and developmental stages, and represent an expanded set of macromolecules associated with the common task of supporting the functionality of each metabolic pathway. As experimentally demonstrated, co-expression analysis can provide a rich approach towards understanding gene function.

Transcriptional and post-transcriptional regulation of ribulose 1,5-bisphosphate carboxylase gene expression in light- and dark-grown amaranth cotyledons

The regulation of expression of the genes encoding the large subunit (LSU) and small subunit (SSU) of ribulose 1,5-bisphosphate carboxylase (RuBPCase) was examined in 1- through 8-day-old, dark-grown (etiolated) and light-grown amaranth cotyledons. RuBPCase specific activity in light-grown cotyledons increased during this 8-day period to a level 15-fold higher than in dark-grown cotyledons. Under both growth conditions, the accumulation of the LSU and SSU polypeptides was not coordinated. Initial detection of the SSU occurred 1 and 2 days after the appearance of the LSU in light- and dark-grown cotyledons, respectively. Furthermore, although the levels of the LSU were similar in both light- and dark-grown seedlings, the amount of the SSU followed clearly the changes in enzyme activity. Synthesis of these two polypeptides was dramatically different in etiolated versus light-grown cotyledons. In light the synthesis of both subunits was first observed on day 2 and continued throughout the growth of the cotyledons. In darkness the rate of synthesis of both subunits was much lower than in light and occurred only as a burst between days 2 and 5 after planting. However, mRNAs for both subunits were present in etiolated cotyledons at similar levels on days 4 through 7 (by Northern analysis) and were functional in vitro, despite their apparent inactivity in vivo after day 5. In addition, since both LSU and SSU mRNA levels were lower in dark- than in light-grown seedlings, our results indicate that both transcriptional and post-transcriptional controls modulate RuBPCase production in developing amaranth cotyledons.

PIN6 is required for nectary auxin response and short stamen development

The PIN family of proteins is best known for its involvement in polar auxin transport and tropic responses. PIN6 (At1g77110) is one of the remaining PIN family members in Arabidopsis thaliana to which a biological function has not yet been ascribed. Here we report that PIN6 is a nectary-enriched gene whose expression level is positively correlated with total nectar production in Arabidopsis, and whose function is required for the proper development of short stamens. PIN6 accumulates in internal membranes consistent with the ER, and multiple lines of evidence demonstrate that PIN6 is required for auxin-dependent responses in nectaries. Wild-type plants expressing auxin-responsive DR5:GFP or DR5:GUS reporters displayed intense signal in lateral nectaries, but pin6 lateral nectaries showed little or no signal for these reporters. Further, exogenous auxin treatment increased nectar production more than tenfold in wild-type plants, but nectar production was not increased in pin6 mutants when treated with auxin. Conversely, the auxin transport inhibitor N-1-naphthylphthalamic acid (NPA) reduced nectar production in wild-type plants by more than twofold, but had no significant effect on pin6 lines. Interestingly, a MYB57 transcription factor mutant, myb57-2, closely phenocopied the loss-of-function mutant pin6-2. However, PIN6 expression was not dependent on MYB57, and RNA-seq analyses of pin6-2 and myb57-2 mutant nectaries showed little overlap in terms of differentially expressed genes. Cumulatively, these results demonstrate that PIN6 is required for proper auxin response and nectary function in Arabidopsis. These results also identify auxin as an important factor in the regulation of nectar production, and implicate short stamens in the maturation of lateral nectaries.

Genetic recombination in plants

Meiotic recombination generates novel allelic arrays on chromosomes. Recent experiments have revealed an extraordinarily nonrandom distribution of recombination breakpoints along the lengths of plant chromosomes; for example, recombination breakpoints often resolve within genic sequences, and thereby generate novel alleles. The mechanism by which recombination breakpoints are determined is an area of active investigation. In addition, recent developments are providing recombination-based technologies for creating targeted alterations in the architecture of plant genomes.

Metabolomics as a Hypothesis-Generating Functional Genomics Tool for the Annotation of Arabidopsis thaliana Genes of "Unknown Function"

Metabolomics is the methodology that identifies and measures global pools of small molecules (of less than about 1,000 Da) of a biological sample, which are collectively called the metabolome. Metabolomics can therefore reveal the metabolic outcome of a genetic or environmental perturbation of a metabolic regulatory network, and thus provide insights into the structure and regulation of that network. Because of the chemical complexity of the metabolome and limitations associated with individual analytical platforms for determining the metabolome, it is currently difficult to capture the complete metabolome of an organism or tissue, which is in contrast to genomics and transcriptomics. This paper describes the analysis of Arabidopsis metabolomics data sets acquired by a consortium that includes five analytical laboratories, bioinformaticists, and biostatisticians, which aims to develop and validate metabolomics as a hypothesis-generating functional genomics tool. The consortium is determining the metabolomes of Arabidopsis T-DNA mutant stocks, grown in standardized controlled environment optimized to minimize environmental impacts on the metabolomes. Metabolomics data were generated with seven analytical platforms, and the combined data is being provided to the research community to formulate initial hypotheses about genes of unknown function (GUFs). A public database (www.PlantMetabolomics.org) has been developed to provide the scientific community with access to the data along with tools to allow for its interactive analysis. Exemplary datasets are discussed to validate the approach, which illustrate how initial hypotheses can be generated from the consortium-produced metabolomics data, integrated with prior knowledge to provide a testable hypothesis concerning the functionality of GUFs.