Yonghui
Dong
and
Asaph
Aharoni
*
Department of Plant Sciences, Weizmann Institute of Science, Rehovot 76100, Israel. E-mail: asaph.aharoni@weizmann.ac.il
First published on 23rd June 2022
Covering: 2017 to 2022
Mass spectrometry imaging (MSI) has become a mature molecular imaging technique that is well-matched for natural product (NP) discovery. Here we present a brief overview of MSI, followed by a thorough discussion of different MSI applications in NP research. This review will mainly focus on the recent progress of MSI in plants and microorganisms as they are the main producers of NPs. Specifically, the opportunity and potential of combining MSI with other imaging modalities and stable isotope labeling are discussed. Throughout, we focus on both the strengths and weaknesses of MSI, with an eye on future improvements that are necessary for the progression of MSI toward routine NP studies. Finally, we discuss new areas of research, future perspectives, and the overall direction that the field may take in the years to come.
Mass spectrometry imaging (MSI) is a well-established analytical tool that allows direct mapping of a wide variety of chemical classes from different biological samples, providing information concerning analyte identity, relative abundance, and spatial distribution. MSI has enjoyed a surge in popularity in the last decade due to its non-targeted and label-free nature. Analytes of interest do not have to be pre-selected or even known prior to MSI analysis, and they can in most cases be detected without any chemical modification or labelling. This is in stark contrast to most histochemical staining techniques, which require the availability of a suitable antibody.4 It is worth noting that multiplexed ion beam imaging (MIBI) is an exception that requires antibodies tagged with monoisotopic metal reporters to image targeted proteins using secondary ion mass spectrometry.5 Most importantly, MSI provides spatio-chemical information which is much more specific as compared to different kinds of microscopic imaging techniques, and much more intuitive similar to colorimetric imaging. These characteristics have made MSI particularly suited for exploratory analysis.
Several key topics such as MSI ionization technologies, sample preparation, and data analysis have been extensively reviewed elsewhere.6–10 This review particularly intends to highlight the emerging applications, overall trends, and recent innovations that expand the scope of MSI in NP studies. In addition, strategies, and opportunities for integrating MSI with other imaging modalities and stable isotope labelling are discussed. Finally, remaining challenges and future promises are presented.
![]() | ||
Fig. 1 Overview of the Mass Spectrometry Imaging (MSI) technique. (A) General workflow of a typical MSI procedure. (1) A thin-section sample is mounted onto a glass slide, coated with or without matrix. By using a laser (MALDI), ion beam (SIMS), or charged solvent stream (DESI), ions are generated at each spot (pixel) which are then detected by a mass spectrometer. (2) Mass spectrum is produced from each spot (pixel). (3) After acquisition, ion intensity of each unique, resolved m/z peak can be extracted from each spot and plotted as a scaled false color heat map-like image, termed MS image or ion image. As exemplified here, single MS images for several different metabolites, such as phosphatidylcholine (PC) 36:4 and naringenin chalcone, are generated to represent their spatial intensity distributions in two different tomato fruit sections. In addition, MS images can be overlaid to compare their distributions. (B) Currently, MSI techniques need to compromise between the “4S-criteria” (i.e. speed, specificity, spatial resolution, and sensitivity) for performance. Priority should be given to each different criteria depending on the analytes of interest, sample tissue size, time constraints, and objective of the study. Panel B was reprinted (adapted) from Schulz et al. (2019)15 with permission. |
Spatial resolution, sensitivity, speed, and specificity are four important aspects (termed 4S-criteria) that must be rigorously considered before conducting an MSI experiment (Fig. 1B). Unfortunately, current MSI approaches need to compromise the 4S-criteria.15 As such, priority should be given to each different criteria depending on the compounds of interest, sample tissue size, time constraints, and objectives of the study. High spatial resolution MSI (e.g. <10 μm) is essential for studies designed to investigate biological differences at cellular level. While improving the spatial resolution (e.g. by the acquisition of smaller pixels) inevitably leads to a significant decrease in throughput and sensitivity as the result of the increased number of pixels and reduced sampling volume (or availability of ions) for each pixel (Fig. 1B). Conversely, low spatial resolution MSI can be acquired in a relatively short time scale and with high sensitivity, but it does not offer fine localization information.16 The development of ultra-high-speed MSI instrumentation, e.g. implementation of a high-repetition-rate laser in MALDI imaging, has allowed to compensate the elevated acquisition time of high resolution MSI. With the continuous improvement in scan speed, acquisition time and throughput will not become a significant issue for high resolution MSI in the future.17,18 High resolving power, accurate mass (e.g. <5 ppm) MSI is vital for specific MS image generation and reliable metabolite identification since it can distinguish targeted metabolites from isobaric ions and complex background noises. Although high resolving power is advantageous, it generally requires long spectrum acquisition time, thus leading to reduced measurement speed (Fig. 1B). High-speed MS imaging is typically used for phenotyping and screening purposes aiming at the analysis of a large number of samples. Therefore, spatial resolution and mass resolution are often sacrificed in favour of throughput. As noted above, designing an ideal MSI experiment is typically application-specific; the balance between the 4S-criteria for a particular experiment is very much dependent on the required outcomes, sample types, and the availability of the MSI instrument. Unlike many novices may think, high spatial resolution MSI is not always necessary to ensure a successful MSI experiment.
Reverse genetics strategies, employed to obtain a phenotype (e.g. changes in metabolite level) by disrupting or modifying a specific gene, are often used to study in vivo gene function. Although changes in metabolites can be quantified by assaying homogenized tissue using LCMS or GCMS, such methods are vulnerable to “bulk dilution effects” which may obscure the gene–metabolite relationship. As MSI enables measuring the spatial dependence of metabolic changes upon gene expression or silencing, it could provide more accurate tissue-specific, gene-to-metabolite correlation, from which significant clues regarding gene function involved in the biosynthesis of related metabolites may be obtained. Dong and co-authors demonstrated that coupling MALDI imaging and various reverse genetics approaches, such as RNA interference (RNAi), virus-induced gene silencing (VIGS), and agroinfiltration, provides a great opportunity in the gene function elucidation.22,23 In one example, they silenced the 2-oxoglutarate dependent dioxygenase GLYCOALKALOID METABOLISM31(GAME31) gene using anthocyanin-guided VIGS (Fig. 2A). The anthocyanin serves as a visual reporter that allowed them to confirm the silencing of the target gene and distinguish between silenced (anthocyanin-free) from non-silenced (anthocyanin-rich) tissue regions (Fig. 2B). Three anthocyanidins, i.e. delphinidin, petunidin, and malvidin, were detected by MALDI imaging, and their distributions were well correlated with the localization of the purple color in both control and GAME31-VIGS tomato fruit (Fig. 2C and D). In control fruit, esculeoside A was found homogeneously localized between silenced and non-silenced sectors, and only trace amounts of α-tomatine were detected (Fig. 2C). In contrast, a significant amount of α-tomatine was detected in GAME31-VIGS tomato fruit and it positively spatially co-occurred with the silenced sector. The distribution of esculeoside A was negatively co-occurred with the silenced sector in the GAME31-VIGS tomato fruit (Fig. 2D), confirming the role of GAME31 in converting α-tomatine to downstream esculeoside A during tomato fruit maturation.24,25 A similar experiment has been reported in studying the PHYTOENE DESATURASE (PDS) gene in maize by coupling MALDI imaging and Raman imaging with VIGS assays.26 Another noteworthy example involves the characterization of two homologous maize genes, GLOSSY2 (Gl2) and Gl2-like, whose specific biochemical functions remain elusive, but both genes are known to affect the accumulation of cuticular lipids.27 In this study, the authors measured and visualized the spatial and temporal patterns of cuticular lipid accumulation on the surface of developing Arabidopsis floral tissues expressing the two genes. MALDI images show that the expression of the maize Gl2 and Gl2-like in the Arabidopsis eceriferum2 (cer2) mutant background could restore the cuticular lipid accumulation on developing flowers surface, suggesting that maize Gl2 and Gl2-like are functional homologs of the Arabidopsis CER2 gene.27
![]() | ||
Fig. 2 Localization of anthocyanins and major steroidal glycoalkaloids (SGAs) in GAME31 Virus-Induced Gene Silencing (VIGS)-silenced ripe tomato fruit using MALDI imaging. (A) Tomato fruits expressing DELILA (Del) and ROSEA1 (Ros1) transcription factors under fruit-specific promoter E8 (E8:Del/Ros1 MT) were subject to VIGS assay using a Tobacco Rattle Virus (TRV) vector carrying fragments of the Del and Ros1 anthocyanin inducer genes and a GAME31 gene fragment. Consequently, the conversion of α-tomatine to its downstream hydroxy- and acetoxy-derivatives was blocked. (B) The non-uniform or patchy accumulation of anthocyanin red pigments distinguishes the non-silenced (marked 1: anthocyanin-rich) and silenced (marked 2: anthocyanin-free) sectors in the same tomato fruit. Scale bar, 100 μm. (C and D) MALDI images of delphinidin (2), petunidin (3), malvidin (4), α-tomatine (5), and esculeoside A (6) in control and GAME31-VIGS tomato fruit, respectively. Control fruit is the E8:Del/Ros tomato fruit infected with the TRV plasmid harbouring merely the Del/Ros1 sequence (without the GAME31 gene). It was used to assess if VIGS by itself could influence SGAs metabolism. As shown in C (5)-(6), the spatial accumulations of α-tomatine and esculeoside A were not correlated with the distribution of anthocyanins, demonstrating that VIGS does not disturb SGAs metabolism. Scale bar, 5 mm. Ac, acetoxy; Glu, glucose; Gal, galactose; Xyl, xylose. Sample sections were coated with 2,5-dihydroxybenzoic acid (DHB) and measured with MALDI imaging in positive ion mode (100 μm raster size). Reprinted (adapted) from Dong et al. (2020)22 with permission. |
Some metabolites accumulate in tissues and cells at or in the vicinity of their sites of biosynthesis. This active-site-directed approach can be used to identify the most representative tissue enriched with targeted metabolites; this tissue is then prioritized to search for the relevant genes. For instance, using MALDI imaging, Heskes et al. found that diterpenoids were exclusively localized to glandular trichomes on the surface of fruit and leaves of Vitex agnus-castus (Lamiaceae). Through analysis of the leaf trichome-specific transcriptome, the authors identified seven candidate genes involved in the diterpenoid biosynthesis.28 In another example, the spatial distribution of different lipid species in jojoba (Simmondsia chinensis) seeds was mapped using MALDI imaging.29 The images revealed heterogeneous distributions of different lipids in the seed tissue (Fig. 3A). For example, wax ester (WE) 42:2 was found enriched in the cotyledon (Fig. 3B) and triacylglycerol (TAG) 62:3 was highly concentrated in the embryonic axis of the seed (Fig. 3C). The overlaid MALDI images of the two lipids clearly demonstrated the differential enrichment of WEs and TAGs in jojoba seeds (Fig. 3D). Transcriptome analysis of the two distinct tissues showed that genes involved in WE synthesis, such as the jojoba fatty acyl-CoA elongase 1 (ScFAE1), fatty acyl-CoA reductase (ScFAR), and wax synthase (ScWS), were all highly expressed and predominant towards the cotyledon; while acyl-CoA: diacylglycerol O-acyltransferase 1 (ScDGAT1), a terminal enzyme in TAG synthesis, was found significantly expressed and biased toward the embryonic axis tissues (Fig. 3E).29
![]() | ||
Fig. 3 Heterogeneous metabolite distributions and gene expression in jojoba seeds. (A) Optical image of a longitude jojoba seed section. Scale bar, 1 mm. (B and C) MALDI images of wax ester (WE) 42:2 and triacylglycerol (TAG) 62:3. (D) Overlaid MALDI images of WE 42:2 (red) and TAG 62:3 (green), demonstrating the differential enrichment of WEs and TAGs in jojoba seeds. (E) Metabolic pathways and corresponding genes involved in WE and TAG biosynthesis in jojoba. The heatmap (blue to magenta color scale) next to each gene name was used to indicate gene expression bias toward embryonic axis tissues (blue) or cotyledon tissues (magenta). Acyl-CoA, acyl coenzyme A; GPAT, glycerol-3-phosphate acyltransferase; PAP, phosphatidic acid phosphatase. Sample sections were coated with 2,5-dihydroxybenzoic acid (DHB) and measured with MALDI imaging in positive ion mode (40 μm raster size). Reprinted (adapted) from Sturtevant et al. (2020)29 with permission. |
Plants and mushrooms are constantly exposed to various biotic and abiotic stresses that result in wounding. Wounding usually induces the expression of specific gene sets to acquire wound (site)-specific metabolites to properly respond to the injury and prevent infection. One prominent example of analyzing wound sites led to the identification of a POLYENE PIGMENT SYNTHASE 1 (PPS1), a previously uncharacterized type of polyketide synthases in the BY1 stereaceous mushroom.30 Upon mechanical wounding of the mushroom mycelium, two antilarval polyenes 18-methyl-19-oxoicosaoctaenoic acid and 20-methyl-21-oxodocosanonaenoic acid were found exclusively synthesized in the wound site, resulting in a yellow coloration of the injured area (Fig. 4A and B). The presence of double bonds and a methyl branch in the two polyenes suggested the existence of a highly reducing polyketide synthase. Homology searches in the BY1 mushroom genome eventually led to the identification of PPS1, and expression analysis showed that PPS1 transcript was massively upregulated following mycelial damage, confirming the role of PPS1 in the polyene biosynthesis.30
![]() | ||
Fig. 4 Identification of polyketide synthase (PPS1) in the damaged mycelium of the BY1 fungus. (A) Optical image of the mycelium of the BY1 fungus 72 h after injury, where polyenes 18-methyl-19-oxoicosaoctaenoic acid and 20-methyl-21-oxodocosanonaenoic acid result in yellow coloration of the injured area. (B) MALDI image of 20-methyl-21-oxodocosanonaenoic acid showing that it was mainly detected in the injured area. Samples were coated with a universal MALDI matrix (1![]() ![]() |
MSI, as a surface analytical technique, is particularly well-positioned for ex vivo enzyme activity characterization. The most promising way for high throughput analysis is through the covalent or non-covalent immobilization of the substrates onto an MS-compatible surface, followed by MSI analysis of the enzymatic reaction. Two variants of this approach are often used. One is to pre-immobilize the enzyme substrate on the surface and then monitor the enzyme activity; the other is to selectively capture the enzyme substrate or product onto the surface after in-solution enzymatic reaction.32,33 An important consideration of this approach is that the surface immobilization should be tolerant to aqueous washing, meanwhile, the products can be readily desorbed from the surface upon MSI analysis to produce adequate signals.32 The Mrksich group proposed a method named imaging self-assembled monolayers (SAM) for MALDI MS (known as iSAMDI-MS) to obtain quantitative, high-throughput maps of reaction progress.34 The SAM are formed on a gold surface through an Au–S bond, which is strong enough to permit repeated aqueous washing, yet weak enough to be desorbed by MALDI laser irradiation. To demonstrate the capability of this approach towards calculating the Michaelis constant (Km) of an enzyme–substrate reaction, the authors used iSAMDI to quantitate 2592 unique experiments performed inflow, where the products were captured onto a SAM that could then be imaged.34 In another study, a MALDI peptide array was fabricated to screen thrombin inhibitors by incorporation of bioactive peptides onto a modified surface of indium tin oxide (ITO) glass slide. The synthesized peptide substrates containing the thrombin cleavage site were tagged with dioctadecanoyl chains through the Diels–Alder reaction to amplify the MALDI signals and meanwhile to avoid the interference caused by MALDI matrix.33 The thrombin inhibitory activities and half-maximal inhibitory concentration (IC50) of proanthocyanidins with different degrees of polymerization (DP) were successfully determined with the fabricated peptide array coupled with MALDI imaging.33
MSI also allows assessing the spatial distribution of enzymatic activity directly from tissue sections. Enzyme histochemistry with MSI, namely “functional mass spectrometry imaging” (fMSI), is a novel field and is expected to resolve the current limitations of enzyme histochemistry. In one pioneering study, Hamilton et al. used fMSI to study the spatial distribution of venom activity across the snake venom gland.35 The venom gland of the brown forest cobra (Naja subfulva) is rich in enzymatically active phospholipases A2 (PLA2) and, by imaging the venom gland sections exposed to two different lipid substrates, the authors were able to map the spatial activity of PLA2 by visualizing both the substrates and the PLA2 products (Fig. 5). So far, fMSI has been mostly performed to investigate high molecular weight (MW) substrates such as peptides and lipids.36,37 While low-MS enzymatic transformation is also amenable to fMSI as demonstrated by Takeo et al. where the endogenous cholinesterase (ChE) activity was accurately determined in Mus musculus and Drosophila melanogaster by using deuterium-labelled acetylcholine and choline as substrates to distinguish from their endogenous metabolites.38 fMSI also has great potential to map multiple enzyme activities in tissue thanks to its untargeted nature.39
![]() | ||
Fig. 5 The use of functional MSI to study phospholipases A2 activity against two different lipid substrates, PC 16:1/16:1 and PC 15:0/18:1-d7, from N.subfulva venom gland. Optical images of venom gland tissue sections in the absence of lipid substrate (A), the addition of the substrate phosphatidylcholine (PC) 16:1/16:1 (D) and a deuterium-labeled substrate PC 15:0/18:1-d7 (G). (B, E and H) MALDI images of the substrate (blue) and their corresponding PLA2 product (yellow) for each section. (C, F and I) Averaged MALDI mass spectra in the absence of lipid substrate and addition of the two substrates. Scale bar, 2 mm. Sample sections were coated with α-cyano-4-hydroxycinnamic acid (CHCA) and measured with MALDI imaging in positive ion mode (80 μm raster size). Reprinted (adapted) from Hamilton et al. (2020)35 with permission. |
Highly localized metabolites with tiny yield may escape the detection from a standard LCMS or GCMS analysis workflow due to the “bulk dilution effect”. This phenomenon has been capitalized in a recent study concerning the identification of virulence factors causing mushroom soft rot disease.41 With the assistance of MALDI imaging, the authors were able to detect a toxin, namely Toxoflavin, which was produced by a bacterial mushroom pathogen Burkholderia gladioli primarily at the infected area of the mushroom (Fig. 6A and B). However, toxoflavin was not detected by LCMS analysis of the entire mushroom extract likely because its final concentration was below the limit of detection following global biomass extraction. By analyzing a ten-fold higher concentrated infected mushroom extract using LCMS the authors detected toxoflavin and confirmed the above assumption. In the same study, a trace amount of the antifungal agent caryoynencin was also detected by LCMS. However, its role in defense against mushroom soft rot disease could not be fully validated with either LCMS or MALDI imaging since caryoynencin is prone to degradation during sample preparation. The presence of highly polarizable units such as four conjugated C–C triple bonds in caryoynencin makes it ideally suitable for detection by Raman spectroscopy. As such, Raman imaging was used to monitor and compare the production of caryoynencin between healthy and infected mushroom areas (Fig. 6C). The result finally confirmed that caryoynencin was synthesized in infected tissues for defense against soft rot disease (Fig. 6D). This excellent example highlights the merits of combining additional orthogonal, chemically, or structurally informative imaging techniques with MSI to obtain more comprehensive evidence on metabolite structure and function. Multimodal imaging is further discussed in Section 4.1.
![]() | ||
Fig. 6 Multimodal identification of virulence factors involved in soft rot disease. (A) Workflow of multimodal imaging of infected mushroom slices. (B) MALDI image of toxoflavin showing that it primarily accumulates in infected mushroom tissue. (C) Density functional theory (DFT) calculated and measured UV Raman spectra of caryoynencin. Measurements were carried out on mushrooms slice: non-infected (black), infected (blue/red). (D) Box plots of the area under the curve of the Raman spectral band region associated with caryoynencin (2120–2220 cm−1) of infected mushroom tissue at 0, 8, 16, 24 h after infection. Sample sections were coated with a universal MALDI matrix (1![]() ![]() |
MSI can be used to investigate the metabolic secretome of various microbial species. It allowed discovering numerous secreted and colony-associated metabolites that are produced by microbes/plants to interact with or respond to neighbouring organisms and environmental stresses.42 For instance, acylsugars are specialized metabolites that are best known to be synthesized in the secretory glandular trichomes of the plant family Solanaceae (e.g. tomato), and evidence suggests that they play important roles in plant chemical defenses against insects and other herbivores.43 Korenblum et al. found that tomato root could also exudate acylsugars when colonized by Bacillus subtilis. MALDI images showed that different acylsugars appeared to accumulate in distinctive and highly localized root domains. For example, acylsucrose S1:5 was localized to the lateral root tips, while acylsucrose S4:19 was mainly distributed along the hairs of the main root.44 Additionally, MSI frequently allows predicting the bioactivity of metabolites based on their dynamic changes in spatial distribution; therefore it can be also used to filter the number of candidates needed to be purified for further characterization. For instance, Bai and co-authors found that the colony shape of a wilt-inducing fungus strain F. oxysporum f. sp. lycopersici 4287 (Fo4287) significantly changed when it was confronted with the strain P.polonicum ACCC31573 (Ppo1); while its colony morphology remained unaffected by the presence of the strain P. polonicum Hy4 (Ppo2) (Fig. 7A).45 Such observation suggested that Ppo1 could release certain metabolites to inhibit the growth of Fo4287. By direct analysing the microbial colony and the surrounding media using MALDI imaging, the authors detected 37 mass features in the negative ion mode. Three of them were predicted to possess antifungal activity based on their spatial distribution patterns using 3 criteria: (1) the candidate is secreted by Ppo1, (2) it could reach the edge of Fo4287 colony, and (3) it is not produced by Ppo2 (Fig. 7B). Subsequent LCMS, NMR, and anti-fungal activity assays of the 3 targeted candidates led to the identification of the novel fructigenine A and B anti-fungal compounds.45 In this example, instead of isolating each compound and testing its bioactivity, MALDI imaging was first applied to select the most promising antifungal candidates by their spatial distribution, therefore largely accelerating the antibiotic identification process.
![]() | ||
Fig. 7 Mapping the localization of candidate antifungal metabolites from Penicillium polonicum. (A) Optical images of the fungi strains grown on potato dextrose agar (PDA) plate. The first lane: a single colony of a wilt-inducing fungus strain F. oxysporum f. sp. lycopersici 4287 (Fo4287) and P. polonicum ACCC31573 (Ppo1), co-culture of Fo4287 and Ppo1 or P. polonicum Hy4 (Ppo2); the second lane: co-culture of Fo4287 and Ppo2 with different distances marked in the images. (B) MALDI images of three candidate m/z signals, i.e. m/z 281.25, m/z 410.25, and m/z 444.04, respectively. (C) Overlaid MALDI images and optical images of the colony shown in red rectangle in (A). Samples were coated with a universal MALDI matrix (1![]() ![]() |
Genome mining, the process of in silico identification of specific functional genes or biosynthetic gene clusters (BGCs) from genome sequence using modern bioinformatics tools, has revealed that many bacteria, fungi, and even higher organisms have far greater potential to produce NPs than what is estimated from classic bioactivity screening.46 This underestimation is partially due to the reason that many BGCs are silent or weakly expressed under standard laboratory conditions, resulting in low production yield which may prevent the detection of NPs by many analytical approaches.45 As such, there is a growing demand for novel genome mining methods and strategies to access these untapped NPs.47 Recently several approaches have been developed for activating silent BGCs, among which high-throughput elicitor screening (HiTES) is an attractive method for efficient detection of elicitors. In this approach, a reporter gene is inserted into a BGC of interest, and small molecule libraries are then screened against the reporter strain to look for suitable transcriptional elicitors. Once a hit is identified via phenotypic changes or assays, the small molecule product is further characterized along with the regulation of the silent BGC. Major drawbacks of this approach include: (i) the BGC must be targeted specifically using genetic information, which limits the throughput of NP discovery, especially for genetically intractable organisms; (ii) it does not provide a direct link of bioactivity to a given cryptic metabolite.48,49 To overcome these shortcomings, Xu et al. have adapted HiTES by replacing genetic reporter assays with MSI as a read-out for NPs production.50 This workflow consists of subjecting the wild-type microorganism to elicitor screening followed by imaging the resulting 500–1000 metabolomes, as a function of each molecule in the library, using LAESI imaging. With this approach, they discovered nine new metabolites in two different bacteria strains. This approach does not require challenging genetic, cloning, or culturing procedures, it can be thus used for both sequenced and un-sequenced bacteria.50 As MSI allows direct readout of the chemical output and metabolic exchange between organism and their environment, it holds great promise in speeding up NPs discovery rate and linking NPs to their biosynthetic gene clusters. This is especially true for organisms that are not genetically accessible.51 For instance, Helfrich and co-authors have isolated over 200 phyllosphere bacteria from Arabidopsis.52 The genomic analysis combined with thousands of binary interaction screens led to the identification of more than 1000 BGCs, with a majority lacking similarity with characterized BGCs. For bacteria harbouring a large number of BGCs and showing a wide inhibition spectrum, MALDI imaging was used to map the spatial distribution of ion peaks around colonies grown on different sensitive indicator strains. As discussed previously that MSI often allows predicting compounds with an antibiotic activity by correlating the shape of the inhibition zone around the colony with its spatial distribution. Metabolites of interest can be thus selected and purified guided by their distribution patterns. This strategy has led to the identification of two novel antibiotics, lysolipid phosphobrevin, and the trans-AT PKS-derived macrobrevin, from one selected strain Brevibacillus sp. Leaf182.52,53
![]() | ||
Fig. 8 The use of unsupervised segmentation to determine molecular phenotypes of four different E. coli colonies based on measured metabolomes. (A) Optical image of TY05, TY06, NHL17, and WT-TesA E. coli colonies. (B, E and F) Segment statistics and mass spectra highlight the primary features contributing to each strain segment. The segmentation results reveal the unique phenotypes of the four strains. The dark blue segment represents the dodecanoic acid producing TY05 strain and red corresponds to NHL17 engineered for octanoic acid production. The yellow and light blue segments denote the control TY06 strain and the broad free fatty acids producing WT-TesA strain. The inset in (F) is scale expanded segment mass spectra. (C and D) MALDI images of dodecanoic acid and octanoic acid, respectively. Samples measured by DESI imaging in positive ion mode with acetonitrile/water (90/10, with 0.5% NH4OH) being spray solvent (50 × 250 μm pixel size). Reprinted (adapted) from Ellis et al. (2021)55 with permission. |
Coupling MSI with different imaging modalities starts with optical microscopy techniques such as classical histology, fluorescence imaging, and immunohistochemistry. Optical microscopy is often performed on the same, or adjacent tissue section to acquire high-resolution images with detailed morphological and structural information. The acquired optical images can be used to determine the region of interest (ROI) for MSI analysis and interpret the MSI data by overlaying MS images with the optical image. The combination of these approaches could lead to the discovery of valuable biochemical information. For instance, Geier et al. have combined fluorescence in situ hybridization (FISH) microscopy and high-resolution atmospheric-pressure MALDI imaging (namely MetaFISH) to study the interactions between a deep-sea mussel and its intracellular symbiotic bacteria. This approach enables direct discovery of specialized metabolites located at the host–symbiont interface.63 Such method holds great promise for studying small scale interactions between the microbiome and the host. When identifying bioactive compounds associated with host–virus interactions (between the alga Emiliania huxleyi and its specific virus), Schleyer et al. found it difficult to accurately pinpoint a viral plaque (area without host cells, lysed by the virus) with naked eyes during the early viral infection.64 As chlorophyll degradation is known to be associated with viral infections, the authors used chlorophyll autofluorescence to identify such plaques. The centre of the plaque which mainly consists of lysed cells can be easily spotted by a significantly reduced chlorophyll signal (Fig. 9). The identified plaques were then subjected to MSI analysis, and a list of known lipid biomarkers involved in viral infection was detected (Fig. 9).64 Wisman et al. constructed a β-glucuronidase (GUS)-expressing strain of Aspergillus oryzae and used it to produce rice koji (a solid culture of Aspergillus oryzae on steamed rice grains). The image of the GUS reporter gene allowed them to clearly observe the fungal growth during koji making. The same grain was subjected to MSI analysis to visualize glucose distribution and correlate it with the degree of mycelial penetration.65
![]() | ||
Fig. 9 Multimodal imaging of known lipid biomarkers associated with host–viral interactions between the alga Emiliania huxleyi and the lytic E. huxleyi virus EhV20153. (A1 and B1) Chlorophyll autofluorescence of the sample (red pseudo-colour) accurately locates the centre of the viral plaque. (A2-8) MALDI images of virus-derived glycosphingolipids (vGSLs) t17:0/h22:0, betaine-like lipid (BLL) 36:6, BLL 38:6, triacylglycerol (TAG) 48:1, diacylglyceryl carboxyhydroxymethylcholine (DGCC) 40:7, glycerophospholipids, such as phosphatidyl-S,S-dimethylpropanethiol (PDPT) 40:7 and phosphatidylcholine (PC) 32:1. Scale bars, 1 mm. (C2-8) Flow-probe-MS imaging of the seven lipid biomarkers. Scale bars: 2 mm. Samples were imaged using both MALDI imaging and Flow-probe-MS imaging. For MALDI imaging, samples were coated with 2,5-dihydroxybenzoic acid (DHB) and measured in positive ion mode (100 μm raster size); for flow-probe-MS imaging, samples were measured in positive ion mode (25 × 500 μm pixel size). Reprinted (adapted) from Schleyer et al. (2019)64 with permission. |
Each MSI technique has its own merits and constraints. For instance, SIMS affords chemical imaging with micrometre to nanometre scale spatial resolution, which is the highest lateral resolution reported of any MSI methods. But it often leads to significant fragmentation of large molecules, which makes metabolite identification challenging. MALDI, by contrast, allows direct detection of intact molecules but the spatial resolution is limited. The spatial resolution for a given MALDI imaging experiment is mostly defined by the laser spot size. Current commercial MALDI imaging systems can now achieve spatial resolutions at 5–20 μm. The development of specialized optical setups, e.g. a transmission geometry optical configuration, has reduced the laser beam to less than 1 μm.66,67 However, other factors also affect the ultimate spatial resolution of a MALDI measurement, such as delocalization during sample preparation, the type of matrix used, the matrix application method and sensitivity of the instrument. As such, typical MALDI experiments are performed at pixel sizes of ∼30–100 μm.16,68,69 By coupling different MSI techniques, one can push the limits on the current “4S-criteria” constraint. For instance, Lanni et al. used a MALDI-guided SIMS imaging approach to study the compositional differences between the biofilms of WT and quorum sensing (QS) mutant Pseudomonas aeruginosa.70 The biofilms were first imaged by low-resolution MALDI imaging. The resulting MALDI images were then used to guide subsequent SIMS imaging at high resolution to examine the selected ROI.70 Multimodal MS imaging is also applied to boost molecular coverage. In a more recent study, Flint et al. used DESI imaging, imaging mass cytometry and laser ablation inductively coupled plasma (LA-ICP) imaging to characterize a 3D cell culture model of an epithelial HCC827 lung adenocarcinoma cell line.71 Combination of the three imaging techniques in this study enables the identification and mapping of the spatial distribution of metabolites, proteins, and elements, thus providing a more comprehensive understanding of the biological processes within the model.
MSI has also been frequently combined with other non-MS-based chemical imaging techniques for additional structural, morphological, biological, and chemical information. For instance, the capability to detect highly diverse molecules benefits from the use of different types of imaging modes, where one modality can detect molecules undetectable by another, as demonstrated in Fig. 6. Coupling MSI with vibrational spectroscopy imaging (VSI) techniques, such as Raman spectroscopy and Fourier Transform Infrared spectroscopy (FT-IR), provide perfect complementary toolsets to study complex biological samples. MSI offers specific molecular information while VSI adds information about the abundance and specificity of the chemical families (e.g. types of chemical bonds and heavy atoms).72 In particular, due to the non-destructive and label-free nature of VSI, the two approaches can be applied to the same sample. Susniak et al. combined MALDI imaging, FT-IR, and Raman spectroscopy to assess the effectiveness of biofertilizers by monitoring spatial changes of and relative abundance of several selected metabolites involved in nitrogen fixation in the pea root nodules.73 Magnetic resonance imaging (MRI) has also been used to complement MSI with detailed structural information. In one example, MRI and DESI imaging were used to study distributions of alkaloids in two structurally distinct regions of maturing areca nuts (i.e. seeds of Areca catechu). MRI images showed that during maturation, the liquid content started solidifying from the peripheral regions of the nut and during the solidification process, brown and white regions started forming. DESI images revealed that arecoline, arecaidine, and guvacoline alkaloids were segregated in the brown region while guvacine was localized in the white region.74 X-ray-based tomography, like micro-computed tomography (micro-CT) and phase-contrast synchrotron radiation based micro-computed tomography (PC-SRμCT), is a non-invasive imaging technique that is able to capture 3D structural information at micrometre-scale resolution. By integrating MALDI imaging with either PC-SRμCT75 or micro-CT,76 Geier et al. have developed a correlative imaging workflow to study the 3D microanatomy, cellular and chemical organization of host-microbe systems.
Significant challenges in terms of imaging data analysis are still pervasive. Most current multimodal workflows collect, process, and visualize the datasets separately. Developing a straightforward approach to integrating all datasets (such as image registration and data fusion) will enable deeper data mining. This is particularly difficult within an imaging regime because the various imaging modalities can have dramatically different spatial resolutions, data structures, and chemical information.61 Having a single hybrid instrument that combines multiple imaging methods would considerably simplify the experimental workflow and reduce the complexity associated with intermediary data processing and visualization steps.77 Zemaitis et al. coupled a 3D printed DESI ion source onto a commercial MALDI MS instrument.78 This approach allows the generation of multiple modes of MS images without the need to exchange ion sources and application of additional external software. In fact, an MS instrument equipped with multimodal MALDI and DESI imaging is already on the market, which delivers multi-layered, information-rich data from a single sample (e.g., Waters™). For multimodal data processing, universal imaging file formats and standard data processing and analysis workflows are in great demand. A unified data format, common to all spectral imaging techniques, would benefit and promote the development of multimodal imaging applications. Iakab et al. developed a software tool to convert Raman datasets to imzML format (the standard format for MSI files), thus allowing analysing Raman images using the same open-source software designed for MSI.79 Nevertheless, multimodal imaging will become more and more common in the future and provide unprecedented insights into NP studies.
One of the most common applications of global SIL in MSI is to determine the elemental composition (EC). The number of labelled atoms can be accurately inferred by calculating the mass difference between labelled and unlabelled peaks, thus effectively reducing the number of candidate ECs calculated by measured m/z value of the unlabelled peak within an accepted mass accuracy window. The use of multi-isotopes in global SIL, such as combining 13C and 15N labelling, can further narrow down the number of possible ECs since the number of both C and N atoms contained in the formula is accurately calculated. In addition, by feeding a specific stable isotope, it is possible to fish out all metabolites containing this element from a biological system. For instance, Kinnel et al. developed a MALDI-isotopic approach for NP discovery.85 In a proof-of-concept study, they analysed the filaments of M. producens JHB cultured with 14N- and 15N-nitrate as nitrogen sources using MALDI. By comparing the mass shift between the isotope cluster, the authors were able to identify a list of N-containing metabolites including a new metabolite cryptomaldamide.85 Due to the lack of retention time information, it is generally difficult to pair a labelled to its corresponding unlabelled peak accurately and efficiently. Consequently, MSI-based global SIL is often combined with LCMS for EC determination. For instance, Nakabayashi et al. integrated 15N labelling, LCMS, MALDI profiling, and MALDI imaging and successfully identified 65 N-containing metabolites in Catharanthus roseus.86 With the same approach, they also identified 45 N-containing monoterpene indole alkaloids in different tissues of Catharanthus roseus.87
Compared to global SIL, tracer-based SIL is more commonly used in MSI for various applications. A labelled precursor is usually administered to the organism to confirm the presence of certain precursor-derived metabolites detected by MSI according to their labelling patterns and spatial distribution.88,89 It can be also applied in an untargeted manner to potentially detect all precursor-derived metabolites. For instance, Feldberg et al. developed a novel approach that integrates in vivo dual tracer-based SIL with LCMS and MALDI imaging (termed DLEMMA-MALDI-imaging) for metabolite identification and metabolic network reconstruction.83 In the study, they used two differentially labelled tyrosine forms as feeding precursors and identified 60 tyrosine-derived metabolites in lemna (Spirodela polyrhiza) (Fig. 10A). Among them, half were detected due to the improved concentration as a result of tyrosine feeding, which highlights the advantage of tracer-based SIL for increasing molecular coverage. Structural isomers have identical m/z values, therefore most of the current MSI approaches are likely to reflect the summed distribution of multiple isomers. In DLEMMA-MALDI-imaging, instead of using the unlabelled m/z value, two labelled m/z values of the metabolites of interest were used to generate ion images, therefore minimizing the isomer interferences, and depicting the true spatial distribution of the metabolite. With this approach, the authors were able to construct a detailed tyrosine-derived metabolic network map in lemna (Fig. 10B and C). In another study, 1H NMR and a multimodal MS imaging approach (i.e., MALDI and NanoSIMS imaging) were applied collectively to examine the impact of severe drought on carbon allocation and metabolic responses within the roots of three tropical rainforest plants. The roots were labelled with position-specific (C1 or C2) 13C-pyruvate, which allows to precisely track the fate of 13C-pyruvate into downstream 13C labelled metabolites.90
![]() | ||
Fig. 10 Detection of tyrosine-derived metabolites in the water-floating plant lemna by using DLEMMA-MALDI-imaging. (A) Schematic representation of the DLEMMA-MALDI-imaging procedure. Three forms of tyrosine (one unlabeled and two differently labeled, i.e., tyrosine-d4 and tyrosine-13C915N1) were used as feeding precursors. After feeding for 24 h, lemna plants were attached to glass slides, desiccated, and subjected to MALDI imaging analysis. (B) Summary of the MALDI imaging results concerning identified tyrosine-derived metabolites. Thirteen metabolites were distributed in the whole plant (circled in red), and seven were primarily localized in roots (circled in blue). The letter ‘A’ denotes that the metabolite was detected in both tyrosine-fed and non-fed sample groups, and the letter ‘B’ means that the metabolite was only detected in the tyrosine-fed sample group. (C) Reconstruction of a putative tyrosine-derived metabolic network in Lemna based on the DLEMMA-MALDI-imaging study. The pathway was constructed based on possible biochemical reactions using 24 tyrosine-derived metabolites. Samples were coated with 2,5-dihydroxybenzoic acid (DHB) and measured with MALDI imaging in positive ion mode (65 μm raster size). Reprinted (adapted) from Feldberg et al. (2018)83 with permission. |
Coupling MSI with SIL also provides the opportunity for spatial flux analysis and kinetics studies. MSI is limited in its ability to determine the production and turnover rates of intermediates as it only captures a snapshot of metabolism at the time of sampling. When combined with SIL, the relative amounts of labelled versus unlabelled metabolites, and the degree of labelling within metabolites, could be used to define the rate of reactions and the fluxes within a given metabolic network. For instance, Romsdahl et al. labelled camelina sativa and pennycress (Thlaspi arvense) embryos with [U-13C]-glucose for 64 h and 120 h to study the spatial differences in carbon flux through phosphatidylcholine (PC).91 MALDI images showed that 13C carbon source was preferentially incorporated into PC molecular species in cotyledons compared to the embryonic axis, which may indicate a greater flux towards the cotyledons of developing embryos. The pattern of 13C-labeling also suggested that there is a greater 13C enrichment in more saturated PC species as compared to more unsaturated PC species, implying that desaturation is a relatively slower process than the incorporation of the acyl chains.91 In another example, SIL (13C and 15N) was coupled to large geometry-SIMS imaging (LG-SIMS) and NanoSIMS imaging to image the C and N2 fixation and cellular distribution of recently fixed C and N at several time points in the filaments of the heterocyst-forming cyanobacterium Anabaena sp. strain PCC 7120.92 A similar approach has been applied by this group to study nutrient kinetics in an algal-bacterial mutualism,93 N2 and C fixation activity of single cells from different diatom-Richelia symbiotic populations.94
The combination of SIL and MSI is still in its infancy. As the sample preparation method, labelling incorporation, and data analysis begin to develop and mature, it is certain to become a rapidly growing area in the near future, providing a window into a thus far hidden dimension of biochemical information.
The future of MSI relies on both innovation in instrumentation development and creativity in analytical methodologies. Coupling MSI with high-resolution, accurate MS, such as Fourier transform ion cyclotron resonance and Orbitrap MS, enables direct determination of EC from isotopic fine structures for small metabolites (<500 Da). The use of on-tissue tandem MS imaging (i.e. MS/MS imaging) significantly improves the confidence in metabolite identification. The development of automated, parallel MS and MS/MS imaging makes it possible to acquire hundreds of high-quality, unique MS/MS spectra without sacrificing spatial resolution and data acquisition time.95,100 As discussed earlier, the use of SIL also improves the confidence in metabolite identification. The development METASPACE (https://metaspace2020.eu/), an online annotation engine for MSI data, has largely facilitated metabolite identification. It performs MS1-based metabolite annotation by searching over ten public metabolomics/lipidomics databases such as HDMB101 and LIPID MAPS,102 and provides information regarding their molecular formulas and corresponding ion adducts, metabolite ion images, values of the metabolite-signal match (MSM) score, and an estimated false discovery rate (FDR) value for confidence-controlled metabolite annotation.103–105
Coupling MSI to more efficient desorption/ionization techniques such as laser-induced post-ionization (e.g. MALDI2) could improve the detection of poorly ionizable metabolites. MALDI2 is one of the most efficient means to enhance ion yields. With this technique, the initial MALDI-particle plume, formed in a fine vacuum environment of a few mbar of N2, is intercepted by the beam of a pulsed PI-laser which gives rise to a significant boost in ion yields for many metabolites.67,106 Chemical derivatization is another important approach to improve ionization efficiency and molecular coverage. In particular, the combination of several different complementary chemical derivatization agents allows the detection of various metabolites in an untargeted manner. To demonstrate this approach, Dueñas et al. (2019) performed on-tissue derivatization using three known chemical reactions: coniferyl aldehyde (CA) for primary amines, Girard's reagent T (GT) for carbonyl groups, and 2-picolylamine (2-PA) for carboxylic acids.96 This multiple on-tissue derivatization strategy was applied to explore the metabolite coverage in roots and leaves of two different maize genotypes (B73 and Mo17), and over six hundred new metabolite features were detected.96
As discussed earlier, MSI is unable to unambiguously resolve distinct isomers. This limitation could be addressed by coupling MSI with ion mobility spectrometry (IMS). IMS is a gas-phase separation technique that separates ions based on the form, size, and charge of the ions and has shown great potential in resolving the isobaric/isomeric metabolites.107 Tandem MS imaging is another widely used strategy to bypass isomeric interference, in which MS/MS imaging is performed to distinguish isomeric species in tissues, and the specific fragment is used to represent the spatial localization of the targeted metabolite.60,108 Multiple reaction monitoring (MRM) imaging (e.g., MRM-DESI imaging) can be regarded as a variant of tandem MS imaging, which has demonstrated improved sensitivity and specificity for the targeted metabolites.109 It is important to note that tandem MS imaging would be practicable only when the isomeric precursors generate different fragment ions. Unfortunately, some isomers are not easily discriminated by tandem mass spectrometry as a stand-alone methodology, such as glycan and lipid isomers.108 For instance, the positioning of carbon–carbon double bonds and stereospecific numbering are key factors to determine lipid functions, but when larger lipids are fragmented, this information is usually lost. Chemical derivatization of carbon–carbon double bonds has been employed to overcome this issue and facilitate specific fragmentation in MSI.110 Alternatively, MSI has been coupled with Paternò–Büchi reaction,111 epoxidation,112 ultraviolet photodissociation,113 ozone-induced dissociation,114 ion–ion reactions,115 and electron-induced dissociation116 to resolve lipid isomers.
The biosynthesis and storage of many NPs are highly compartmentalized. For instance, in plants, terpenes are produced in the disk cells and stored in the secretory cavity in glandular trichomes, with terpene precursors geranyl diphosphate and monoterpene formation occurring primarily in plastids, and farnesyl diphosphate and sesquiterpenes being synthesized predominantly in the cytosol.117 Measuring metabolic heterogeneity at the cellular or sub-cellular level thus has the potential to discover unforeseen details masked by averaging populations of cells and contribute to a deeper understanding of the metabolism and functions of NPs. For this reason, single-cell MSI is increasingly appreciated and various single-cell MSI techniques have been developed. Currently, SIMS is the predominant technique for single-cell MSI due to its sub-micron spatial resolution and surface sensitivity. MIBI, a reporter-based method using heavy metal tags coupled with SIMS, has allowed for targeted imaging analysis of proteins at the subcellular level.5,118 The development of optical fiber-based laser ablation electrospray ionization mass spectrometry (f-LAESI-MS),119–121 Nano-DESI,122 and laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS)123,124 have enabled MSI analysis of single cells under ambient conditions. Recent methodological and technological advances in MALDI imaging have improved its spatial resolution to 1–10 μm, making MALDI imaging at cellular and subcellular possible.125–127 However, sensitivity and speed are two major challenges (the 4S-criteria constraint), which need to be addressed to allow wide-scale application of single-cell MSI. Nevertheless, along with the further methodological and instrumental development, single-cell MSI holds significant promise in deepening our understanding of metabolic heterogeneity across cells.
With spatial resolution already approaching sub-micrometre levels and 3D MSI becoming more commonplace, the number of spectra acquired per experiment has increased exponentially. This data complexity gives rise to memory and computational challenges, namely “the curse of dimensionality”.128 Data analysis challenges are further exacerbated when integrating MSI with other imaging modalities. More powerful MSI-specific algorithms are in high demand for more efficient and effective data mining, visualization, image fusion, and interpretation. The reproducibility of MSI data analysis is another significant issue that needs to be addressed in the near future. The use of open-source software has done a remarkable job addressing this challenge but complete reporting of software versions and parameters is required, which is not yet routine in the current MSI studies.129 Improving sample preparation and workflows will also improve the quality and reproducibility of collected data. Data sharing is another critical aspect in MSI community because it offers opportunities for collaboration and promotes reproducibility. Previously, MSI data were archived in metabolomics and proteomics repositories, such as MetaboLights130 and PRIDE.131 The advent of OpenMSI132 and its extension BASTet133 have provided an important mean to share and reproduce MSI analyses. The construction and popularization of METASPACE has further accelerated to achieve data sharing, effective use, and verification of MSI result.104,105 It has become the largest MSI repository, which contains over 6000 MSI datasets. Users could easily share their MSI datasets and annotation result using a URL from METASPACE. The development of such tool represents an important direction towards establishing repositories for storing, sharing, and comparing MSI results.
This journal is © The Royal Society of Chemistry 2022 |