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.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2782-3) contains supplementary material, which is available to authorized users.

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.

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.

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.

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 B₂) 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.