A Topological Map of the Compartmentalized Arabidopsis thaliana Leaf Metabolome

Abstract

Background

The extensive subcellular compartmentalization of metabolites and metabolism in eukaryotic cells is widely acknowledged and represents a key factor of metabolic activity and functionality. In striking contrast, the knowledge of actual compartmental distribution of metabolites from experimental studies is surprisingly low. However, a precise knowledge of, possibly all, metabolites and their subcellular distributions remains a key prerequisite for the understanding of any cellular function.

Methodology/Principal Findings

Here we describe results for the subcellular distribution of 1,117 polar and 2,804 lipophilic mass spectrometric features associated to known and unknown compounds from leaves of the model plant Arabidopsis thaliana. Using an optimized non-aqueous fractionation protocol in conjunction with GC/MS- and LC/MS-based metabolite profiling, 81.5% of the metabolic data could be associated to one of three subcellular compartments: the cytosol (including the mitochondria), vacuole, or plastids. Statistical analysis using a marker-‘free’ approach revealed that 18.5% of these metabolites show intermediate distributions, which can either be explained by transport processes or by additional subcellular compartments.

Conclusion/Significance

Next to a functional and conceptual workflow for the efficient, highly resolved metabolite analysis of the fractionated Arabidopsis thaliana leaf metabolome, a detailed survey of the subcellular distribution of several compounds, in the graphical format of a topological map, is provided. This complex data set therefore does not only contain a rich repository of metabolic information, but due to thorough validation and testing by statistical methods, represents an initial step in the analysis of metabolite dynamics and fluxes within and between subcellular compartments.

About the Authors:

Stephan Krueger

Contributed equally to this work with: Stephan Krueger, Patrick Giavalisco, Dirk Steinhauser

Affiliation: Botanical Institute, University of Cologne, Cologne, Germany

Patrick Giavalisco

Contributed equally to this work with: Stephan Krueger, Patrick Giavalisco, Dirk Steinhauser

Affiliation: Max Planck Institute of Molecular Plant Physiology, Potsdam-Golm, Germany

Leonard Krall

Affiliation: Max Planck Institute of Molecular Plant Physiology, Potsdam-Golm, Germany

Marie-Caroline Steinhauser

Affiliation: Max Planck Institute of Molecular Plant Physiology, Potsdam-Golm, Germany

Dirk Büssis

Affiliation: GABI Managing Office, c/o Max Planck Institute of Molecular Plant Physiology, Potsdam-Golm, Germany

Bjoern Usadel

Affiliation: Max Planck Institute of Molecular Plant Physiology, Potsdam-Golm, Germany

Ulf-Ingo Flügge

Affiliation: Botanical Institute, University of Cologne, Cologne, Germany

Alisdair R. Fernie

Affiliation: Max Planck Institute of Molecular Plant Physiology, Potsdam-Golm, Germany

Lothar Willmitzer

Affiliation: Max Planck Institute of Molecular Plant Physiology, Potsdam-Golm, Germany

Dirk Steinhauser

Contributed equally to this work with: Stephan Krueger, Patrick Giavalisco, Dirk Steinhauser

* E-mail: [email protected]

Affiliation: Max Planck Institute of Molecular Plant Physiology, Potsdam-Golm, Germany

Introduction

The partitioning of cellular functions and metabolism into subcellular compartments is a fundamental feature of all eukaryotic cells. Subcellular compartments are usually delineated by a lipid bilayer to maintain compartment integrity and specific microenvironments. Though physically and biochemically distinct, these compartments and their metabolic contents are interlinked by inter-compartmentally transported metabolites [1], [2], [3]. This translocation, as well as the turnover of metabolites, can be exceptionally fast [2], [4], making the reliable determination of metabolites in subcellular compartments challenging. Consequently, the development of methods and strategies to determine the metabolic composition of these compartments is required to gain a comprehensive understanding of the cellular biochemistry.

While subcellular distributions have been determined for a limited number of metabolites using genetically encoded metabolic sensors [5], [6] or direct mass imaging methods on surface tissues [7], the number of approaches devoted towards deciphering subcellular distributions of multiple metabolites is rather limited. The main challenge using destructive approaches is that in order to prevent the leakage of metabolites out of organelles the analysis needs to be performed under anhydrous conditions thus rendering subcellular metabolite analyses strikingly different from e.g. organelle oriented proteomic studies [8].

Non-aqueous fractionation (NAF) is a powerful technique to separate subcellular compartments, and their molecular compositions, under conditions where biological activities are completely arrested due to rapid freezing and dehydration of the sample material [9], [10], [11]. Cellular constituents in proximity to each other aggregate to small particles during lyophilization of the ground sample material. These particles, mainly fragments of cellular compartments, are then separated by their composition-dependent density using equilibrium centrifugation in a non-aqueous gradient. Using compartment-specific marker abundances throughout collected gradient fractions, compartment enrichment and compartmental separation can be assessed. As well, subcellular metabolite distributions can be calculated, usually by applying a two- or three-compartmental calculation strategy [10], [11].

While NAF was first applied to study animal nuclei [9] and mitochondrial high-energy phosphates in mammalian cells [12], [13], [14], later on it has been used mostly in plant sciences in order to determine the partitioning of photosynthetic assimilates in leaves [10], [11], [15], [16], [17], storage organs [18], [19], [20], rose petals [21] or for analysis of specific pathways [22].

In the past decade, technological breakthroughs in mass spectrometry (MS) and nuclear magnetic resonance spectrometry (NMR) [23] have paved the way for comprehensive analyses of an organism's metabolic composition [24], [25]. Even though NMR provides advantages for quantitative and structural metabolomics [26], LC/MS- and GC/MS-based metabolite profiling have become the methods of choice for a general overview of cellular metabolism due to their high throughput, compound coverage, and sensitivity [27]. Despite the increasing use of MS-based metabolite profiling, it has only been combined with NAF in a limited number of studies, basically to unravel the subcellular location of primary metabolites in soybean leaves and potato tubers by means of targeted analyses [20], [28], [29].

Here, we describe the subcellular distribution of a broad range of polar and lipophilic compounds in leaves of the model plant Arabidopsis thaliana obtained using three orthogonal MS-based analytical approaches, namely GC-TOF/MS for primary and LC-FT/MS analyses for lipids and semipolar, secondary metabolites. The provided data, which can be regarded as a resource documenting a metabolomic survey of a compartmentally separated leaf, clearly distinguishes the cytosol, the plastids, and the vacuole from one another. Using statistical approaches we were able to demonstrate the robustness of our analyses, assign chemical compounds to the resolved compartments, and to validate our results using structurally annotated (known) metabolites. We further demonstrate that the localizations of several known metabolites and structurally undetermined compounds (unknowns) are difficult to unambiguously explain on the basis of three compartments due to either unresolved compartments, or the interconnections of subcellular metabolic networks.

Results and Discussion

Non-aqueous fractionation of Arabidopsis leaves allows clear separation of three subcellular compartments

In order to analye the subcellular compartmentalization of the plant metabolome, NAF was performed on three independent replicates of pooled Arabidopsis leaves from soil grown plants harvested three hours after the onset of light using an optimized NAF protocol [22] (Text S1). NAF separates fragments of subcellular compartments and organelles in a continuous density gradient. Due to the variable composition-dependent density of the fragments, their segregation reflects continuous compartmental distributions throughout the gradient [11]. To unambiguously assign a specific compartment to these distributions, abundances of compartment-specific markers within the six collected gradient fractions were determined. These marker distributions, which must be sufficiently distinct from each other, were then used to evaluate the compartmental enrichment and separation of distinct organelles or subcellular spaces (Figure 1). Nitrate, as vacuolar marker [22], showed a clear enrichment in the densest fraction F1 with 40.1±2.1% (as mean ± SD) which is in agreement with the vacuolar H⁺-ATPase abundance (Figure 1). The cytosolic marker UGPase [30] was relatively equally distributed across the gradients with abundances ranging from 12.3±2.6% to 23.3±4.2%, showing slight increases (F1: 18.9±1.5%; F6: 23.3±4.2%) in the most distant fractions (Figure 1). Contrarily to nitrate, the plastidic marker NADP-GAPDH [11] was clearly enriched in the lightest fraction F6, with 66.8±3.4%, which is in agreement with the abundance of the light harvesting complex (LHC) (Figure 1). Citrate synthase, used as a mitochondrial marker [31] was detected throughout the gradients and revealed a similar distribution as observed for the cytosol, but with decreased abundance in fraction F1 (9.3±2%) and an enrichment in fraction F6 (34.8±9.7%; Figure 1).

[Figure omitted. See PDF.]

Figure 1. Distribution of compartment-specific markers in non-aqueous gradients from Arabidopsis thaliana leaves.

(A) The distribution of vacuolar (nitrate), cytosolic (UGPase), mitochondrial (citrate synthase), and plastidic (GAPDH) markers are shown as the average of three independent gradients. The mean values and standard deviations of marker enzyme activities (UGPase, citrate synthase, GAPDH) or relative concentrations (nitrate) in each fraction are depicted as percentage from total (scaled data). Significant differences (P_BH<0.05, Benjamini-Hochberg corrected) using t-test within each fraction compared to the cytosolic (cyt.) or mitochondrial (mit.) marker are shown as green colored boxes below the graph. Grey boxes illustrate uncorrected significant (P<0.05) differences. (B) Western blots detecting LHC (plastidic) and vacuolar H⁺-ATPase (vacuolar) membrane proteins in each fraction are shown for one representative gradient to confirm the distribution of the plastidic and vacuolar compartment throughout the gradients. The pixel intensities quantified using ImageJ (http://rsb.info.nih.gov/ij) are drawn as bar diagrams.

https://doi.org/10.1371/journal.pone.0017806.g001

Despite the clear intermediate distribution between the cytosolic and plastidic compartment, the mitochondrial marker revealed a relatively large standard deviation and was not enriched in any fraction as compared to the other markers (Figure 1). Therefore, and in agreement with previous reports [e.g. 10], [20], [22], the mitochondrial compartment was not, even though a clear trend could be observed, considered to be unambiguously delineated from the cytosolic compartment. However, with the broad separation of the other markers we were clearly able to obtain an excellent separation of the vacuolar, the cytosolic, and the plastidic compartments by non-aqueous fractionation of Arabidopsis leaf material.

Non-aqueous fractionation produces consistent fraction separation

A total of 18 fractions, resulting from three independent gradients comprised of six fractions each, were subjected to the three MS platforms for polar and lipophilic metabolite analysis (Data S1, S2, S3). In the following, the MS data refers to the analytical approach applied rather than the exact chemical properties of the detected metabolites.

Using GC/MS, 203 analytes, comprising 88 unique metabolites, were consistently identified with 93 (45.8%) and 110 (54.2%) analytes of known and unknown chemical structure (Data S4). High-resolution LC/MS analyses of lipophilic and secondary metabolites resulted in the consistent monitoring of 2,804 and 910 mass spectrometric features (afterwards analytes) comprising 726 and 461 non-redundant peaks (T/S clusters, Text S1). Database searches revealed that 88 (3.1%) and 31 (3.4%) analytes of lipophilic and secondary metabolite profiling represented known metabolites and further 362 (12.9%) and 224 (24.6%) produced database hits with single or multiple potential chemical structures (Data S4).

In order to test whether the individual metabolome data are consistent among the independent gradients and discriminative with respect to fraction separation, principal component (PCA) and hierarchical cluster (HCA) analyses were performed on the scaled and, for the PCA, additionally log₂-transformed metabolite data (Figure 2). Non-parametric ANOVA using the Mantel test [32] supported separation of the six fraction groups despite low matrix correlations of r = 0.43, r = 0.28, and r = 0.45 (P<0.001) for primary, lipophilic, and secondary metabolite data, respectively. Sequential expanding of the fraction grouping (Figure S1) showed significant (P<0.001, r = 0.97) differences between the F6 and the other fractions for lipophilic metabolites (Figure 2E). Primary (r = 0.84) and secondary (r = 0.89) metabolite data statistically (P<0.001) support the distance of the plastidic (F6) and vacuolar (F1) enriched fractions from the remaining ones, even though relatively high matrix correlations (r = 0.63 and r = 0.66, P<0.001) are observed if the two further clusters, comprising the intermediate-dense fractions F2-F3 and F4-F5, are not merged (Figure S1). Mantel tests between sample distance matrices, to determine the overall similarity in terms of similar fraction separation underlying the different metabolome data, showed a very high correlation of r = 0.92 (P<0.001) between primary and secondary metabolite data. However, both primary and secondary metabolite data revealed significant but lower correlations with r = 0.82 (P<0.001) and r = 0.63 (P<0.005), respectively, when compared to the lipophilic metabolite data.

[Figure omitted. See PDF.]

Figure 2. Principal component (PCA; A–C) and hierarchical cluster (HCA; D–F) analyses of metabolite data.

Both PCA and HCA plots demonstrate a good separation of the six fractions from each other independent of the three major compound classes, (A, D) primary, (B, E) lipophilic and (C, F) secondary metabolites. PCA and HCA were performed on scaled data. For PCA, data were additionally log₂-transformed; HCA analyses are based on Euclidean distances among samples. Identical gradient fractions are encoded with the same color and shape as depicted in the graph legend. The 95% confidence ellipse is drawn as a grey dotted line on the basis of the mean, standard deviation, and correlation of the three independent gradient replicates per fraction. To aid interpretation of HCA graphs (D–F) same fractions are identically color-coded (see PCA legend) at the bottom sidebar. The unbiased cluster P-values, calculated using multiscale bootstrap resampling, are depicted as red-colored numbers at each node. The fraction groups explaining the highest variance of data and revealing a good cohesion within and separation among assigned fractions are depicted at the bottom of the HCA plots. Fraction groups, evaluated using resampling and gap statistic, were assembled by sequential merging of neighboring sample clusters with fractions assigned using membership majority voting (Figure S1). All graphs clearly support that the plastidic enriched, lowest density fraction F6 (group A) and to a lesser extent the vacuolar enriched, most dense fraction F1 are separated from the intermediately-dense fractions F2 to F5, even though for primary and secondary metabolite data two less-well separated fraction groups (F2–F3 and F4–F5) can be assumed (solid and dotted line; Figure S1).

https://doi.org/10.1371/journal.pone.0017806.g002

In essence, the three metabolite data sets showed consistency among the data derived from the independent gradients and supported, visually (Figure 2) and statistically (see above; Figure S1), the separation of compartmental-enriched fractions. The plastidic (F6) and to a lesser extent the vacuolar enriched fraction (F1) are distinct from the majority of the intermediate-dense fractions (F2–5). Gap statistics suggested overall less well-separated clusters (Figure S1D–F), likely because there is a continuous distribution of compartments and their metabolite content throughout the gradients (Figure 1).

Non-aqueous fractionation results in robust compartmental fractionation

As NAF in combination with the MS-based analysis is a complex procedure, various error sources can affect the compartmental separation and downstream estimation of subcellular distributions. To evaluate the consistency and robustness of NAF-derived data further markers were measured or selected from our MS data. Starch was used as additional plastidic marker, as it is synthesized and stored as semi-crystalline granules in plastids during the day [33]. Digalactosyldiacylglycerol (DGDG), a group of galactolipids with high abundances in both envelope and thylakoid membranes [34] were further utilized as plastidic markers. Many classes of secondary metabolites like glucosinolates and flavonoids are reported to be commonly stored in the vacuole (or vacuolar inclusions) of several different plant species [e.g. 35], [36], [37], [38], [39], [40], [41], [42] and thus, represent ideal vacuolar markers. Therefore, based on a targeted analysis we selected a number of glucosinolates and flavonoids/sinapate esters (afterwards for simplicity called flavonoids) (cf. Data S4) reported to be found in Arabidopsis (KNApSAcK database and references therein [43]; [44]). As additional markers for the cytosol, triacylglycerides [45] and glyceroceramids, a class of lipids within the sphingolipid group localized to the plasma membrane and to a lesser extent to the tonoplast [46], were used. The results for all nine marker distributions are shown in Figure 3 and demonstrate a high reproducibility of the marker distribution between the three gradients. Likewise, the between-gradient variation of markers designating the same compartment are relatively small with coefficients of variation of 29±8.5%, 19.1±5.8%, and 21.8±7.7% for the plastidic, cytosolic, and vacuolar compartment, respectively (all as mean ± SD; Figure S2A).

[Figure omitted. See PDF.]

Figure 3. Side-by-side bar plots of marker distribution representing the same subcellular compartment throughout the three independent gradients.

For all graphs scaled data were used. Marker names are provided as graph headers including the compartmental representation, i.e. cpl. - plastidic, vac. - vacuolar, and cyt. – cytosolic compartment. The bars are colored as follows: gradient 1 - black, gradient 2 – white, and gradient 3 – grey. For DGDG (cpl.), GlcCer (cyt.), TAG (cyt.), glucosinolates (vac.), and flavonoids (vac.) the abundance per fraction is based on the robust average of multiple analytes representing the individual compound class (Data S4).

https://doi.org/10.1371/journal.pone.0017806.g003

Use of multiple markers results in robust compartmental designation and assignment

As described above, in addition to the three markers used to assign subcellular compartments we took advantage of the fact that within our metabolite data several compounds could be assigned to a specific compartment in an unambiguous way. The availability of these additional markers allowed us to rigorously test the reproducibility of the fractionation procedure and to assess the magnitude of cohesion within and separation between the three delineated compartments.

Classical multidimensional scaling (CMD; Figure 4) and HCA on normalized Manhattan distances (Figure S2A) among markers and gradients clearly demonstrate the separation of the three considered compartments (Figures 4 and S2). The individual clusters reveal high silhouette information with 0.71±0.05, 0.60±0.06, and 0.61±0.05 for the plastidic, cytosolic, and vacuolar compartment, respectively, with a cluster-solution average of 0.64±0.07 (all as mean ± SD). Thus, a high cohesion within and separation among the clusters is observed, which is supported by gap statistic (Figure S2B) and non-parametric ANOVA using Mantel test (P<0.001, r = 0.77). The spread of compartmental clusters, estimated as the clusterwise average of their normalized Manhattan distances, within- and between-gradients is very similar (data not shown), possessing low between-gradient cluster spreads of 0.12±0.05, 0.11±0.03, and 0.11±0.03 (all as mean ± SD) for the plastidic, cytosolic and vacuolar compartment, respectively. Interestingly, the plastidic compartment revealed the largest cluster diameter based on the maximum clusterwise normalized Manhattan distance (0.16±0.02; 0.24) compared to the cytosolic (0.13±0.03; 0.17) and vacuolar (0.14±0.01; 0.19) compartment (within-gradient diameter as mean ± SD followed by between-gradient diameter).

[Figure omitted. See PDF.]

Figure 4. Consistency and robustness of compartmental separation within and between gradients.

Manhattan distances among markers for each gradient were converted into a principal coordinates (PCo) space using classical multidimensional scaling (CMD) for the three independent and 729 non-redundant combinations of randomly assembled gradients. Shapes colored in magenta show the data points for the three independent gradients (G1, G2, and G3 as depicted in the figure) of each selected marker. Data points from simulated gradients are depicted as circles with coloration according to the individual compartments: green – chloroplasts, blue – cytosol, and grey – vacuole. For each compartment the 95% confidence ellipse is drawn as a dashed grey line on the basis of the mean, standard deviation, and correlation of the 729 non-redundant gradient combinations. The mitochondrial compartment (yellow circles) shows an overlapping distribution with the cytosol where the majority of data are in between the plastidic and cytosolic compartments. The principal coordinates 1 and 2 explain together in average about 96.7% of the total variance of the underlying distance matrices. For each of the three resolved compartments the distribution of average abundances including their standard deviations throughout the three independent gradients is depicted as bar plots.

https://doi.org/10.1371/journal.pone.0017806.g004

For further robustness evaluation, fractions were systematically assembled into all possible non-redundant artificial combinations (simulated gradients) and the computed normalized Manhattan distances subjected to CMD analyses. Markers representing the same compartment (Figure 4) are in close proximity to each other and do not reveal large variations within their distribution in principal coordinates space. With respect to compartmental clusters, the 95% confidence ellipses and gap statistic (data not shown) clearly support their separation as neither they, nor any data point overlap with any other cluster (Figure 4).

As initially observed, the mitochondrial marker citrate synthase is clearly distributed between the plastid and cytosol with overlap towards the cytosol (Figure 4). When including the mitochondrial compartment the silhouette information of the cluster solution drops (0.5±0.2) with clusterwise values of 0.55±0.10, −0.01±0.10, 0.50±0.14, and 0.61±0.05 for the plastidic, mitochondrial, cytosolic, and vacuolar compartment, respectively (all as mean ± SD). The mitochondrial compartment itself revealed a large between-gradient cluster spread of 0.20±0.003 and diameter of 0.20 compared to the other compartments (see above; Figure 4). In addition, a fourth cluster (mitochondrial) is not fully supported by gap statistic (Figure S2B).

As we lack further unambiguous markers to clearly distinguish the mitochondrial compartment, we therefore prefer to not consider the contribution of this compartment from our NAF gradients. Interestingly, the problem of unequivocally separating mitochondria was already described in previous NAF studies [10], [20], [22]. Reasons for this might be the small size and dispersed localization of mitochondria in the plant cell, as they exist as a population of physically discrete organelles [47]. They are also highly motile within the cell, associating with specific compartments, as seen under stress conditions [48], through association with the actin cytoskeleton [49]. In consequence, the cytosolic compartment must be considered in a broader sense as it represents metabolites with clear cytosolic and/or possible mitochondrial localization.

Different computational strategies to estimate subcellular distributions result in similar downstream findings and demonstrate statistical robustness

As the compartment-specific markers support an enrichment and robust separation of the three compartments, they enabled a marker-based determination of subcellular metabolite distributions for these compartments (Data S4).

Two different computational algorithms were used, namely the non-negative least square (NNLS) algorithm by Lawson and Hanson [50] and the best fit algorithm (BFA) by Riens [11]. Essentially, both approaches solve a system of linear equations defined by the subcellular compartments (designated by marker abundances throughout the gradient) to find the relative assignment of a metabolite to the compartments by minimizing the discrepancy between the measured and computationally estimated (fitted) fraction abundances of this metabolite. NNLS is based on an active-set approach and seeks linear least square solutions that are also non-negative by minimizing the Euclidean distance [50], which corresponds to the square root of the residual sum of squares (Eq. 1). BFA is based on a heuristic approach and tests all possible subcellular distributions using 1% intervals, i.e. (1^st) vacuole 100%, cytosol 0%, plastid 0%, (2^nd) vacuole 99%, cytosol 1%, plastid 0% and so forth, by minimizing the Q-value (Eq. 2), the Euclidean distance divided by the number of fractions – 1 [11]. Whereas BFA solutions add up to 100% across the considered compartments, the NNLS solutions, as no other constraints other than non-negative values are given, can sum to above or below 100% (Data S4).

Application of both BFA and NNLS revealed that the average of the estimated subcellular distributions from the three independent gradients are very similar, characterized by a mean difference of about −0.1±2.3% (as mean ± SD) where 95% of the BFA-to-NNLS differences lay in a range of −5% to 4% with tailing at about ≥100% due to the abovementioned algorithms differences (Figure 5A). Comparison of the BFA solutions derived from the independent and simulated gradient data revealed overall small differences of the averaged subcellular distributions with a mean difference of about 0±1.2% (as mean ± SD) where 95% of the differences lay in a range of −3% to 2% (Figure 5B). A similar behavior was observed when using the NNLS solutions (cf. Data S4).

[Figure omitted. See PDF.]

Figure 5. Diagnostic plots showing the differences in estimated subcellular distributions using (A) BFA and NNLS algorithm on the three independent gradients, and (B) using BFA algorithm on the three independent and 729 simulated gradients.

The difference versus average plot (left) shows the differences (D) in dependence of the averages (A) of estimated subcellular distributions in a compartment (C, estimated as percentage) between two computational strategies (S): D_i = C_{i [S1]} – C_{i [S2]}, and A_i = 0.5 x (C_{i [S1]} + C_{i [S2]}). The corresponding computational strategies depicted in the plots are for (A) S1 = BFA and S2 = NNLS solutions, and for (B) S1 = 3 (GRD) and S2 = 729 (SIM) gradients. The histogram plot (right) shows the distribution of the M values in 1% intervals with blue solid lines indicating the 2.5% and 97.5% quantiles (95% range of observed differences). For all comparisons the average of estimated subcellular distributions (based on 3 or 729 gradients) are used.

https://doi.org/10.1371/journal.pone.0017806.g005

As both BFA and NNLS resulted in similar estimates of subcellular distribution, and computation on simulated gradients revealed the overall statistical robustness, we used the BFA solutions estimated on the three independent gradient data for all further analyses.

A three-compartmental distribution strategy sufficiently explains the majority of observed analyte distributions

As criteria for a best fit, the mathematically related measures Euclidean distance, residual sum of squares, or the Q-value are usually considered (see above and Eq. 1–3). A tight fit is reflected in small values of the outputs from these equations, however ‘small’ is difficult to define. Therefore we used the normalized Manhattan distance (Eq. 3), the sum of absolute differences between the measured and fitted data, as it ranges from 0 to 100% (or 0.0 to 1.0 if expressed relative) on scaled NAF data. Thus, it describes the total percentage discrepancy (TPD) between the model and the measurements. Subcellular distributions were considered as insufficiently explained (‘unexplained’, cf. Figure 6) if both the average TPD exceeded 10% and the TPD from individual gradients exceeded 10% in ≥50% of the cases. The 10% cutoff was chosen because the difference of individual markers to their respective compartment-specific average was 6.9±1.8% (as mean ± SD) with a maximum of 10.3%.

[Figure omitted. See PDF.]

Figure 6. Manually constructed classification tree for the partitioning of analytes into compartments and intermediate units based on their estimated subcellular distribution as well as compartmental abundance and variability.

The classification tree was used to assign an analyte into a group defining its subcellular distribution type (type) and an associated mode (mode) which represents one of the resolved compartments or the overlap between them. Assignments are based on the mean and standard deviations of the subcellular distribution, estimated using the best fit algorithm (BFA), for each analyte based on the three independent gradient data. Analytes with insufficiently explained subcellular distributions according to the selected compartment-specific markers are accounted as unexplained and are not further considered in this tree. Analytes revealing sufficiently explained fits are accounted as ‘shared’ if the minimum of the percentage value (min  =  mean - SD) in the most abundant compartment is overlapping with the maximum percentage value (max  =  mean + SD) of any other considered compartment. The corresponding mode is defined by the overlapping compartments regarding the most abundant compartment. Analytes are considered ‘specific’ with the mode according to the estimated most abundant compartment, if the minimum of the most abundant compartment is ≥75% and the values of all other (low abundant) compartments are negligible, i.e. ≤10%. If the minimum in the most abundant compartment is larger than the sum of the maxima among the other compartment and ≥66.67% (2/3 compartments), analytes are accounted as ‘dominantly’ distributed in the respective compartment. Analytes are accounted as ‘enriched’ in a compartment, if the value of the most abundant fraction is ≥50% than the sum of the other compartment. If none of this decisions result in an assignment the analyte is considered as being shared but with enrichment in a particular compartment (‘shared*’).

https://doi.org/10.1371/journal.pone.0017806.g006

Using BFA the subcellular distributions of 3,198 (81.5%) out of all 3,922 analytes are considered as sufficiently explained by the averaged compartment-specific markers using a three-compartmental calculation strategy, considering the vacuolar, cytosolic, and plastidic compartments (Table 1; Data S4 and S5). Consequently, the subcellular distributions of 724 (18.5%) analytes are insufficiently explained. Classification based on the compartment-specificity of the three nearest neighboring markers using k-nearest neighbor algorithm (kNN, with k = 3) facilitated the assignment of 487 (67.3%), 174 (24.0%), and 63 (8.7%) of these insufficiently explained analytes into the cytosolic, plastidic, and vacuolar compartments, respectively (Table 1). Interestingly, the mitochondrial marker citrate synthase was considered as insufficiently explained but, as mentioned before, it showed a more cytosol-like distribution and therefore was assigned to the cytosol (Data S4). Similarly, sucrose, a metabolite which is synthesized in the cytosol and transported into sink organs via the phloem, was assigned to the cytosol (Data S4). Earlier conducted NAF studies [10], [11] have already indicated that the observed sucrose distribution could not clearly be ascribed to the cytosolic, plastidic, or vacuolar compartment, most likely due to the greatly higher amounts present in sieve tubes [11].

[Figure omitted. See PDF.]

Table 1. Overview of compartmental assignment results.

https://doi.org/10.1371/journal.pone.0017806.t001

Therefore, analytes revealing insufficiently explained subcellular distributions may indicate the presence and influence of unconsidered compartments, as the compartment-specific markers used do not encompasses their distribution and thus are inadequate to precisely explain the observed distributions.

Guilt by association – a three-compartmental distribution strategy facilitates compartmental classification of metabolites

Several metabolites are known to localize to more than one specific compartment; similarly, in our analysis we have observed compounds present in more than one compartment. Explanations for these observations can be the occurrence of biochemical pathways in multiple compartments as well as transport of compounds between compartments. In order to account for these situations, analytes were classified into specific compartments or intermediate assignments based on their compartmental abundance using a classification tree (Figure 6). Five classes were considered:

1. ‘specific’,

2. ‘dominant’ where, for both, the analyte pool sizes are located dominantly but to different degrees within a designated compartment,

3. ‘enriched’ where the pool size in the most abundant compartment is roughly higher than the sum of the others,

4. ‘shared’ between compartments,

5. shared with enrichment in a specific compartment (‘shared*’).

These classification results (Table 1; Figure S3) have been visualized as a topological map of the compartmentalized metabolome (Figure 7; for single analytes see Data S6). In detail, 82 GC/MS analytes (40.4%) were classified as specific or dominant and 47 (23.2%) as shared between compartments. Of the specific or dominant class, 48.8% of analytes were assigned to the cytosol, 22% localized to the plastids, and 29.3% were localized to the vacuole. Of the 488 (53.6%) specific and dominantly assigned analytes derived from secondary metabolism, 48% were assigned to the vacuole, 42% to the cytosol, and only 10% were localized to the plastid. 1,657 (59.1%) out of all lipophilic analytes displayed specific and dominant subcellular distributions. Of them, the majority, 63%, were assigned to the cytosol and 36.8% to the plastids. Lipophilic compounds showing specific or dominant pool sizes in the vacuole are negligible (0.2%).

[Figure omitted. See PDF.]

Figure 7. A topological map of the compartmentalized Arabidopsis thaliana leaf metabolome for (A) all and (B) selected analytes.

The classification for the partitioning of analytes into compartments and intermediate units is based on the best fit - estimated subcellular distribution as well as compartmental abundance and variability for each analyte (for details see Figure 6). The topological map (cf. Data S6 for single analytes) of the classification results for (A) all and (B) selected metabolites is visualized in principal coordinates (PCo) space on the basis of averaged Manhattan distances among analytes for the three independent gradients. To aid interpretation, analytes of the classes ‘specific’ and ‘dominant’ were both assigned into the respective compartment and color-coded accordingly: green – chloroplast, blue – cytosol, and grey – vacuole. Analytes assigned as being shared between two compartments are color coded as depicted in the figure. With exception of analytes with insufficiently explained (unexplained) subcellular distributions, all other analytes not belonging to one of the above-mentioned classes are defined as ‘others’.

https://doi.org/10.1371/journal.pone.0017806.g007

Overall the compartmental assignment varies regarding the major compound classes (Table 1; Table S1, S2, S3). Many lipophilic metabolites can be found localized to plastids and the cytosol, as both encompass large internal membrane systems. As well, plastids are the site of plant fatty acids synthesis [34], [51]. In contrast, many secondary metabolites are dominant or even specific for the vacuole and cytosol, reflecting their synthesis and storage location as has been supported by protein localization studies [e.g. 36], [41], [52]. Analytes of primary metabolism revealed the largest diversity regarding compartmental class assignments as these compounds are crucial constituents for many pathways. They are localized to each compartment, with 29.6% analytes revealing shared pool sizes.

Literature confirmation of selected metabolites demonstrates robustness, relevance, and facilitates hypothesis deduction

As described above we followed a comprehensive approach to assign as many analytes into specific compartments as possible. To the best of our knowledge, this study, with respect to its comprehensiveness, is the first of its kind. Therefore we decided to validate the data by linking it to prior knowledge.

In our study, most of the amino acids were highly abundant in the cytosol and chloroplasts (Table S1), which is in agreement with results obtained for leaves of other plant species [11], [53]. Proline, which is synthesized in chloroplasts and the cytosol of mesophyll cells [54], is dominantly plastidic localized (67±6%; Table S1). This localization fits its function as ROS scavenger and singlet oxygen quencher during photosynthesis [55]. Methylglucopyranoside (MeG), a secondary metabolite synthesized by direct transfer of methanol onto glucose in the cytosol of Geum montanum, is rapidly transported into the vacuole where it accumulates to more than 95% [56]. The high vacuolar abundance of MeG (95±9%; Table S1) indicates that MeG metabolism in Arabidopsis might be similar, at least in terms of storage. Recently, it was shown that myo-inositol accumulates in the cytosol and not in the vacuole of Mesembryanthemum crystallinum [57] supporting its cytosolic localization (100±1%) in Arabidopsis (Table S1). Phytol is released during chlorophyll degradation by a chloroplast-located pheophytinase [58]. The free phytol residue is redirected into chloroplast lipid metabolism [59] which would support an abundant plastidic pool (90±10%) as shown (Table S1).

Malate and fumarate were localized mainly in the cytosol and did not accumulate in the vacuole (Table S1) which is in contrast to reports indicating a large vacuolar pool [10]. However, in C3 plants malate accumulates during the day with a maximum at the end of the light period, only being transported into the vacuole after reaching a threshold concentration [60]. A further observation concerns the predominant aliphatic glucosinolate in Arabidopsis, 4-methylsulfinylbutyl glucosinolate (glucoraphanin), which revealed a dominant pool size within the vacuole (88±19%; Table S3). Glucoraphanin can be hydrolyzed by myrosinase into 5-methylsufinylpentylnitrile, which was dominantly localized in the cytosol (80±18%). In Arabidopsis myrosinase is localized in the vacuole of idioblastic cells of the phloem parenchyma [61], whereas glucosinolates are commonly reported to be stored in the vacuole, indicating that substrate and enzyme are likely not co-localized in the same cell [36], [61], [62]. The detection of glucoraphanin in the vacuole and the degradation product 5-methylsufinylpentylnitrile in the cytosol (Figure 7B) therefore provides evidence for the transport and catabolism of glucosinolates under physiological conditions that does not involve tissue disruption by herbivore attacks. Even though little is known about glucosinolate catabolism in plants, their concentrations can significantly vary in leaves during diurnal cycle [63], [64] or specific glucosinolates can be degraded during developmental processes [65].

Based on the within-compartment distance of the three markers, the plastidic compartment seemed to be resolved to a higher resolution than the others. Whereas starch is stored within the plastidial stroma [33], the galactolipids, MGDG and DGDG (Table S2), are found in both envelope and thylakoid membranes [34]. Surprisingly, NADP-GAPDH, an enzyme found within the stroma, is clearly deviant from both (Figure 4), and also showed a very similar distribution and close proximity to chlorophyll (Figure 7B; Data S6) and the light harvesting complex (Figure 1B). Studies in spinach [66] and Synechocystis [67] have provided circumstantial evidences that the Calvin cycle multienzyme complex seems to be bound to thylakoid membranes and thus may indicate a partial separation of the thylakoid and the stroma of plastids under our NAF conditions.

Recurring distribution patterns throughout the gradients suggests the existence and contribution of previously unconsidered compartments

As described above, the vast majority (81.5%) of analytes could be assigned to one of the five classes. In contrast, the subcellular distributions for another 724 analytes (18.5%) could not be sufficiently estimated as the compartment-specific markers did not encompass their distribution (Table 1; Data S4). These include aspartate, asparagine, glutamate, glutamine, serine as well as the mitochondrial marker citrate synthase amongst others. Considering that the distribution for these metabolites resembles, to some extent, the situation for the mitochondrial compartment, which also could not be unambiguously delineated, therefore we speculated that unresolved or unconsidered compartments may contribute to the recurring and unexplained distribution patterns. To test this hypothesis on the 724 analytes we tried to identify analyte groups characterized by similar, yet distant and reproducible distributions using a marker-‘free’-based classification by k-medoids clustering, allowing only the cytosolic compartment to be partitioned into different clusters without being assigned to another compartment.

This resulted in the identification of seven clusters of which two are represented by the cytosolic compartment (Figure S4A; Data S4). Out of the 724 analytes with insufficiently explained distributions 339 (46.8%) were assigned into one of the two cytosolic, 125 (17.3%) into the plastidic, and 26 (3.6%) into the vacuolar cluster. Further 234 (32.3%) analytes were assigned into three novel clusters (Figure 8). Mantel tests performed as non-parametric ANOVA revealed significant (P<0.001) and intermediate matrix correlations of r = 0.51±0.02 (as mean ± SD). The same approach, but restricted to the clusters reflecting the three considered compartments, resulted in a higher average matrix correlation of r = 0.82±0.01. Using gap static both biologically-driven cluster solutions are supported, however it indicated that the seven considered clusters are less well-separated (Figure S4C). Despite a certain degree of cluster overlap, line plots displayed robust intermediate cluster distributions in between the delineated subcellular compartments (Figure 8).

[Figure omitted. See PDF.]

Figure 8. Scatter and distribution plots of analytes with compartment-specific and unresolved subcellular distributions.

Analytes with compartment-specific distributions were identified using a classification tree based assignment (Figure 6). Analytes with insufficiently explained subcellular distribution were grouped according to k-medoids clustering (Figure S4A). For visualization, Manhattan distances among analytes for each of the three independent gradients were averaged and then converted into a principal coordinates (PCo) space. Analytes were color-coded according to their cluster membership. For each identified cluster the distribution of members throughout the gradients (grey lines) and the robust average distribution including standard deviations (black lines) are depicted as line plots. Despite the overlap of the intermediate clusters with the resolved compartments, recurring and stable distribution patterns have been observed.

https://doi.org/10.1371/journal.pone.0017806.g008

Specifically, the largest intermediate cluster ‘cpl-cyt’ (n = 94) showed similarity to the robust consensus distribution of the plastidic compartment, displaying increased abundance in fraction F5 and a decreased abundance in fraction F6 (Figure 8). Amongst others aspartate, glutamate, asparagine, and the mitochondrial marker citrate synthase are assigned into this cluster. Interestingly, it was shown that the glutamine synthase GLN2 is targeted between both chloroplasts and mitochondria and facilitates ammonium recovery by transferring ammonium to glutamate during photorespiration [68]. Aspartate aminotransferase activity in mitochondria indicates that aspartate, as its substrate, is also present [69]. Together with the mitochondrial marker citrate synthase this intermediate cluster may represent metabolites captured in transport between the plastids and mitochondria but as well as the cytosol, as serine, involved in photorespiration, is assigned into one of the cytosolic clusters (Data S4).

The cluster ‘vac-cyt’ comprises 80 mainly unknown analytes (Data S4) and has similarity to the robust consensus distribution of the vacuolar compartment. The abundances in the densest fractions F1 and F2 are similar, whereas for the vacuolar compartment the abundance in fraction F1 is about 2-fold higher compared to F2 (Figure 8).

The smallest cluster ‘cyt-vac’ comprises 60 members of which 59 are unknown secondary metabolites (Data S4). It strongly overlaps with the cytosolic and the ‘vac-cyt’ clusters, and shows the highest abundances in the fractions F2 and to a lesser extent F3 (Figure 8). Interestingly, most of these analytes are relatively large (average m/z 640) and have a relatively late retention time (53 with RTs greater than 14 min), indicating that these compounds could be very hydrophobic. At this point it might appear speculative to hypothesize about the provenience of these compounds since many reasonable explanations seem possible, still it is tempting to propose that this unusual cluster with specific distribution (as for the cluster ‘vac-cyt’) could be a derivative of the highly heterogeneous vacuole [70]. Another likely explanation could be that we are capturing some vesicles channeled between compartments [71], [72] or that we simply see an unconsidered compartment like the endoplasmatic reticulum (ER). The later would be supported by the structurally annotated metabolite 4-hydroxybenzoate, a precursor for the synthesis of the antimicrobial metabolite shikonin [73], [74] as well as an intermediate in ubiquinone biosynthesis. In both cases the biosynthetic reactions involving 4-hydroxybenzoate are localized in the ER and Golgi apparatus [75], [76] or in small vesicle derived from the ER [73], [77], [78]. A targeted proteomic or immunological approach towards the enzymes involved in these reactions might strengthen or dismiss this hypothesis.

Nevertheless, despite the identification of robust recurrent distribution patterns (Figure 8), the observed distributions are generally not distinctive enough when compared to the defined subcellular compartments. However, when this approach was applied on all analytes (Figures S4B), the intermediate cluster ‘cpl-cyt’ was supported (cf. Figures S4D), demonstrating that the observed intermediate distributions can be robustly identified. Even though a further subcellular compartment cannot be unambiguously delineated, the subcellular distributions of analytes with sufficiently explained distributions assigned into this cluster might be partially overestimated as this cluster comprises the mitochondrial marker and therefore metabolites shared between the mitochondria and plastids/cytosol (see above; Data S4).

Concluding remarks

By using an untargeted metabolic approach in combination with the development of an advanced method for critical analysis of NAF-derived metabolic data, we have gathered a comprehensive description of a compartmentalized (with regard to the cytosol (including the mitochondria), chloroplast, and vacuole) metabolome of an eukaryotic organism. The resultant comprehensive metabolic map of Arabidopsis leaves provides a resource that can serve as a basis to identify constraints and key processes as targets for biotechnology or for systems-biology driven research.

A precise understanding of how metabolites are synthesized, stored, and transported is critical for a better understanding of subcellular biochemical networks which will be important in biotechnological applications, as well as providing a basis to refine metabolic models by considering the subcellular localization of dominant pool sizes. This fact is of particular importance for plant energy metabolism which is closely linked with the plant plastid, mitochondria, and cytosol. In frame with this it will be of interest to sufficiently delineate not only the mitochondria from the cytosol but also to uncover novel subcellular distributions. While marker-‘free’ reconstructions showed the contribution of unconsidered compartments in our data, an unambiguous designation and biological description for these compartments could not be achieved as they are mainly comprised of structurally unknown analytes. Currently, this represents one of the main limitations in NAF studies, as even the subcellular localization of structurally identified (known) metabolites are often not described in literature and even then their localization might still be variable. Therefore it is clear that a comprehensive framework of markers needs to be established to align and assemble metabolites based on the measurement of known, unambiguously localizable molecules. For this purpose it will be necessary to include, along with the metabolic data, more protein analyses. These could be either provided using more antibody-based assays or by performing proteomic measurements on the gradient fractions. Nevertheless, having developed the presented metabolomics resource we have also laid the groundwork needed in order to perform and analyze more complex experiments, such as a time course or changing environmental conditions.

With the biological validation of the dataset, and the promise in the future to be able to name some of the unknowns, this topographical map can aid in the discovery of novel transporters, biosynthesis enzymes, and generate hypotheses for undiscovered pathways. As NAF and the whole metabolomics platform are applicable to any eukaryotic organism, the provided optimized protocol (Text S1) for Arabidopsis and statistical workflow should be adaptable to many other organisms.

Materials and Methods

Plant growth

All wild-type Arabidopsis thaliana Col-0 plants were grown on soil for two weeks under short day conditions (8 h light) before being transferred for three weeks to long day conditions (16 h light) with 140 µmol m⁻² s⁻¹ photon flux density and a temperature of 21°C at 50% relative humidity. A total of 4–8 g pooled plant leaf material from individual plants was harvested at the beginning of the light period (about 3 h after light switched on), snap-frozen in liquid nitrogen, and stored at −80°C until use.

Non-aqueous fractionation

For determination of subcellular metabolite levels, cellular compartments were separated using density gradient centrifugation under non-aqueous conditions according to the methods for leaf material [11] with optimized conditions [22] (Text S1). Frozen Arabidopsis leaf material was homogenized using a ball mill, pre-cooled in liquid N₂ to avoid thawing, instead of using a mortar as mortar-ground material was insufficiently filtered through a 20 µm nylon net (used instead of quartz wool (data not shown)). The gradient volume, composed of the non-polar solvents tetrachlorethylene/heptane, was increased from 12 to 28 mL using a much smaller linear density ρ from 1.43 g cm⁻³ to 1.62 g cm⁻³. Most of the sample material was focused in the middle fractions with exception of the plastidic compartment enriched within the top fractions (data not shown). By testing several centrifugation velocities and durations, equilibrium distribution was already achieved at 5,000 g and 50 min instead of 25,000 g and 180 min [cf. 11], shortening the exposure time of sample material to the non-aqueous solvents.

SDS–PAGE and Western blotting

SDS–PAGE and Western blotting were conducted as described [79]. Western blots were blocked with skimmed milk and probed with polyclonal primary antibody against the light harvesting complex (LHC) from Pisum sativum or the subunit E of the vacuolar type H⁺-ATPase (V-ATPase; Abcam plc, Cambridge, UK). Anti-rabbit horse radish peroxidase-conjugated secondary antibodies were used to detect primary antibodies. All blots were developed using ECL Western blotting kit (GE Healthcare, Munich, Germany).

Enzyme and metabolite assays

Enzyme assay extracts were prepared according to Geigenberger and Stitt [80]. NADP-dependent glyceraldehyde-3-phosphate dehydrogenase (GAPDH, EC 1.2.1.12) was measured as described by Stitt et al. [81]. Uridine diphosphate (UDP)-glucose-pyrophosphorylase (UGPase, EC 2.7.7.9) was assayed according to Zrenner et al. [82]. Citrate synthase (EC 2.3.3.1) activity was determined as described [83]. Chlorophyll was extracted twice with 80% (v/v) and once with 50% (v/v) hot ethanol (30 min, 95°C) and determined as outlined by Arnon [84]. Starch was measured from the remaining pellet of ethanolic extracts according to Hendriks et al. [85]. Nitrate was analyzed by enzymatic reaction as described [86].

Metabolite profiling

For GC-TOF/MS analyses, dried fraction aliquots were extracted with cold 10∶3∶1 (v/v/v) methanol:chloroform:water (MCW) solution and two extract aliquots (100 µL, 150 µL) were derivatized and analyzed as described [87] with m/z acquisition of 85–750. The established GC/MS protocol allows quantification of sugars, sugar alcohols, organic and amino acids, ascorbate and some lipophilic compounds [87], [88], [89].

For LC/MS analyses lipophilic and secondary metabolites were extracted from dried fraction aliquots with cold 2.5∶1∶1 (v/v/v) MCW solution under shaking and sonication. After phase separation, aliquots of the upper, aqueous phase and lower, organic phase were dried and resuspended in ddH₂0 (secondary metabolites) or 50∶20∶25 (v/v/v) isopropanol/hexane/water (lipids). Extraction and derivatization of individual soluble thiols (cystein, γ-glutamylcysteine, glutathione) were performed as described [90]. UPLC separation of soluble thiols, secondary, and lipophilic metabolites were performed on a Waters Acquity UPLC system (Waters, Mildford, MA, USA) equipped with a BEH C₁₈ (thiols), a HSS T3 C₁₈ (secondary metabolites), or a BEH C₈ (lipids) reversed phase column (Waters) coupled to a Fourier Transform Ion Cyclotron Resonance Mass Spectrometer (thiols) or an Exactive Orbitrap (secondary and lipophilic metabolites) (both Thermo Fisher Scientific, Bremen, Germany). Mass spectra were recorded in full scan, positive ion mode with m/z acquisition of 100–1500 and 200–600 using 25,000 and 50,000 ppm resolution for soluble thiols and secondary or lipophilic metabolites, respectively (Text S1).

MS data analyses

GC/MS data were processed and aligned as described [87] using a curated library of authentic standards and unknown Arabidopsis compounds comprising 1,032 unique spectral entries (Krall et al., in prep.). The aligned data with 413 found library entries, were evaluated and curated (Text S1). The filtered raw GC/MS data comprises 40 samples and 203 curated analytes with 1 (0.01%) missing value. All GC/MS data were expressed relative to U-¹³C-sorbitol and extract replicates averaged after TIC normalization (Data S1).

High-resolution MS data were aligned or peaks extracted using GeneData (v5.3.7, Basel, Swizerland) and Xcalibur (v2.06, Thermo). Aligned FT-MS data, comprising 16,262 and 53,785 time-m/z features (afterwards analytes) of lipophilic and secondary metabolites, were filtered for consistently found analytes (Text S1). These resultant peak lists were then searched against KEGG [91] and KNApSAcK [43] for secondary metabolites using an in-house developed database search tool (GoBioSpace, Hummel et al., unpublished) while the lipid data was searched against an in-house compiled lipid database (Giavalisco et al., submitted). These filtered and uncurated data were derived from 20 samples comprising 2,804 and 910 analytes with 1,125 (2%) and 457 (2.5%) missing values for lipophilic and secondary metabolites, respectively (Data S2–S3). These analytes were annotated onto three levels: unknown, if no database hit could be assigned; match if an unverified database hit was assigned; and known for orthogonally validated database hits. The validation of known metabolites does not include the use of authentic reference standards, but instead relies on previously described compounds for Arabidopsis, the use of validated fragmentation patterns, and mass shifts of ¹³C, ¹⁵N, and ³⁴S isotope labeled Arabidopsis thaliana samples (Giavalisco et al., submitted). In order to estimate the number of potential non-redundant analytes within the FT-MS data, a correlative approach similar as described [92] was conducted by defining time/similarity (T/S) clusters (Text S1).

The individual MS data were assembled into a joint data set including metabolites measured by targeted MS approaches (thiols) and metabolic assays (chlorophyll, starch) (Data S4).

Statistical analyses and visualization

All statistical analyses were performed if not otherwise stated according to Sokal and Rohlf [32] using R 2.9.1.

Metabolite data were normalized to adjust for sample amount variations using the total ion count within and among gradients (Text S1). Analyte abundances were expressed as percentage from total (scaled data). Missing values were imputed by principal component analyses (PCA) [93]. Outliers, extreme deviations from the respective fraction means, were detected by a boxplot approach and replaced with the corresponding fraction mean to promote extraction of biological relevant and robust information (Text S1). The processed, i.e. normalized, imputed, and outlier-removed data are provided as supplemental data (Data S1, S2, S3).

Robust consensus distributions throughout gradients were computed using Tukey's biweight. The t-test was performed two-sided with equal or unequal variance determined using F-test. P-values were adjusted by Benjamini-Hochberg correction (P_BH) [94] to control the false discovery rate. Mantel tests were performed as Pearson's matrix correlations (r) between distance matrices or as non-parametric ANOVA. HCA using average linkage clustering were performed on Euclidean (Eq. 1) or Manhattan distances. P-values for cluster nodes were computed with R's pvclust [95]. Classical multidimensional scaling (CMD; [96]) on normalized Manhattan distances (Eq. 3) among analytes was used to reflect distances as points in principal coordinates space. This approach was used to visualize and assess the proximity of a metabolite (or compartment) to the delineated compartments. Gap statistic was performed to estimate the number of clusters [97].

To estimate the robustness of downstream results, fractions were randomly assembled into a total of 729 (726 random +3 original combination) non-redundant artificial gradients and analyses were repeated.

Compartmental distribution and assignment

Subcellular metabolite distributions were computed using the BestFit command line tool (available upon request) by a three-compartmental distribution strategy utilizing the best fit (BFA) [11] and non-negative least square (NNLS) [50] algorithm. The abundances of all markers delineating the same compartment were averaged for each gradient separately prior to computation.

Analytes were assigned onto the three resolved subcellular compartments using a k-nearest neighbor (kNN) approach [98] with k = 3 nearest neighbors (estimated using cross-evaluation) on normalized Manhattan distances (Data S4). Refined compartmental assignments (Data S4) were performed using a classification tree based on observed subcellular distribution (Figure 6) and marker-‘free’ by means of robust k-medoids clustering (PAM, partitioning around medoids). The number of clusters (k) was determined by allowing only the cytosolic compartment (represented by three compartment-specific markers) to be partitioned into different clusters without being assigned onto another compartment. The validity of identified cluster numbers was evaluated using gap statistics. Non-parametric ANOVA by means of Mantel test was performed on 5 randomly selected cluster members for each cluster and repeated 999 times.

Equations

(Eq. 1) Euclidean distance d_E

RSS  =  residual sum of squares

(Eq. 2) Q-value

(Eq. 3) Normalized Manhattan distance d_m (on scaled data)(3a)(3b)

with (Manhattan distance)

Supporting Information

[Figure omitted. See PDF.]

Figure S1.

Box plots illustrating (A–C) the silhouette information and matrix correlation of assembled fraction group solutions as well as (D–F) the gap statistics to estimate the number of sample clusters on the basis of (A, D) primary, (B, E) lipophilic and (C, F) secondary metabolite data.

https://doi.org/10.1371/journal.pone.0017806.s001

(DOC)

Figure S2.

(A) Heatmap and cluster distribution of selected markers representing the three resolved subcellular compartments and (B) gap curves to estimate the number of marker clusters.

https://doi.org/10.1371/journal.pone.0017806.s002

(DOC)

Figure S3.

Venn diagrams of compartmental assignments of analytes separated according to the major compound classes (A) primary, (B) lipophilic, and (C) secondary metabolite data.

https://doi.org/10.1371/journal.pone.0017806.s003

(DOC)

Figure S4.

Graphical visualization of (A, B) classification results and (C, D) gap curves based on k-medoids clustering regarding (A, C) analytes with insufficiently explained (unexplained) subcellular distributions and (B, D) all analytes.

https://doi.org/10.1371/journal.pone.0017806.s004

(DOC)

Table S1.

Subcellular metabolite distribution and assignment results for selected major compound classes of primary metabolic compounds.

https://doi.org/10.1371/journal.pone.0017806.s005

(DOC)

Table S2.

Subcellular metabolite distribution and assignment results for selected major compound classes of lipophilic metabolic compounds.

https://doi.org/10.1371/journal.pone.0017806.s006

(DOC)

Table S3.

Subcellular metabolite distribution and assignment results for selected major compound classes of secondary metabolic compounds.

https://doi.org/10.1371/journal.pone.0017806.s007

(DOC)

Text S1.

Supplemental extended methods and detailed supplemental data description. Further details of non-aqueous fractionation, mass-spectrometry based metabolome analyses as well as data and statistical analyses are provided.

https://doi.org/10.1371/journal.pone.0017806.s008

(PDF)

Data S1.

Raw and processed GC-TOF/MS data of primary metabolites.

https://doi.org/10.1371/journal.pone.0017806.s009

(XLS)

Data S2.

Raw and processed UPLC-FT/MS data of lipophilic metabolites.

https://doi.org/10.1371/journal.pone.0017806.s010

(XLS)

Data S3.

Raw and processed UPLC-FT/MS data of secondary metabolites.

https://doi.org/10.1371/journal.pone.0017806.s011

(XLS)

Data S4.

Fused metabolome data set covering analyte annotations as well as results of estimated subcellular distributions and compartmental assignments.

https://doi.org/10.1371/journal.pone.0017806.s012

(XLS)

Data S5.

Distribution of measured and fitted fraction abundances of analytes across the gradient based on three independent gradient data.

https://doi.org/10.1371/journal.pone.0017806.s013

(PDF)

Data S6.

Scatter plots of analytes and compartment-specific markers in the principal coordinates space for visual assessment of subcellular location.

https://doi.org/10.1371/journal.pone.0017806.s014

(PDF)

Acknowledgments

We would like to thank Aenne Eckardt, Gudrun Wolter, and Antje Bolze for excellent technical assistance regarding GC/MS and LC/MS analyses. We acknowledge Sebastian Klie and Dr. Takayuki Tohge for their critical reading of this manuscript. Furthermore, the comments from the two anonymous referees are gratefully acknowledged.

Author Contributions

Conceived and designed the experiments: SK PG DS. Performed the experiments: SK PG LK M-CS DS. Analyzed the data: SK PG LK DS. Contributed reagents/materials/analysis tools: BU DB U-IF ARF LW. Wrote the paper: SK PG LK U-IF ARF LW DS. Guided the study: SK PG LK U-IF ARF LW DS.

Citation: Krueger S, Giavalisco P, Krall L, Steinhauser M-C, Büssis D, Usadel B, et al. (2011) A Topological Map of the Compartmentalized Arabidopsis thaliana Leaf Metabolome. PLoS ONE 6(3): e17806. https://doi.org/10.1371/journal.pone.0017806

References

1. Linka N, Weber AP (2010) Intracellular metabolite transporters in plants. Mol Plant 3: 21–53.N. LinkaAP Weber2010Intracellular metabolite transporters in plants.Mol Plant32153

2. Weber AP, Fischer K (2007) Making the connections–the crucial role of metabolite transporters at the interface between chloroplast and cytosol. FEBS Lett 581: 2215–2222.AP WeberK. Fischer2007Making the connections–the crucial role of metabolite transporters at the interface between chloroplast and cytosol.FEBS Lett58122152222

3. Martinoia E, Maeshima M, Neuhaus HE (2007) Vacuolar transporters and their essential role in plant metabolism. J Exp Bot 58: 83–102.E. MartinoiaM. MaeshimaHE Neuhaus2007Vacuolar transporters and their essential role in plant metabolism.J Exp Bot5883102

4. Stitt M, Fernie AR (2003) From measurements of metabolites to metabolomics: an ‘on the fly’ perspective illustrated by recent studies of carbon-nitrogen interactions. Curr Opin Biotechnol 14: 136–144.M. StittAR Fernie2003From measurements of metabolites to metabolomics: an ‘on the fly’ perspective illustrated by recent studies of carbon-nitrogen interactions.Curr Opin Biotechnol14136144

5. Deuschle K, Chaudhuri B, Okumoto S, Lager I, Lalonde S, et al. (2006) Rapid metabolism of glucose detected with FRET glucose nanosensors in epidermal cells and intact roots of Arabidopsis RNA-silencing mutants. Plant Cell 18: 2314–2325.K. DeuschleB. ChaudhuriS. OkumotoI. LagerS. Lalonde2006Rapid metabolism of glucose detected with FRET glucose nanosensors in epidermal cells and intact roots of Arabidopsis RNA-silencing mutants.Plant Cell1823142325

6. Gutscher M, Pauleau AL, Marty L, Brach T, Wabnitz GH, et al. (2008) Real-time imaging of the intracellular glutathione redox potential. Nat Methods 5: 553–559.M. GutscherAL PauleauL. MartyT. BrachGH Wabnitz2008Real-time imaging of the intracellular glutathione redox potential.Nat Methods5553559

7. Wiseman JM, Ifa DR, Zhu Y, Kissinger CB, Manicke NE, et al. (2008) Desorption electrospray ionization mass spectrometry: Imaging drugs and metabolites in tissues. Proc Natl Acad Sci U S A 105: 18120–18125.JM WisemanDR IfaY. ZhuCB KissingerNE Manicke2008Desorption electrospray ionization mass spectrometry: Imaging drugs and metabolites in tissues.Proc Natl Acad Sci U S A1051812018125

8. Aebersold R, Mann M (2003) Mass spectrometry-based proteomics. Nature 422: 198–207.R. AebersoldM. Mann2003Mass spectrometry-based proteomics.Nature422198207

9. Behrens M (1932) Untersuchungen an isolierten Zell- und Gewebsbestandteilen. I. Mitteilung: Isolierung von Zellkernen des Kalbsherzmuskels. Zeits physiol Chem 209: 59–74.M. Behrens1932Untersuchungen an isolierten Zell- und Gewebsbestandteilen. I. Mitteilung: Isolierung von Zellkernen des Kalbsherzmuskels.Zeits physiol Chem2095974

10. Gerhardt R, Heldt HW (1984) Measurement of Subcellular Metabolite Levels in Leaves by Fractionation of Freeze-Stopped Material in Nonaqueous Media. Plant Physiol 75: 542–547.R. GerhardtHW Heldt1984Measurement of Subcellular Metabolite Levels in Leaves by Fractionation of Freeze-Stopped Material in Nonaqueous Media.Plant Physiol75542547

11. Riens B, Lohaus G, Heineke D, Heldt HW (1991) Amino Acid and Sucrose Content Determined in the Cytosolic, Chloroplastic, and Vacuolar Compartments and in the Phloem Sap of Spinach Leaves. Plant Physiol 97: 227–233.B. RiensG. LohausD. HeinekeHW Heldt1991Amino Acid and Sucrose Content Determined in the Cytosolic, Chloroplastic, and Vacuolar Compartments and in the Phloem Sap of Spinach Leaves.Plant Physiol97227233

12. Elbers R, Heldt HW, Schmucker P, Soboll S, Wiese H (1974) Measurement of the ATP/ADP ratio in mitochondria and in the extramitochondrial compartment by fractionation of freeze-stopped liver tissue in non-aqueous media. Hoppe Seylers Z Physiol Chem 355: 378–393.R. ElbersHW HeldtP. SchmuckerS. SobollH. Wiese1974Measurement of the ATP/ADP ratio in mitochondria and in the extramitochondrial compartment by fractionation of freeze-stopped liver tissue in non-aqueous media.Hoppe Seylers Z Physiol Chem355378393

13. Soboll S, Akerboom TP, Schwenke WD, Haase R, Sies H (1980) Mitochondrial and cytosolic ATP/ADP ratios in isolated hepatocytes. A comparison of the digitonin method and the non-aqueous fractionation procedure. Biochem J 192: 951–954.S. SobollTP AkerboomWD SchwenkeR. HaaseH. Sies1980Mitochondrial and cytosolic ATP/ADP ratios in isolated hepatocytes. A comparison of the digitonin method and the non-aqueous fractionation procedure.Biochem J192951954

14. Rauch U, Schulze K, Witzenbichler B, Schultheiss HP (1994) Alteration of the cytosolic-mitochondrial distribution of high-energy phosphates during global myocardial ischemia may contribute to early contractile failure. Circ Res 75: 760–769.U. RauchK. SchulzeB. WitzenbichlerHP Schultheiss1994Alteration of the cytosolic-mitochondrial distribution of high-energy phosphates during global myocardial ischemia may contribute to early contractile failure.Circ Res75760769

15. Gerhardt R, Stitt M, Heldt HW (1987) Subcellular Metabolite Levels in Spinach Leaves: Regulation of Sucrose Synthesis during Diurnal Alterations in Photosynthetic Partitioning. Plant Physiol 83: 399–407.R. GerhardtM. StittHW Heldt1987Subcellular Metabolite Levels in Spinach Leaves: Regulation of Sucrose Synthesis during Diurnal Alterations in Photosynthetic Partitioning.Plant Physiol83399407

16. Winter H, Robinson DG, Heldt HW (1993) Subcellular volumes and metabolite concentrations in barley leaves. Planta 191: 180–190.H. WinterDG RobinsonHW Heldt1993Subcellular volumes and metabolite concentrations in barley leaves.Planta191180190

17. Fettke J, Eckermann N, Tiessen A, Geigenberger P, Steup M (2005) Identification, subcellular localization and biochemical characterization of water-soluble heteroglycans (SHG) in leaves of Arabidopsis thaliana L.: distinct SHG reside in the cytosol and in the apoplast. Plant J 43: 568–585.J. FettkeN. EckermannA. TiessenP. GeigenbergerM. Steup2005Identification, subcellular localization and biochemical characterization of water-soluble heteroglycans (SHG) in leaves of Arabidopsis thaliana L.: distinct SHG reside in the cytosol and in the apoplast.Plant J43568585

18. Heineke D, Riens B, Grosse H, Hoferichter P, Peter U, et al. (1991) Redox Transfer across the Inner Chloroplast Envelope Membrane. Plant Physiol 95: 1131–1137.D. HeinekeB. RiensH. GrosseP. HoferichterU. Peter1991Redox Transfer across the Inner Chloroplast Envelope Membrane.Plant Physiol9511311137

19. Shannon JC, Pien FM, Cao H, Liu KC (1998) Brittle-1, an adenylate translocator, facilitates transfer of extraplastidial synthesized ADP–glucose into amyloplasts of maize endosperms. Plant Physiol 117: 1235–1252.JC ShannonFM PienH. CaoKC Liu1998Brittle-1, an adenylate translocator, facilitates transfer of extraplastidial synthesized ADP–glucose into amyloplasts of maize endosperms.Plant Physiol11712351252

20. Farré EM, Tiessen A, Roessner U, Geigenberger P, Trethewey RN, et al. (2001) Analysis of the compartmentation of glycolytic intermediates, nucleotides, sugars, organic acids, amino acids, and sugar alcohols in potato tubers using a nonaqueous fractionation method. Plant Physiol 127: 685–700.EM FarréA. TiessenU. RoessnerP. GeigenbergerRN Trethewey2001Analysis of the compartmentation of glycolytic intermediates, nucleotides, sugars, organic acids, amino acids, and sugar alcohols in potato tubers using a nonaqueous fractionation method.Plant Physiol127685700

21. Yamada K, Norikoshi R, Suzuki K, Imanishi H, Ichimura K (2009) Determination of subcellular concentrations of soluble carbohydrates in rose petals during opening by nonaqueous fractionation method combined with infiltration-centrifugation method. Planta 230: 1115–1127.K. YamadaR. NorikoshiK. SuzukiH. ImanishiK. Ichimura2009Determination of subcellular concentrations of soluble carbohydrates in rose petals during opening by nonaqueous fractionation method combined with infiltration-centrifugation method.Planta23011151127

22. Krueger S, Niehl A, Lopez Martin MC, Steinhauser D, Donath A, et al. (2009) Analysis of cytosolic and plastidic serine acetyltransferase mutants and subcellular metabolite distributions suggests interplay of the cellular compartments for cysteine biosynthesis in Arabidopsis. Plant Cell Environ 32: 349–367.S. KruegerA. NiehlMC Lopez MartinD. SteinhauserA. Donath2009Analysis of cytosolic and plastidic serine acetyltransferase mutants and subcellular metabolite distributions suggests interplay of the cellular compartments for cysteine biosynthesis in Arabidopsis.Plant Cell Environ32349367

23. Pan Z, Raftery D (2007) Comparing and combining NMR spectroscopy and mass spectrometry in metabolomics. Anal Bioanal Chem 387: 525–527.Z. PanD. Raftery2007Comparing and combining NMR spectroscopy and mass spectrometry in metabolomics.Anal Bioanal Chem387525527

24. Fiehn O (2001) Combining genomics, metabolome analysis, and biochemical modelling to understand metabolic networks. Comp Funct Genomics 2: 155–168.O. Fiehn2001Combining genomics, metabolome analysis, and biochemical modelling to understand metabolic networks.Comp Funct Genomics2155168

25. Fernie AR, Trethewey RN, Krotzky AJ, Willmitzer L (2004) Metabolite profiling: from diagnostics to systems biology. Nat Rev Mol Cell Biol 5: 763–769.AR FernieRN TretheweyAJ KrotzkyL. Willmitzer2004Metabolite profiling: from diagnostics to systems biology.Nat Rev Mol Cell Biol5763769

26. Eisenreich W, Bacher A (2007) Advances of high-resolution NMR techniques in the structural and metabolic analysis of plant biochemistry. Phytochemistry 68: 2799–2815.W. EisenreichA. Bacher2007Advances of high-resolution NMR techniques in the structural and metabolic analysis of plant biochemistry.Phytochemistry6827992815

27. Dettmer K, Aronov PA, Hammock BD (2007) Mass spectrometry-based metabolomics. Mass Spectrom Rev 26: 51–78.K. DettmerPA AronovBD Hammock2007Mass spectrometry-based metabolomics.Mass Spectrom Rev265178

28. Farré EM, Fernie AR, Willmitzer L (2008) Analysis of subcellular metabolite levels of potato tubers (Solanum tuberosum) displaying alterations in cellular or extracellular sucrose metabolism. Metabolomics 4: 161–170.EM FarréAR FernieL. Willmitzer2008Analysis of subcellular metabolite levels of potato tubers (Solanum tuberosum) displaying alterations in cellular or extracellular sucrose metabolism.Metabolomics4161170

29. Benkeblia N, Shinano T, Osaki M (2007) Metabolite profiling and assessment of metabolome compartmentation of soybean leaves using non-aqueous fractionation and GC-MS analysis. Metabolomics 3: 297–305.N. BenkebliaT. ShinanoM. Osaki2007Metabolite profiling and assessment of metabolome compartmentation of soybean leaves using non-aqueous fractionation and GC-MS analysis.Metabolomics3297305

30. Oparka K, Viola R, Wright K, Prior D (1992) Sugar transport and metabolism in the potato tuber. In: Pollock CJ, Farrar JF, Gordon AJ, editors. Carbon Partitioning within and between Organisms. Oxford: BIOS Scientific Publishers. pp. 91–114.K. OparkaR. ViolaK. WrightD. Prior1992Sugar transport and metabolism in the potato tuber.CJ PollockJF FarrarAJ GordonCarbon Partitioning within and between OrganismsOxfordBIOS Scientific Publishers91114

31. Stitt M, Lilley RM, Heldt HW (1982) Adenine Nucleotide Levels in the Cytosol, Chloroplasts, and Mitochondria of Wheat Leaf Protoplasts. Plant Physiol 70: 971–977.M. StittRM LilleyHW Heldt1982Adenine Nucleotide Levels in the Cytosol, Chloroplasts, and Mitochondria of Wheat Leaf Protoplasts.Plant Physiol70971977

32. Sokal RR, Rohlf FJ (1995) Biometry: The principles and practice of statistics in biological research. New York: W.H. Freeman and Company. 337 p.RR SokalFJ Rohlf1995Biometry: The principles and practice of statistics in biological research.New YorkW.H. Freeman and Company337

33. Zeeman SC, Kossmann J, Smith AM (2010) Starch: its metabolism, evolution, and biotechnological modification in plants. Annu Rev Plant Biol 61: 209–234.SC ZeemanJ. KossmannAM Smith2010Starch: its metabolism, evolution, and biotechnological modification in plants.Annu Rev Plant Biol61209234

34. Dörmann P, Benning C (2002) Galactolipids rule in seed plants. Trends Plant Sci 7: 112–118.P. DörmannC. Benning2002Galactolipids rule in seed plants.Trends Plant Sci7112118

35. Marrs KA, Alfenito MR, Lloyd AM, Walbot V (1995) A glutathione S-transferase involved in vacuolar transfer encoded by the maize gene Bronze-2. Nature 375: 397–400.KA MarrsMR AlfenitoAM LloydV. Walbot1995A glutathione S-transferase involved in vacuolar transfer encoded by the maize gene Bronze-2.Nature375397400

36. Kelly PJ, Bones A, Rossiter JT (1998) Sub-cellular immunolocalization of the glucosinolate sinigrin in seedlings of Brassica juncea. Planta 206: 370–377.PJ KellyA. BonesJT Rossiter1998Sub-cellular immunolocalization of the glucosinolate sinigrin in seedlings of Brassica juncea.Planta206370377

37. Debeaujon I, Peeters AJ, Leon-Kloosterziel KM, Koornneef M (2001) The TRANSPARENT TESTA12 gene of Arabidopsis encodes a multidrug secondary transporter-like protein required for flavonoid sequestration in vacuoles of the seed coat endothelium. Plant Cell 13: 853–871.I. DebeaujonAJ PeetersKM Leon-KloosterzielM. Koornneef2001The TRANSPARENT TESTA12 gene of Arabidopsis encodes a multidrug secondary transporter-like protein required for flavonoid sequestration in vacuoles of the seed coat endothelium.Plant Cell13853871

38. Lambrix V, Reichelt M, Mitchell-Olds T, Kliebenstein DJ, Gershenzon J (2001) The Arabidopsis epithiospecifier protein promotes the hydrolysis of glucosinolates to nitriles and influences Trichoplusia ni herbivory. Plant Cell 13: 2793–2807.V. LambrixM. ReicheltT. Mitchell-OldsDJ KliebensteinJ. Gershenzon2001The Arabidopsis epithiospecifier protein promotes the hydrolysis of glucosinolates to nitriles and influences Trichoplusia ni herbivory.Plant Cell1327932807

39. Yazaki K (2005) Transporters of secondary metabolites. Curr Opin Plant Biol 8: 301–307.K. Yazaki2005Transporters of secondary metabolites.Curr Opin Plant Biol8301307

40. Burow M, Rice M, Hause B, Gershenzon J, Wittstock U (2007) Cell- and tissue-specific localization and regulation of the epithiospecifier protein in Arabidopsis thaliana. Plant Mol Biol 64: 173–185.M. BurowM. RiceB. HauseJ. GershenzonU. Wittstock2007Cell- and tissue-specific localization and regulation of the epithiospecifier protein in Arabidopsis thaliana.Plant Mol Biol64173185

41. Marinova K, Pourcel L, Weder B, Schwarz M, Barron D, et al. (2007) The Arabidopsis MATE transporter TT12 acts as a vacuolar flavonoid/H+ -antiporter active in proanthocyanidin-accumulating cells of the seed coat. Plant Cell 19: 2023–2038.K. MarinovaL. PourcelB. WederM. SchwarzD. Barron2007The Arabidopsis MATE transporter TT12 acts as a vacuolar flavonoid/H+ -antiporter active in proanthocyanidin-accumulating cells of the seed coat.Plant Cell1920232038

42. Zhao J, Dixon RA (2010) The ‘ins’ and ‘outs’ of flavonoid transport. Trends Plant Sci 15: 72–80.J. ZhaoRA Dixon2010The ‘ins’ and ‘outs’ of flavonoid transport.Trends Plant Sci157280

43. Shinbo Y, Nakamura Y, Altaf-Ul-Amin U, Asahi H, Kurokawa K, et al. (2006) KNApSAcK: A comprehensive species-metabolite relationship database. In: Saito K, Dixon RA, Willmitzer L, editors. Plant Metabolomics (Biotechnology in Agriculture and Forestry) Heidelberg. pp. 165–184.Y. ShinboY. NakamuraU. Altaf-Ul-AminH. AsahiK. Kurokawa2006KNApSAcK: A comprehensive species-metabolite relationship database.K. SaitoRA DixonL. WillmitzerPlant Metabolomics (Biotechnology in Agriculture and Forestry) Heidelberg165184

44. Nakabayashi R, Kusano M, Kobayashi M, Tohge T, Yonekura-Sakakibara K, et al. (2009) Metabolomics-oriented isolation and structure elucidation of 37 compounds including two anthocyanins from Arabidopsis thaliana. Phytochemistry 70: 1017–1029.R. NakabayashiM. KusanoM. KobayashiT. TohgeK. Yonekura-Sakakibara2009Metabolomics-oriented isolation and structure elucidation of 37 compounds including two anthocyanins from Arabidopsis thaliana.Phytochemistry7010171029

45. Lin W, Oliver DJ (2008) Role of triacylglycerols in leaves. Plant Sci 175: 233–237.W. LinDJ Oliver2008Role of triacylglycerols in leaves.Plant Sci175233237

46. Lynch DV, Dunn TM (2004) An introduction to plant sphingolipids and a review of recent advances in understanding their metabolism and function. New Phyt 161: 677–702.DV LynchTM Dunn2004An introduction to plant sphingolipids and a review of recent advances in understanding their metabolism and function.New Phyt161677702

47. Logan DC (2010) Mitochondrial fusion, division and positioning in plants. Biochem Soc Trans 38: 789–795.DC Logan2010Mitochondrial fusion, division and positioning in plants.Biochem Soc Trans38789795

48. Noctor G, De Paepe R, Foyer CH (2007) Mitochondrial redox biology and homeostasis in plants. Trends Plant Sci 12: 125–134.G. NoctorR. De PaepeCH Foyer2007Mitochondrial redox biology and homeostasis in plants.Trends Plant Sci12125134

49. Doniwa Y, Arimura S, Tsutsumi N (2007) Mitochondria use actin filaments as rails for fast translocation in Arabidopsis and tobacco cells. Plant Biotechnol J 24: 441–447.Y. DoniwaS. ArimuraN. Tsutsumi2007Mitochondria use actin filaments as rails for fast translocation in Arabidopsis and tobacco cells.Plant Biotechnol J24441447

50. Lawson CL, Hanson RJ (1995) Solving Least Squares Problems. Classics in Applied Mathematics. Philadelphia: SIAM. CL LawsonRJ Hanson1995Solving Least Squares Problems. Classics in Applied Mathematics.PhiladelphiaSIAM

51. Ohlrogge J, Browse J (1995) Lipid biosynthesis. Plant Cell 7: 957–970.J. OhlroggeJ. Browse1995Lipid biosynthesis.Plant Cell7957970

52. Ono E, Hatayama M, Isono Y, Sato T, Watanabe R, et al. (2006) Localization of a flavonoid biosynthetic polyphenol oxidase in vacuoles. Plant J 45: 133–143.E. OnoM. HatayamaY. IsonoT. SatoR. Watanabe2006Localization of a flavonoid biosynthetic polyphenol oxidase in vacuoles.Plant J45133143

53. Winter H, Lohaus G, Heldt HW (1992) Phloem Transport of Amino Acids in Relation to their Cytosolic Levels in Barley Leaves. Plant Physiol 99: 996–1004.H. WinterG. LohausHW Heldt1992Phloem Transport of Amino Acids in Relation to their Cytosolic Levels in Barley Leaves.Plant Physiol999961004

54. Lehmann S, Funck D, Szabados L, Rentsch D (2010) Proline metabolism and transport in plant development. Amino Acids. S. LehmannD. FunckL. SzabadosD. Rentsch2010Proline metabolism and transport in plant development.Amino Acids

55. Szekely G, Abraham E, Cseplo A, Rigo G, Zsigmond L, et al. (2008) Duplicated P5CS genes of Arabidopsis play distinct roles in stress regulation and developmental control of proline biosynthesis. Plant J 53: 11–28.G. SzekelyE. AbrahamA. CseploG. RigoL. Zsigmond2008Duplicated P5CS genes of Arabidopsis play distinct roles in stress regulation and developmental control of proline biosynthesis.Plant J531128

56. Aubert S, Choler P, Pratt J, Douzet R, Gout E, et al. (2004) Methyl-beta-D-glucopyranoside in higher plants: accumulation and intracellular localization in Geum montanum L. leaves and in model systems studied by 13C nuclear magnetic resonance. J Exp Bot 55: 2179–2189.S. AubertP. CholerJ. PrattR. DouzetE. Gout2004Methyl-beta-D-glucopyranoside in higher plants: accumulation and intracellular localization in Geum montanum L. leaves and in model systems studied by 13C nuclear magnetic resonance.J Exp Bot5521792189

57. Schneider S, Beyhl D, Hedrich R, Sauer N (2008) Functional and physiological characterization of Arabidopsis INOSITOL TRANSPORTER1, a novel tonoplast-localized transporter for myo-inositol. Plant Cell 20: 1073–1087.S. SchneiderD. BeyhlR. HedrichN. Sauer2008Functional and physiological characterization of Arabidopsis INOSITOL TRANSPORTER1, a novel tonoplast-localized transporter for myo-inositol.Plant Cell2010731087

58. Schelbert S, Aubry S, Burla B, Agne B, Kessler F, et al. (2009) Pheophytin pheophorbide hydrolase (pheophytinase) is involved in chlorophyll breakdown during leaf senescence in Arabidopsis. Plant Cell 21: 767–785.S. SchelbertS. AubryB. BurlaB. AgneF. Kessler2009Pheophytin pheophorbide hydrolase (pheophytinase) is involved in chlorophyll breakdown during leaf senescence in Arabidopsis.Plant Cell21767785

59. Ischebeck T, Zbierzak AM, Kanwischer M, Dormann P (2006) A salvage pathway for phytol metabolism in Arabidopsis. J Biol Chem 281: 2470–2477.T. IschebeckAM ZbierzakM. KanwischerP. Dormann2006A salvage pathway for phytol metabolism in Arabidopsis.J Biol Chem28124702477

60. Martinoia E, Rentsch D (1994) Malate compartmentation: responses to a complex metabolism. Annu Rev Plant Physiol Plant Mol Biol 45: 447–467.E. MartinoiaD. Rentsch1994Malate compartmentation: responses to a complex metabolism.Annu Rev Plant Physiol Plant Mol Biol45447467

61. Andreasson E, Bolt Jorgensen L, Hoglund AS, Rask L, Meijer J (2001) Different myrosinase and idioblast distribution in Arabidopsis and Brassica napus. Plant Physiol 127: 1750–1763.E. AndreassonL. Bolt JorgensenAS HoglundL. RaskJ. Meijer2001Different myrosinase and idioblast distribution in Arabidopsis and Brassica napus.Plant Physiol12717501763

62. Koroleva OA, Davies A, Deeken R, Thorpe MR, Tomos AD, et al. (2000) Identification of a new glucosinolate-rich cell type in Arabidopsis flower stalk. Plant Physiol 124: 599–608.OA KorolevaA. DaviesR. DeekenMR ThorpeAD Tomos2000Identification of a new glucosinolate-rich cell type in Arabidopsis flower stalk.Plant Physiol124599608

63. Rosa EAS, Heaney RK, Rego FC, Fenwick GR (1994) The Variation of Glucosinolate Concentration during a Single Day in Young Plants of Brassica oleracea var acephala and capitata. J Sci Food Agric 66: 457–463.EAS RosaRK HeaneyFC RegoGR Fenwick1994The Variation of Glucosinolate Concentration during a Single Day in Young Plants of Brassica oleracea var acephala and capitata.J Sci Food Agric66457463

64. Rosa EAS, Heaney RK, Portas CAM, Fenwick GR (1996) Changes in Glucosinolate Concentrations in Brassica Crops (B oleracea and B napus) Throughout Growing Seasons. J Sci Food Agric 71: 237–244.EAS RosaRK HeaneyCAM PortasGR Fenwick1996Changes in Glucosinolate Concentrations in Brassica Crops (B oleracea and B napus) Throughout Growing Seasons.J Sci Food Agric71237244

65. Brown PD, Tokuhisa JG, Reichelt M, Gershenzon J (2003) Variation of glucosinolate accumulation among different organs and developmental stages of Arabidopsis thaliana. Phytochemistry 62: 471–481.PD BrownJG TokuhisaM. ReicheltJ. Gershenzon2003Variation of glucosinolate accumulation among different organs and developmental stages of Arabidopsis thaliana.Phytochemistry62471481

66. Süss KH, Arkona C, Manteuffel R, Adler K (1993) Calvin cycle multienzyme complexes are bound to chloroplast thylakoid membranes of higher plants in situ. Proc Natl Acad Sci U S A 90: 5514–5518.KH SüssC. ArkonaR. ManteuffelK. Adler1993Calvin cycle multienzyme complexes are bound to chloroplast thylakoid membranes of higher plants in situ.Proc Natl Acad Sci U S A9055145518

67. Agarwal R, Ortleb S, Sainis JK, Melzer M (2009) Immunoelectron microscopy for locating calvin cycle enzymes in the thylakoids of synechocystis 6803. Mol Plant 2: 32–42.R. AgarwalS. OrtlebJK SainisM. Melzer2009Immunoelectron microscopy for locating calvin cycle enzymes in the thylakoids of synechocystis 6803.Mol Plant23242

68. Taira M, Valtersson U, Burkhardt B, Ludwig RA (2004) Arabidopsis thaliana GLN2-encoded glutamine synthetase is dual targeted to leaf mitochondria and chloroplasts. Plant Cell 16: 2048–2058.M. TairaU. ValterssonB. BurkhardtRA Ludwig2004Arabidopsis thaliana GLN2-encoded glutamine synthetase is dual targeted to leaf mitochondria and chloroplasts.Plant Cell1620482058

69. Schultz CJ, Coruzzi GM (1995) The aspartate aminotransferase gene family of Arabidopsis encodes isoenzymes localized to three distinct subcellular compartments. Plant J 7: 61–75.CJ SchultzGM Coruzzi1995The aspartate aminotransferase gene family of Arabidopsis encodes isoenzymes localized to three distinct subcellular compartments.Plant J76175

70. Paris N, Stanley CM, Jones RL, Rogers JC (1996) Plant cells contain two functionally distinct vacuolar compartments. Cell 85: 563–572.N. ParisCM StanleyRL JonesJC Rogers1996Plant cells contain two functionally distinct vacuolar compartments.Cell85563572

71. Chanda A, Roze LV, Kang S, Artymovich KA, Hicks GR, et al. (2009) A key role for vesicles in fungal secondary metabolism. Proc Natl Acad Sci U S A 106: 19533–19538.A. ChandaLV RozeS. KangKA ArtymovichGR Hicks2009A key role for vesicles in fungal secondary metabolism.Proc Natl Acad Sci U S A1061953319538

72. Echeverria E (2000) Vesicle-mediated solute transport between the vacuole and the plasma membrane. Plant Physiol 123: 1217–1226.E. Echeverria2000Vesicle-mediated solute transport between the vacuole and the plasma membrane.Plant Physiol12312171226

73. Schulze-Lefert P (2004) Knocking on the heaven's wall: pathogenesis of and resistance to biotrophic fungi at the cell wall. Curr Opin Plant Biol 7: 377–383.P. Schulze-Lefert2004Knocking on the heaven's wall: pathogenesis of and resistance to biotrophic fungi at the cell wall.Curr Opin Plant Biol7377383

74. Sircar D, Mitra A (2008) Evidence for p-hydroxybenzoate formation involving enzymatic phenylpropanoid side-chain cleavage in hairy roots of Daucus carota. J Plant Physiol 165: 407–414.D. SircarA. Mitra2008Evidence for p-hydroxybenzoate formation involving enzymatic phenylpropanoid side-chain cleavage in hairy roots of Daucus carota.J Plant Physiol165407414

75. Swiezewska E, Dallner G, Andersson B, Ernster L (1993) Biosynthesis of ubiquinone and plastoquinone in the endoplasmic reticulum-Golgi membranes of spinach leaves. J Biol Chem 268: 1494–1499.E. SwiezewskaG. DallnerB. AnderssonL. Ernster1993Biosynthesis of ubiquinone and plastoquinone in the endoplasmic reticulum-Golgi membranes of spinach leaves.J Biol Chem26814941499

76. Ohara K, Kokado Y, Yamamoto H, Sato F, Yazaki K (2004) Engineering of ubiquinone biosynthesis using the yeast coq2 gene confers oxidative stress tolerance in transgenic tobacco. Plant J 40: 734–743.K. OharaY. KokadoH. YamamotoF. SatoK. Yazaki2004Engineering of ubiquinone biosynthesis using the yeast coq2 gene confers oxidative stress tolerance in transgenic tobacco.Plant J40734743

77. Yamaga Y, Nakanishi K, Fukui H, Tabata M (1993) Intracellular localization of p-hydroxybenzoate geranyltransferase, a key enzyme involved in shikonin biosynthesis. Phytochemistry 32: 633–636.Y. YamagaK. NakanishiH. FukuiM. Tabata1993Intracellular localization of p-hydroxybenzoate geranyltransferase, a key enzyme involved in shikonin biosynthesis.Phytochemistry32633636

78. Yazaki K, Kunihisa M, Fujisaki T, Sato F (2002) Geranyl diphosphate:4-hydroxybenzoate geranyltransferase from Lithospermum erythrorhizon. Cloning and characterization of a ket enzyme in shikonin biosynthesis. J Biol Chem 277: 6240–6246.K. YazakiM. KunihisaT. FujisakiF. Sato2002Geranyl diphosphate:4-hydroxybenzoate geranyltransferase from Lithospermum erythrorhizon. Cloning and characterization of a ket enzyme in shikonin biosynthesis.J Biol Chem27762406246

79. Sambrook J, Fritsch EF, Maniatis T (1989) Molecular Cloning: A Laboratory Manual. New York: Cold Spring Harbor Laboratory Press. J. SambrookEF FritschT. Maniatis1989Molecular Cloning: A Laboratory Manual.New YorkCold Spring Harbor Laboratory Press

80. Geigenberger P, Stitt M (1993) Sucrose synthase catalyses a readily reversible reaction in vivo in developing potato tubers and other plant tissues. Planta 189: 329–339.P. GeigenbergerM. Stitt1993Sucrose synthase catalyses a readily reversible reaction in vivo in developing potato tubers and other plant tissues.Planta189329339

81. Stitt M, Wirtz W, Heldt HW (1983) Regulation of Sucrose Synthesis by Cytoplasmic Fructosebisphosphatase and Sucrose Phosphate Synthase during Photosynthesis in Varying Light and Carbon Dioxide. Plant Physiol 72: 767–774.M. StittW. WirtzHW Heldt1983Regulation of Sucrose Synthesis by Cytoplasmic Fructosebisphosphatase and Sucrose Phosphate Synthase during Photosynthesis in Varying Light and Carbon Dioxide.Plant Physiol72767774

82. Zrenner R, Willmitzer L, Sonnewald U (1993) Analysis of the expression of potato uridinediphosphate-glucose pyrophosphorylase and its inhibition by antisense RNA. Planta 190: 247–252.R. ZrennerL. WillmitzerU. Sonnewald1993Analysis of the expression of potato uridinediphosphate-glucose pyrophosphorylase and its inhibition by antisense RNA.Planta190247252

83. Bergmeyer HU (1987) Methods of Enzymatic Analysis. Weinheim: VCH. HU Bergmeyer1987Methods of Enzymatic Analysis.WeinheimVCH

84. Arnon DI, Hoagland DR (1939) A Comparison of Water Culture and Soil as Media for Crop Production. Science 89: 512–514.DI ArnonDR Hoagland1939A Comparison of Water Culture and Soil as Media for Crop Production.Science89512514

85. Hendriks JH, Kolbe A, Gibon Y, Stitt M, Geigenberger P (2003) ADP-glucose pyrophosphorylase is activated by posttranslational redox-modification in response to light and to sugars in leaves of Arabidopsis and other plant species. Plant Physiol 133: 838–849.JH HendriksA. KolbeY. GibonM. StittP. Geigenberger2003ADP-glucose pyrophosphorylase is activated by posttranslational redox-modification in response to light and to sugars in leaves of Arabidopsis and other plant species.Plant Physiol133838849

86. Cross JM, von Korff M, Altmann T, Bartzetko L, Sulpice R, et al. (2006) Variation of enzyme activities and metabolite levels in 24 Arabidopsis accessions growing in carbon-limited conditions. Plant Physiol 142: 1574–1588.JM CrossM. von KorffT. AltmannL. BartzetkoR. Sulpice2006Variation of enzyme activities and metabolite levels in 24 Arabidopsis accessions growing in carbon-limited conditions.Plant Physiol14215741588

87. Krall L, Huege J, Catchpole G, Steinhauser D, Willmitzer L (2009) Assessment of sampling strategies for gas chromatography-mass spectrometry (GC-MS) based metabolomics of cyanobacteria. J Chromatogr B Analyt Technol Biomed Life Sci 877: 2952–2960.L. KrallJ. HuegeG. CatchpoleD. SteinhauserL. Willmitzer2009Assessment of sampling strategies for gas chromatography-mass spectrometry (GC-MS) based metabolomics of cyanobacteria.J Chromatogr B Analyt Technol Biomed Life Sci87729522960

88. Roessner U, Luedemann A, Brust D, Fiehn O, Linke T, et al. (2001) Metabolic profiling allows comprehensive phenotyping of genetically or environmentally modified plant systems. Plant Cell 13: 11–29.U. RoessnerA. LuedemannD. BrustO. FiehnT. Linke2001Metabolic profiling allows comprehensive phenotyping of genetically or environmentally modified plant systems.Plant Cell131129

89. Lisec J, Schauer N, Kopka J, Willmitzer L, Fernie AR (2006) Gas chromatography mass spectrometry-based metabolite profiling in plants. Nat Protoc 1: 387–396.J. LisecN. SchauerJ. KopkaL. WillmitzerAR Fernie2006Gas chromatography mass spectrometry-based metabolite profiling in plants.Nat Protoc1387396

90. Hell R, Bergmann L (1990) γ-Glutamylcysteine synthetase in higher plants: catalytic properties and subcellular localization. Planta 180: 603–612.R. HellL. Bergmann1990γ-Glutamylcysteine synthetase in higher plants: catalytic properties and subcellular localization.Planta180603612

91. Kanehisa M, Goto S (2000) KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res 28: 27–30.M. KanehisaS. Goto2000KEGG: kyoto encyclopedia of genes and genomes.Nucleic Acids Res282730

92. Werner E, Croixmarie V, Umbdenstock T, Ezan E, Chaminade P, et al. (2008) Mass spectrometry-based metabolomics: accelerating the characterization of discriminating signals by combining statistical correlations and ultrahigh resolution. Anal Chem 80: 4918–4932.E. WernerV. CroixmarieT. UmbdenstockE. EzanP. Chaminade2008Mass spectrometry-based metabolomics: accelerating the characterization of discriminating signals by combining statistical correlations and ultrahigh resolution.Anal Chem8049184932

93. Stacklies W, Redestig H, Scholz M, Walther D, Selbig J (2007) pcaMethods–a bioconductor package providing PCA methods for incomplete data. Bioinformatics 23: 1164–1167.W. StackliesH. RedestigM. ScholzD. WaltherJ. Selbig2007pcaMethods–a bioconductor package providing PCA methods for incomplete data.Bioinformatics2311641167

94. Benjamini Y, Hochberg Y (1995) Controlling the False Discovery Rate - a Practical and Powerful Approach to Multiple Testing. Journal of the Royal Statistical Society Series B-Methodological 57: 289–300.Y. BenjaminiY. Hochberg1995Controlling the False Discovery Rate - a Practical and Powerful Approach to Multiple Testing.Journal of the Royal Statistical Society Series B-Methodological57289300

95. Suzuki R, Shimodaira H (2006) Pvclust: an R package for assessing the uncertainty in hierarchical clustering. Bioinformatics 22: 1540–1542.R. SuzukiH. Shimodaira2006Pvclust: an R package for assessing the uncertainty in hierarchical clustering.Bioinformatics2215401542

96. Cox TF, Cox MAA (1994) Multidimensional Scaling. Monographs on Statistics and Applied Probability. Boca Raton: Chapman and Hall/CRC. TF CoxMAA Cox1994Multidimensional Scaling. Monographs on Statistics and Applied Probability.Boca RatonChapman and Hall/CRC

97. Tibshirani R, Walther G, Hastie T (2001) Estimating the Number of Clusters in a Dataset via the Gap Statistic. J R Stat Soc Ser B 63: 411–423.R. TibshiraniG. WaltherT. Hastie2001Estimating the Number of Clusters in a Dataset via the Gap Statistic.J R Stat Soc Ser B63411423

98. Ripley BD (1996) Pattern Recognition and Neural Networks. Cambridge: Cambridge University Press. BD Ripley1996Pattern Recognition and Neural Networks.CambridgeCambridge University Press