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.

@en

Multimodal imaging by matrix-assisted laser desorption ionisation mass spectrometry imaging (MALDI MSI) and microscopy holds potential for understanding pathological mechanisms by mapping molecular signatures from the tissue microenvironment to specific cell populations. However, existing software solutions for MALDI MSI data analysis are incomplete, require programming skills and contain laborious manual steps, hindering broadly applicable, reproducible, and high-throughput analysis to generate impactful biological discoveries. Here, we present msiFlow, an accessible open-source, platform-independent and vendor-neutral software for end-to-end, high-throughput, transparent and reproducible analysis of multimodal imaging data. msiFlow integrates all necessary steps from raw data import to analytical visualisation along with state-of-the-art and self-developed algorithms into automated workflows. Using msiFlow, we unravel the molecular heterogeneity of leukocytes in infected tissues by spatial regulation of ether-linked phospholipids containing arachidonic acid. We anticipate that msiFlow will facilitate the broad applicability of MSI in multimodal imaging to uncover context-dependent cellular regulations in disease states.

@en

While probabilistic constellation shaping (PCS) enables rate and reach adaption with finer granularity [1] (Cho and Winzer, 2009), it imposes signal processing challenges at the receiver. Since the distribution of PCS-quadrature amplitude modulation (QAM) signals tends to be Gaussian, conventional blind polarization demultiplexing algorithms are not suitable for them [2] (Johnson et al., 1998). It is known that independently and identically distributed (iid) Gaussian signals, when mixed, cannot be recovered/separated from their mixture. For PCS-QAM signals, there are algorithms such as [3] and [4] Dris et al. (2019) and Athuraliya et al. (2004) which are designed by extending conventional blind algorithms used for uniform QAM signals. In these algorithms, an initialization point is obtained by processing only a part of the mixed signal, which have non-Gaussian statistics. In this article, we propose an alternative method wherein we add temporal correlations at the transmitter, which are subsequently exploited at the receiver in order to separate the polarizations. We will refer to the proposed method as frequency domain (FD) joint diagonalization (JD) probability aware-multi modulus algorithm (pr-MMA), and it is suited to channels with moderate polarization mode dispersion (PMD) effects. Furthermore, we extend our previously proposed JD-MMA [5] (Bajaj et al., 2022) by replacing the standard MMA with a pr-MMA, improving its performance. Both FDJD-pr-MMA and JD-pr-MMA are evaluated for a diverse range of PCS (entropy $\mathcal {H}$) of 64-QAM over a first-order PMD channel that is simulated in a proof-of-concept setup. A MMA initialized with a memoryless constant modulus algorithm (CMA) is used as a benchmark. We show that at a differential group delay (DGD) of 10% of symbol period T$_{\text{symb}}$ and 18 dB SNR/pol., JD-pr-MMA successfully demultiplexes the PCS signals, while CMA-MMA fails drastically. Furthermore, we demonstrate that the newly proposed FDJD-pr-MMA is robust against moderate PMD effects by evaluating it over a DGD of up to 40% of T$_{\text{symb}}$. Our results show that the proposed FDJD-pr-MMA successfully equalizes PMD channels with a DGD up to 20% of T$_{\text{symb}}$.

@en

Correction to: Nature Cell Biologyhttps://doi.org/10.1038/s41556-023-01194-w. Published online 19 July 2023. In the version of this article originally published, the name of Tianyang Xu was misspelled as Tiangyang Xu. The name has been corrected in the HTML and PDF versions of the article.

@en

The lack of standardization in antibody validation remains a major contributor to irreproducibility of human research. To address this, we have applied a standardized approach to validate a panel of antibodies to identify 18 major cell types and 5 extracellular matrix compartments in the human kidney by immunofluorescence (IF) microscopy. We have used these to generate an organ mapping antibody panel for two-dimensional (2-D) and three-dimensional (3-D) cyclical IF (CyCIF) to provide a more detailed method for evaluating tissue segmentation and volumes using a larger panel of markers than would normally be possible using standard fluorescence microscopy. CyCIF also makes it possible to perform multiplexed IF microscopy of whole slide images, which is a distinct advantage over other multiplexed imaging technologies that are applicable to limited fields of view. This enables a broader view of cell distributions across larger anatomical regions, allowing a better chance to capture localized regions of dysfunction in diseased tissues. These methods are broadly accessible to any laboratory with a fluorescence microscope, enabling spatial cellular phenotyping in normal and disease states. We also provide a detailed solution for image alignment between CyCIF cycles that can be used by investigators to perform these studies without programming experience using open-sourced software. This ability to perform multiplexed imaging without specialized instrumentation or computational skills opens the door to integration with more highly dimensional molecular imaging modalities such as spatial transcriptomics and imaging mass spectrometry, enabling the discovery of molecular markers of specific cell types, and how these are altered in disease.NEW & NOTEWORTHY We describe here validation criteria used to define on organ mapping panel of antibodies that can be used to define 18 cell types and five extracellular matrix compartments using cyclical immunofluorescence (CyCIF) microscopy. As CyCIF does not require specialized instrumentation, and image registration required to assemble CyCIF images can be performed by any laboratory without specialized computational skills, this technology is accessible to any laboratory with access to a fluorescence microscope and digital scanner.

@en

There is a need to define regions of gene activation or repression that control human kidney cells in states of health, injury, and repair to understand the molecular pathogenesis of kidney disease and design therapeutic strategies. Comprehensive integration of gene expression with epigenetic features that define regulatory elements remains a significant challenge. We measure dual single nucleus RNA expression and chromatin accessibility, DNA methylation, and H3K27ac, H3K4me1, H3K4me3, and H3K27me3 histone modifications to decipher the chromatin landscape and gene regulation of the kidney in reference and adaptive injury states. We establish a spatially-anchored epigenomic atlas to define the kidney’s active, silent, and regulatory accessible chromatin regions across the genome. Using this atlas, we note distinct control of adaptive injury in different epithelial cell types. A proximal tubule cell transcription factor network of ELF3, KLF6, and KLF10 regulates the transition between health and injury, while in thick ascending limb cells this transition is regulated by NR2F1. Further, combined perturbation of ELF3, KLF6, and KLF10 distinguishes two adaptive proximal tubular cell subtypes, one of which manifested a repair trajectory after knockout. This atlas will serve as a foundation to facilitate targeted cell-specific therapeutics by reprogramming gene regulatory networks.@en

Purpose: Mutations in the gene that encodes the enzyme acid sphingomyelinase (ASMase) are associated with Niemann-Pick disease, a lysosomal storage disorder. Mice that lack ASMase (ASMase-/-) exhibit age-related retinal degeneration and large increases in accumulation of lipofuscin in the retinal pigment epithelium (RPE). We examined which lipid species accumulate in the retina and the RPE of ASMase-/- mice and whether the retinal degeneration is associated with impaired photoreceptor metabolism and retinyl chromophore processing. Methods: NADPH availability and all-trans retinol formation after rhodopsin bleaching were measured in isolated single rod photoreceptors with fluorescence imaging; sphingolipid levels in retinas and RPEs were measured with LC/MS; relative abundances of different lipid species in different retinal layers were measured with MALDI imaging mass spectrometry. Results: There was no detectable difference in the kinetics of all-trans retinol formation or the NADPH-generating capacity between ASMase-/- and wild-type mice. Sphingomyelin levels were much higher in the retinas and RPEs of ASMase-/- animals compared to wild type, but there were no significant differences for ceramides. There was a large increase in the abundance of bis(monoacylglycero)phosphates (BMPs) in ASMase-/- mice, indicative of lysosomal dysfunction, but no substantial changes were detected for the bis-retinoid A2E. Conclusions: Lysosomal dysfunction and retinal degeneration in ASMase-/- mice are not associated with defects in rod photoreceptor metabolism that affect all-trans retinol formation and availability of NADPH. Lysosomal dysfunction in ASMase-/- mice is not associated with bis-retinoid A2E accumulation.

@en

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.

@en

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.

@en

Macroscopic x-ray fluorescence imaging spectroscopy (MA-XRF) and reflectance imaging spectroscopy (RIS) are important tools in the analysis of cultural heritage objects, both for conservation and art historical research purposes. The elemental and molecular distributions provided by MA-XRF and RIS respectively, are particularly useful for the identification and mapping of pigments in easel paintings. While MA-XRF has relatively established data processing methods based on modeling of the underlying physics, RIS data cannot be modeled with sufficient precision and its processing has considerable room for improvements. This work seeks to improve RIS data processing workflows in the short wavelength infrared range (SWIR, 1000–2500 nm) with a novel method that fits Gaussian profiles to pigment-specific absorption features, and we compare its performance to MA-XRF for the task of semi-quantitative pigment mapping, evaluating their limits of detection (LODs) and the matrix effects that affect their signals. Two pigments are considered in this work, lead white and blue verditer, which are mapped in SWIR RIS using the first overtone of -OH stretching of their primary compounds, hydrocerussite (Pb₃(CO₃)₂(OH)₂) and azurite (Cu₃(CO₃)₂(OH)₂), at 1447 and 1497 nm respectively, and in MA-XRF using the Pb-L and Cu-K fluorescence signals. The methods are evaluated using two sets of custom-prepared paint samples, as well as a 16th-century painting, discussing the identification, mapping, and semi-quantitative analysis of the considered pigments. We found SWIR RIS to be a pigment-specific method with a longer linear range but inferior LODs and penetration depth when compared to MA-XRF, the latter is often not capable of discriminating between different pigments with identical elemental markers. We furthermore present a novel color scale that allows the simultaneous visualization of signals above and below a confidence limit.

@en

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.

@en

The Human BioMolecular Atlas Program (HuBMAP) aims to create a multi-scale spatial atlas of the healthy human body at single-cell resolution by applying advanced technologies and disseminating resources to the community. As the HuBMAP moves past its first phase, creating ontologies, protocols and pipelines, this Perspective introduces the production phase: the generation of reference spatial maps of functional tissue units across many organs from diverse populations and the creation of mapping tools and infrastructure to advance biomedical research.

@en

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.

@en

Gangliosides are acidic glycosphingolipids, containing ceramide moieties and oligosaccharide chains with one or more sialic acid residue(s) and are highly diverse isomeric structures with distinct biological roles. Matrix-assisted laser desorption/ionization imaging mass spectrometry (MALDI IMS) enables the untargeted spatial analysis of gangliosides, among other biomolecules, directly from tissue sections. Integrating trapped ion mobility spectrometry with MALDI IMS allows for the analysis of isomeric lipid structures in situ. Here, we demonstrate the gas-phase separation and identification of disialoganglioside isomers GD1a and GD1b that differ in the position of a sialic acid residue, in multiple samples, including a standard mixture of both isomers, a biological extract, and directly from thin tissue sections. The unique spatial distributions of GD1a/b (d36:1) and GD1a/b (d38:1) isomers were determined in rat hippocampus and spinal cord tissue sections, demonstrating the ability to structurally characterize and spatially map gangliosides based on both the carbohydrate chain and ceramide moieties.

@en

We propose a novel method for blind polarization-demultiplexing of probabilistically shaped signals for coherent receivers. The method is capable of separating signals with (quasi) Gaussian distributions by exploiting temporal correlations added to the transmit signals. The proposed method is evaluated in challenging mixing scenarios.

@en

Introduction: Although Staphylococcus aureus is the leading cause of biofilm-related infections, the lipidomic distributions within these biofilms is poorly understood. Here, lipidomic mapping of S. aureus biofilm cross-sections was performed to investigate heterogeneity between horizontal biofilm layers. Methods: S. aureus biofilms were grown statically, embedded in a mixture of carboxymethylcellulose/gelatin, and prepared for downstream matrix-assisted laser desorption/ionization imaging mass spectrometry (MALDI IMS). Trapped ion mobility spectrometry (TIMS) was also applied prior to mass analysis. Results: Implementation of TIMS led to a ∼ threefold increase in the number of lipid species detected. Washing biofilm samples with ammonium formate (150 mM) increased signal intensity for some bacterial lipids by as much as tenfold, with minimal disruption of the biofilm structure. MALDI TIMS IMS revealed that most lipids localize primarily to a single biofilm layer, and species from the same lipid class such as cardiolipins CL(57:0) – CL(66:0) display starkly different localizations, exhibiting between 1.5 and 6.3-fold intensity differences between layers (n = 3, p < 0.03). No horizontal layers were observed within biofilms grown anaerobically, and lipids were distributed homogenously. Conclusions: High spatial resolution analysis of S. aureus biofilm cross-sections by MALDI TIMS IMS revealed stark lipidomic heterogeneity between horizontal S. aureus biofilm layers demonstrating that each layer was molecularly distinct. Finally, this workflow uncovered an absence of layers in biofilms grown under anaerobic conditions, possibly indicating that oxygen contributes to the observed heterogeneity under aerobic conditions. Future applications of this workflow to study spatially localized molecular responses to antimicrobials could provide new therapeutic strategies.

@en

Bacterial pathogens have evolved virulence factors to colonize, replicate, and disseminate within the vertebrate host. Although there is an expanding body of literature describing how bacterial pathogens regulate their virulence repertoire in response to environmental signals, it is challenging to directly visualize virulence response within the host tissue microenvironment. Multimodal imaging approaches enable visualization of host-pathogen molecular interactions. Here we demonstrate multimodal integration of high spatial resolution imaging mass spectrometry and microscopy to visualize Staphylococcus aureus envelope modifications within infected murine and human tissues. Data-driven image fusion of fluorescent bacterial reporters and matrix-assisted laser desorption/ionization Fourier transform ion cyclotron resonance imaging mass spectrometry uncovered S. aureus lysyl-phosphatidylglycerol lipids, localizing to select bacterial communities within infected tissue. Absence of lysyl-phosphatidylglycerols is associated with decreased pathogenicity during vertebrate colonization as these lipids provide protection against the innate immune system. The presence of distinct staphylococcal lysyl-phosphatidylglycerol distributions within murine and human infections suggests a heterogeneous, spatially oriented microbial response to host defenses.

@en

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.

@en

We demonstrate the utility of combining silicon nanopost arrays (NAPA) and trapped ion mobility imaging mass spectrometry (TIMS IMS) for high spatial resolution and specificity mapping of neutral lipid classes in tissue. Ionization of neutral lipid species such as triglycerides (TGs), cholestryl esters (CEs), and hexosylceramides (HexCers) from biological tissues has remained a challenge for imaging applications. NAPA, a matrix-free laser desorption ionization substrate, provides enhanced ionization efficiency for the above-mentioned neutral lipid species, providing complementary lipid coverage to matrix-assisted laser desorption ionization (MALDI). The combination of NAPA and TIMS IMS enables imaging of neutral lipid species at 20 μm spatial resolution while also increasing molecular coverage greater than 2-fold using gas-phase ion mobility separations. This is a significant improvement with respect to sensitivity, specificity, and spatial resolution compared to previously reported imaging studies using NAPA alone. Improved specificity for neutral lipid analysis using TIMS IMS was shown using rat kidney tissue to separate TGs, CEs, HexCers, and phospholipids into distinct ion mobility trendlines. Further, this technology allowed for the separation of isomeric species, including mobility resolved isomers of Cer(d42:2) (m/z 686.585) with distinct spatial localizations measured in rat kidney tissue section.

@en

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.

@en