Lipidomics on nanoneedles
We assessed the suitability of nanoneedles for DESI-MSI lipidomics, utilizing nanoneedle arrays with 2 μm pitch, 4.00 μm peak, 50 nm tip diameter, 600 nm base diameter, 77.2 m2 g−1 floor space, 50% porosity and seven.9 nm common pore diameter, on 8 × 8 mm silicon substrates (Supplementary Fig. 1). The evaluation of porcine mind lipid extract revealed that nanoneedles and flat surfaces supply comparable sensitivity for essentially the most consultant lipids, in keeping with established DESI limits of detection (LODs)19 (Prolonged Information Fig. 1a,b and Supplementary Desk 1). Nanoneedle substrates sometimes displayed barely increased LODs than flat silicon, suggesting that nanotextured surfaces can dampen the desorption–ionization course of. Nevertheless, the relative peak depth correlated properly between nanoneedles and flat substrates, indicating that the evaluation on nanoneedles precisely recapitulated the lipid composition of the combination (Prolonged Information Fig. 1c and Supplementary Fig. 2). Optimization of DESI-MSI analytical parameters20 for nanoneedles considerably elevated the variety of detectable options by an element of two.47, from 408 ± 83 to 1,007 ± 141 (Prolonged Information Fig. 1d–f).
Lipidomics imaging of molecular replicas
We then established the flexibility of nanoneedles to generate a molecular reproduction by imprinting frozen mouse mind (Fig. 2a,b). Computational modelling indicated that, throughout imprinting, a number of micrometres of tissue transiently thaw, permitting assortment of biomolecules (Supplementary Fig. 3). The ensuing reproduction fashioned a uniform molecular adsorbate on the nanoneedles, whereas utilizing a flat substrate yielded poor adsorbate uniformity, presenting areas with out protection and areas of adsorbate accumulation (Fig. 2c). The nanoneedle reproduction absolutely reproduced the unique tissue and macroscopically retained its morphology, whereas the flat reproduction reproduced solely lower than half of the tissue, missing clear morphology (Fig. second,e).
a, Schematics of the nanoneedle imprinting system. Created with BioRender.com. b, Images illustrating the system approaching (high) and imprinting (backside) nanoneedles on tissue. c, Scanning electron microscopy photographs of the molecular adsorbate on flat and nanoneedle substrates after imprinting. d, Images of the imprinted tissue and fluorescence microscopy photographs of the corresponding molecular reproduction. e, Quantification of replica-to-tissue space ratio. Information are imply ± s.e.m. (N = 3 unbiased organic replicates). Statistics: two-tailed unpaired t-test. f, Quantification of RNA, lipids and proteins from nanoneedle replicas with completely different floor chemistries (as-etched, APTES, oxidized) and a 10-µm mind tissue part. ata are imply ± s.e.m. (N = 3 unbiased organic replicates). Statistics: two-way ANOVA with a post-hoc Tukey’s multiple-comparisons take a look at. g, DESI-MSI maps of gray matter (m/z 834.53, blue) and white matter (m/z 888.62, inexperienced) lipids from a tissue slice and molecular replicas generated utilizing porous silicon nanoneedles (pSi nN), strong silicon nanoneedles (strong nN), porous silicon nanoneedles over a porous silicon layer (pSi nN + layer) and a porous silicon layer (pSi layer). h, Quantification of the variety of options as a operate of the boldness threshold set for peak identification in samples from g. Information are imply ± s.e.m. (N = 3 unbiased organic replicates). i, The function depend at a peak threshold of two.0 (95% confidence) for samples from g. Information are imply ± s.e.m. (N = 3 unbiased organic replicates). j, The correlation matrix of spectra from three tissue sections and their corresponding replicas from g. ok, The relative floor space of molecular replicas to unique tissue for samples from g. Information are imply ± s.e.m. (N = 3 unbiased organic replicates). Statistics: extraordinary one-way ANOVA with post-hoc Tukey’s multiple-comparisons take a look at. l, A correlation matrix evaluating three pSi nN replicas with the unique tissue part. m, A correlation scatter plot of lipids peak intensities between the tissue part and the primary reproduction. n, A correlation scatter plot of lipid peak intensities between the primary and second reproduction.
We then evaluated the composition of the molecular adsorbate on replicas generated from nanoneedles as-etched, O2-plasma oxidized and (3-aminopropyl)-triethoxysilane (APTES) functionalized, assessing the function of floor polarity and cost on biomolecular assortment effectivity21,22. Oxidation and APTES respectively generated destructive and constructive fees, offering extra polar surfaces than the as-etched porous silicon. For all of the functionalizations thought-about, the quantity of collected RNA, proteins and lipids in contrast favourably with people who may very well be extracted from the reference tissue (Fig. 2f). RNA assortment ranged between 210 ng and 350 ng, whereas protein assortment ranged between 2.5 μg and 4.3 μg. Lipid assortment was extra environment friendly on APTES and oxidized nanoneedles in contrast with as-etched ones, ranging between 200 μg and 600 μg. These information indicated that plasma oxidized nanoneedles had been appropriate for spatial lipidomics, as a result of their lipid assortment effectivity of 410 ± 95 μg was properly inside the vary by which a lot of options had been detectable (Fig. 2f) and in contrast favourably with a 10-μm tissue part (210 ± 23 μg), recognized to be imageable by DESI-MSI. The usage of oxidized nanoneedles additionally mitigated potential issues relating to interference from the desorption of polymerized APTES.
We additionally in contrast the standard of the reproduction generated utilizing porous silicon nanoneedles (pSi nN) with these obtained utilizing strong silicon nanoneedles (strong nN), porous silicon nanoneedles over a steady porous silicon layer (pSi nN + layer) and a steady porous silicon layer (pSi layer)23,24,25 (Supplementary Fig. 4). pSi nN displayed the most important variety of options and the least variability, adopted by the strong nN, whereas the pSi nN + layer and pSi layer yielded much less options with increased variability (Fig. 2g–i). The pSi nN additionally outperformed the opposite substrates of their spectral correlation to the unique tissue, carefully matched by strong nN (Fig. 2j), and most carefully reproduced the floor space of the unique tissue (Fig. 2g,ok).
The imprinting was reproducible, with three replicas from the identical floor preserving morphological options and displaying robust correlation with the reference part and with one another throughout the entire vary from low- to high-abundance lipids (Fig. 2l–n and Supplementary Fig. 5).
These information indicated {that a} molecular reproduction generated utilizing porous silicon nanoneedles may precisely and robustly reproduce the lipidomic profile of a tissue. Nevertheless, the presence of a steady porous layer lowered the effectivity of the desorption–ionization course of, favouring using substrates with nanoneedles solely. As well as, porous nanoneedles marginally outperformed strong ones, producing replicas that extra carefully represented the unique tissue by way of space, variety of options and correlation.
Spatial lipidomics of murine glioma replicas
We in contrast the standard of lipidomic imaging on molecular replicas with its proximal tissue part5. Hierarchical cluster evaluation (HCA) revealed that the highest two clusters of each part and reproduction appeared to align with the areas of gray and white matter inside the mouse mind structure (Fig. 3a). The white matter cluster aligned with mind areas containing extremely myelinated fibres, matching the expression sample of myelin genes akin to myelin fundamental protein in areas such because the corpus callosum. In the meantime the gray matter cluster corresponded to areas with increased cell physique densities and fewer myelinated fibres, such because the cortex4.
a, HCA maps of the highest two clusters from the DESI-MSI information of a murine mind part and its molecular reproduction, similar to gray matter (blue) and white matter (yellow). b, Common spectra displaying normalized counts versus m/z, for the highest two clusters related to gray matter (GM) and white matter (WM) in part (S) and reproduction (R). c, Volcano plots of the differential lipid distribution between gray and white matter, highlighting the function of m/z 834.53 (PS 40:6) and 888.62 (SHexCer 42:2;O2) as gray and white matter biomarkers in each part and reproduction. d, An overlay map of TIC-normalized depth for m/z 834.53 (blue) and 888.62 (yellow). e, HCA maps of the highest three clusters from the DESI-MSI dataset of a tumour-bearing mouse mind part and reproduction, similar to gray matter (GM, blue), white matter (WM, yellow), and tumour (T, purple). White bins point out areas used for lipid abundance evaluation. f, Common spectra displaying normalized counts versus m/z, for the three clusters related to gray matter (GM), white matter (WM) and tumour (T) in part (S) and reproduction (R). g, Volcano plots evaluating lipid distribution between GM and T, highlighting essentially the most consultant lipid peaks in each part and reproduction. h, An overlay map of cumulative TIC-normalized depth for peaks marking gray matter (m/z 834.53 and m/z 600.51, blue), white matter (m/z 844.64, m/z 860.64 and m/z 888.62, yellow) and tumour (m/z 682.59, m/z 716.52 and m/z 736.65, purple) in part and reproduction. The white bins point out areas used for lipid abundance evaluation throughout GM, WM and T. i,j, Heatmaps displaying the relative abundance of 33 lipid peaks (from evaluation in g) alongside the foremost axis of the white field in e and h for the part (i) and reproduction (j). Prime and backside strips above every heatmap present the cluster project and relative lipid depth for corresponding pixels from e and h.
The lipid composition of those two clusters confirmed this morphological classification (Fig. 3b,c). The upper abundance of m/z 888.62 (sulfated hexosylceramide (SHexCer) 42:2;O2) and decrease abundance of m/z 834.53 (phosphatidylserine (PS) 40:6) in white matter had been in keeping with established DESI-MSI mind biomarkers26. Mapping the relative depth of those two peaks throughout the pattern precisely delineated white and gray matter areas in each the reproduction and the reference part (Fig. 3d). The congruent distribution of those lipid species throughout reproduction and part confirmed the feasibility of producing a lipidomic profile map from a nanoneedle reproduction of the mouse mind. Minor discrepancies in mind morphology had been attributed to the inherent anatomical variations between adjoining tissue sections.
We then used the molecular reproduction to characterize a cancerous lesion inside a mouse glioma5. The highest three HCA clusters categorized white matter, gray matter and tumour areas with comparable effectiveness in part and reproduction (Fig. 3e). When analysing a dataset together with all spectra throughout tissue and reproduction, HCA equally discriminated tumour, white matter and gray matter, highlighting the correspondence of the lipid profile on reproduction and part. Principal element evaluation of intercluster variance additional confirmed that each the reproduction and the part displayed comparable attribute options, enabling correct mapping of tumour margins and transitions between tissue sorts (Supplementary Fig. 6).
Comparative evaluation of the typical mass spectra for every cluster allowed the lipid compositions of white matter, gray matter and tumour to be distinguished and essentially the most related differentially plentiful species to be recognized (Fig. 3f,g and Supplementary Fig. 7). Once more, clusters related to white matter in each the part and reproduction exhibited a distinguished peak at m/z 888.62 (SHexCer 42:2;O2), whereas the gray matter clusters confirmed a distinguished peak at m/z 834.53 (PS 40:6). The tumour cluster confirmed a rise within the abundance of m/z 736.65 (ceramide (Cer) 46:1;O4), m/z 682.59 (Cer 42:0;O4) and m/z 716.52 (phosphatidylethanolamine (PE) 34:1). The depth map for these attribute species successfully discriminated tumour, white matter and gray matter within the reproduction in addition to within the part (Fig. 3h), demonstrating the reproduction’s capacity to precisely reproduce the lipidomic map of the unique tissue, even within the presence of a cancerous lesion. The evaluation of two further, distinct mind samples confirmed the reproducibility of our strategy (Supplementary Fig. 8).
We assessed the spatial distribution of the 33 most consultant species within the part and reproduction by evaluating their abundance alongside the foremost axis (size) of the white field (Fig. 3h and Supplementary Fig. 9). The heatmaps from each part and reproduction revealed a transparent and congruent sample in lipid abundance throughout tumour, gray matter and white matter (Fig. 3i,j and Supplementary Fig. 9). The tumour–gray matter transition occurred at 3.1 mm from the left edge. The tumour was wealthy in lipids at m/z 736.65 (Cer 46:1;O4), m/z 682.59 (Cer 42:0;O4) and m/z 740.51 (phosphatidylcholine (PC) 30:0), whereas gray matter confirmed excessive ranges of m/z 600.51 (Cer 36:1;O2) and m/z 834.53 (PS 40:6). At 4 mm, a 0.3 mm white matter layer was enriched in m/z 844.64 (PS O-41:1), m/z 860.64 (PS 41:0) and m/z 888.62 (SHexCer 42:2;O2) adopted by gray matter (4.3–4.7 mm) with a lipid profile in keeping with the earlier gray matter area.
These information point out {that a} molecular reproduction can map the biomolecular structure of a tumour-bearing tissue. The standard of the reproduction info permits unsupervised classification of particular person pixels with comparable constancy to the unique pattern enabling mapping lipid abundance with comparable high quality to the tissue.
Spatial lipidomics of HG reproduction
We then captured the molecular heterogeneity of a HG biopsy. Histopathological evaluation recognized three distinct areas of wholesome mind infiltrated by tumour (I), bulk tumour (T) and necrosis (N) (Fig. 4a). HCA may distinguish these three areas and establish their lipid composition within the reproduction and the part (Fig. 4b–e). Mapping the distribution of the height at m/z 600.51 (Cer 36:1;O2) characterised the infiltrated tumour, whereas m/z 888.62 (SHexCer 42:2;O2) distinguished the majority tumour and m/z 680.54 (Cer 43:0;O3) preferentially related to areas of necrosis (Fig. 4f).
a, A bright-field picture of the haematoxylin and eosin (H&E)-stained part of HG biopsy, displaying infiltrated tumour (I), bulk tumour (T) and necrosis (N). b,c, HCA maps of the highest three clusters from DESI-MSI photographs of the part (b) and reproduction (c). d, Volcano plots displaying the differential lipid abundance evaluation throughout I, T and N. constantly figuring out m/z 600.51 (Cer 36:1;O2) for I, m/z 888.62 (SHexCer 42:2;O2) for T and m/z 680.54 (Cer 43:0;O3) for N in part and reproduction. e, A heatmap of the ten most consultant lipids defining I, T and N in each part and reproduction. f, An overlay map of TIC-normalized depth for m/z 600.51 (Cer 36:1;O2, purple), 888.62 (SHexCer 42:2;O2, blue) and 680.54 (Cer 43:0;O3, yellow).
Machine studying inference of HG grade
To judge whether or not replicas may match the predictive energy of unique tissues, we collected DESI-MSI maps of molecular replicas (HG-r) and reference sections (HG-s) from 23 HG samples with recognized prognoses (Supplementary Desk 2). Further sections (HG 11, 12) and replicas (HG 6, 18) had been collected for intrapatient and interreplica validations, for a complete of 25 sections and 25 replicas.
We noticed a excessive correlation between replicas and their reference tissue part, indicating a preserved molecular profile (Fig. 5a). The correlation with the reference part ranked first with a likelihood of 0.85, and ranked within the high 5 with a likelihood of 0.96, properly above the random expectations of 0.04 and 0.2, respectively (Fig. 5b). Replicas from the identical part preserved a standard molecular signature, displaying increased correlation with their reference part than with another (HG 6, 18) (Fig. 5c).
a, A heatmap of correlation values between 25 tissue sections and 25 replicas from 23 biopsy samples. Samples are labelled HG plus a sequential quantity; a number of sections or replicas from the identical biopsy embrace a trailing underscore and quantity. Sections are marked with ‘s’ and replicas with ‘r’. The purple bins point out matched part–reproduction pairs. b, A depend plot of correlation rank between matching tissue sections and molecular replicas. c, Pearson correlation evaluation of the molecular replicas HG6_1-r, HG6_2-r, HG18_1-r and HG18_2-r, with all 25 sections. The arrows mark the reference HG6-s (blue) and HG18-s (orange) part. d, Receiver operator traits for the logistic regression classifier utilized to the HG dataset. Information are imply ± s.d. from 100 random-seed classifications. e, Receiver operator traits for logistic regression cross-inference: mannequin educated on sections, examined on replicas and vice versa. Information are imply ± s.d. from 100 random-seed classifications. TPR, true constructive charge; FPR, false constructive charge. f,g, SHAP worth distributions of top-ranked options within the reproduction (f) and the part mannequin (g). Optimistic SHAP values point out contribution to high-grade predictions. h, TIC-normalized ion maps of chosen replicas (R) and tissue sections (S) displaying the highest 4 peaks from f, alongside H&E-stained morphology.
We subsequent assessed illness grade prediction utilizing sections and replicas evaluating 4 classification fashions: logistic regression, determination timber, XGBoost27 and LightGBM28. All fashions carried out higher than a random mannequin for each sections and replicas (Fig. 5d and Supplementary Fig. 10), with illness grade assignments displaying a excessive diploma of concordance. Logistic regression demonstrated superior efficiency (Fig. 5d and Supplementary Fig. 10) and thus was thought-about for additional evaluation. The classification efficiency was comparable for sections and replicas with a imply (±s.d.) space below the curve (AUC) for the receiver working attribute (ROC) of 0.71 ± 0.05 and 0.75 ± 0.02, respectively (Fig. 5d and Supplementary Desk 3). Bootstrapping statistical validation confirmed that the noticed classification efficiency was not as a result of probability (Supplementary Fig. 11). One of the best-performing fashions agreed on 25 out of the 27 part–reproduction pairs, additional indicating that the fashions captured comparable discriminative info from the spectral information (Supplementary Desk 2). Cross-model prediction evaluation making use of the fashions educated on sections to foretell the replicas achieved an AUC of 0.79 ± 0.01, whereas fashions educated on replicas to foretell the sections achieved an AUC of 0.70 ± 0.05, additional supporting the similarity of their molecular profiles (Fig. 5e and Supplementary Desk 3). Shapley additive explanations (SHAP)29 confirmed that part and reproduction shared three of the highest 5 most predictive peaks, whereas 9 out of the highest ten displayed the identical correlation with lesion grade (Fig. 5f,g). Among the many high 4 reproduction peaks, two (m/z 790.55 (PE 40:6) and m/z 834.53 (PS 40:6)) correlated with low-grade tumours and two (m/z 885.55 (phosphatidylinositol (PI) 38:4) and m/z 788.55 (PS 36:1)) with excessive grade. Their depth maps throughout consultant tumours aligned with anticipated tumour grades (Fig. 5h). This evaluation additionally captured the spatial heterogeneity of tumours, figuring out the decrease part of HG 25—histologically infiltrated wholesome mind—as low grade, enriched within the m/z values of 790.55 and 834.53, and the sides of HG 23 as excessive grade, enriched in m/z 788.55 and depleted in low-grade markers. These outcomes established that the molecular reproduction generated utilizing nanoneedles can change a tissue part for illness state classification.
Molecular replicas of residing tissues
The non-destructive molecular assortment enabled by nanoneedles provides a singular alternative to pattern residing tissue longitudinally, permitting one to trace its spatiotemporal lipidomic evolution. We examined this capacity by monitoring the response of reside glioma-bearing tissue slices to temozolomide (TMZ) remedy. We collected longitudinal molecular replicas on day 0 (D0T1, handled; D0T0, management) and day 5 (D5T1, handled; D5T0, management) whereas exposing the samples to TMZ or dimethylsulfoxide from day 3 (Fig. 6a).
a, A schematic of the longitudinal evaluation workflow. Created with BioRender.com. b, Images (left) of mind 1 (B1) displaying handled (T1) and untreated (T0) tissue slices at day 0 (D0) and day 5 (D5), with corresponding DESI-MSI maps (proper) of gray matter (m/z 600.51, blue), white matter (m/z 888.62, yellow) and tumour (m/z 682.59, purple) markers from molecular replicas of the photographed acute tissue slice. c, Quantification of gray matter, white matter and tumour areas from DESI-MSI molecular replicas information throughout time and remedy in three mind slices. Statistics: extraordinary two-way ANOVA with post-hoc Šidák multiple-comparisons take a look at. n.s., not important. d, A volcano plot of lipid modifications over time in untreated brains, highlighting markers for white matter (yellow), gray matter (blue) and tumour (purple) and considerably elevated (brown) or decreased (purple) species. N = 3 mind slices. e, A heatmap of the differentially plentiful species (orange and purple) proven in d. f–h, Violin plots of gray matter (f), white matter (g) and tumour (h) marker abundance in mind 1 throughout time and remedy. Information: strong line: median; dashed strains: higher and decrease quartiles. i, A heatmap of z scores for treatment-dependent important lipids. j, DESI-MSI maps of handled and untreated molecular replicas at day 0 and day 5 for the lipid at m/z 692.45 (PS 29:0, purple) within the tumour alongside gray matter (m/z 600.51, Cer 36:1;O2, blue) and white matter (m/z 888.62, SHexCer 42:2;O2, yellow) markers. ok, Violin plots of m/z 692.45 (PS 29:0) abundance. Information: strong line: median; dashed strains: higher and decrease quartiles.
We first assessed the feasibility of producing molecular replicas from residing tissue with minimal perturbation. Nanoneedle imprinting didn’t affect mind slice viability instantly (Prolonged Information Fig. 2a,b) or over the 5 days of tradition (Prolonged Information Fig. 2c,d), whereas the TMZ remedy selectively lowered the viability of cells within the tumour area by day 5 (Prolonged Information Fig. 2e–g). The day 0 and day 5 replicas captured all the tissue slice space (Prolonged Information Fig. 2h).
Spatiotemporal lipidomics of molecular replicas
We used these longitudinal molecular replicas to check the spatiotemporal variations within the lipidomic profile in response to TMZ remedy. The DESI imaging of the molecular replicas at each timepoints confirmed well-defined mind areas (Fig. 6b and Supplementary Fig. 12a,c). Tumour, gray matter and white matter had been distinguished by the identical lipids recognized for the frozen sections (Figs. 3h and 6b). The morphology of those areas aligned properly with the anticipated construction of the tumour-bearing mind and carefully matched the mind construction noticed within the slices, which developed over the 5 days in tradition. Over time, and no matter remedy, the dimensions of the gray matter barely elevated and the white matter decreased, in keeping with the anticipated axonal degeneration in cultured mind slices30,31 (Fig. 6c). Against this, the tumour dimension barely elevated in untreated samples whereas reducing considerably in response to TMZ remedy, supporting the statement that the remedy selectively focused tumour cells (Fig. 6c and Prolonged Information Fig. 2g).
We then assessed the temporal evolution of lipid composition by evaluating lipid abundance between untreated samples at day 0 (D0T0) and day 5 (D5T0) throughout the three brains. Among the many 629 generally recognized lipids on day 0, 21% elevated in abundance over time, whereas 79% decreased with 31 of them changing into undetectable (Fig. 6d). Amongst these, 27 lipids confirmed important modifications, with 8 rising and 19 reducing in abundance (Fig. 6e). The white matter and tumour markers didn’t range considerably, whereas the gray matter marker abundance elevated over time (Fig. 6f–h and Supplementary Fig. 12b,d). The info point out an general lower within the abundance and variety of species as a consequence of culturing the mind slices. This evaluation captures mind tissue remodelling throughout culturing and divulges modifications in lipid composition and spatial distribution inside particular person slices, insights made doable via longitudinal sampling enabled by nanoneedles.
We then analysed the impact of remedy particularly inside the tumour. Total, the relative abundance of 42% of lipids elevated with remedy, whereas 58% decreased. Statistical evaluation recognized 11 important lipids (Fig. 6i), with 8 reducing and three rising in abundance with TMZ remedy. Essentially the most important, m/z 692.45 (PS 29:0), selectively disappeared in handled tumours throughout all samples (Fig. 6j,ok). Different lipid species confirmed distinct, treatment-dependent abundance patterns, highlighting the facility of longitudinal molecular reproduction evaluation to detect lipid-specific and specimen-specific remedy responses (Supplementary Fig. 13). Phosphatidylserines (PS29:0, PS 42:8, PS O-41:0, PS 38:4, PS 36:1) and phosphatidylinositol (PI O-39:1) featured prominently among the many species reducing in response to TMZ (Supplementary Desk 1). Phosphatidylserine upregulation is a central function of glioblastoma multiforme32 that’s modulated by TMZ remedy33, together with via ferroptosis34,35,36. The PI3K–AKT is a ubiquitously upregulated pathway in glioblastoma37, which is leveraged by TMZ38.
These information present that TMZ remedy alters tumour lipid composition and distribution, notably by selectively lowering particular metabolites. Utilizing longitudinal information from molecular replicas, we characterised the response of the tumour and recognized goal lipids beforehand linked to glioma and TMZ. This evaluation provides invaluable insights into metabolic pathways affected by the remedy and potential biomarkers of therapeutic efficacy.