Walid M. Abdelmoula

Automatic Registration of Mass Spectrometry Imaging Data Sets to the Allen Brain Atlas

Mass spectrometry imaging holds great potential for understanding the molecular basis of neurological disease. Several key studies have demonstrated its ability to uncover disease-related biomolecular changes in rodent models of disease, even if highly localized or invisible to established histological methods. The high analytical reproducibility necessary for the biomedical application of mass spectrometry imaging means it is widely developed in mass spectrometry laboratories. However, many lack the expertise to correctly annotate the complex anatomy of brain tissue, or have the capacity to analyze the number of animals required in preclinical studies, especially considering the significant variability in sizes of brain regions. To address this issue, we have developed a pipeline to automatically map mass spectrometry imaging data sets of mouse brains to the Allen Brain Reference Atlas, which contains publically available data combining gene expression with brain anatomical locations. Our pipeline enables facile and rapid interanimal comparisons by first testing if each animals tissue section was sampled at a similar location and enabling the extraction of the biomolecular signatures from specific brain regions.

Large-Scale Mass Spectrometry Imaging Investigation of Consequences of Cortical Spreading Depression in a Transgenic Mouse Model of Migraine

AbstractCortical spreading depression (CSD) is the electrophysiological correlate of migraine aura. Transgenic mice carrying the R192Q missense mutation in the Cacna1a gene, which in patients causes familial hemiplegic migraine type 1 (FHM1), exhibit increased propensity to CSD. Herein, mass spectrometry imaging (MSI) was applied for the first time to an animal cohort of transgenic and wild type mice to study the biomolecular changes following CSD in the brain. Ninety-six coronal brain sections from 32 mice were analyzed by MALDI-MSI. All MSI datasets were registered to the Allen Brain Atlas reference atlas of the mouse brain so that the molecular signatures of distinct brain regions could be compared. A number of metabolites and peptides showed substantial changes in the brain associated with CSD. Among those, different mass spectral features showed significant (t-test, P < 0.05) changes in the cortex, 146 and 377 Da, and in the thalamus, 1820 and 1834 Da, of the CSD-affected hemisphere of FHM1 R192Q mice. Our findings reveal CSD- and genotype-specific molecular changes in the brain of FHM1 transgenic mice that may further our understanding about the role of CSD in migraine pathophysiology. The results also demonstrate the utility of aligning MSI datasets to a common reference atlas for large-scale MSI investigations. Graphical Abstractᅟ

Data-driven identification of prognostic tumor subpopulations using spatially mapped t-SNE of mass spectrometry imaging data

Significance Mass spectrometry imaging provides untargeted spatiomolecular information necessary to uncover molecular intratumor heterogeneity. The challenge has been to identify those tumor subpopulations that drive patient outcomes within the highly complex datasets (hyperdimensional data, intratumor heterogeneity, and patient variation). Here we report an automatic, unbiased pipeline to nonlinearly map the hyperdimensional data into a 3D space, and identify molecularly distinct, clinically relevant tumor subpopulations. We demonstrate this pipeline’s ability to uncover subpopulations statistically associated with patient survival in primary tumors of gastric cancer and with metastasis in primary tumors of breast cancer. The identification of tumor subpopulations that adversely affect patient outcomes is essential for a more targeted investigation into how tumors develop detrimental phenotypes, as well as for personalized therapy. Mass spectrometry imaging has demonstrated the ability to uncover molecular intratumor heterogeneity. The challenge has been to conduct an objective analysis of the resulting data to identify those tumor subpopulations that affect patient outcome. Here we introduce spatially mapped t-distributed stochastic neighbor embedding (t-SNE), a nonlinear visualization of the data that is able to better resolve the biomolecular intratumor heterogeneity. In an unbiased manner, t-SNE can uncover tumor subpopulations that are statistically linked to patient survival in gastric cancer and metastasis status in primary tumors of breast cancer.

Automatic Generic Registration of Mass Spectrometry Imaging Data to Histology Using Nonlinear Stochastic Embedding

The combination of mass spectrometry imaging and histology has proven a powerful approach for obtaining molecular signatures from specific cells/tissues of interest, whether to identify biomolecular changes associated with specific histopathological entities or to determine the amount of a drug in specific organs/compartments. Currently there is no software that is able to explicitly register mass spectrometry imaging data spanning different ionization techniques or mass analyzers. Accordingly, the full capabilities of mass spectrometry imaging are at present underexploited. Here we present a fully automated generic approach for registering mass spectrometry imaging data to histology and demonstrate its capabilities for multiple mass analyzers, multiple ionization sources, and multiple tissue types.

Histology-Guided High-Resolution Matrix-Assisted Laser Desorption Ionization Mass Spectrometry Imaging.

Mass spectrometry imaging (MSI) is widely used for clinical research because when combined with histopathological analysis the molecular signatures of specific cells/regions can be extracted from the often-complex histologies of pathological tissues. The ability of MSI to stratify patients according to disease, prognosis, and response is directly attributable to this cellular specificity. MSI developments are increasingly focused on further improving specificity, through higher spatial resolution to better localize the signals or higher mass resolution to better resolve molecular ions. Higher spatial/mass resolution leads to increased data size and longer data acquisition times. For clinical applications, which analyze large series of patient tissues, this poses a challenge to keep data load and acquisition time manageable. Here we report a new tool to perform histology guided MSI; instead of analyzing large parts of each tissue section the histology from adjacent tissue sections is used to focus the analysis on the areas of interest, e.g., comparable cell types in different patient tissues, thereby minimizing data acquisition time and data load. The histology tissue section is annotated and then automatically registered to the MSI-prepared tissue section; the registration transformation is then applied to the annotations, enabling them to be used to define the MSI measurement regions. Using a series of formalin-fixed, paraffin-embedded human myxoid liposarcoma tissues, we demonstrate an 80% reduction of data load and acquisition time, thereby enabling high resolution (mass or spatial) to be more readily applied to clinical research. The software is freely available for download.

Segmentation of Choroidal Neovascularization in Fundus Fluorescein Angiograms

Choroidal neovascularization (CNV) is a common manifestation of age-related macular degeneration (AMD). It is characterized by the growth of abnormal blood vessels in the choroidal layer causing blurring and deterioration of the vision. In late stages, these abnormal vessels can rupture the retinal layers causing complete loss of vision at the affected regions. Determining the CNV size and type in fluorescein angiograms is required for proper treatment and prognosis of the disease. Computer-aided methods for CNV segmentation is needed not only to reduce the burden of manual segmentation but also to reduce inter- and intraobserver variability. In this paper, we present a framework for segmenting CNV lesions based on parametric modeling of the intensity variation in fundus fluorescein angiograms. First, a novel model is proposed to describe the temporal intensity variation at each pixel in image sequences acquired by fluorescein angiography. The set of model parameters at each pixel are used to segment the image into regions of homogeneous parameters. Preliminary results on datasets from 21 patients with Wet-AMD show the potential of the method to segment CNV lesions in close agreement with the manual segmentation.

Postmortem MRI and histology demonstrate differential iron accumulation and cortical myelin organization in early- and late-onset Alzheimer's disease

Previous MRI studies reported cortical iron accumulation in early-onset (EOAD) compared to late-onset (LOAD) Alzheimer disease patients. However, the pattern and origin of iron accumulation is poorly understood. This study investigated the histopathological correlates of MRI contrast in both EOAD and LOAD. T2*-weighted MRI was performed on postmortem frontal cortex of controls, EOAD, and LOAD. Images were ordinally scored using predefined criteria followed by histology. Nonlinear histology-MRI registration was used to calculate pixel-wise spatial correlations based on the signal intensity. EOAD and LOAD were distinguishable based on 7T MRI from controls and from each other. Histology-MRI correlation analysis of the pixel intensities showed that the MRI contrast is best explained by increased iron accumulation and changes in cortical myelin, whereas amyloid and tau showed less spatial correspondence with T2*-weighted MRI. Neuropathologically, subtypes of Alzheimers disease showed different patterns of iron accumulation and cortical myelin changes independent of amyloid and tau that may be detected by high-field susceptibility-based MRI.

Whole-brain microscopy meets in vivo neuroimaging: techniques, benefits, and limitations

Magnetic resonance imaging, positron emission tomography, and optical imaging have emerged as key tools to understand brain function and neurological disorders in preclinical mouse models. They offer the unique advantage of monitoring individual structural and functional changes over time. What remained unsolved until recently was to generate whole-brain microscopy data which can be correlated to the 3D in vivo neuroimaging data. Conventional histological sections are inappropriate especially for neuronal tracing or the unbiased screening for molecular targets through the whole brain. As part of the European Society for Molecular Imaging (ESMI) meeting 2016 in Utrecht, the Netherlands, we addressed this issue in the Molecular Neuroimaging study group meeting. Presentations covered new brain clearing methods, light sheet microscopes for large samples, and automatic registration of microscopy to in vivo imaging data. In this article, we summarize the discussion; give an overview of the novel techniques; and discuss the practical needs, benefits, and limitations.

Interactive Visual Exploration of 3D Mass Spectrometry Imaging Data Using Hierarchical Stochastic Neighbor Embedding Reveals Spatiomolecular Structures at Full Data Resolution

Technological advances in mass spectrometry imaging (MSI) have contributed to growing interest in 3D MSI. However, the large size of 3D MSI data sets has made their efficient analysis and visualization and the identification of informative molecular patterns computationally challenging. Hierarchical stochastic neighbor embedding (HSNE), a nonlinear dimensionality reduction technique that aims at finding hierarchical and multiscale representations of large data sets, is a recent development that enables the analysis of millions of data points, with manageable time and memory complexities. We demonstrate that HSNE can be used to analyze large 3D MSI data sets at full mass spectral and spatial resolution. To benchmark the technique as well as demonstrate its broad applicability, we have analyzed a number of publicly available 3D MSI data sets, recorded from various biological systems and spanning different mass-spectrometry ionization techniques. We demonstrate that HSNE is able to rapidly identify regions of interest within these large high-dimensionality data sets as well as aid the identification of molecular ions that characterize these regions of interest; furthermore, through clearly separating measurement artifacts, the HSNE analysis exhibits a degree of robustness to measurement batch effects, spatially correlated noise, and mass spectral misalignment.

Quantitative assessment of age-related macular degeneration using parametric modeling of the leakage transfer function: Preliminary results

Age-related macular degeneration (AMD) is a major cause of blindness and visual impairment in older adults. The wet form of the disease is characterized by abnormal blood vessels forming a choroidal neovascular membrane (CNV), that result in destruction of normal architecture of the retina. Current evaluation and follow up of wet AMD include subjective evaluation of Fluorescein Angiograms (FA) to determine the activity of the lesion and monitor the progression or regression of the disease. However, this subjective evaluation prevents accurate monitoring of the disease progression or regression in response to a pharmacologic agent. In this work, we present a method that allows objective assessment of the activity of a CNV lesion which can be statistically compared across different patient and time points. The method is based on a hypothesis that the discrepancy in the time-intensity signals among the diseased and normal retinal areas are due to an implicit transfer function whose parameters can be used to characterize the retina. The method begins with parametric modeling of the temporal variation of the lesion and background intensities. Then, the values of the model parameters are used to evaluate the change in the activity of the disease. Preliminary results on five datasets show that the calculated parameters are highly correlated with the Visual Acuity (VA) of the patients.