From Models to Crop Species: Caveats and Solutions for Translational Metabolomics

Although plant metabolomics is largely carried out on Arabidopsis it is essentially genome-independent, and thus potentially applicable to a wide range of species. However, transfer between species, or even between different tissues of the same species, is not facile. This is because the reliability of protocols for harvesting, handling and analysis depends on the biological features and chemical composition of the plant tissue. In parallel with the diversification of model species it is important to establish good handling and analytic practice, in order to augment computational comparisons between tissues and species. Liquid chromatography–mass spectrometry (LC–MS)-based metabolomics is one of the powerful approaches for metabolite profiling. By using a combination of different extraction methods, separation columns, and ion detection, a very wide range of metabolites can be analyzed. However, its application requires careful attention to exclude potential pitfalls, including artifactual changes in metabolite levels during sample preparation under variations of light or temperature and analytic errors due to ion suppression. Here we provide case studies with two different LC–MS-based metabolomics platforms and four species (Arabidopsis thaliana, Chlamydomonas reinhardtii, Solanum lycopersicum, and Oryza sativa) that illustrate how such dangers can be detected and circumvented.

Initially, the use of metabolic profiling in plants, as indeed in all species, was restricted to diagnostic approaches in which the obtained profiles were used as markers for a range of biological conditions (Sauter et al., 1988;Meyer et al., 2007;Semel et al., 2007;Carmo-Silva et al., 2009;Scherling et al., 2009;Widodo Patterson et al., 2009). Although such studies remain highly important, particularly in medical research (Nicholson and Wilson, 2003;Griffin and Nicholls, 2006), more sophisticated uses of metabolic profiling have recently been developed, including identifying regulated enzymes and exploring the regulatory structure of pathways (Tiessen et al., 2002;Arrivault et al., 2009), searching for unexpected effects of genetic manipulation (Catchpole et al., 2005), screening wild species for beneficial chemical composition El-Lithy et al., 2005), gaining a more comprehensive view of metabolic regulation and as part of integrative analyses for the systemic response of environmental genetic perturbations (Hirai et al., 2004Fukushima et al., 2009;Sulpice et al., 2009;Trenkamp et al., 2009). In addition to these uses, metabolomics is proving to be a powerful tool for gene functional annotation in plants. There are now several examples of Arabidopsis genes that have been identified with the help of metabolomic approaches including MYB transcription factors (Hirai et al., 2007;Stracke et al., 2007), O-methyltransferase , glycosyltransferases (Tohge et al., 2005;Yonekura-Sakakibara et al., 2007, acyltransferases (Luo et al., 2007), UDP-rhamnose synthase (Yonekura-Sakakibara et al., 2008), and pyrophosphorylase (Okazaki et al., 2009) with the approach being equally effective in other species (Aharoni et al., 2000;Goossens et al., 2003;Achnine et al., 2005;Yamazaki et al., 2008).
One advantage that metabolomics has over transcriptomics [with the exception of next-generation sequencing tools, see Detlef Weigels recent review (Schneeberger and Weigel, 2011)] and proteomics is that it is essentially genome-independent (Stitt and Fernie, 2003) and as such can be applied to a species whose genome has not been sequenced as easily as those whose has. This "democratization" of biology allows in depth functional analyses of many species for which a complete and fully annotated genome is not yet available (Schneeberger and Weigel, 2011). Despite this fact caution needs to be taken when adopting a method set up for one tissue of one species to analyze another tissue of that species or even another species. This is especially so for metabolite profiling. The plant kingdom contains an incredibly rich chemical diversity (St-Pierre and De Luca, 2000). It is obvious that this chemical diversity poses a large challenge and stimulates research in developing new and increasingly powerful approaches to separate, detect, and identify metabolites. However, it also raises important challenges for experimental design, sample handling, and validation of analytic procedures. This is because tissue composition affects the reliability with which a particular metabolite can be reliably extracted and analyzed. This problem is particularly acute when using liquid chromatography-mass spectrometry (LC-MS) due to the so called ion suppression effects wherein the composition of the extract affects the efficiency of ionization of some of its constituent analytes . That said, a number of relatively simple control tests, in combination with the growing number of chemoinformatic tools for metabolomics (Tohge and Fernie, 2009;Bais et al., 2010;Carroll et al., 2010;Cottret et al., 2010;Xia and Wishart, 2010), should at least ameliorate this phenomenon and hence facilitate high-quality translational metabolomics.
Driven by the increasing diversification of plant research away from the principle model species A. thaliana we present here case studies in which methods developed for this species are assessed for use in determining metabolite levels either in the unicellular algae C. reinhardtii or in the crop species rice and tomato. For the former we assessed the analysis of primary metabolism using an LC-MS/MS method developed to deliver validated measurements of the levels of Calvin-Benson cycle intermediates, organic acids, nucleotide-sugars, and nucleotides in Arabidopsis rosettes (Arrivault et al., 2009). Given that information documenting the transfer of gas chromatography-mass spectrometry (GC-MS)based methods of analysis of primary metabolites has already been extensively supplied for potato and tomato (Roessner et al., 2001;Roessner-Tunali et al., 2003), we chose crop species to focus our studies on secondary metabolism. The two LC-MS-based methods applied in this study complement standard and well-established GC-MS methods by greatly increasing the range of metabolites that can be analyzed.
Here some examples of how can be performed using two different LC-MS-based metabolomics platforms, on one algal and two crop species and A. thaliana are shown. The combined results illustrate important experimental controls which should be implemented alongside computation algorithms in order to successfully adapt protocols that have been established for another biological system. This also applies to other LC-MS-based methods (Okazaki et al., 2009;Kanno et al., 2010). We additionally discuss how such studies could be used in conjuncture with novel tools for combined sequence comparison and co-expression analysis (Mutwil et al., 2011, and Ruprecht et al., this issue) in order to improve gene functional predictions from Arabidopsis to crop species.

CELL CULTURE AND EXTRACTION PROCEDURES
Chlamydomonas reinhardtii strain CC-1690 wild type mt+ was acquired from the Chlamydomonas Genetics Center (Duke University, Durham, NC, USA). Single colonies were used to inoculate the growth media containing 5 mM Hepes, 1 mM K-phosphates, Beijerinck salts (final concentrations of 7.5 mM NH 4 Cl, 0.34 mM CaCl 2 , 0.41 mM MgSO 4 ) and trace salt solution (final concentrations of 184 μM H 3 BO 3 , 77 μM ZnSO 4 , 26 μM MnCl 2 , 18 μM FeSO 4 , 7 μM CoCl 2 , 6 μM CuSO 4 , 1 μM (NH 4 ) 6 Mo 7 O 24 ; Harris, 1989;May et al., 2008;Kempa et al., 2009) at 25˚C under constant illumination with 400 μmol photon m −2 s −1 and continuous shaking. The amount of NH 4 Cl was reduced to 4 mM for the experiment shown in Figure 1 to reduce the impact of ion suppression. Before harvesting, cells were grown to a density of 3 × 10 6 cells ml −1 and dark-adapted for a minimum of 20 min before transferring 1 ml of cells to a cuvette and exposed them to 660 μmol photon m −2 s −1 under continuous stirring. Before illumination and at different time points after illumination the suspension was quenched by vigorously adding 2 ml of −70˚C methanol (70%). The entire mix was then lyophilized to dryness at −80˚C and extracted at 4˚C by a chloroform:methanol:water (1:2:5 [v/v]) mixture. Water fractions of three subsequent washes were collected, concentrated by lyophilization at −80˚C and filtered before metabolite measurement (80 μl of extracted culture in 100 μl sample measured) by ion pair (reverse-phase) chromatography triple quadrupole MS (IPC-MS/MS) detection.

PLANT MATERIALS AND EXTRACTION PROCEDURES
Arabidopsis thaliana ecotype Col-0 and S. lycopersicum (M82) were grown in soil in a controlled environmental chamber (16 h light/8 h dark photoperiod; 21˚C at 145 μmol photon m −2 s −1 and 25˚C at 500 μmol photon m −2 s −1 , respectively). A. thaliana used for recovery test was grown in general 1/2 MS agar plate in controlled plant growth chamber (16 h light/8 h dark photoperiod; 21˚C). Rice (O. sativa, Nipponbare) seeds were pre-germinated in tap water at 28˚C for 10 days. Plantlets were transferred to a controlled growth chamber with 12 h day length at 700 μmol photon m −2 s −1 . Plant material was frozen in liquid nitrogen, ground into powder and stored at −80˚C until use.

ION PAIR (REVERSE-PHASE) CHROMATOGRAPHY TRIPLE QUADRUPOLE MS (IPC-MS/MS)
Primary metabolite analysis by IPC-MS/MS was carried out on a Dionex HPLC system (Sunnyvale, CA, USA) coupled to a Finnigan TSQ Quantum Discovery MS-Q3 (Thermo Fisher Scientific, Waltham, USA) equipped with an electrospray ionization (ESI) interface. It was operated as described in Arrivault et al. (2009). Chromatographic separation was obtained at 35˚C by a multi-step gradient with online-degassed eluent A (10 mM tributylamine aqueous solution, adjusted to pH 4.95 with 15 mM acetic acid) and eluent B (methanol) applied to a Gemini (C18) 150 mm × 2.00 mm inner diameter 5 μm, 110 Å particle column (Phenomenex, Aschaffenburg, Germany). The MS-Q3 device was operated in the negative ion scanning mode with selected reaction monitoring (SRM). The MS-parameters for each compound are documented in the Supplementary Data of Arrivault et al. (2009). Calibration curves using authentic standards were used to calculate absolute amounts of metabolites in algal samples. Data were processed using LC-quan 2.5.6 SP1 software.

HARVESTING -OBTAINING REPRESENTATIVE MATERIAL AND AVOIDING HANDLING-INDUCED CHANGES
Expression of genes and activity of enzymes associated with photosynthesis, respiration, and energy metabolism are rapidly affected by changes in environmental conditions. Transcriptional and metabolic regulation by the circadian clock has been defined (Harmer et al., 2000;Gibon et al., 2006;Fukushima et al., 2009;Kerwin et al., 2011). Many metabolites showed marked diurnal changes. Problems related to variation in clock and diurnal rhythms can be circumvented by harvesting plants at the same time in the 24-h cycle. They can also be affected by shorter term fluctuations. The exact timing of harvesting and avoidance of perturbation of metabolism during harvesting, by for instance shading of leaves or changes in the oxygen tension are therefore critical (see Geigenberger et al., 2000). Additionally, rapid and complete quenching of metabolic activity is crucial to ensure faithful measurement of the intracellular metabolite content (discussed in details below).
To highlight the rapid metabolic response following light treatment, metabolite profiling was performed on the model organism C. reinhardtii (Figure 1). Quenching was performed by rapid mixing of the algal suspension in the light with an excess of very cold (−70˚C) methanol to instantaneously freeze the cells. Metabolite profiling by IPC-MS/MS analysis of short term light treatment was performed in dark-adapted material and after 0.25, 4.5, and 7.5 min illumination. Ribulose-1,5-bisphosphate (RuBP), which is a major metabolite in the Calvin-Benson cycle, was already significantly elevated by a light treatment of 15 s ( Figure 1A). The levels of malate, the late step in TCA cycle, also displayed significant increases upon illumination, whilst the level of the amino acid aspartate was not altered. Figure 1B illustrates why avoiding perturbations of the conditions by the mean of rapid harvesting is critical for metabolite profiling. Darkening significantly and almost instantaneously influences operation of the photosystems and the delivery of ATP and NADPH. The levels of RuBP, sedoheptulose-1,7-bisphosphate (SBP), fructose-1,6-bisphosphate (FBP), ADP-glucose (ADPG), and isocitrate were significantly decreased, while adenosine 5 -diphosphate (ADP) and adenosine monophosphate (AMP) increased within 2 s of darkening. Thus, harvesting protocols that lead to even very brief decrease or increase in the light intensity preceding quenching or during the quenching process will lead to erroneous estimates for the levels of metabolites. An identical problem arises in higher plants. This is due to the simple fact that the fluxes in the Calvin-Benson cycle are so high that many of the metabolites in the cycle as well as ATP and NADPH have short turnover times of 1 s or less (Arrivault et al., 2009).
This problem is of course especially critical for processes like photosynthesis, where fluxes are very fast and metabolite pools are small and turn over very quickly. However, it illustrates the more general points that (i) all available information about the turnover times of the metabolites-of-interest should be collected, evaluated, and used to design an appropriate harvesting and quenching protocol and (ii) that this protocol should be validated by checking if slowing down or speeding up the harvesting process modifies the levels of metabolites that are found in the harvested material.

QUENCHING OF ENZYMATIC ACTIVITIES AND DIFFERENCES BETWEEN PLANT SPECIES AND CHEMICAL PROPERTIES
Quenching of metabolic activity is not only essential to stop metabolic turnover in the running pathways, but also to inhibit other enzymatic activities that can destroy the metabolite after tissue disruption. An old but still instructive example of this is the precautions needed to determine pyrophosphate levels in plants (Weiner et al., 1987). Pyrophosphatase activity in leaves is so high that it can hydrolyse all the pyrophosphate in a leaf extract in <0.05 s. In an intact tissue, the vast majority of the pyrophosphatase activity is in the plastids whilst the pyrophosphate is in the cytosol. As soon as the tissue is disrupted the pyrophosphatase comes into contact with and destroys the pyrophosphate. To measure pyrophosphate it is therefore essential that enzymatic activity is completely stopped by rapid quenching and remains totally inactive during all subsequent stages in sample handlings as a fraction of a percent would be enough to destroy all the pyrophosphate within few seconds. Such problems can be routinely identified by recovery experiments in which representative amounts of authentic standards are added to the plant material before extraction, and it is checked that the added standard can be quantitatively detected in the final extract .
Whilst secondary metabolites do not display such rapid responses to changes in the environment as those observed for primary metabolism (see an example in Kusano et al., 2011), they, like primary metabolites, are highly susceptible to degradation by enzymes that come in contact with them after tissue disruption. For example glucosinolates are converted into isothiocyanates by myrosinase in Arabidopsis (Tierens et al., 2001;Barth and Jander, 2006). Given that degradative enzymes typically remain potent subsequent to freezing in liquid nitrogen when the extract is thawed, particular care must be taken during the extraction procedure. To illustrate this point, three extraction procedures were conducted in addition of our original extraction procedure (extraction a), using frozen powder of plant materials. To test if breakdown enzymes were definitively inactivated during this extraction procedure, extracts were incubated at 37˚C for 1 h (extraction b). A pre-incubation of sample material at 37˚C for 1 h prior extraction was conducted to test the extent of metabolite degradation upon thawing (extraction d). Secondary metabolite extraction is routinely performed at liquid nitrogen temperature, so to test if another temperature would affect metabolic composition, addition of buffer on frozen sample was performed at ice temperature (extraction c). All extractions were performed using pre-cooled extraction buffer (10˚C). Metabolite breakdown was assessed in A. thaliana leaves, O. sativa leaves, and S. lycopersicum fruits by mean of LC-MS. Total ion chromatograms and relative peak areas of selected metabolites are presented in Figures 2 and 3, respectively.  (2010), (b) extracts obtained by (a) method were incubated at 37˚C for 1 h, (c) extraction buffer was immediately added to frozen material on ice, (d) frozen sample was incubated at 37˚C for 1 h before extraction. All extractions were performed using pre-cooled extraction buffer (10˚C). Arrows show peaks which were newly or not detected in treated samples.

Frontiers in Plant Science | Plant Physiology
For almost all plant species, an incubation of the material after addition of extraction buffer at 37˚C for 1 h had no observable consequences on the total ion chromatograms and on the relative peak areas of selected metabolites [Figures 2 and 3, respectively. Comparison between (a) and (b)]. This result indicates that the major secondary metabolites are not broken down at 37˚C, provided the tissue has been taken up in the extraction buffer. However, metabolite profiling of samples extracted after 1 h (pre-incubation of the disrupted tissue at 37˚C) revealed that samples were significantly changed in some compound species (pointed by arrows in Figure 2). A more detailed analysis revealed that within a plant species and between various plant species the metabolite classes were differently affected [ Figure 3, comparison between (a) and (d)]. In general non-pigmented flavonoids such as flavonol glycoside in Arabidopsis leaves and S. lycopersicum fruits, and glycoflavone in O. sativa leaves were stable. By contrast, red-pigmented flavonoids namely anthocyanin derivatives in Arabidopsis were significantly decreased following the 37˚C preincubation. Furthermore, phenylpropanoids in Arabidopsis such as sinapoyl-derivatives were broken down by pre-incubation. The Brassica species specific secondary metabolites glucosinolates are well-known compounds which can be broken down by myrosinase (Tierens et al., 2001;Barth and Jander, 2006). The glucosinolate levels in extracts with 37˚C incubation before extraction were not detected [ Figure 3A, comparison between (a) and (d)]. Despite the breakdown of phenylpropanoids by enzyme activity in Arabidopsis, phenylpropanoids in tomato fruit such as chlorogenic acid related compounds were generally unaffected by enzymes. With the exception of esculeoside related compounds, levels of glycoalkaloids which are the major alkaloid in tomato fruits were not significantly changed ( Figure 3B). These data show that, whatever the plant species, a pre-incubation of the plant material at 37˚C prior to extraction leads generally to various levels of secondary metabolite breakdown due to the presence of active enzymes when plant material is thawed out. Addition of extraction buffer to frozen material not at liquid nitrogen temperature (extraction c) led to a significant decrease in the levels of glucosinolates and anthocyanin derivatives ( Figure 3A). This shows that the temperature during addition of the extraction buffer is also an important factor.
These results taken together illustrate that the tissue extraction should be carried out in the proper way with attention being taken to empirically optimize the extraction method for each and every new tissue measured. The effect of 20 min sonication was also evaluated in the same manner, but no differences were observed (data not shown). That said it is important to note that sonication should only be performed if this can be managed without an increase in temperature.

ION SUPPRESSION EFFECTS CAUSED BY GROWTH MEDIA
Although LC-MS analysis is a highly sensitive technique, ion suppression is a general problem of LC-MS analytical platforms due to altered ESI of a target ion by a contamination (Ikonomou et al., 1990;Kebarle and Tang, 1993;Buhrman et al., 1996;Matuszewski et al., 1998Matuszewski et al., , 2003King et al., 2000;Fernie et al., 2004). It is actually not a single event but a range of response-reducing phenomena which should be avoided as much as possible. While there is, however, no universal solution to this problem, understanding difference between samples and assessing the effects of ion suppression affords greater confidence in the accuracy of the results.
An example of ion suppression caused by growth conditions of Chlamydomonas is shown in Figure 4; Table 1. This example illustrates how the composition of the growth media can dramatically affect the reliability of metabolite analyses in Chlamydomonas. Following quenching of the algal suspension by mixing with an excess of cold methanol (see above), we took the entire suspension for analysis. This was necessary because some metabolites leak www.frontiersin.org FIGURE 3 | Comparison of secondary metabolite levels between different extraction methods. Relative peak area of the result described in Figure 2 of (A) Arabidopsis leaves, (B) tomato fruits, and (C) rice leaves, are shown. Y -axis indicates percentage calculated by average of peak area of treatment (a). (a-d) indicate different extraction methods as described in Figure 2 and Section "Materials and Methods." Blue related colors indicate flavonoid related compounds (flavonol and glycoflavone). Red color indicates red-pigmented flavonoid related compounds (anthocyanin). Yellow, phenylpropanoid related compounds (sinapoyl-derivatives and chlorogenate-related compound); green, glucosinolate related compounds in Arabidopsis (methionine-derived aliphatic glucosinolate and tryptophan-related aromatic glucosinolate); purple, indicates glycoalkaloid; gray, putative sterol derivative; black, unknown compounds. The m/z value detected in negative ion detection is indicated between parentheses. Average of three experimental replicates ±SD. Significant differences by t -test (p < 0.05) are marked by asterisks.
out of the cells into the methanol-water mix and are therefore lost when the quenched cells are harvested by centrifugation (data not shown, see also Krall et al., 2009). This means that metabolites from the cells must be analyzed in a matrix that contains methanol and all the components of the suspension medium. Because the Chlamydomonas cells are quite diluted, components of the growth

FIGURE 4 | Example of ion suppression caused by growth media. (A-C)
A standard mix containing all measured metabolites was mixed with individual components of the growth media to identify which component(s) of the growth media most severely influence ionization during electrospray ionization (ESI). (A) The standard mix was mixed with the whole growth medium (dark gray) and independently with the medium components Hepes, K-phosphates, trace salts, and Beijerinck salts (light gray). (B) In a second experiment, the individual components of the Beijerinck salts were tested for ion suppression and enhancement. Standards were mixed to all Beijerinck salts (dark gray) and independently to its components CaCl 2 , MgSO 4 , and NH 4 Cl (light gray). (C) A third experiment investigated if the anion or the cation was responsible for ion suppression caused by ammonium chloride. Standards were mixed to NH 4 Cl (dark gray), NH 4 HCO 3 , or NaCl (light gray). The data is shown as box plots of the average values (calculated from three technical replicates) for 24 metabolites. Significant outliers (p < 0.05) are identified in the figure panels. For a complete list of percentage of ion suppression and enhancement for all metabolites see Table A1 in Appendix.
medium are present in rather large amounts compared to metabolites in the cells. Unfortunately, some components of the growth medium lead to ion suppression.
We were alerted to ion suppression by three routine checks. First, comparison of the spectrum of metabolites with those expected from earlier studies of metabolites during photosynthesis in Arabidopsis showed low levels of several metabolites. Second, we checked whether the signal for each metabolite shows a strictly linear relationship to the amount of extract applied. In the case of ion suppression, the estimated levels of many metabolites decreased www.frontiersin.org

CC1690 cells grown in growth media and quenched by methanol (column V).
Values represent recovery ratios for all standards (n = 3, SD of raw data <18%).
when more samples were applied. Third, we checked the recovery of authentic standards added to the extract and found it was very low. Attempts to analyze high concentrations of sample resulted in an almost complete suppression of all metabolite signals, including those of spiked standards (Figure A1 in Appendix). A five-time dilution allowed an average of 87-93% recovery of the spiked standards (Table 1; Figure 4). However, there was still a residual ion suppression, and this varied from metabolite to metabolite (Table 1). This obviously still prevents reliable and comparable analysis of metabolite levels. We therefore carried out a further series of experiments to identify the major sources of ion suppression, in order to modify the growth medium and circumvent this problem. To show that ion suppression was caused by components of the growth medium and not the biological sample itself or the methanol, we first mixed known amounts of standards with either methanol and/or the growth media compared to the same known amounts of standards mixed with a sample containing Chlamydomonas cells (Table 1). For more than half of the metabolites, <15% of the signal was suppressed by the media, but for many other metabolites including dihydroxyacetone-phosphate (DHAP), fructose-6phosphate (F6P), glycerate, NAD + , ribose-5-phosphate (R5P), sedoheptulose-7-phosphate (S7P), UDP-glucose (UDPG), and xylulose-5-phosphate/ribulose-5-phosphate (X5P/Ru5P) the signals were reduced by 15-50%. For G6P ion suppression caused a decrease of intensities by >50%. Therefore, for all these metabolites the absolute values have to be treated with extreme caution. Further, small changes in the extent of ion suppression can lead to changes in the relative signals for the various metabolites. It almost goes without saying that mixing ISs of these metabolites to each sample would allow a much more precise determination of their absolute amounts.
To minimize such errors, we systematically investigated which of the salts in the medium could contribute to the loss of signal due to ion suppression or ion enhancement during ionization by ESI. The growth media used in this study consisted of Hepes, K-phosphates, Beijerinck salts, and trace salts (for details see Materials and Methods). A sequence of experiments was performed to unravel the effects of the individual salts from this growth media (Figure 4, for details see Table A1 in Appendix). Hepes, K-phosphate, and the trace salt solution had only minor ion suppression effects ( Figure 4A, for details see Table A1 in Appendix). Hepes caused significant ion suppression of G6P but less than the Beijerinck salts. K-phosphates caused weak but significant ion induction of ADP, and FBP, RuBP, SBP, aconitate, and isocitrate. The trace salt solution caused overestimation of aspartate due to ion enhancement. However, most of the ion suppression observed due to the growth media in the sample could be attributed to the Beijerinck salts present in the media (Figure 4A). In a second experiment, ammonium chloride, one component of the Beijerinck salts, was shown to be responsible for the major part of the residual ion suppression whereas MgSO 4 and CaCl 2 had minor effects (Figure 4B). In a third experiment, the chloride anion was found to be the main reason why ammonium chloride causes ion suppression ( Figure 4C). From the 24 metabolites routinely measured with this method, the signal of 10 was significantly suppressed in the presence of the Beijerinck salt (Figure A1 in Appendix). With the exception of UDPG which was found to be suppressed by MgSO 4 , the chloride anion in the growth media was found to be responsible for ion suppression. Thus, simply replacing the chloride anion with bicarbonate greatly decreased ion suppression ( Figure 4C).
For subsequent measurements a growth medium with lower ammonium chloride concentration was used (Figure 1). Alternatively, a medium in which ammonium chloride is replaced with ammonium bicarbonate could be used or, as mentioned earlier ISs for each metabolite could be added to the sample to assure accurate metabolite measurements. This example shows that different degrees of ion suppression, both with respect of individual metabolites and with respect to different samples, can be generated by differences in growth conditions and culture components. These results imply that erroneous results will also be obtained if such changes occur as a result of growing algae in different conditions, or are generated in time as a result of the algae using nutrients.

Frontiers in Plant Science | Plant Physiology
The specific issue with growth medium components does not arise with higher plants. Nevertheless, this example serves as a warning that differential ion accumulation in plant tissues might affect ion suppression. In attempt to circumvent this problem the total ion chromatogram should be carefully checked in areas which appear to be strongly affected. If strong ion suppression is observed, both dilution and recovery tests in which standard compounds are added to the extracts should be performed (see Fernie and Keurentjes, 2011). More generally, such problems can be identified by routine checks that the signal is linear with the amount of applied extract and (where available) that authentic samples can be quantitatively recovered after addition to the extract.

ION SUPPRESSION EFFECTS CAUSED BY DIFFERENT TISSUE TYPES
Metabolite composition varies between plant species and also between different tissues of the same plant. It is therefore expected to observe different levels of ion suppression within these samples. To evaluate this, an IS mixture (sinigrin, isovitexin, and saponarin) was added to the same volume (ratio 1/1) of extracts from Arabidopsis leaves, tomato fruits, and rice leaves, respectively. This was performed with four different known concentrations of the IS mixture with three experimental replicates of the step of mixing solutions. As control, the 50% diluted original standard mixture was also analyzed without being mixed with plant extract. Peak areas for each IS were determined and are presented in Figure 5. For all IS compounds, the strongest ion suppression was observed when they were added to rice leaves, followed by Arabidopsis leaves, and the lowest was seen for tomato fruit extracts. As expected, the ion suppression in a mixture of Arabidopsis leaf and tomato fruit extracts (ratio 1/1) was intermediate to what was observed in the corresponding independent extracts.
Plant material used for metabolic determination is often a mixture of tissues. For example, plant seedling is a mixture of hypocotyl and root, or fruit samples are a mixture of pericarp, seed, and peel. Ion suppression caused by differences between tissue types, in comparative analysis between mutants, transgenic, time course, and stress treatment is relatively minor. But in case of comparison between different plant species (as shown above with Arabidopsis and rice leaves), wild accessions which have important phenotypical differences, or mutants with strong phenotypic differences, such problems due to differential ion suppression could easily arise. For example, seedling samples that differ in the relative amount of shoot/root or hypocotyl/root might be particularly susceptible to differential ion suppression. The same problem is raised in the case of harvesting flower samples with a varying ratio of flower/sepal and/or (flower/pedicel).
To evaluate this problem, recombination analyses (Fernie and Keurentjes, 2011) were evaluated using Arabidopsis samples containing various ratio of leaves and roots or flower extracts, focusing on the general secondary metabolites which were detected in both tissues. The percentage recovery was simply estimated for evaluation using theoretical concentration of extracts mixture, [(level in leaves × A%) + (level in roots (or flowers) × B%)]/100] ( Table 2). Mixture of different tissue types results in >100% recovery (i.e., less ion suppression) for some peaks (e.g., IS), presumably because the area which is strongly suppressed differs between leaf and root extracts. Increase of recovery was observed in some peaks

FIGURE 5 | Example of ion suppression caused by plant extracts.
Standard solutions (20, 10, 5, 1 μg ml −1 ) of three standard compounds, (A) sinigrin, (B) saponarin, and (C) isovitexin, were mixed to the same volume of plant extracts (0.2 mg FW μl −1 ) of Arabidopsis leaves, tomato fruits, rice leaves and mixture of Arabidopsis leaves and tomato fruits (ratio 1/1). Average and ±SD were calculated from three experimental replicates.
(<117%) at an equivalent mixture of leaf and root extracts. On the other hand, recovery ratio in mixture of leaves and flowers showed more significant variance (43 ∼ 122%). This experiment highlights the value of preliminary analyses in order to check for ion suppression. It is furthermore useful to understand the range of variance in instances in which the value of IS is unstable between samples.

CHEMICAL DIVERSITY AND PEAK ANNOTATION USING CROSS SPECIES COMPARISON
Many metabolite databases for LC-MS are available, such as MASSBANK (Horai et al., 2010), METLIN , MS2T , KNApSAcK (Shinbo et al., 2006), and Flavonoid Viewer (Arita and Suwa, 2008). These greatly aid in the prediction and annotation of detected peaks . That said, technical improvement of peak identification and annotation still represents a major hurdle for LC-MS-based www.frontiersin.org metabolite profiling platforms. The identification of secondary metabolites is obstructed by the insufficient availability of standard compounds. It is impossible to comprehensively purchase standard substances since the diversity of their chemical structure is far too large. Moreover, complex compounds are largely unavailable commercially and those that are available are often prohibitively expensive. Furthermore, LC-MS studies are complicated by the fact that the levels of secondary metabolites are highly divergent between different organs, growth conditions, and species (Petersen, 2007;Hanhineva et al., 2008;Matsuda et al., 2010). For these reasons, peak identification is generally performed by the use of combinatorial strategies whereby the literature information is taken alongside available compounds in an attempt to identify specific peaks (see for example Tohge et al., 2005;Giavalisco et al., 2009). Given the recent explosion of genome information afforded firstly by microarray analyses and more recently by nextgeneration sequencing (review of (Schneeberger and Weigel, 2011), further tools for translational biology are becoming available. One such example, PlaNet, was described recently by Mutwil et al. (2011). Following this approach gene sequences can be connected between plant species on the basis of BLAST homology searches and then the positions in co-expression networks can be ascertained and finally it is possible to link unknown genes to annotated metabolic genes. As such this approach holds great promise both for gene functional annotation and via use of mutant plants in the annotation of unknown metabolites Mutwil et al., 2011). It demonstrated the utility of this approach by identifying candidate genes of the general and species specific flavonoid pathways. It is likely that integrating metabolomics data on all the species currently in PlaNet will greatly aid this process and is certainly a research avenue that should be pursued in the near future.

CONCLUSION
Whilst applying a method established for another species is likely not to be overly problematic for screening purposes and for a first insight into the metabolome of an organism, the examples presented here demonstrate that when more precise information is required considerable effort should be put into establishing both the qualitative and quantitative reliability of any LC-MS-based metabolic profiling method. As evidenced by the ion suppression (and ion enhancement) examples particular care must be taken with this issue as well as in ensuring that the extraction procedure is appropriate for the tissue under study. Once these important controls have been adhered to a wide array of computational resources are available (Tohge and Fernie, 2009), which will greatly aid in translational research. Given that the trend Frontiers in Plant Science | Plant Physiology in plant science research is to move away from the single model species of A. thaliana, such tools will become increasingly important. However, it is prudent to note that uncritical use of such tools without adequate controls of the type demonstrated here may well result in inaccurate representations of the metabolome. The best way to approach a new tissue, species, or even a dramatic mutant/transgenic line is to adopt both experimental and computational approaches to ensure the highest possible data quality.

ACKNOWLEDGMENTS
We thank Dr. Yozo Okazaki in RIKEN PSC for useful discussion, and Dr. Mark-Aurel Schöttler at the Max-Planck-Institute of Molecular Plant Physiology (MPIMP) for kindly providing instrumentation for Chlamydomonas illumination and harvesting. We thank Prof. Dr. Martin Steup in University of Potsdam and Dr. Wagner L. Araújo at MPIMP for expert comments. This work is partially funded by the German Federal Ministry of Education and Research by the FORSYS BMBF grant (GoFORSYS) and the Max-Planck Society (MPG). Research activity of Takayuki Tohge was supported by the Alexander von Humboldt Foundation. Tabea Mettler was supported by the International Max-Planck Research School (IMPRS). Adam James Carroll was supported through a grant to the Australian Research Council Centre of Excellence in Plant Energy Biology. FIGURE A1 | Sample dilution led to decreased ion suppression. Illuminated Chlamydomonas reinhardtii CC-1690 suspension was quenched, harvested, and extracted as in Figure 1. Four different dilutions of the extract (diluted with water) were then spiked with a reference standard mix and compared to the standard mix on its own to assess the extent of ion suppression. The fivefold dilution was selected for experiments shown in Figures 1 and 4; Table 1 as for several metabolites the 10-fold dilution was below detection limit. Note that the average of ion suppression was 77% for the fivefold dilution in this experiment. For experiment shown in Figures 1 and 4; Table 1, the average ion suppression for the fivefold dilution was <13%. This improvement was vey likely achieved by using a quadrupole (Finning TSQ Quantum Ultra) with a larger ion transfer tube diameter. The plot shows the mean ion suppression for 24 metabolites, each measured with two technical replicates, ±SD. www.frontiersin.org