MetaCyc and AraCyc. Metabolic pathway databases for plant research

MetaCyc (http://metacyc.org) contains experimentally determined biochemical pathways to be used as a reference database for metabolism. In conjunction with the Pathway Tools software, MetaCyc can be used to computationally predict the metabolic pathway complement of an annotated genome. To increase the breadth of pathways and enzymes, more than 60 plant-specific pathways have been added or updated in MetaCyc recently. In contrast to MetaCyc, which contains metabolic data for a wide range of organisms, AraCyc is a species-specific database containing only enzymes and pathways found in the model plant Arabidopsis (Arabidopsis thaliana). AraCyc (http://arabidopsis.org/tools/aracyc/) was the first computationally predicted plant metabolism database derived from MetaCyc. Since its initial computational build, AraCyc has been under continued curation to enhance data quality and to increase breadth of pathway coverage. Twenty-eight pathways have been manually curated from the literature recently. Pathway predictions in AraCyc have also been recently updated with the latest functional annotations of Arabidopsis genes that use controlled vocabulary and literature evidence. AraCyc currently features 1,418 unique genes mapped onto 204 pathways with 1,156 literature citations. The Omics Viewer, a user data visualization and analysis tool, allows a list of genes, enzymes, or metabolites with experimental values to be painted on a diagram of the full pathway map of AraCyc. Other recent enhancements to both MetaCyc and AraCyc include implementation of an evidence ontology, which has been used to provide information on data quality, expansion of the secondary metabolism node of the pathway ontology to accommodate curation of secondary metabolic pathways, and enhancement of the cellular component ontology for storing and displaying enzyme and pathway locations within subcellular compartments.     

Network constrained clustering for gene microarray data

Many bioinformatics problems can be tackled from a fresh angle offered by the network perspective. Directly inspired by metabolic network structural studies, we propose an improved gene clustering approach for inferring gene signaling pathways from gene microarray data. Based on the construction of co-expression networks that consists of both significantly linear and non-linear gene associations together with controlled biological and statistical significance, our approach tends to group functionally related genes into tight clusters despite their expression dissimilarities. We illustrate our approach and compare it to the traditional clustering approaches on a yeast galactose metabolism dataset and a retinal gene expression dataset. Our approach greatly outperforms the traditional approach in rediscovering the relatively well known galactose metabolism pathway in yeast and in clustering genes of the photoreceptor differentiation pathway. AVAILABILITY: The clustering method has been implemented in an R package "GeneNT" that is freely available from: http://www.cran.org.     

Arabidopsis glucosyltransferase UGT74B1 functions in glucosinolate biosynthesis and auxin homeostasis

Glucosinolates are a class of secondary metabolites with important roles in plant defense and human nutrition. Here, we characterize a putative UDP-glucose:thiohydroximate S-glucosyltransferase, UGT74B1, to determine its role in the Arabidopsis glucosinolate pathway. Biochemical analyses demonstrate that recombinant UGT74B1 specifically glucosylates the thiohydroximate functional group. Low Km values for phenylacetothiohydroximic acid (approximately 6 microm) and UDP-glucose (approximately 50 microm) strongly suggest that thiohydroximates are in vivo substrates of UGT74B1. Insertional loss-of-function ugt74b1 mutants exhibit significantly decreased, but not abolished, glucosinolate accumulation. In addition, ugt74b1 mutants display phenotypes reminiscent of auxin overproduction, such as epinastic cotyledons, elongated hypocotyls in light-grown plants, excess adventitious rooting and incomplete leaf vascularization. Indeed, during early plant development, mutant ugt74b1 seedlings accumulate nearly threefold more indole-3-acetic acid than the wild type. Other phenotypes, however, such as chlorosis along the leaf veins, are likely caused by thiohydroximate toxicity. Analysis of UGT74B1 promoter activity during plant development reveals expression patterns consistent with glucosinolate metabolism and induction by auxin treatment. The results are discussed in the context of known mutations in glucosinolate pathway genes and their effects on auxin homeostasis. Taken together, our work provides complementary in vitro and in vivo evidence for a primary role of UGT74B1 in glucosinolate biosynthesis.     

Prediction of operon-like gene clusters in the Arabidopsis thaliana genome based on co-expression analysis of neighboring genes

Operon-like arrangements of genes occur in eukaryotes ranging from yeasts and filamentous fungi to nematodes, plants, and mammals. In plants, several examples of operon-like gene clusters involved in metabolic pathways have recently been characterized, e.g. the cyclic hydroxamic acid pathways in maize, the avenacin biosynthesis gene clusters in oat, the thalianol pathway in Arabidopsis thaliana, and the diterpenoid momilactone cluster in rice. Such operon-like gene clusters are defined by their co-regulation or neighboring positions within immediate vicinity of chromosomal regions. A comprehensive analysis of the expression of neighboring genes therefore accounts a crucial step to reveal the complete set of operon-like gene clusters within a genome. Genome-wide prediction of operon-like gene clusters should contribute to functional annotation efforts and provide novel insight into evolutionary aspects acquiring certain biological functions as well. We predicted co-expressed gene clusters by comparing the Pearson correlation coefficient of neighboring genes and randomly selected gene pairs, based on a statistical method that takes false discovery rate (FDR) into consideration for 1469 microarray gene expression datasets of A. thaliana. We estimated that A. thaliana contains 100 operon-like gene clusters in total. We predicted 34 statistically significant gene clusters consisting of 3 to 22 genes each, based on a stringent FDR threshold of 0.1. Functional relationships among genes in individual clusters were estimated by sequence similarity and functional annotation of genes. Duplicated gene pairs (determined based on BLAST with a cutoff of E<10(-5)) are included in 27 clusters. Five clusters are associated with metabolism, containing P450 genes restricted to the Brassica family and predicted to be involved in secondary metabolism. Operon-like clusters tend to include genes encoding bio-machinery associated with ribosomes, the ubiquitin/proteasome system, secondary metabolic pathways, lipid and fatty-acid metabolism, and the lipid transfer system.     

Glucose conjugation of anthranilate by the Arabidopsis UGT74F2 glucosyltransferase is required for tryptophan mutant blue fluorescence

Plant mutants with defects in intermediate enzymes of the tryptophan biosynthetic pathway often display a blue fluorescent phenotype. This phenotype results from the accumulation of the fluorescent tryptophan precursor anthranilate, the bulk of which is found in a glucose-conjugated form. To elucidate factors that control fluorescent tryptophan metabolites, we conducted a genetic screen for suppressors of blue fluorescence in the Arabidopsis trp1-100 mutant, which has a defect in the second enzymatic step of the tryptophan pathway. This screen yielded loss-of-function mutations in the UDP-glucosyltransferase gene UGT74F2. The bacterially expressed UGT74F2 enzyme catalyzed a conjugation reaction, with free anthranilate and UDP-glucose as substrates, that yielded the same fluorescent glucose ester compound as extracted from the trp1-100 mutant. These results indicate that sugar conjugation of anthranilate by UGT74F2 allows its stable accumulation in plant tissues. A highly related Arabidopsis enzyme UGT74F1 could also catalyze this reaction in vitro and could complement the ugt74F2 mutation when overexpressed in vivo. However, the UGT74F1 gene is expressed at a lower level than the UGT74F2 gene. Therefore, even though UGT74F1 and UGT74F2 have redundant conjugating activities toward anthranilate, UGT74F2 is the major source of this activity in the plant.     

THE EVOLUTION OF GENES IN BRANCHED METABOLIC PATHWAYS

Simulation models of the evolution of genes in a branched metabolic pathway subject to stabilizing selection on flux are described and analyzed. The models are based either on metabolic control theory (MCT), with the assumption that enzymes are far from saturation, or on Michaelis–Menten kinetics, which allows for saturation and near saturation. Several predictions emerge from the models: (1) flux control evolves to be concentrated at pathway branch points, including the first enzyme in the pathway. (2) When flux is far from its optimum, adaptive substitutions occur disproportionately often in branching enzymes. (3) When flux is near its optimum, adaptive substitutions occur disproportionately often in nonbranching enzymes. (4) Slightly deleterious substitutions occur disproportionately often in nonbranching enzymes. (5) In terms of both flux control and patterns of substitution, pathway branches are similar to those predicted for linear pathways. These predictions provide null hypotheses for empirical examination of the evolution of genes in metabolic pathways.     

Comparative analysis of Arabidopsis UGT74 glucosyltransferases reveals a special role of UGT74C1 in glucosinolate biosynthesis

The study of glucosinolates and their regulation has provided a powerful framework for the exploration of fundamental questions about the function, evolution, and ecological significance of plant natural products, but uncertainties about their metabolism remain. Previous work has identified one thiohydroximate S‐glucosyltransferase, UGT74B1, with an important role in the core pathway, but also made clear that this enzyme functions redundantly and cannot be the sole UDP‐glucose dependent glucosyltransferase (UGT) in glucosinolate synthesis. Here, we present the results of a nearly comprehensive in vitro activity screen of recombinant Arabidopsis Family 1 UGTs, which implicate other members of the UGT74 clade as candidate glucosinolate biosynthetic enzymes. Systematic genetic analysis of this clade indicates that UGT74C1 plays a special role in the synthesis of aliphatic glucosinolates, a conclusion strongly supported by phylogenetic and gene expression analyses. Finally, the ability of UGT74C1 to complement phenotypes and chemotypes of the ugt74b1‐2 knockout mutant and to express thiohydroximate UGT activity in planta provides conclusive evidence for UGT74C1 being an accessory enzyme in glucosinolate biosynthesis with a potential function during plant adaptation to environmental challenge.     

Engineering of red cells of <Emphasis Type="Italic">Arabidopsis thaliana</Emphasis> and comparative genome-wide gene expression analysis of red cells versus wild-type cells

We report metabolic engineering of Arabidopsis red cells and genome-wide gene expression analysis associated with anthocyanin biosynthesis and other metabolic pathways between red cells and wild-type (WT) cells. Red cells of A. thaliana were engineered for the first time from the leaves of production of anthocyanin pigment 1-Dominant (pap1-D). These red cells produced seven anthocyanin molecules including a new one that was characterized by LC–MS analysis. Wild-type cells established as a control did not produce anthocyanins. A genome-wide microarray analysis revealed that nearly 66 and 65% of genes in the genome were expressed in the red cells and wild-type cells, respectively. In comparison with the WT cells, 3.2% of expressed genes in the red cells were differentially expressed. The expression levels of 14 genes involved in the biosynthetic pathway of anthocyanin were significantly higher in the red cells than in the WT cells. Microarray and RT-PCR analyses demonstrated that the TTG1–GL3/TT8–PAP1 complex regulated the biosynthesis of anthocyanins. Furthermore, most of the genes with significant differential expression levels in the red cells versus the WT cells were characterized with diverse biochemical functions, many of which were mapped to different metabolic pathways (e.g., ribosomal protein biosynthesis, photosynthesis, glycolysis, glyoxylate metabolism, and plant secondary metabolisms) or organelles (e.g., chloroplast). We suggest that the difference in gene expression profiles between the two cell lines likely results from cell types, the overexpression of PAP1, and the high metabolic flux toward anthocyanins.     

Secondary metabolism, inventive evolution and biochemical diversity--a review.

Evidence now being obtained through nucleotide (nt) sequence analysis supports the concept that secondary metabolism has arisen by modification of existing primary metabolic reactions. Although amino acid sequence identity deduced from nt sequences of genes encoding proteins from related primary and secondary metabolic pathways is sufficient to indicate a common ancestry, the match is often better when genes in different rather than in the same species are compared. The information so far available suggests that gene transfer between organisms has been an important factor in the evolution of secondary metabolism. Many secondary pathways may be of relatively ancient origin and they may have arisen only infrequently. Much subsequent elaboration of the pathways has probably taken place after their acquisition by other species and so has been influenced by a variety of selective conditions. The characteristic diversity of secondary metabolites and their functions can be accounted for by the random manner in which the pathways initially evolved and have subsequently been exploited.     

Reciprocal regulation ofArabidopsis UGT78D2 and BANYULS is critical for regulation of the metabolic flux of anthocyanidins to condensed tannins in developing seed coats

Anthocyanins are major color pigments in plants. Their biosynthetic pathways are well established, and the majority of these biosynthetic enzymes have been identified in model plants such asArabidopsis, maize, and petunia. One exception inArabidopsis is UDP-glucose:flavonoid 3-O-glucosyltransferase (UF3GT). This enzyme is known as Bronze-1 (Bz1 ) in maize, where it converts anthocyanidins to anthocyanins. Phylogenetic sequence analysis of theArabidopsis thaliana UDP-glycosyltransferase (UGT) family previously indicated that UGT78D1, UGT78D2, and UGT78D3 cluster together with UF3GTs from other species. Here, we report thatUGT78D2 encodes a cytosolic UGT that is functionally consistent with maize Bz-1. Biochemically, UGT78D2 catalyzes the glucosylation of both flavonols and anthocyanidins at the 3-OH position. A T-DNA-insertedugt78d2 mutant accumulates very little anthocyanin and lacks 3-O-glucosylated quercetin. Expression analysis indicated thatUGT78D2, in opposite toBANYULS, is highly expressed in anthocyanin-accumulating seedlings but repressed in condensed tannin-accumulating seed coats. This suggests that the reciprocal regulation of these two genes is important in directing the metabolic flux to either anthocyanins or condensed tannins. Consistent with this, the ectopic expression of UGT78D2 produces purple-colored seed coats due to the accumulation of anthocyanins. Taken together, our data indicate thatUGT78D2 encodes an enzyme equivalent to maize Bz1, and that the reciprocal regulation of UGT78D2 and BANYULS is critical for the regulation of metabolic flux of anthocyanidins inArabidopsis.     

Genome organization in Arabidopsis thaliana: a survey for genes involved in isoprenoid and chlorophyll metabolism

The isoprenoid biosynthetic pathway provides intermediates for the synthesis of a multitude of natural products which serve numerous biochemical functions in plants: sterols (isoprenoids with a C30 backbone) are essential components of membranes; carotenoids (C40) and chlorophylls (which contain a C20 isoprenoid side-chain) act as photosynthetic pigments; plastoquinone, phylloquinone and ubiquinone (all of which contain long isoprenoid side-chains) participate in electron transport chains; gibberellins (C20), brassinosteroids (C30) and abscisic acid (C15) are phytohormones derived from isoprenoid intermediates; prenylation of proteins (with C15 or C20 isoprenoid moieties) may mediate subcellular targeting and regulation of activity; and several monoterpenes (C10), sesquiterpenes (C15) and diterpenes (C20) have been demonstrated to be involved in plant defense. Here we present a comprehensive analysis of genes coding for enzymes involved in the metabolism of isoprenoid-derived compounds in Arabidopsis thaliana. By combining homology and sequence motif searches with knowledge regarding the phylogenetic distribution of pathways of isoprenoid metabolism across species, candidate genes for these pathways in A. thaliana were obtained. A detailed analysis of the vicinity of chromosome loci for genes of isoprenoid metabolism in A. thaliana provided evidence for the clustering of genes involved in common pathways. Multiple sequence alignments were used to estimate the number of genes in gene families and sequence relationship trees were utilized to classify their individual members. The integration of all these datasets allows the generation of a knowledge-based metabolic map of isoprenoid metabolic pathways in A. thaliana and provides a substantial improvement of the currently available gene annotation.     

PathAligner

MOTIVATION: Analysis of metabolic pathways is a central topic in understanding the relationship between genotype and phenotype. The rapid accumulation of biological data provides the possibility of studying metabolic pathways at both the genomic and the metabolic levels. Retrieving metabolic pathways from current biological data sources, reconstructing metabolic pathways from rudimentary pathway components, and aligning metabolic pathways with each other are major tasks. Our motivation was to develop a conceptual framework and computational system that allows the retrieval of metabolic pathway information and the processing of alignments to reveal the similarities between metabolic pathways. RESULTS: PathAligner extracts metabolic information from biological databases via the Internet and builds metabolic pathways with data sources of genes, sequences, enzymes, metabolites etc. It provides an easy-to-use interface to retrieve, display and manipulate metabolic information. PathAligner also provides an alignment method to compare the similarity between metabolic pathways. AVAILABILITY: PathAligner is available at http://bibiserv.techfak.uni-bielefeld.de/pathaligner.     

Analysis of metabolic networks using a pathway distance metric through linear programming

The solution of the shortest path problem in biochemical systems constitutes an important step for studies of their evolution. In this paper, a linear programming (LP) algorithm for calculating minimal pathway distances in metabolic networks is studied. Minimal pathway distances are identified as the smallest number of metabolic steps separating two enzymes in metabolic pathways. The algorithm deals effectively with circularity and reaction directionality. The applicability of the algorithm is illustrated by calculating the minimal pathway distances for Escherichia coli small molecule metabolism enzymes, and then considering their correlations with genome distance (distance separating two genes on a chromosome) and enzyme function (as characterised by enzyme commission number). The results illustrate the effectiveness of the LP model. In addition, the data confirm that propinquity of genes on the genome implies similarity in function (as determined by co-involvement in the same region of the metabolic network), but suggest that no correlation exists between pathway distance and enzyme function. These findings offer insight into the probable mechanism of pathway evolution.     

Expression profiling of <Emphasis Type="Italic">Streptomyces peucetius</Emphasis> metabolic genes using DNA microarray analysis

Doxorubicin (DXR) and daunorubicin (DNR) are anthracycline antibiotics produced by Streptomyces peucetius and widely used as cancer chemotherapeutic agents. To improve their productivity, regulation of DXR/DNR synthesis genes as well as central metabolic pathway genes must be understood more clearly. So far, studies have focused on DXR/DNR gene regulation. To investigate the correlation between the central metabolic pathway genes and DXR/DNR productivity, we selected 265 genes involved in glycolysis, fermentation, the citric acid cycle, butanoate metabolism, etc., and searched for their sequences in the S. peucetius genome by comparing gene sequences to those of Streptomyces coelicolor. The homologous genes were amplified by PCR and arrayed on glass microarray slides. Gene expression was monitored under two different growth media conditions, R2YE and NDYE. Genes involved in the production of malonyl-CoA and propionyl-CoA, the main precursors for doxorubicin synthesis, were mainly upregulated in NDYE media. Genes related to acetyl-CoA and the urea cycle were also upregulated. These changes in gene expression were confirmed by real-time RT-PCR.     

Identification and functional analysis of light-responsive unique genes and gene family members in rice

Functional redundancy limits detailed analysis of genes in many organisms. Here, we report a method to efficiently overcome this obstacle by combining gene expression data with analysis of gene-indexed mutants. Using a rice NSF45K oligo-microarray to compare 2-week-old light- and dark-grown rice leaf tissue, we identified 365 genes that showed significant 8-fold or greater induction in the light relative to dark conditions. We then screened collections of rice T-DNA insertional mutants to identify rice lines with mutations in the strongly light-induced genes. From this analysis, we identified 74 different lines comprising two independent mutant lines for each of 37 light-induced genes. This list was further refined by mining gene expression data to exclude genes that had potential functional redundancy due to co-expressed family members (12 genes) and genes that had inconsistent light responses across other publicly available microarray datasets (five genes). We next characterized the phenotypes of rice lines carrying mutations in ten of the remaining candidate genes and then carried out co-expression analysis associated with these genes. This analysis effectively provided candidate functions for two genes of previously unknown function and for one gene not directly linked to the tested biochemical pathways. These data demonstrate the efficiency of combining gene family-based expression profiles with analyses of insertional mutants to identify novel genes and their functions, even among members of multi-gene families.     

Functional co-expression of xenobiotic metabolizing enzymes, rat cytochrome P450 1A1 and UDP-glucuronosyltransferase 1A6, in yeast microsomes

Xenobiotic Phase I and Phase II reactions in hepatocytes occur sequentially and cooperatively during the metabolism of various chemical compounds including drugs. In order to investigate the sequential metabolism of 7-ethoxycoumarin (7EC) as model substrate in vitro, xenobiotic metabolizing enzymes, rat cytochrome P450 1A1 (P450 1A1) and UDP-glucuronosyltransferase 1A6 (UGT1A6) were co-expressed in Saccharomyces cerevisiae AH22. Rat P450 1A1 and yeast NADPH-P450 reductase were expressed on a multicopy plasmid (pGYR1) in the yeast. Rat UGT1A6 cDNA with a yeast alcohol dehydrogenase I promoter and terminator was integrated into yeast chromosomal DNA to achieve the stable expression. Co-expression of P450 1A1 and UGT1A6 in yeast microsomes was confirmed by immunoblot analysis. Protease treatment of the microsomes showed the correct topological orientation of UGT to the membranes. The metabolism of 7EC to 7-hydroxycoumarin (7HC) and its glucuronide in yeast microsomes was analyzed by reverse phase HPLC. In a co-expression system containing 7EC, NADPH and UDP-glucuronic acid, glucuronide formation was detected after a lag phase, following the accumulation of 7HC. In the case of P450 1A1 and UGT1A6, efficient coupling of hydroxylation and glucuronidation in 7EC metabolism was not observed in the co-expression system. This P450 and UGT co-expression system in yeast allows the sequential biotransformation of xenobiotics to be simulated in vitro.