Improved Mass Calibration in MALDI MSI Using Neural Network-Based Recalibration.

Matrix-assisted laser desorption/ionization mass spectrometry imaging (MALDI MSI) is a powerful imaging method for generating molecular maps of biological samples and has numerous applications in biomedical research. A key challenge in MALDI MSI is to reliably map observed mass peaks to theoretical masses, which can be difficult due to mass shifts that occur during the measurement process. In this paper, we propose MassShiftNet, a novel self-supervised learning framework for mass recalibration. We train a neural network on a data dependent and specifically augmented training data set to directly estimate and correct the mass shift in the observed spectra. In our evaluation, we show that this method is both able to reduce the absolute mass error and to increase the relative mass alignment between peptide MSI spectra acquired from FFPE-fixated tissue using a MALDI time-of-flight (TOF) instrument.

Spatial probabilistic mapping of metabolite ensembles in mass spectrometry imaging.

Mass spectrometry imaging vows to enable simultaneous spatially resolved investigation of hundreds of metabolites in tissues, but it primarily relies on traditional ion images for non-data-driven metabolite visualization and analysis. The rendering and interpretation of ion images neither considers nonlinearities in the resolving power of mass spectrometers nor does it yet evaluate the statistical significance of differential spatial metabolite abundance. Here, we outline the computational framework moleculaR ( https://github.com/CeMOS-Mannheim/moleculaR ) that is expected to improve signal reliability by data-dependent Gaussian-weighting of ion intensities and that introduces probabilistic molecular mapping of statistically significant nonrandom patterns of relative spatial abundance of metabolites-of-interest in tissue. moleculaR also enables cross-tissue statistical comparisons and collective molecular projections of entire biomolecular ensembles followed by their spatial statistical significance evaluation on a single tissue plane. It thereby fosters the spatially resolved investigation of ion milieus, lipid remodeling pathways, or complex scores like the adenylate energy charge within the same image.

Multimodal Lung Cancer Subtyping Using Deep Learning Neural Networks on Whole Slide Tissue Images and MALDI MSI.

Artificial intelligence (AI) has shown potential for facilitating the detection and classification of tumors. In patients with non-small cell lung cancer, distinguishing between the most common subtypes, adenocarcinoma (ADC) and squamous cell carcinoma (SqCC), is crucial for the development of an effective treatment plan. This task, however, may still present challenges in clinical routine. We propose a two-modality, AI-based classification algorithm to detect and subtype tumor areas, which combines information from matrix-assisted laser desorption/ionization (MALDI) mass spectrometry imaging (MSI) data and digital microscopy whole slide images (WSIs) of lung tissue sections. The method consists of first detecting areas with high tumor cell content by performing a segmentation of the hematoxylin and eosin-stained (H&E-stained) WSIs, and subsequently classifying the tumor areas based on the corresponding MALDI MSI data. We trained the algorithm on six tissue microarrays (TMAs) with tumor samples from N = 232 patients and used 14 additional whole sections for validation and model selection. Classification accuracy was evaluated on a test dataset with another 16 whole sections. The algorithm accurately detected and classified tumor areas, yielding a test accuracy of 94.7% on spectrum level, and correctly classified 15 of 16 test sections. When an additional quality control criterion was introduced, a 100% test accuracy was achieved on sections that passed the quality control (14 of 16). The presented method provides a step further towards the inclusion of AI and MALDI MSI data into clinical routine and has the potential to reduce the pathologist's work load. A careful analysis of the results revealed specific challenges to be considered when training neural networks on data from lung cancer tissue.

Multicenter Evaluation of Tissue Classification by Matrix-Assisted Laser Desorption/Ionization Mass Spectrometry Imaging.

Many studies have demonstrated that tissue phenotyping (tissue typing) based on mass spectrometric imaging data is possible; however, comprehensive studies assessing variation and classifier transferability are largely lacking. This study evaluated the generalization of tissue classification based on Matrix Assisted Laser Desorption/Ionization (MALDI) mass spectrometric imaging (MSI) across measurements performed at different sites. Sections of a tissue microarray (TMA) consisting of different formalin-fixed and paraffin-embedded (FFPE) human tissue samples from different tumor entities (leiomyoma, seminoma, mantle cell lymphoma, melanoma, breast cancer, and squamous cell carcinoma of the lung) were prepared and measured by MALDI-MSI at different sites using a standard protocol (SOP). Technical variation was deliberately introduced on two separate measurements via a different sample preparation protocol and a MALDI Time of Flight mass spectrometer that was not tuned to optimal performance. Using standard data preprocessing, a classification accuracy of 91.4% per pixel was achieved for intrasite classifications. When applying a leave-one-site-out cross-validation strategy, accuracy per pixel over sites was 78.6% for the SOP-compliant data sets and as low as 36.1% for the mistuned instrument data set. Data preprocessing designed to remove technical variation while retaining biological information substantially increased classification accuracy for all data sets with SOP-compliant data sets improved to 94.3%. In particular, classification accuracy of the mistuned instrument data set improved to 81.3% and from 67.0% to 87.8% per pixel for the non-SOP-compliant data set. We demonstrate that MALDI-MSI-based tissue classification is possible across sites when applying histological annotation and an optimized data preprocessing pipeline to improve generalization of classifications over technical variation and increasing overall robustness.

Robust subtyping of non-small cell lung cancer whole sections through MALDI mass spectrometry imaging.

Subtyping of the most common non-small cell lung cancer (NSCLC) tumor types adenocarcinoma (ADC) and squamous cell carcinoma (SqCC) is still a challenge in the clinical routine and a correct diagnosis is crucial for an adequate therapy selection. Matrix-assisted laser desorption/ionization (MALDI) mass spectrometry imaging (MSI) has shown potential for NSCLC subtyping but is subject to strong technical variability and has only been applied to tissue samples assembled in tissue microarrays (TMAs). To our knowledge, a successful transfer of a classifier from TMAs to whole sections, which are generated in the standard clinical routine, has not been presented in the literature as of yet. We introduce a classification algorithm using extensive preprocessing and a classifier (either a neural network or a linear discriminant analysis (LDA)) to robustly classify whole sections of ADC and SqCC lung tissue. The classifiers were trained on TMAs and validated and tested on whole sections. Vital for a successful application on whole sections is the extensive preprocessing and the use of whole sections for hyperparameter selection. The classification system with the neural network/LDA results in 99.0%/98.3% test accuracy on spectra level and 100.0%/100.0% test accuracy on whole section level, respectively, and, therefore, provides a powerful tool to support the pathologist's decision making process. The presented method is a step further towards a clinical application of MALDI MSI and artificial intelligence for subtyping of NSCLC tissue sections.

Deeply Supervised UNet for Semantic Segmentation to Assist Dermatopathological Assessment of Basal Cell Carcinoma.

Accurate and fast assessment of resection margins is an essential part of a dermatopathologist's clinical routine. In this work, we successfully develop a deep learning method to assist the dermatopathologists by marking critical regions that have a high probability of exhibiting pathological features in whole slide images (WSI). We focus on detecting basal cell carcinoma (BCC) through semantic segmentation using several models based on the UNet architecture. The study includes 650 WSI with 3443 tissue sections in total. Two clinical dermatopathologists annotated the data, marking tumor tissues' exact location on 100 WSI. The rest of the data, with ground-truth sectionwise labels, are used to further validate and test the models. We analyze two different encoders for the first part of the UNet network and two additional training strategies: (a) deep supervision, (b) linear combination of decoder outputs, and obtain some interpretations about what the network's decoder does in each case. The best model achieves over 96%, accuracy, sensitivity, and specificity on the Test set.

Cross-Normalization of MALDI Mass Spectrometry Imaging Data Improves Site-to-Site Reproducibility.

Matrix-assisted laser desorption/ionization mass spectrometry imaging (MALDI MSI) is an established tool for the investigation of formalin-fixed paraffin-embedded (FFPE) tissue samples and shows a high potential for applications in clinical research and histopathological tissue classification. However, the applicability of this method to serial clinical and pharmacological studies is often hampered by inevitable technical variation and limited reproducibility. We present a novel spectral cross-normalization algorithm that differs from the existing normalization methods in two aspects: (a) it is based on estimating the full statistical distribution of spectral intensities and (b) it involves applying a non-linear, mass-dependent intensity transformation to align this distribution with a reference distribution. This method is combined with a model-driven resampling step that is specifically designed for data from MALDI imaging of tryptic peptides. This method was performed on two sets of tissue samples: a single human teratoma sample and a collection of five tissue microarrays (TMAs) of breast and ovarian tumor tissue samples (N = 241 patients). The MALDI MSI data was acquired in two labs using multiple protocols, allowing us to investigate different inter-lab and cross-protocol scenarios, thus covering a wide range of technical variations. Our results suggest that the proposed cross-normalization significantly reduces such batch effects not only in inter-sample and inter-lab comparisons but also in cross-protocol scenarios. This demonstrates the feasibility of cross-normalization and joint data analysis even under conditions where preparation and acquisition protocols themselves are subject to variation.

Classification of Inflammatory Bowel Disease from Formalin-Fixed, Paraffin-Embedded Tissue Biopsies via Imaging Mass Spectrometry.

PURPOSE: Discrimination between ulcerative colitis (UC) and Crohn's disease (CD) by histologic features alone can be challenging and often leads to inaccurate initial diagnoses in inflammatory bowel disease (IBD) patients. This is mostly due to an overlap of clinical and histologic features. However, exact diagnosis is not only important for patient treatment but it also has a socioeconomic impact. It is therefore important to develop and improve diagnostic tools complementing traditional histomorphological approaches. EXPERIMENTAL DESIGN: In this retrospective proof-of-concept study, the utilization of MALDI imaging is explored in combination with multi-variate data analysis methods to classify formalin-fixed, paraffin-embedded (FFPE) colon biopsies from UC (87 biopsies, 14 patients), CD (71 biopsies, 14 patients), and normal colonic (21 biopsies, 14 patients) tissues. RESULTS: The proposed method results in an overall balanced accuracy of 85.7% on patient and of 80.4% on sample level, thus demonstrating that the assessment of IBD from FFPE tissue specimens via MALDI imaging is feasible. CONCLUSIONS AND CLINICAL RELEVANCE: The results emphasize the high potential of this method to distinguish IBD subtypes in FFPE tissue sections, which is a prerequisite for further investigations in retrospective multicenter studies, as well as for a future implementation into clinical routine.

Peak detection for MALDI mass spectrometry imaging data using sparse frame multipliers.

MALDI mass spectrometry imaging (MALDI MSI) is a spatially resolved analytical tool for biological tissue analysis by measuring mass-to-charge ratios of ionized molecules. With increasing spatial and mass resolution of MALDI MSI data, appropriate data analysis and interpretation is getting more and more challenging. A reliable separation of important peaks from noise (aka peak detection) is a prerequisite for many subsequent processing steps and should be as accurate as possible. We propose a novel peak detection algorithm based on sparse frame multipliers, which can be applied to raw MALDI MSI data without prior preprocessing. The accuracy is evaluated on a simulated data set in comparison with state-of-the-art algorithms. These results also show the proposed method's robustness to baseline and noise effects. In addition, the method is evaluated on real MALDI-TOF data sets, whereby spatial information can be included in the peak picking process. SIGNIFICANCE: The field of proteomics, in particular MALDI Imaging, encompasses huge amounts of data. The processing and preprocessing of this data in order to segment or classify spatial structures of certain peptides or isotope patterns can hence be cumbersome and includes several independent processing steps. In this work, we propose a simple peak-picking algorithm to quickly analyze large raw MALDI Imaging data sets, which has a better sensitivity than current state-of-the-art algorithms. Further, it is possible to get an overall overview of the entire data set showing the most significant and spatially localized peptide structures and, hence, contributes all data driven evaluation of MALDI Imaging data.

Using the Chemical Noise Background in MALDI Mass Spectrometry Imaging for Mass Alignment and Calibration.

Matrix-assisted laser desorption/ionization mass spectrometry imaging (MALDI MSI) is an established tool for the investigation of formalin fixed paraffin embedded (FFPE) tissue samples and shows a high potential for applications in clinical research and histopathological diagnosis. The applicability and accuracy of this method, however, heavily depends on the quality of the acquired data, and in particular mass misalignment in axial time-of-flight (TOF) MSI continues to be a serious issue. We present a mass alignment and recalibration method that is specifically designed to operate on MALDI peptide imaging data. The proposed method exploits statistical properties of the characteristic chemical noise background observed in peptide imaging experiments. By comparing these properties to a theoretical peptide mass model, the effective mass shift of each spectrum is estimated and corrected. The method was evaluated on a cohort of 31 FFPE tissue samples, pursuing a statistical validation approach to estimate both the reduction of relative misalignment, as well as the increase in absolute mass accuracy. Our results suggest that a relative mass precision of approximately 5 ppm and an absolute accuracy of approximately 20 ppm are achievable using our method.

Supervised non-negative matrix factorization methods for MALDI imaging applications.

MOTIVATION: Non-negative matrix factorization (NMF) is a common tool for obtaining low-rank approximations of non-negative data matrices and has been widely used in machine learning, e.g. for supporting feature extraction in high-dimensional classification tasks. In its classical form, NMF is an unsupervised method, i.e. the class labels of the training data are not used when computing the NMF. However, incorporating the classification labels into the NMF algorithms allows to specifically guide them toward the extraction of data patterns relevant for discriminating the respective classes. This approach is particularly suited for the analysis of mass spectrometry imaging (MSI) data in clinical applications, such as tumor typing and classification, which are among the most challenging tasks in pathology. Thus, we investigate algorithms for extracting tumor-specific spectral patterns from MSI data by NMF methods. RESULTS: In this article, we incorporate a priori class labels into the NMF cost functional by adding appropriate supervised penalty terms. Numerical experiments on a MALDI imaging dataset confirm that the novel supervised NMF methods lead to significantly better classification accuracy and stability as compared with other standard approaches. AVAILABILITY AND IMPLEMENTATON: https://gitlab.informatik.uni-bremen.de/digipath/Supervised_NMF_Methods_for_MALDI.git. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Targeted Feature Extraction in MALDI Mass Spectrometry Imaging to Discriminate Proteomic Profiles of Breast and Ovarian Cancer.

PURPOSE: To develop a mass spectrometry imaging (MSI) based workflow for extracting m/z values related to putative protein biomarkers and using these for reliable tumor classification. EXPERIMENTAL DESIGN: Given a list of putative breast and ovarian cancer biomarker proteins, a set of related m/z values are extracted from heterogeneous MSI datasets derived from formalin-fixed paraffin-embedded tissue material. Based on these features, a linear discriminant analysis classification model is trained to discriminate the two tumor types. RESULTS: It is shown that the discriminative power of classification models based on the extracted features is increased compared to the automatic training approach, especially when classifiers are applied to spectral data acquired under different conditions (instrument, preparation, laboratory). CONCLUSIONS AND CLINICAL RELEVANCE: Robust classification models not confounded by technical variation between MSI measurements are obtained. This supports the assumption that the classification of the respective tumor types is based on biological rather than technical differences, and that the selected features are related to the proteomic profiles of the tumor types under consideration.

Deep learning for tumor classification in imaging mass spectrometry.

Motivation: Tumor classification using imaging mass spectrometry (IMS) data has a high potential for future applications in pathology. Due to the complexity and size of the data, automated feature extraction and classification steps are required to fully process the data. Since mass spectra exhibit certain structural similarities to image data, deep learning may offer a promising strategy for classification of IMS data as it has been successfully applied to image classification. Results: Methodologically, we propose an adapted architecture based on deep convolutional networks to handle the characteristics of mass spectrometry data, as well as a strategy to interpret the learned model in the spectral domain based on a sensitivity analysis. The proposed methods are evaluated on two algorithmically challenging tumor classification tasks and compared to a baseline approach. Competitiveness of the proposed methods is shown on both tasks by studying the performance via cross-validation. Moreover, the learned models are analyzed by the proposed sensitivity analysis revealing biologically plausible effects as well as confounding factors of the considered tasks. Thus, this study may serve as a starting point for further development of deep learning approaches in IMS classification tasks. Availability and implementation: https://gitlab.informatik.uni-bremen.de/digipath/Deep_Learning_for_Tumor_Classification_in_IMS. Contact: jbehrmann@uni-bremen.de or christianetmann@uni-bremen.de. Supplementary information: Supplementary data are available at Bioinformatics online.

A new classification method for MALDI imaging mass spectrometry data acquired on formalin-fixed paraffin-embedded tissue samples.

Matrix-assisted laser desorption/ionization imaging mass spectrometry (MALDI IMS) shows a high potential for applications in histopathological diagnosis, and in particular for supporting tumor typing and subtyping. The development of such applications requires the extraction of spectral fingerprints that are relevant for the given tissue and the identification of biomarkers associated with these spectral patterns. We propose a novel data analysis method based on the extraction of characteristic spectral patterns (CSPs) that allow automated generation of classification models for spectral data. Formalin-fixed paraffin embedded (FFPE) tissue samples from N=445 patients assembled on 12 tissue microarrays were analyzed. The method was applied to discriminate primary lung and pancreatic cancer, as well as adenocarcinoma and squamous cell carcinoma of the lung. A classification accuracy of 100% and 82.8%, resp., could be achieved on core level, assessed by cross-validation. The method outperformed the more conventional classification method based on the extraction of individual m/z values in the first application, while achieving a comparable accuracy in the second. LC-MS/MS peptide identification demonstrated that the spectral features present in selected CSPs correspond to peptides relevant for the respective classification. This article is part of a Special Issue entitled: MALDI Imaging, edited by Dr. Corinna Henkel and Prof. Peter Hoffmann.

An approach to optimize sample preparation for MALDI imaging MS of FFPE sections using fractional factorial design of experiments.

A standardized workflow for matrix-assisted laser desorption/ionization imaging mass spectrometry (MALDI imaging MS) is a prerequisite for the routine use of this promising technology in clinical applications. We present an approach to develop standard operating procedures for MALDI imaging MS sample preparation of formalin-fixed and paraffin-embedded (FFPE) tissue sections based on a novel quantitative measure of dataset quality. To cover many parts of the complex workflow and simultaneously test several parameters, experiments were planned according to a fractional factorial design of experiments (DoE). The effect of ten different experiment parameters was investigated in two distinct DoE sets, each consisting of eight experiments. FFPE rat brain sections were used as standard material because of low biological variance. The mean peak intensity and a recently proposed spatial complexity measure were calculated for a list of 26 predefined peptides obtained by in silico digestion of five different proteins and served as quality criteria. A five-way analysis of variance (ANOVA) was applied on the final scores to retrieve a ranking of experiment parameters with increasing impact on data variance. Graphical abstract MALDI imaging experiments were planned according to fractional factorial design of experiments for the parameters under study. Selected peptide images were evaluated by the chosen quality metric (structure and intensity for a given peak list), and the calculated values were used as an input for the ANOVA. The parameters with the highest impact on the quality were deduced and SOPs recommended.

New techniques in CT angiography.

Computed tomographic (CT) angiography has been improved significantly with the introduction of four- to 64-section spiral CT scanners, which offer rapid acquisition of isotropic data sets. A variety of techniques have been proposed for postprocessing of the resulting images. The most widely used techniques are multiplanar reformation (MPR), thin-slab maximum intensity projection, and volume rendering. Sophisticated segmentation algorithms, vessel analysis tools based on a centerline approach, and automatic lumen boundary definition are emerging techniques; bone removal with thresholding or subtraction algorithms has been introduced. These techniques increasingly provide a quality of vessel analysis comparable to that achieved with intraarterial three-dimensional rotational angiography. Neurovascular applications for these various image postprocessing methods include steno-occlusive disease, dural sinus thrombosis, vascular malformations, and cerebral aneurysms. However, one should keep in mind the potential pitfalls of these techniques and always double-check the final results with source or MPR imaging.

CT angiography of peripheral arterial bypass grafts: Accuracy and time-effectiveness of quantitative image analysis with an automated software tool.

RATIONALE AND OBJECTIVES: Qualitative analysis of computed tomography (CT) angiography data often is limited by intra- and interobserver variability. The purpose of this study was to evaluate the time-effectiveness and accuracy of a quantitative CT angiography data analysis using automated software in comparison with qualitative axial and coronal CT image reading in patients with peripheral bypass grafts. MATERIALS AND METHODS: Twenty-eight patients with 33 saphenous bypass grafts underwent 4-channel (n = 21) and 16-channel (n = 7) CT angiography. Two readers evaluated in consensus the CT data qualitatively on axial and coronal reconstructions and with the software regarding the presence of graft stenoses, aneurysmal changes, and arteriovenous fistulas. The time for data analysis was taken and the accuracy was compared with the results from digital subtraction angiography (DSA). RESULTS: No significant difference was present between data analysis time using axial and coronal CT images (4.9 +/- 1.5 minutes) and when using the software tool (5.5 +/- 1.4 minutes). Good (kappa = 0.652) to excellent (kappa = 1.000) intermodality agreement was present between qualitative and quantitative CT analysis regarding graft-related abnormalities. Sensitivity and specificity for diagnosing stenoses, aneurysms, and fistula did not differ significantly (P > .025) between qualitative CT image reading and the automated software tool. CONCLUSIONS: CT angiography analysis of peripheral bypass grafts using an automated software tool is similar regarding time-effectiveness and accuracy when compared with qualitative CT data analysis on axial and coronal images. It may assist in determining the significance of an abnormality and can yield objective morphometric data of vessel calibers.

New vessel analysis tool for morphometric quantification and visualization of vessels in CT and MR imaging data sets.

Image processing algorithms and a prototypical research software tool have been developed for visualization and quantitative analysis of vessels in data sets from computed tomography and magnetic resonance imaging. The software is based on a sequence of processing steps, which are as follows: (a) vessel segmentation based on a region growing algorithm, (b) interactive "premasking" to optionally exclude interfering structures close to the vessels of interest, (c) distance transform-based skeletonization, (d) multiplanar reformation orthogonal to the vessel path, (e) identification of the lumen boundary on the orthogonal cross-section images, and (f) morphometric measurements. The development of the algorithmic components and the application user interface has been carried out in close cooperation with clinical users to achieve a high degree of usability and flexible support of work flow. The software has been successfully applied to the intracranial arteries, carotid arteries, and abdominal and thoracic aorta, as well as the renal, coronary, and peripheral arteries.