Bioinformatics Advance Access originally published online on March 7, 2007
Bioinformatics 2007 23(11):1418-1423; doi:10.1093/bioinformatics/btm040
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
MedicCyc: a biochemical pathway database for Medicago truncatula
The Samuel Roberts Noble Foundation, Inc., Ardmore, OK 73401, USA
*To whom correspondence should be addressed.
| ABSTRACT |
|---|
|
|
|---|
Motivation: There is an imperative need to integrate functional genomics data to obtain a more comprehensive systems-biology view of the results. We believe that this is best achieved through the visualization of data within the biological context of metabolic pathways. Accordingly, metabolic pathway reconstruction was used to predict the metabolic composition for Medicago truncatula and these pathways were engineered to enable the correlated visualization of integrated functional genomics data.
Results: Metabolic pathway reconstruction was used to generate a pathway database for M. truncatula (MedicCyc), which currently features more than 250 pathways with related genes, enzymes and metabolites. MedicCyc was assembled from more than 225 000 M. truncatula ESTs (MtGI Release 8.0) and available genomic sequences using the Pathway Tools software and the MetaCyc database. The predicted pathways in MedicCyc were verified through comparison with other plant databases such as AraCyc and RiceCyc. The comparison with other plant databases provided crucial information concerning enzymes still missing from the ongoing, but currently incomplete M. truncatula genome sequencing project. MedicCyc was further manually curated to remove non-plant pathways, and Medicago-specific pathways including isoflavonoid, lignin and triterpene saponin biosynthesis were modified or added based upon available literature and in-house expertise. Additional metabolites identified in metabolic profiling experiments were also used for pathway predictions. Once the metabolic reconstruction was completed, MedicCyc was engineered to visualize M. truncatula functional genomics datasets within the biological context of metabolic pathways.
Availability: freely accessible at http://www.noble.org/MedicCyc/
Contact: lwsumner{at}noble.org
| 1 INTRODUCTION |
|---|
|
|
|---|
A comprehensive and accurate understanding of both primary and secondary metabolism is fundamental to the knowledge of plant biology. Although primary metabolism is one of the best reconstructed networks within biological systems, our understanding of secondary metabolism and its metabolic regulation is still quite limited. Improved functional genomic technologies such as transcriptomics, proteomics and metabolomics offer stimulating opportunities for a more thorough understanding of metabolism and we currently utilize an integrated functional genomics approach to study the model legume Medicago truncatula (Achnine et al., 2005; Broeckling et al., 2005; Lei et al., 2005; Sumner et al., 2003; Suzuki et al., 2005). It is presumed that the model legume M. truncatula will be the fourth plant genome to be fully sequenced (Town, 2006; Young et al., 2005), and as of January 2007,
190 Mbp of euchromatin sequence has been acquired (http://compbio.dfci.harvard.edu/tgi/cgi-bin/tgi/gimain.pl?gudb=medicago) representing approximately 25 000 predicted genes. At the same time, the collection of ESTs (http://www.tigr.org/tigr-scripts/tgi/T_index.cgi?species=medicago) has increased to 18 238 non-redundant singleton ESTs and 30 000 non-redundant tentative consensus (TC) sequences (MtGI Release 8.0). With the continually mounting sequence of M. truncatula, release of the fully sequenced Arabidopsis (Arabidopsis Genome Initiative, 2000), rice (International Rice Genome Sequencing Project, 2005), and black cottonwood (Tuskan et al., 2006) genomes, and the initiation of many other plant sequencing projects; there is a fast growing need for a better understanding of the functional aspects of sequenced genomes. This demand can be best met through an integrated functional genomics approach, which typically includes a combination of transcriptomics, proteomics and/or metabolomics. Genomic and expressed tag sequencing of M. truncatula (Young et al., 2005) have provided resources for the assimilation of custom and commercial cDNA and oligonucleotide high-density microarray chip sets for transcriptome analysis (Barnett et al., 2004; Covitz et al., 1998; Kuster et al., 2004). Similarly, M. truncatula nucleotide sequencing and advances in mass spectrometry have enabled large-scale profiling of proteins (Gallardo et al., 2003; Imin et al., 2004; Lei et al., 2005; Mathesius et al., 2001; Watson et al., 2003). The development of metabolomics is more recent (Fiehn, 2002; Sumner et al., 2003; Weckwerth, 2003), and the number and variety of metabolites that can be routinely measured have remarkably increased over the last several years (Broeckling et al., 2005; Chen et al., 2003; Huhman and Sumner, 2002; Huhman et al., 2005; Liu et al., 2003). Biochemical pathway maps composed of genes, proteins and metabolites (Zhang et al., 2005) are powerful tools around which one can compile the biological context of functional genomics datasets. Currently, several legume databases exist and include The Legume Information System (http://www.comparative-legumes.org/) (Gonzales et al., 2005), The Medicago truncatula Consortium Database (MtDB, http://www.medicago.org/MtDB/) (Lamblin et al., 2003) and the Medicago Gene Index (http://compbio.dfci.harvard.edu/tgi/cgi-bin/tgi/gimain.pl?gudb=medicago)) (Cannon et al., 2005). Of these, only the TIGR database provides linkages of tentative consensus sequences to metabolic pathways, and to the best of our knowledge, none of the databases listed are currently capable of visualizing functional genomics data (e.g. expression data obtained by microarray analysis) within the context of biological pathways. However, such databases do exist for other species. For example, the visualization tool MapMan has been developed for the display of metabolite and transcript data onto metabolic pathways of Arabidopsis and tomato (Thimm et al., 2004; Urbanczyk-Wochniak et al., 2006; Usadel et al., 2005). Similarly, a commercial software has been modified to visualize Arabidopsis transcript and metabolite data (Lange and Ghassemian, 2005). To the best of our knowledge, there is no dedicated tool or publicly accessible database designed for accommodating or visualizing M. truncatula functional genomics data. As a result, we have developed MedicCyc, a biochemical pathway database for M. truncatula which can be used as a reference for M. truncatula as well as other legume species and is freely accessible at http://www.noble.org/MedicCyc/.
MedicCyc was constructed using MetaCyc, a database composed of all experimentally determined biochemical pathways for small molecule metabolism (Caspi et al., 2006; Keseler et al., 2005; Krieger et al., 2004), and the related Pathway Tools software (Karp et al., 2002; Paley and Karp, 2002). MetaCyc consists of reactions, enzymes, metabolites, feedback regulation information and genes that encode enzymes accumulated from a large number of species to enable the prediction of pathways in new species. The current version of MetaCyc (10.0) contains approximately 700 pathways from more than 600 species, ranging from microbes to plants and human (http://metacyc.org/). The Pathway Tools software allows automated generation of a pathways database through functional assignment of genes and manual editing of pathways using a graphical user interface. Currently, more than 90% of pathways integrated into MetaCyc are manually curated and include literature citation and species information (Caspi et al., 2006). To date only two pathway genome databases (PGDBs) have been created for plants. The first was AraCyc for Arabidopsis thaliana (http://www.arabidopsis.org/tools/aracyc/) (Mueller et al., 2003; Zhang et al., 2005) and more recently, the Gramene Community (a Resource for Comparative Grass Genomics) reconstructed the biochemical pathways for rice and created RiceCyc—a web-based tool for viewing gene annotations mapped to various biochemical pathways in two model plant species, Arabidopsis and rice respectively (for more details visit http://www.gramene.org/pathway/). The utilization of a common Cyc database format provides a consistent and more equitable platform for the comparison of reconstructed plant pathways between Medicago and other available PlantCyc databases. This is possible through the utilization of a Pathway Tools comparative module. Pathway comparisons between plant species can also reveal current gaps in the knowledge of Medicago metabolism.
MedicCyc, a curated and integrated pathway database for M. truncatula, was developed to enable further study of metabolism and functional genomics in M. truncatula. The goals of this project were: (i) to use pathway reconstruction for predicting the metabolic composition of M. truncatula; (ii) to provide a platform for visualization of the integrated functional genomics datasets. Furthermore, the long-term goals are: (iii) to enable the prediction and annotation of unknown genes through comparison with other sequenced plant genomes and detected metabolites; and (iv) to provide an additional resource for the metabolic engineering of legume metabolism and the enhancements of desired traits.
| 2 METHODS |
|---|
|
|
|---|
2.1. MedicCyc construction using transcript annotation
MedicCyc was constructed using metabolic annotated transcripts selected from The Affymetrix GeneChip® Medicago Genome Array. The Medicago Genome Array was designed specifically to monitor gene expression in Madicago truncatula, Medicago sativa and additionally the symbiotic organism Sinorhizobium meliloti (http://www.affymetrix.com/products/arrays/specific/medicago.affx). Sequence information for this array was selected from diverse data sources including the TIGR M. truncatula gene index (MtGI Release 8.0), automated gene predictions from The International Medicago Genome Annotation Group (IMGAG, http://www.medicago.org/genome/IMGAG.php) (Town, 2006), gene predictions from the S. meliloti genome, and M.sativa EST information released by TIGR. The array consists of more than 32 167 Medicago truncatula EST/mRNA-based and chloroplast gene-based probe sets as well as 18 733 M. truncatula IMGAG and phase 2/3 BAC prediction-based probe sets, which are believed to cover the majority of the transcribed part of the M. truncatula genome. It has been predicted that
25% of the Arabidopsis genome encodes metabolic enzymes and the value may be similar or even higher for legumes (Bell et al., 2001). Accordingly, 4400 transcripts were pre-selected from Medicago Genome Array based upon the corresponding annotations that suggest that these genes are involved in metabolism. Due to the current and incomplete stage of M. truncatula genomic sequencing, genes and corresponding transcripts have not been completely associated with chromosomes. Therefore, MedicCyc transcripts were associated as arbitrary groups which do not represent the typical chromosomal grouping utilized by Pathologic. After completion of the genomic sequencing, all transcripts used to build MedicCyc will be linked and grouped to specific chromosomes. The analog annotations were formatted into a PathoLogic-specific set according to the documentation for Pathway Tools (Keseler et al., 2005) and used for the first MedicCyc database construction. The annotations were edited manually to remove extraneous words and characters that could interfere with the enzyme name-matching software, which was analogous to the process used for construction of AraCyc (Mueller et al., 2003). Enzymes labeled as putative or similar to were also included in the dataset. Additionally, the selected transcripts were coupled with corresponding EC (Enzyme Commission) numbers, which were used to map the transcripts onto pathways and to avoid problems with enzyme synonyms. The initial MedicCyc database was generated using the PathoLogic Pathway Prediction module of Pathway Tools ver. 10.0. The initial Medicago pathways (MedicCyc ver. beta) were inferred using MetaCyc as a reference database of metabolic pathways, and later MedicCyc ver. 1.0 was curated using AraCyc and RiceCyc PGDBs as co-reference databases.
2.2 MedicCyc construction using previously detected metabolites
To evaluate the completeness of the M. truncatula metabolic reconstruction, we compared the compounds in MedicCyc with compounds previously detected and identified in metabolic profiling experiments (Broeckling et al., 2005; Chen et al., 2003; Huhman and Sumner, 2002; Huhman et al., 2005; Liu et al., 2003). Forty-two compounds reported in these studies were not found in the initial version of MedicCyc (Table 1). Missing compounds included secondary metabolites recently reported in M. truncatula, e.g. different isoflavonoids (Deavours et al., 2006; Farag et al., 2007). The absent metabolites were manually inserted into the database and most of them were used to reconstruct new metabolic pathways, e.g. triterpene saponin biosynthesis. Continuing metabolic profiling experiments will further extend these pathways as our knowledge of secondary metabolite biosynthesis and degradation is still far from complete.
|
| 3 RESULTS |
|---|
|
|
|---|
3.1 MedicCyc statistics
After completing the initial automated construction of MedicCyc (MedicCyc ver. beta), PathoLogic recognized 3418 enzymes for which the defined functions were known, and another 707 gene products identified as putative enzymes. The putative enzymes consisted of generic annotations such as MAP kinase and the precise functions of these enzymes are still unknown. In total, the initial MedicCyc ver. beta database contained 408 pathways (Table 2), which were comprised of 1678 enzymatic reactions and 1232 compounds. It was anticipated that the automated construction would include non-plant pathways because the MetaCyc reference database includes a large number of non-plant species. Accordingly, 126 pathways that could not be confirmed in plants were deleted from the first Pathways Tools analysis.
|
Seventeen pathways were manually reconstructed (Table 3) to compensate for incomplete pathways retrieved from MetaCyc (e.g. superpathway of phenylalanine, tyrosine and tryptophan biosynthesis I). Several pathway fragments were observed in the beta version of MedicCyc due to the lack of Medicago annotations. These included mainly primary metabolic pathways, which are presumed to be conserved across plant species based upon comparisons of the predicted Medicago pathways with pathways from Arabidopsis and rice, i.e. AraCyc and RiceCyc, respectively. Accordingly, the partially reconstructed Medicago pathways noted in Table 3 were extended by copying complete pathways directly from MetaCyc. A few selected pathways related to secondary metabolism were altered using specific knowledge about Medicago metabolism or based upon additional data supporting the presence of novel pathways in Medicago, e.g. detection and identification of the pathway metabolites. Specific new pathway variants included pathways related to isoflavonoid biosynthesis, lignin biosynthesis or anthocyanin biosynthesis.
|
The final, curated version of MediCyc (1.0) consists of 254 pathways (Table 2), 3426 enzymes and 1215 compounds. Most of the reconstructed pathways are grouped into biosynthesis classes, which are similar to Arabidopsis and rice (Table 4). The number of enzymatic reactions in different MedicCyc pathways varied. However, a different number of Medicago annotated genes were automatically implemented into single pathways, depending upon their biosynthesis class. The most complete reconstructed pathways belong to Energy metabolism, which is analogous to the situation described for Arabidopsis (Mueller et al., 2003). Other pathways such as fatty acids and lipids biosynthesis or amines and polyamines degradation contain only single Medicago genes (e.g. cutin biosynthesis, choline biosynthesis or ureide degradation). It is anticipated that further genomic sequencing and metabolic profiling experiments will provide additional information for verifying MedicCyc completeness as we probe deeper into the metabolome of M. truncatula and as more sensitive detection techniques become available.
|
Biochemical pathway reconstruction provides not only information about predicted pathways, but also about selected reactions for which no corresponding enzymes have been assigned yet. These reactions are referred to as pathway holes, and are detected by the Pathway Tools software (Keseler et al., 2005). So far 587 pathways holes were identified in MedicCyc. The list of pathway holes (all missing enzymes) is available on the MedicCyc web page (http://www.noble.org/MedicCyc/). It is presumed that a significant portion of the missing enzymes will be revealed as the M. truncatula genomic sequencing project moves closer towards completion.
3.2 Functional genomics data visualization
It is useful to visualize functional genomics data within a biochemical context for the interpretation of large-scale omics data. Thus, reconstructed Medicago pathways can also be used in data analysis via the Omics Viewer, a software module of Pathways Tools. The Omics Viewer—can visualize data from gene expression, protein expression, gene family analysis or metabolite profiling experiments within the biological context of metabolic networks (Zhang et al., 2005). Each reaction (represented as a connection line between the compounds) can be color coded in accordance with the expression level of the gene and the protein catalyzing the reaction (Fig. 1). As we used expressed sequence information from the Medicago Affymetrix Chipset, MedicCyc was engineered to visualize experimental transcript data specifically from this platform. The results obtained using alternative mircroarray platforms can also be accommodated assuming each transcript element is represented with a unique nomenclature (TIGR ID or IMGAG ID). The MedicCyc Omics Viewer was also designed to display metabolic profiling data. Each detected metabolite, which was previously implemented into MedicCyc database, can be color coded in accordance with its recorded abundance. Correlated metabolite and expression data can then be viewed in the form of a general overview map (Fig. 1A) or as in the expanded single pathway view (Fig. 1B and C). The visualization viewer is available through the MedicCyc web site which enables visualization of user's own data.
|
| 4 DISCUSSION AND CONCLUSION |
|---|
|
|
|---|
The aims of this project were to use pathway reconstruction to predict the metabolic composition of M. truncatula and further exploit this information as a platform for the visualization of integrated functional genomics datasets. This was achieved through the construction of MedicCyc. This database was manually curated and the final curated version of MedicCyc (1.0) consists of 254 pathways, 3426 enzymes and 1215 compounds (Table 2). A significant number of enzymes were identified as still missing from the currently incomplete M. truncatula genome sequencing project through comparisons with AraCyc. As the M. truncatula genome sequencing project moves toward completion, MedicCyc will be updated. The MedicCyc database is currently supported by The Noble Foundation and additional funding for continued curation and improvements will be sought. Continued curation will be necessary to address new transcript annotations and new metabolic pathways related to Medicago metabolism which become available during on-going sequencing and annotation of the M. truncatula genome. Additional future efforts will focus on the inclusion of subcellular localizations for specific known enzymatic isoforms and metabolites. Continual metabolic profiling experiments are expected to provide additional information concerning novel metabolic content in M. truncatula and will be added as available. It is anticipated that this resource will serve as a repository for our current and future understanding of Medicago metabolism, provide a fundamental tool for the visualization of functional genomics datasets, serve as a contributor to larger databases (MetaCyc, MtDB, LIS, etc.), facilitate comparative studies of pathways across species, enable the prediction and annotation of unknown genes, and facilitate metabolic engineering of legume metabolism for the enhancement of desired traits. We strongly encourage input from the scientific community through the incorporated tool available online.
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
We thank Ms Keisha Wood for assistance with the database curation and Mr Roland Stolfa for technical assistance in establishing the MedicCyc database. We also thank Drs Patrick Zao and Xinbin Dai for consultation, and we express our appreciation to Dr Yuhong Tang for providing unpublished Affymetrix microarray data. This work was financially supported by The Samuel Roberts Noble Foundation. Funding to pay the Open Access publication charges was provided by The Samuel Roberts Noble Foundation.
Conflict of Interest: none declared.
| FOOTNOTES |
|---|
Associate Editor: Alex Bateman
Received on October 26, 2006; revised on January 29, 2007; accepted on February 1, 2007
| REFERENCES |
|---|
|
|
|---|
Achnine L, et al. Genomics-based selection and functional characterization of triterpene glycosyltransferases from the model legume. Medicago truncatula. Plant J. (2005) 41:875–887.[CrossRef][Web of Science][Medline]
Arabidopsis Genome Initiative. Analysis of the genome sequence of the flowering plant Arabidopsis thaliana. Nature (2000) 408:796–815.[CrossRef][Medline]
Barnett MJ, et al. A dual-genome Symbiosis Chip for coordinate study of signal exchange and development in a prokaryote-host interaction. PNAS (2004) 101:16636–16641.
Bell CJ, et al. The Medicago Genome Initiative: a model legume database. Nucleic Acids Res. (2001) 29:114–117.
Broeckling CD, et al. Metabolic profiling of Medicago truncatula cell cultures reveals the effects of biotic and abiotic elicitors on metabolism. J. Exp. Bot. (2005) 56:323–336.
Cannon SB, et al. Databases and information integration for the Medicago truncatula genome and transcriptome. Plant Physiol. (2005) 138:38–46.
Caspi R, et al. MetaCyc: a multiorganism database of metabolic pathways and enzymes. Nucleic Acids Res. (2006) 34:D511–D516.
Chen F, et al. Profiling phenolic metabolites in transgenic alfalfa modified in lignin biosynthesis. Phytochemistry (2003) 64:1013–1021.[CrossRef][Web of Science][Medline]
Covitz PA, et al. Expressed sequence tags from a root-hair-enriched Medicago truncatula cDNA library. Plant Physiol. (1998) 117:1325–1332.
Deavours BE, et al. Functional analysis of members of the isoflavone and isoflavanone O-methyltransferase enzyme families from the model legume Medicago truncatula. Plant Mol. Biol. (2006) 67:2262–2268.
Farag MA, et al. Metabolic profiling and systematic identification of flavonoids and isoflavonoids in roots and cell suspension cultures of Medicago truncatula using HPLC-UV-ESI-MS and GC-MS. Phytochemistry (2007) 68:342–354.[CrossRef][Web of Science][Medline]
Fiehn O. Metabolomics-the link between genotypes and phenotypes. Plant Mol. Biol. (2002) 48:155–171.[CrossRef][Web of Science][Medline]
Gallardo K, et al. Proteomics of Medicago truncatula seed development establishes the time frame of diverse metabolic processes related to reserve accumulation. Plant Physiol. (2003) 133:664–682.
Gonzales MD, et al. The Legume Information System (LIS): an integrated information resource for comparative legume biology. Nucleic Acids Res. (2005) 33:D660–D665.
Huhman DV, Sumner LW. Metabolic profiling of saponins in Medicago sativa and Medicago truncatula using HPLC coupled to an electrospray ion-trap mass spectrometer. Phytochemistry (2002) 59:347–360.[CrossRef][Web of Science][Medline]
Huhman DV, et al. Quantification of saponins in aerial and subterranean tissues of Medicago truncatula. J. Agric Food Chem. (2005) 53:1914–1920.[CrossRef][Web of Science][Medline]
Imin N, et al. Proteome reference maps of Medicago truncatula embryogenic cell cultures generated from single protoplasts. Proteomics (2004) 4:1883–1896.[CrossRef][Web of Science][Medline]
International Rice Genome Sequencing Project. The map-based sequence of the rice genome. Nature (2005) 436:793–800.[CrossRef][Medline]
Karp PD, et al. The Pathway Tools software. Bioinformatic (2002) 18:S225–S232.[Abstract]
Keseler IM, et al. EcoCyc: a comprehensive database resource for Escherichia coli. Nucleic Acids Res. (2005) 33:D334–D337.
Krieger CJ, et al. MetaCyc: a multiorganism database of metabolic pathways and enzymes. Nucleic Acids Res. (2004) 32:D438–D442.
Kuster H, et al. Construction and validation of cDNA-based Mt6k-RIT macro- and microarrays to explore root endosymbioses in the model legume Medicago truncatula. J. Biotechnol. (2004) 108:95–113.[CrossRef][Web of Science][Medline]
Lamblin A-FJ, et al. MtDB: a database for personalized data mining of the model legume Medicago truncatula transcriptome. Nucleic Acids Res. (2003) 31:196–201.
Lange B, Ghassemian M. Comprehensive post-genomic data analysis approaches integrating biochemical pathway maps. Phytochemistry (2005) 66:413–451.[CrossRef][Web of Science][Medline]
Lei Z, et al. A two-dimensional electrophoresis proteomic reference map and systematic identification of 1367 proteins from a cell suspension culture of the model legume Medicago truncatula. Mol. Cell Proteomics (2005) 4:1812–1825.
Liu C, et al. Regiospecific hydroxylation of isoflavones by cytochrome p450 81E enzymes from Medicago truncatula. Plant J. (2003) 36:471–484.[CrossRef][Web of Science][Medline]
Mathesius U, et al. Establishment of a root proteome reference map for the model legume Medicago truncatula using the expressed sequence tag database for peptide mass fingerprinting. Proteomics (2001) 1:1424–1440.[CrossRef][Web of Science][Medline]
Mueller LA, et al. AraCyc: a biochemical pathway database for arabidopsis. Plant Physiol. (2003) 132:453–460.
Paley SM, Karp PD. Evaluation of computational metabolic-pathway predictions for Helicobacter pylori. Bioinformatics (2002) 18:715–724.
Sumner LW, et al. Plant metabolomics: large-scale phytochemistry in the functional genomics era. Phytochemistry (2003) 62:817–836.[CrossRef][Web of Science][Medline]
Suzuki H, et al. Methyl jasmonate and yeast elicitor induce differential transcriptional and metabolic re-programming in cell suspension cultures of the model legume Medicago truncatula. Planta (2005) 220:696–707.[CrossRef][Web of Science][Medline]
Thimm O, et al. MapMan: a user-driven tool to display genomics data sets onto diagrams of metabolic pathways and other biological processes. Plant J. (2004) 37:914–939.[CrossRef][Web of Science][Medline]
Town CD. Annotating the genome of Medicago truncatula. Curr. Opin. Plant Biol. (2006) 9:122–127.[CrossRef][Web of Science][Medline]
Tuskan GA, et al. The genome of black cottonwood, Populus trichocarpa (Torr. & Gray). Science (2006) 313:1596–1604.
Urbanczyk-Wochniak E, et al. Conversion of MapMan to allow the analysis of transcript data from Solanaceous species: effects of genetic and environmental alterations in energy metabolism in the leaf. Plant Mol. Biol. (2006) 60:773–792.[CrossRef][Web of Science][Medline]
Usadel B, et al. Extension of the visualization tool MapMan to allow statistical analysis of arrays, display of corresponding genes, and comparison with known responses. Plant Physiol. (2005) 138:1195–1204.
Watson BS, et al. Mapping the proteome of Barrel Medic (Medicago truncatula). Plant Physiol (2003) 131:1104–1123.
Weckwerth W. Metabolomics in systems biology. Annu. Rev. Plant Biol. (2003) 54:669–689.[CrossRef][Medline]
Young ND, et al. Sequencing the genespaces of Medicago truncatula and Lotus japonicus. Plant Physiol. (2005) 137:1174–1181.
Zhang P, et al. MetaCyc and AraCyc. Metabolic pathway databases for plant research. Plant Physiol. (2005) 138:27–37.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
