Mass spectrometry imaging of lipids: untargeted consensus spectra reveal spatial distributions in Niemann-Pick disease type C1[S]

Mass spectrometry imaging (MSI) is a tool to rapidly map the spatial location of analytes without the need for tagging or a reporter system. Niemann-Pick disease type C1 (NPC1) is a neurodegenerative, lysosomal storage disorder characterized by accumulation of unesterified cholesterol and sphingolipids in the endo-lysosomal system. Here, we use MSI to visualize lipids including cholesterol in cerebellar brain tissue from the NPC1 symptomatic mouse model and unaffected controls. To complement the imaging studies, a data-processing pipeline was developed to generate consensus mass spectra, thereby using both technical and biological image replicates to assess differences. The consensus spectra are used to determine true differences in lipid relative abundance; lipid distributions can be determined in an unbiased fashion without prior knowledge of location. We show the cerebellar distribution of gangliosides GM1, GM2, and GM3, including variants of lipid chain length. We also performed MALDI-MSI of cholesterol. Further analysis of lobules IV/V and X of the cerebellum gangliosides indicates regional differences. The specificity achieved highlights the power of MSI, and this new workflow demonstrates a universal approach for addressing reproducibility in imaging experiments applied to NPC1.

Additionally, lyso-sphingolipids in plasma and amniotic fluid have been quantified using targeted MS methods (33). Although whole tissue lysate lipidomics studies are invaluable, as they are sensitive, specific and can be used to obtain precise quantities, any information regarding spatial distribution is lost.
In the context of lysosomal storage disorders, few studies have reported lipid imaging using MS. Lipid mapping has been performed in a mouse model of Hunter syndrome, a carbohydrate metabolism disorder in which spatial distributions of the gangliosides GM2 and GM3 were obtained alongside immunohistochemistry (34). A method for enhancing neutral lipid mapping was shown in a mouse model of Fabry disease (35). Other lysosomal disorders with imaging studies include Farber (36), Gaucher (37,38), and Sandhoff disease (39,40). Recently, the I1061T NPC1 point mutant model was used to develop an infrared spectroscopy-MALDI imaging workflow. In this approach initial Fourier-transform infrared imaging data is used to guide MALDI imaging to allow for tissue annotation. The statistical component, including a t-test based feature and peak-picking extraction, was performed; however, this approach required specialized tools to obtain segmentation information (41). It is also important to note that while an NPC1 model was used, this study did not look at cerebellar changes given the coronal sectioning performed nor did the authors provide insight into lipid mapping related to NPC1 disease. These reports reveal the importance of spatially localizing lipids, particularly in NPC1, which has been far more troublesome for antibody-based approaches and provide an insight into disease pathology.
Given the complexity and size of MSI experiments, reproducibility is a paramount consideration in experimental design. For example, the goal of many MSI tissue profiling studies is to identify disease markers relative to their location in the sample. In order to have significance of a specific mass-to-charge ratio, any finding must be spatially, analytically, and biologically reproducible. To avoid the pursuit of "beautiful noise," various techniques in preprocessing MSI data have been implemented (42,43) and reviewed in greater detail (44). Pre-processing strategies apply mathematical and logical operations to raw spectral data such that replicates can be combined, compared, and analyzed ostensibly with minimized contributions from biological, instrumental, and random variance. The most common of these approaches is to normalize the total ion current to generate the image. The difficulty with such an approach is that it only averages signal variance and does not take into account the overall spectral integrity. It thereby assumes that each spectrum collected over a given area has a comparable number of features. This assumption may be acceptable when dealing with homogenous results of different chemical samples. However, imaging of biological replicates must be considered as a set of separate experiments, and thus even multiple acquisitions of serial sections can lead to great mass spectral variability. Eijkel et al. (45) developed an algorithm to correlate secondary ion MS and MALDI imaging data from the same tissue. The authors attempted to address the data reduction of imaging experiments by conducting spectral binning to unit resolution, which is inappropriate to do when using high-resolution mass analyzers. They performed principal component analysis and canonical correlation analysis to determine unique features in their datasets, but replicate analysis was not addressed in the workflow, thus potentially leading to an underestimation of features (45). We argue that simply publishing MSI data as representative images rather than constructs of replicate spectra underestimates the true variance between analyses, thereby potentially misleading the reader as to the significance of features. Such an approach furthermore excludes rather than incorporates the majority of data generated in the investigation.
In the present study, we sought to develop an algorithm to evaluate both biological and technical replicate imaging experiments and then identify altered lipids and map their spatial distribution using a mouse model of NPC1. We imaged the cerebellum of 7-week-old null BALB/c-Npc nih , hereafter Npc1 / (46), and control littermates, using MALDI-MSI. Relative to the disease progression, the 7-week time point displays classical phenotypes of ataxia and tremors that recapitulate the human disease. The newly developed algorithm was used to evaluate Npc1 +/+ and Npc1 / MSI spectra datasets of the cerebellum to validate reproducibility between replicates. This method provides an unbiased process to screen for unique and biologically interesting features that would not be detected from a single tissue image alone and also provides a way to form representative mass spectra for a dataset of replicates.

MATERIALS AND METHODS
All reagents were used as supplied unless otherwise noted. Purified water was obtained via a Barnstead GenPure (Thermo Fisher Scientific). All reagents were obtained from Sigma-Aldrich unless noted. Hematoxylin was purchased from Ricca Chemical Co.

Experimental model
All experiments were performed in accordance with University of Illinois at Chicago IACUC-approved protocols. Balb/c npc nih (Npc1 +/ ) mice were obtained from Jackson Laboratories (RRID:IMSR JAX:003092), and a breeding colony was maintained in our laboratory. Genotyping was performed using polymerase chain reaction as previously reported (17). The primer sequences used for genotyping were: (FWD8F) 5′-GGTGCTGGACAGC-CAAGTA-3′ and (REVINTR3) 5′-GATGGTCTGTTCTCCCATG-3′. At 7 weeks of age, control (Npc1 +/+ ) and null mutant (Npc1 / ) mice were euthanized via CO 2 asphyxiation followed by decapitation. Whole brain tissue was dissected and immediately frozen in dry ice to maintain spatial integrity and stored at 80°C.

H&E staining
Tissue staining was performed using rehydration and dehydration steps. Briefly, tissue sections on microscope slides were immersed in 95% ethanol solution for 30

Sample preparation
Frozen, intact brain was cryosectioned at 20°C using a Microm HM 525 (Thermo Fisher Scientific). Serial tissue sections of 16 m thickness were thaw-mounted directly on MALDI stainless steel target plates and optical microscope slides for H&E staining and then were placed under vacuum for 10 min, then stored at 80°C until further analysis. The MALDI stainless-steel plate was dried under vacuum for 10 min to remove residual moisture. The plate containing the tissue section was then submerged in 50 mM ammonium formate solution (Alfa Aesar) for 20 s (47). Matrix application was conducted using an artist-model airbrush to produce an even coating. The matrix 9-aminoacridine was prepared at 10 mg/ml concentration in acetone (Thermo Fisher Scientific). In an effort to have reproducible matrix application, the density of matrix applied was calculated by weighing the target plates before and after the application with the desired amount of matrix being 0.15-0.25 mg/cm 2 . The analysis of matrix deposition is provided in supplemental Fig. S1. The mass spectrum for the on-tissue analysis (supplemental Fig. S1A-D) is provided for each replicate to evaluate m/z 612.3, a matrix ion of 1,5-diaminonapthalene (48). An example image of a matrix ion is provided to demonstrate reproducible matrix application (supplemental Fig.  S1E). Slight differences are observed from the tissue compared with the stainless-steel plate, which is confirmed by recent reports of ion suppression in MALDI-MSI experiments (49). Biological replicate analyses (n = 5 Npc1 +/+ , n = 6 Npc1 / ) were conducted with n = 3-4 technical replicates (serial sections) collected for each animal. Matrix sublimation of dihydroxybenzoic acid (DHB; Sigma-Aldrich) was conducted for cholesterol imaging using a homemade sublimation apparatus. DHB matrix was dissolved in acetone (30 mg/ml) and transferred into the bottom of the sublimation apparatus, then evaporated using a stream of nitrogen gas. The preweighed target plate was then attached to the flat bottom of the apparatus using copper tape. The sand bath temperature was set to 120°C, and the trap was filled with ice slush. The internal vacuum was kept at 80 mTorr before starting the sublimation which lasts for 2 min. The target plate was weighed again and measured to obtain the density. Initial lipid assignments were made by comparing accurate mass measurements with the LIPID MAPS database (http://www.lipidmaps.org). Second, tandem MS data was acquired to yield signature fragment ions of each lipid species.

MALDI-MSI
MSI was performed using a model 4800 Plus MALDI TOF/ TOF Analyzer (Sciex) equipped with a 200 Hz Nd-YAG pulsed laser (355 nm). The approximate laser spot size was 90 µm. Data were acquired in both the positive and negative ion reflectron modes as separate experiments. Red phosphorus was used for external mass calibration as previously described (50). The number of laser shots per pixel was set at 255 and the raster distance between each pixel was set to 50 m using the 4800 Imaging Tool (version 3.2; www.maldi-msi.org, M. Stoeckli). Data processing for the MS ion images was conducted using TissueView software (Version 1.1; Sciex) or MSiReader (51). Regions-of-interest were determined based on ion intensities compared with the MALDI plate which was then compared with a representative H&E-stained section. Data analysis was carried out using the R-Programming Language (https://www.r-project.org/about.html).

LC/MS/MS
Cerebellar tissue from three Npc1 +/+ and Npc1 / animals at 7 weeks of age were homogenized in 1X PBS. Protein concentration was measured using the Pierce BCA Protein Assay (Thermo Scientific). Extraction of cerebellar lipids from aliquots of each homogenized tissue was carried out using a modified chloroform/ methanol method to extract gangliosides originally developed by Folch et. al (52). The upper aqueous layer from the chloroform/ methanol protocol was collected, and the lipid extracts were resuspended in 9:1 (v/v) methanol:toluene with volumes normalized according to the protein concentration. Lipids were separated using an Agilent EclipsePlus C18 RRHD C18 column (100 × 2.1 mm, 1.8 m) coupled to an Agilent 1290 Infinity II system. The column was maintained at 50°C at a flow rate of 250 l/min. The mobile phases consisted of solvent A, 60:40 (v/v) acetonitrile:water with 10 mM ammonium acetate, and solvent B, 90:10 (v/v) isopropanol:acetonitrile with 10 mM ammonium acetate. MS analysis was carried out in "auto MS/MS mode" (8MS/s and 6 MS/MS/s) in triplicate analysis with an Agilent 6545 quadrupole TOF mass spectrometer in negative ion mode.

RESULTS AND DISCUSSION
The goal of this study was to address variability in MALDI-MSI experiments while incorporating technical and biological measurements to obtain strong biomarker candidates in the neurodegenerative disease NPC1. To carry this out, we used the Npc1 null mouse model, which is characterized by accumulation of unesterified cholesterol and glycosphingolipids in the late endosomal/lysosomal systems as well as displaying progressive neuron loss prominent in cerebellar Purkinje neurons. The mechanism detailing the interplay of cholesterol and sphingolipid accumulation and neuronal death is still under investigation. Therefore, this study demonstrates a useful tool for all MSI studies, while providing additional information about NPC1 disease specifically.
The general experimental workflow is similar to a traditional MSI, except that this study included technical and biological replicates with these replicates being used to develop the data analysis pipeline. The strategy here addresses the challenges in obtaining spatially mapped analytes in multiple datasets while employing statistically robust methods to evaluate such data. To accomplish this, we performed MSI of lipid species in cerebellar tissue from the Npc1-null mouse model. This model is ideal for such a development because known accumulating lipid species, such as gangliosides, may serve as internal controls while simultaneously allowing the discovery of new potential lipid biomarkers.
Four experiments were carried out to generate molecular maps of lipids in NPC1. First, positive-ion mode was used to detect basic lipids, such as phosphatidylcholines. Second, negative-ion mode was used to obtain data for acidic lipids such as sulfatides, phosphatidylinositols, phosphatidylglycerols, and phosphatidylethanolamines. Third, negativemode data was collected at a higher mass-to-charge ratio range to obtain data for gangliosides, and, finally, we evaluated cholesterol. As observed in Fig. 1 (Top), gangliosides accumulate in the mutant animal tissues analyzed compared with controls and regional specificity was observed. Interestingly, we discern primary accumulation of several gangliosides in lobule X of the cerebellum for the Npc1 mutant animal tissue. It is important to note that images are shown as a total count scale for the specific ion of interest.
An added advantage of using the MSI technique over traditional staining methodologies, such as antibodies, is the specificity achieved for detection of different fatty acid compositions of the same molecular family. Here, we are able to differentiate and simultaneously detect multiple chain length forms of GM1, GM2, and GM3 (Fig. 1). Based on accurate mass of each molecule, they most likely contain stearic or arachidic acids as their fatty acyls. The spatial distributions of the possible GM2 and GM3 stearic forms are different. GM2 (stearic and arachidic forms) reside mainly in regions that appear to be the granular layer of the lobule, whereas GM3 (stearic acid) is found in what we propose is the Purkinje layer. This supports previous findings that GM2 and GM3 are not sequestered together and may have differing roles involving different biochemical mechanisms (28,53,54). Interestingly, the two GM3 species (steric and arachidonic acid) appear to localize differently within lobule X. Although gangliosides have been shown to be elevated in NPC1 disease, these new data indicate the specific distribution at a symptomatic time point in the neurodegeneration progression. Interestingly, lobule X typically does not show loss of Purkinje neurons at this time point in the disease progression in this mouse model; therefore, to observe this strong accumulation is of significant interest. As opposed to high spatial resolution microscopy studies that provide single cell or subcellular imaging, this approach addresses a means to look across the entire tissue for trends. Furthermore, the molecular specificity obtained regarding fatty acid composition is new given that antibody-based staining has only reported total GM1, GM2, or GM3 for NPC1 (28). To provide insight into biological and technical variability, we evaluated images from three serial sections obtained from the same animal. Figure 2 displays several ganglioside replicate images in both mutant and control animals. It should be noted that these images were obtained by manually investigating the lipids of interest, using MS spectra and then generating maps for each genotype and replicate. Therefore, initial assignments are made by accurate mass measurements, followed by ontissue tandem MS experiments. The exact location of double bonds and length of the chains is inferred from database notation and cannot be precisely determined via this method. As a secondary verification, we then evaluated relative differences in ganglioside abundance between control and mutant animals using LC/MS/MS in which elevated amounts were confirmed in cerebellar tissue (supplemental Fig. S2). That is to say that increased levels of GM1, GM2, and GM3 were observed in NPC1 mutant cerebellar via both methods. Fig. 1. Representative MSI images of gangliosides in the cerebellum acquired from Npc1 / (upper) and Npc1 +/+ (lower) mice at 7 weeks of age with corresponding H&E-stained sections (left). Cerebellar lobules noted I and X. Raster distance between pixels is set at 50 m. Scale bar indicates maximum intensity for each ion observed in both the control and mutant images. In all cases, the intensity of the ganglioside was higher in the mutant. The negative ion mass to charge for each ganglioside species is included. Fig. 2. MALDI-MS images of major gangliosides found to be increased in 7-week-old Npc1 / (A) and Npc1 +/+ (B) cerebellum with serial sections for each mouse at a pitch of 50 m. Substantial accumulation is detected for GM3 (d36:1) (m/z 1,179.8) in the Purkinje layer, specifically in lobule X, whereas the accumulation of GM2 (d36:1) (m/z 1,382.8) and GM2 (d38:1) (m/z 1,410.8) is determined to be localized in lobule X. GM1 (m/z 1,544.9) does not reveal localization. On the other hand, GM1 (m/z 1,572.9) does not appear to be accumulating or localizing in the KO mouse compared with the WT. Triplicate analysis shows reproducible ion intensities and lipid distributions, and images were normalized by the TIC of each section using MSiReader.

2450
Journal of Lipid Research Volume 59, 2018

Unbiased detection of cerebellar lipid changes
A challenge with analyzing MSI datasets is having the ability to analyze multiple replicates and identifying analytes with differing abundance or spatial distribution in an unbiased and timely manner. To address this challenge, we developed an algorithmic-based method for clustering replicates, peak selection, and normalization. The following scheme depicts the algorithm executed: Centroided spectra → Peak selection → Normalize → Assess Similarity → Generate Consensus → Peak match to image. The three steps "Normalize," "Assess Similarity," and "Generate Consensus" are embedded in the program used for quality threshold (QT) clustering described below and in previous work evaluating MALDI linear and tandem MS spectra (55,56). For the purposes of this work, however, we have developed the peak selection step and modified the basis on which consensus spectra are compared.
MALDI-MSI spectra tend to have great variation in both absolute intensity and in the baseline noise level due to variations in matrix application and amounts of material present in any given pixel (57). Although every attempt is in place to minimize these effects, recognizing the possibility for heterogeneous sample preparation is crucial. In order to analyze such spectra, we devised an approach based on an analysis of the frequency distribution of peak intensities in a spectrum (Fig. 3). That is, plotting peak intensities independent of the mass-to-charge-ratio at which they occur yields an intensity frequency distribution that can be processed subsequently to yield the ions of most significance in a spectrum. Manual inspection of a number of spectra suggested that incorporating the 40 most intense peaks into a "filtered spectrum" yielded the ions of most interest in that spectrum. Fig. 3A, B are examples of high and low noise baseline in a dataset, respectively. The intensity frequency distributions of these two spectra are shown below each of them (Fig. 3C, D), where the vertical red lines indicate the cut-off for the 40 most intense peaks, i.e., those peaks appearing to the right of the vertical red line. Note that for the spectrum in Fig. 3A, the origin of the intensity frequency distribution is 30 (shown in Fig. 3C), whereas the intensity frequency distribution for the spectrum in Fig. 3B has an origin of one. The number of peaks to be accepted, i.e., the position of the vertical red line, could be modified as required for particular analyses.
The QT clustering algorithm normalizes spectra to unit vectors and then forms clusters of peaks within the multiple spectra that depend on a sensitivity threshold related to the mass measurement precision. Following the clustering, a consensus spectrum for the dataset is formed, and each replicate is compared with the consensus using a dot product analysis. The dot products of each replicate to the consensus are used as a similarity metric; statistical confidence intervals can be calculated and used to evaluate the quality of each replicate spectrum with the aim of including only the most robust replicates into the consensus (55). Outlier spectra can be removed, and the algorithm can be rerun until only the most robust replicates are used to form the consensus spectrum. The overall algorithm workflow in the context of the present study is represented in Fig. 4. Multiple animals were studied for both genotypes (n = 5 and 6 for control and mutant, respectively) and technical replicates (n = 3 or 4 serial sections) for each animal were collected. Each box represents the number of total spectra collected for the MALDI-MSI ion mode of interest. Upon correlation analysis by dot product, and considering the 95% confidence interval, the final column of boxes depicts the number of replicates that passed and were used to generate the final consensus spectrum for the dataset. The resulting consensus spectra for each dataset are a normalized spectrum that represents all replicates for a given condition. Figure 5A-C shows the consensus spectra for each of the three MS datasets collected with a focus on lipids in control and Npc1 mutant animals. Comparison of the spectra in Fig. 5A, showing the ganglioside mass region, we observed differences in normalized intensity for m/z 1,383 and 1,411 corresponding to the two GM2 species previously discussed (see Fig. 1). Tandem MS spectra for each are provided in the supplemental data (supplemental Figs. S3 and S4), where the presence of a sialic acid moiety (m/z 290) Fig. 3. Example of selecting ions from spectra for each sample for the newly developed algorithm. Centroid mass spectrum of a typical noisy spectrum (A) and clean spectrum (B). Frequency distribution analysis of intensities for the corresponding spectra collected (C, D) to illustrate the approach to remove noise. The red lines show implementation of an algorithm to select the 40 most intense peaks from each spectrum. Note that the use of 40 peaks was determined for these dataset and can be changed for any given analysis. AU, arbitrary units. confirms the identity as ganglioside. Unfortunately, the exact lipid chain length and degree of saturation cannot be determined by the on-tissue tandem MS spectra. Additionally, isobaric interference in GM2(d38:1) (supplemental Fig.  S4) shows a fragment ion at m/z 673.5. A common limitation in most MALDI-MSI experiments is the lack of separation of isobaric or near-isobaric species prior to detection as well as issues owing to instrumental limitations mainly in the timed ion selector. Interestingly, the GM1 species are not drastically elevated by this metric, thus demonstrating the power of the normalization strategy and unbiased detection tool, although manual inspection suggests a slight increase in the NPC1 mutant. In fact, the ratios of normalized peak intensity for GM1(d36:1) and GM(38:1) are 1.05 and 0.98, respectively, for the mutant over the control. Compared with LC/MS/MS data, these ratios are slightly decreased, thereby suggesting a conservative approach to identifying differences. To supplement the lipid assignments initially based upon accurate mass measurement, on tissue tandem MS measurements were collected for GM1(d36:1) at m/z 1,544.8 and GM1(d38:1) at m/z 1,572.8 (supplemental Figs. S5 and S6, respectively). The GM3 lipids assigned also included accurate mass measurement and tandem MS spectra (supplemental Figs. S7 and S8). Fig. 5B, C display the consensus spectra for the low mass negative ion mode and the positive ion mode in which significant differences were not observed. It is important to note at this point that the positive ion mode spectra were subjected to an additional step of data processing to remove all ions less than m/z 705, as these mass-to-chargeratio values represent "chemical noise" that we associated with the matrix background. A summary of the consensus results is provided in Fig. 5D for all MS modes collected.
Mean values for each of the consensus spectra and the dot product analysis between the two genotypes is included. Notably, the dot product analysis of control versus Npc1 mutant cerebellar imaging datasets indicated differences with a value of 0.739 for the negative-ion-mode ganglioside region; however, notable spectral differences were not observed in the negative ion mode (0.986) and positive ion mode (0.994) by the dot product analysis. Although this tool drives the further investigation of differences, manual interrogation is still needed to map and evaluate each image.

Total unesterified cholesterol mapping in the cerebellum
Given that NPC is a cholesterol storage disease, we utilized MALDI-MSI to map unesterified cholesterol in the cerebellum of control and Npc1 mutant mice (Fig. 6). In control animals, a slightly more abundant localization of cholesterol is observed in the white matter region, whereas this is less pronounced in the Npc1 mutants (Fig. 6A). This observation is consistent with reduced myelin in NPC1 (58,59). Owing to the spatial resolution limits of our instrument, we are unable to localize intra-lysosomal storage in the Npc1 mutant animals. Using our consensus algorithm (Fig. 6B), we evaluated the normalized levels of cholesterol across the cerebellum. The consensus spectra reveal similar levels of cholesterol from biological replicates (n = 3 each), which included technical analyses (n = 7 for Npc1 +/+ and n = 8 Npc1 / ). The indistinguishable difference in cholesterol levels agrees with previous work in which liver, kidney, spleen, and brain tissue from the Npc1 null mouse model, which were analyzed by GC/MS for cholesterol and oxysterol levels. Compared with controls, the brain did not have elevated levels of cholesterol despite altered levels in other tissues (60).

Spatially resolved lipids in lobules
The age of the animals in the present study represents a symptomatic point in the NPC1 disease progression. With that knowledge, and based upon our preliminary data, we further analyzed the dataset to investigate changes we observed in only lobule X of the cerebellum. The original analysis considered data across the entire cerebellum, where at the 7-week time point, significant cerebellar degeneration has occurred (19,61). Therefore, when observing differences by manual mapping, spatial, and quantitative differences for the GM3 species, as an example, were not found in the consensus spectra comparison. We attribute this to the possibility that in the earlier lobules, Purkinje neurons are known to degenerate, and, therefore, the accumulation in the later lobules is less pronounced by this method. To address this observation, further refinement of the approach was carried out. First, rather than analyzing the entire cerebellar tissue, only lobule X was analyzed ( Fig. 7A-C). This is appropriate, given that this region of the cerebellum is preserved in the late stages of NPC1 disease and may provide insight into the neurodegeneration mechanism of NPC1 (19). Second, the analysis included the 60 most intense ions, which was determined after manual evaluation of the spectra. Similar differences were Fig. 4. General workflow for collection and processing of multiple biological and technical replicates in MALDI-MSI datasets. Biological replicates were processed by investigating n = 6 (control) and n = 5 (mutant) animals. Cerebellar tissue was imaged using MALDI-MSI, and spectra were collected in three experiments to detect GM, negative, and positive lipid ions. Using technical replicates from each animal, as indicated by "n," the newly developed alignment, normalization, and consensus spectra algorithm was utilized. To evaluate spectral reproducibility, a dot product analysis was carried out and evaluated based on correlation. Outlier spectra were removed, and those that "passed" are indicated by having a correlation above the 95% CI.

2452
Journal of Lipid Research Volume 59, 2018 observed in the consensus spectra for lobule X, as was observed in the whole cerebellar tissue analysis.
Based upon our imaging and lobule X analysis, we sought to determine whether the differences in gangliosides were only present in lobule X. We then evaluated the ganglioside mass range in cerebellar lobule IV/V to evaluate differences compared with lobule X in control and NPC1 mutant animals. This region of the cerebellum was affected Consensus spectra and summary of MALDI-MSI replicate analysis in NPC1. Consensus mass spectra are provided for the GM mass region (A), negative ion (B), and positive ion (C), respectively. Data for WT is indicated in blue (below) and for the mutant in red (above). D: Data table shows average correlation coefficient between replicates for each dataset and the mean value of the 95% CI. The dot product analysis of the consensus spectra is also provided between the control and disease datasets (far right column). Dot products >0.9 are indicative of high similarity and, therefore, do not suggest spectral differences between control and mutant datasets. +/+ (right) mice at 7 weeks of age. Spatially, MSI images reveal an overall similar distribution between the genotypes however cholesterol may be slightly higher in the white matter of control animals. Scale bar represents 0-9,000 counts. B: Consensus spectra alignment following normalization of the top most intense 20 ions is provided for the cholesterol mass range for the whole cerebellum. The consensus spectra reveals similar normalized abundance of cholesterol for mutant (top, red) and WT (bottom, blue) tissue. Three animals of each genotype were analyzed, which included n = 8 and 7 technical replicates for mutant and WT, respectively. earlier in the disease in terms of Purkinje neuron loss. Because the abundance differences in lobule X were so pronounced, we sought to evaluate a region of the cerebellum that may have differences, yet those differences may not be quite as obvious. Fig. 7D displays the consensus spectrum for the ganglioside mass range of lobule IV/V in control and Npc1 mutant cerebella. Ions corresponding to examples of ganglioside molecules GM1, GM2, and GM3 are denoted with stars. Consistent with our analysis in both the complete cerebellum and lobule X alone, elevation of gangliosides was observed in the mutant animal compared with control. From these data, we can conclude that the spectrum-based QT clustering normalization tool provides unbiased evaluating of biological and technical replicate imaging experiments to reveal alterations in disease models.

CONCLUSIONS
Using MALDI-MSI, lipid mapping, and visualization of different gangliosides species including fatty acidspecific composition of GM3, GM2, and GM1 in the cerebellum of a 7-week-old Npc1 / mouse model was performed. Using this label-free technique, we were able to detect variants of different fatty acid chain length and determine the localizations of gangliosides. The use of unbiased peak selection and QT clustering for calculation of consensus spectra was able to validate spectral-expected features expressed differentially between the mutant and control datasets. Proofof-concept experiments included MALDI-MSI analysis of cholesterol in the cerebellum. Validation by MS/MS and LC/MS analysis confirmed alterations spatially across the entire cerebellum and within specific lobule regions. With these tools set in place, continued studies on can be performed to evaluate the lipid changes during the entire disease progression of NPC1. Fig. 7. Lobule analysis in NPC1. Consensus spectra alignment following normalization of the top most intense 60 ions is provided for the GM mass range (A), negative ion mode (B), and positive ion mode (C) with a focus on lobule X, which does not show Purkinje cell loss in NPC1 disease. Significant spectral intensity differences are observed in the ganglioside region. D: Consensus spectrum obtained from evaluating the ganglioside mass range for lobules IV/V in control and Npc1 mutant cerebella. Example gangliosides presented in Fig. 1 are denoted with stars including m/z 1,179.8; 1,382.9; and 1,573.0 representing variants of GM3, GM2, and GM1, respectively. Similar to whole cerebellar analysis, these sphingolipids are altered in Npc1 mutant lobule IV/V compared with control.