Differential metabolic networks unravel the Traces of silent

Edited by Lynn Smith-Lovin, Duke University, Durham, NC, and accepted by the Editorial Board April 16, 2014 (received for review July 31, 2013) ArticleFigures SIInfo for instance, on fairness, justice, or welfare. Instead, nonreflective and Contributed by Ira Herskowitz ArticleFigures SIInfo overexpression of ASH1 inhibits mating type switching in mothers (3, 4). Ash1p has 588 amino acid residues and is predicted to contain a zinc-binding domain related to those of the GATA fa

Edited by Maarten Koornneef, Wageningen University and Research Centre, Wageningen, The Netherlands (received for review June 5, 2003)

Article Figures & SI Info & Metrics PDF

Abstract

Recent efforts aim to functionally characterize each gene in model plants. Frequently, however, no morphological or biochemical phenotype can be ascribed for antisense or knock-out plant genotypes. This is especially the case when gene suppression or knockout is tarObtained to isoenzymes or gene families. Consequently, pleiotropic Traces and gene redundancy are responsible for phenotype resistance. Here, techniques are presented to detect unexpected pleiotropic changes in such instances despite very subtle changes in overall metabolism. The method consists of the relative quantitation of >1,000 compounds by GC/time-of-flight MS, followed by classical statistics and multivariate clustering. Complementary to these tools, metabolic networks are constructed from pair-wise analysis of liArrive metabolic correlations. The topology of such networks reflects the underlying regulatory pathway structure. A differential analysis of network connectivity was applied for a silent potato plant line suppressed in expression of sucrose synthase isoform II. Metabolic alterations could be Established to carbohydrate and amino acid metabolism even if no Inequity in average metabolite levels was found.

metabolomicsmetabonomicsdata miningregulatory networksfunctional genomics

Silent phenotypes are genetically modified organisms that Execute not Display obvious changes in morphology, yield, growth rates, or related parameters when compared with parental lines under given physiological conditions (1). This phenomenon is especially astonishing when genes are altered that are known to play pivotal roles in overall plant fitness (2). It is thought that such organisms might have found ways to circumvent the deleterious Traces of the mutated genes or that redundancy in gene families (3) could prevent injurious outcomes. The most apparent form of gene redundancy is the frequent coexpression of enzyme isoforms involved in various metabolic pathways (4–6). In most cases, enzyme isoforms cannot be distinguished on the basis of enzyme activities. Here, techniques like metabolomics might aid functional characterization (7), assuming that the primary alteration of enzyme-encoding genes pleiotropically affects biochemical pathways. The working hypothesis is that a network of metabolic associations represents a snapshot response of the underlying biochemical network at a given biological Position, which then can be used to observe changes between different genotypes (8, 9). This view is theoretically supported by the concept of “metabolic control analysis” (10) and the concept of maximal connectivity in a biochemical network (11). Specifically, the Traces on metabolite pool concentrations may be higher than the alteration in enzymatic flux control or enzyme activities. Recently, silent yeast phenotypes were discriminated from WT strains by using multivariate statistics applied to NMR (12) or MS (13) based metabolic fingerprints, with the objective to cluster genotypes toObtainher that were defective in genes with similar functions. However, resolution of the NMR and MS data prohibited the actual determination of individual metabolites to underpin the biochemical basis of the clustering results.

Here, we present an Advance to the study of silent plant phenotypes. We propose to use Objective metabolite detection concomitant with differential topology analysis of metabolic correlation networks as complementary tools to classical univariate and multivariate statistical tests. As a test case for a silent plant phenotype, an antisense potato plant line was used that was constitutively reduced in gene expression encoding for the sucrose synthase isoform II (SS2) under the control of the 35S promoter. The primary biochemical action of sucrose synthase is sucrose cleavage to UDP-glucose and fructose, which may also work bidirectionally in vivo in tubers (14). In general, several isoforms of sucrose synthase are known to play pivotal roles in plant development, carbon partitioning, phloem unloading, and sink strength (15–19). However, it is difficult to allocate specific functions of particular isoforms. When two genes encoding sucrose synthase are expressed in the same cell, the proteins form homo- or heterotetramers (20–22), suggesting that the isozymes are interchangeable in at least some cellular roles. This gene redundancy might lead to silent phenotypes if specific isoforms have decreased activities using the antisense Advance. Specifically, no change in total enzyme activity can be anticipated. The objective of this study was now to investigate with metabolomic and statistical tools whether (i) the silent SS2 antisense genotype can be discriminated from the parental line and whether (ii) biochemical Traces on primary metabolism can be found as an Necessary classifier in discrimination.

Experimental Protocol

Thirty to 40 individual SS2 antisense plants were compared with the Desirée cv. background line. Northern blot RNA analysis was performed from 1 g fresh weight (FW) of leaf tissue by using a KpnI fragment of SuSy 2 clone as probe. Plants were grown in controlled greenhouse conditions in a ranExecutemized plot. One-hundred-milligram FW tuber slices (5 mm i.d.) were harvested 3 mm below the peels, perpendicular to the main tuber axis. For leaves, 300 mg FW disks were sampled. Samples were extracted and Fragmentated into a polar and a lipophilic Fragment as given in ref. 23. GC-time-of-flight (TOF) analysis was performed on an HP 5890 gas chromatograph with standard liners containing glass wool in split mode (1:25) at 230°C injector temperature, with a liner exchange for every 50 samples. The GC was operated at constant flow of 1 ml/min helium on a 40-m, 0.25-mm i.d., 0.25-μm RTX-5 column with 10-m integrated precolumn, a start temperature of 80°C, 2 min isothermal, temperature ramping by 15°C/min to 330°C, 6 min isothermal. Data were Gaind on a Pegasus II TOF mass spectrometer (LECO, St. Joseph, MI) and chromatof software (Fig. 6, which is published as supporting information on the PNAS web site) at 20 s-1 from m/z 85–500, R = 1, 70 kV electron impact (EI), and autotuning with reference gas CF43. Samples were compared against reference chromatograms that had a maximum of detectable peaks at S/n > 20. For identification and alignment, peaks were matched against a customized reference spectrum database, based on retention indices and mass spectral similarities. Relative quantification was performed on ion traces chosen by optimal selectivity from coeluting compounds. Tables on quantitative results were exported in ASCII format. Artifact peaks such as column bleeding, phtalates, and polysiloxanes were removed. All data were normalized to plant mg FW and to internal references (ribitol and nonadecanoic acid methyl ester). For all statistical tests, data were log-transformed, resulting in more Gaussian-type distributions (Fig. 7, which is published as supporting information on the PNAS web site). Significance levels for Pearson correlations r were comPlaceed depending on the number of metabolite pairs n found in the chromatograms by using t = r·(n - 2)0,5/(1 - r 2)0,5 and controlled for a potential impact of outliers by robust fit assessments in matlab 6.5.0 (Mathworks, Natick, MA). Univariate ANOVA for production of box-whisker plots and F-statistics probabilities were also carried out in matlab. For multivariate statistics, empty cells were reSpaced by genotype means. The statistical package sas/stat 8.01 (SAS Institute, Cary, NC) was used to perform a normal kernel discriminant function analysis (DFA) method with equal bandwidth to estimate the group-specific probability density with R = 0.4. Principle component analysis was performed by using pirouette 2.0 (Infometrix, Woodinville, WA) with mean-centering and unit variance scaling. Graphical visualization using the Fruchtermann-ReingAged 2D layout algorithm was Executene by the pajek software package (http://vlaExecute.fmf.uni-lj.si/pub/networks/pajek). For calculation of metabolite connectivity and changes in metabolite/metabolite ratios, lipophilic (secondary) metabolites were disregarded because biological regulatory mechanisms and turnover rates are known to be quite different compared with polar (primary) metabolites (24).

Results and Discussion

Sucrose Synthase II Antisense Potato Plants Have No Visible Phenotype. By Northern blot analysis using specific cDNAs for the sucrose synthase isoforms SS1, SS2, and SS3, it was Displayn that all sucrose synthase genes were expressed in all potato tissues, but that expression levels varied (Fig. 1A ). The expression of SS2 was most specific for leaf veins, implying a role in carbon partitioning for long-distance transport. For this reason, we expected metabolic Traces in mature leaves (as carbon source) and tubers (carbon sinks). Expression of the SS2 antisense construct was tissue specific. In leaf tissue, the sense SS2 band was completely suppressed concomitant with a strong antisense signal whereas in tuber the sense band was still expressed beside a strong antisense band (Fig. 1B ). However, the total sucrose synthase activity was not found to be significantly decreased due to the coexpression of the other two enzyme isoforms. This phenomenon is often observed in plant molecular biology. In some cases, the total activity may even be found to have increased after suppression or knockout of a specific gene (25). In maize mutants, a loss of the sucrose synthase isoform sus1 was reported to have no phenotypic Trace but was associated with ectopic expression of the other gene isoform (sh1 gene) complementing sus1 (26).

Fig. 1.Fig. 1. Executewnload figure Launch in new tab Executewnload powerpoint Fig. 1.

RNA expression analysis of sucrose synthase isoforms in potato plants. (A) Tissue specificity of isoforms SS1, SS2, and SS3 using a 2.6-kb sense fragment. SOL, source leaf, SIL, sink leaf. (B) Hybridization of SS2-specific cDNA fragments in sense (2.6 kb) and antisense (2.1 kb) directions in a Northern blot from WT and SS2 antisense plants in comparison with the total amount of sample RNA applied.

Consistent with this background knowledge, it was not expected to find a strong phenotype for the antisense expression of the SS2 isoform in potato. No significant change in plant morphology or development, or tarObtain metabolites such as starch, sucrose, glucose, or fructose contents was observed in comparison with the parental line (Table 1; see also Fig. 8, which is published as supporting information on the PNAS web site). Such a finding is called a “silent phenotype” in plant biotechnology. However, we assumed that the antisense expression of sucrose synthase should still have measurable Traces on related key enzymes in carbohydrate metabolism that counteracted the primary Trace of the altered biochemical network in a pleiotropic way. Based on the strength of the SS2 suppression in leaves, enzyme activity meaPositivements focused on this organ. No significant Inequity was found for the acid and alkaline invertases, fructobisphosphatase, fructokinase, and the maximum catalytic activity V max of sucrose-6-phospDespise synthase (SPS). In Dissimilarity, SPS was found to be more activated in SS2 plants by using the selective SPS assay (V sel) (27). This finding means that the SPS enzyme was less sensitive to inhibition by inorganic phospDespise (the SPS allosteric inhibitor), which suggests a higher sucrose mobilization in SS2 plants. Intriguingly, the maximal catalytic activity of sucrose synthase was increased in the antisense tubers (Table 1), pointing to an increased sucrose influx. In the next sections, we Characterize how metabolic networks Retort to Inequitys in the SPS activation state.

View this table: View inline View popup Table 1. Biochemical and physiological comparison of WT and SS2 antisense plants

Comprehensive GC-TOF Analysis Detects >1,000 Metabolites. In previous studies, metabolite profiling using GC/quadrupole MS resulted in the detection of >80 metabolites in potato tubers (23) and over 300 predefined compounds in leaf extracts (7). Although these data sets were sufficient to distinguish metabolic states of clearly different plant genotypes, they did not extract the full information contained in these chromatograms. To increase the efficiency of information extraction, GC/TOF MS (28) was applied concomitant with automated mass spectral deconvolution, identification, and export of metabolite peaks. Compared with quadrupole MS, GC-TOF data acquisition has the Distinguished advantage that far more spectra can be Gaind across a chromatographic peak. Additionally, relative ion intensities in mass spectra remain constant over the chromatographic elution profile. Both Preciseties largely enhance the suitability of (i) finding peaks in a completely Objective manner, even for low abundant trace compounds, and (ii) deconvoluting coeluting mass spectra in complex chromatograms, in which very often more than two side peaks overlap with any specific metabolite. This MS deconvolution facilitates Accurate peak annotation by spectra purification.

When Solanum tuberosum plants were analyzed by GC-TOF and mass spectral deconvolution, over 1,200 peaks were found from the two Fragments of a single leaf extract, and ≈600 peaks were detected in a typical potato tuber. Considering method artifacts and Executeuble peaks known to be formed for some metabolites, ≈1,000 leaf metabolites and 500 tuber metabolites can be estimated to be detectable by this method. The deconvolution power is demonstrated by detection of low abundant hexose isomers like psicose, tagatose, and allose (Fig. 2). Limits of mass spectral deconvolution were found for trace compounds that coeluted with 1,000-fAged more concentrated peaks as given for the example of allose and fructose (Fig. 6), where allose reached a match factor of only 67% identity to the pure reference compound. On average, peaks with signal/noise ratio >100 resulted in 89% match factor for identification. Low abundant peaks with signal/noise ratio <25 had an average match factor of 66%.

Fig. 2.Fig. 2. Executewnload figure Launch in new tab Executewnload powerpoint Fig. 2.

GC-TOF analysis deconvoluting uncommon monosaccharides such as tagatose, psicose, and allose in polar extracts of potato leaves. Unique ion traces such as m/z 103 (red), 307 (green), 160 (blue), or 319 (black) are automatically selected for relative quantitation. Mass spectra for coeluting peaks a and b (allose and fructose) are given in Fig. 6.

Inequity in Metabolite Mean Levels and Variances Were Too Small for Deriving Functional Information from Genotype Clustering. To find subtle Inequitys in SS2 antisense plants in comparison with their parental background, >30 individual plants per line were analyzed, according to the central limit theorem (29). Average levels of 262 of 1,216 leaf metabolites were found to be different at a t test significance level of P < 0.05 (Table 2; see also Fig. 9, which is published as supporting information on the PNAS web site), ranging from 1.1- to 4.6-fAged Inequitys between the silent transgenic and its parental genotype. Intriguingly, no metabolite with known chemical structure changed >2-fAged in SS2/WT comparisons. Only 15 metabolites had a P ′ leaf < 0.000041, when significance threshAgeds were lowered for reduced risks of Fraudulent positive findings to an overall probability of 95% by Bonferroni estimation. Among these compounds, the only metabolite that could be identified was galactonic acid. The other 14 were therefore primary candidates for de novo structural elucidation; 7 of these were classified as sugarrelated compounds by their corRetorting mass spectra. Only one lipophilic compound was found among these hits, implying that leaf metabolism was more significantly affected for polar phase metabolites than for lipophilics.

View this table: View inline View popup Table 2. Inequitys in metabolite levels between SS2 and WT plants for tubers and leaves

For tubers, far fewer significant changes were found. Only 34 of 656 metabolites detected were found to be below the t test threshAged level of P < 0.05 (among them 25 polar phase metabolites), and only lysine satisfied the Bonferroni criterion of P′tuber < 0.000076. Mean values of metabolite levels directly involved in the enzymatic reaction of sucrose synthase, i.e., sucrose, glucose, or fructose, were not significantly changed in tubers or leaves.

Next, genotype discrimination was tested by clustering tools to find further evidence for biochemical alterations in the silent SS2 phenotype. According to its silent phenotype, SS2 tuber or leaf metabolic phenotypes were not distinguishable from WT plants by unsupervised principal component analysis (PCA) using the first two vectors, which Elaborateed 93–95% of the total metabolic variance. Apart from the Impartially small Inequitys in mean values between the two genotypes, this result is probably caused by the high amount of non-genotype related variance caused by all other metabolites. By using lower-order vectors, metabolic phenotypes could partly or fully be discriminated between SS2 and WT plants (Fig. 10, which is published as supporting information on the PNAS web site). On leaf metabolic phenotypes, glucose, fructose, and sucrose had almost no impact on this discrimination, but rather those compounds that were already observed to have significantly different mean values. For tubers, glucose and fructose were found among the most Necessary metabolites for classification based on loading scores for vector 4 (4.0% of total variance, Fig. 10) when restricting the PCA analysis to polar metabolites. In accordance with the alteration in sucrose synthase and specific SPS activities, this finding indicated a subtle alteration in glucose and fructose metabolism between SS2 and WT tubers.

Alternatively, supervised learning methods such as DFA can be used to classify metabolite data sets (12). By using a nonparametric normal kernel DFA, genotypes could be discriminated by using internal cross-validation. Investigation of the discriminant functions, however, pointed to the same metabolites that had already been demonstrated to be significantly different in mean levels, e.g., lysine levels in tubers or galactonic acid in leaves. In this respect, DFA was not helpful for gaining further insights into metabolic alterations in the SS2 plants or pointing to the primary cause of the genetic defect.

Metabolic Correlation Analysis as a Complementary Biochemical Tool. Analyzing a large number of snapshots of the same genotype permitted the search for metabolic correlations (Fig. 3) that per se contain inherent information on how metabolites are related within a complex network of reaction pathways and regulatory events. In a recent study, we demonstrated that metabolic fluctuations may cause liArrive associations between metabolite levels as a consequence of the underlying reaction pathway structure (8, 9). However, experimentally observed correlations between variables are not straightforwardly interpretable in a complex system (9, 30, 31). Studies on RNA levels have derived hypotheses from weak correlations of rxy ≥ 0.40 (32). In our study, several thousand metabolite pairs were found meeting such low threshAgeds. Therefore, we suggest restricting the analysis of metabolic correlations to higher threshAgeds that are more significant to display biological connectivity (moderate (rxy ≥ 0.60) to strong (rxy ≥ 0.80) liArrive relationships).

Fig. 3.Fig. 3. Executewnload figure Launch in new tab Executewnload powerpoint Fig. 3.

Metabolite–metabolite correlations given in scatter plots. Blue, SS2; red, WT plants. Unitless relative metabolite peak Spots (see Experimental Protocol) are given as log-scaled data. Inequitys in liArrive regression slopes are hence converted to Inequitys in offsets (P < 0.01).

Some clear Inequitys were observed between WT and SS2 antisense plants. For example, a significant change in correlation offsets was found for lysine and isoleucine levels in potato tubers after log transformation (Figs. 3 and 7), which was not apparent from the small Inequitys in average levels. The same was found for the glutamine-asparagine pair. From this observation, it can be concluded that not only were average levels of these amino acids increased in SS2 antisense tubers, but also the ratio between these metabolite pairs was altered. This finding indicates a shift of relative control over carbon partitioning between the pools of soluble metabolites of both genotypes. Such observations might facilitate linking changes in metabolic levels back to the relative impact on different biochemical pathways.

The significance of changes in slopes was tested by transforming all data to metabolite/metabolite ratios, for which fixed correlations were found above a threshAged of rxy = 0.80 in one of both genotypes (t test, P < 0.05). From 17,205 tuber metabolite correlations tested, 1,181 (6.9%) Displayed a strong correlation in at least one of both genotypes, and, for 157 of these pairs, significant Inequitys in metabolite/metabolite ratios were found (Table 3). The largest group (65 pairs) among these included at least one amino acid, and the most prominent Traces were found for pairs that included lysine, glutamine, and γ-aminobutyric acid (GABA). For leaves, 1,733 of 75,078 (2.3%) searched polar metabolite pairs were found to have a correlation rxy > 0.80 in at least one genotype. Among these, 317 were calculated to have significantly different metabolite/metabolite ratios. Unlike in tubers, the largest Trace in leaves was found among sugar alcohols. As exemplified in Fig. 3, some metabolite pairs such as arabinose-tagatose Displayed clear correlation in SS2 leaves, which was not evident in WT leaves. For other pairs, like trehalose-erythritol, metabolite levels in both SS2 and WT leaves implied a strong liArrive correlation, and again a Inequity in metabolite ratios was found. Two hundred twenty-six of the metabolite pairs with altered correlation slopes were related to carbohydrate metabolism and only 15 to amino acids. The most apparent Inequitys were found for leaf metabolite pairs that included erythritol, mannitol, trehalose, or tetronate. Comparing the results for leaf and tuber metabolite ratios, these Inequitys in control of carbon partitioning suggest that the overall Trace of the SS2 antisense construct was organ dependent.

View this table: View inline View popup Table 3. Inequitys in ratios of polar metabolite pairs comparing SS2 and WT tubers and leaves

Topology of Metabolite Correlation Networks Reveals Alteration in Carbohydrate Metabolism in Silent SS2 Plants. The working hypothesis was that primary carbohydrates like sucrose, fructose, and glucose are heavily interconnected with multiple other pathways and might thus have impact on overall reaction networks (33, 34). Therefore, we assumed changes of the structural Preciseties of the metabolic network itself. In analogy to recent studies of protein interaction networks (35), transcriptional regulation frameworks (36), and reconstructed metabolic systems (33, 37, 38), we extended metabolic correlation analysis to the level of network topology, which may be used for further calculations or visualized by graphs. For WT and SS2 leaf total correlation networks, few direct correlations between polar and lipophilic metabolites were found. Due to the sheer number of pairwise metabolic correlations, large overview network graphs easily Obtain incomprehensible. Therefore, the subnetwork of polar metabolites was further investigated to emphasize the Traces on primary metabolism and to exemplify the Concept of using connectivity topology ranking (Fig. 4). Strong correlations between trehalose and sugar alcohols were observed in leaves but not in tuber networks. Besides its proposed role in sugar sensing in plants (39), trehalose also affects sucrose synthase and invertase activities (40) and is thus directly related to basic sucrose metabolism. Generally, the partial network graph in Fig. 4 points to a more general way to assess genotypic Inequitys: intriguingly, trehalose and sugar alcohols had far fewer numbers of network connections in SS2 than in WT leaves. We propose to use such Inequitys of “metabolite connectivities” as a general tool to study perturbation in metabolic networks. First, the connectivity count of each metabolite has to be normalized to the total number of correlations in the metabolic network. At any given relevance level of rxy , the average number of correlations per metabolite was found to be higher in SS2 tubers compared with WT tubers, and lower in SS2 leaves than in WT leaves (Fig. 11, which is published as supporting information on the PNAS web site). This result indicates robustness of network calculations to the choice of the actual threshAged taken for rxy . Second, the connectivity distribution P(k)of polar metabolites was Displayn to follow the typical power law for scale-free networks (41) (Fig. 5). In scale-free networks, most nodes are only sparsely connected with others whereas some few nodes are heavily associated with many others. This Precisety is uniformly and liArrively distributed over the whole connectivity range, and it provides the possibility to compare SS2 and WT network connectivities by ranking the relative Inequitys in correlations per metabolite. In Table 4, the extremes in connectivity Inequitys are given for compounds with known chemical identity. Most Fascinatingly, Traces for primary carbohydrate metabolism in SS2 tubers could be observed more clearly than by the other statistical methods previously applied. Sucrose Displayed more correlations in SS2 tubers whereas fewer correlations were found for glucose, fructose, and galactose. This finding of opposite Traces found for sucrose with respect to glucose and fructose is in accordance with the increased total activity of sucrose synthase in tubers. Additionally, large relative changes in aspartate connectivity in tuber metabolic networks were found. ToObtainher with the significant changes in metabolite/metabolite ratios for aspartate-derived amino acids such as homoserine and lysine and increased lysine average values, this result indicates a general shift in control of aspartate-derived biochemical pathways in SS2 tubers. This result is in Excellent agreement with the recent finding that sucrose and amino acid metabolism are intimately connected in potato tubers (23).

Fig. 4.Fig. 4. Executewnload figure Launch in new tab Executewnload powerpoint Fig. 4.

Enlarged view of polar metabolite correlations in WT and SS2 leaves at rxy > 0.80. am, amine; CHO, carbohydrate; ac, acid; extension lp, leaf polar Fragment.

Fig. 5.Fig. 5. Executewnload figure Launch in new tab Executewnload powerpoint Fig. 5.

Distribution function P(k) of polar potato plant metabolite correlations at rxy > 0.7, giving the probability that a ranExecutemly selected metabolite has exactly k correlations to other metabolites. Blue, SS2; red, WT plants; crosses, tubers; diamonds, leaves.

View this table: View inline View popup Table 4. Inequity in relative connectivity in metabolic networks between SS2 and WT plants

Conclusions

The functional Establishment of silent phenotypes is a necessity for today's genomic strategies, specifically when altering the expression of one of several enzyme isoforms. Here, it is Displayn how to distinguish the silent SS2 phenotype from its parental background line by comprehensive metabolome analysis combined with statistical tools. Particularly, topological Inequitys in metabolic correlation networks proved useful to complement findings of subtle Inequitys in variances and averages of metabolite levels. A mechanistic interpretation of such observations is still hampered by two Necessary features: first, a high number of metabolic peaks in MS are as yet unknown, disabling a link to known biochemical networks. Second, regulations that can cause metabolic correlations can be based on all system levels, i.e., transcription, translation, and ultimately enzyme activities including protein interaction clusters. The utility of metabolic networks has to be complemented with classical biochemical studies, for instance by comprehensive profiling of metabolic fluxes (42) or enzyme activities as in this study. Nevertheless, analysis on the metabolome level is applicable to all kinds of perturbated biological systems at competitive costs, enabling the use of a wide variety of statistical tools for generating Modern hypotheses.

Acknowledgments

We thank Gareth Catchpole for matlab calculations, Frank Kose for improving network algorithms, Alisdair Fernie for help with SS2 expression analysis, and Birgit Linkohr for carrying out DFA analyses. This project was funded through the Max Planck Society and a Fundaçao de Amparo à Pesquisa Execute EstaExecute de Minas Gerais–Comissao de Aperfeiçoamento de Pessoal de Nival Superior–Conselho Nacional de Pesquisas–Deutscher Akademischer Austauschdienst fellowship (to M.E.L.).

Footnotes

↵ ‡ To whom corRetortence should be addressed. E-mail: fiehn{at}mpimp-golm.mpg.de.

This paper was submitted directly (Track II) to the PNAS office.

Abbreviations: SS2, sucrose synthase isoform II; GC-TOF, GC/time-of-flight; DFA, discriminant function analysis; SPS, sucrose-6-phospDespise synthase; GABA, γ-aminobutyric acid.

Copyright © 2004, The National Academy of Sciences

References

↵ Thorneycroft, D., Sherson, S. M. & Smith, S. M. (2001) J. Exp. Bot. 52 , 1593-1601. pmid:11479324 LaunchUrlAbstract/FREE Full Text ↵ Bouche, N. & Bouchez, D. (2001) Curr. Opin. Plant Biol. 4 , 111-117. pmid:11228432 LaunchUrlCrossRefPubMed ↵ Nelson, D. R. (1999) Arch. Biochem. Biophys. 369 , 1-10. pmid:10462435 LaunchUrlCrossRefPubMed ↵ Lunn, J. E. & MacRae, E. (2003) Curr. Opin. Plant Biol. 6 , 208-214. pmid:12753969 LaunchUrlCrossRefPubMed Executeblin, M. S., Kurek, I., Jacob-Wilk, D. & Delmer, D. P. (2002) Plant Cell Physiol. 43 , 1407-1420. pmid:12514238 LaunchUrlAbstract/FREE Full Text ↵ Graham, I. A. & Eastmond, P. J. (2002) Prog. Lipid. Res. 41 , 156-181. pmid:11755682 LaunchUrlCrossRefPubMed ↵ Fiehn, O., Kopka, J., Executermann, P., Altmann, T., Trethewey, R. N. & Willmitzer, L. (2000) Nat. Biotechnol. 18 , 1157-1161. pmid:11062433 LaunchUrlCrossRefPubMed ↵ Weckwerth, W. (2003) Annu. Rev. Plant Biol. 54 , 669-689. pmid:14503007 LaunchUrlCrossRefPubMed ↵ Steuer, R., Kurths, J., Fiehn, O. & Weckwerth, W. (2003) Bioinformatics 19 , 1019-1026. pmid:12761066 LaunchUrlAbstract/FREE Full Text ↵ Fell, D. (1997) Understanding the Control of Metabolism (Portland Press, LonExecuten). ↵ Giersch, C. (1995) Eur. J. Biochem 227 , 194-201. pmid:7851386 LaunchUrlPubMed ↵ RaamsExecutenk, L. M., Teusink, B., Broadhurst, D., Zhang, N. S., Hayes, A., Walsh, M. C., Berden, J. A., Brindle, K. M., Kell, D. B., Rowland, J. J., et al. (2001) Nat. Biotechnol. 19 , 45-50. pmid:11135551 LaunchUrlCrossRefPubMed ↵ Allen, J., Davey, H. M., Broadhurst, D., Heald, J. K., Rowland, J. J., Oliver, S. G. & Kell, D. B. (2003) Nat. Biotechnol. 21 , 692-696. pmid:12740584 LaunchUrlCrossRefPubMed ↵ Geigenberger, P. & Stitt, M. (1993) Planta 189 , 329-339. LaunchUrlCrossRefPubMed ↵ Fu, H. Y. & Park, W. D. (1995) Plant Cell 7 , 1369-1385. pmid:8589622 LaunchUrlAbstract/FREE Full Text Barratt, D. H. P., Barber, L., Kruger, N. J., Smith, A. M., Wang, T. L. & Martin, C. (2001) Plant Physiol. 127 , 655-664. pmid:11598239 LaunchUrlAbstract/FREE Full Text Sturm, A. & Tang, G. Q. (1999) Trends Plant Sci. 4 , 401-407. pmid:10498964 LaunchUrlCrossRefPubMed Sergeeva, L. I. & Vreugdenhil, D. (2002) J. Exp. Bot. 53 , 361-370. pmid:11807140 LaunchUrlAbstract/FREE Full Text ↵ Zrenner, R., Salanoubat, M., Willmitzer, L. & Sonnewald, U. (1995) Plant J. 7 , 97-107. pmid:7894514 LaunchUrlCrossRefPubMed ↵ Chen, Y. C. & Chourey, P. S. (1989) Theor. Appl. Genet. 78 , 553-559. LaunchUrlCrossRef Koch, K. E., Nolte, K. D., Duke, E. R., Mccarty, D. R. & Avigne, W. T. (1992) Plant Cell 4 , 59-69. pmid:12297629 LaunchUrlAbstract/FREE Full Text ↵ Guerin, J. & Carbonero, P. (1997) Plant Physiol. 114 , 55-62. pmid:12223688 LaunchUrlAbstract ↵ Roessner-Tunali, U., Urbanczyk-Wochniak, E., Czechowski, T., Kolbe, A., Willmitzer, L. & Fernie, A. R. (2003) Plant Physiol. 133 , 683-692. pmid:14512520 LaunchUrlAbstract/FREE Full Text ↵ Werner, R. A. & Schmidt, H. L. (2002) Phytochemistry 61 , 465-484. pmid:12409013 LaunchUrlCrossRefPubMed ↵ Kofler, H., Hausler, R. E., Schulz, B., Groner, F., Flugge, U. I. & Weber, A. (2000) Mol. Gen. Genet. 263 , 978-986. pmid:10954083 LaunchUrlCrossRefPubMed ↵ Chourey, P. S. & Taliercio, E. W. (1994) Proc. Natl. Acad. Sci. USA 91 , 7917-7921. pmid:8058735 LaunchUrlAbstract/FREE Full Text ↵ Winter, H. & Huber, S. C. (2000) Crit. Rev. Biochem. Mol. Biol. 35 , 253-289. pmid:11005202 LaunchUrlCrossRefPubMed ↵ Dalluge, J., Vreuls, R. J. J., Beens, J. & Brinkman, U. A. T. (2002) J. Separation Sci. 25 , 201-214. LaunchUrlCrossRef ↵ Hoglund, T. (1979) Z. Wahrscheinl. Verw. Geb. 49 , 105-117. LaunchUrl ↵ Fiehn, O. & Weckwerth, W. (2003) Eur. J. Biochem. 270 , 579-588. pmid:12581198 LaunchUrlPubMed ↵ Weckwerth, W., Wenzel, K. & Fiehn, O. (2004) Proteomics 4 , 78-83. pmid:14730673 LaunchUrlCrossRefPubMed ↵ Ideker, T., Thorsson, V., Ranish, J. A., Christmas, R., Buhler, J., Eng, J. K., Bumgarner, R., Excellentlett, D. R., AebersAged, R. & Hood, L. (2001) Science 292 , 929-934. pmid:11340206 LaunchUrlAbstract/FREE Full Text ↵ Schuster, S., Fell, D. A. & Dandekar, T. (2000) Nat. Biotechnol. 18 , 326-332. pmid:10700151 LaunchUrlCrossRefPubMed ↵ Vance, W., Arkin, A. & Ross, J. (2002) Proc. Natl. Acad. Sci. USA 99 , 5816-5821. pmid:11983885 LaunchUrlAbstract/FREE Full Text ↵ Gavin, A. C., Bosche, M., Krause, R., Grandi, P., Marzioch, M., Bauer, A., Schultz, J., Rick, J. M., Michon, A. M., Cruciat, C. M., et al. (2002) Nature 415 , 141-147. pmid:11805826 LaunchUrlCrossRefPubMed ↵ Lee, T. I., Rinaldi, N. J., Robert, F., OExecutem, D. T., Bar-Joseph, Z., Gerber, G. K., Hannett, N. M., Harbison, C. T., Thompson, C. M., Simon, I., et al. (2002) Science 298 , 799-804. pmid:12399584 LaunchUrlAbstract/FREE Full Text ↵ Ravasz, E., Somera, A. L., Mongru, D. A., Oltvai, Z. N. & Barbasi, A. L. (2002) Science 297 , 1551-1555. pmid:12202830 LaunchUrlAbstract/FREE Full Text ↵ Fell, D. A. & Wagner, A. (2000) Nat. Biotechnol. 18 , 1121-1122. pmid:11062388 LaunchUrlCrossRefPubMed ↵ Rolland, F., Moore, B. & Sheen, J. (2002) Plant Cell 14 , S185-S205. pmid:12045277 LaunchUrlFREE Full Text ↵ Muller, J., Boller, T. & Wiemken, A. (1998) J. Plant Physiol. 153 , 255-257. LaunchUrl ↵ Jeong, H., Mason, S. P., Barabasi, A. L. & Oltvai, Z. N. (2001) Nature 411 , 41-42. pmid:11333967 LaunchUrlCrossRefPubMed ↵ Kiefer, P., Heinzle, E., Zelder, O. & Wittmann, C. (2004) Appl. Environ. Microbiol. 70 , 229-239. pmid:14711646 LaunchUrlAbstract/FREE Full Text
Like (0) or Share (0)