A mathematical comparison of non-negative matrix factorization related methods with practical implications for the analysis of mass spectrometry imaging data

Adipose tissue composition determines its computed tomography radiodensity

OBJECTIVES

Adipose tissue radiodensity in computed tomography (CT) performed before surgeries can predict surgical difficulty. Despite its clinical importance, little is known about what influences radiodensity. This study combines desorption electrospray ionization mass spectrometry imaging (DESI-MSI) and electrospray ionization (ESI) with machine learning to unveil how chemical composition of adipose tissue determines its radiodensity.

METHODS

Patients in the study underwent abdominal surgeries. Before surgery, CT radiodensity of fat near operated sites was measured. Fifty-three fat samples were collected and analyzed by DESI-MSI, ESI, and histology, and then sorted by radiodensity, demographic parameters, and adipocyte size. A non-negative matrix factorization (NMF) algorithm was developed to differentiate between high and low radiodensities.

RESULTS

No associations between radiodensity and patient age, gender, weight, height, or fat origin were found. Body mass index showed negative correlation with radiodensity. A substantial difference in chemical composition between adipose tissues of high and low radiodensities was observed. More radiodense tissues exhibited greater abundance of high molecular weight species, such as phospholipids of various types, ceramides, cholesterol esters and diglycerides, and about 70% smaller adipocyte size. Less radiodense tissue showed high abundance of short acyl-tail fatty acids.

CONCLUSIONS

This study unveils the connection between abdominal adipose tissue radiodensity and its chemical composition. Because the radiodensity of the fat around the surgical site is associated with surgical difficulty, it is important to understand how adipose tissue composition affects this parameter. We conclude that fat tissue with a higher content of various phospholipids and waxy lipids is more CT radiodense.

CLINICAL RELEVANCE STATEMENT

This study establishes the connection between the CT radiodensity of adipose tissue and its chemical composition. Clinicians may use this information for preoperative planning of surgical procedures, potentially modifying their surgical approach (for example, performing partial nephrectomy openly rather than laparoscopically).

KEY POINTS

• Adipose tissue radiodensity values in computed tomography images taken prior to the surgery can potentially predict surgery difficulty. • Fifty-three human specimens were analyzed by advanced mass spectrometry, molecular imaging, and machine learning to establish the key features that determine Hounsfield units' values of adipose tissue. • The findings of this research will enable clinicians to better prepare for surgical procedures and select operative strategies.

Spatiochemical Characterization of the Pancreas Using Mass Spectrometry Imaging and Topological Data Analysis

Mass Spectrometry Imaging (MSI) is a technique used to identify the spatial distribution of molecules in tissues. An MSI experiment results in large amounts of high dimensional data, so efficient computational methods are needed to analyze the output. Topological Data Analysis (TDA) has proven to be effective in all kinds of applications. TDA focuses on the topology of the data in high dimensional space. Looking at the shape in a high dimensional data set can lead to new or different insights. In this work, we investigate the use of Mapper, a form of TDA, applied on MSI data. Mapper is used to find data clusters inside two healthy mouse pancreas data sets. The results are compared to previous work using UMAP for MSI data analysis on the same data sets. This work finds that the proposed technique discovers the same clusters in the data as UMAP and is also able to uncover new clusters, such as an additional ring structure inside the pancreatic islets and a better defined cluster containing blood vessels. The technique can be used for a large variety of data types and sizes and can be optimized for specific applications. It is also computationally similar to UMAP for clustering. Mapper is a very interesting method, especially its use in biomedical applications.

Automated Library Generation and Serendipity Quantification Enables Diverse Discovery in Coordination Chemistry

Library generation experiments are a key part of the discovery of new materials, methods, and models in chemistry, but the question of how to generate high quality libraries to enable discovery is nontrivial. Herein, we use coordination chemistry to demonstrate the automation of many of the workflows used for library generation in automated hardware including the Chemputer. First, we explore the target-oriented synthesis of three influential coordination complexes, to validate key synthetic operations in our system; second, the generation of focused libraries in chemical and process space; and third, the development of a new workflow for prospecting library formation. This involved Bayesian optimization using a Gaussian process as surrogate model combined with a metric for novelty (or serendipity) quantification based on mass spectrometry data. In this way, we show directed exploration of a process space toward those areas with rarer observations and build a picture of the diversity in product distributions present across the space. We show that this effectively "engineers" serendipity into our search through the unexpected appearance of acetic anhydride, formed in situ, and solvent degradation products as ligands in an isolable series of three Co(III) anhydride complexes.

Emerging Computational Methods in Mass Spectrometry Imaging

Mass spectrometry imaging (MSI) is a powerful analytical technique that generates maps of hundreds of molecules in biological samples with high sensitivity and molecular specificity. Advanced MSI platforms with capability of high-spatial resolution and high-throughput acquisition generate vast amount of data, which necessitates the development of computational tools for MSI data analysis. In addition, computation-driven MSI experiments have recently emerged as enabling technologies for further improving the MSI capabilities with little or no hardware modification. This review provides a critical summary of computational methods and resources developed for MSI data analysis and interpretation along with computational approaches for improving throughput and molecular coverage in MSI experiments. This review is focused on the recently developed artificial intelligence methods and provides an outlook for a future paradigm shift in MSI with transformative computational methods.

Applications of multivariate analysis and unsupervised machine learning to ToF-SIMS images of organic, bioorganic, and biological systems

Time-of-flight secondary ion mass spectrometry (ToF-SIMS) imaging offers a powerful, label-free method for exploring organic, bioorganic, and biological systems. The technique is capable of very high spatial resolution, while also producing an enormous amount of information about the chemical and molecular composition of a surface. However, this information is inherently complex, making interpretation and analysis of the vast amount of data produced by a single ToF-SIMS experiment a considerable challenge. Much research over the past few decades has focused on the application and development of multivariate analysis (MVA) and machine learning (ML) techniques that find meaningful patterns and relationships in these datasets. Here, we review the unsupervised algorithms—that is, algorithms that do not require ground truth labels—that have been applied to ToF-SIMS images, as well as other algorithms and approaches that have been used in the broader family of mass spectrometry imaging (MSI) techniques. We first give a nontechnical overview of several commonly used classes of unsupervised algorithms, such as matrix factorization, clustering, and nonlinear dimensionality reduction. We then review the application of unsupervised algorithms to various organic, bioorganic, and biological systems including cells and tissues, organic films, residues and coatings, and spatially structured systems such as polymer microarrays. We then cover several novel algorithms employed for other MSI techniques that have received little attention from ToF-SIMS imaging researchers. We conclude with a brief outline of potential future directions for the application of MVA and ML algorithms to ToF-SIMS images.