Skip to main content

Main menu

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
    • Front Matter Portal
    • Journal Club
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses
  • Submit
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Accessibility Statement
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian

User menu

  • Log in
  • My Cart

Search

  • Advanced search
Home
Home
  • Log in
  • My Cart

Advanced Search

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
    • Front Matter Portal
    • Journal Club
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses
  • Submit
Research Article

The metabolic signature related to high plant growth rate in Arabidopsis thaliana

Rhonda C. Meyer, Matthias Steinfath, Jan Lisec, Martina Becher, Hanna Witucka-Wall, Ottó Törjék, Oliver Fiehn, Änne Eckardt, Lothar Willmitzer, Joachim Selbig, and Thomas Altmann
  1. Departments of *Genetics and
  2. ‡Bioinformatics, Institute of Biochemistry and Biology, University of Potsdam, Karl-Liebknecht-Strasse 24–25, 14476 Potsdam, Germany; and
  3. §Max-Planck-Institute of Molecular Plant Physiology, Am Muehlenberg 1, 14476 Golm, Germany

See allHide authors and affiliations

PNAS March 13, 2007 104 (11) 4759-4764; https://doi.org/10.1073/pnas.0609709104
Rhonda C. Meyer
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: rcmeyer@uni-potsdam.de
Matthias Steinfath
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jan Lisec
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Martina Becher
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Hanna Witucka-Wall
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Ottó Törjék
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Oliver Fiehn
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Änne Eckardt
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Lothar Willmitzer
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Joachim Selbig
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Thomas Altmann
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  1. Edited by Maarten Koornneef, Wageningen University and Research Centre, Wageningen, The Netherlands, and approved January 3, 2007 (received for review November 1, 2006)

  • Article
  • Figures & SI
  • Info & Metrics
  • PDF
Loading

Abstract

The decline of available fossil fuel reserves has triggered world-wide efforts to develop alternative energy sources based on plant biomass. Detailed knowledge of the relations of metabolism and biomass accumulation can be expected to yield powerful novel tools to accelerate and enhance energy plant breeding programs. We used metabolic profiling in the model Arabidopsis to study the relation between biomass and metabolic composition using a recombinant inbred line (RIL) population. A highly significant canonical correlation (0.73) was observed, revealing a close link between biomass and a specific combination of metabolites. Dividing the entire data set into training and test sets resulted in a median correlation between predicted and true biomass of 0.58. The demonstrated high predictive power of metabolic composition for biomass features this composite measure as an excellent biomarker and opens new opportunities to enhance plant breeding specifically in the context of renewable resources.

  • biomass
  • canonical correlation
  • metabolic profiling
  • recombinant inbred line population
  • biomarker

Multicellular organisms have to optimize the use of available resources to fit their needs in terms of energy, biosynthetic building blocks, and reserves. Green plants unlike animals produce their own organic compounds. Their ability to grow thus solely depends on their own photosynthetic and metabolic capacity. Biomass accumulation in the vegetative growth phase of a plant can therefore be regarded as the ultimate expression of its metabolic performance. Plants function as integrated systems, in which metabolic and developmental pathways draw on common resource pools and respond to changes in environmental energy and resource supplies (1). The distribution of metabolites between growth, production of defense compounds and storage compounds therefore has to be very tightly regulated. Growth and the concomitant drain of metabolites into cellular components has to be adjusted to the metabolic capacity of the system, i.e., the ability to supply sufficient amounts of organic compounds. This regulation is demonstrated by several observations of growth depression upon reduction of primary metabolism such as sucrose synthesis (2, 3). Growth ceases upon severe starvation caused by an extended dark period and is reinitiated only after a lag period of several hours after relief from the starvation by reillumination (4). Recent observations of the roles of the DELLA proteins in plants indicate that plant growth is limited to a submaximum level to enable plants to cope with unfavorable conditions (5). Thus, growth rate has to be adjusted to the metabolic status of a plant that needs to be translated into an appropriate response. This interaction between metabolism and the growth regulatory mechanisms may operate in two ways: either a high supply of metabolites triggers growth, or growth drains metabolites to a minimum tolerable level upon which growth is restricted. Metabolites may exert control on growth either by acting as substrates for the synthesis of cellular components, that become limiting under conditions of maximum tolerable growth, or by acting as signals that are sensed leading to subsequent changes in growth. Sugars such as glucose and sucrose have been shown to act as metabolic signals and to be involved in the control of plant growth and development (6). Trehalose-6-phosphate has recently been shown to be involved in signaling of the plant sugar status and in control of growth and development (7, 8).

Metabolic profiling is a mass-spectrometry (MS)- or NMR-based technology for an unbiased analysis of the metabolome of a given biological system with a high diagnostic power (9). Thus, in case of e.g., yeast or plants metabolic analysis allows to distinguish between different genotypes, developmental status or environmental conditions (10–12). In the case of humans, metabolomic approaches allow us to predict the response of individuals to drugs opening aspects of personalized drug treatments (13). In addition, single or a small number of metabolites can be extracted from metabolic profiling studies that have the potential to be developed into rapidly accessible biomarkers (14).

Based on the above considerations between the metabolic status of a plant system and growth and the proven high diagnostic power of metabolic profiling approaches, we decided to test whether biomass of a plant is correlated with and can thus be predicted by its metabolic composition. To this end we took advantage of a recombinant inbred line (RIL) population of Arabidopsis thaliana derived from a cross between the Arabidopsis lines Col-0 and C24 (15), which in previous studies showed strong transgressive segregation for biomass (16). RILs represent permanent segregating populations of homozygous lines, which allow to reduce the environmental variance in replicated experiments (17). The extensive biochemical variation in Arabidopsis is largely under genetic control (11). Therefore, the use of such a population for an exploratory analysis of relations between growth and metabolite levels is particularly advantageous (over e.g., using environmental perturbations to modulate growth and metabolism) as it offers the opportunity to identify the genetic determinants of all studied traits in addition to the determination of correlations.

As shown below, when applying multivariate analysis to the combined data sets of biomasses and metabolic profiles, a statistically highly significant correlation between metabolic composition and biomass was obtained. We believe this result to be of high relevance for our basic understanding of plant growth and metabolism and to have obvious implications for breeding of high plant biomass producers, an aspect which in recent years has become of increasing importance regarding renewable resources as energy supply (18, 19). It furthermore provides precedence for the utility of molecular profiling data to extract biomarkers with high predictive power for a complex trait.

Results

Biomass and Metabolite Profile Determination of Col-0/C24 RILs.

The combined analysis of biomass and metabolic profile was performed on a total of 1,144 genotypes. Of these lines, 429 genotypes were derived from a RIL population from the reciprocal crosses Col-0 × C24 (228 lines) and C24 × Col-0 (201 lines) and 715 lines were derived from crosses of the RILs to parents Col-0 and C24. All plants were grown under controlled conditions in six replicated experiments. Plants were harvested 15 days after sowing and used for shoot biomass determination or were pooled and frozen for metabolite profiling by gas-chromatography/mass-spectrometry (GC/MS). The distribution of mean biomass within the population clearly shows transgressive segregation (Fig. 1). We detected no significant differences in biomass (t test, P = 0.238) between the two subpopulations, and therefore treated the RILs as one population in subsequent analyses.

Fig. 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 1.

Distribution of shoot biomass in the recombinant inbred line (RIL) population. Shown is the mean biomass (milligrams per plant) estimated by REML. The arrow indicates the biomass determined for the parental lines C24 (1.265 mg per plant) and Col-0 (1.254 mg per plant). The histogram of the shoot biomass of the RIL crosses to the parents is shown in Inset.

From the metabolic profiling data we took only those metabolites into account, which were detected reproducibly in at least 85% of the samples analyzed. Major groups among these metabolites are organic acids, sugars, sugar phosphates, polyols, amines and amino acids. Concentrations could be determined for a set of 181 compounds, 84 of which were assigned a chemical structure by comparison with a library (20, 21). The remaining compounds were classified into chemical groups by using representative masses.

Canonical Correlation Reveals a Close Link Between Biomass and a Specific Combination of Metabolites.

In a first approach distributions of single metabolites were queried for their predictive power with respect to the biomass distribution by calculating pairwise correlations between all 181 measured metabolite levels and biomass [supporting information (SI) Table 2]. Because a normal distribution cannot be assumed for all variables rank correlation was used as a robust estimation of the correlation coefficient. The highest absolute correlation found was for a carbohydrate, which yielded a value of 0.266. Although the correlation is statistically highly significant (P value of 5.17 × 10−20), it can only explain 7.07% of the variance. Other significantly correlated compounds are ethanolamine (0.238; P = 3.87 × 10−16), fructose-6-phosphate (−0.177; P = 1.65 × 10−9, glutamine (−0.177, P = 1.81 × 10−9), glucose-6-phosphate (−0.175; P = 2.44 × 10−9 and citric acid (−0.175; P = 2.80 × 10−9). Their individual contribution to the explained variance is smaller than 5.64%.

In the second approach we applied multivariate tools to analyze the relationships between the two large groups of metabolite and biomass variables. Canonical correlation analysis (CCA) is a multivariate technique often used in psychological, climate and ecological studies to quantify the associations between two separate data sets measured on the same experimental units (22–25). In contrast to the aforementioned pairwise correlation analysis, CCA yielded a much stronger correlation of 0.73. This value corresponds to 53.29% of variance explained by the linear combination of metabolites, almost 10 times more than explained by any individual metabolite. To test the significance of this result, the biomass vector was permutated 50,000 times. At this point the maximum correlation did not increase significantly with additional permutations. This maximum value is 0.46. The distance between the median of the random correlations and the estimated value amounts to 17 standard deviations (Fig. 2A), which for normal distributions corresponds to a P value of 4.1 × 10−65 strongly suggesting that the model is statistically highly significant.

Fig. 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 2.

Significance (A) and predictive power (B) of the multiplicative model. (A) Histogram of canonical correlations between the metabolite matrix and random permutations of the biomass vector. The vertical line on the right corresponds to the canonical correlation between the actual biomass vector (DW) and the metabolite matrix (X)R. The distance to the median of the random correlations amounts to 17 standard deviations. (B) Prediction of the biomass by the metabolite matrix. Shown is one representative example of 20 repeats in the cross-validation. Size of the training set was 1,086, the 58 data points of the test set are displayed. The straight line represents the exact prediction.

Predictive Power of Metabolic Composition for Biomass.

In a final step we wanted to test the predictive power of metabolite composition for biomass. To this end, we decided to apply the partial least square (PLS) approach, because CCA yields the maximum correlation and thus an upper limit for the true correlation, but is notoriously inferior to other methods, especially PLS, for cross-validation (26). (compare also SI Text). Thus, the metabolite matrix and biomass vector were divided into training and test sets. The PLS coefficients estimated in the training set explaining 90% of the variance of the training data were used to predict the biomass in the test set. This procedure was repeated 20 times. For a size of the training set of 1086 genotypes we obtained a median correlation between predicted and true biomass of 0.58 in the remaining 58 genotypes (representing the test set) confirming a strong predictive power of metabolic composition for biomass (Fig. 2B). To evaluate the significance the same permutation as for CCA was applied. For each of the 500 permutations a cross-validation was performed. The median of the corresponding correlations was −0.001 ± 0.052, thus, using the same assumption as above, we estimate a P value of 3.4 × 10−29.

Metabolites Most Relevant for Biomass Accumulation.

As a next step in our analysis we extracted the metabolites most relevant for biomass accumulation by their correlation to the canonical variate (27). The first 44 metabolites with significant correlations are listed in Table 1 and displayed on biochemical pathways (28) in Fig. 3. Strongly represented are central metabolism derived compounds such as glucose-6-phosphate and fructose-6-phosphate, members of the tricarboxylic acid (TCA) cycle such as succinate, citrate and malate, members of the membrane/phospholipid biosynthesis such as glycerol-3-phosphate, ethanolamine and sinapine, or sucrose. A list of all relevant metabolites is given in SI Table 3.

View this table:
  • View inline
  • View popup
Table 1.

List of top 44 signature metabolites ranked according to the strength of the canonical correlation.

Fig. 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 3.

Representation of the most important metabolites known by structure according to CCA on biochemical pathways. This representation of metabolism indicates all known metabolites we analyzed by using GC/MS that could be annotated in MapMan (28). Red color visualizes metabolites which are high ranked in CCA (positions 1–44), with ranking according to the color-coded scale bar.

Discussion

We took advantage of an Arabidopsis thaliana RIL population for a parallel and integrative analysis of vegetative biomass accumulation and metabolic composition to answer the question whether or not biomass can be described as a function of metabolic composition.

As outlined in the Results section, pairwise correlation analysis of biomass and single metabolites could explain a maximum of 7% of the total variance observed in biomass. These data strongly suggest that there is no single “magic” compound detectable, which could explain the biomass variance in a satisfying way. In contrast, canonical correlation analysis yielded a highly significant (the estimated P value based on permutations is lower than 10−64) canonical correlation of 0.73 (compare Fig. 2A). Furthermore, in cross-validations a median correlation of 0.58 between the predicted and the observed biomass was observed (compare Fig. 2B).

This result demonstrates that a combination of the levels of a large number of metabolites rather than few individual metabolites show a close correlation with growth. It indicates that variation in growth coincides with characteristic combinatorial changes of metabolite levels, whereas individual metabolites may fluctuate largely independently of alterations in growth. To exclude the possibility that the strong correlation between biomass and metabolic composition is simply due to coincidental overlap of quantitative trait loci (QTL) for biomass and metabolites, we performed a QTL analysis on the RIL data set (429 lines) and detected a total of 157 QTL for 84 metabolites and six QTL for biomass (data not shown). Of the latter only two colocate with significantly more metabolite QTL than expected by random, thus making this explanation highly unlikely.

Inspection of the metabolites highly ranked in CCA and thus representing the main drivers of the correlation shows that central metabolism derived metabolites are strongly represented. Of high relevance are the three metabolic intermediates of the hexose phosphate pool, fructose-6-phosphate, glucose-6-phosphate, and glucose-1-phosphate, which link carbon flow from photosynthesis and starch and sucrose metabolism with cell wall formation, the oxidative pentose phosphate pathway (it provides substrates for nucleic acid synthesis and for lignin, polyphenol and amino acid synthesis) and glycolysis. Members of the TCA cycle such as succinate, citrate, and malate are highly ranked. This finding underpins the central importance of this pathway which together with reactions of the glycolysis pathway and the oxidative phosphorylation constitutes a key process delivering carbon skeletons, reduction equivalents, and energy for the vast majority of biochemical pathways. Also highly ranked is sucrose, the major transport form of carbon from source to sink tissue and which is central to the export from the sources and the import to the sinks. It thus represents the interface between carbohydrate production and utilization at the whole plant level Other metabolites such as glycerol-3-phosphate or ethanolamine play a major role in membrane/phospholipid biosynthesis. The anti-oxidant ascorbic acid (vitamin C) has been implicated in cell division (29) and plant growth regulation by means of its role as enzyme cofactor (30). Glutamine as a central metabolite in nitrogen assimilation and the major primary donor of reduced nitrogen is also found amongst the most important metabolites. This observation is contrasted by the fact that nearly all other amino acids analyzed are of rather low contribution based on the CCA. Further highly ranking metabolites can be assigned to general stress metabolites such as sinapine as the major phenylpropanoid in Brassicaceae, ornithine, the polyamines putrescine and spermidine, and trehalose. Thus, a link between the metabolites ranked high in the CCA and biomass accumulation is plausible because central metabolism and stress response are of utmost importance to plant growth, and thus biomass.

Another noteworthy observation is that the canonical variate determined by means of a multiplicative model resulted in closer correlations between the predicted and the observed biomass values than by means of an additive model (data not shown). It indicates that the involved metabolites act synergistically rather than additively which is very plausible as the aforementioned closely interlinked pathways of carbon metabolism are required for different cellular components that all are crucial for growth/biomass formation. The strong reciprocal interrelation between nitrogen and carbon assimilation would also strongly argue for synergistic and not additive effects between key metabolites representing these classes of biochemical compounds as observed in our case. Similar arguments can be made for e.g., ethanolamine synthesized via serine as a major constituent of membranes or sinapic acid as the major phenylpropanoid component in Arabidopsis.

A surprising observation from our data are the occurrence of both positive and negative correlations between metabolites and biomass. The large majority of known metabolites displaying a negative correlation to the biomass vector are the aforementioned intermediates of central metabolic pathways including sucrose, glucose- and fructose-6-phosphate, the TCA cycle members citric acid, succinate or malic acid, as well as the amino acids glutamine and phenylalanine. On the other hand, amongst the positively correlated metabolites are a large fraction of unknown chemical structure as well as some metabolites discussed in stress response such as nicotinic acid (31) or putrescine (32), or the stress metabolite trehalose discussed in connection with drought resistance (33). A negative correlation suggests that pool sizes of these metabolites are reduced to low levels when strong growth occurs. It is conceivable that this process involves mostly metabolites providing the major building blocks for growth such as the central metabolites mentioned. In conclusion, this observation would suggest that growth drives metabolism and not vice versa. This finding would indicate that high growth rates cause a depletion of central metabolite pools rather than growth being enhanced through increased supply of substrates for the synthesis of cellular components. A similar conclusion of metabolism driven by growth has been derived from a study of the relationship between tomato fruit size and metabolites (34). In this scenario, the positively correlated metabolites could play a role in plant defense against abiotic and biotic stress and it is comprehensible that higher concentrations of these metabolites would coincide with better armed plants. For both groups of substances, however, the relation with growth may be nonlinear. On the one hand, the reduction of central metabolite levels below a certain minimum necessary to sustain high flux rates may result in growth limitation and thus a breakdown of a linear negative relationship. Similarly, a positive effect on growth because of elevated stress tolerance may be achieved in a certain range of stress metabolite levels above which no further beneficial or even detrimental effects may occur. As the procedures applied here determine linear correlations, it is not unexpected that no tighter relationships (stronger correlations) were detected. A complementary hypothesis regards metabolites not primarily as chemicals for growth and defense but rather as signals. Under this assumption positively correlated metabolites are positive signals regulating plant growth and the contrary would be true for negatively correlated metabolites. In the context of signal molecules the large number of positively correlated compounds of as yet unknown structure is worth noting and stresses the need for identification of their chemical nature. They might constitute unusual products of metabolic side reactions that are derived from primary metabolites generated for signaling purposes and which can move to sites of perception without further conversion along the major metabolic reactions or transport pathways. Further studies querying some testable predictions from such models (e.g., the presence of receptors/sensors or the elicitation of specific responses in case of signaling metabolites) are needed to validate these models.

Conclusion

Using an Arabidopsis thaliana RIL population and conducting a combined analysis of biomass and metabolite profiles allowed the prediction of biomass as a function of metabolic composition providing a direct proof for the hypothesis that metabolic composition is related to biomass and thus growth. The observations made here further extend this hypothesis toward the notion that major global changes in metabolism are the result of variation in growth rather than vice versa. In addition to fostering our basic understanding, these data are of immediate potential for a number of applied purposes. The possibility to predict biomass on the basis of the metabolic signature of a plant presents a first precedence for the use of metabolite profiles as biomarkers with high predictive power and could potentially revolutionize the selection and thus breeding process for biomass producers such as trees that are cultivated for decades before harvest. Identification of highly productive genotypes already at an early growth stage would result in enormous time and cost-savings. In the light of reduced availability of fossil fuels and increasing reliance on bio-derived energy, the importance of such an opportunity can hardly be overestimated.

Materials and Methods

Creation of Recombinant Inbred Line (RIL) Population.

Two reciprocal sets of RILs were developed from a cross between the two Arabidopsis thaliana accessions C24 and Col-0 as described elsewhere (15). The population consisted of 228 Col-0 × C24 F8 and 201 C24 × Col-0 F8 individual lines.

Plant Cultivation.

The RILs were planted in a split plot design with 54 incomplete blocks and four replicates, repeated six times. Plants were grown in 1:1 mixture of GS 90 soil and vermiculite in 96-well-trays. Six plants of the same line were grown per well. Seeds were germinated in a growth chamber at 6°C for 2 days before transfer to a long-day regime (16 h fluorescent light [120 μmol·m−2·s−1] at 20°C and 60% relative humidity/8 h dark at 18°C and 75% relative humidity). To avoid position effects, trays were rotated around the growth chamber every 2 days.

Shoot Dry Biomass.

Shoot dry biomass was determined 15 days after sowing. Plants from the same well were harvested together and placed in a vacuum oven at 80°C for 48 h. Dry biomass was recorded by using an analysis balance. Mean shoot dry biomass in mg/plant per plant was estimated by using the linear mixed model (35) G + E:E·G + E·GC + E·GC·T where E is experiment, G is genotype, GC is growth chamber, T is tray (REML procedure in Genstat). Biomass in the two subpopulations was compared with a two-sided t test. We detected no significant differences in biomass (P = 0.238) between the two subpopulations, and treated the RILs as one population in subsequent analyses.

Metabolite Data.

Sample Preparation, Measurement, and Data Processing.

Samples for the analysis of metabolic composition were collected together with the material for dry biomass analysis at 15 days after sowing. Harvested material (shoot and leaf) was cooled below −80°C immediately and kept at this temperature until further processing. Derivatization, GC/MS analysis, and data processing were done as described elsewhere (36). All 181 metabolic signatures that have been evaluated within this experiment are listed in SI Table 4. The GC/MS spectra of evaluated metabolites that are unknown with respect to their actual chemical formula but can be repeatedly found in Arabidopsis are available in SI Table 4.

The extracted metabolite data consist of unique mass intensity values for each referenced compound and measurement respectively. These raw data were normalized and otherwise directly used for analysis. This method allows between sample comparisons but no quantitative statements about single metabolites.

Normalization.

Metabolite data were normalized by dividing each raw value by the median of all measurements of a day for one metabolite.

Missing Value Estimation.

For the canonical correlation analysis (CCA) missing value replacement is necessary. The 6% missing values in the metabolite matrix were imputed with a self-organizing map (SOM) algorithm (37). The mean square error was estimated by the comparison of known values with those calculated from the SOM algorithm. The coefficient of variation (root mean square error divided by the mean) was 0.3.

Integrated Analysis of Phenotypic and Metabolite Data.

Linear models for the relation between metabolite profile and biomass.

The relation between biomass and metabolite profile was measured by simple Spearman correlation between the dry biomass and relative abundances of all metabolites, and by a more complex multiplicative model. The first corresponds to the following model, referred to as model 1: Embedded Image The second model can be described by: Embedded Image B denotes the biomass, x the relative metabolite abundance and c the corresponding constants for all i metabolites.

Multivariate linear analysis.

Canonical correlation analysis (CCA) calculates the highest possible correlation between linear combinations of the columns from two matrices with the same number of rows. If the second matrix has only one column, this procedure corresponds to a ordinary least square (OLS) regression. The correlation thus found is called canonical correlation, the corresponding linear combination canonical variate. The mathematical foundation is described in the literature (23, 37). The R function cancor was used to calculate the canonical correlation between metabolites and biomass. For cross-validation a partial least square (PLS) regression was performed. This method (38) seeks to maximize the covariance instead of the correlation between the matrices. To carry out the procedure the R function plsr was used. These functions are publicly available (www.r-project.org). All procedures were applied after missing value estimation followed by normalization of the metabolic matrix. To test the robustness of the selection of the signature metabolites, we applied the following procedure: with 90% of the 1,144 genotypes chosen at random the canonical variate was calculated and the important metabolites selected as described above. Selected metabolites, which were not in the original list of 44 metabolites, were regarded as false. This procedure was repeated 100 times. We obtained a median “false positive rate” of 0.048 (±0.034).

Acknowledgments

We thank Hans-Peter Piepho for the development of the mixed model for the estimation of mean biomass values and Christine Hackett for useful discussions. We are grateful to Wolfram Stacklies who performed the missing value estimation. We thank Anke Kalkbrenner, Cindy Marona, Melanie Teltow, and Monique Zeh for excellent technical assistance, and Katrin Seehaus and Dirk Zerning for plant cultivation. This work was supported by Deutsche Forschungsgemeinschaft (DFG) Grants AL 387/6-1 and -6-2 (to T.A. and R.C.M.), European Community Grant QLG2-CT-2001-01097 (to T.A.), DFG Grants FI 842/2-1 and WI 550/3-2 (to O.F. and L.W.), and the Max-Planck-Society.

Footnotes

  • ↵†To whom correspondence should be addressed. E-mail: rcmeyer{at}uni-potsdam.de
  • Author contributions: R.C.M., M.S., and J.L. contributed equally to this work; R.C.M., M.S., O.F., L.W., J.S., and T.A. designed research; R.C.M., M.S., J.L., M.B., H.W.-W., O.T., and Ä.E. performed research; J.L. contributed new reagents/analytic tools; R.C.M., M.S., J.L., and J.S. analyzed data; and R.C.M., M.S., J.L., L.W., J.S., and T.A. wrote the paper.

  • The authors declare no conflict of interest.

  • This article is a PNAS direct submission.

  • This article contains supporting information online at www.pnas.org/cgi/content/full/0609709104/DC1.

Abbreviations

RIL,
recombinant inbred line;
CCA,
canonical correlation analysis;
PLS,
partial least square;
TCA,
tricarboxylic acid;
QTL,
quantitative trait loci.
  • Received November 1, 2006.
  • © 2007 by The National Academy of Sciences of the USA

Freely available online through the PNAS open access option.

References

  1. ↵
    1. Tonsor SJ,
    2. Alonso-Blanco C,
    3. Koornneef M
    (2005) Plant Cell Environ 28:2–20.
    OpenUrlCrossRef
  2. ↵
    1. Chen S,
    2. Hajirezaei M,
    3. Peisker M,
    4. Tschiersch H,
    5. Sonnewald U,
    6. Börnke F
    (2005) Planta 221:479–492.
    OpenUrlCrossRefPubMed
  3. ↵
    1. Fernie AR,
    2. Tauberger E,
    3. Lytovchenko A,
    4. Roessner U,
    5. Willmitzer L,
    6. Trethewey RN
    (2002) Planta 214:510–520.
    OpenUrlCrossRefPubMed
  4. ↵
    1. Gibon Y,
    2. Bläsing OE,
    3. Hannemann J,
    4. Carillo P,
    5. Hohne M,
    6. Hendriks JHM,
    7. Palacios N,
    8. Cross J,
    9. Selbig J,
    10. Stitt M
    (2004) Plant Cell 16:3304–3325.
    OpenUrlAbstract/FREE Full Text
  5. ↵
    1. Achard P,
    2. Cheng H,
    3. De Grauwe L,
    4. Decat J,
    5. Schoutteten H,
    6. Moritz T,
    7. Van der Straeten D,
    8. Peng JR,
    9. Harberd NP
    (2006) Science 311:91–94.
    OpenUrlAbstract/FREE Full Text
  6. ↵
    1. Gibson S
    (2005) Curr Opin Plant Biol 8:93–102.
    OpenUrlCrossRefPubMed
  7. ↵
    1. Kolbe A,
    2. Tiessen A,
    3. Schlüpmann H,
    4. Paul M,
    5. Ulrich S,
    6. Geigenberger P
    (2005) Proc Natl Acad Sci USA 102:11118–11123.
    OpenUrlAbstract/FREE Full Text
  8. ↵
    1. Schlüpmann H,
    2. Pellny T,
    3. Van Dijken A,
    4. Smeekens S,
    5. Paul M
    (2003) Proc Natl Acad Sci USA 100:6849–6854.
    OpenUrlAbstract/FREE Full Text
  9. ↵
    1. Fernie AR,
    2. Trethewey RN,
    3. Krotzky AJ,
    4. Willmitzer L
    (2004) Nat Rev Mol Cell Biol 5:763–769.
    OpenUrlCrossRefPubMed
  10. ↵
    1. Allen J,
    2. Davey HM,
    3. Broadhurst D,
    4. Heald JK,
    5. Rowland JJ,
    6. Oliver SG,
    7. Kell DB
    (2003) Nat Biotechnol 21:692–696.
    OpenUrlCrossRefPubMed
  11. ↵
    1. Keurentjes JJB,
    2. Fu J,
    3. de Vos CHR,
    4. Lommen A,
    5. Hall RD,
    6. Bino RJ,
    7. van der Plas LHW,
    8. Jansen RC,
    9. Vreugdenhil D,
    10. Koornneef M
    (2006) Nat Genet 38:842–849.
    OpenUrlCrossRefPubMed
  12. ↵
    1. Tarpley L,
    2. Duran AL,
    3. Kebrom TH,
    4. Sumner LW
    (2005) BMC Plant Biol 5:8.
    OpenUrlCrossRefPubMed
  13. ↵
    1. Clayton TA,
    2. Lindon JC,
    3. Cloarec O,
    4. Antti H,
    5. Charuel C,
    6. Hanton G,
    7. Provost JP,
    8. Le Net JL,
    9. Baker D,
    10. Walley RJ,
    11. et al.
    (2006) Nature 440:1073–1077.
    OpenUrlCrossRefPubMed
  14. ↵
    1. Lindon J,
    2. Holmes E,
    3. Bollard M,
    4. Stanley E,
    5. Nicholson J
    (2004) Biomarkers 9:1–31.
    OpenUrlCrossRefPubMed
  15. ↵
    1. Törjék O,
    2. Witucka-Wall H,
    3. Meyer RC,
    4. von Korff M,
    5. Kusterer B,
    6. Rautengarten C,
    7. Altmann T
    (2006) Theor Appl Genet 113:1551–1561.
    OpenUrlCrossRefPubMed
  16. ↵
    1. Meyer RC,
    2. Törjék O,
    3. Becher M,
    4. Altmann T
    (2004) Plant Physiol 134:1813–1823.
    OpenUrlAbstract/FREE Full Text
  17. ↵
    1. Alonso-Blanco C,
    2. Koornneef M,
    3. Stam P
    (1998) Methods Mol Biol 82:137–146.
    OpenUrlPubMed
  18. ↵
    1. Somerville C
    (2006) Science 312:1277.
    OpenUrlAbstract/FREE Full Text
  19. ↵
    1. Schubert C
    (2006) Nat Biotechnol 24:777–784.
    OpenUrlCrossRefPubMed
  20. ↵
    1. Kopka J,
    2. Schauer N,
    3. Krüger S,
    4. Birkemeyer C,
    5. Usadel B,
    6. Bergmüller E,
    7. Dörmann P,
    8. Gibon Y,
    9. Stitt M,
    10. Willmitzer L,
    11. et al.
    (2005) Bioinformatics 21:1635–1638.
    OpenUrlAbstract/FREE Full Text
  21. ↵
    1. Schauer N,
    2. Steinhauser D,
    3. Strelkov S,
    4. Schomburg D,
    5. Allison G,
    6. Moritz T,
    7. Lundgren K,
    8. Roessner-Tunali U,
    9. Forbes MG,
    10. Willmitzer L,
    11. et al.
    (2005) FEBS Lett 579:1332–1337.
    OpenUrlCrossRefPubMed
  22. ↵
    1. Gittins R
    (1985) Canonical Analysis: A Review with Applications in Ecology (Springer, Berlin).
  23. ↵
    1. Hotelling H
    (1935) J Educ Psychol 26:139–143.
    OpenUrlCrossRef
  24. ↵
    1. Laudadio T,
    2. Pels P,
    3. De Lathauwer L,
    4. van Hecke P,
    5. van Huffel S
    (2005) Magnet Reson Med 54:1519–1529.
    OpenUrlCrossRef
  25. ↵
    1. Zwiers FW,
    2. von Storch H
    (2004) Int J Climatol 24:665–680.
    OpenUrlCrossRef
  26. ↵
    1. Frank IE,
    2. Friedman JH
    (1993) Technometrics 35:109–135.
    OpenUrlCrossRef
  27. ↵
    1. South-East Swedish Breast Cancer Study Group,
    2. Razavi AR,
    3. Gill H,
    4. Stal O,
    5. Sundquist M,
    6. Thorstenson S,
    7. Ahlfeldt H,
    8. Shahsavar N
    (2005) BMC Med Inform Decis Mak 5:29.
    OpenUrlCrossRefPubMed
  28. ↵
    1. Thimm O,
    2. Bläsing O,
    3. Gibon Y,
    4. Nagel A,
    5. Meyer S,
    6. Krüger P,
    7. Selbig J,
    8. Müller LA,
    9. Rhee SY,
    10. Stitt M
    (2004) Plant J 37:914–939.
    OpenUrlCrossRefPubMed
  29. ↵
    1. Liso R,
    2. Calabrese G,
    3. Bitonti M,
    4. Arrigoni O
    (1984) Exp Cell Res 150:314–320.
    OpenUrlCrossRefPubMed
  30. ↵
    1. Smirnoff N
    (2000) Curr Opin Plant Biol 3:229–235.
    OpenUrlCrossRefPubMed
  31. ↵
    1. Hageman GJ,
    2. Stierum RH
    (2001) Mutat Res Fund Mol M 475:45–56.
    OpenUrlCrossRef
  32. ↵
    1. Tkachenko AG,
    2. Pshenichnov MR,
    3. Nesterova LY
    (2001) Microbiology 70:422–428.
    OpenUrlCrossRef
  33. ↵
    1. Garg AK,
    2. Kim JK,
    3. Owens TG,
    4. Ranwala AP,
    5. Do CY,
    6. Kochian LV,
    7. Wu RJ
    (2002) Proc Natl Acad Sci USA 99:15898–15903.
    OpenUrlAbstract/FREE Full Text
  34. ↵
    1. Schauer N,
    2. Semel Y,
    3. Roessner U,
    4. Gur A,
    5. Balbo I,
    6. Carrari F,
    7. Pleban T,
    8. Perez-Melis A,
    9. Bruedigam C,
    10. Kopka J,
    11. et al.
    (2006) Nat Biotechnol 24:447–454.
    OpenUrlCrossRefPubMed
  35. ↵
    1. Piepho HP,
    2. Büchse A,
    3. Emrich K
    (2003) J Agron Crop Sci 189:310–322.
    OpenUrlCrossRef
  36. ↵
    1. Lisec J,
    2. Schauer N,
    3. Kopka J,
    4. Willmitzer L,
    5. Fernie AR
    (2006) Nat Protoc 1:387–396.
    OpenUrlCrossRefPubMed
  37. ↵
    1. Kuss M,
    2. Graepel T
    (2003) Technical Report 108 (Max Planck Inst Biol Cybern, Tübingen, Germany).
  38. ↵
    1. Wold H
    (1975) Soft Modelling by Latent Variables (Academic, London).
PreviousNext
Back to top
Article Alerts
Email Article

Thank you for your interest in spreading the word on PNAS.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
The metabolic signature related to high plant growth rate in Arabidopsis thaliana
(Your Name) has sent you a message from PNAS
(Your Name) thought you would like to see the PNAS web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
The metabolic signature related to high plant growth rate in Arabidopsis thaliana
Rhonda C. Meyer, Matthias Steinfath, Jan Lisec, Martina Becher, Hanna Witucka-Wall, Ottó Törjék, Oliver Fiehn, Änne Eckardt, Lothar Willmitzer, Joachim Selbig, Thomas Altmann
Proceedings of the National Academy of Sciences Mar 2007, 104 (11) 4759-4764; DOI: 10.1073/pnas.0609709104

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Request Permissions
Share
The metabolic signature related to high plant growth rate in Arabidopsis thaliana
Rhonda C. Meyer, Matthias Steinfath, Jan Lisec, Martina Becher, Hanna Witucka-Wall, Ottó Törjék, Oliver Fiehn, Änne Eckardt, Lothar Willmitzer, Joachim Selbig, Thomas Altmann
Proceedings of the National Academy of Sciences Mar 2007, 104 (11) 4759-4764; DOI: 10.1073/pnas.0609709104
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Mendeley logo Mendeley
Proceedings of the National Academy of Sciences: 104 (11)
Table of Contents

Submit

Sign up for Article Alerts

Jump to section

  • Article
    • Abstract
    • Results
    • Discussion
    • Conclusion
    • Materials and Methods
    • Acknowledgments
    • Footnotes
    • Abbreviations
    • References
  • Figures & SI
  • Info & Metrics
  • PDF

You May Also be Interested in

Smoke emanates from Japan’s Fukushima nuclear power plant a few days after tsunami damage
Core Concept: Muography offers a new way to see inside a multitude of objects
Muons penetrate much further than X-rays, they do essentially zero damage, and they are provided for free by the cosmos.
Image credit: Science Source/Digital Globe.
Water from a faucet fills a glass.
News Feature: How “forever chemicals” might impair the immune system
Researchers are exploring whether these ubiquitous fluorinated molecules might worsen infections or hamper vaccine effectiveness.
Image credit: Shutterstock/Dmitry Naumov.
Venus flytrap captures a fly.
Journal Club: Venus flytrap mechanism could shed light on how plants sense touch
One protein seems to play a key role in touch sensitivity for flytraps and other meat-eating plants.
Image credit: Shutterstock/Kuttelvaserova Stuchelova.
Illustration of groups of people chatting
Exploring the length of human conversations
Adam Mastroianni and Daniel Gilbert explore why conversations almost never end when people want them to.
Listen
Past PodcastsSubscribe
Horse fossil
Mounted horseback riding in ancient China
A study uncovers early evidence of equestrianism in ancient China.
Image credit: Jian Ma.

Similar Articles

Site Logo
Powered by HighWire
  • Submit Manuscript
  • Twitter
  • Facebook
  • RSS Feeds
  • Email Alerts

Articles

  • Current Issue
  • Special Feature Articles – Most Recent
  • List of Issues

PNAS Portals

  • Anthropology
  • Chemistry
  • Classics
  • Front Matter
  • Physics
  • Sustainability Science
  • Teaching Resources

Information

  • Authors
  • Editorial Board
  • Reviewers
  • Subscribers
  • Librarians
  • Press
  • Cozzarelli Prize
  • Site Map
  • PNAS Updates
  • FAQs
  • Accessibility Statement
  • Rights & Permissions
  • About
  • Contact

Feedback    Privacy/Legal

Copyright © 2021 National Academy of Sciences. Online ISSN 1091-6490