Spatially targeted proteomics analyzes the proteome of specific cell types and functional regions within tissue. While spatial context is often essential to understanding biological processes, interpreting sub-region-specific protein profiles can pose a challenge due to the high-dimensional nature of the data. Here, we develop a multivariate approach for rapid exploration of differential protein profiles acquired from distinct tissue regions and apply it to analyze a published spatially targeted proteomics data set collected from Staphylococcus aureus-infected murine kidney, 4 and 10 days postinfection. The data analysis process rapidly filters high-dimensional proteomic data to reveal relevant differentiating species among hundreds to thousands of measured molecules. We employ principal component analysis (PCA) for dimensionality reduction of protein profiles measured by microliquid extraction surface analysis mass spectrometry. Subsequently, k-means clustering of the PCA-processed data groups samples by chemical similarity. Cluster center interpretation revealed a subset of proteins that differentiate between spatial regions of infection over two time points. These proteins appear involved in tricarboxylic acid metabolomic pathways, calcium-dependent processes, and cytoskeletal organization. Gene ontology analysis further uncovered relationships to tissue damage/repair and calcium-related defense mechanisms. Applying our analysis in infectious disease highlighted differential proteomic changes across abscess regions over time, reflecting the dynamic nature of host-pathogen interactions.

@en

Imaging mass spectrometry (IMS) provides untargeted, highly multiplexed maps of molecular distributions in tissue. Ion images are routinely presented as heatmaps and can be overlaid onto complementary microscopy images that provide greater context. However, heatmaps use transparency blending to visualize both images, obscuring subtle quantitative differences and distribution gradients. Here, we developed a contour mapping approach that combines information from IMS ion intensity distributions with that of stained microscopy. As a case study, we applied this approach to imaging data from Staphylococcus aureus-infected murine kidney. In a univariate, or single molecular species, use-case of the contour map representation of IMS data, certain lipids colocalizing with regions of infection were selected using Pearson’s correlation coefficient. Contour maps of these lipids overlaid with stained microscopy showed enhanced visualization of lipid distributions and spatial gradients in and around the bacterial abscess as compared to traditional heatmaps. The full IMS data set comprising hundreds of individual ion images was then grouped into a smaller subset of representative patterns using non-negative matrix factorization (NMF). Contour maps of these multivariate NMF images revealed distinct molecular profiles of the major abscesses and surrounding immune response. This contour mapping workflow also enabled a molecular visualization of the transition zone at the host-pathogen interface, providing potential clues about the spatial molecular dynamics beyond what histological staining alone provides. In summary, we developed a new IMS-based contour mapping approach to augment classical stained microscopy images, providing an enhanced and more interpretable visualization of IMS-microscopy multimodal molecular imaging data sets.

@en

Gangliosides are acidic glycosphingolipids, containing ceramide moieties and oligosaccharide chains with one or more sialic acid residue(s) and are highly diverse isomeric structures with distinct biological roles. Matrix-assisted laser desorption/ionization imaging mass spectrometry (MALDI IMS) enables the untargeted spatial analysis of gangliosides, among other biomolecules, directly from tissue sections. Integrating trapped ion mobility spectrometry with MALDI IMS allows for the analysis of isomeric lipid structures in situ. Here, we demonstrate the gas-phase separation and identification of disialoganglioside isomers GD1a and GD1b that differ in the position of a sialic acid residue, in multiple samples, including a standard mixture of both isomers, a biological extract, and directly from thin tissue sections. The unique spatial distributions of GD1a/b (d36:1) and GD1a/b (d38:1) isomers were determined in rat hippocampus and spinal cord tissue sections, demonstrating the ability to structurally characterize and spatially map gangliosides based on both the carbohydrate chain and ceramide moieties.

@en

Introduction: Although Staphylococcus aureus is the leading cause of biofilm-related infections, the lipidomic distributions within these biofilms is poorly understood. Here, lipidomic mapping of S. aureus biofilm cross-sections was performed to investigate heterogeneity between horizontal biofilm layers. Methods: S. aureus biofilms were grown statically, embedded in a mixture of carboxymethylcellulose/gelatin, and prepared for downstream matrix-assisted laser desorption/ionization imaging mass spectrometry (MALDI IMS). Trapped ion mobility spectrometry (TIMS) was also applied prior to mass analysis. Results: Implementation of TIMS led to a ∼ threefold increase in the number of lipid species detected. Washing biofilm samples with ammonium formate (150 mM) increased signal intensity for some bacterial lipids by as much as tenfold, with minimal disruption of the biofilm structure. MALDI TIMS IMS revealed that most lipids localize primarily to a single biofilm layer, and species from the same lipid class such as cardiolipins CL(57:0) – CL(66:0) display starkly different localizations, exhibiting between 1.5 and 6.3-fold intensity differences between layers (n = 3, p < 0.03). No horizontal layers were observed within biofilms grown anaerobically, and lipids were distributed homogenously. Conclusions: High spatial resolution analysis of S. aureus biofilm cross-sections by MALDI TIMS IMS revealed stark lipidomic heterogeneity between horizontal S. aureus biofilm layers demonstrating that each layer was molecularly distinct. Finally, this workflow uncovered an absence of layers in biofilms grown under anaerobic conditions, possibly indicating that oxygen contributes to the observed heterogeneity under aerobic conditions. Future applications of this workflow to study spatially localized molecular responses to antimicrobials could provide new therapeutic strategies.

@en

Bacterial pathogens have evolved virulence factors to colonize, replicate, and disseminate within the vertebrate host. Although there is an expanding body of literature describing how bacterial pathogens regulate their virulence repertoire in response to environmental signals, it is challenging to directly visualize virulence response within the host tissue microenvironment. Multimodal imaging approaches enable visualization of host-pathogen molecular interactions. Here we demonstrate multimodal integration of high spatial resolution imaging mass spectrometry and microscopy to visualize Staphylococcus aureus envelope modifications within infected murine and human tissues. Data-driven image fusion of fluorescent bacterial reporters and matrix-assisted laser desorption/ionization Fourier transform ion cyclotron resonance imaging mass spectrometry uncovered S. aureus lysyl-phosphatidylglycerol lipids, localizing to select bacterial communities within infected tissue. Absence of lysyl-phosphatidylglycerols is associated with decreased pathogenicity during vertebrate colonization as these lipids provide protection against the innate immune system. The presence of distinct staphylococcal lysyl-phosphatidylglycerol distributions within murine and human infections suggests a heterogeneous, spatially oriented microbial response to host defenses.

@en

The search for molecular species that are differentially expressed between biological states is an important step towards discovering promising biomarker candidates. In imaging mass spectrometry (IMS), performing this search manually is often impractical due to the large size and high-dimensionality of IMS datasets. Instead, we propose an interpretable machine learning workflow that automatically identifies biomarker candidates by their mass-to-charge ratios, and that quantitatively estimates their relevance to recognizing a given biological class using Shapley additive explanations (SHAP). The task of biomarker candidate discovery is translated into a feature ranking problem: given a classification model that assigns pixels to different biological classes on the basis of their mass spectra, the molecular species that the model uses as features are ranked in descending order of relative predictive importance such that the top-ranking features have a higher likelihood of being useful biomarkers. Besides providing the user with an experiment-wide measure of a molecular species' biomarker potential, our workflow delivers spatially localized explanations of the classification model's decision-making process in the form of a novel representation called SHAP maps. SHAP maps deliver insight into the spatial specificity of biomarker candidates by highlighting in which regions of the tissue sample each feature provides discriminative information and in which regions it does not. SHAP maps also enable one to determine whether the relationship between a biomarker candidate and a biological state of interest is correlative or anticorrelative. Our automated approach to estimating a molecular species' potential for characterizing a user-provided biological class, combined with the untargeted and multiplexed nature of IMS, allows for the rapid screening of thousands of molecular species and the obtention of a broader biomarker candidate shortlist than would be possible through targeted manual assessment. Our biomarker candidate discovery workflow is demonstrated on mouse-pup and rat kidney case studies.

@en

We demonstrate the utility of combining silicon nanopost arrays (NAPA) and trapped ion mobility imaging mass spectrometry (TIMS IMS) for high spatial resolution and specificity mapping of neutral lipid classes in tissue. Ionization of neutral lipid species such as triglycerides (TGs), cholestryl esters (CEs), and hexosylceramides (HexCers) from biological tissues has remained a challenge for imaging applications. NAPA, a matrix-free laser desorption ionization substrate, provides enhanced ionization efficiency for the above-mentioned neutral lipid species, providing complementary lipid coverage to matrix-assisted laser desorption ionization (MALDI). The combination of NAPA and TIMS IMS enables imaging of neutral lipid species at 20 μm spatial resolution while also increasing molecular coverage greater than 2-fold using gas-phase ion mobility separations. This is a significant improvement with respect to sensitivity, specificity, and spatial resolution compared to previously reported imaging studies using NAPA alone. Improved specificity for neutral lipid analysis using TIMS IMS was shown using rat kidney tissue to separate TGs, CEs, HexCers, and phospholipids into distinct ion mobility trendlines. Further, this technology allowed for the separation of isomeric species, including mobility resolved isomers of Cer(d42:2) (m/z 686.585) with distinct spatial localizations measured in rat kidney tissue section.

@en

Low molecular weight metabolites are essential for defining the molecular phenotypes of cells. However, spatial metabolomics tools often lack the sensitivity, specify, and spatial resolution to provide comprehensive descriptions of these species in tissue. MALDI imaging mass spectrometry (IMS) of low molecular weight ions is particularly challenging as MALDI matrix clusters are often nominally isobaric with multiple metabolite ions, requiring high resolving power instrumentation or derivatization to circumvent this issue. An alternative to this is to perform ion mobility separation before ion detection, enabling the visualization of metabolites without the interference of matrix ions. Additional difficulties surrounding low weight metabolite visualization include high resolution imaging, while maintaining sufficient ion numbers for broad and representative analysis of the tissue chemical complement. Here, we use MALDI timsTOF IMS to image low molecular weight metabolites at higher spatial resolution than most metabolite MALDI IMS experiments (20 μm) while maintaining broad coverage within the human kidney. We demonstrate that trapped ion mobility spectrometry (TIMS) can resolve matrix peaks from metabolite signal and separate both isobaric and isomeric metabolites with different distributions within the kidney. The added ion mobility data dimension dramatically increased the peak capacity for spatial metabolomics experiments. Through this improved sensitivity, we have found >40 low molecular weight metabolites in human kidney tissue, such as argininic acid, acetylcarnitine, and choline that localize to the cortex, medulla, and renal pelvis, respectively. Future work will involve further exploring metabolomic profiles of human kidneys as a function of age, sex, and race.

@en

Identifying the spatial distributions of biomolecules in tissue is crucial for understanding integrated function. Imaging mass spectrometry (IMS) allows simultaneous mapping of thousands of biosynthetic products such as lipids but has needed a means of identifying specific cell-types or functional states to correlate with molecular localization. We report, here, advances starting from identity marking with a genetically encoded fluorophore. The fluorescence emission data were integrated with IMS data through multimodal image processing with advanced registration techniques and data-driven image fusion. In an unbiased analysis of spleens, this integrated technology enabled identification of ether lipid species preferentially enriched in germinal centers. We propose that this use of genetic marking for microanatomical regions of interest can be paired with molecular information from IMS for any tissue, cell-type, or activity state for which fluorescence is driven by a gene-tracking allele and ultimately with outputs of other means of spatial mapping.

@en

Metals play an essential role in biological systems and are required as structural or catalytic co-factors in many proteins. Disruption of the homeostatic control and/or spatial distributions of metals can lead to disease. Imaging technologies have been developed to visualize elemental distributions across a biological sample. Measurement of elemental distributions by imaging mass spectrometry and imaging X-ray fluorescence are increasingly employed with technologies that can assess histological features and molecular compositions. Data from several modalities can be interrogated as multimodal images to correlate morphological, elemental, and molecular properties. Elemental and molecular distributions have also been axially resolved to achieve three-dimensional volumes, dramatically increasing the biological information. In this review, we provide an overview of recent developments in the field of metal imaging with an emphasis on multimodal studies in two and three dimensions. We specifically highlight studies that present technological advancements and biological applications of how metal homeostasis affects human health.

@en

Lipids are a structurally diverse class of molecules with important biological functions including cellular signaling and energy storage. Matrix-assisted laser desorption/ionization (MALDI) imaging mass spectrometry (IMS) allows for direct mapping of biomolecules in tissues. Fully characterizing the structural diversity of lipids remains a challenge due to the presence of isobaric and isomeric species, which greatly complicates data interpretation when only m/z information is available. Integrating ion mobility separations aids in deconvoluting these complex mixtures and addressing the challenges of lipid IMS. Here, we demonstrate that a MALDI quadrupole time-of-flight (Q-TOF) mass spectrometer with trapped ion mobility spectrometry (TIMS) enables a >250% increase in the peak capacity during IMS experiments. MALDI TIMS-MS separation of lipid isomer standards, including sn backbone isomers, acyl chain isomers, and double-bond position and stereoisomers, is demonstrated. As a proof of concept, in situ separation and imaging of lipid isomers with distinct spatial distributions were performed using tissue sections from a whole-body mouse pup.

@en

Imaging mass spectrometry (IMS) is a rapidly advancing molecular imaging modality that can map the spatial distribution of molecules with high chemical specificity. IMS does not require prior tagging of molecular targets and is able to measure a large number of ions concurrently in a single experiment. While this makes it particularly suited for exploratory analysis, the large amount and high-dimensional nature of data generated by IMS techniques make automated computational analysis indispensable. Research into computational methods for IMS data has touched upon different aspects, including spectral preprocessing, data formats, dimensionality reduction, spatial registration, sample classification, differential analysis between IMS experiments, and data-driven fusion methods to extract patterns corroborated by both IMS and other imaging modalities. In this work, we review unsupervised machine learning methods for exploratory analysis of IMS data, with particular focus on (a) factorization, (b) clustering, and (c) manifold learning. To provide a view across the various IMS modalities, we have attempted to include examples from a range of approaches including matrix assisted laser desorption/ionization, desorption electrospray ionization, and secondary ion mass spectrometry-based IMS. This review aims to be an entry point for both (i) analytical chemists and mass spectrometry experts who want to explore computational techniques; and (ii) computer scientists and data mining specialists who want to enter the IMS field.

@en

Imaging mass spectrometry (IMS) enables the spatially targeted molecular assessment of biological tissues at cellular resolutions. New developments and technologies are essential for uncovering the molecular drivers of native physiological function and disease. Instrumentation must maximize spatial resolution, throughput, sensitivity, and specificity, because tissue imaging experiments consist of thousands to millions of pixels. Here, we report the development and application of a matrix-assisted laser desorption/ionization (MALDI) trapped ion-mobility spectrometry (TIMS) imaging platform. This prototype MALDI timsTOF instrument is capable of 10 μm spatial resolutions and 20 pixels/s throughput molecular imaging. The MALDI source utilizes a Bruker SmartBeam 3-D laser system that can generate a square burn pattern of <10 × 10 μm at the sample surface. General image performance was assessed using murine kidney and brain tissues and demonstrate that high-spatial-resolution imaging data can be generated rapidly with mass measurement errors <5 ppm and ∼40 000 resolving power. Initial TIMS-based imaging experiments were performed on whole-body mouse pup tissue demonstrating the separation of closely isobaric [PC(32:0) + Na]⁺ and [PC(34:3) + H]⁺ (3 mDa mass difference) in the gas phase. We have shown that the MALDI timsTOF platform can maintain reasonable data acquisition rates (>2 pixels/s) while providing the specificity necessary to differentiate components in complex mixtures of lipid adducts. The combination of high-spatial-resolution and throughput imaging capabilities with high-performance TIMS separations provides a uniquely tunable platform to address many challenges associated with advanced molecular imaging applications. ©

@en

The correlation of imaging mass spectrometry (IMS) with histopathology can help relate novel molecular findings obtained through IMS to the well-characterized and validated histopathology knowledge base. The quality of correlation between these two modalities is limited by the quality of the spatial mapping that is obtained by registration of the two image types. In this work, we develop novel workflows for MALDI IMS-to-microscopy data registration and analysis using nondestructive IMS-compatible wide field autofluorescence (AF) microscopy combined with computational image registration. First, a substantially automated procedure for high-accuracy registration between IMS and microscopy data of the same section is described that explicitly links the MALDI laser ablation pattern imaged by microscopy to its corresponding IMS pixel. Subsequent examination of the registered data allows for high-confidence colocalization of image features between the two modalities, down to single-cell scales within tissue. Building on this IMS-microscopy spatial mapping, we furthermore demonstrate the automated spatial correlation between IMS measurements from serial sections. This AF-registration-driven inter-section analysis, using a combination of nonlinear AF-to-AF and IMS-to-AF image registrations, can be applied to tissue sections that are prepared and imaged with different sample preparations (e.g., lipids vs proteins) and/or that are measured using different spatial resolutions. Importantly, all registrations, whether within a single section or across serial sections, are entirely independent of the IMS intensity signal content and thus unbiased by it.

@en

Diseases are characterized by distinct changes in tissue molecular distribution. Molecular analysis of intact tissues traditionally requires preexisting knowledge of, and reagents for, the targets of interest. Conversely, label-free discovery of disease-Associated tissue analytes requires destructive processing for downstream identification platforms. Tissue-based analyses therefore sacrifice discovery to gain spatial distribution of known targets or sacrifice tissue architecture for discovery of unknown targets. To overcome these obstacles, we developed a multimodality imaging platform for discovery-based molecular histology. We apply this platform to a model of disseminated infection triggered by the pathogen Staphylococcus aureus, leading to the discovery of infection-Associated alterations in the distribution and abundance of proteins and elements in tissue in mice. These data provide an unbiased, three-dimensional analysis of how disease affects the molecular architecture of complex tissues, enable culture-free diagnosis of infection through imaging-based detection of bacterial and host analytes, and reveal molecular heterogeneity at the host-pathogen interface.

@en

The molecular identification of species of interest is an important part of an imaging mass spectrometry (IMS) experiment. The high resolution accurate mass capabilities of Fourier transform ion cyclotron resonance mass spectrometry (FT-ICR MS) have recently been shown to facilitate the identification of proteins in matrix-assisted laser desorption/ionization (MALDI) IMS. However, these experiments are typically limited to proteins giving rise to ions of relatively low m/z due to difficulties transmitting and measuring large molecular weight ions of low charge states. Herein we have modified the source gas manifold of a commercial MALDI FT-ICR MS to regulate the gas flow and pressure to maximize the transmission of large m/z protein ions through the ion funnel region of the instrument. By minimizing the contribution of off-axis gas disruption to ion focusing and maximizing the effective potential wall confining the ions through pressure optimization, the signal-to-noise ratios (S/N) of most protein species were improved by roughly 1 order of magnitude compared to normal source conditions. These modifications enabled the detection of protein standards up to m/z 24 000 and the detection of proteins from tissue up to m/z 22 000 with good S/N, roughly doubling the mass range for which high quality protein ion images from rat brain and kidney tissue could be produced. Due to the long time-domain transients (>4 s) required to isotopically resolve high m/z proteins, we have used these data as part of an FT-ICR IMS-microscopy data-driven image fusion workflow to produce estimated protein images with both high mass and high spatial resolutions.

@en

Histology-directed imaging mass spectrometry (IMS) is a spatially targeted IMS acquisition method informed by expert annotation that provides rapid molecular characterization of select tissue structures. The expert annotations are usually determined on digital whole slide images of histological stains where the staining preparation is incompatible with optimal IMS preparation, necessitating serial sections: one for annotation, one for IMS. Registration is then used to align staining annotations onto the IMS tissue section. Herein, we report a next-generation histology-directed platform implementing IMS-compatible autofluorescence (AF) microscopy taken prior to any staining or IMS. The platform enables two histology-directed workflows, one that improves the registration process between two separate tissue sections using automated, computational monomodal AF-to-AF microscopy image registration, and a registration-free approach that utilizes AF directly to identify ROIs and acquire IMS on the same section. The registration approach is fully automated and delivers state of the art accuracy in histology-directed workflows for transfer of annotations (∼3-10 μm based on 4 organs from 2 species) while the direct AF approach is registration-free, allowing targeting of the finest structures visible by AF microscopy. We demonstrate the platform in biologically relevant case studies of liver stage malaria and human kidney disease with spatially targeted acquisition of sparsely distributed (composing less than one tenth of 1% of the tissue section area) malaria infected mouse hepatocytes and glomeruli in the human kidney case study.

@en

Imaging mass spectrometry (IMS) is a molecular imaging technology that can measure thousands of biomolecules concurrently without prior tagging, making it particularly suitable for exploratory research. However, the data size and dimensionality often makes thorough extraction of relevant information impractical. To help guide and accelerate IMS data analysis, we recently developed a framework that integrates IMS measurements with anatomical atlases, opening up opportunities for anatomy-driven exploration of IMS data. One example is the automated anatomical interpretation of ion images, where empirically measured ion distributions are automatically decomposed into their underlying anatomical structures. While offering significant potential, IMS-atlas integration has thus far been restricted to the Allen Mouse Brain Atlas (AMBA) and mouse brain samples. Here, we expand the applicability of this framework by extending towards new animal species and a new set of anatomical atlases retrieved from the Scalable Brain Atlas (SBA). Furthermore, as many SBA atlases are based on magnetic resonance imaging (MRI) data, a new registration pipeline was developed that enables direct non-rigid IMS-to-MRI registration. These developments are demonstrated on protein-focused FTICR IMS measurements from coronal brain sections of a Parkinson's disease (PD) rat model. The measurements are integrated with an MRI-based rat brain atlas from the SBA. The new rat-focused IMS-atlas integration is used to perform automated anatomical interpretation and to find differential ions between healthy and diseased tissue. IMS-atlas integration can serve as an important accelerator in IMS data exploration, and with these new developments it can now be applied to a wider variety of animal species and modalities. This article is part of a Special Issue entitled: MALDI Imaging, edited by Dr. Corinna Henkel and Prof. Peter Hoffmann.

@en

Lung cancer is the leading cause of cancer death. Beyond first line treatment, few therapeutic options are available, particularly for squamous cell carcinoma (SCC). Here, we have explored the phospholipidomes of 30 human SCCs and found that they almost invariably (in 96.7% of cases) contain phospholipids with longer acyl chains compared to matched normal tissues. This trait was confirmed using in situ 2D-imaging MS on tissue sections and by phospholipidomics of tumor and normal lung tissue of the L-Ikka^KA/KA mouse model of lung SCC. In both human and mouse, the increase in acyl chain length in cancer tissue was accompanied by significant changes in the expression of acyl chain elongases (ELOVLs). Functional screening of differentially expressed ELOVLs by selective gene knockdown in SCC cell lines followed by phospholipidomics revealed ELOVL6 as the main elongation enzyme responsible for acyl chain elongation in cancer cells. Interestingly, inhibition of ELOVL6 drastically reduced colony formation of multiple SCC cell lines in vitro and significantly attenuated their growth as xenografts in vivo in mouse models. These findings identify acyl chain elongation as one of the most common traits of lung SCC discovered so far and pinpoint ELOVL6 as a novel potential target for cancer intervention.

@en

Neurofibromatosis type 1 (NF1) is a common neurogenetic disorder, in which affected individuals develop tumors of the nervous system. Children with NF1 are particularly prone to brain tumors (gliomas) involving the optic pathway that can result in impaired vision. Since tumor formation and expansion requires a cooperative tumor microenvironment, it is important to identify the cellular and acellular components associated with glioma development and growth. In this study, we used 3-D matrix assisted laser desorption ionization imaging mass spectrometry (MALDI IMS) to measure the distributions of multiple molecular species throughout optic nerve tissue in mice with and without glioma, and to explore their spatial relationships within the 3-D volume of the optic nerve and chiasm. 3-D IMS studies often involve extensive workflows due to the high volume of sections required to generate high quality 3-D images. Herein, we present a workflow for 3-D data acquisition and volume reconstruction using mouse optic nerve tissue. The resulting 3-D IMS data yield both molecular similarities and differences between glioma-bearing and wild-type (WT) tissues, including protein distributions localizing to different anatomical subregions. Biological significance The current work addresses a number of challenges in 3-D MALDI IMS, driven by the small size of the mouse optic nerve and the need to maintain consistency across multiple 2-D IMS experiments. The 3-D IMS data yield both molecular similarities and differences between glioma-bearing and wild-type (WT) tissues, including protein distributions localizing to different anatomical subregions, which could then be targeted for identification and related back to the biology observed in gliomas of the optic nerve.@en