decreasing CV. Sphingolipids populated the
top CV ranking positions, confirming that the
sphingolipid pathway (fig. S1F) is subjected to
high cell-to-cell variability (Fig. 2B).
To test whether the sphingolipid metabolic
pathway is coordinately modulated or if differ-
ent sphingolipid subsets are controlled in-
dependently in different cells, we created a
pairwise lipid-lipid correlation (PearsonÕsR)
matrix (Fig. 2C). Although phospholipid spe-
cies did not form biochemically meaningful
cliques, sphingolipids were clustered in groupsconsisting of compounds bearing the same head
group (i.e., OH with Cers, hexose with HexCers,
trihexose with Gb3s, andN-acetyl-hexose-
trihexose with Gb4s) but with different Cer
backbones (mostly 34:1, 40:1, 42:1, and 42:2)
(fig. S1F). This suggests that specific enzymaticCapolupoet al.,Science 376 , eabh1623 (2022) 15 April 2022 2 of 12
Fig. 1. Single-pixel MALDI-
MSI analysis on dHFs.
(A) Schematic drawing of
MALDI-MSI workflow.
Cells were fixed, matrix
was deposited, and MALDI-
MSI was performed by
rasterizing the laser across
a selected area. For each
spot, a mass spectrum was
collected and mass images
were obtained for each
ion by plottingm/zintensity
at the correspondingx
andycoordinates (left panel).
For peak identification,
total lipid extracts were
analyzed by ESI-LC/MS
(right panel). Lipids identi-
fied by ESI-LC/MS were
then compared with
the ones obtained by
MALDI-MSI. (B) Ion images
(50mm^2 /pixel; 354 ×
218 pixels) of selected lipids
recorded in positive-ion
mode. Insets show individual
cells images at higher mag-
nification. [PC(34:1)+Na]+,
phosphatidylcholine with
acyl chains consisting
of 34 carbon atoms and
one double bond complexed
with sodium; [SM(34:1)+Na]+,
SM with a backbone of
34 carbon atoms and
one double bond complexed
with sodium; [PC O-(36:4)
+H]+, phosphatidylcholine
plasmalogen with acyl
chains consisting of
36 carbon atoms and
four double bonds
complexed with hydrogen;
TIC, total ion current. Scale
bar, 500mm. (C) Images
displaying at each location
the PCA coordinate of each
pixel. PC1, PC5, and PC6
values are displayed using
a divergent color map;
positive coordinates are
shown in red and negative in
blue. Insets show individual
cell images at higher magnification. (D) Bar plots showing the contribution of the top 10 lipids with higher (red) and lower (blue) loadings for
each PC. Miniatures in the upper left corner show the entire distribution.
A BCmin 0 Max−0.150.00.15−0.20.00.2−0.20.00.2D10
45.56
175504.26486
1 124.542371084
.54922SM 41:1SM 40:1SM 38:1SM4 1:2
CerP 46:6SM4 2:1
SM40:2
PC 32:2
SM 36:0PA 34:1PC P-36:5PE O-40:^6
Gb4 40:1
Gb4 42:1Gb4 34:1
Gb4 42:20.20.2-0.2Loadings PC6Loadings PC5Loadings PC10.1-0.1-0.213
44.6423310
35.517861040.5441912
69.59146
1080.57622
Gb4 42:1Gb4 42:2
Gb43 4:1Gb4 40 :1
PA O - 37:0PC 33:2
Gb3 42:2Gb334:1
800.55
65
2PE O-^4 0:
61 002.57985
Gb3 40:1
SM 36:0
1160.7928
9Gb342:110
3 8.58432
1037.5808
1053.^5
5438
1213.58968 1 077.57352
69
7.526231173.59751045.^56175CerP36:^3
CerP46:6
PA
O
-37:0PC 40:6
Cer 34 :1
PC32:2
Cer 42 :1
PE34
:1
PE 38:5PC
39:^5
PE
3 8:2
PC-O 34:1PC #1 PC Coordinate
PC #5
PC #6 PC #6
Lipid SpeciesMin MaxSpectrum Peak Intensity
Normalized to TIC[PC O-(36:4)+H]+
m/z 768.5917[SM(34:1)+Na]+
m/z 725.5561[PC(34:1)+Na]+
m/z 760.585Matrix depositionLaserIntensity m/zPixelsMALDI-MSILipid #1 Lipid #2 Lipid #205...Optical ImageESI-LC/MSLipid callingLipid #205RTIntensityLipid #2PO
OOOOOHON
OPO
NHOOOH
HON
OLipid #1HOOOH
OHOH
ONH
HO OTAG
Cer
GlcCerPCPSPELBPA PGPAPI SMPC #5PC #1RESEARCH | RESEARCH ARTICLE