Testing for presence of known and unknown molecules in imaging mass spectrometry

Benchmarking spatial clustering methods with spatially resolved transcriptomics data

Spatial clustering, which shares an analogy with single-cell clustering, has expanded the scope of tissue physiology studies from cell-centroid to structure-centroid with spatially resolved transcriptomics (SRT) data. Computational methods have undergone remarkable development in recent years, but a comprehensive benchmark study is still lacking. Here we present a benchmark study of 13 computational methods on 34 SRT data (7 datasets). The performance was evaluated on the basis of accuracy, spatial continuity, marker genes detection, scalability, and robustness. We found existing methods were complementary in terms of their performance and functionality, and we provide guidance for selecting appropriate methods for given scenarios. On testing additional 22 challenging datasets, we identified challenges in identifying noncontinuous spatial domains and limitations of existing methods, highlighting their inadequacies in handling recent large-scale tasks. Furthermore, with 145 simulated data, we examined the robustness of these methods against four different factors, and assessed the impact of pre- and postprocessing approaches. Our study offers a comprehensive evaluation of existing spatial clustering methods with SRT data, paving the way for future advancements in this rapidly evolving field.

NECTAR: A New Algorithm for Characterizing and Correcting Noise in QToF-Mass Spectrometry Imaging Data

A typical mass spectrometry imaging experiment yields a very high number of detected peaks, many of which are noise and thus unwanted. To select only peaks of interest, data preprocessing tasks are applied to raw data. A statistical study to characterize three types of noise in MSI QToF data (random, chemical, and background noise) is presented through NECTAR, a new NoisE CorrecTion AlgoRithm. Random noise is confirmed to be dominant at lower m/z values (∼50-400 Da) while systematic chemical noise dominates at higher m/z values (>400 Da). A statistical approach is presented to demonstrate that chemical noise can be corrected to reduce its presence by a factor of ∼3. Reducing this effect helps to determine a more reliable baseline in the spectrum and therefore a more reliable noise level. Peaks are classified according to their spatial S/N on the single ion images, and background noise is thus removed from the list of peaks of interest. This new algorithm was applied to MALDI and DESI QToF data generated from the analysis of a mouse pancreatic tissue section to demonstrate its applicability and ability to filter out these types of noise in a relevant data set. PCA and t-SNE multivariate analysis reviews of the top 4000 peaks and the final 744 and 299 denoised peak list for MALDI and DESI, respectively, suggests an effective removal of uninformative peaks and proper selection of relevant peaks.

MALDI MS Imaging of Chickpea Seeds (Cicer arietinum) and Crab's Eye Vine (Abrus precatorius) after Tryptic Digestion Allows Spatially Resolved Identification of Plant Proteins

Matrix-assisted laser desorption/ionization mass spectrometry (MALDI MS) imaging following in situ enzymatic digestion is a versatile analytical method for the untargeted investigation of protein distributions, which has rarely been used for plants so far. The present study describes a workflow for in situ tryptic digestion of plant seed tissue for MALDI MS imaging. Substantial modifications to the sample preparation procedure for mammalian tissues were necessary to cater to the specific properties of plant materials. For the first time, distributions of tryptic peptides were successfully visualized in plant tissue using MS imaging with accurate mass detection. Sixteen proteins were visualized and identified in chickpea seeds showing different distribution patterns, e.g., in the cotyledons, radicle, or testa. All tryptic peptides were detected with a mass resolution higher than 60,000 as well as a mass accuracy better than 1.5 ppm root-mean-square error and were matched to results from complementary liquid chromatography-MS/MS (LC-MS/MS) data. The developed method was also applied to crab's eye vine seeds for targeted MS imaging of the toxic protein abrin, showing the presence of abrin-a in all compartments. Abrin (59 kDa), as well as the majority of proteins visualized in chickpeas, was larger than 50 kDa and would thus not be readily accessible by top-down MS imaging. Since antibodies for plant proteins are often not readily available, in situ digestion MS imaging provides unique information, as it makes the distribution and identification of larger proteins in plant tissues accessible in an untargeted manner. This opens up new possibilities in the field of plant science as well as to assess the nutritional quality and/or safety of crops.

Supervised topological data analysis for MALDI mass spectrometry imaging applications

BACKGROUND

Matrix-assisted laser desorption/ionization mass spectrometry imaging (MALDI MSI) displays significant potential for applications in cancer research, especially in tumor typing and subtyping. Lung cancer is the primary cause of tumor-related deaths, where the most lethal entities are adenocarcinoma (ADC) and squamous cell carcinoma (SqCC). Distinguishing between these two common subtypes is crucial for therapy decisions and successful patient management.

RESULTS

We propose a new algebraic topological framework, which obtains intrinsic information from MALDI data and transforms it to reflect topological persistence. Our framework offers two main advantages. Firstly, topological persistence aids in distinguishing the signal from noise. Secondly, it compresses the MALDI data, saving storage space and optimizes computational time for subsequent classification tasks. We present an algorithm that efficiently implements our topological framework, relying on a single tuning parameter. Afterwards, logistic regression and random forest classifiers are employed on the extracted persistence features, thereby accomplishing an automated tumor (sub-)typing process. To demonstrate the competitiveness of our proposed framework, we conduct experiments on a real-world MALDI dataset using cross-validation. Furthermore, we showcase the effectiveness of the single denoising parameter by evaluating its performance on synthetic MALDI images with varying levels of noise.

CONCLUSION

Our empirical experiments demonstrate that the proposed algebraic topological framework successfully captures and leverages the intrinsic spectral information from MALDI data, leading to competitive results in classifying lung cancer subtypes. Moreover, the framework's ability to be fine-tuned for denoising highlights its versatility and potential for enhancing data analysis in MALDI applications.

Multiple multi-sample testing under arbitrary covariance dependency

Modern high-throughput biomedical devices routinely produce data on a large scale, and the analysis of high-dimensional datasets has become commonplace in biomedical studies. However, given thousands or tens of thousands of measured variables in these datasets, extracting meaningful features poses a challenge. In this article, we propose a procedure to evaluate the strength of the associations between a nominal (categorical) response variable and multiple features simultaneously. Specifically, we propose a framework of large-scale multiple testing under arbitrary correlation dependency among test statistics. First, marginal multinomial regressions are performed for each feature individually. Second, we use an approach of multiple marginal models for each baseline-category pair to establish asymptotic joint normality of the stacked vector of the marginal multinomial regression coefficients. Third, we estimate the (limiting) covariance matrix between the estimated coefficients from all marginal models. Finally, our approach approximates the realized false discovery proportion of a thresholding procedure for the marginal p-values for each baseline-category logit pair. The proposed approach offers a sensible trade-off between the expected numbers of true and false findings. Furthermore, we demonstrate a practical application of the method on hyperspectral imaging data. This dataset is obtained by a matrix-assisted laser desorption/ionization (MALDI) instrument. MALDI demonstrates tremendous potential for clinical diagnosis, particularly for cancer research. In our application, the nominal response categories represent cancer (sub-)types.

Systems Immunology Approaches to Metabolism

Over the last decade, immunometabolism has emerged as a novel interdisciplinary field of research and yielded significant fundamental insights into the regulation of immune responses. Multiple classical approaches to interrogate immunometabolism, including bulk metabolic profiling and analysis of metabolic regulator expression, paved the way to appreciating the physiological complexity of immunometabolic regulation in vivo. Studying immunometabolism at the systems level raised the need to transition towards the next-generation technology for metabolic profiling and analysis. Spatially resolved metabolic imaging and computational algorithms for multi-modal data integration are new approaches to connecting metabolism and immunity. In this review, we discuss recent studies that highlight the complex physiological interplay between immune responses and metabolism and give an overview of technological developments that bear the promise of capturing this complexity most directly and comprehensively.

Spatial segmentation of mass spectrometry imaging data featuring selected principal components

Spatial segmentation aims to find homogeneous/heterogeneous subgroups of spectra or ion images in mass spectrometry imaging (MSI) data. The maps it generated inform researchers of vital characteristics of the data and thus provide the basis for strategizing further biological analysis. Dimensional reduction and clustering are two basic steps of segmentation. Due to the variations in the quality, resolution, density of spectral information, and sizes, not all datasets could be segmented ideally with combinations of different dimensional reduction and clustering algorithms. Here, we proposed a segmentation pipeline that utilized pattern compression by principal component analysis (PCA) and represented by principal components. Instead of preprocessed or raw MSI data, normalized principal components were used for the segmentation process. Multiple datasets of rat brains and mouse kidneys were tested, and the proposed segmentation pipeline presented the obvious advantage of easy-to-use and can be readily intergraded with other existing innovative pipelines.

Spatially aware dimension reduction for spatial transcriptomics

Spatial transcriptomics are a collection of genomic technologies that have enabled transcriptomic profiling on tissues with spatial localization information. Analyzing spatial transcriptomic data is computationally challenging, as the data collected from various spatial transcriptomic technologies are often noisy and display substantial spatial correlation across tissue locations. Here, we develop a spatially-aware dimension reduction method, SpatialPCA, that can extract a low dimensional representation of the spatial transcriptomics data with biological signal and preserved spatial correlation structure, thus unlocking many existing computational tools previously developed in single-cell RNAseq studies for tailored analysis of spatial transcriptomics. We illustrate the benefits of SpatialPCA for spatial domain detection and explores its utility for trajectory inference on the tissue and for high-resolution spatial map construction. In the real data applications, SpatialPCA identifies key molecular and immunological signatures in a detected tumor surrounding microenvironment, including a tertiary lymphoid structure that shapes the gradual transcriptomic transition during tumorigenesis and metastasis. In addition, SpatialPCA detects the past neuronal developmental history that underlies the current transcriptomic landscape across tissue locations in the cortex.

What are we imaging? Software tools and experimental strategies for annotation and identification of small molecules in mass spectrometry imaging

Mass spectrometry imaging (MSI) has become a widespread analytical technique to perform nonlabeled spatial molecular identification. The Achilles' heel of MSI is the annotation and identification of molecular species due to intrinsic limitations of the technique (lack of chromatographic separation and the difficulty to apply tandem MS). Successful strategies to perform annotation and identification combine extra analytical steps, like using orthogonal analytical techniques to identify compounds; with algorithms that integrate the spectral and spatial information. In this review, we discuss different experimental strategies and bioinformatics tools to annotate and identify compounds in MSI experiments. We target strategies and tools for small molecule applications, such as lipidomics and metabolomics. First, we explain how sample preparation and the acquisition process influences annotation and identification, from sample preservation to the use of orthogonal techniques. Then, we review twelve software tools for annotation and identification in MSI. Finally, we offer perspectives on two current needs of the MSI community: the adaptation of guidelines for communicating confidence levels in identifications; and the creation of a standard format to store and exchange annotations and identifications in MSI.

Quantifying Spatial Heterogeneity of Tumor-Infiltrating Lymphocytes to Predict Survival of Individual Cancer Patients

Tumor-infiltrating lymphocytes (TILs), identified on HE-stained histopathological images in the cancer area, are indicators of the adaptive immune response against cancers and play a major role in personalized cancer immunotherapy. Recent works indicate that the spatial organization of TILs may be prognostic of disease-specific survival and recurrence. However, there are a limited number of methods that were proposed and tested in analyses of the spatial structure of TILs. In this work, we evaluated 14 different spatial measures, including the one developed for other omics data, on 10,532 TIL maps from 23 cancer types in terms of reproducibility, uniqueness, and impact on patient survival. For each spatial measure, 16 different scenarios for the definition of prognostic factor were tested. We found no difference in survival prediction when TIL maps were stored as binary images or continuous TIL probability scores. When spatial measures were discretized into a low and high category, a higher correlation with survival was observed. Three measures with the highest cancer prognosis capability were spatial autocorrelation, GLCM M1, and closeness centrality. Most of the tested measures could be further tuned to increase prediction performance.

MALDI-MSI Towards Multimodal Imaging: Challenges and Perspectives

Multimodal imaging is a powerful strategy for combining information from multiple images. It involves several fields in the acquisition, processing and interpretation of images. As multimodal imaging is a vast subject area with various combinations of imaging techniques, it has been extensively reviewed. Here we focus on Matrix-assisted Laser Desorption Ionization Mass Spectrometry Imaging (MALDI-MSI) coupling other imaging modalities in multimodal approaches. While MALDI-MS images convey a substantial amount of chemical information, they are not readily informative about the morphological nature of the tissue. By providing a supplementary modality, MALDI-MS images can be more informative and better reflect the nature of the tissue. In this mini review, we emphasize the analytical and computational strategies to address multimodal MALDI-MSI.

A versatile deep learning architecture for classification and label-free prediction of hyperspectral images

Hyperspectral imaging is a technique that provides rich chemical or compositional information not regularly available to traditional imaging modalities such as intensity imaging or color imaging based on the reflection, transmission, or emission of light. Analysis of hyperspectral imaging often relies on machine learning methods to extract information. Here, we present a new flexible architecture, the U-within-U-Net, that can perform classification, segmentation, and prediction of orthogonal imaging modalities on a variety of hyperspectral imaging techniques. Specifically, we demonstrate feature segmentation and classification on the Indian Pines hyperspectral dataset and simultaneous location prediction of multiple drugs in mass spectrometry imaging of rat liver tissue. We further demonstrate label-free fluorescence image prediction from hyperspectral stimulated Raman scattering microscopy images. The applicability of the U-within-U-Net architecture on diverse datasets with widely varying input and output dimensions and data sources suggest that it has great potential in advancing the use of hyperspectral imaging across many different application areas ranging from remote sensing, to medical imaging, to microscopy.

Data Filtering and Its Prioritization in Pipelines for Spatial Segmentation of Mass Spectrometry Imaging

Mass spectrometry imaging (MSI) could provide vast amounts of data at the temporal-spatial scale in heterogeneous biological specimens, which challenges us to segment accurately suborgans/microregions from complex MSI data. Several pipelines had been proposed for MSI spatial segmentation in the past decade. More importantly, data filtering was found to be an efficient procedure to improve the outcomes of MSI segmentation pipelines. It is not clear, however, how the filtering procedure affects the MSI segmentation. An improved pipeline was established by elaborating the filtering prioritization and filtering algorithm. Lipidomic-characteristic-based MSI data of a whole-body mouse fetus was used to evaluate the established pipeline on localization of the physiological position of suborgans by comparing with three commonly used pipelines and commercial SCiLS Lab software. Two structural measurements were used to quantify the performances of the pipelines including the percentage of abnormal edge pixel (PAEP) and CHAOS. Our results demonstrated that the established pipeline outperformed the other pipelines in visual inspection, spatial consistence, time-cost, and robustness analysis. For example, the dorsal pallium (isocortex) and hippocampal formation (Hpf) regions, midbrain, cerebellum, and brainstem on the mouse brain were annotated and located by the established pipeline. As a generic pipeline, the established pipeline could help with the accurate assessment and screening of drug/chemical-induced targeted organs and exploration of the progression and molecular mechanisms of diseases. The filter-based strategy is expected to become a critical component in the standard operating procedure of MSI data sets.

Esmraldi: efficient methods for the fusion of mass spectrometry and magnetic resonance images

Background

Mass spectrometry imaging (MSI) is a family of acquisition techniques producing images of the distribution of molecules in a sample, without any prior tagging of the molecules. This makes it a very interesting technique for exploratory research. However, the images are difficult to analyze because the enclosed data has high dimensionality, and their content does not necessarily reflect the shape of the object of interest. Conversely, magnetic resonance imaging (MRI) scans reflect the anatomy of the tissue. MRI also provides complementary information to MSI, such as the content and distribution of water.

Results

We propose a new workflow to merge the information from 2D MALDI–MSI and MRI images. Our workflow can be applied to large MSI datasets in a limited amount of time. Moreover, the workflow is fully automated and based on deterministic methods which ensures the reproducibility of the results. Our methods were evaluated and compared with state-of-the-art methods. Results show that the images are combined precisely and in a time-efficient manner.

Conclusion

Our workflow reveals molecules which co-localize with water in biological images. It can be applied on any MSI and MRI datasets which satisfy a few conditions: same regions of the shape enclosed in the images and similar intensity distributions.

Implications of Peak Selection in the Interpretation of Unsupervised Mass Spectrometry Imaging Data Analyses

Mass spectrometry imaging can produce large amounts of complex spectral and spatial data. Such data sets are often analyzed with unsupervised machine learning approaches, which aim at reducing their complexity and facilitating their interpretation. However, choices made during data processing can impact the overall interpretation of these analyses. This work investigates the impact of the choices made at the peak selection step, which often occurs early in the data processing pipeline. The discussion is done in terms of visualization and interpretation of the results of two commonly used unsupervised approaches: t-distributed stochastic neighbor embedding and k-means clustering, which differ in nature and complexity. Criteria considered for peak selection include those based on hypotheses (exemplified herein in the analysis of metabolic alterations in genetically engineered mouse models of human colorectal cancer), particular molecular classes, and ion intensity. The results suggest that the choices made at the peak selection step have a significant impact in the visual interpretation of the results of either dimensionality reduction or clustering techniques and consequently in any downstream analysis that relies on these. Of particular significance, the results of this work show that while using the most abundant ions can result in interesting structure-related segmentation patterns that correlate well with histological features, using a smaller number of ions specifically selected based on prior knowledge about the biochemistry of the tissues under investigation can result in an easier-to-interpret, potentially more valuable, hypothesis-confirming result. Findings presented will help researchers understand and better utilize unsupervised machine learning approaches to mine high-dimensionality data.

Spatial Metabolomics and Imaging Mass Spectrometry in the Age of Artificial Intelligence

Spatial metabolomics is an emerging field of omics research that has enabled localizing metabolites, lipids, and drugs in tissue sections, a feat considered impossible just two decades ago. Spatial metabolomics and its enabling technology-imaging mass spectrometry-generate big hyper-spectral imaging data that have motivated the development of tailored computational methods at the intersection of computational metabolomics and image analysis. Experimental and computational developments have recently opened doors to applications of spatial metabolomics in life sciences and biomedicine. At the same time, these advances have coincided with a rapid evolution in machine learning, deep learning, and artificial intelligence, which are transforming our everyday life and promise to revolutionize biology and healthcare. Here, we introduce spatial metabolomics through the eyes of a computational scientist, review the outstanding challenges, provide a look into the future, and discuss opportunities granted by the ongoing convergence of human and artificial intelligence.

Spatial metabolomics of in situ host-microbe interactions at the micrometre scale

Spatial metabolomics describes the location and chemistry of small molecules involved in metabolic phenotypes, defence molecules and chemical interactions in natural communities. Most current techniques are unable to spatially link the genotype and metabolic phenotype of microorganisms in situ at a scale relevant to microbial interactions. Here, we present a spatial metabolomics pipeline (metaFISH) that combines fluorescence in situ hybridization (FISH) microscopy and high-resolution atmospheric-pressure matrix-assisted laser desorption/ionization mass spectrometry to image host-microbe symbioses and their metabolic interactions. The metaFISH pipeline aligns and integrates metabolite and fluorescent images at the micrometre scale to provide a spatial assignment of host and symbiont metabolites on the same tissue section. To illustrate the advantages of metaFISH, we mapped the spatial metabolome of a deep-sea mussel and its intracellular symbiotic bacteria at the scale of individual epithelial host cells. Our analytical pipeline revealed metabolic adaptations of the epithelial cells to the intracellular symbionts and variation in metabolic phenotypes within a single symbiont 16S rRNA phylotype, and enabled the discovery of specialized metabolites from the host-microbe interface. metaFISH provides a culture-independent approach to link metabolic phenotypes to community members in situ and is a powerful tool for microbiologists across fields.

LipostarMSI: Comprehensive, Vendor-Neutral Software for Visualization, Data Analysis, and Automated Molecular Identification in Mass Spectrometry Imaging

Mass Spectrometry Imaging (MSI) is an established and powerful MS technique that enables molecular mapping of tissues and cells finding widespread applications in academic, medical, and pharmaceutical industries. As both the applications and MSI technology have undergone rapid growth and improvement, the challenges associated both with analyzing large datasets and identifying the many detected molecular species have become apparent. The lack of readily available and comprehensive software covering all necessary data analysis steps has further compounded this challenge. To address this issue we developed LipostarMSI, comprehensive and vendor-neutral software for targeted and untargeted MSI data analysis. Through user-friendly implementation of image visualization and co-registration, univariate and multivariate image and spectral analysis, and for the first time, advanced lipid, metabolite, and drug metabolite (MetID) automated identification, LipostarMSI effectively streamlines biochemical interpretation of the data. Here, we introduce LipostarMSI and case studies demonstrating the versatility and many capabilities of the software.

SPUTNIK: an R package for filtering of spatially related peaks in mass spectrometry imaging data

Summary

SPUTNIK is an R package consisting of a series of tools to filter mass spectrometry imaging peaks characterized by a noisy or unlikely spatial distribution. SPUTNIK can produce mass spectrometry imaging datasets characterized by a smaller but more informative set of peaks, reduce the complexity of subsequent multi-variate analysis and increase the interpretability of the statistical results.

Availability and implementation

SPUTNIK is freely available online from CRAN repository and at https://github.com/paoloinglese/SPUTNIK. The package is distributed under the GNU General Public License version 3 and is accompanied by example files and data.

Supplementary information

Supplementary data are available at Bioinformatics online.

Clusterwise Peak Detection and Filtering Based on Spatial Distribution To Efficiently Mine Mass Spectrometry Imaging Data

Mass spectrometry imaging (MSI) has the potential to reveal the localization of thousands of biomolecules such as metabolites and lipids in tissue sections. The increase in both mass and spatial resolution of today’s instruments brings on considerable challenges in terms of data processing; accurately extracting meaningful signals from the large data sets generated by MSI without losing information that could be clinically relevant is one of the most fundamental tasks of analysis software. Ion images of the biomolecules are generated by visualizing their intensities in 2-D space using mass spectra collected across the tissue section. The intensities are often calculated by summing each compound’s signal between predefined sets of borders (bins) in the m/z dimension. This approach, however, can result in mixed signals from different compounds in the same bin or splitting the signal from one compound between two adjacent bins, leading to low quality ion images. To remedy this problem, we propose a novel data processing approach. Our approach consists of a sensitive peak detection method able to discover both faint and localized signals by utilizing clusterwise kernel density estimates (KDEs) of peak distributions. We show that our method can recall more ground-truth molecules, molecule fragments, and isotopes than existing methods based on binning. Furthermore, it automatically detects previously reported molecular ions of lipids, including those close in m/z, in an experimental data set.

Data-Driven Rescoring of Metabolite Annotations Significantly Improves Sensitivity

When analyzing mass spectrometry imaging data sets, assigning a molecule to each of the thousands of generated images is a very complex task. Recent efforts have taken lessons from (tandem) mass spectrometry proteomics and applied them to imaging mass spectrometry metabolomics, with good results. Our goal is to go a step further in this direction and apply a well established, data-driven method to improve the results obtained from an annotation engine. By using a data-driven rescoring strategy, we are able to consistently improve the sensitivity of the annotation engine while maintaining control of statistics like estimated rate of false discoveries. All the code necessary to run a search and extract the additional features can be found at https://github.com/anasilviacs/sm-engine and to rescore the results from a search in https://github.com/anasilviacs/rescore-metabolites .

Characterization of microbial mixtures by mass spectrometry

MS applications in microbiology have increased significantly in the past 10 years, due in part to the proliferation of regulator-approved commercial MALDI MS platforms for rapid identification of clinical infections. In parallel, with the expansion of MS technologies in the "omics" fields, novel MS-based research efforts to characterize organismal as well as environmental microbiomes have emerged. Successful characterization of microorganisms found in complex mixtures of other organisms remains a major challenge for researchers and clinicians alike. Here, we review recent MS advances toward addressing that challenge. These include sample preparation methods and protocols, and established, for example, MALDI, as well as newer, for example, atmospheric pressure ionization (API) techniques. MALDI mass spectra of intact cells contain predominantly information on the highly expressed house-keeping proteins used as biomarkers. The API methods are applicable for small biomolecule analysis, for example, phospholipids and lipopeptides, and facilitate species differentiation. MS hardware and techniques, for example, tandem MS, including diverse ion source/mass analyzer combinations are discussed. Relevant examples for microbial mixture characterization utilizing these combinations are provided. Chemometrics and bioinformatics methods and algorithms, including those applied to large scale MS data acquisition in microbial metaproteomics and MS imaging of biofilms, are highlighted. Select MS applications for polymicrobial culture analysis in environmental and clinical microbiology are reviewed as well.

BASIS: High-performance bioinformatics platform for processing of large-scale mass spectrometry imaging data in chemically augmented histology

Mass Spectrometry Imaging (MSI) holds significant promise in augmenting digital histopathologic analysis by generating highly robust big data about the metabolic, lipidomic and proteomic molecular content of the samples. In the process, a vast quantity of unrefined data, that can amount to several hundred gigabytes per tissue section, is produced. Managing, analysing and interpreting this data is a significant challenge and represents a major barrier to the translational application of MSI. Existing data analysis solutions for MSI rely on a set of heterogeneous bioinformatics packages that are not scalable for the reproducible processing of large-scale (hundreds to thousands) biological sample sets. Here, we present a computational platform (pyBASIS) capable of optimized and scalable processing of MSI data for improved information recovery and comparative analysis across tissue specimens using machine learning and related pattern recognition approaches. The proposed solution also provides a means of seamlessly integrating experimental laboratory data with downstream bioinformatics interpretation/analyses, resulting in a truly integrated system for translational MSI.

Automated Morphological and Morphometric Analysis of Mass Spectrometry Imaging Data: Application to Biomarker Discovery

Mass spectrometry imaging datasets are mostly analyzed in terms of average intensity in regions of interest. However, biological tissues have different morphologies with several sizes, shapes, and structures. The important biological information, contained in this highly heterogeneous cellular organization, could be hidden by analyzing the average intensities. Finding an analytical process of morphology would help to find such information, describe tissue model, and support identification of biomarkers. This study describes an informatics approach for the extraction and identification of mass spectrometry image features and its application to sample analysis and modeling. For the proof of concept, two different tissue types (healthy kidney and CT-26 xenograft tumor tissues) were imaged and analyzed. A mouse kidney model and tumor model were generated using morphometric - number of objects and total surface - information. The morphometric information was used to identify m/z that have a heterogeneous distribution. It seems to be a worthwhile pursuit as clonal heterogeneity in a tumor is of clinical relevance. This study provides a new approach to find biomarker or support tissue classification with more information. Graphical Abstract ᅟ.

FDR-controlled metabolite annotation for high-resolution imaging mass spectrometry

High-mass-resolution imaging mass spectrometry promises to localize hundreds of metabolites in tissues, cell cultures, and agar plates with cellular resolution, but it is hampered by the lack of bioinformatics tools for automated metabolite identification. We report pySM, a framework for false discovery rate (FDR)-controlled metabolite annotation at the level of the molecular sum formula, for high-mass-resolution imaging mass spectrometry (https://github.com/alexandrovteam/pySM). We introduce a metabolite-signal match score and a target-decoy FDR estimate for spatial metabolomics.

Probabilistic Segmentation of Mass Spectrometry (MS) Images Helps Select Important Ions and Characterize Confidence in the Resulting Segments

Mass spectrometry imaging is a powerful tool for investigating the spatial distribution of chemical compounds in a biological sample such as tissue. Two common goals of these experiments are unsupervised segmentation of images into newly discovered homogeneous segments and supervised classification of images into predefined classes. In both cases, the important secondary goals are to characterize the uncertainty associated with the segmentation and with the classification and to characterize the spectral features that define each segment or class. Recent analysis methods have focused on the spatial structure of the data to improve results. However, they either do not address these secondary goals or do this with separate post hoc procedures.

We introduce spatial shrunken centroids, a statistical model-based framework for both supervised classification and unsupervised segmentation. It takes as input sets of previously detected, aligned, quantified, and normalized spectral features and expresses both spatial and multivariate nature of the data using probabilistic modeling. It selects informative subsets of spectral features that define each unsupervised segment or supervised class and quantifies and visualizes the uncertainty in spatial segmentations and in tissue classification. In the unsupervised setting, it also guides the choice of an appropriate number of segments. We demonstrate the usefulness of this framework in a supervised human renal cell carcinoma experimental dataset and several unsupervised experimental datasets, including a pig fetus cross-section, three rodent brains, and a controlled image with known ground truth. This framework is available for use within the open-source R package Cardinal as part of a full pipeline for the processing, visualization, and statistical analysis of mass spectrometry imaging experiments.

Using collective expert judgements to evaluate quality measures of mass spectrometry images

Motivation: Imaging mass spectrometry (IMS) is a maturating technique of molecular imaging. Confidence in the reproducible quality of IMS data is essential for its integration into routine use. However, the predominant method for assessing quality is visual examination, a time consuming, unstandardized and non-scalable approach. So far, the problem of assessing the quality has only been marginally addressed and existing measures do not account for the spatial information of IMS data. Importantly, no approach exists for unbiased evaluation of potential quality measures.

Results: We propose a novel approach for evaluating potential measures by creating a gold-standard set using collective expert judgements upon which we evaluated image-based measures. To produce a gold standard, we engaged 80 IMS experts, each to rate the relative quality between 52 pairs of ion images from MALDI-TOF IMS datasets of rat brain coronal sections. Experts’ optional feedback on their expertise, the task and the survey showed that (i) they had diverse backgrounds and sufficient expertise, (ii) the task was properly understood, and (iii) the survey was comprehensible. A moderate inter-rater agreement was achieved with Krippendorff’s alpha of 0.5. A gold-standard set of 634 pairs of images with accompanying ratings was constructed and showed a high agreement of 0.85. Eight families of potential measures with a range of parameters and statistical descriptors, giving 143 in total, were evaluated. Both signal-to-noise and spatial chaos-based measures performed highly with a correlation of 0.7 to 0.9 with the gold standard ratings. Moreover, we showed that a composite measure with the linear coefficients (trained on the gold standard with regularized least squares optimization and lasso) showed a strong linear correlation of 0.94 and an accuracy of 0.98 in predicting which image in a pair was of higher quality.

Availability and implementation: The anonymized data collected from the survey and the Matlab source code for data processing can be found at: https://github.com/alexandrovteam/IMS_quality.

Contact:theodore.alexandrov@embl.de

An efficient TOF-SIMS image analysis with spatial correlation and alternating non-negativity-constrained least squares

MOTIVATION

Advances in analytical instrumentation towards acquiring high-resolution images of mass spectrometry constantly demand efficient approaches for data analysis. This is particularly true of time-of-flight secondary ion mass spectrometry imaging where recent advances enable acquisition of high-resolution data in multiple dimensions. In many applications, the distribution of different species from a sampled surface is spatially continuous in nature and a model that incorporates the spatial correlation across the surface would be preferable to estimations at discrete spatial locations. A key challenge here is the capability to analyse the high-resolution multidimensional data to extract relevant information reliably and efficiently.

RESULTS

We propose a framework based on alternating non-negativity-constrained least squares which accounts for the spatial correlation across the sample surface. The proposed method also decouples the computational complexity of the estimation procedure from the image resolution, which significantly reduces the processing time. We evaluate the performance of the algorithm with biochemical image datasets generated from mixture of metabolites.