Spatially Resolved Top-Down Proteomics of Tissue Sections Based on a Microfluidic Nanodroplet Sample Preparation Platform

Conventional proteomic approaches measure the averaged signal from mixed cell populations or bulk tissues, leading to the dilution of signals arising from subpopulations of cells that might serve as important biomarkers. Recent developments in bottom-up proteomics have enabled spatial mapping of cellular heterogeneity in tissue microenvironments. However, bottom-up proteomics cannot unambiguously define and quantify proteoforms, which are intact (i.e., functional) forms of proteins capturing genetic variations, alternatively spliced transcripts and posttranslational modifications. Herein, we described a spatially resolved top-down proteomics (TDP) platform for proteoform identification and quantitation directly from tissue sections. The spatial TDP platform consisted of a nanodroplet processing in one pot for trace samples-based sample preparation system and an laser capture microdissection-based cell isolation system. We improved the nanodroplet processing in one pot for trace samples sample preparation by adding benzonase in the extraction buffer to enhance the coverage of nucleus proteins. Using ∼200 cultured cells as test samples, this approach increased total proteoform identifications from 493 to 700; with newly identified proteoforms primarily corresponding to nuclear proteins. To demonstrate the spatial TDP platform in tissue samples, we analyzed laser capture microdissection-isolated tissue voxels from rat brain cortex and hypothalamus regions. We quantified 509 proteoforms within the union of top-down mass spectrometry-based proteoform identification and characterization and TDPortal identifications to match with features from protein mass extractor. Several proteoforms corresponding to the same gene exhibited mixed abundance profiles between two tissue regions, suggesting potential posttranslational modification-specific spatial distributions. The spatial TDP workflow has prospects for biomarker discovery at proteoform level from small tissue sections.

Proteome Mapping of the Human Pancreatic Islet Microenvironment Reveals Endocrine-Exocrine Signaling Sphere of Influence

The need for a clinically accessible method with the ability to match protein activity within heterogeneous tissues is currently unmet by existing technologies. Our proteomics sample preparation platform, named microPOTS (Microdroplet Processing in One pot for Trace Samples), can be used to measure relative protein abundance in micron-scale samples alongside the spatial location of each measurement, thereby tying biologically interesting proteins and pathways to distinct regions. However, given the smaller pixel/voxel number and amount of tissue measured, standard mass spectrometric analysis pipelines have proven inadequate. Here we describe how existing computational approaches can be adapted to focus on the specific biological questions asked in spatial proteomics experiments. We apply this approach to present an unbiased characterization of the human islet microenvironment comprising the entire complex array of cell types involved while maintaining spatial information and the degree of the islet's sphere of influence. We identify specific functional activity unique to the pancreatic islet cells and demonstrate how far their signature can be detected in the adjacent tissue. Our results show that we can distinguish pancreatic islet cells from the neighboring exocrine tissue environment, recapitulate known biological functions of islet cells, and identify a spatial gradient in the expression of RNA processing proteins within the islet microenvironment.

Coupling microdroplet-based sample preparation, multiplexed isobaric labeling, and nanoflow peptide fractionation for deep proteome profiling of tissue microenvironment

There is increasing interest in developing in-depth proteomic approaches for mapping tissue heterogeneity at a cell-type-specific level to better understand and predict the function of complex biological systems, such as human organs. Existing spatially resolved proteomics technologies cannot provide deep proteome coverages due to limited sensitivity and poor sample recovery. Herein, we seamlessly combined laser capture microdissection with a low-volume sample processing technology that includes a microfluidic device named microPOTS (Microdroplet Processing in One pot for Trace Samples), the multiplexed isobaric labelling, and a nanoflow peptide fractionation approach. The integrated workflow allowed to maximize proteome coverage of laser-isolated tissue samples containing nanogram proteins. We demonstrated the deep spatial proteomics can quantify more than 5,000 unique proteins from a small-sized human pancreatic tissue pixel (∼60,000 µm2) and reveal unique islet microenvironments.

The Impact of the Mass Analyzer and Tissue Section Thickness on Spatial N-Glycomics with MALDI-MSI

We compared matrix-assisted laser desorption/ionization mass spectrometry imaging (MALDI-MSI) spatial N-glycomics data sets from Fourier-transform ion cyclotron resonance (FTICR) and orthogonal accelerated time-of-flight (timsTOF) mass spectrometers of FFPE preserved human kidney samples. We also tested different tissue section thicknesses. In these analyses, we assessed the impact of the mass analyzer and tissue section thickness on N-glycan coverage, sensitivity, and histological alignment. Our results indicate negligible differences in N-glycan coverage between the two mass analyzers, where N-glycan annotation numbers remained consistent and were highly reproducible. The timsTOF-MS analyses demonstrated significant advantages with higher duty cycles and better lateral resolution, allowing for finer spatial resolution without compromising signal integrity. Specifically, timsTOF was able to generate detailed MALDI-MS images at 20 μm step size, accurately identifying N-glycan Hex:5 HexNAc:5 dHex:1 as a tubular-specific marker without observable delocalization. Despite minor annotation discrepancies, where only three species detected by FTICR were not detected by using timsTOF, and a few false-positive annotations from the timsTOF analysis attributed to lower mass resolving power, the overall consistency between the instruments was high. Importantly, tissue section thickness did not affect analysis sensitivity in the timsTOF analyses, with the average glycan signal intensity remaining stable between 7 and 2 μm sections. These findings demonstrate that 2 μm thick tissue slices can be effectively used in spatial N-glycomics workflows, maintaining sensitivity while enhancing confidence in pathohistological evaluations.

From single cell to spatial multi-omics: unveiling molecular mechanisms in dynamic and heterogeneous systems

Single-cell multi-omics and spatial technology have been widely applied to biomedical studies and recently to environmental studies. The cell size detected by single-cell omics ranges from ∼2 µm (e.g., Bacillus subtilis) to ∼120 µm (e.g., human oocytes). Simultaneous detection of single-cell multi-omics is available to human and plant tissues while limited to microbial samples. Spatial technology enables mapping the detected biomolecules in situ. The recent advances in Matrix-Assisted Laser Desorption/Ionization-Mass Spectrometry Imaging and Micro/Nanodroplet Processing in One Pot for Trace Samples for the first time allow the application of spatial multi-omics in highly heterogeneous environmental samples composed of plants, fungi, and bacteria. We envision that these technologies will continue to advance our understanding of unique cell types, their developmental trajectory, and the intercellular signaling and interaction within biological samples.

Spatial top-down proteomics for the functional characterization of human kidney

BACKGROUND

The Human Proteome Project has credibly detected nearly 93% of the roughly 20,000 proteins which are predicted by the human genome. However, the proteome is enigmatic, where alterations in amino acid sequences from polymorphisms and alternative splicing, errors in translation, and post-translational modifications result in a proteome depth estimated at several million unique proteoforms. Recently mass spectrometry has been demonstrated in several landmark efforts mapping the human proteoform landscape in bulk analyses. Herein, we developed an integrated workflow for characterizing proteoforms from human tissue in a spatially resolved manner by coupling laser capture microdissection, nanoliter-scale sample preparation, and mass spectrometry imaging.

RESULTS

Using healthy human kidney sections as the case study, we focused our analyses on the major functional tissue units including glomeruli, tubules, and medullary rays. After laser capture microdissection, these isolated functional tissue units were processed with microPOTS (microdroplet processing in one-pot for trace samples) for sensitive top-down proteomics measurement. This provided a quantitative database of 616 proteoforms that was further leveraged as a library for mass spectrometry imaging with near-cellular spatial resolution over the entire section. Notably, several mitochondrial proteoforms were found to be differentially abundant between glomeruli and convoluted tubules, and further spatial contextualization was provided by mass spectrometry imaging confirming unique differences identified by microPOTS, and further expanding the field-of-view for unique distributions such as enhanced abundance of a truncated form (1-74) of ubiquitin within cortical regions.

CONCLUSIONS

We developed an integrated workflow to directly identify proteoforms and reveal their spatial distributions. Of the 20 differentially abundant proteoforms identified as discriminate between tubules and glomeruli by microPOTS, the vast majority of tubular proteoforms were of mitochondrial origin (8 of 10) while discriminate proteoforms in glomeruli were primarily hemoglobin subunits (9 of 10). These trends were also identified within ion images demonstrating spatially resolved characterization of proteoforms that has the potential to reshape discovery-based proteomics because the proteoforms are the ultimate effector of cellular functions. Applications of this technology have the potential to unravel etiology and pathophysiology of disease states, informing on biologically active proteoforms, which remodel the proteomic landscape in chronic and acute disorders.

Protein N-glycans in Healthy and Sclerotic Glomeruli in Diabetic Kidney Disease

BACKGROUND

Diabetes is expected to directly impact renal glycosylation, yet to date, there has not been a comprehensive evaluation of alterations in N-glycan composition in the glomeruli of patients with diabetic kidney disease (DKD).

METHODS

We used untargeted mass spectrometry imaging to identify N-glycan structures in healthy and sclerotic glomeruli in FFPE sections from needle biopsies of five patients with DKD and three healthy kidney samples. Regional proteomics was performed on glomeruli from additional biopsies from the same patients to compare the abundances of enzymes involved in glycosylation. Secondary analysis of single nuclei transcriptomics (snRNAseq) data was used to inform on transcript levels of glycosylation machinery in different cell types and states.

RESULTS

We detected 120 N-glycans, and among them identified twelve of these protein post-translated modifications that were significantly increased in glomeruli. All glomeruli-specific N-glycans contained an N-acetyllactosamine (LacNAc) epitope. Five N-glycan structures were highly discriminant between sclerotic and healthy glomeruli. Sclerotic glomeruli had an additional set of glycans lacking fucose linked to their core, and they did not show tetra-antennary structures that are common in healthy glomeruli. Orthogonal omics analyses revealed lower protein abundance and lower gene expression involved in synthesizing fucosylated and branched N-glycans in sclerotic podocytes. In snRNAseq and regional proteomics analyses, we observed that genes and/or proteins involved in sialylation and LacNAc synthesis were also downregulated in DKD glomeruli, but this alteration remained undetectable by our spatial N-glycomics assay.

CONCLUSIONS

Integrative spatial glycomics, proteomics, and transcriptomics revealed protein N-glycosylation characteristic of sclerotic glomeruli in DKD.

Coupling Microdroplet-Based Sample Preparation, Multiplexed Isobaric Labeling, and Nanoflow Peptide Fractionation for Deep Proteome Profiling of the Tissue Microenvironment

There is increasing interest in developing in-depth proteomic approaches for mapping tissue heterogeneity in a cell-type-specific manner to better understand and predict the function of complex biological systems such as human organs. Existing spatially resolved proteomics technologies cannot provide deep proteome coverage due to limited sensitivity and poor sample recovery. Herein, we seamlessly combined laser capture microdissection with a low-volume sample processing technology that includes a microfluidic device named microPOTS (microdroplet processing in one pot for trace samples), multiplexed isobaric labeling, and a nanoflow peptide fractionation approach. The integrated workflow allowed us to maximize proteome coverage of laser-isolated tissue samples containing nanogram levels of proteins. We demonstrated that the deep spatial proteomics platform can quantify more than 5000 unique proteins from a small-sized human pancreatic tissue pixel (∼60,000 μm2) and differentiate unique protein abundance patterns in pancreas. Furthermore, the use of the microPOTS chip eliminated the requirement for advanced microfabrication capabilities and specialized nanoliter liquid handling equipment, making it more accessible to proteomic laboratories.

Spatial top-down proteomics for the functional characterization of human kidney

Background

The Human Proteome Project has credibly detected nearly 93% of the roughly 20,000 proteins which are predicted by the human genome. However, the proteome is enigmatic, where alterations in amino acid sequences from polymorphisms and alternative splicing, errors in translation, and post-translational modifications result in a proteome depth estimated at several million unique proteoforms. Recently mass spectrometry has been demonstrated in several landmark efforts mapping the human proteoform landscape in bulk analyses. Herein, we developed an integrated workflow for characterizing proteoforms from human tissue in a spatially resolved manner by coupling laser capture microdissection, nanoliter-scale sample preparation, and mass spectrometry imaging.

Results

Using healthy human kidney sections as the case study, we focused our analyses on the major functional tissue units including glomeruli, tubules, and medullary rays. After laser capture microdissection, these isolated functional tissue units were processed with microPOTS (microdroplet processing in one-pot for trace samples) for sensitive top-down proteomics measurement. This provided a quantitative database of 616 proteoforms that was further leveraged as a library for mass spectrometry imaging with near-cellular spatial resolution over the entire section. Notably, several mitochondrial proteoforms were found to be differentially abundant between glomeruli and convoluted tubules, and further spatial contextualization was provided by mass spectrometry imaging confirming unique differences identified by microPOTS, and further expanding the field-of-view for unique distributions such as enhanced abundance of a truncated form (1-74) of ubiquitin within cortical regions.

Conclusions

We developed an integrated workflow to directly identify proteoforms and reveal their spatial distributions. Where of the 20 differentially abundant proteoforms identified as discriminate between tubules and glomeruli by microPOTS, the vast majority of tubular proteoforms were of mitochondrial origin (8 of 10) where discriminate proteoforms in glomeruli were primarily hemoglobin subunits (9 of 10). These trends were also identified within ion images demonstrating spatially resolved characterization of proteoforms that has the potential to reshape discovery-based proteomics because the proteoforms are the ultimate effector of cellular functions. Applications of this technology have the potential to unravel etiology and pathophysiology of disease states, informing on biologically active proteoforms, which remodel the proteomic landscape in chronic and acute disorders.

Untargeted Spatial Metabolomics and Spatial Proteomics on the Same Tissue Section

An increasing number of spatial multiomic workflows have recently been developed. Some of these approaches have leveraged initial mass spectrometry imaging (MSI)-based spatial metabolomics to inform the region of interest (ROI) selection for downstream spatial proteomics. However, these workflows have been limited by varied substrate requirements between modalities or have required analyzing serial sections (i.e., one section per modality). To mitigate these issues, we present a new multiomic workflow that uses desorption electrospray ionization (DESI)-MSI to identify representative spatial metabolite patterns on-tissue prior to spatial proteomic analyses on the same tissue section. This workflow is demonstrated here with a model mammalian tissue (coronal rat brain section) mounted on a poly(ethylene naphthalate)-membrane slide. Initial DESI-MSI resulted in 160 annotations (SwissLipids) within the METASPACE platform (≤20% false discovery rate). A segmentation map from the annotated ion images informed the downstream ROI selection for spatial proteomics characterization from the same sample. The unspecific substrate requirements and minimal sample disruption inherent to DESI-MSI allowed for an optimized, downstream spatial proteomics assay, resulting in 3888 ± 240 to 4717 ± 48 proteins being confidently directed per ROI (200 μm × 200 μm). Finally, we demonstrate the integration of multiomic information, where we found ceramide localization to be correlated with SMPD3 abundance (ceramide synthesis protein), and we also utilized protein abundance to resolve metabolite isomeric ambiguity. Overall, the integration of DESI-MSI into the multiomic workflow allows for complementary spatial- and molecular-level information to be achieved from optimized implementations of each MS assay inherent to the workflow itself.

Advanced multi-modal mass spectrometry imaging reveals functional differences of placental villous compartments at microscale resolution

The placenta is a complex and heterogeneous organ that links the mother and fetus, playing a crucial role in nourishing and protecting the fetus throughout pregnancy. Integrative spatial multi-omics approaches can provide a systems-level understanding of molecular changes underlying the mechanisms leading to the histological variations of the placenta during healthy pregnancy and pregnancy complications. Herein, we advance our metabolome-informed proteome imaging (MIPI) workflow to include lipidomic imaging, while also expanding the molecular coverage of metabolomic imaging by incorporating on-tissue chemical derivatization (OTCD). The improved MIPI workflow advances biomedical investigations by leveraging state-of-the-art molecular imaging technologies. Lipidome imaging identifies molecular differences between two morphologically distinct compartments of a placental villous functional unit, syncytiotrophoblast (STB) and villous core. Next, our advanced metabolome imaging maps villous functional units with enriched metabolomic activities related to steroid and lipid metabolism, outlining distinct molecular distributions across morphologically different villous compartments. Complementary proteome imaging on these villous functional units reveals a plethora of fatty acid- and steroid-related enzymes uniquely distributed in STB and villous core compartments. Integration across our advanced MIPI imaging modalities enables the reconstruction of active biological pathways of molecular synthesis and maternal-fetal signaling across morphologically distinct placental villous compartments with micrometer-scale resolution.

Changes in maternal blood and placental lipidomic profile in obesity and gestational diabetes: Evidence for sexual dimorphism

Introduction

Obesity and gestational diabetes (GDM) are associated with adverse pregnancy outcomes and program the offspring for cardiometabolic disease in a sexually dimorphic manner. The placenta transfers lipids to the fetus and uses these substrates to support its own metabolism impacting the amount of substrate available to the growing fetus.

Methods

We collected maternal plasma and placental villous tissue following elective cesarean section at term from women who were lean (pre-pregnancy BMI 18.5-24.9), obese (BMI>30) and type A2 GDM (matched to obese BMI) with male or female fetus (n=4 each group). Lipids were extracted and fatty acid composition of different lipid classes were analyzed by LC-MS/MS analysis. Significant changes in GDM vs obese, GDM vs lean, and obese vs lean were determined using t-test with a Tukey correction set at p<0.05.

Results

In placental samples 436 lipids were identified, among which 85 showed significant changes. Of note only in male placentas significant decreases in C22:6 - docosahexaenoic acid (DHA) in phosphatidylcholine (PC) and triglyceride lipid species were seen when comparing tissue from GDM women to lean. In maternal plasma we observed no effect of obesity. GDM or fetal sex.

Conclusion

This is the first study assessing fatty acid composition of lipids in matched maternal plasma and placental tissue from lean, obese, and GDM women stratified by fetal sex. It highlights how GDM affects distribution of fatty acids in lipid classes changes in a sexually dimorphic manner in the placenta.

Sexual dimorphism in lipidomic changes in maternal blood and placenta associated with obesity and gestational diabetes: A discovery study

INTRODUCTION

The placenta uses lipids and other nutrients to support its own metabolism hence impacting the type and amount of these substrates available to the growing fetus. Maternal obesity and gestational diabetes (GDM) can disrupt placental lipid metabolism and thus lead to altered fetal growth contributing to adverse pregnancy outcomes and developmentally programing the offspring for disease in later life. Understanding obesity and GDM driven changes in placental lipid metabolism is thus important.

METHODS

We collected maternal plasma and placental villous tissue following elective cesarean section at term from women who were lean (pre-pregnancy BMI 18.5-24.9), obese (BMI>30) or obese with type A2 GDM n = 8 each group (4 male and 4 female placentas). Fatty acid composition of different lipid classes was analyzed by LC-MS/MS analysis. Significant changes in GDM vs obese, GDM vs lean, and obese vs lean were determined in both a fetal sex-dependent and independent manner.

RESULTS

In placenta 436 lipids were identified, among which 85 showed significant changes. We report significant changes in placental triglyceride, phosphatidylcholine, and phosphatidylinositol lipids containing essential fatty acids- DHA and AA in GDM, with male placentas driving these changes. In maternal plasma, 284 lipids were identified with 14 showing significant changes, but we observed no changes based on fetal sex.

DISCUSSION

Maternal obesity and GDM impact placental lipid composition in a sexually dimorphic manner. The alteration in specific lipid classes can impact cellular energetics and placental function.

Mapping microhabitats of lignocellulose decomposition by a microbial consortium

The leaf-cutter ant fungal garden ecosystem is a naturally evolved model system for efficient plant biomass degradation. Degradation processes mediated by the symbiotic fungus Leucoagaricus gongylophorus are difficult to characterize due to dynamic metabolisms and spatial complexity of the system. Herein, we performed microscale imaging across 12-µm-thick adjacent sections of Atta cephalotes fungal gardens and applied a metabolome-informed proteome imaging approach to map lignin degradation. This approach combines two spatial multiomics mass spectrometry modalities that enabled us to visualize colocalized metabolites and proteins across and through the fungal garden. Spatially profiled metabolites revealed an accumulation of lignin-related products, outlining morphologically unique lignin microhabitats. Metaproteomic analyses of these microhabitats revealed carbohydrate-degrading enzymes, indicating a prominent fungal role in lignocellulose decomposition. Integration of metabolome-informed proteome imaging data provides a comprehensive view of underlying biological pathways to inform our understanding of metabolic fungal pathways in plant matter degradation within the micrometer-scale environment.

Spatial Mapping of Plant N-Glycosylation Cellular Heterogeneity Inside Soybean Root Nodules Provided Insights Into Legume-Rhizobia Symbiosis

Although ubiquitously present, information on the function of complex N-glycan posttranslational modification in plants is very limited and is often neglected. In this work, we adopted an enzyme-assisted matrix-assisted laser desorption/ionization mass spectrometry imaging strategy to visualize the distribution and identity of N-glycans in soybean root nodules at a cellular resolution. We additionally performed proteomics analysis to probe the potential correlation to proteome changes during symbiotic rhizobia-legume interactions. Our ion images reveal that intense N-glycosylation occurs in the sclerenchyma layer, and inside the infected cells within the infection zone, while morphological structures such as the cortex, uninfected cells, and cells that form the attachment with the root are fewer N-glycosylated. Notably, we observed different N-glycan profiles between soybean root nodules infected with wild-type rhizobia and those infected with mutant rhizobia incapable of efficiently fixing atmospheric nitrogen. The majority of complex N-glycan structures, particularly those with characteristic Lewis-a epitopes, are more abundant in the mutant nodules. Our proteomic results revealed that these glycans likely originated from proteins that maintain the redox balance crucial for proper nitrogen fixation, but also from enzymes involved in N-glycan and phenylpropanoid biosynthesis. These findings indicate the possible involvement of Lewis-a glycans in these critical pathways during legume-rhizobia symbiosis.