Mass Spectrometry Based Imaging of Labile Glucosides in Plants

Mass spectrometry based imaging is a powerful tool to investigate the spatial distribution of a broad range of metabolites across a variety of sample types. The recent developments in instrumentation and computing capabilities have increased the mass range, sensitivity and resolution and rendered sample preparation the limiting step for further improvements. Sample preparation involves sectioning and mounting followed by selection and application of matrix. In plant tissues, labile small molecules and specialized metabolites are subject to degradation upon mechanical disruption of plant tissues. In this study, the benefits of cryo-sectioning, stabilization of fragile tissues and optimal application of the matrix to improve the results from MALDI mass spectrometry imaging (MSI) is investigated with hydroxynitrile glucosides as the main experimental system. Denatured albumin proved an excellent agent for stabilizing fragile tissues such as Lotus japonicus leaves. In stem cross sections of Manihot esculenta, maintaining the samples frozen throughout the sectioning process and preparation of the samples by freeze drying enhanced the obtained signal intensity by twofold to fourfold. Deposition of the matrix by sublimation improved the spatial information obtained compared to spray. The imaging demonstrated that the cyanogenic glucosides (CNglcs) were localized in the vascular tissues in old stems of M. esculenta and in the periderm and vascular tissues of tubers. In MALDI mass spectrometry, the imaged compounds are solely identified by their m/z ratio. L. japonicus MG20 and the mutant cyd1 that is devoid of hydroxynitrile glucosides were used as negative controls to verify the assignment of the observed masses to linamarin, lotaustralin, and linamarin acid.

Introduction

Plants synthesize a diverse range of specialized metabolites to fend off herbivores and pests and to communicate with the environment (Møller, 2010; Neilson et al., 2013). In addition, some specialized metabolites like cyanogenic glucosides (CNglcs) may serve as storage and transport forms of reduced carbon and nitrogen used to fine tune primary metabolism (Jenrich et al., 2007; Picmanova et al., 2015; Nielsen et al., 2016; Bjarnholt et al., 2018). Understanding these interactions is key to developing robust crop plants for the future.

Mass spectrometry imaging (MSI) is a powerful technique to guide elucidation of the functional properties of specialized metabolites in nature by enabling determination of their precise cellular localization based on their molecular mass (Shroff et al., 2008; Hamm et al., 2010; Li B. et al., 2011; Sarsby et al., 2012; Li et al., 2013; Becker et al., 2014; Bjarnholt et al., 2014; Hölscher et al., 2014; Boughton et al., 2015; Shroff et al., 2015; Jarvis et al., 2017). Recent technological advances have provided increased sensitivity and improved mass and spatial resolution of their localization (Li et al., 2013; Boughton et al., 2015). Current, approaches to MSI for the spatial measurement of low mass natural products and their metabolites are dominated by Matrix Assisted Laser Desorption Ionisation (MALDI) (Boughton et al., 2015; Boughton and Hamilton, 2017; Weber et al., 2017). For small molecule analysis in complex matrices, the inability to distinguish between the presence of near-isobaric compounds with similar mass-to-charge (m/z) ratios provides a significant challenge for mass spectrometry (MS) based identification. Ultra-high mass resolution by MALDI – Fourier Transform Ion Cyclotron Resonance-MS (FTMS) provides significant advantages, using both high mass accuracy to identify the molecular formula (typically less than 2 ppm mass error) and ultra-high mass resolution to resolve metabolites with similar m/z (utilizing resolving powers >100,000).

Mapping of labile bio-active natural products to specific cells and tissue-types requires new instrumental and sample preparation methods. This is because the fixation and preparation processes may result in loss or structural degradation of the metabolites, either removing them entirely or decreasing their relative concentrations in the tissue or altering their cellular localisation by diffusion. Hydrolysis of glycosylated natural products during preparation for MALDI-MSI is attributed to the presence of glycosidases which during sectioning and processing are brought into contact with their otherwise compartmentalized substrates as demonstrated for many two component defense systems like those based on CNglcs (Morant et al., 2008). To stabilize tissues between the inherent time difference between sampling and preparation, tissues are generally flash-frozen and stored to retain the structure and distribution of endogenous metabolites prior to analysis. Frozen tissues require cryo-sectioning to generate thin tissue sections which are then generally thaw-mounted directly to a target suitable for analysis. For frozen plant samples containing CNglcs, the cryo-sectioning and freeze-thaw steps disrupt the cellular structures and results in mixing of the CNglcs with β-glucosidase enzymes leading to hydrolysis. This lowers signal intensity and impairs the spatial distribution analyses.

Here we present a novel broad based technique to embed, prepare, section, and mount plant tissue suitable for LDI and MALDI-MSI. The technique circumvents hydrolysis and diffusion of labile glycosides. The method was developed using Manihot esculenta (cassava) with the aim of being able to accurately monitor the distribution of CNglcs in different tissues at ultra-high mass resolution using FTMS. M. esculenta produces the two mono-CNglcs, linamarin (1) and lotaustralin (2), a number of cyanogenic di-glcs and apiosides as well as structurally related amides, acids (3) and anitriles (Picmanova et al., 2015; Figure 1), with the latter three being intermediates in the hypothesized endogenous turn-over pathways of CNglcs (Nielsen et al., 2016; Bjarnholt et al., 2018). A second CNglc producing species, Lotus japonicus, was used to verify the developed technique. This species produces both linamarin and lotaustralin and the non-cyanogenic hydroxynitrile glucosides rhodiocyanosides A (4) and D (5) (Bjarnholt and Møller, 2008). A non-cyanogenic L. japonicus mutant (cyd1) harboring a non-functional CYP79A1 enzyme was also included in the present study and provided an experimental negative control system (Takos et al., 2011). Whereas linamarin is the main CNglc present in cassava (Nartey, 1968; Koch et al., 1992; Lykkesfeldt and Møller, 1994), lotaustralin is the main CNglc in L. japonicus (Forslund et al., 2004).

FIGURE 1

FIGURE 1. Structures of hydroxynitrile glucosides and their endogenous turn-over products in Manihot esculenta (Me) and Lotus japonicus (Lj).

Cyanogenic glycosides (CNglcs) are specialized metabolites widely distributed throughout the plant kingdom (Gleadow and Møller, 2014). They are α-hydroxynitrile glycosides and as noted above form part of a two-component defense system detonated by a specific β-glucosidase that releases a hydrogen cyanide bomb upon disruption of the cellular structure, e.g., as mediated by a feeding herbivore (Figure 2; Morant et al., 2008). CNglcs are formed from amino acids which in the general pathway are converted into α-hydroxynitriles (cyanohydrins) by the action of two membrane bound cytochrome P450s and then via the action of glycosyltransferases converted into mono and di-glycosides (Conn, 1969).

FIGURE 2

FIGURE 2. Hydrolysis of the cyanogenic glucoside (CNglc) linamarin in M. esculenta and L. japonicus as a result of tissue disruption and the action of the β-glucosidase linamarase to form the α-hydroxynitrile, acetone cyanohydrin. Acetone cyanohydrin is dissociated into hydrogen cyanide and acetone by the action of an α-hydroxynitrile lyase. At alkaline pH, the latter process proceeds non-enzymatically.

Once considered to function only as components of plant chemical defense systems, it has become apparent that CNglcs contribute to enhance plant (Møller, 2010; Gleadow and Møller, 2014). They show species-specific chemical profiles and patterns of localization, both of which can change during plant ontogeny (Halkier and Møller, 1989; Adewusi, 1990; Koch et al., 1992; Nambisan and Sundaresan, 1994; Nielsen et al., 2016). Little is known about how plants balance biosynthesis, transport, storage, and turnover of CNglcs and how this is influenced by changing environmental conditions. Localisation of CNglcs and related metabolites with high spatial resolution within a range of tissues may be able to help shed light on some of these areas.

Results

Sectioning and Embedding Method

The preparation of plant tissues for mass spectrometric imaging required modifications when compared to standard procedures used with animal cells (Norris and Caprioli, 2013; Spengler, 2015). The rigid structure of the plant cell wall, the high intracellular water content and presence of large vacuoles render frozen sections of most plant tissues delicate to handle and prone to fracturing. To bypass this physical property, sections of a thickness above 35 μm were prepared to reduce breakage. For M. esculenta stem and tuber samples, section thicknesses of 40 and 60 μm, respectively, offered the best overall analytical details. In the case of L. japonicus leaves, a tissue thickness of 40 μm provided the best sections for analysis. Sections were mounted onto glass slides using either a standard ‘freeze-thaw’ approach or a ‘freeze-dry’ approach where frozen sections were transferred onto carbon double sided tape prior to freeze-drying.

Due to the relatively thin and delicate structural nature of leaves, it is not possible to section intact frozen leaves without prior imbedding. Sturdier plant tissues, such as the stems and tubers of M. esculenta plants used in this study may be sectioned without imbedding. In this study, we exploited the use of denatured albumin as an inexpensive and readily available embedding medium for fine tissue structures. The leaf material was sandwiched between two pieces of denatured albumin, then gently frozen over a surface of liquid nitrogen and sectioned. This approach offered leaf tissue sections suitable for MALDI-MSI analysis due to the very low background from the albumin, ease of handling and good sectioning properties (Kellersberger et al., 2011, 2012). Albumin contains very few small molecules that could interfere with any MALDI analysis, being made up of water (88–90%), protein (10–12%), and carbohydrate (0.24%) and containing almost no lipid.

Matrix Deposition Strategies Affect Distribution of Compounds

Matrix selection and deposition is an essential part of MALDI-MSI. 2,5-Dihydroxybenzoic acid (DHB) was chosen as matrix since it has been shown to work well with low molecular mass molecules such as nucleotides, peptides, lipids, and saccharides in positive ion mode (Zaima et al., 2010). Two deposition technologies were tested: Spray (wet) and sublimation (dry).

Sections were first mounted using a freeze-dry method and analyzed by MALDI-MSI. Initially, we focussed on measuring the distribution of linamarin reported as the most abundant CNglc in the M. esculenta tissue using wide mass range profiling (Nartey, 1968; Lykkesfeldt and Møller, 1994; Selmar, 1994). Linamarin was predominantly observed as the potassium salt [M+K]⁺ at m/z 286.06873 (calc. 286.06875, 0.05 ppm mass error) with a relatively low response compared with other endogenous metabolites (data not shown). The presence of the potassium salt agrees with previous findings (Kebarle, 2000; Kristensen et al., 2005). In the vacuole, the concentration of potassium ions may reach 500 mM (de Lacerda et al., 2003; Kanto et al., 2012). To increase sensitivity for linamarin, a method using Continuous Accumulation of Selected Ions (CASI) scan across a narrow m/z range was employed to filter highly abundant ions at different m/z. Vast differences were observed in linamarin response and localisation when comparing sections with DHB applied by either spray or sublimation (Figure 3). Specifically, results from the sprayed section showed a higher degree of tissue disruption and poor resolution of linamarin distribution, including very low signal (if any) in comparison to the signal obtained in equivalent areas of tissues when the matrix was applied by sublimation (Figure 3).

FIGURE 3

FIGURE 3. Distribution of linamarin ([M+K]+ m/z 286.06873) in 40 μm sections of a M. esculenta stem excised approximately 5 cm from the apex and covered with 2,5-dihydroxybenzoic acid (DHB) either by TM spray (A,B) and sublimation (C,D). Data was not normalized, Scale represents 0–100% signal intensity. Bar: 500 μm.

Sample Mounting Techniques Affect Metabolite Stability, Intensity, and Distribution

In M. esculenta, hydrolysis of the CNglcs linamarin and lotaustralin is catalyzed by linamarase, which is a highly stable β-glucosidase (Elias et al., 1997). To determine the stability and dislocation of CNglcs during the section mounting process, two different mounting techniques were compared, ‘freeze-thaw’ and ‘freeze-dry’ using sequential sections from the same tissue sample. For the mounting techniques, a standard freeze-thaw mounting to glass slide followed by vacuum drying was compared to mounting a frozen section to carbon double-sided tape followed by freeze-drying. In both treatments, sublimation with the matrix DHB was carried out before CNglc analysis. The experimental material used was M. esculenta stem tissues from young parts of the M. esculenta plant excised 5 cm from the apex and from older parts excised 5 cm above the soil surface. Images of freeze-dried vs freeze-thaw mounted sections on both young and older stem tissues were collected in the same analytical instrument sequence to compare the localisation of linamarin across the stem using the two different mounting techniques (Figure 4). Results showed that the standard freeze-thaw mounting approach led to an overall twofold to fourfold decrease in linamarin ion intensity when compared to the freeze-dried carbon tape based method in similar areas of tissue (Figure 5). Overall, the results indicate that a significant amount of the CNglc content is lost using the freeze-thaw mounting technique compromising to some extent the ability to accurately assess the distribution of linamarin in the tissue section.

FIGURE 4

FIGURE 4. The linamarin ([M+K]+ m/z 286.06873) signal obtained in 40 μm cross sections of M. esculenta stem excised 5 cm below apex (A,E) or 5 cm above soil surfaced (C,G) mounted either by freeze-dry (A–D) or freeze-thaw methodology (E–H). Scale represents 0–100% normalized signal intensity.

FIGURE 5

FIGURE 5. Box and dot plot of normalized linamarin ([M+K]+ m/z 286.06873) signal intensities in young and older 40 μm cross sections of a cassava stem mounted either by the freeze-thaw method (FT) or freeze-dry method (FD). Red dots represent outliers.

Distribution and Cyanogenic Compounds in M. esculenta

Manihot esculenta stem sections from the upper and lower part of the plant and from tubers proved to be useful proof of concept samples for visualizing the differential distribution of linamarin, lotaustralin [M+K]⁺ m/z 300.08435 (calc. 300.08439, 0.15 ppm error) and linamarin acid [M+K]⁺ m/z 305.06326 (calc. 305.06332, 2.16 ppm error). Linamarin and lotaustralin were observed to be present and evenly distributed throughout all tissue types of the young stem material except that the concentrations were somewhat reduced in the central core. In the sections from the older part of the stem, linamarin and lotaustralin were restricted to the first cell rows of the cortex, to cells around the phloem and to the parenchymatous cells surrounding the secondary xylem in the center of the stem (Figure 6). For tubers, linamarin and lotaustralin was localized to the periderm and vascular ring. Linamarin acid was found predominantly in the epidermis of stems and was detectable in the entire tuber but predominantly in the periderm. Cross-sections and longitudinal sections of old stems showed the same spatial distribution.

FIGURE 6

FIGURE 6. Distribution of linamarin ([M+K]+ m/z 286.06873) (A–D), lotaustralin ([M+K]+ m/z 300.08435) (E–H), and linamarin acid ([M+K]+ m/z 305.06326) (I–L) in 40 μm sections of young (5 cm from apex) and old (5 cm from soil surface) M. esculenta stem section and in 60 μm sections of tuber, all mounted using the freeze-dry method and sublimated with 2,5-DHB. Bars indicate 1 mm for all stem sections and 2 mm for all tuber sections. Scale represents 0–100% normalized signal intensity.

Testing of the Sample Preparation Technique and Observed Masses in Lotus japonicus

Once the stable MALDI-MSI preparation and sampling method, utilizing the denatured albumin embedding and freeze-dry method was established for M. esculenta, the techniques were tested on L. japonicus tissue. In the wild-type, CYP79D3 catalyzes the first committed step in CNglc biosynthesis in the leaves (Forslund et al., 2004). The cyd1 mutant of L. japonicus does not produce CNglcs because it contains a non-functional CYP79D3 (Takos et al., 2010). L. japonicus MG20 wildtype plants and cyd1 mutants were analyzed for ions expected to be the potassium adducts of lotaustralin, linamarin, and rhodiocyanoside A/D (observed [M+K]⁺ m/z 298.06944, calc. 298.06874, 2.33 ppm error). The signals for these adducts were observed in wild-type plants but were absent in the cyd1 mutant, thus linking these ions to the in vivo presence of their respective metabolites (Figures 7, 8). Linamarin, lotaustralin, and rhodiocyanoside A/D appeared to localize to the epidermal layers on both sides of the leaf. To further understand the role of the acid forms of CNglcs, the localisation of linamarin acid and lotaustralin acid were examined. An ion corresponding to the potassium adduct of linamarin acid [M+K]⁺ m/z 305.06519 (calc. 305.0639, 4.23 ppm mass error) was observed in both wild-type and cyd1. The respective signal for the potassium adduct of lotaustralin acid was not observed.

FIGURE 7

FIGURE 7. Presence of linamarin ([M+K]+ m/z 286.06873), lotaustralin ([M+K]+ m/z 300.08435), rhodiocyanoside A/D ([M+K]+ m/z 298.06944), and linamarin acid ([M+K]+ m/z 305.06326) in 40 μm cross sections of Lotus japonicus MG20 (A,C,E,G,I) and cyd1 (B,D,F,H,J) leaves mounted using the freeze-dry method and sublimated with 2,5-DHB. Scale represents 0–100% normalized signal intensity.

FIGURE 8

FIGURE 8. Box and dot plot of normalized signal intensities of linamarin ([M+K]⁺ m/z 286.06873), lotaustralin ([M+K]⁺ m/z 300.08435), rhodiocyanoside A/D ([M+K]⁺ m/z 298.06944), and linamarin acid ([M+K]⁺ m/z 305.06326) in 40 μm cross sections of L. japonicus MG20 and cyd1 leaves mounted using the freeze-dry method and sublimated with 2,5-DHB. Red dots represent outliers.

Discussion

The classically assigned physiological function of CNglcs is their action as defense compounds through cyanogenesis. The defense response is activated when otherwise compartmentalized CNglcs and hydrolytic enzymes are brought into contact by tissue disruption (Morant et al., 2008). More recently CNglcs have been attributed functions as transport and storage forms of reduced nitrogen and as ROS scavengers (Møller, 2010; Picmanova et al., 2015; Nielsen et al., 2016; Bjarnholt et al., 2018), resulting in the transient formation of the corresponding amides and carboxylic acids. Spatial analysis by MSI may guide the understanding of the orchestration of such functions at the tissue and cellular level and result in the discovery of hitherto unrecognized new functionalities.

The separation of CNglcs and hydrolytic enzymes is partly disrupted upon tissue sectioning. In an aqueous cellular environment, the β-glycosidases are active and may by simple diffusion processes gain contact with and hydrolyse CNglcs present. A recent report on MSI of CNglcs in flax (Linum usitatissimum) seed noted the potential for rapid hydrolysis of CNglcs but no measures to prevent this were provided (Dalisay et al., 2015). We observed that the mounting technique and method for application of matrix had a significant impact on signal intensity and localisation of CNglcs in M. esculenta tissues and we suggest keeping tissues frozen during all steps from sampling until matrix application (Figure 9).

FIGURE 9

FIGURE 9. Workflow diagram of preparation of tissue containing labile compounds for mass spectrometry based imaging.

Denatured albumin has previously been used in a number of different applications including quantitation directly from images (Kellersberger et al., 2011, 2012). A typical histological approach involves embedment of tissues in an Optimal Cutting Temperature (OCT) compound prior to freezing. OCT compounds have been shown to smear across the surface of tissues during cryo-sectioning and to significantly suppress analyte signal (Schwartz et al., 2003). Other embedding media including agar (Marques et al., 2014), gelatine (Ye et al., 2013; Horn et al., 2014; Korte and Lee, 2014), and carboxymethylcellulose (Yoshimura et al., 2012; Bencivenni et al., 2014) are available and have been employed for a variety of different tissue types. In general, aqueous gels or solutions have a tendency to spread or wet surrounding tissue slices when thawed, potentially reducing the spatial information obtained and inducing the hydrolysis of labile compounds. We found that handling the tissues while frozen and the inclusion of a freeze-dry step circumvented these negative effects. Inclusion of a freeze-dry step is recommended and would mitigate wetting if using a gel based embedding medium. More recently, methods using frozen gel embedding mediums and cryotape transfer have been demonstrated in barley and quinoa seeds but require specific cryotape (Gorzolka et al., 2016; Jarvis et al., 2017). The denatured albumin embedding, section and freeze-dry method described here presents a simple and low-cost alternative to the sample handling procedures described above.

In a single scan of a tissue sample, full range unbiased MSI offers an unprecedented account of the metabolic status with spatial resolution reaching the cellular level. In the event of low signal intensity for a desired analyte, the advanced features of the Bruker SolariX system enable targeted CASI experiments and trapping, increasing the sensitivity for individual metabolites. The CASI based method applies a quadrupole filter thus excluding competing ions from the detection cell. When coupled to high mass resolution, the approach provides excellent analyte signal intensity and capability to measure distribution of specific metabolites of interest. This was documented in the studies of the M. esculenta tissue samples, where the CASI based methods provided the best result by excluding competing ions from the detection cell increasing the sensitivity for the specific metabolites in question greatly, particularly when testing different mounting methods.

MALDI-MSI results for M. esculenta tubers demonstrated that linamarin and lotaustralin are localized in the periderm and vascular ring tissue (Figure 6). This pattern corresponds to earlier findings using DESI-MSI on fresh sections and orthogonal LC-MS on laser dissected tissues in which these CNglcs were also found to be localized in these regions (Li et al., 2013). Extending these previous finding, we show for the first time the localisation of the CNglcs linamarin and lotaustralin in stems and the presence of the endogenous turn-over product linamarin acid in stems and tubers of M. esculenta. Localisation of CNglc pathway transcripts and enzymes in M. esculenta stems have shown the pathway to localize to vascular tissue, cortex, and around the laticifers (Kannangara et al., 2011). The same localisation patterns were found for linamarin and lotaustralin in the data presented here (Figure 6). Linamarin acid is predominant in the sections obtained from the older part of the stem, which corresponds to the observations reported by Picmanova et al. (2015). Other intermediates of the proposed CNglc turnover pathway were not found but are generally present in very low concentrations (Picmanova et al., 2015), potentially below the limits of detection within this experiment. Glucosinolates—sulfur containing defense compounds found in Brassicaceae species—share similar initial biosynthetic steps and mode of bioactivation (Bak et al., 1998; Ishida et al., 2014). Localisation studies of glucosinolate biosynthesis and sulfur in stems of Arabidopsis thaliana linked both to specific cells in the vascular tissues (Koroleva et al., 2000; Li J. et al., 2011). The similar localisation patterns for glucosinolates and CNglcs in the vascular tissue might indicate they share even more functional parallels. The localisation of linamarin and lotaustralin to the first row of cortex cells and the vascular tissue in the older parts of the stem could offer protection against both generalist herbivores and phloem feeders.

In MSI, compounds are observed as their respective charged ions and are then first putatively annotated by the accurate mass match and the isotopic peak profile. Due to the possibility of many metabolites possessing a similar mass, an orthogonal technique, such as tandem-MS directly off tissue, LC-MS, NMR or the use of prior information are required to confirm identity (Boughton et al., 2015; McDonnell et al., 2015; Dias et al., 2016). The presence of linamarin, lotaustralin and linamarin acid have previously been identified using various orthogonal techniques in both M. esculenta and L. japonicus (including rhodiocyanosides A/D) (Forslund et al., 2004; Picmanova et al., 2015; Nielsen et al., 2016; Bjarnholt et al., 2018), but their spatial distributions have not been visualized using MALDI-MSI. In L. japonicus, both a cyanogenic wild-type and an acyanogenic mutant were available to test the sample preparation technique developed with M. esculenta and to ensure that the observed masses correspond to the CNglcs of interest and their related metabolites. A simple comparison between the two types can provide a clear statement (presence or absence) whether the observed mass indeed correspond to the metabolite of interest. One peculiar observation was the presence of linamarin acid in both wild-type and the acyanogenic cyd1 mutant. This indicates that linamarin acid may be formed by an additional route not related to endogenous turn-over of linamarin.

The presented method provides a broadly applicable preparation technique that enables imaging of plant metabolites using MALDI-MSI. For labile metabolites acted upon by hydrolysing enzymes it was found essential to halt enzymatic action by freeze-drying. Likewise, dry sublimation based matrix deposition significantly increased the metabolite signal above the signal obtained using wet spray based matrix deposition techniques. High mass accuracy, complementary LC-MS data and the availability and use of biological control samples like cyanogenic and acyanogenic cultivars provide independent coupling of m/z values observed to the known presence or absence of the labile metabolites. The advanced features of the MALDI-FT-ICR-MSI allowed targeted analysis of plant tissue metabolites with low abundancy, which would otherwise be difficult to detect. MSI is a potent new technology but cannot stand alone in structural identification of the constituents present in the samples analyzed. Control images obtained with reference compounds added to the dried tissue sample offers an alternative to the use of mutants. Appropriate sample preparation methodologies need to be adhered to in order to gain trustworthy results.

Experimental Procedures

Plant Material

Manihot esculenta (cassava) cultivar MAus7 (kindly provided by Professor Roslyn Gleadow, Monash University, Melbourne) was propagated from stem cuttings and grown in a greenhouse at a 30°C/16 h day and 20°C/8 h night regime. L. japonicus MG20 and cyd1 mutant seeds were provided by Dr. Daniela Lai and the plants grown in a greenhouse at 21°C day/16 h and 18°C night/8 h regime.

Tissue Preparation

Manihot esculenta stem (‘young’ tissue: 5 cm from apex, ‘old’ tissue: 5 cm above soil with stem diameters of approximately 4 and 8 mm) and tuber (with approximately diameter of 12 mm) tissues were harvested directly from 3 to 5 months old plants and immediately frozen in a 50 mL Falcon^® tube placed in liquid nitrogen. Commercial chicken eggs were boiled and peeled. The denatured egg white composed of albumin was collected. Suitably sized pieces were then cut to size for embedding the leaf sections. L. japonicus leaf samples were harvested directly from the plant and embedded in a 2 cm cube of denatured albumin (boiled egg-white). The albumin block was wrapped in cellophane and frozen gently by placing it directly above a liquid nitrogen surface. Tissues and the frozen blocks were mounted to a chuck using Tissue-Tek^® Optimal Cutting Temperature (O.C.T.) compound. Sections were cut at 35–60 μm thickness making sure no OCT compound came into contact with the sectioned tissue. For freeze-thaw mounting, sections were transferred onto pre-chilled (kept in the cryomicrotome at -20°C) Menzel–Gläser Superfrost Ultra-Plus Glass slides then gently warmed from the rear of the slide by pressing a finger until the glass was warm enough for the section to thaw and attach to the slide. For freeze-dried mounting, sections were transferred to pre-chilled glass slides with mounted carbon double-sided tape (ProSciTech, Thuringowa, QLD, Australia) and gently pressed onto the carbon double-sided tape using a fine paintbrush. All slides prepared with sections were kept frozen in a 50 mL Falcon^® tube placed in liquid nitrogen and were freeze-dried overnight. Mounted sections of both techniques were then stored in a vacuum-desiccator at room temperature.

Matrix Deposition

To compare the different matrix deposition strategies, mounted tissue sections were covered with 2,5-DHB (Sigma-Aldrich) using either wet (spray) or dry (sublimation). Wet matrix deposition was carried out using a HTX TM-Sprayer (HTX Industries, Chapel Hill, NC, United States) fitted with a Shimadzu LC20-AD HPLC pump (Shimadzu Australia, Rydalmere, NSW, Australia) using the following settings: Flow rate: 0.1 mL/min; Gas flow rate: 10 L/min; Nozzle temperature: 30°C; Spray conditions: 8 passes, 900 mm/min, spacing of 2 mm with alternate passes at a 90° offset and repeat passes set to an offset of 1 mm. DHB was dissolved in ethanol (100%) at a concentration of 10 mg/mL providing a matrix concentration of xx mg/mm². For dry application, DHB was sublimed onto tissue sections using a custom built sublimation apparatus at temperatures of 130–140°C at vacuum pressures less than 0.1 mBar for a period of 7–8 min to generate a matrix coverage of 0.3 ± 0.1 mg/cm².

Mass Spectrometer

For spatial mass spectrometric analysis, a Bruker (Bruker Daltonik, Bremen, Germany) SolariX XR 7 Tesla Hybrid ESI/MALDI-FT-ICR-MS was used. The instrument was operated in the positive ion mode using optimized instrumental settings. For full scan experiments, a mass range from m/z 100–800 was employed with analysis instrument set to broadband mode with a time domain for acquisition of 2M providing an estimated resolving power of 130,000 at m/z 400. For targeted analysis, CASI was used with a center mass set to m/z 286.1 (linamarin [M+K]⁺), m/z 300.1 (lotaustralin [M+K]⁺) and m/z 305.1 (linamarin acid [M+K]⁺) with 1 Da window chosen using a total of 1 ICR cell fills. The laser was set to 28 or 30% power using the minimum spot size resulting in ablation spots of approximately 35–40 μm in diameter. A total of 1,500–2,000 shots were fired per spectrum at a frequency of 2 kHz within a 50 μm × 50 μm array. Optical images of tissue sections were acquired using an Epson Photosmart 4480 flatbed scanner using a minimum setting of 4,800 d.p.i.

Data Analysis

Acquired mass spectrometry data were analyzed using Compass FlexImaging 4.1 (Bruker) and were normalized to Total Ion Chromatogram (TIC) and hierarchical cluster analysis conducted using a peak tolerance of 2 ppm mass error. Individual spectra were analyzed using Data Analysis 4.2 with peak lists generated using signal to noise threshold = 4 and 0.15% base peak height threshold. For quantification purposes similar areas of tissue (cell type and total area) were selected by generating specific regions of interest (ROI) using FlexImaging software. The ROI were exported as a spectrum list and imported into Data Analysis where an average spectrum was generated for each specific ROI. The average signal intensity of the respective m/z was then compared. Data were imported into SCiLS lab software with the respective ROIs and Box-dot plots were generated.

Author Contributions

FBS, BB, BLM, and KJ planned the project and the experiments. FBS, AH, KJ, and BB conceptualized and developed the methods. FBS collected and prepared the samples for analysis. FBS and BB analyzed the samples, interpreted the results, and drafted the manuscript. BB and DT did additional data analysis. FBS, BB, KJ, AH, and BLM commented on the manuscript for final version.

Funding

This research was supported by the VILLUM Research Center for Plant Plasticity; by the Center for Synthetic Biology (UCPH Excellence Programme for Interdisciplinary Research); and by an ERC Advanced Grant to BLM (ERC-2012-ADG_20120314). Imaging mass spectrometry was conducted at Metabolomics Australia located at the School of BioSciences, The University of Melbourne, Australia, which is an NCRIS initiative under Bioplatforms Australia Pty Ltd.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

Adewusi, S. R. A. (1990). Turnover of dhurrin in green sorghum seedlings. Plant Physiol. 94, 1219–1224. doi: 10.1104/pp.94.3.1219

CrossRef Full Text | Google Scholar

Bak, S., Kahn, R. A., Nielsen, H. L., Møller, B. L., and Halkier, B. A. (1998). Cloning of three A-type cytochromes p450, CYP71E1, CYP98, and CYP99 from Sorghum bicolor (L.) Moench by a PCR approach and identification by expression in Escherichia coli of CYP71E1 as a multifunctional cytochrome p450 in the biosynthesis of the cyanogenic glucoside dhurrin. Plant Mol. Biol. 36, 393–405. doi: 10.1023/A:1005915507497

PubMed Abstract | CrossRef Full Text | Google Scholar

Becker, L., Carre, V., Poutaraud, A., Merdinoglu, D., and Chaimbault, P. (2014). MALDI mass spectrometry imaging for the simultaneous location of resveratrol, pterostilbene and viniferins on grapevine leaves. Molecules 19, 10587–10600. doi: 10.3390/molecules190710587

PubMed Abstract | CrossRef Full Text | Google Scholar

Bencivenni, M., Faccini, A., Zecchi, R., Boscaro, F., Moneti, G., Dossena, A., et al. (2014). Electrospray MS and MALDI imaging show that non-specific lipid-transfer proteins (LTPs) in tomato are present as several isoforms and are concentrated in seeds. J. Mass Spectrom. 49, 1264–1271. doi: 10.1002/jms.3454

PubMed Abstract | CrossRef Full Text | Google Scholar

Bjarnholt, N., Li, B., D’Alvise, J., and Janfelt, C. (2014). Mass spectrometry imaging of plant metabolites - principles and possibilities. Nat. Prod. Rep. 31, 818–837. doi: 10.1039/C3NP70100J

PubMed Abstract | CrossRef Full Text | Google Scholar

Bjarnholt, N., and Møller, B. L. (2008). Hydroxynitrile glucosides. Phytochemistry 69, 1947–1961. doi: 10.1016/j.phytochem.2008.04.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Bjarnholt, N., Neilson, E. H., Crocoll, C., Jørgensen, K., Motawie, M., Olsen, C. E., et al. (2018). Glutathione transferases catalyze recycling of auto-toxic cyanogenic glucosides in sorghum. Plant J. 94, 1109–1125. doi: 10.1111/tpj.13923

PubMed Abstract | CrossRef Full Text | Google Scholar

Boughton, B. A., and Hamilton, B. (2017). “Spatial Metabolite profiling by matrix-assisted laser desorption ionization mass spectrometry imaging,” in Metabolomics: from Fundamentals to Clinical Applications, ed. A. Sussulini (Switzerland: Springer International Publishing), 291–321.

Google Scholar

Boughton, B. A., Thinagaran, D., Sarabia, D., Bacic, A., and Roessner, U. (2015). Mass spectrometry imaging for plant biology: a review. Phytochem. Rev. 15, 445–488. doi: 10.1007/s11101-015-9440-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Conn, E. E. (1969). Cyanogenic glycosides. J. Agric. Food Chem. 17, 519–526. doi: 10.1021/jf60163a014

CrossRef Full Text | Google Scholar

Dalisay, D. S., Kim, K. W., Lee, C., Yang, H., Rubel, O., Bowen, B. P., et al. (2015). Dirigent protein-mediated lignan and cyanogenic glucoside formation in flax seed: integrated omics and MALDI mass spectrometry imaging. J. Nat. Prod. 78, 1231–1242. doi: 10.1021/acs.jnatprod.5b00023

PubMed Abstract | CrossRef Full Text | Google Scholar

de Lacerda, C. F., Cambraia, J., Oliva, M. A., Ruiz, H. A., and Prisco, J. T. N. (2003). Solute accumulation and distribution during shoot and leaf development in two sorghum genotypes under salt stress. Environ. Exp. Bot. 49, 107–120. doi: 10.1016/S0098-8472(02)00064-3

CrossRef Full Text | Google Scholar

Dias, D., Jones, O., Beale, D., Boughton, B., Benheim, D., Kouremenos, K., et al. (2016). Current and future perspectives on the structural identification of small molecules in biological systems. Metabolites 6:46. doi: 10.3390/metabo6040046

PubMed Abstract | CrossRef Full Text | Google Scholar

Elias, M., Nambisan, B., and Sudhakaran, P. R. (1997). Isoforms of linamarase in cassava (Manihot esculenta). J. Plant Biochem. Biotechnol. 6, 63–67. doi: 10.1007/bf03263012

CrossRef Full Text | Google Scholar

Forslund, K., Morant, M., Jørgensen, B., Olsen, C. E., Asamizu, E., Sato, S., et al. (2004). Biosynthesis of the nitrile glucosides rhodiocyanoside A and D and the cyanogenic glucosides lotaustralin and linamarin in Lotus japonicus. Plant Physiol. 135, 71–84. doi: 10.1104/pp.103.038059

PubMed Abstract | CrossRef Full Text | Google Scholar

Gleadow, R. M., and Møller, B. L. (2014). “Cyanogenic glycosides: synthesis, physiology, and phenotypic plasticity,” in Annual Review of Plant Biology, Vol. 65, ed. S. S. Merchant (Palo Alto, CA: Annual Reviews), 155–185.

Google Scholar

Gorzolka, K., Kölling, J., Nattkemper, T. W., and Niehaus, K. (2016). Spatio-temporal metabolite profiling of the barley germination process by MALDI MS imaging. PLoS One 11:e0150208. doi: 10.1371/journal.pone.0150208

PubMed Abstract | CrossRef Full Text | Google Scholar

Halkier, B. A., and Møller, B. L. (1989). Biosynthesis of the cyanogenic glucoside dhurrin in seedlings of Sorghum bicolor (L.) Moench and partial purification of the enzyme system involved. Plant Physiol. 90, 1552–1559. doi: 10.1104/pp.90.4.1552

PubMed Abstract | CrossRef Full Text | Google Scholar

Hamm, G., Carre, V., Poutaraud, A., Maunit, B., Frache, G., Merdinoglu, D., et al. (2010). Determination and imaging of metabolites from Vitis vinifera leaves by laser desorption/ionisation time-of-flight mass spectrometry. Rapid Commun. Mass Spectrom. 24, 335–342. doi: 10.1002/rcm.4395

PubMed Abstract | CrossRef Full Text | Google Scholar

Hölscher, D., Dhakshinamoorthy, S., Alexandrov, T., Becker, M., Bretschneider, T., Buerkert, A., et al. (2014). Phenalenone-type phytoalexins mediate resistance of banana plants (Musa spp.) to the burrowing nematode Radopholus similis. Proc. Natl. Acad. Sci. U.S.A. 111, 105–110. doi: 10.1073/pnas.1314168110

PubMed Abstract | CrossRef Full Text | Google Scholar

Horn, P. J., Sturtevant, D., and Chapman, K. D. (2014). Modified oleic cottonseeds show altered content, composition and tissue-specific distribution of triacylglycerol molecular species. Biochimie 96, 28–36. doi: 10.1016/j.biochi.2013.08.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Ishida, M., Hara, M., Fukino, N., Kakizaki, T., and Morimitsu, Y. (2014). Glucosinolate metabolism, functionality and breeding for the improvement of Brassicaceae vegetables. Breed. Sci. 64, 48–59. doi: 10.1270/jsbbs.64.48

PubMed Abstract | CrossRef Full Text | Google Scholar

Jarvis, D. E., Ho, Y. S., Lightfoot, D. J., Schmöckel, S. M., Li, B., Borm, T. J. A., et al. (2017). The genome of Chenopodium quinoa. Nature 542, 307–312. doi: 10.1038/nature21370

PubMed Abstract | CrossRef Full Text | Google Scholar

Jenrich, R., Trompetter, I., Bak, S., Olsen, C. E., Moller, B. L., and Piotrowski, M. (2007). Evolution of heteromeric nitrilase complexes in Poaceae with new functions in nitrile metabolism. Proc. Natl. Acad. Sci. U.S.A. 104, 18848–18853. doi: 10.1073/pnas.0709315104

PubMed Abstract | CrossRef Full Text | Google Scholar

Kannangara, R., Motawia, M. S., Hansen, N. K. K., Paquette, S. M., Olsen, C. E., Møller, B. L., et al. (2011). Characterization and expression profile of two UDP-glucosyltransferases, UGT85K4 and UGT85K5, catalyzing the last step in cyanogenic glucoside biosynthesis in cassava. Plant J. 68, 287–301. doi: 10.1111/j.1365-313X.2011.04695.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Kanto, U., Jutamanee, K., Osotsapar, Y., and Jattupornpong, S. (2012). Effect of swine manure extract on leaf nitrogen concentration, chlorophyll content, total potassium in plant parts and starch content in fresh tuber yield of cassava. J. Plant Nutr. 35, 688–703. doi: 10.1080/01904167.2012.653074

CrossRef Full Text | Google Scholar

Kebarle, P. (2000). A brief overview of the present status of the mechanisms involved in electrospray mass spectrometry. J. Mass Spectrom. 35, 804–817. doi: 10.1002/1096-9888(200007)35:7<804::AID-JMS22>3.0.CO;2-Q

PubMed Abstract | CrossRef Full Text | Google Scholar

Kellersberger, K. A., Cornett, D. S., Becker, M., Schuerenberg, M., Belau, E., Groseclose, M. R., et al. (2011). Advancing the IMS Toolkit: A Model for the Optimization of Sample Preparation in Drug Imaging. Billerica, MA: BRUKER.

Kellersberger, K. A., Groseclose, M. R., Wagner, D. S., and Castellino, S. (2012). The Sausage Omelette: Breakfast of (IMS) Champions or a Novel Quantitation Tool for Small Molecule Imaging? Billerica, MA: BRUKER.

Koch, B., Nielsen, V. S., Halkier, B. A., Olsen, C. E., and Møller, B. L. (1992). The biosynthesis of cyanogenic glucosides in seedlings of cassava (Manihot esculenta Crantz). Arch. Biochem. Biophys. 292, 141–150. doi: 10.1016/0003-9861(92)90062-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Koroleva, O. A., Davies, A., Deeken, R., Thorpe, M. R., Tomos, A. D., and Hedrich, R. (2000). Identification of a new glucosinolate-rich cell type in Arabidopsis flower stalk. Plant Physiol. 124, 599–608. doi: 10.1104/pp.124.2.599

PubMed Abstract | CrossRef Full Text | Google Scholar

Korte, A. R., and Lee, Y. J. (2014). MALDI-MS analysis and imaging of small molecule metabolites with 1,5-diaminonaphthalene (DAN). J. Mass Spectrom. 49, 737–741. doi: 10.1002/jms.3400

PubMed Abstract | CrossRef Full Text | Google Scholar

Kristensen, C., Morant, M., Olsen, C. E., Ekstrøm, C. T., Galbraith, D. W., Møller, B. L., et al. (2005). Metabolic engineering of dhurrin in transgenic Arabidopsis plants with marginal inadvertent effects on the metabolome and transcriptome. Proc. Natl. Acad. Sci. U.S.A. 102, 1779–1784. doi: 10.1073/pnas.0409233102

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, B., Bjarnholt, N., Hansen, S. H., and Janfelt, C. (2011). Characterization of barley leaf tissue using direct and indirect desorption electrospray ionization imaging mass spectrometry. J. Mass Spectrom. 46, 1241–1246. doi: 10.1002/jms.2010

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, J., Kristiansen, K. A., Hansen, B. G., and Halkier, B. A. (2011). Cellular and subcellular localization of flavin-monooxygenases involved in glucosinolate biosynthesis. J. Exp. Bot. 62, 1337–1346. doi: 10.1093/jxb/erq369

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, B., Knudsen, C., Hansen, N. K., Jørgensen, K., Kannangara, R., Bak, S., et al. (2013). Visualizing metabolite distribution and enzymatic conversion in plant tissues by desorption electrospray ionization mass spectrometry imaging. Plant J. 74, 1059–1071. doi: 10.1111/tpj.12183

PubMed Abstract | CrossRef Full Text | Google Scholar

Lykkesfeldt, J., and Møller, B. L. (1994). Cyanogenic glycosides in cassava, Manihot esculenta Crantz. Acta Chem. Scand. 48, 178–180. doi: 10.3891/acta.chem.scand.48-0178

CrossRef Full Text | Google Scholar

Marques, J. V., Dalisay, D. S., Yang, H., Lee, C., Davin, L. B., and Lewis, N. G. (2014). A multi-omics strategy resolves the elusive nature of alkaloids in Podophyllum species. Mol. Biosyst. 10, 2838–2849. doi: 10.1039/c4mb00403e

PubMed Abstract | CrossRef Full Text | Google Scholar

McDonnell, L. A., Römpp, A., Balluff, B., Heeren, R. M. A., Albar, J. P., Andrén, P. E., et al. (2015). Discussion point: reporting guidelines for mass spectrometry imaging. Anal. Bioanal. Chem. 407, 2035–2045. doi: 10.1007/s00216-014-8322-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Møller, B. L. (2010). Functional diversifications of cyanogenic glucosides. Curr. Opin. Plant Biol. 13, 338–347. doi: 10.1016/j.pbi.2010.01.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Morant, A. V., Jørgensen, K., Jørgensen, C., Paquette, S. M., Sanchez-Perez, R., Møller, B. L., et al. (2008). b-glucosidases as detonators of plant chemical defense. Phytochemistry 69, 1795–1813. doi: 10.1016/j.phytochem.2008.03.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Nambisan, B., and Sundaresan, S. (1994). Distribution of linamarin and its metabolizing enzymes in cassava tissues. J. Sci. Food Agric. 66, 503–507. doi: 10.1002/jsfa.2740660413

CrossRef Full Text | Google Scholar

Nartey, F. (1968). Studies on cassava Manihot utilissima Pohl -.I. Cyanogenesis - biosynthesis of linamarin and lotaustralin in etiolated seedlings. Phytochemistry 7:1307. doi: 10.1016/S0031-9422(00)85629-0

CrossRef Full Text | Google Scholar

Neilson, E. H., Goodger, J. Q. D., Woodrow, I. E., and Møller, B. L. (2013). Plant chemical defense: at what cost? Trends Plant Sci. 18, 250–258. doi: 10.1016/j.tplants.2013.01.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Nielsen, L. J., Stuart, P., Pičmanová, M., Rasmussen, S., Olsen, C. E., Harholt, J., et al. (2016). Dhurrin metabolism in the developing grain of Sorghum bicolor (L.) Moench investigated by metabolite profiling and novel clustering analyses of time-resolved transcriptomic data. BMC Genomics 17:1021. doi: 10.1186/s12864-016-3360-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Norris, J. L., and Caprioli, R. M. (2013). Analysis of tissue specimens by matrix-assisted laser desorption/ionization imaging mass spectrometry in biological and clinical research. Chem. Rev. 113, 2309–2342. doi: 10.1021/cr3004295

PubMed Abstract | CrossRef Full Text | Google Scholar

Picmanova, M., Neilson, E. H., Motawia, M. S., Olsen, C. E., Agerbirk, N., Gray, C. J., et al. (2015). A recycling pathway for cyanogenic glycosides evidenced by the comparative metabolic profiling in three cyanogenic plant species. Biochem. J. 469, 375–389. doi: 10.1042/BJ20150390

PubMed Abstract | CrossRef Full Text | Google Scholar

Sarsby, J., Towers, M. W., Stain, C., Cramer, R., and Koroleva, O. A. (2012). Mass spectrometry imaging of glucosinolates in Arabidopsis flowers and siliques. Phytochemistry 77, 110–118. doi: 10.1016/j.phytochem.2012.01.026

PubMed Abstract | CrossRef Full Text | Google Scholar

Schwartz, S. A., Reyzer, M. L., and Caprioli, R. M. (2003). Direct tissue analysis using matrix-assisted laser desorption/ionization mass spectrometry: practical aspects of sample preparation. J. Mass Spectrom. 38, 699–708. doi: 10.1002/jms.505

PubMed Abstract | CrossRef Full Text | Google Scholar

Selmar, D. (1994). Translocation of cyanogenic glucosides in cassava. Acta Hortic. 375, 61–67. doi: 10.17660/ActaHortic.1994.375.3

CrossRef Full Text | Google Scholar

Shroff, R., Schramm, K., Jeschke, V., Nemes, P., Vertes, A., Gershenzon, J., et al. (2015). Quantification of plant surface metabolites by matrix-assisted laser desorption-ionization mass spectrometry imaging: glucosinolates on Arabidopsis thaliana leaves. Plant J. 81, 961–972. doi: 10.1111/tpj.12760

PubMed Abstract | CrossRef Full Text | Google Scholar

Shroff, R., Vergara, F., Muck, A., Svatos, A., and Gershenzon, J. (2008). Nonuniform distribution of glucosinolates in Arabidopsis thaliana leaves has important consequences for plant defense. Proc. Natl. Acad. Sci. U.S.A. 105, 6196–6201. doi: 10.1073/pnas.0711730105

PubMed Abstract | CrossRef Full Text | Google Scholar

Spengler, B. (2015). Mass spectrometry imaging of biomolecular information. Anal. Chem. 87, 64–82. doi: 10.1021/ac504543v

PubMed Abstract | CrossRef Full Text | Google Scholar

Takos, A., Lai, D., Mikkelsen, L., Abou Hachem, M., Shelton, D., Motawia, M. S., et al. (2010). Genetic screening identifies cyanogenesis-deficient mutants of Lotus japonicus and reveals enzymatic specificity in hydroxynitrile glucoside metabolism. Plant Cell 22, 1605–1619. doi: 10.1105/tpc.109.073502

PubMed Abstract | CrossRef Full Text | Google Scholar

Takos, A. M., Knudsen, C., Lai, D., Kannangara, R., Mikkelsen, L., Motawia, M. S., et al. (2011). Genomic clustering of cyanogenic glucoside biosynthetic genes aids their identification in Lotus japonicus and suggests the repeated evolution of this chemical defence pathway. Plant J. 68, 273–286. doi: 10.1111/j.1365-313X.2011.04685.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Weber, R. J. M., Lawson, T. N., Salek, R. M., Ebbels, T. M. D., Glen, R. C., Goodacre, R., et al. (2017). Computational tools and workflows in metabolomics: an international survey highlights the opportunity for harmonisation through Galaxy. Metabolomics 13:12. doi: 10.1007/s11306-016-1147-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Ye, H., Gemperline, E., Venkateshwaran, M., Chen, R., Delaux, P. M., Howes-Podoll, M., et al. (2013). MALDI mass spectrometry-assisted molecular imaging of metabolites during nitrogen fixation in the Medicago truncatula-Sinorhizobium meliloti symbiosis. Plant J. 75, 130–145. doi: 10.1111/tpj.12191

PubMed Abstract | CrossRef Full Text | Google Scholar

Yoshimura, Y., Zaima, N., Moriyama, T., and Kawamura, Y. (2012). Different localization patterns of anthocyanin species in the pericarp of black rice revealed by imaging mass spectrometry. PLoS One 7:e31285. doi: 10.1371/journal.pone.0031285

PubMed Abstract | CrossRef Full Text | Google Scholar

Zaima, N., Hayasaka, T., Goto-Inoue, N., and Setou, M. (2010). Matrix-assisted laser desorption/ionization imaging mass spectrometry. Int. J. Mol. Sci. 11, 5040–5055. doi: 10.3390/ijms11125040

PubMed Abstract | CrossRef Full Text | Google Scholar