Quantifying biological samples using Linear Poisson Independent Component Analysis for MALDI-ToF mass spectra

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.

Mass Spectrometry Imaging for Spatial Chemical Profiling of Vegetative Parts of Plants

The detection of chemical species and understanding their respective localisations in tissues have important implications in plant science. The conventional methods for imaging spatial localisation of chemical species are often restricted by the number of species that can be identified and is mostly done in a targeted manner. Mass spectrometry imaging combines the ability of traditional mass spectrometry to detect numerous chemical species in a sample with their spatial localisation information by analysing the specimen in a 2D manner. This article details the popular mass spectrometry imaging methodologies which are widely pursued along with their respective sample preparation and the data analysis methods that are commonly used. We also review the advancements through the years in the usage of the technique for the spatial profiling of endogenous metabolites, detection of xenobiotic agrochemicals and disease detection in plants. As an actively pursued area of research, we also address the hurdles in the analysis of plant tissues, the future scopes and an integrated approach to analyse samples combining different mass spectrometry imaging methods to obtain the most information from a sample of interest.

Habitat Imaging of Tumors Enables High Confidence Sub-Regional Assessment of Response to Therapy

Imaging biomarkers are used in therapy development to identify and quantify therapeutic response. In oncology, use of MRI, PET and other imaging methods can be complicated by spatially complex and heterogeneous tumor micro-environments, non-Gaussian data and small sample sizes. Linear Poisson Modelling (LPM) enables analysis of complex data that is quantitative and can operate in small data domains. We performed experiments in 5 mouse models to evaluate the ability of LPM to identify responding tumor habitats across a range of radiation and targeted drug therapies. We tested if LPM could identify differential biological response rates. We calculated the theoretical sample size constraints for applying LPM to new data. We then performed a co-clinical trial using small data to test if LPM could detect multiple therapeutics with both improved power and reduced animal numbers compared to conventional t-test approaches. Our data showed that LPM greatly increased the amount of information extracted from diffusion-weighted imaging, compared to cohort t-tests. LPM distinguished biological response rates between Calu6 tumors treated with 3 different therapies and between Calu6 tumors and 4 other xenograft models treated with radiotherapy. A simulated co-clinical trial using real data detected high precision per-tumor treatment effects in as few as 3 mice per cohort, with p-values as low as 1 in 10,000. These findings provide a route to simultaneously improve the information derived from preclinical imaging while reducing and refining the use of animals in cancer research.

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

RATIONALE

Non-negative matrix factorization (NMF) has been used extensively for the analysis of mass spectrometry imaging (MSI) data, visualizing simultaneously the spatial and spectral distributions present in a slice of tissue. The statistical framework offers two related NMF methods: probabilistic latent semantic analysis (PLSA) and latent Dirichlet allocation (LDA), which is a generative model. This work offers a mathematical comparison between NMF, PLSA, and LDA, and includes a detailed evaluation of Kullback-Leibler NMF (KL-NMF) for MSI for the first time. We will inspect the results for MSI data analysis as these different mathematical approaches impose different characteristics on the data and the resulting decomposition.

METHODS

The four methods (NMF, KL-NMF, PLSA, and LDA) are compared on seven different samples: three originated from mice pancreas and four from human-lymph-node tissues, all obtained using matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS).

RESULTS

Where matrix factorization methods are often used for the analysis of MSI data, we find that each method has different implications on the exactness and interpretability of the results. We have discovered promising results using KL-NMF, which has only rarely been used for MSI so far, improving both NMF and PLSA, and have shown that the hitherto stated equivalent KL-NMF and PLSA algorithms do differ in the case of MSI data analysis. LDA, assumed to be the better method in the field of text mining, is shown to be outperformed by PLSA in the setting of MALDI-MSI. Additionally, the molecular results of the human-lymph-node data have been thoroughly analyzed for better assessment of the methods under investigation.

CONCLUSIONS

We present an in-depth comparison of multiple NMF-related factorization methods for MSI. We aim to provide fellow researchers in the field of MSI a clear understanding of the mathematical implications using each of these analytical techniques, which might affect the exactness and interpretation of the results.

Quantitative Mass Spectrometry Imaging of Biological Systems

Mass spectrometry imaging (MSI) is a powerful, label-free technique that provides detailed maps of hundreds of molecules in complex samples with high sensitivity and subcellular spatial resolution. Accurate quantification in MSI relies on a detailed understanding of matrix effects associated with the ionization process along with evaluation of the extraction efficiency and mass-dependent ion losses occurring in the analysis step. We present a critical summary of approaches developed for quantitative MSI of metabolites, lipids, and proteins in biological tissues and discuss their current and future applications.

A reformulation of pLSA for uncertainty estimation and hypothesis testing in bio-imaging

MOTIVATION

Probabilistic latent semantic analysis (pLSA) is commonly applied to describe mass spectra (MS) images. However, the method does not provide certain outputs necessary for the quantitative scientific interpretation of data. In particular, it lacks assessment of statistical uncertainty and the ability to perform hypothesis testing. We show how linear Poisson modelling advances pLSA, giving covariances on model parameters and supporting χ2 testing for the presence/absence of MS signal components. As an example, this is useful for the identification of pathology in MALDI biological samples. We also show potential wider applicability, beyond MS, using magnetic resonance imaging (MRI) data from colorectal xenograft models.

RESULTS

Simulations and MALDI spectra of a stroke-damaged rat brain show MS signals from pathological tissue can be quantified. MRI diffusion data of control and radiotherapy-treated tumours further show high sensitivity hypothesis testing for treatment effects. Successful χ2 and degrees-of-freedom are computed, allowing null-hypothesis thresholding at high levels of confidence.

AVAILABILITY AND IMPLEMENTATION

Open-source image analysis software available from TINA Vision, www.tina-vision.net.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A new method for the high-precision assessment of tumor changes in response to treatment

Motivation

Imaging demonstrates that preclinical and human tumors are heterogeneous, i.e. a single tumor can exhibit multiple regions that behave differently during both development and also in response to treatment. The large variations observed in control group, tumors can obscure detection of significant therapeutic effects due to the ambiguity in attributing causes of change. This can hinder development of effective therapies due to limitations in experimental design rather than due to therapeutic failure. An improved method to model biological variation and heterogeneity in imaging signals is described. Specifically, linear Poisson modeling (LPM) evaluates changes in apparent diffusion co-efficient between baseline and 72 h after radiotherapy, in two xenograft models of colorectal cancer. The statistical significance of measured changes is compared to those attainable using a conventional t-test analysis on basic apparent diffusion co-efficient distribution parameters.

Results

When LPMs were applied to treated tumors, the LPMs detected highly significant changes. The analyses were significant for all tumors, equating to a gain in power of 4-fold (i.e. equivalent to having a sample size 16 times larger), compared with the conventional approach. In contrast, highly significant changes are only detected at a cohort level using t-tests, restricting their potential use within personalized medicine and increasing the number of animals required during testing. Furthermore, LPM enabled the relative volumes of responding and non-responding tissue to be estimated for each xenograft model. Leave-one-out analysis of the treated xenografts provided quality control and identified potential outliers, raising confidence in LPM data at clinically relevant sample sizes.

Availability and implementation

TINA Vision open source software is available from www.tina-vision.net.

Supplementary information

Supplementary data are available at Bioinformatics online.