Evaluation and comparison of unsupervised methods for the extraction of spatial patterns from mass spectrometry imaging data (MSI).

For the extraction of spatially important regions from mass spectrometry imaging (MSI) data, different clustering methods have been proposed. These clustering methods are based on certain assumptions and use different criteria to assign pixels into different classes. For high-dimensional MSI data, the curse of dimensionality also limits the performance of clustering methods which are usually overcome by pre-processing the data using dimension reduction techniques. In summary, the extraction of spatial patterns from MSI data can be done using different unsupervised methods, but the robust evaluation of clustering results is what is still missing. In this study, we have performed multiple simulations on synthetic and real MSI data to validate the performance of unsupervised methods. The synthetic data were simulated mimicking important spatial and statistical properties of real MSI data. Our simulation results confirmed that K-means clustering with correlation distance and Gaussian Mixture Modeling clustering methods give optimal performance in most of the scenarios. The clustering methods give efficient results together with dimension reduction techniques. From all the dimension techniques considered here, the best results were obtained with the minimum noise fraction (MNF) transform. The results were confirmed on both synthetic and real MSI data. However, for successful implementation of MNF transform the MSI data requires to be of limited dimensions.

A methodological approach to correlate tumor heterogeneity with drug distribution profile in mass spectrometry imaging data.

BACKGROUND: Drug mass spectrometry imaging (MSI) data contain knowledge about drug and several other molecular ions present in a biological sample. However, a proper approach to fully explore the potential of such type of data is still missing. Therefore, a computational pipeline that combines different spatial and non-spatial methods is proposed to link the observed drug distribution profile with tumor heterogeneity in solid tumor. Our data analysis steps include pre-processing of MSI data, cluster analysis, drug local indicators of spatial association (LISA) map, and ions selection. RESULTS: The number of clusters identified from different tumor tissues. The spatial homogeneity of the individual cluster was measured using a modified version of our drug homogeneity method. The clustered image and drug LISA map were simultaneously analyzed to link identified clusters with observed drug distribution profile. Finally, ions selection was performed using the spatially aware method. CONCLUSIONS: In this paper, we have shown an approach to correlate the drug distribution with spatial heterogeneity in untargeted MSI data. Our approach is freely available in an R package 'CorrDrugTumorMSI'.

Past-in-the-Future. Peak detection improves targeted mass spectrometry imaging.

Mass spectrometry imaging is a valuable tool for visualizing the localization of drugs in tissues, a critical issue especially in cancer pharmacology where treatment failure may depend on poor drug distribution within the tumours. Proper preprocessing procedures are mandatory to obtain quantitative data of drug distribution in tumours, even at low intensity, through reliable ion peak identification and integration. We propose a simple preprocessing and quantification pipeline. This pipeline was designed starting from classical peak integration methods, developed when "microcomputers" became available for chromatography, now applied to MSI. This pre-processing approach is based on a novel method using the fixed mass difference between the analyte and its 5 d derivatives to set up a mass range gate. We demonstrate the use of this pipeline for the evaluating the distribution of the anticancer drug paclitaxel in tumour sections. The procedure takes advantage of a simple peak analysis and allows to quantify the drug concentration in each pixel with a limit of detection below 0.1 pmol mm-2 or 10 µg g-1. Quantitative images of paclitaxel distribution in different tumour models were obtained and average paclitaxel concentrations were compared with HPLC measures in the same specimens, showing <20% difference. The scripts are developed in Python and available through GitHub, at github.com/FrancescaFalcetta/Imaging_of_drugs_distribution_and_quantifications.git.

Drug-Homogeneity Index in Mass-Spectrometry Imaging.

Enhancing drug penetration in solid tumors is an interesting clinical issue of considerable importance. In preclinical research, mass-spectrometry imaging is a promising technique for visualizing drug distribution in tumors under different treatment conditions and its application in this field is rapidly increasing. However, in view of the huge variability among MSI data sets, drug homogeneity is usually manually assessed by an expert, and this approach is biased by interobserver variability and lacks reproducibility. We propose a new texture-based feature, the drug-homogeneity index (DHI), which provides an objective, automated measure of drug homogeneity in MSI data. A simulation study on synthetic data sets showed that previously known texture features do not give an accurate picture of intratumor drug-distribution patterns and are easily influenced by the tumor-tissue morphology. The DHI has been used to study the distribution profile of the anticancer drug paclitaxel in various xenograft models, which were either pretreated or not pretreated with antiangiogenesis compounds. The conclusion is that drug homogeneity is better in the pretreated condition, which is in agreement with previous experimental findings published by our group. This study shows that DHI could be useful in preclinical studies as a new parameter for the evaluation of protocols for better drug penetration.