==== Front Int J Mol Sci Int J Mol Sci ijms International Journal of Molecular Sciences 1422-0067 MDPI 33256100 10.3390/ijms21238987 ijms-21-08987 Review Plant Single-Cell Metabolomics—Challenges and Perspectives https://orcid.org/0000-0002-7200-8808de Souza Leonardo Perez 1* https://orcid.org/0000-0003-1359-7611Borghi Monica 2 Fernie Alisdair 1* 1 Max Planck Institute of Molecular Plant Physiology, Am Müehlenberg 1, Golm, 14476 Potsdam, Germany 2 Department of Biology, Utah State University, 1435 Old Main Hill, Logan, UT 84322, USA; monica.borghi@usu.edu * Correspondence: lperez@mpimp-golm.mpg.de (L.P.d.S.); fernie@mpimp-golm.mpg.de (A.F.) 26 11 2020 12 2020 21 23 898702 11 2020 25 11 2020 © 2020 by the authors.2020Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).Omics approaches for investigating biological systems were introduced in the mid-1990s and quickly consolidated to become a fundamental pillar of modern biology. The idea of measuring the whole complement of genes, transcripts, proteins, and metabolites has since become widespread and routinely adopted in the pursuit of an infinity of scientific questions. Incremental improvements over technical aspects such as sampling, sensitivity, cost, and throughput pushed even further the boundaries of what these techniques can achieve. In this context, single-cell genomics and transcriptomics quickly became a well-established tool to answer fundamental questions challenging to assess at a whole tissue level. Following a similar trend as the original development of these techniques, proteomics alternatives for single-cell exploration have become more accessible and reliable, whilst metabolomics lag behind the rest. This review summarizes state-of-the-art technologies for spatially resolved metabolomics analysis, as well as the challenges hindering the achievement of sensu stricto metabolome coverage at the single-cell level. Furthermore, we discuss several essential contributions to understanding plant single-cell metabolism, finishing with our opinion on near-future developments and relevant scientific questions that will hopefully be tackled by incorporating these new exciting technologies. cell type specific metabolismmetabolomicssingle-cellmass spectrometry imaging ==== Body 1. Introduction The advent of genomics immediately followed by similar conceptual frameworks to investigate transcriptomes, proteomes, and metabolomes represented a paradigm shift in biological systems investigation. The appealing idea of holistically assessing such systems has translated into rapid developments for systems biology. Researchers can now investigate multiple processes simultaneously, revealing essential mechanisms involved in regulating development and responses to the environment. For practical reasons, such techniques have been mainly applied to bulk samples consisting of a large number of cells for which results correspond to populations’ averages (Figure 1). In such experiments, the stochasticity of biological processes leading to cell heterogeneity is often considered not to be biologically relevant. Indeed, this is often the case, and for many applications such as characterizing mutants of central metabolic pathways [1,2] and identifying genes involved in the production of specialized metabolites [3,4,5], the use of averages is undoubtedly suitable. However, cell heterogeneity has been shown to play important biological roles in many situations for which averaging would mask relevant mechanistic insights [6]. In plants, several works highlighted the importance of cell-specific metabolism in regulating essential physiological processes such as the metabolism of the shoot apical meristem [7], the regulation of stomatal closure by guard cells and subsidiary cells [8,9], C4 metabolism [10,11,12], and the evolution of specialized metabolism [13]. However, most of these studies involve cell-specific labor-intensive protocols for cell isolation or reporter lines targeting few metabolites. True metabolomics at the cellular level remains a daunting task due to innumerable challenges in measuring metabolites. 2. Technical Challenges Current coverage of the metabolome still lags far behind genomics, transcriptomics, and proteomics because of the technical limitations imposed by the nature of metabolites. DNA, RNA, and proteins exhibit high regularity as they are constituted by a set of repeating unities, namely nucleotides and amino acids. As a consequence, these classes of molecules have characteristic physicochemical properties that are similar between them. Metabolites, on the other hand, exhibit much broader physicochemical diversity hindering their global analysis by a single technique. The broadest coverages of the metabolome achievable by date rely heavily upon the high sensitivity of mass spectrometry techniques hyphenated to efficient separation provided by gas and liquid chromatography. Although current technological advances provide considerable resolution in benchmark equipment such as orbitraps and quadrupole time-of-flight mass spectrometers (QTOFs), the combination of these two techniques is still essential for overcoming matrix effects providing maximum metabolome coverage. Following the trend, metabolomics once again lags its predecessors in the pursuit of single-cell systems biology. Single-cell genomics and transcriptomics saw rapid popularization in the last years [14], followed more recently by proteomics [15,16,17]. Here, in addition to the aforementioned technical hurdles, sensitivity also imposes a challenge for single-cell metabolomics. DNA and RNA analysis presents a significant technical advantage as the genetic material can be amplified, yielding considerably more sensitive detection over proteins and metabolites. Recent developments in proteomics have explored alternatives such as fluorescent tags providing a comprehensive increase in sensitivity. On the other hand, metabolites cannot be amplified, and their broad dynamic range of concentrations has a considerable impact on the observable metabolome. Furthermore, improving detection through derivatization reactions is also complicated by their wide chemical diversity and a higher propensity to structural modifications compared to bulkier proteins. Moreover, the minute concentrations and volumes of material represent an issue for using classical platforms relying on chromatographic separation. We can classify attempts to achieve cellular resolution metabolomics in three main groups: those that attempt at isolating enough material of a specific cell type to perform the analysis on platforms used for regular metabolomics, which we will refer to as single-cell-type metabolomics as coined by Reference [18]; those based on micromanipulation of single cells [19,20]; and those based on mass spectrometry imaging (MSI) [21,22]. In the next sections, we briefly summarize some of the main vantages and disadvantages of the different approaches (Figure 2). 3. Single-Cell and Single-Cell-Type Metabolomics In an ideal scenario, direct extraction of the inner content of a cell, or the cell as a whole, followed by metabolite profiling, represents the optimal procedure for preserving the natural cellular environment and assessing individual cellular heterogeneity across an organism. However, mass spectrometry performed within such low volumes and concentrations is generally limited to detecting only a small set of compounds. Moreover, performing chromatography with such material is even more challenging, and most platforms skip this method altogether. The lack of chromatographic separation results in increased matrix effects such as ion suppression, simply put, the signal reduction due to ionization interference between species simultaneously reaching the ionization source, therefore, negatively affecting the detection of most analytes. One of the few platforms established for such analysis, “Live-MS” performs single-cell metabolite profiling by sucking out the cell content under video-microscopy observations with the help of a metal-coated microcapillary such as a nanospray tip. The sample is further transferred into a mass spectrometer via a nano-electrospray ionization plume [19,23]. A more viable alternative from the analytical point of view is to sample many specific cells before the metabolomics experiment in single-cell-type experiments [18]. The main advantage is the possibility of using traditional LC/GC-MS-based platforms providing high throughput, optimal sensitivity, and coverage due to the chromatographic separation. The main limitations are imposed by the chosen cell sampling technique [24,25,26]. Some particularly exciting works include applications of laser microdissection (LMD)-based techniques such as laser microdissection and pressure catapulting (LMPC) and laser capture microdissection (LCM) [18], as well as fluorescence-activated cell sorting (FACS) [27]. LMD-based techniques are a great option as they preserve contextual information from spatial cell distribution. However, they are significantly limited in terms of throughput. LMD is a labor-intensive technique requiring an experienced operator to harvest the cells [18]. FACS, on the other hand, provides a high throughput alternative to isolating specific cells. However, the necessity to obtain single-cell suspensions is far from trivial, considerably affecting the metabolome [28]. Moreover, the inherent introduction of perturbations due to cell manipulation by all these techniques is particularly troublesome when considering the rapid changes of the metabolome with the turnover time of some metabolites being fractions of a second [29]. Despite these limitations, the recent improvement in data processing capacity and machine learning algorithms brings exciting advances to fill some of these gaps. A great example has recently been shown using image analysis algorithms, machine-learning, and high-throughput microscopy to recognize individual cells in suspensions or tissue and automatically guide extraction through LCM or micromanipulation in the so-called computer-assisted microscopy isolation (CAMI) [30]. Similarly, exciting improvements have also been developed for FACS [31]. However, the issues related to obtaining cell suspensions for this technique are likely hard to overcome and particularly challenging for plant sciences, as discussed below. 4. Mass Spectrometry Imaging (MSI) MSI is a general term encompassing multiple technologies capable of providing spatially resolved ionization of samples for mass spectrometry-based metabolite profiling [21,22]. The multiple techniques essentially provide different tradeoffs related to sample preparation, the lateral resolution of the ionization spot, degree of fragmentation, and ionization range (m/z). We briefly describe here some of the most common ionization platforms that we believe cover an attractive complementary space of features, namely matrix-assisted laser desorption/ionization (MALDI) [32,33], secondary ion mass spectrometry (SIMS) [34], desorption electrospray ionization (DESI) [35,36], and laser-ablation electrospray ionization (LAESI) [37] (Figure 3). MALDI is the most popular ionization method adaptable to MSI [38]. In MALDI, a matrix applied to the sample is excited by a laser; this energy is further transferred to the sample resulting in the ionization event [32,33]. It is particularly good at ionizing large molecules above 500 m/z, often suffering from matrix interference signals below this mass range [39]. Several groups have developed extensive work involving MALDI’s application as a platform for MSI with multiple applications into the analysis of plant samples [22]. Despite limited biologically relevant insights, these works tackle some of the main challenges in achieving comprehensive spatially resolved metabolomics, including sample preparation, the lateral resolution of ionization, and multiplex data acquisition. Preparation for MALDI usually comprises cryo-sectioning and lyophilizing a frozen sample embedded in some media before applying the matrix by either a sprayer or solvent-free sublimation [39]. These methods offer an advantage over cell isolation in terms of metabolome integrity whilst also preserving the relative localization of cells and allowing them to assess the intercellular space [40]. However, the process still lacks significant improvements in throughput. The choice of method for matrix deposition and its composition are particularly important factors in MALDI ionization. Comparing traditional spray and solvent-free sublimation methods as an example show that the former may promote metabolite delocalization, an issue amended by the latter method in the detriment of other metabolites not being ionized [41]. Moreover, matrix crystalline structure is a relevant factor limiting lateral resolution [42]. That said, matrix optimization is an active field in technological developments for MALDI imaging applications [43]. Several works have described matrix optimization for specific compound classes [44], as well as exciting approaches to expand the coverage based on derivatizations [45] and post-ionization strategies [46]. Another recent trend involves using nanoparticles instead of organic matrixes, and it shows promising results for ionizing the smaller range of metabolites and providing increased spatial resolution [47,48]. MALDI’s lateral resolution is usually in the range of 50–10 µm, even though some reports manage to achieve numbers as low as 2–5 µm in customized systems [40,49]. Factors limiting resolution again include the matrix structure and also qualitative aspects of the laser. UV lasers provide higher resolutions of up to 10 µm. However, they have several disadvantages compared to IR lasers, such as limitations in matrix absorption [21]. SIMS is an alternative to MALDI that relies on ion beams instead of a laser to ionize the samples [50]. Such a mechanism results in a more fragmented ionization and removes the necessity of any matrix and limitations due to laser’s diffraction limit, thus providing higher reproducibility and resolution below 2 µm [50]. Moreover, SIMS allows for the acquisition of 3D imaging through the use of dual beans. Indeed, all these advantages have been recently combined in a commercial system that includes the ultra-high resolution of orbitrap analyzers [51]. Despite the advantage of more straightforward sample preparation, SIMS-based platforms’ limitation is the need for samples to be ionized under a high vacuum. A few works try overcoming such limitations, for instance, via the use of cryogenic orbiSIMS to evaluate semi-volatile organic compounds that would otherwise be vaporized before ionization [52]. Nevertheless, DESI and LAESI offer promising alternatives for direct ionization of samples with minimal treatment. In DESI, a solvent stream originated from an electrospray probe is directed at an angle toward the sample at ambient pressure, propelling secondary ions to the analyzer [53]. One of the biggest limitations of DESI is its comparatively low resolution in the order of 100 µm [50]. Finally, LAESI combines laser ablation followed by post-ionization via an electrospray. A typical resolution is in the order of 200–300 µm; however, it can reach better resolution than DESI with the additional advantage of ionizing through multiple layers of tissue [21,54,55,56]. Despite its many advantages concerning in situ analysis, MSI platforms offer significant challenges regarding data analysis [57]. As a technique in its infancy, data processing standards, such as normalization, are still lacking. Indeed, only relatively recently, an open cross-platform data format was developed [58]. Quantification is also challenging, with few works providing absolute metabolite levels [59]. 5. Spatially Resolved Metabolomics in Plants: Current Status, Challenges, and Future Prospects The tremendous metabolic diversity that evolved in plants at the level of tissues and organs makes single-cell metabolomics a suitable tool for investigations targeting these cell-specific chemical signatures. To resolve metabolic diversity at a tissue level, spatial single-cell mass spectrometry, performed alone or in combination with single-cell mass spectrometry, has primarily been applied to plant science (Table 1). In Catharanthus roseus, for example, the combination of MSI and single-cell MS provided evidence of a developmentally driven process that segregates branches of the terpenoid indole alkaloid (TIA) biosynthetic pathway into specific anatomical structures. As idioblasts and laticifers differentiate while leaves grow and expand, single-cell MS detected the appearance of new metabolic intermediates so that an initial draft of the TIA pathway could be written in its completeness once leaves fully developed [60,61]. From an evolutionary perspective, this is a remarkable discovery as it provides an additional example of the parallel evolution between biochemical processes and anatomical structures, which often occurs in plants. In a similarly elegant experiment, Livingston and colleagues used a combination of different techniques, which included measurements of trichome intrinsic fluorescence and microcapillary-assisted metabolite extraction followed by GC-MS, UHPLC-MS/MS, and RNA sequencing analyses, to lay out the developmental trajectories of Cannabis sativa trichomes from sessile to stalked and the parallel changes occurring in the composition of their metabolites [62]. Questions concerning color pattern formation in flowers have also been answered via MSI. For instance, a recent study revealed that the deep-blue color of the nectar guides of Viola cornuta petals is due to the colocalization of the anthocyanin violanin and numerous colorless flavonol 3-O-glycosides [63]. By surrounding violanin, flavonol molecules prevent self-stacking and the consequent shift in the spectrum of light absorbance [64]. Additionally, they protect the chromophore of violanin from hydration, hence inhibiting the formation of colorless chalcones. It is not accidental that the great majority of these studies focused on tissues that accumulate specialized metabolites in a large abundance and can be relatively easily accessed, such as glandular trichomes [62,65,66,67], laticifers [60,68], and floral petals [69,70,71]. Indeed, as metabolites physiologically accumulate in these organs and structures, their concentration is already optimized to detect an MS signal of sufficient quality for the molecular identification of compounds. In all other cases, to achieve a proper concentration, metabolites must be extracted from a hundred thousand identical cells. For reasons that we explain below, harvesting such a large number of cells from plant tissues is extremely labor-intensive, as it emerges when comparisons with similar systems utilized in animal studies are made. Animal cell lines established after cell disaggregation from tissues followed by subculturing usually maintain similar physiological and biochemical characteristics as their organ of origin [72]. As such, immortal animal cell lines, for example, HeLa cells, have been successfully utilized to investigate metabolic responses to drugs and growth regulators [57]. Indeed, as metabolic changes in animal cell cultures mirror changes in intact organs, single-cell metabolomics is a powerful system to predict metabolic trajectories induced by medical treatments [73]. Conversely, plant liquid and solid (callus) cultures are made of cells in an undifferentiated status maintained with a balanced ratio of auxins and cytokinins [74]. As plant cell cultures are phenotypically and biochemically very distant from their differentiated counterparts, methods other than culturing must be adopted to collect a large number of cells of a specific lineage. Fluorescence-activated cell sorting (FACS) has been successfully employed to collect a large amount of GFP-tagged lines from plant tissues [27], but FACS applied to samples destined for the analysis of metabolites is a very challenging procedure. As reagents for cell protoplasting are potential contaminants of the MS detector and metabolites are prone to fast degradation, well-established methods that are used to collect RNA from fluorescent-tagged and sorted cells need further adjustments when applied to single-cell metabolomics [27]. In addition, the spatial distribution of differentiated cells in plant tissues and the prospect of obtaining protoplasts from these cells are tremendous limitations to the pursuit of harvesting cells of a single type in an amount that is sufficient for metabolite analysis. In roots, the continuous development and radial organization of layered tissues make single-cell collection via FACS relatively easily attainable, as well as facilitating the interpretation of imaging at a reasonable lateral resolution [75,76,77]. Conversely, organs that at maturity show a high degree of anatomical complexity, for example, flowers, are not equally suitable for such analyses. Not to mention that the process of protoplasting, which removes the cell wall—an intrinsic component of all plant cells— washes away compounds that cells secrete and deposit in the extracellular space. These compounds often have relevant physiological functions, for example, phenylalanine derivatives which confer protection to fungal pathogens [78]. Therefore, the removal of the cell wall can make data interpretation difficult, as it may weaken the link between chemical phenotypes and physiological functions. These challenges are at least partially resolved when tissues undergoing single-cell analyses are already composed of a large number of identical cells. Thus, for seeds and grains where cells with well-defined chemistry spatially cluster to form seed coat, embryo, and endosperm, MSI has largely been utilized to resolve in situ localization of metabolites. For example, in the oil-seed crops Camelina sativa and Brassica napus, as well as in Arabidopsis, MSI helped to determine the distribution of lipids in the embryo of wild-type and transgenic lines [79,80,81,82,83]. In barley and wheat, the spatial distribution of sugars and proteins between the endosperm and aleurone layer has been the main object of study [84,85,86]. Furthermore, in inbreds of maize, amino acids, sugar alcohols, organic acids, phospholipids, and triacylglycerols were observed within the embryo and radicle [87]. Unicellular structures and unicellular organisms such as pollen grains, algae, and microalgae (diatoms) represent another exception. For instance, metabolomics of pollen grains, which are unicellular haploid male gametophytes, has been performed with the most disparate array of techniques. As a result of these investigations, the molecular structure and composition of sporopollenin have recently been unraveled [88,89]. Sporopollenin is an extraordinarily inert and resistant polymer, the acquisition of which by land plants represents a focal adaptation to life outside water. The spectacular inertness of sporopollenin toward the most disparate analytical techniques made the search for its molecular structure hard to obtain, and at the same time, very desirable given the multitude of promising applications in the fields of material engineering and nanotechnology. Progress has also been made in the analysis of lipids, proteins, and the mechanisms of accumulation of flavonoid glycosides on the surface of pollen grains [90,91,92,93], as well as the metabolic processes underlying pollen germination and pollen tube elongation [94]. However, given the complexity of these chemical signals, their physiological function is not yet fully understood [95]. Finally, we briefly mention here that while collecting abundant pollen from male microsporangia of gymnosperms is usually easy endeavors, autogamous angiosperms generally produce a tiny amount of pollen, for which the collection of whole anthers is a necessary step. In algal research, live single-cell metabolomics helped elucidate the metabolic rearrangements occurring in response to environmental perturbations such as low nutrient and variation in light regimes [96,97,98]. Phytoplankton, which is primarily composed of microalgae and minor amounts of protists and bacteria, contributes to global biogeochemical cycles of carbon, nitrogen, phosphorus, and silicate. Therefore, gaining an understanding of the physiological status of phytoplankton cells holds excellent promises for environmental research. Besides, chemotyping of microalgae via pipelines that utilize live single-cell MS is currently exploited for taxonomic identification [99]. Despite the broad applications in environmental research, initial studies on algal metabolomics mostly focused on the model organism Chlamydomonas reinhardtii because of applications in the biotechnology industry and biofuel production [100]. In the recent past, metabolomic and transcriptomic approaches have been extensively utilized to investigate the responses to external determinants of algal growth such as temperature, light intensity, salinity, and nutrient availability [101,102], while today’s research mostly shifted toward functional genomic studies that aim at understanding the genetic mechanisms of this metabolic plasticity [103]. In plants, the application of single-cell MS to functional genomic studies has so far been scant, although initial studies that employed known Arabidopsis thaliana mutants as a proof of concept showed great promises [104,105,106]. More recently, the combination of direct infusion metabolomics and MSI was used to characterize the signaling pathway of feronia mutants, revealing an interesting phenotype associated with high levels of oxylipin arabidopsides, and suggesting chloroplastic localization [107]. Other new prospects of single-cell metabolomic applications pertain to research on plant–pathogen interactions where MSI is currently being utilized to analyze plant metabolites synthesized in response to pathogens’ infection. Here, initial studies that used to visualize metabolites present on the plant surface have further expanded to include analysis of metabolites that accumulate deeper in plant tissues, which can be seen after tissue fracturing and sectioning [108]. A recent study has shown this new approach’s relevance when it analyzed transcriptome and metabolome responses of susceptible and resistant soybean cultivars to aphid infestation. As aphids are insects that feed on phloem sap, plant metabolites conferring resistance to aphids are expected to be found in the phloem. Conversely, MSI revealed the accumulation of isoflavones in mesophyll and epidermal cells, suggesting a role for these compounds in the non-phloem defense response induced by feeding [109]. MSI has also been used to investigate the distribution of glucosinolates across Arabidopsis leaves and the response of lepidopterans oviposition to the detected metabolites’ concentrations [59]. In plant–bacteria symbiotic associations, MSI has been employed to study metabolite distribution in roots and nodules of wild-type and mutant genotypes of Medicago truncatula [110,111], and more recently of soybean [112,113]. Still, single-cell metabolomics has found very little application in the field of plant developmental research where the combination of functional genomics and metabolomics holds the promise to pave the way toward a better understanding of how and to what extent anatomy and metabolism are mutually coordinated. This is at least partially due to challenges related to sample preparation, as the presence of abundant water, cell wall, and cuticles make the process laborious in plants [114]. Additionally, the resolution and annotation of metabolites of the central pathway, which accumulate in lower abundance than specialized metabolites (see above), represent an additional challenge. Similarly, phytohormones that play crucial roles at the cellular level in the development and environmental responses are of great interest but challenging to detect even with traditional methods [115]. A work using the “Live-MS” platform [116] to investigate the response of two phytohormones, ABA and JA-Ile, have shown promising results being able to detect some of the expected changes. However, their results also point to limitations of the technique which still suffers from high variability. Finally, it is worth mentioning that plants primarily utilize hexoses and various sugar polymers for storage, transport, and organ to organ communication, of which MS annotation is not always easily attainable. ijms-21-08987-t001_Table 1Table 1 Summary of spatially resolved plant metabolomics works applying mass spectrometry-based platforms described in this review. Species Technique Cell-Type/Tissue Compounds Reference Arabidopsis FACS Roots Multiple [27] Arabidopsis MALDI Leaves Glucosinolates [60] Catharanthus roseus MALDI and Live-MS Laticifers and idioblasts from leaves TIA [61] Catharanthus roseus MALDI and Live-MS Laticifers, idioblast, parenchyma, and epidermal cells from stems TIA [62] Viola cornuta MALDI Petals Flavonoids [64] Rauvolfia tetraphylla DESI Stem, leaves, root, and fruits Indole alkaloids [69] Hypericum perforatum DESI Petals and leaves Hyperforin [71] Datura stramonium DESI Petals and leaves Sugars, atropine, and scopolamine [71] Maize MALDI Roots Amino acids [76] Maize MALDI Roots Lipids, sugars, and benzoxazinoid [77] Glycyrrhiza glabra MALDI Roots Flavonoids and triterpenoids [78] Camelina sativa MALDI Seed Lipids [80] Camelina sativa MALDI Seed Lipids [81] Camelina sativa MALDI Seed Lipids [82] Brassica napus MALDI Seed Lipids [83] Arabidopsis MALDI Seed Lipids [84] Barley MALDI Germinating seeds Multiple [85] Maize MALDI Germinating seeds Multiple [88] Lycopodium clavatum SIMS and MALDI Polen Sporopollenin [90] Poa alpina MALDI Polen Multiple [91] Arabidopsis MALDI Leaves Oxylipins [108] Rice MALDI Leaves Multiple [109] Soybean MALDI Leaves Multiple [109] Soybean MALDI Leaves Isoflavones [110] Medicago truncatula MALDI Root nodules Multiple [111] Medicago truncatula MALDI Root nodules Multiple [112] Soybean MALDI Root nodules Multiple [113] Soybean LAESI Root nodules Multiple [114] Vicia faba Live-MS Leaves Phytohormones [117] 6. Conclusions and Future Perspectives There are just over 40 different types of cells described in plant tissues [117]. As most metabolomics experiments capture data of whole tissues, our knowledge is largely biased toward prevailing cells such as mesophyll cells in leaves [118] and endosperm in seeds [119,120]. However, several works highlight the striking differences in cell-specific metabolism and the impact that less recurrent cell types have in regulating and integrating crucial physiological processes, including transpiration and photosynthesis [121,122]. Moreover, assessing metabolic heterogeneity across cells belonging to a tissue has the potential to unravel unforeseen details masked by averaging such populations of cells, thereby contributing to a deeper understanding of metabolic regulation [6]. Techniques for measuring single-cell metabolites have recently gone through considerable improvements providing exciting insights into metabolic compartmentalization. Nevertheless, some of the metabolomics most outstanding achievements rely on high throughput and comprehensive metabolome coverage. Both parameters are still considerably limited in the current single-cell and spatially resolved platforms. The many advantages of single-cell profiling described here represent an enormous potential when applied to large throughput experiments. Single-cell transcriptomics of different human tissues has recently been utilized to identify Quantitative Trait Loci (QTL) associated with expression and splicing variants (eQTLs and sQTLs, respectively) affected by the background genetic variation of different individuals [123]. A similar approach to plant tissues has not yet been adopted. However, it represents a significant potential if applied to large populations to understand, among others, the effect of environmental perturbations at a single-cell level. Improvements in various aspects of mass spectrometry aspects, particularly resolution and sensitivity, have been instrumental in facilitating the measurement of the spatial distribution of metabolites through single-cell and MSI platforms. The introduction and broad adoption of other technologies into metabolomics platforms, such as nanoLC and ion mobility, are likely to play important roles in further reducing issues concerning limited sample and sensitivity in single-cell metabolomics and matrix effects in MSI, respectively. Better integration of current technologies with other imaging platforms such as microscopy also offers a promising way to improve experiments throughput and information [73]. Finally, as these technologies mature, we can foresee their adoption to even the most challenging applications in current tissue level metabolomics. A recurrent question that has proven essential to characterize metabolism is the definition of metabolic fluxes rather than a simple description of relative metabolite levels as routinely performed [124]. Methods for integrating multi-omics of single cells are also an exciting boundary to be crossed [125]. We can anticipate considerable hurdles for generating such datasets. However, this could represent an outstanding means of reducing experimental complexity while improving the statistical power of systems biology studies. Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. Author Contributions L.P.d.S. and M.B. surveyed the literature, all authors contributed with writing and revising the manuscript. All authors have read and agreed to the published version of the manuscript. Funding This research was funded by the European Commission, Horizon 2020 Framework Programme, H2020 European Institute of Innovation and Technology, grant numbers: FPA No. 664620, No 739582, SGA-CSA No 664621. Conflicts of Interest The authors declare no conflict of interest. Abbreviations DESI Desorption electrospray ionization FACS Fluorescence-activated cell sorting GC Gas chromatography LAESI Laser-ablation electrospray ionization LC Liquid chromatography LCM Laser capture microdissection LMD Laser microdissection LMPC Laser microdissection and pressure catapulting MALDI Matrix-assisted laser desorption/ionization MS Mass spectrometry MSI Mass spectrometry imaging QTOF Quadrupole time-of-flight mass spectrometer SIMS Secondary ion mass spectrometry TIA Terpenoid indole alkaloid UHPLC Ultra-high-performance liquid chromatography UV Ultraviolet Figure 1 The effect of cell heterogeneity, exemplified by different cell types in a leaf, in metabolomics analysis. (A) Diagram of a leaf cross-section depicting typical cell types. (B) Hypothetical values for a metabolite X differentially accumulated in multiple cell types highlighting the averaging effect of pooling cells together in a traditional metabolomics experiment. Trichomes (TC), epidermal cells (EP), mesophyll cells (MS), bundle sheath cells (BS), and guard cells (GC). Figure 2 Overview of experimental steps and data structure from the different approaches for cell-specific metabolomics. Figure 3 Schematic representation of the different ionization strategies used for mass spectrometry imaging (MSI). (A) MALDI, (B) secondary ion mass spectrometry (SIMS), (C) desorption electrospray ionization (DESI), (D) laser-ablation electrospray ionization (LAESI). ==== Refs References 1. Daloso D.M. Müller K. Obata T. Florian A. Tohge T. Bottcher A. Riondet C. Bariat L. Carrari F. Nunes-Nesi A. Thioredoxin, a master regulator of the tricarboxylic acid cycle in plant mitochondria Proc. Natl. Acad. Sci. USA 2015 112 E1392 E1400 10.1073/pnas.1424840112 25646482 2. Araújo W.L. Ishizaki K. Nunes-Nesi A. Larson T.R. Tohge T. Krahnert I. Witt S. Obata T. Schauer N. Graham I.A. Identification of the 2-Hydroxyglutarate and Isovaleryl-CoA Dehydrogenases as Alternative Electron Donors Linking Lysine Catabolism to the Electron Transport Chain of Arabidopsis Mitochondria Plant Cell 2010 22 1549 1563 10.1105/tpc.110.075630 20501910 3. Tohge T. Wendenburg R. Ishihara H. Nakabayashi R. Watanabe M. Sulpice R. Hoefgen R. Takayama H. Saito K. Stitt M. Characterization of a recently evolved flavonol-phenylacyltransferase gene provides signatures of natural light selection in Brassicaceae Nat. Commun. 2016 7 12399 10.1038/ncomms12399 27545969 4. Perez de Souza L. Garbowicz K. Brotman Y. Tohge T. Fernie A.R. The Acetate Pathway Supports Flavonoid and Lipid Biosynthesis in Arabidopsis Plant Physiol. 2020 182 857 869 10.1104/pp.19.00683 31719153 5. Cárdenas P.D. Sonawane P.D. Heinig U. Jozwiak A. Panda S. Abebie B. Kazachkova Y. Pliner M. Unger T. Wolf D. Pathways to defense metabolites and evading fruit bitterness in genus Solanum evolved through 2-oxoglutarate-dependent dioxygenases Nat. Commun. 2019 10 5169 10.1038/s41467-019-13211-4 31727889 6. Altschuler S.J. Wu L.F. Cellular Heterogeneity: Do Differences Make a Difference? Cell 2010 141 559 563 10.1016/j.cell.2010.04.033 20478246 7. Fleming A. Metabolic aspects of organogenesis in the shoot apical meristem J. Exp. Bot. 2006 57 1863 1870 10.1093/jxb/erj178 16690623 8. Medeiros D.B. Perez Souza L. Antunes W.C. Araújo W.L. Daloso D.M. Fernie A.R. Sucrose breakdown within guard cells provides substrates for glycolysis and glutamine biosynthesis during light-induced stomatal opening Plant J. 2018 94 583 594 10.1111/tpj.13889 29543357 9. Wang H. Yan S. Xin H. Huang W. Zhang H. Teng S. Yu Y.-C. Fernie A.R. Lu X. Li P. A Subsidiary Cell-Localized Glucose Transporter Promotes Stomatal Conductance and Photosynthesis Plant Cell 2019 31 1328 1343 10.1105/tpc.18.00736 30996077 10. Arrivault S. Obata T. Szecówka M. Mengin V. Guenther M. Hoehne M. Fernie A.R. Stitt M. Metabolite pools and carbon flow during C4 photosynthesis in maize: 13CO2 labeling kinetics and cell type fractionation J. Exp. Bot. 2016 68 283 298 10.1093/jxb/erw414 27834209 11. Islam M.M. Al-Siyabi A. Saha R. Obata T. Dissecting metabolic flux in C4 plants: Experimental and theoretical approaches Phytochem. Rev. 2018 17 1253 1274 10.1007/s11101-018-9579-8 12. Weissmann S. Ma F. Furuyama K. Gierse J. Berg H. Shao Y. Taniguchi M. Allen D.K. Brutnell T.P. Interactions of C4 Subtype Metabolic Activities and Transport in Maize Are Revealed through the Characterization of DCT2 Mutants Plant Cell 2016 28 466 484 10.1105/tpc.15.00497 26813621 13. Schenck C.A. Last R.L. Location, location! cellular relocalization primes specialized metabolic diversification FEBS J. 2020 287 1359 1368 10.1111/febs.15097 31623016 14. Luo C. Fernie A.R. Yan J. Single-Cell Genomics and Epigenomics: Technologies and Applications in Plants Trends Plant Sci. 2020 25 1030 1040 10.1016/j.tplants.2020.04.016 32532595 15. Labib M. Kelley S.O. Single-cell analysis targeting the proteome Nat. Rev. Chem. 2020 4 143 158 10.1038/s41570-020-0162-7 16. Marx V. A dream of single-cell proteomics Nat. Methods 2019 16 809 812 10.1038/s41592-019-0540-6 31406385 17. Slavov N. Unpicking the proteome in single cells Science 2020 367 512 513 10.1126/science.aaz6695 32001644 18. Misra B.B. Assmann S.M. Chen S. Plant single-cell and single-cell-type metabolomics Trends Plant Sci. 2014 19 637 646 10.1016/j.tplants.2014.05.005 24946988 19. Fujii T. Matsuda S. Tejedor M.L. Esaki T. Sakane I. Mizuno H. Tsuyama N. Masujima T. Direct metabolomics for plant cells by live single-cell mass spectrometry Nat. Protoc. 2015 10 1445 1456 10.1038/nprot.2015.084 26313480 20. Yang Y. Huang Y. Wu J. Liu N. Deng J. Luan T. Single-cell analysis by ambient mass spectrometry TrAC Trends Anal. Chem. 2017 90 14 26 10.1016/j.trac.2017.02.009 21. Bjarnholt N. Li B. D’Alvise J. Janfelt C. Mass spectrometry imaging of plant metabolites–principles and possibilities Nat. Prod. Rep. 2014 31 818 837 10.1039/C3NP70100J 24452137 22. Hansen R.L. Lee Y.J. High-Spatial Resolution Mass Spectrometry Imaging: Toward Single Cell Metabolomics in Plant Tissues Chem. Rec. 2018 18 65 77 10.1002/tcr.201700027 28685965 23. Masujima T. Live Single-cell Mass Spectrometry Anal. Sci. 2009 25 953 960 10.2116/analsci.25.953 19667470 24. Gross A. Schoendube J. Zimmermann S. Steeb M. Zengerle R. Koltay P. Technologies for Single-Cell Isolation Int. J. Mol. Sci. 2015 16 16897 16919 10.3390/ijms160816897 26213926 25. Zhang X. Wei X. Wei Y. Chen M. Wang J. The up-to-date strategies for the isolation and manipulation of single cells Talanta 2020 218 121147 10.1016/j.talanta.2020.121147 32797903 26. Couvillion S.P. Zhu Y. Nagy G. Adkins J.N. Ansong C. Renslow R.S. Piehowski P.D. Ibrahim Y.M. Kelly R.T. Metz T.O. New mass spectrometry technologies contributing towards comprehensive and high throughput omics analyses of single cells Analyst 2019 144 794 807 10.1039/C8AN01574K 30507980 27. Moussaieff A. Rogachev I. Brodsky L. Malitsky S. Toal T.W. Belcher H. Yativ M. Brady S.M. Benfey P.N. Aharoni A. High-resolution metabolic mapping of cell types in plant roots Proc. Natl. Acad. Sci. USA 2013 110 E1232 E1241 10.1073/pnas.1302019110 23476065 28. Reichard A. Asosingh K. Best Practices for Preparing a Single Cell Suspension from Solid Tissues for Flow Cytometry Cytom. Part A 2019 95 219 226 10.1002/cyto.a.23690 29. Arrivault S. Guenther M. Ivakov A. Feil R. Vosloh D. Van Dongen J.T. Sulpice R. Stitt M. Use of reverse-phase liquid chromatography, linked to tandem mass spectrometry, to profile the Calvin cycle and other metabolic intermediates in Arabidopsis rosettes at different carbon dioxide concentrations Plant J. 2009 59 826 839 10.1111/j.1365-313X.2009.03902.x 19453453 30. Brasko C. Smith K. Molnar C. Farago N. Hegedus L. Balind A. Balassa T. Szkalisity A. Sukosd F. Kocsis K. Intelligent image-based in situ single-cell isolation Nat. Commun. 2018 9 226 10.1038/s41467-017-02628-4 29335532 31. Isozaki A. Mikami H. Hiramatsu K. Sakuma S. Kasai Y. Iino T. Yamano T. Yasumoto A. Oguchi Y. Suzuki N. A practical guide to intelligent image-activated cell sorting Nat. Protoc. 2019 14 2370 2415 10.1038/s41596-019-0183-1 31278398 32. Karas M. Bahr U. Matrix-Assisted Laser Desorption-Ionization (MALDI) Mass Spectrometry: Principles and Applications Selected Topics in Mass Spectrometry in the Biomolecular Sciences Caprioli R.M. Malorni A. Sindona G. Springer Netherlands Dordrecht, The Netherlands 1997 33 53 33. Schwamborn K. Kriegsmann M. Weichert W. MALDI imaging mass spectrometry—From bench to bedside Biochim. Biophys. Acta (BBA) Proteins Proteom. 2017 1865 776 783 10.1016/j.bbapap.2016.10.014 34. Walker A.V. Secondary Ion Mass Spectrometry Encyclopedia of Spectroscopy and Spectrometry 3rd ed. Lindon J.C. Tranter G.E. Koppenaal D.W. Academic Press Oxford, UK 2017 44 49 35. Hollerbach A. Ayrton S. Jarmusch A. Graham Cooks R. Desorption Electrospray Ionization: Methodology and Applications Encyclopedia of Spectroscopy and Spectrometry 3rd ed. Lindon J.C. Tranter G.E. Koppenaal D.W. Academic Press Oxford, UK 2017 401 408 36. Cooks R.G. Ouyang Z. Takats Z. Wiseman J.M. Ambient Mass Spectrometry Science 2006 311 1566 1570 10.1126/science.1119426 16543450 37. Wu C. Dill A.L. Eberlin L.S. Cooks R.G. Ifa D.R. Mass spectrometry imaging under ambient conditions Mass Spectrom. Rev. 2013 32 218 243 10.1002/mas.21360 22996621 38. Boughton B.A. Hamilton B. Spatial Metabolite Profiling by Matrix-Assisted Laser Desorption Ionization Mass Spectrometry Imaging Metabolomics: From Fundamentals to Clinical Applications Sussulini A. Springer International Publishing Cham, Switzerland 2017 291 321 39. Sturtevant D. Lee Y.-J. Chapman K.D. Matrix assisted laser desorption/ionization-mass spectrometry imaging (MALDI-MSI) for direct visualization of plant metabolites in situ Curr. Opin. Biotechnol. 2016 37 53 60 10.1016/j.copbio.2015.10.004 26613199 40. Korte A.R. Yandeau-Nelson M.D. Nikolau B.J. Lee Y.J. Subcellular-level resolution MALDI-MS imaging of maize leaf metabolites by MALDI-linear ion trap-Orbitrap mass spectrometer Anal. Bioanal. Chem. 2015 407 2301 2309 10.1007/s00216-015-8460-5 25618761 41. Gemperline E. Rawson S. Li L. Optimization and Comparison of Multiple MALDI Matrix Application Methods for Small Molecule Mass Spectrometric Imaging Anal. Chem. 2014 86 10030 10035 10.1021/ac5028534 25331774 42. Li S. Zhang Y. Liu J.a. Han J. Guan M. Yang H. Lin Y. Xiong S. Zhao Z. Electrospray deposition device used to precisely control the matrix crystal to improve the performance of MALDI MSI Sci. Rep. 2016 6 37903 10.1038/srep37903 27885266 43. Calvano C.D. Monopoli A. Cataldi T.R.I. Palmisano F. MALDI matrices for low molecular weight compounds: An endless story? Anal. Bioanal. Chem. 2018 410 4015 4038 10.1007/s00216-018-1014-x 29682685 44. Korte A.R. Lee Y.J. MALDI-MS analysis and imaging of small molecule metabolites with 1,5-diaminonaphthalene (DAN) J. Mass Spectrom. 2014 49 737 741 10.1002/jms.3400 25044901 45. Dueñas M.E. Larson E.A. Lee Y.J. Toward Mass Spectrometry Imaging in the Metabolomics Scale: Increasing Metabolic Coverage Through Multiple On-Tissue Chemical Modifications Front. Plant Sci. 2019 10 10.3389/fpls.2019.00860 46. Soltwisch J. Kettling H. Vens-Cappell S. Wiegelmann M. Müthing J. Dreisewerd K. Mass spectrometry imaging with laser-induced postionization Science 2015 348 211 215 10.1126/science.aaa1051 25745064 47. Yagnik G.B. Hansen R.L. Korte A.R. Reichert M.D. Vela J. Lee Y.J. Large Scale Nanoparticle Screening for Small Molecule Analysis in Laser Desorption Ionization Mass Spectrometry Anal. Chem. 2016 88 8926 8930 10.1021/acs.analchem.6b02732 27573492 48. Hansen R.L. Dueñas M.E. Lee Y.J. Sputter-Coated Metal Screening for Small Molecule Analysis and High-Spatial Resolution Imaging in Laser Desorption Ionization Mass Spectrometry J. Am. Soc. Mass Spectrom. 2019 30 299 308 10.1007/s13361-018-2081-0 30341582 49. Römpp A. Schäfer K.C. Guenther S. Wang Z. Köstler M. Leisner A. Paschke C. Schramm T. Spengler B. High-resolution atmospheric pressure infrared laser desorption/ionization mass spectrometry imaging of biological tissue Anal. Bioanal. Chem. 2013 405 6959 6968 10.1007/s00216-013-7180-y 23877173 50. Touboul D. Brunelle A. What more can TOF-SIMS bring than other MS imaging methods? Bioanalysis 2016 8 367 369 10.4155/bio.16.11 26891858 51. Passarelli M.K. Pirkl A. Moellers R. Grinfeld D. Kollmer F. Havelund R. Newman C.F. Marshall P.S. Arlinghaus H. Alexander M.R. The 3D OrbiSIMS—Label-free metabolic imaging with subcellular lateral resolution and high mass-resolving power Nat. Methods 2017 14 1175 1183 10.1038/nmeth.4504 29131162 52. Newell C.L. Vorng J.-L. MacRae J.I. Gilmore I.S. Gould A.P. Cryogenic OrbiSIMS Localizes Semi-Volatile Molecules in Biological Tissues Angew. Chem. Int. Ed. 2020 59 18194 18200 10.1002/anie.202006881 53. Wiseman J.M. Ifa D.R. Song Q. Cooks R.G. Tissue Imaging at Atmospheric Pressure Using Desorption Electrospray Ionization (DESI) Mass Spectrometry Angew. Chem. Int. Ed. 2006 45 7188 7192 10.1002/anie.200602449 54. Kulkarni P. Wilschut R.A. Verhoeven K.J.F. van der Putten W.H. Garbeva P. LAESI mass spectrometry imaging as a tool to differentiate the root metabolome of native and range-expanding plant species Planta 2018 248 1515 1523 10.1007/s00425-018-2989-4 30140978 55. Samarah L.Z. Khattar R. Tran T.H. Stopka S.A. Brantner C.A. Parlanti P. Veličković D. Shaw J.B. Agtuca B.J. Stacey G. Single-Cell Metabolic Profiling: Metabolite Formulas from Isotopic Fine Structures in Heterogeneous Plant Cell Populations Anal. Chem. 2020 92 7289 7298 10.1021/acs.analchem.0c00936 32314907 56. Etalo D.W. De Vos R.C.H. Joosten M.H.A.J. Hall R.D. Spatially Resolved Plant Metabolomics: Some Potentials and Limitations of Laser-Ablation Electrospray Ionization Mass Spectrometry Metabolite Imaging Plant Physiol. 2015 169 1424 1435 10.1104/pp.15.01176 26392264 57. Alexandrov T. Spatial Metabolomics and Imaging Mass Spectrometry in the Age of Artificial Intelligence Annu. Rev. Biomed. Data Sci. 2020 3 61 87 10.1146/annurev-biodatasci-011420-031537 58. Schramm T. Hester Z. Klinkert I. Both J.-P. Heeren R.M.A. Brunelle A. Laprévote O. Desbenoit N. Robbe M.-F. Stoeckli M. imzML—A common data format for the flexible exchange and processing of mass spectrometry imaging data J. Proteom. 2012 75 5106 5110 10.1016/j.jprot.2012.07.026 22842151 59. Shroff R. Schramm K. Jeschke V. Nemes P. Vertes A. Gershenzon J. Svatoš A. Quantification of plant surface metabolites by matrix-assisted laser desorption–ionization mass spectrometry imaging: Glucosinolates on Arabidopsis thaliana leaves Plant J. 2015 81 961 972 10.1111/tpj.12760 25600688 60. Yamamoto K. Takahashi K. Caputi L. Mizuno H. Rodriguez-Lopez C.E. Iwasaki T. Ishizaki K. Fukaki H. Ohnishi M. Yamazaki M. The complexity of intercellular localisation of alkaloids revealed by single-cell metabolomics New Phytol. 2019 224 848 859 10.1111/nph.16138 31436868 61. Yamamoto K. Takahashi K. Mizuno H. Anegawa A. Ishizaki K. Fukaki H. Ohnishi M. Yamazaki M. Masujima T. Mimura T. Cell-specific localization of alkaloids in Catharanthus roseus stem tissue measured with Imaging MS and Single-cell MS Proc. Natl. Acad. Sci. USA 2016 113 3891 3896 10.1073/pnas.1521959113 27001858 62. Livingston S.J. Quilichini T.D. Booth J.K. Wong D.C.J. Rensing K.H. Laflamme-Yonkman J. Castellarin S.D. Bohlmann J. Page J.E. Samuels A.L. Cannabis glandular trichomes alter morphology and metabolite content during flower maturation Plant J. 2020 101 37 56 10.1111/tpj.14516 31469934 63. Sugahara K. Kitao K. Watanabe T. Yamagaki T. Imaging Mass Spectrometry Analysis of Flavonoids in Blue Viola Petals and Their Enclosure Effects on Violanin during Color Expression Anal. Chem. 2019 91 896 902 10.1021/acs.analchem.8b03815 30521315 64. Goto T. Kondo T. Structure and Molecular Stacking of Anthocyanins—Flower Color Variation Angew. Chem. Int. Ed. 1991 30 17 33 10.1002/anie.199100171 65. Bertea C.M. Freije J.R. van der Woude H. Verstappen F.W.A. Perk L. Marquez V. De Kraker J.W. Posthumus M.A. Jansen B.J.M. de Groot A. Identification of Intermediates and Enzymes Involved in the Early Steps of Artemisinin Biosynthesis in Artemisia annua Planta Med. 2005 71 40 47 10.1055/s-2005-837749 15678372 66. Frerigmann H. Böttcher C. Baatout D. Gigolashvili T. Glucosinolates are produced in trichomes of Arabidopsis thaliana Front. Plant Sci. 2012 3 10.3389/fpls.2012.00242 23115560 67. Schilmiller A. Shi F. Kim J. Charbonneau A.L. Holmes D. Daniel Jones A. Last R.L. Mass spectrometry screening reveals widespread diversity in trichome specialized metabolites of tomato chromosomal substitution lines Plant J. 2010 62 391 403 10.1111/j.1365-313X.2010.04154.x 20113441 68. Mohana Kumara P. Uma Shaanker R. Pradeep T. UPLC and ESI-MS analysis of metabolites of Rauvolfia tetraphylla L. and their spatial localization using desorption electrospray ionization (DESI) mass spectrometric imaging Phytochemistry 2019 159 20 29 10.1016/j.phytochem.2018.11.009 30562679 69. Shi F. Flanigan P.M. Archer J.J. Levis R.J. Ambient Molecular Analysis of Biological Tissue Using Low-Energy, Femtosecond Laser Vaporization and Nanospray Postionization Mass Spectrometry J. Am. Soc. Mass Spectrom. 2016 27 542 551 10.1007/s13361-015-1302-z 26667178 70. Thunig J. Hansen S.H. Janfelt C. Analysis of Secondary Plant Metabolites by Indirect Desorption Electrospray Ionization Imaging Mass Spectrometry Anal. Chem. 2011 83 3256 3259 10.1021/ac2004967 21473636 71. Hölscher D. Shroff R. Knop K. Gottschaldt M. Crecelius A. Schneider B. Heckel D.G. Schubert U.S. Svatoš A. Matrix-free UV-laser desorption/ionization (LDI) mass spectrometric imaging at the single-cell level: Distribution of secondary metabolites of Arabidopsis thaliana and Hypericum species Plant J. 2009 60 907 918 10.1111/j.1365-313X.2009.04012.x 19732382 72. Bhatia S. Naved T. Sardana S. Introduction to Pharmaceutical Biotechnology, Volume 3 Animal Tissue Culture and Biopharmaceuticals IOP Publishing Bristol, UK 2019 73. Rappez L. Stadler M. Triana S. Phapale P. Heikenwalder M. Alexandrov T. Spatial single-cell profiling of intracellular metabolomes in situ bioRxiv 2019 510222 10.1101/510222 74. Fehér A. Callus, Dedifferentiation, Totipotency, Somatic Embryogenesis: What These Terms Mean in the Era of Molecular Plant Biology? Front. Plant Sci. 2019 10 10.3389/fpls.2019.00536 75. O’Neill K.C. Lee Y.J. Visualizing Genotypic and Developmental Differences of Free Amino Acids in Maize Roots With Mass Spectrometry Imaging Front. Plant Sci. 2020 11 10.3389/fpls.2020.00639 76. Feenstra A.D. Dueñas M.E. Lee Y.J. Five Micron High Resolution MALDI Mass Spectrometry Imaging with Simple, Interchangeable, Multi-Resolution Optical System J. Am. Soc. Mass Spectrom. 2017 28 434 442 10.1007/s13361-016-1577-8 28050871 77. Li B. Bhandari D.R. Janfelt C. Römpp A. Spengler B. Natural products in Glycyrrhiza glabra (licorice) rhizome imaged at the cellular level by atmospheric pressure matrix-assisted laser desorption/ionization tandem mass spectrometry imaging Plant J. 2014 80 161 171 10.1111/tpj.12608 25040821 78. Doppler M. Kluger B. Bueschl C. Steiner B. Buerstmayr H. Lemmens M. Krska R. Adam G. Schuhmacher R. Stable Isotope-Assisted Plant Metabolomics: Investigation of Phenylalanine-Related Metabolic Response in Wheat Upon Treatment With the Fusarium Virulence Factor Deoxynivalenol Front. Plant Sci. 2019 10 10.3389/fpls.2019.01137 31736983 79. Horn P.J. Silva J.E. Anderson D. Fuchs J. Borisjuk L. Nazarenus T.J. Shulaev V. Cahoon E.B. Chapman K.D. Imaging heterogeneity of membrane and storage lipids in transgenic Camelina sativa seeds with altered fatty acid profiles Plant J. 2013 76 138 150 10.1111/tpj.12278 23808562 80. Marmon S. Sturtevant D. Herrfurth C. Chapman K. Stymne S. Feussner I. Two Acyltransferases Contribute Differently to Linolenic Acid Levels in Seed Oil Plant Physiol. 2017 173 2081 2095 10.1104/pp.16.01865 28235891 81. Usher S. Han L. Haslam R.P. Michaelson L.V. Sturtevant D. Aziz M. Chapman K.D. Sayanova O. Napier J.A. Tailoring seed oil composition in the real world: Optimising omega-3 long chain polyunsaturated fatty acid accumulation in transgenic Camelina sativa Sci. Rep. 2017 7 6570 10.1038/s41598-017-06838-0 28747792 82. Woodfield H.K. Sturtevant D. Borisjuk L. Munz E. Guschina I.A. Chapman K. Harwood J.L. Spatial and Temporal Mapping of Key Lipid Species in Brassica napus Seeds Plant Physiol. 2017 173 1998 2009 10.1104/pp.16.01705 28188274 83. Sturtevant D. Dueñas M.E. Lee Y.-J. Chapman K.D. Three-dimensional visualization of membrane phospholipid distributions in Arabidopsis thaliana seeds: A spatial perspective of molecular heterogeneity Biochim. Biophys. Acta (BBA) Mol. Cell Biol. Lipids 2017 1862 268 281 10.1016/j.bbalip.2016.11.012 84. Gorzolka K. Kölling J. Nattkemper T.W. Niehaus K. Spatio-Temporal Metabolite Profiling of the Barley Germination Process by MALDI MS Imaging PLoS ONE 2016 11 e0150208 10.1371/journal.pone.0150208 26938880 85. Rolletschek H. Melkus G. Grafahrend-Belau E. Fuchs J. Heinzel N. Schreiber F. Jakob P.M. Borisjuk L. Combined Noninvasive Imaging and Modeling Approaches Reveal Metabolic Compartmentation in the Barley Endosperm Plant Cell 2011 23 3041 3054 10.1105/tpc.111.087015 21856793 86. Shewry P.R. Wan Y. Hawkesford M.J. Tosi P. Spatial distribution of functional components in the starchy endosperm of wheat grains J. Cereal Sci. 2020 91 102869 10.1016/j.jcs.2019.102869 32089586 87. Feenstra A.D. Alexander L.E. Song Z. Korte A.R. Yandeau-Nelson M.D. Nikolau B.J. Lee Y.J. Spatial Mapping and Profiling of Metabolite Distributions during Germination Plant Physiol. 2017 174 2532 2548 10.1104/pp.17.00652 28634228 88. Li F.-S. Phyo P. Jacobowitz J. Hong M. Weng J.-K. The molecular structure of plant sporopollenin Nat. Plants 2019 5 41 46 10.1038/s41477-018-0330-7 30559416 89. Mikhael A. Jurcic K. Schneider C. Karr D. Fisher G.L. Fridgen T.D. Diego-Taboada A. Georghiou P.E. Mackenzie G. Banoub J. Demystifying and unravelling the molecular structure of the biopolymer sporopollenin Front. Plant Sci. 2020 34 e8740 10.1002/rcm.8740 32003875 90. Diehn S. Zimmermann B. Tafintseva V. Seifert S. Bağcıoğlu M. Ohlson M. Weidner S. Fjellheim S. Kohler A. Kneipp J. Combining Chemical Information From Grass Pollen in Multimodal Characterization Plant Cell 2020 10 10.3389/fpls.2019.01788 32082348 91. Grunewald S. Marillonnet S. Hause G. Haferkamp I. Neuhaus H.E. Veß A. Hollemann T. Vogt T. The Tapetal Major Facilitator NPF2.8 Is Required for Accumulation of Flavonol Glycosides on the Pollen Surface in Arabidopsis thaliana Plant Cell 2020 32 1727 1748 10.1105/tpc.19.00801 32156687 92. Kenđel A. Zimmermann B. Chemical Analysis of Pollen by FT-Raman and FTIR Spectroscopies Front. Plant Sci. 2020 11 10.3389/fpls.2020.00352 93. Wan X. Wu S. Li Z. An X. Tian Y. Lipid Metabolism: Critical Roles in Male Fertility and Other Aspects of Reproductive Development in Plants Mol. Plant 2020 13 955 983 10.1016/j.molp.2020.05.009 32434071 94. Selinski J. Scheibe R. Pollen tube growth: Where does the energy come from? Plant Signal. Behav. 2014 9 e977200 10.4161/15592324.2014.977200 25482752 95. Borghi M. Fernie A.R. Outstanding questions in flower metabolism Plant J. 2020 103 1275 1288 10.1111/tpj.14814 32410253 96. Baumeister T.U.H. Vallet M. Kaftan F. Svatoš A. Pohnert G. Live Single-Cell Metabolomics With Matrix-Free Laser/Desorption Ionization Mass Spectrometry to Address Microalgal Physiology Front. Plant Sci. 2019 10 10.3389/fpls.2019.00172 30833957 97. Jaschinski T. Helfrich E.J.N. Bock C. Wolfram S. Svatoš A. Hertweck C. Pohnert G. Matrix-free single-cell LDI-MS investigations of the diatoms Coscinodiscus granii and Thalassiosira pseudonana J. Mass Spectrom. 2014 49 136 144 10.1002/jms.3316 24677306 98. Sun M. Yang Z. Wawrik B. Metabolomic Fingerprints of Individual Algal Cells Using the Single-Probe Mass Spectrometry Technique Front. Plant Sci. 2018 9 10.3389/fpls.2018.00571 29760716 99. Baumeister T.U.H. Vallet M. Kaftan F. Guillou L. Svatoš A. Pohnert G. Identification to species level of live single microalgal cells from plankton samples with matrix-free laser/desorption ionization mass spectrometry Metabolomics 2020 16 28 10.1007/s11306-020-1646-7 32090296 100. Bölling C. Fiehn O. Metabolite Profiling of Chlamydomonas reinhardtii under Nutrient Deprivation Plant Physiol. 2005 139 1995 2005 10.1104/pp.105.071589 16306140 101. Lee D.Y. Park J.-J. Barupal D.K. Fiehn O. System response of metabolic networks in Chlamydomonas reinhardtii to total available ammonium Mol. Cell Proteom. 2012 11 973 988 10.1074/mcp.M111.016733 102. Schreiber F. Ackermann M. Environmental drivers of metabolic heterogeneity in clonal microbial populations Curr. Opin. Biotechnol. 2020 62 202 211 10.1016/j.copbio.2019.11.018 31874388 103. Krismer J. Tamminen M. Fontana S. Zenobi R. Narwani A. Single-cell mass spectrometry reveals the importance of genetic diversity and plasticity for phenotypic variation in nitrogen-limited Chlamydomonas ISME J. 2017 11 988 998 10.1038/ismej.2016.167 27935588 104. Cha S. Song Z. Nikolau B.J. Yeung E.S. Direct Profiling and Imaging of Epicuticular Waxes on Arabidopsis thaliana by Laser Desorption/Ionization Mass Spectrometry Using Silver Colloid as a Matrix Anal. Chem. 2009 81 2991 3000 10.1021/ac802615r 19290666 105. Günl M. Neumetzler L. Kraemer F. de Souza A. Schultink A. Pena M. York W.S. Pauly M. AXY8 Encodes an α-Fucosidase, Underscoring the Importance of Apoplastic Metabolism on the Fine Structure of Arabidopsis Cell Wall Polysaccharides Plant Cell 2011 23 4025 4040 10.1105/tpc.111.089193 22080600 106. Korte A.R. Song Z. Nikolau B.J. Lee Y.J. Mass spectrometric imaging as a high-spatial resolution tool for functional genomics: Tissue-specific gene expression of TT7 inferred from heterogeneous distribution of metabolites in Arabidopsis flowers Anal. Methods 2012 4 474 481 10.1039/C2AY05618F 107. Hansen R.L. Guo H. Yin Y. Lee Y.J. FERONIA mutation induces high levels of chloroplast-localized Arabidopsides which are involved in root growth Plant J. 2019 97 341 351 10.1111/tpj.14123 30300943 108. Klein A.T. Yagnik G.B. Hohenstein J.D. Ji Z. Zi J. Reichert M.D. MacIntosh G.C. Yang B. Peters R.J. Vela J. Investigation of the Chemical Interface in the Soybean–Aphid and Rice–Bacteria Interactions Using MALDI-Mass Spectrometry Imaging Anal. Chem. 2015 87 5294 5301 10.1021/acs.analchem.5b00459 25914940 109. Hohenstein J.D. Studham M.E. Klein A. Kovinich N. Barry K. Lee Y.-J. MacIntosh G.C. Transcriptional and Chemical Changes in Soybean Leaves in Response to Long-Term Aphid Colonization Front. Plant Sci. 2019 10 10.3389/fpls.2019.00310 30930925 110. Gemperline E. Jayaraman D. Maeda J. Ané J.-M. Li L. Multifaceted investigation of metabolites during nitrogen fixation in Medicago via high resolution MALDI-MS imaging and ESI-MS J. Am. Soc. Mass Spectrom. 2015 26 149 158 10.1007/s13361-014-1010-0 25323862 111. Ye H. Gemperline E. Venkateshwaran M. Chen R. Delaux P.-M. Howes-Podoll M. Ané J.-M. Li L. MALDI mass spectrometry-assisted molecular imaging of metabolites during nitrogen fixation in the Medicago truncatula–Sinorhizobium meliloti symbiosis Plant J. 2013 75 130 145 10.1111/tpj.12191 23551619 112. Veličković D. Agtuca B.J. Stopka S.A. Vertes A. Koppenaal D.W. Paša-Tolić L. Stacey G. Anderton C.R. Observed metabolic asymmetry within soybean root nodules reflects unexpected complexity in rhizobacteria-legume metabolite exchange ISME J. 2018 12 2335 2338 10.1038/s41396-018-0188-8 29899508 113. Agtuca B.J. Stopka S.A. Evans S. Samarah L. Liu Y. Xu D. Stacey M.G. Koppenaal D.W. Paša-Tolić L. Anderton C.R. Metabolomic profiling of wild-type and mutant soybean root nodules using laser-ablation electrospray ionization mass spectrometry reveals altered metabolism Plant J. 2020 103 1937 1958 10.1111/tpj.14815 32410239 114. Boughton B.A. Thinagaran D. Sarabia D. Bacic A. Roessner U. Mass spectrometry imaging for plant biology: A review Phytochem. Rev. 2016 15 445 488 10.1007/s11101-015-9440-2 27340381 115. Tarkowská D. Novák O. Oklestkova J. Strnad M. The determination of 22 natural brassinosteroids in a minute sample of plant tissue by UHPLC–ESI–MS/MS Anal. Bioanal. Chem. 2016 408 6799 6812 10.1007/s00216-016-9807-2 27531032 116. Shimizu T. Miyakawa S. Esaki T. Mizuno H. Masujima T. Koshiba T. Seo M. Live Single-Cell Plant Hormone Analysis by Video-Mass Spectrometry Plant Cell Physiol. 2015 56 1287 1296 10.1093/pcp/pcv042 25759328 117. Martin C. Bhatt K. Baumann K. Shaping in plant cells Curr. Opin. Plant Biol. 2001 4 540 549 10.1016/S1369-5266(00)00213-2 11641071 118. Gago J. Daloso D.d.M. Figueroa C.M. Flexas J. Fernie A.R. Nikoloski Z. Relationships of Leaf Net Photosynthesis, Stomatal Conductance, and Mesophyll Conductance to Primary Metabolism: A Multispecies Meta-Analysis Approach Plant Physiol. 2016 171 265 279 10.1104/pp.15.01660 26977088 119. Wen W. Li D. Li X. Gao Y. Li W. Li H. Liu J. Liu H. Chen W. Luo J. Metabolome-based genome-wide association study of maize kernel leads to novel biochemical insights Nat. Commun. 2014 5 3438 10.1038/ncomms4438 24633423 120. Hu C. Shi J. Quan S. Cui B. Kleessen S. Nikoloski Z. Tohge T. Alexander D. Guo L. Lin H. Metabolic variation between japonica and indica rice cultivars as revealed by non-targeted metabolomics Sci. Rep. 2014 4 5067 10.1038/srep05067 24861081 121. Robaina-Estévez S. Daloso D.M. Zhang Y. Fernie A.R. Nikoloski Z. Resolving the central metabolism of Arabidopsis guard cells Sci. Rep. 2017 7 8307 10.1038/s41598-017-07132-9 28814793 122. Töpfer N. Braam T. Shameer S. Ratcliffe R.G. Sweetlove L.J. Alternative CAM Modes Provide Environment-Specific Water-Saving Benefits in a Leaf Metabolic Model Plant Cell 2020 10.1105/tpc.20.00132 123. Kim-Hellmuth S. Aguet F. Oliva M. Muñoz-Aguirre M. Kasela S. Wucher V. Castel S.E. Hamel A.R. Viñuela A. Roberts A.L. Cell type–specific genetic regulation of gene expression across human tissues Science 2020 369 eaaz8528 10.1126/science.aaz8528 32913075 124. Shrestha B. Ten Major Future Challenges in Single-Cell Metabolomics Single Cell Metabolism: Methods and Protocols Shrestha B. Springer New York, NY, USA 2020 219 223 125. Macaulay I.C. Ponting C.P. Voet T. Single-Cell Multiomics: Multiple Measurements from Single Cells Trends Genet. 2017 33 155 168 10.1016/j.tig.2016.12.003 28089370