Development of a complete advanced computational workflow for high-resolution LDI-MS metabolomics imaging data processing and visualization

La imatge per espectrometria de masses (MSI) mapeja la distribució espacial de les molècules en una mostra. Això permet extreure informació Metabolòmica espacialment corralada d'una secció de teixit. MSI no s'usa àmpliament en la metabolòmica espacial a causa de diverses limitacions relacionades amb les matrius MALDI, incloent la generació d'ions que interfereixen en el rang de masses més baix i la difusió lateral dels compostos. Hem desenvolupat un flux de treball que millora l'adquisició de metabòlits en un instrument MALDI utilitzant un "sputtering" per dipositar una nano-capa d'Au directament sobre el teixit. Això minimitza la interferència dels senyals del "background" alhora que permet resolucions espacials molt altes. S'ha desenvolupat un paquet R per a la visualització d'imatges i processament de les dades MSI, tot això mitjançant una implementació optimitzada per a la gestió de la memòria i la programació concurrent. A més, el programari desenvolupat inclou també un algoritme per a l'alineament de masses que millora la precisió de massa.La imagen por espectrometría de masas (MSI) mapea la distribución espacial de las moléculas en una muestra. Esto permite extraer información metabolòmica espacialmente corralada de una sección de tejido. MSI no se usa ampliamente en la metabolòmica espacial debido a varias limitaciones relacionadas con las matrices MALDI, incluyendo la generación de iones que interfieren en el rango de masas más bajo y la difusión lateral de los compuestos. Hemos desarrollado un flujo de trabajo que mejora la adquisición de metabolitos en un instrumento MALDI utilizando un “sputtering” para depositar una nano-capa de Au directamente sobre el tejido. Esto minimiza la interferencia de las señales del “background” a la vez que permite resoluciones espaciales muy altas. Se ha desarrollado un paquete R para la visualización de imágenes y procesado de los datos MSI, todo ello mediante una implementación optimizada para la gestión de la memoria y la programación concurrente. Además, el software desarrollado incluye también un algoritmo para el alineamiento de masas que mejora la precisión de masa.Mass spectrometry imaging (MSI) maps the spatial distributions of molecules in a sample. This allows extracting spatially-correlated metabolomics information from tissue sections. MSI is not widely used in spatial metabolomics due to several limitations related with MALDI matrices, including the generation of interfering ions and in the low mass range and the lateral compound delocalization. We developed a workflow to improve the acquisition of metabolites using a MALDI instrument. We sputter an Au nano-layer directly onto the tissue section enabling the acquisition of metabolites with minimal interference of background signals and ultra-high spatial resolution. We developed an R package for image visualization and MSI data processing, which is optimized to manage datasets larger than computer’s memory using a mutli-threaded implementation. Moreover, our software includes a label-free mass alignment algorithm for mass accuracy enhancement

Unsupervised neural networks as a support tool for pathology diagnosis in MALDI-MSI experiments:A case study on thyroid biopsies

Artificial intelligence is getting a foothold in medicine for disease screening and diagnosis. While typical machine learning methods require large labeled datasets for training and validation, their application is limited in clinical fields since ground truth information can hardly be obtained on a sizeable cohort of patients. Unsupervised neural networks - such as Self-Organizing Maps (SOMs) - represent an alternative approach to identifying hidden patterns in biomedical data. Here we investigate the feasibility of SOMs for the identification of malignant and non-malignant regions in liquid biopsies of thyroid nodules, on a patient-specific basis. MALDI-ToF (Matrix Assisted Laser Desorption Ionization -Time of Flight) mass spectrometry-imaging (MSI) was used to measure the spectral profile of bioptic samples. SOMs were then applied for the analysis of MALDI-MSI data of individual patients' samples, also testing various pre-processing and agglomerative clustering methods to investigate their impact on SOMs' discrimination efficacy. The final clustering was compared against the sample's probability to be malignant, hyperplastic or related to Hashimoto thyroiditis as quantified by multinomial regression with LASSO. Our results show that SOMs are effective in separating the areas of a sample containing benign cells from those containing malignant cells. Moreover, they allow to overlap the different areas of cytological glass slides with the corresponding proteomic profile image, and inspect the specific weight of every cellular component in bioptic samples. We envision that this approach could represent an effective means to assist pathologists in diagnostic tasks, avoiding the need to manually annotate cytological images and the effort in creating labeled datasets

Classification of cancer cell lines using matrix-assisted laser desorption/ionization time-of-flight mass spectrometry and statistical analysis

Over the past decade, matrix-assisted laser desorption/ionization time‑of‑flight mass spectrometry (MALDI‑TOF MS) has been established as a valuable platform for microbial identification, and it is also frequently applied in biology and clinical studies to identify new markers expressed in pathological conditions. The aim of the present study was to assess the potential of using this approach for the classification of cancer cell lines as a quantifiable method for the proteomic profiling of cellular organelles. Intact protein extracts isolated from different tumor cell lines (human and murine) were analyzed using MALDI‑TOF MS and the obtained mass lists were processed using principle component analysis (PCA) within Bruker Biotyper® software. Furthermore, reference spectra were created for each cell line and were used for classification. Based on the intact protein profiles, we were able to differentiate and classify six cancer cell lines: two murine melanoma (B16‑F0 and B164A5), one human melanoma (A375), two human breast carcinoma (MCF7 and MDA‑MB‑231) and one human liver carcinoma (HepG2). The cell lines were classified according to cancer type and the species they originated from, as well as by their metastatic potential, offering the possibility to differentiate non‑invasive from invasive cells. The obtained results pave the way for developing a broad‑based strategy for the identification and classification of cancer cell

Differentiating tumor heterogeneity in formalin‐fixed paraffin‐embedded (FFPE) prostate adenocarcinoma tissues using principal component analysis of matrix‐assisted laser desorption/ionization imaging mass spectral data

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/135165/1/rcm7776.pdfhttp://deepblue.lib.umich.edu/bitstream/2027.42/135165/2/rcm7776_am.pd

Cytomolecular Classification of Thyroid Nodules Using Fine-Needle Washes Aspiration Biopsies

Fine-needle aspiration biopsies (FNA) represent the gold standard to exclude the malignant nature of thyroid nodules. After cytomorphology, 20-30% of cases are deemed "indeterminate for malignancy" and undergo surgery. However, after thyroidectomy, 70-80% of these nodules are benign. The identification of tools for improving FNA's diagnostic performances is explored by matrix-assisted laser-desorption ionization mass spectrometry imaging (MALDI-MSI). A clinical study was conducted in order to build a classification model for the characterization of thyroid nodules on a large cohort of 240 samples, showing that MALDI-MSI can be effective in separating areas with benign/malignant cells. The model had optimal performances in the internal validation set (n = 70), with 100.0% (95% CI = 83.2-100.0%) sensitivity and 96.0% (95% CI = 86.3-99.5%) specificity. The external validation (n = 170) showed a specificity of 82.9% (95% CI = 74.3-89.5%) and a sensitivity of 43.1% (95% CI = 30.9-56.0%). The performance of the model was hampered in the presence of poor and/or noisy spectra. Consequently, restricting the evaluation to the subset of FNAs with adequate cellularity, sensitivity improved up to 76.5% (95% CI = 58.8-89.3). Results also suggest the putative role of MALDI-MSI in routine clinical triage, with a three levels diagnostic classification that accounts for an indeterminate gray zone of nodules requiring a strict follow-up

A framework for the statistical analysis of mass spectrometry imaging experiments

Mass spectrometry (MS) imaging is a powerful investigation technique for a wide range of biological applications such as molecular histology of tissue, whole body sections, and bacterial films , and biomedical applications such as cancer diagnosis. MS imaging visualizes the spatial distribution of molecular ions in a sample by repeatedly collecting mass spectra across its surface, resulting in complex, high-dimensional imaging datasets. Two of the primary goals of statistical analysis of MS imaging experiments are classification (for supervised experiments), i.e. assigning pixels to pre-defined classes based on their spectral profiles, and segmentation (for unsupervised experiments), i.e. assigning pixels to newly discovered segments with relatively homogenous and distinct spectral profiles. To accomplish these goals, this research provides both statistical methods and statistical computing tools. First, we propose a novel spatial shrunken centroids framework for performing classification and segmentation of MS imaging experiments with feature selection. Spatial shrunken centroids combines spatial smoothing with statistical regularization in a model-based framework appropriate for both supervised and unsupervised settings. Second, we provide Cardinal, a free and open-source R package for processing, visualization, and statistical analysis of MS imaging experiments. Cardinal is the first R package designed specifically for MS imaging, and the first software for MS imaging that focuses specifically on experiments and statistical analysis. In addition to providing tools for statistical analysis, it also provides infrastructure to enable other statisticians to more easily develop new methods for MS imaging experiments. Lastly, to enable scalability of Cardinal to larger-than-memory datasets, we provide matter, a free and open-source R package for statistical computing with structured datasets-on-disk, such as MS imaging data files. Together, spatial shrunken centroids, Cardinal, and matter aim to allow scalable statistical analysis for high-resolution, high-throughput MS imaging experiments

Intact cell mass spectrometry coupled with machine learning reveals minute changes induced by single gene silencing

Intact (whole) cell MALDI TOF mass spectrometry is a commonly used tool in clinical microbiology for several decades. Recently it was introduced to analysis of eukaryotic cells, including cancer and stem cells. Besides targeted metabolomic and proteomic applications, the intact cell MALDI TOF mass spectrometry provides a sufficient sensitivity and specificity to discriminate cell types, isogenous cell lines or even the metabolic states. This makes the intact cell MALDI TOF mass spectrometry a promising tool for quality control in advanced cell cultures with a potential to reveal batch-to-batch variation, aberrant clones, or unwanted shifts in cell phenotype. However, cellular alterations induced by change in expression of a single gene has not been addressed by intact cell mass spectrometry yet. In this work we used a well-characterized human ovarian cancer cell line SKOV3 with silenced expression of a tumor suppressor candidate 3 gene (TUSC3). TUSC3 is involved in co-translational N-glycosylation of proteins with well-known global impact on cell phenotype. Altogether, this experimental design represents a highly suitable model for optimization of intact cell mass spectrometry and analysis of spectral data. Here we investigated five machine learning algorithms (k-nearest neighbors, decision tree, random forest, partial least squares discrimination, and artificial neural network) and optimized their performance either in pure populations or in two-component mixtures composed of cells with normal or silenced expression of TUSC3. All five algorithms reached accuracy over 90 % and were able to reveal even subtle changes in mass spectra corresponding to alterations of TUSC3 expression. In summary, we demonstrate that spectral fingerprints generated by intact cell MALDI-TOF mass spectrometry coupled to a machine learning classifier can reveal minute changes induced by alteration of a single gene, and therefore contribute to the portfolio of quality control applications in routine cell and tissue cultures

A mass spectrometry imaging based approach for prognosis prediction in UICC stage I/II colon cancer

Simple Summary Tumor treatment is heavily dictated by the tumor progression status. However, in colon cancer, it is difficult to predict disease progression in the early stages. In this study, we have employed a proteomic analysis using matrix-assisted laser desorption/ionization mass spectrometry imaging (MALDI-MSI). MALDI-MSI is a technique that measures the molecular content of (tumor) tissue. We analyzed tumor samples of 276 patients. If the patients developed distant metastasis, they were considered to have a more aggressive tumor type than the patients that did not. In this comparative study, we have developed bioinformatics methods that can predict the tendency of tumor progression and advance a couple of molecules that could be used as prognostic markers of colon cancer. The prediction of tumor progression can help to choose a more adequate treatment for each individual patient. Abstract Currently, pathological evaluation of stage I/II colon cancer, following the Union Internationale Contre Le Cancer (UICC) guidelines, is insufficient to identify patients that would benefit from adjuvant treatment. In our study, we analyzed tissue samples from 276 patients with colon cancer utilizing mass spectrometry imaging. Two distinct approaches are herein presented for data processing and analysis. In one approach, four different machine learning algorithms were applied to predict the tendency to develop metastasis, which yielded accuracies over 90% for three of the models. In the other approach, 1007 m/z features were evaluated with regards to their prognostic capabilities, yielding two m/z features as promising prognostic markers. One feature was identified as a fragment from collagen (collagen 3A1), hinting that a higher collagen content within the tumor is associated with poorer outcomes. Identification of proteins that reflect changes in the tumor and its microenvironment could give a very much-needed prediction of a patient’s prognosis, and subsequently assist in the choice of a more adequate treatment