A Core Metabolome Response of Maize Leaves Subjected to Long-Duration Abiotic Stresses.

Abiotic stresses reduce crop growth and yield in part by disrupting metabolic homeostasis and triggering responses that change the metabolome. Experiments designed to understand the mechanisms underlying these metabolomic responses have usually not used agriculturally relevant stress regimes. We therefore subjected maize plants to drought, salt, or heat stresses that mimic field conditions and analyzed leaf responses at metabolome and transcriptome levels. Shared features of stress metabolomes included synthesis of raffinose, a compatible solute implicated in tolerance to dehydration. In addition, a marked accumulation of amino acids including proline, arginine, and γ-aminobutyrate combined with depletion of key glycolysis and tricarboxylic acid cycle intermediates indicated a shift in balance of carbon and nitrogen metabolism in stressed leaves. Involvement of the γ-aminobutyrate shunt in this process is consistent with its previously proposed role as a workaround for stress-induced thiamin-deficiency. Although convergent metabolome shifts were correlated with gene expression changes in affected pathways, patterns of differential gene regulation induced by the three stresses indicated distinct signaling mechanisms highlighting the plasticity of plant metabolic responses to abiotic stress.

Corrigendum: Divisions of labor in the thiamin biosynthetic pathway among organs of maize.

[This corrects the article on p. 370 in vol. 5, PMID: 25136345.].

Arabidopsis TH2 Encodes the Orphan Enzyme Thiamin Monophosphate Phosphatase.

To synthesize the cofactor thiamin diphosphate (ThDP), plants must first hydrolyze thiamin monophosphate (ThMP) to thiamin, but dedicated enzymes for this hydrolysis step were unknown and widely doubted to exist. The classical thiamin-requiring th2-1 mutation in Arabidopsis thaliana was shown to reduce ThDP levels by half and to increase ThMP levels 5-fold, implying that the THIAMIN REQUIRING2 (TH2) gene product could be a dedicated ThMP phosphatase. Genomic and transcriptomic data indicated that TH2 corresponds to At5g32470, encoding a HAD (haloacid dehalogenase) family phosphatase fused to a TenA (thiamin salvage) family protein. Like the th2-1 mutant, an insertional mutant of At5g32470 accumulated ThMP, and the thiamin requirement of the th2-1 mutant was complemented by wild-type At5g32470 Complementation tests in Escherichia coli and enzyme assays with recombinant proteins confirmed that At5g32470 and its maize (Zea mays) orthologs GRMZM2G148896 and GRMZM2G078283 are ThMP-selective phosphatases whose activity resides in the HAD domain and that the At5g32470 TenA domain has the expected thiamin salvage activity. In vitro and in vivo experiments showed that alternative translation start sites direct the At5g32470 protein to the cytosol and potentially also to mitochondria. Our findings establish that plants have a dedicated ThMP phosphatase and indicate that modest (50%) ThDP depletion can produce severe deficiency symptoms.

Systematic identification and analysis of frequent gene fusion events in metabolic pathways.

BACKGROUND: Gene fusions are the most powerful type of in silico-derived functional associations. However, many fusion compilations were made when <100 genomes were available, and algorithms for identifying fusions need updating to handle the current avalanche of sequenced genomes. The availability of a large fusion dataset would help probe functional associations and enable systematic analysis of where and why fusion events occur. RESULTS: Here we present a systematic analysis of fusions in prokaryotes. We manually generated two training sets: (i) 121 fusions in the model organism Escherichia coli; (ii) 131 fusions found in B vitamin metabolism. These sets were used to develop a fusion prediction algorithm that captured the training set fusions with only 7 % false negatives and 50 % false positives, a substantial improvement over existing approaches. This algorithm was then applied to identify 3.8 million potential fusions across 11,473 genomes. The results of the analysis are available in a searchable database at http://modelseed.org/projects/fusions/ . A functional analysis identified 3,000 reactions associated with frequent fusion events and revealed areas of metabolism where fusions are particularly prevalent. CONCLUSIONS: Customary definitions of fusions were shown to be ambiguous, and a stricter one was proposed. Exploring the genes participating in fusion events showed that they most commonly encode transporters, regulators, and metabolic enzymes. The major rationales for fusions between metabolic genes appear to be overcoming pathway bottlenecks, avoiding toxicity, controlling competing pathways, and facilitating expression and assembly of protein complexes. Finally, our fusion dataset provides powerful clues to decipher the biological activities of domains of unknown function.

Experimental and Metabolic Modeling Evidence for a Folate-Cleaving Side-Activity of Ketopantoate Hydroxymethyltransferase (PanB).

Tetrahydrofolate (THF) and its one-carbon derivatives, collectively termed folates, are essential cofactors, but are inherently unstable. While it is clear that chemical oxidation can cleave folates or damage their pterin precursors, very little is known about enzymatic damage to these molecules or about whether the folate biosynthesis pathway responds adaptively to damage to its end-products. The presence of a duplication of the gene encoding the folate biosynthesis enzyme 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase (FolK) in many sequenced bacterial genomes combined with a strong chromosomal clustering of the folK gene with panB, encoding the 5,10-methylene-THF-dependent enzyme ketopantoate hydroxymethyltransferase, led us to infer that PanB has a side activity that cleaves 5,10-methylene-THF, yielding a pterin product that is recycled by FolK. Genetic and metabolic analyses of Escherichia coli strains showed that overexpression of PanB leads to accumulation of the likely folate cleavage product 6-hydroxymethylpterin and other pterins in cells and medium, and-unexpectedly-to a 46% increase in total folate content. In silico modeling of the folate biosynthesis pathway showed that these observations are consistent with the in vivo cleavage of 5,10-methylene-THF by a side-activity of PanB, with FolK-mediated recycling of the pterin cleavage product, and with regulation of folate biosynthesis by folates or their damage products.

Bacterial and plant HAD enzymes catalyse a missing phosphatase step in thiamin diphosphate biosynthesis.

The penultimate step of thiamin diphosphate (ThDP) synthesis in plants and many bacteria is dephosphorylation of thiamin monophosphate (ThMP). Non-specific phosphatases have been thought to mediate this step and no genes encoding specific ThMP phosphatases (ThMPases) are known. Comparative genomic analysis uncovered bacterial haloacid dehalogenase (HAD) phosphatase family genes (from subfamilies IA and IB) that cluster on the chromosome with, or are fused to, thiamin synthesis genes and are thus candidates for the missing phosphatase (ThMPase). Three typical candidates (from Anaerotruncus colihominis, Dorea longicatena and Syntrophomonas wolfei) were shown to have efficient in vivo ThMPase activity by expressing them in an Escherichia coli strain engineered to require an active ThMPase for growth. In vitro assays confirmed that these candidates all preferred ThMP to any of 45 other phosphate ester substrates tested. An Arabidopsis thaliana ThMPase homologue (At4g29530) of unknown function whose expression pattern and compartmentation fit with a role in ThDP synthesis was shown to have in vivo ThMPase activity in E. coli and to prefer ThMP to any other substrate tested. However, insertional inactivation of the At4g29530 gene did not affect growth or the levels of thiamin or its phosphates, indicating that Arabidopsis has at least one other ThMPase gene. The Zea mays orthologue of At4g29530 (GRMZM2G035134) was also shown to have ThMPase activity. These data identify HAD genes specifying the elusive ThMPase activity, indicate that ThMPases are substrate-specific rather than general phosphatases and suggest that different evolutionary lineages have recruited ThMPases independently from different branches of the HAD family.

Functional Diversity of Haloacid Dehalogenase Superfamily Phosphatases from Saccharomyces cerevisiae: BIOCHEMICAL, STRUCTURAL, AND EVOLUTIONARY INSIGHTS.

The haloacid dehalogenase (HAD)-like enzymes comprise a large superfamily of phosphohydrolases present in all organisms. The Saccharomyces cerevisiae genome encodes at least 19 soluble HADs, including 10 uncharacterized proteins. Here, we biochemically characterized 13 yeast phosphatases from the HAD superfamily, which includes both specific and promiscuous enzymes active against various phosphorylated metabolites and peptides with several HADs implicated in detoxification of phosphorylated compounds and pseudouridine. The crystal structures of four yeast HADs provided insight into their active sites, whereas the structure of the YKR070W dimer in complex with substrate revealed a composite substrate-binding site. Although the S. cerevisiae and Escherichia coli HADs share low sequence similarities, the comparison of their substrate profiles revealed seven phosphatases with common preferred substrates. The cluster of secondary substrates supporting significant activity of both S. cerevisiae and E. coli HADs includes 28 common metabolites that appear to represent the pool of potential activities for the evolution of novel HAD phosphatases. Evolution of novel substrate specificities of HAD phosphatases shows no strict correlation with sequence divergence. Thus, evolution of the HAD superfamily combines the conservation of the overall substrate pool and the substrate profiles of some enzymes with remarkable biochemical and structural flexibility of other superfamily members.

A directed-overflow and damage-control N-glycosidase in riboflavin biosynthesis.

Plants and bacteria synthesize the essential human micronutrient riboflavin (vitamin B2) via the same multi-step pathway. The early intermediates of this pathway are notoriously reactive and may be overproduced in vivo because riboflavin biosynthesis enzymes lack feedback controls. In the present paper, we demonstrate disposal of riboflavin intermediates by COG3236 (DUF1768), a protein of previously unknown function that is fused to two different riboflavin pathway enzymes in plants and bacteria (RIBR and RibA respectively). We present cheminformatic, biochemical, genetic and genomic evidence to show that: (i) plant and bacterial COG3236 proteins cleave the N-glycosidic bond of the first two intermediates of riboflavin biosynthesis, yielding relatively innocuous products; (ii) certain COG3236 proteins are in a multi-enzyme riboflavin biosynthesis complex that gives them privileged access to riboflavin intermediates; and (iii) COG3236 action in Arabidopsis thaliana and Escherichia coli helps maintain flavin levels. COG3236 proteins thus illustrate two emerging principles in chemical biology: directed overflow metabolism, in which excess flux is diverted out of a pathway, and the pre-emption of damage from reactive metabolites.

Divisions of labor in the thiamin biosynthetic pathway among organs of maize.

The B vitamin thiamin is essential for central metabolism in all cellular organisms including plants. While plants synthesize thiamin de novo, organs vary widely in their capacities for thiamin synthesis. We use a transcriptomics approach to appraise the distribution of de novo synthesis and thiamin salvage pathways among organs of maize. We identify at least six developmental contexts in which metabolically active, non-photosynthetic organs exhibit low expression of one or both branches of the de novo thiamin biosynthetic pathway indicating a dependence on inter-cellular transport of thiamin and/or thiamin precursors. Neither the thiazole (THI4) nor pyrimidine (THIC) branches of the pathway are expressed in developing pollen implying a dependence on import of thiamin from surrounding floral and inflorescence organs. Consistent with that hypothesis, organs of the male inflorescence and flowers are shown to have high relative expression of the thiamin biosynthetic pathway and comparatively high thiamin contents. By contrast, divergent patterns of THIC and THI4 expression occur in the shoot apical meristem, embyro sac, embryo, endosperm, and root-tips suggesting that these sink organs acquire significant amounts of thiamin via salvage pathways. In the root and shoot meristems, expression of THIC in the absence of THI4 indicates a capacity for thiamin synthesis via salvage of thiazole, whereas the opposite pattern obtains in embryo and endosperm implying that seed storage organs are poised for pyrimidine salvage. Finally, stable isotope labeling experiments set an upper limit on the rate of de novo thiamin biosynthesis in maize leaf explants. Overall, the observed patterns of thiamin biosynthetic gene expression mirror the strategies for thiamin acquisition that have evolved in bacteria.

High-throughput comparison, functional annotation, and metabolic modeling of plant genomes using the PlantSEED resource.

The increasing number of sequenced plant genomes is placing new demands on the methods applied to analyze, annotate, and model these genomes. Today's annotation pipelines result in inconsistent gene assignments that complicate comparative analyses and prevent efficient construction of metabolic models. To overcome these problems, we have developed the PlantSEED, an integrated, metabolism-centric database to support subsystems-based annotation and metabolic model reconstruction for plant genomes. PlantSEED combines SEED subsystems technology, first developed for microbial genomes, with refined protein families and biochemical data to assign fully consistent functional annotations to orthologous genes, particularly those encoding primary metabolic pathways. Seamless integration with its parent, the prokaryotic SEED database, makes PlantSEED a unique environment for cross-kingdom comparative analysis of plant and bacterial genomes. The consistent annotations imposed by PlantSEED permit rapid reconstruction and modeling of primary metabolism for all plant genomes in the database. This feature opens the unique possibility of model-based assessment of the completeness and accuracy of gene annotation and thus allows computational identification of genes and pathways that are restricted to certain genomes or need better curation. We demonstrate the PlantSEED system by producing consistent annotations for 10 reference genomes. We also produce a functioning metabolic model for each genome, gapfilling to identify missing annotations and proposing gene candidates for missing annotations. Models are built around an extended biomass composition representing the most comprehensive published to date. To our knowledge, our models are the first to be published for seven of the genomes analyzed.

A cross-kingdom Nudix enzyme that pre-empts damage in thiamin metabolism.

Genes specifying the thiamin monophosphate phosphatase and adenylated thiazole diphosphatase steps in fungal and plant thiamin biosynthesis remain unknown, as do genes for ThDP (thiamin diphosphate) hydrolysis in thiamin metabolism. A distinctive Nudix domain fused to Tnr3 (thiamin diphosphokinase) in Schizosaccharomyces pombe was evaluated as a candidate for these functions. Comparative genomic analysis predicted a role in thiamin metabolism, not biosynthesis, because free-standing homologues of this Nudix domain occur not only in fungi and plants, but also in proteobacteria (whose thiamin biosynthesis pathway has no adenylated thiazole or thiamin monophosphate hydrolysis steps) and animals (which do not make thiamin). Supporting this prediction, recombinant Tnr3 and its Saccharomyces cerevisiae, Arabidopsis and maize Nudix homologues lacked thiamin monophosphate phosphatase activity, but were active against ThDP, and up to 60-fold more active against diphosphates of the toxic thiamin degradation products oxy- and oxo-thiamin. Deleting the S. cerevisiae Nudix gene (YJR142W) lowered oxythiamin resistance, overexpressing it raised resistance, and expressing its plant or bacterial counterparts restored resistance to the YJR142W deletant. By converting the diphosphates of damaged forms of thiamin into monophosphates, the Tnr3 Nudix domain and its homologues can pre-empt the misincorporation of damaged diphosphates into ThDP-dependent enzymes, and the resulting toxicity.

CathaCyc, a metabolic pathway database built from Catharanthus roseus RNA-Seq data.

The medicinal plant Madagascar periwinkle (Catharanthus roseus) synthesizes numerous terpenoid indole alkaloids (TIAs), such as the anticancer drugs vinblastine and vincristine. The TIA pathway operates in a complex metabolic network that steers plant growth and survival. Pathway databases and metabolic networks reconstructed from 'omics' sequence data can help to discover missing enzymes, study metabolic pathway evolution and, ultimately, engineer metabolic pathways. To date, such databases have mainly been built for model plant species with sequenced genomes. Although genome sequence data are not available for most medicinal plant species, next-generation sequencing is now extensively employed to create comprehensive medicinal plant transcriptome sequence resources. Here we report on the construction of CathaCyc, a detailed metabolic pathway database, from C. roseus RNA-Seq data sets. CathaCyc (version 1.0) contains 390 pathways with 1,347 assigned enzymes and spans primary and secondary metabolism. Curation of the pathways linked with the synthesis of TIAs and triterpenoids, their primary metabolic precursors, and their elicitors, the jasmonate hormones, demonstrated that RNA-Seq resources are suitable for the construction of pathway databases. CathaCyc is accessible online (http://www.cathacyc.org) and offers a range of tools for the visualization and analysis of metabolic networks and 'omics' data. Overlay with expression data from publicly available RNA-Seq resources demonstrated that two well-characterized C. roseus terpenoid pathways, those of TIAs and triterpenoids, are subject to distinct regulation by both developmental and environmental cues. We anticipate that databases such as CathaCyc will become key to the study and exploitation of the metabolism of medicinal plants.

Identification and characterization of the missing pyrimidine reductase in the plant riboflavin biosynthesis pathway.

Riboflavin (vitamin B2) is the precursor of the flavin coenzymes flavin mononucleotide and flavin adenine dinucleotide. In Escherichia coli and other bacteria, sequential deamination and reduction steps in riboflavin biosynthesis are catalyzed by RibD, a bifunctional protein with distinct pyrimidine deaminase and reductase domains. Plants have two diverged RibD homologs, PyrD and PyrR; PyrR proteins have an extra carboxyl-terminal domain (COG3236) of unknown function. Arabidopsis (Arabidopsis thaliana) PyrD (encoded by At4g20960) is known to be a monofunctional pyrimidine deaminase, but no pyrimidine reductase has been identified. Bioinformatic analyses indicated that plant PyrR proteins have a catalytically competent reductase domain but lack essential zinc-binding residues in the deaminase domain, and that the Arabidopsis PyrR gene (At3g47390) is coexpressed with riboflavin synthesis genes. These observations imply that PyrR is a pyrimidine reductase without deaminase activity. Consistent with this inference, Arabidopsis or maize (Zea mays) PyrR (At3g47390 or GRMZM2G090068) restored riboflavin prototrophy to an E. coli ribD deletant strain when coexpressed with the corresponding PyrD protein (At4g20960 or GRMZM2G320099) but not when expressed alone; the COG3236 domain was unnecessary for complementing activity. Furthermore, recombinant maize PyrR mediated NAD(P)H-dependent pyrimidine reduction in vitro. Import assays with pea (Pisum sativum) chloroplasts showed that PyrR and PyrD are taken up and proteolytically processed. Ablation of the maize PyrR gene caused early seed lethality. These data argue that PyrR is the missing plant pyrimidine reductase, that it is plastid localized, and that it is essential. The role of the COG3236 domain remains mysterious; no evidence was obtained for the possibility that it catalyzes the dephosphorylation that follows pyrimidine reduction.

Evidence that the folate-dependent proteins YgfZ and MnmEG have opposing effects on growth and on activity of the iron-sulfur enzyme MiaB.

The folate-dependent protein YgfZ of Escherichia coli participates in the synthesis and repair of iron-sulfur (Fe-S) clusters; it belongs to a family of enzymes that use folate to capture formaldehyde units. Ablation of ygfZ is known to reduce growth, to increase sensitivity to oxidative stress, and to lower the activities of MiaB and other Fe-S enzymes. It has been reported that the growth phenotype can be suppressed by disrupting the tRNA modification gene mnmE. We first confirmed the latter observation using deletions in a simpler, more defined genetic background. We then showed that deleting mnmE substantially restores MiaB activity in ygfZ deletant cells and that overexpressing MnmE with its partner MnmG exacerbates the growth and MiaB activity phenotypes of the ygfZ deletant. MnmE, with MnmG, normally mediates a folate-dependent transfer of a formaldehyde unit to tRNA, and the MnmEG-mediated effects on the phenotypes of the ΔygfZ mutant apparently require folate, as evidenced by the effect of eliminating all folates by deleting folE. The expression of YgfZ was unaffected by deleting mnmE or overexpressing MnmEG or by folate status. Since formaldehyde transfer is a potential link between MnmEG and YgfZ, we inactivated formaldehyde detoxification by deleting frmA. This deletion had little effect on growth or MiaB activity in the ΔygfZ strain in the presence of formaldehyde, making it unlikely that formaldehyde alone connects the actions of MnmEG and YgfZ. A more plausible explanation is that MnmEG erroneously transfers a folate-bound formaldehyde unit to MiaB and that YgfZ reverses this.

Mutational analysis of YgfZ, a folate-dependent protein implicated in iron/sulphur cluster metabolism.

Proteins of the YgfZ family occur in all domains of life and are characterized by the conserved dodecapeptide motif KGC[Y/F]-x-GQE-x(3) -[R/K]. YgfZ proteins are known to participate in assembly or repair of iron/sulphur clusters, and to require folate for biological activity, but their mechanism of action is unknown. To assess the importance of individual residues in the conserved motif, Escherichia coli Ygf Z was expressed from a plasmid in a ΔygfZ strain and subjected to alanine-scanning mutagenesis. The impacts on YgfZ functionality were evaluated by assays of growth and of the in vivo activity of the iron/sulphur enzyme MiaB, which modifies tRNA. By these criteria, the motif's tyrosine residue (Y229) had a detectable influence but only the cysteine residue (C228) was critical, for only the C228A mutant failed to complement the growth and MiaB activity phenotypes of the ΔygfZ strain. Immunoblots confirmed that the latter result was not simply because of a low level of the C228A mutant protein. Collectively, these data demonstrate a pivotal role for the Ygf Z motif's cysteine residue and a subsidiary one for the adjacent tyrosine, and help formulate a hypothesis about the folate requirement of Ygf Z proteins.