Glomeruli filter blood through the coordination of podocytes, mesangial cells, fenestrated endothelial cells, and the glomerular basement membrane. Cellular changes, such as podocyte loss, are associated with pathologies like diabetic kidney disease. However, little is known regarding the in situ molecular profiles of specific cell types and how these profiles change with disease. Matrix-assisted laser desorption/ionization imaging mass spectrometry (MALDI IMS) is well-suited for untargeted tissue mapping of a wide range of molecular classes. Importantly, additional imaging modalities can be integrated with MALDI IMS to associate these biomolecular distributions to specific cell types. Here, we integrated workflow combining MALDI IMS and multiplexed immunofluorescence (MxIF) microscopy. High spatial resolution MALDI IMS (5 μm) was used to determine lipid distributions within human glomeruli from a normal portion of fresh-frozen kidney cancer nephrectomy tissue revealing intra-glomerular lipid heterogeneity. Mass spectrometric data were linked to specific glomerular cell types and substructures through new methods that enable MxIF microscopy to be performed on the same tissue section following MALDI IMS, without sacrificing signal quality from either modality. Machine learning approaches were combined enabling cell type segmentation and identification based on MxIF data. This was followed by mining of cell type or cluster-associated MALDI IMS signatures using classification and interpretable machine learning. This allowed automated discovery of spatially specific molecular markers for glomerular cell types and substructures as well as lipids correlated to deep and superficial glomeruli. Overall, our work establishes a toolbox for probing molecular signatures of glomerular cell types and substructures within tissue microenvironments providing a framework applicable to other kidney tissue features and organ systems.

Introduction: Age related macular degeneration (AMD) causes legal blindness worldwide, with few therapeutic targets in early disease and no treatments for 80% of cases. Extracellular deposits, including drusen and subretinal drusenoid deposits (SDD; also called reticular pseudodrusen), disrupt cone and rod photoreceptor functions and strongly confer risk for advanced disease. Due to the differential cholesterol composition of drusen and SDD, lipid transfer and cycling between photoreceptors and support cells are candidate dysregulated pathways leading to deposit formation. The current study explores this hypothesis through a comprehensive lipid compositional analysis of SDD. Methods: Histology and transmission electron microscopy were used to characterize the morphology of SDD. Highly sensitive tools of imaging mass spectrometry (IMS) and nano liquid chromatography tandem mass spectrometry (nLC-MS/MS) in positive and negative ion modes were used to spatially map and identify SDD lipids, respectively. An interpretable supervised machine learning approach was utilized to compare the lipid composition of SDD to regions of uninvolved retina across 1873 IMS features and to automatically discern candidate markers for SDD. Immunohistochemistry (IHC) was used to localize secretory phospholipase A2 group 5 (PLA2G5). Results: Among the 1873 detected features in IMS data, three lipid classes, including lysophosphatidylcholine (LysoPC), lysophosphatidylethanolamine (LysoPE) and lysophosphatidic acid (LysoPA) were observed nearly exclusively in SDD while presumed precursors, including phosphatidylcholine (PC), phosphatidylethanolamine (PE) and phosphatidic acid (PA) lipids were detected in SDD and adjacent photoreceptor outer segments. Molecular signals specific to SDD were found in central retina and elsewhere. IHC results indicated abundant PLA2G5 in photoreceptors and retinal pigment epithelium (RPE). Discussion: The abundance of lysolipids in SDD implicates lipid remodeling or degradation in deposit formation, consistent with ultrastructural evidence of electron dense lipid-containing structures distinct from photoreceptor outer segment disks and immunolocalization of secretory PLA2G5 in photoreceptors and RPE. Further studies are required to understand the role of lipid signals observed in and around SDD.

Imaging mass spectrometry (IMS) provides untargeted, highly multiplexed maps of molecular distributions in tissue. Ion images are routinely presented as heatmaps and can be overlaid onto complementary microscopy images that provide greater context. However, heatmaps use transparency blending to visualize both images, obscuring subtle quantitative differences and distribution gradients. Here, we developed a contour mapping approach that combines information from IMS ion intensity distributions with that of stained microscopy. As a case study, we applied this approach to imaging data from Staphylococcus aureus-infected murine kidney. In a univariate, or single molecular species, use-case of the contour map representation of IMS data, certain lipids colocalizing with regions of infection were selected using Pearson’s correlation coefficient. Contour maps of these lipids overlaid with stained microscopy showed enhanced visualization of lipid distributions and spatial gradients in and around the bacterial abscess as compared to traditional heatmaps. The full IMS data set comprising hundreds of individual ion images was then grouped into a smaller subset of representative patterns using non-negative matrix factorization (NMF). Contour maps of these multivariate NMF images revealed distinct molecular profiles of the major abscesses and surrounding immune response. This contour mapping workflow also enabled a molecular visualization of the transition zone at the host-pathogen interface, providing potential clues about the spatial molecular dynamics beyond what histological staining alone provides. In summary, we developed a new IMS-based contour mapping approach to augment classical stained microscopy images, providing an enhanced and more interpretable visualization of IMS-microscopy multimodal molecular imaging data sets.

Spatially targeted proteomics analyzes the proteome of specific cell types and functional regions within tissue. While spatial context is often essential to understanding biological processes, interpreting sub-region-specific protein profiles can pose a challenge due to the high-dimensional nature of the data. Here, we develop a multivariate approach for rapid exploration of differential protein profiles acquired from distinct tissue regions and apply it to analyze a published spatially targeted proteomics data set collected from Staphylococcus aureus-infected murine kidney, 4 and 10 days postinfection. The data analysis process rapidly filters high-dimensional proteomic data to reveal relevant differentiating species among hundreds to thousands of measured molecules. We employ principal component analysis (PCA) for dimensionality reduction of protein profiles measured by microliquid extraction surface analysis mass spectrometry. Subsequently, k-means clustering of the PCA-processed data groups samples by chemical similarity. Cluster center interpretation revealed a subset of proteins that differentiate between spatial regions of infection over two time points. These proteins appear involved in tricarboxylic acid metabolomic pathways, calcium-dependent processes, and cytoskeletal organization. Gene ontology analysis further uncovered relationships to tissue damage/repair and calcium-related defense mechanisms. Applying our analysis in infectious disease highlighted differential proteomic changes across abscess regions over time, reflecting the dynamic nature of host-pathogen interactions.

The glomerulus is a multicellular functional tissue unit (FTU) of the nephron that is responsible for blood filtration. Each glomerulus contains multiple substructures and cell types that are crucial for their function. To understand normal aging and disease in kidneys, methods for high spatial resolution molecular imaging within these FTUs across whole slide images is required. Here we demonstrate a workflow using microscopy-driven selected sampling to enable 5 μm pixel size matrix-assisted laser desorption/ionization imaging mass spectrometry (MALDI IMS) of all glomeruli within whole slide human kidney tissues. Such high spatial resolution imaging entails large numbers of pixels, increasing the data acquisition times. Automating FTU-specific tissue sampling enables high-resolution analysis of critical tissue structures, while concurrently maintaining throughput. Glomeruli were automatically segmented using coregistered autofluorescence microscopy data, and these segmentations were translated into MALDI IMS measurement regions. This allowed high-throughput acquisition of 268 glomeruli from a single whole slide human kidney tissue section. Unsupervised machine learning methods were used to discover molecular profiles of glomerular subregions and differentiate between healthy and diseased glomeruli. Average spectra for each glomerulus were analyzed using Uniform Manifold Approximation and Projection (UMAP) and k-means clustering, yielding 7 distinct groups of differentiated healthy and diseased glomeruli. Pixel-wise k-means clustering was applied to all glomeruli, showing unique molecular profiles localized to subregions within each glomerulus. Automated microscopy-driven, FTU-targeted acquisition for high spatial resolution molecular imaging maintains high-throughput and enables rapid assessment of whole slide images at cellular resolution and identification of tissue features associated with normal aging and disease.

Here, we describe the preservation and preparation of human kidney tissue for interrogation by histopathology, imaging mass spectrometry, and multiplexed immunofluorescence. Custom image registration and integration techniques are used to create cellular and molecular atlases of this organ system. Through careful optimization, we ensure high-quality and reproducible datasets suitable for cross-patient comparisons that are essential to understanding human health and disease. Moreover, each of these steps can be adapted to other organ systems or diseases, enabling additional atlas efforts.

The search for molecular species that are differentially expressed between biological states is an important step towards discovering promising biomarker candidates. In imaging mass spectrometry (IMS), performing this search manually is often impractical due to the large size and high-dimensionality of IMS datasets. Instead, we propose an interpretable machine learning workflow that automatically identifies biomarker candidates by their mass-to-charge ratios, and that quantitatively estimates their relevance to recognizing a given biological class using Shapley additive explanations (SHAP). The task of biomarker candidate discovery is translated into a feature ranking problem: given a classification model that assigns pixels to different biological classes on the basis of their mass spectra, the molecular species that the model uses as features are ranked in descending order of relative predictive importance such that the top-ranking features have a higher likelihood of being useful biomarkers. Besides providing the user with an experiment-wide measure of a molecular species' biomarker potential, our workflow delivers spatially localized explanations of the classification model's decision-making process in the form of a novel representation called SHAP maps. SHAP maps deliver insight into the spatial specificity of biomarker candidates by highlighting in which regions of the tissue sample each feature provides discriminative information and in which regions it does not. SHAP maps also enable one to determine whether the relationship between a biomarker candidate and a biological state of interest is correlative or anticorrelative. Our automated approach to estimating a molecular species' potential for characterizing a user-provided biological class, combined with the untargeted and multiplexed nature of IMS, allows for the rapid screening of thousands of molecular species and the obtention of a broader biomarker candidate shortlist than would be possible through targeted manual assessment. Our biomarker candidate discovery workflow is demonstrated on mouse-pup and rat kidney case studies.

Identifying the spatial distributions of biomolecules in tissue is crucial for understanding integrated function. Imaging mass spectrometry (IMS) allows simultaneous mapping of thousands of biosynthetic products such as lipids but has needed a means of identifying specific cell-types or functional states to correlate with molecular localization. We report, here, advances starting from identity marking with a genetically encoded fluorophore. The fluorescence emission data were integrated with IMS data through multimodal image processing with advanced registration techniques and data-driven image fusion. In an unbiased analysis of spleens, this integrated technology enabled identification of ether lipid species preferentially enriched in germinal centers. We propose that this use of genetic marking for microanatomical regions of interest can be paired with molecular information from IMS for any tissue, cell-type, or activity state for which fluorescence is driven by a gene-tracking allele and ultimately with outputs of other means of spatial mapping.

Histology-directed imaging mass spectrometry (IMS) is a spatially targeted IMS acquisition method informed by expert annotation that provides rapid molecular characterization of select tissue structures. The expert annotations are usually determined on digital whole slide images of histological stains where the staining preparation is incompatible with optimal IMS preparation, necessitating serial sections: one for annotation, one for IMS. Registration is then used to align staining annotations onto the IMS tissue section. Herein, we report a next-generation histology-directed platform implementing IMS-compatible autofluorescence (AF) microscopy taken prior to any staining or IMS. The platform enables two histology-directed workflows, one that improves the registration process between two separate tissue sections using automated, computational monomodal AF-to-AF microscopy image registration, and a registration-free approach that utilizes AF directly to identify ROIs and acquire IMS on the same section. The registration approach is fully automated and delivers state of the art accuracy in histology-directed workflows for transfer of annotations (∼3-10 μm based on 4 organs from 2 species) while the direct AF approach is registration-free, allowing targeting of the finest structures visible by AF microscopy. We demonstrate the platform in biologically relevant case studies of liver stage malaria and human kidney disease with spatially targeted acquisition of sparsely distributed (composing less than one tenth of 1% of the tissue section area) malaria infected mouse hepatocytes and glomeruli in the human kidney case study.

The correlation of imaging mass spectrometry (IMS) with histopathology can help relate novel molecular findings obtained through IMS to the well-characterized and validated histopathology knowledge base. The quality of correlation between these two modalities is limited by the quality of the spatial mapping that is obtained by registration of the two image types. In this work, we develop novel workflows for MALDI IMS-to-microscopy data registration and analysis using nondestructive IMS-compatible wide field autofluorescence (AF) microscopy combined with computational image registration. First, a substantially automated procedure for high-accuracy registration between IMS and microscopy data of the same section is described that explicitly links the MALDI laser ablation pattern imaged by microscopy to its corresponding IMS pixel. Subsequent examination of the registered data allows for high-confidence colocalization of image features between the two modalities, down to single-cell scales within tissue. Building on this IMS-microscopy spatial mapping, we furthermore demonstrate the automated spatial correlation between IMS measurements from serial sections. This AF-registration-driven inter-section analysis, using a combination of nonlinear AF-to-AF and IMS-to-AF image registrations, can be applied to tissue sections that are prepared and imaged with different sample preparations (e.g., lipids vs proteins) and/or that are measured using different spatial resolutions. Importantly, all registrations, whether within a single section or across serial sections, are entirely independent of the IMS intensity signal content and thus unbiased by it.