Title: Interpretable dimensionality reduction and classification of mass spectrometry imaging data in a visceral pain model via non-negative matrix factorization
Mass spectrometry imaging (MSI) is a powerful scientific tool for understanding the spatial distribution of biochemical compounds in tissue structures. In this paper, we introduce three novel approaches in MSI data processing to perform the tasks of data augmentation, feature ranking, and image registration. We use these approaches in conjunction with non-negative matrix factorization (NMF) to resolve two of the biggest challenges in MSI data analysis, namely: 1) the large file sizes and associated computational resource requirements and 2) the complexity of interpreting the very high dimensional raw spectral data. There are many dimensionality reduction techniques that address the first challenge but do not necessarily result in readily interpretable features, leaving the second challenge unaddressed. We demonstrate that NMF is an effective dimensionality reduction algorithm that reduces the size of MSI datasets by three orders of magnitude with limited loss of information, yielding spatial and spectral components with meaningful correlation to tissue structure that may be used directly for subsequent data analysis without the need for additional clustering steps. This analysis is demonstrated on an MSI dataset from female Sprague-Dawley rats for an animal model of comorbid visceral pain hypersensitivity (CPH). We find that high-dimensional MSI data (∼ 100,000 ions per pixel) can be reduced to 20 spectral NMF components with < 20% loss in reconstruction accuracy. The resulting spatial NMF components are reproducible and correlate well with H&E-stained tissue images. These components may also be used to generate images with enhanced specificity for different tissue types. Small patches of NMF data (i.e., 20 spatial NMF components over 20 × 20 pixels) provide an accuracy of ∼ 87% in classifying CPH vs naïve control subjects. This paper presents the novel data processing methodologies that were used to produce these results, encompassing novel data processing pipelines for data augmentation to support training for classification, ranking of features according to their contribution to classification, and image registration to enhance tissue-specific imaging.  more » « less
Award ID(s):
2020624
PAR ID:
10565622
Author(s) / Creator(s):
; ; ; ; ; ; ;
Editor(s):
Yousefi, Bardia
Publisher / Repository:
PLOS
Date Published:
Journal Name:
PLOS ONE
Volume:
19
Issue:
10
ISSN:
1932-6203
Page Range / eLocation ID:
e0300526
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Laser ablation inductively-coupled plasma mass spectrometry (LA-ICP-MS) imaging and matrix assisted laser desorption ionization mass spectrometry imaging (MALDI-MSI) are complementary methods that measure distributions of elements and biomolecules in tissue sections. Quantitative correlations of the information provided by these two imaging modalities requires that the datasets be registered in the same coordinate system, allowing for pixel-by-pixel comparisons. We describe here a computational workflow written in Python that accomplishes this registration, even for adjacent tissue sections, with accuracies within ±50 μm. The value of this registration process is demonstrated by correlating images of tissue sections from mice injected with gold nanomaterial drug delivery systems. Quantitative correlations of the nanomaterial delivery vehicle, as detected by LA-ICP-MS imaging, with biochemical changes, as detected by MALDI-MSI, provide deeper insight into how nanomaterial delivery systems influence lipid biochemistry in tissues. Moreover, the registration process allows the more precise images associated with LA-ICP-MS imaging to be leveraged to achieve improved segmentation in MALDI-MS images, resulting in the identification of lipids that are most associated with different sub-organ regions in tissues. 
    more » « less
  2. Spatial transcriptomics enables spatially resolved gene expression analysis, but accompanying histology images are often degraded by fiducial markers and background regions, hindering interpretation. To address this, we introduce Vispro, an end-to-end automated image processing tool optimized for 10× Visium data. Vispro includes modules for fiducial marker detection, image restoration, tissue region detection, and segmentation of disconnected tissue areas. By enhancing image quality, Vispro improves the accuracy and performance of downstream analyses, including tissue and cell segmentation, image registration, gene expression imputation guided by histological context, and spatial domain detection. 
    more » « less
  3. Motivation: Mass spectrometry imaging (MSI) characterizes the molecular composition of tissues at spatial resolution, and has a strong potential for distinguishing tissue types, or disease states. This can be achieved by supervised classification, which takes as input MSI spectra, and assigns class labels to subtissue locations. Unfortunately, developing such classifiers is hindered by the limited availability of training sets with subtissue labels as the ground truth. Subtissue labeling is prohibitively expensive, and only rough annotations of the entire tissues are typically available. Classifiers trained on data with approximate labels have sub-optimal performance. Results: To alleviate this challenge, we contribute a semi-supervised approach mi-CNN. mi-CNN implements multiple instance learning with a convolutional neural network (CNN). The multiple instance aspect enables weak supervision from tissue-level annotations when classifying subtissue locations. The convolutional architecture of the CNN captures contextual dependencies between the spectral features. Evaluations on simulated and experimental datasets demonstrated that mi-CNN improved the subtissue classification as compared to traditional classifiers. We propose mi-CNN as an important step towards accurate subtissue classification in MSI, enabling rapid distinction between tissue types and disease states. 
    more » « less
  4. Abstract We introduce a new method based on nonnegative matrix factorization, Neural NMF, for detecting latent hierarchical structure in data. Datasets with hierarchical structure arise in a wide variety of fields, such as document classification, image processing, and bioinformatics. Neural NMF recursively applies NMF in layers to discover overarching topics encompassing the lower-level features. We derive a backpropagation optimization scheme that allows us to frame hierarchical NMF as a neural network. We test Neural NMF on a synthetic hierarchical dataset, the 20 Newsgroups dataset, and the MyLymeData symptoms dataset. Numerical results demonstrate that Neural NMF outperforms other hierarchical NMF methods on these data sets and offers better learned hierarchical structure and interpretability of topics. 
    more » « less
  5. Tugizimana, Fidele (Ed.)
    Research and innovation in metabolomics tools to measure metabolite accumulation within plants have led to important discoveries with respect to the improvement of plant stress tolerance, development, and crop yield. Traditional metabolomics analyses have commonly utilized gas chromatography–mass spectrometry and liquid chromatography–mass spectrometry, but these methods are often performed without regard for the spatial locations of metabolites within tissues. Methods for mass spectral imaging (MSI) have recently been developed to detect and spatially resolve metabolite accumulation and are rapidly being adopted on a wider scale. Since 2010, the number of publications incorporating mass spectral imaging has grown from approximately 80 articles to over 378 on a yearly basis, constituting an increase of at least 350% during this time frame. Spatially resolved metabolite accumulation data provides unique insights into the function and regulation of plant biochemical pathways. Mass spectral imaging is commonly paired with desorption ionization technologies, including matrix-assisted laser desorption ionization (MALDI) and desorption electrospray ionization (DESI), to generate accurate, spatially resolved metabolomics data from prepared tissue segments. Here, we describe the most recent advancements in sample preparation methods, mass spectral imaging technologies, and data processing tools that have been developed to address the limits of MSI technology. Additionally, we summarize recent applications of MSI technologies in plant metabolomics and discuss potential avenues for future research advancements within the plant biology community through the use of these technologies. 
    more » « less