空间代谢组学：单细胞空间代谢流分析新方法

生物体内的代谢物和脂质不仅是细胞的关键组成模块，它们在信号传导、表观基因组调控、免疫、炎症和癌症发展中同样具有重要作用和意义。代谢组学分析是我们了解、评估生物体、器官和细胞状态的重要方式。而单细胞技术通过展示组织内部甚至单克隆细胞之间的细胞异质性，将生物学研究推进至新维度。质谱成像（MSI）技术可以从样品中创建特定化合物的图像，这些图像是由样品表面获得的数千个质谱生成的。每个记录的质谱都会为图像贡献一个像素，而每个质谱中的峰都可以生成一个图像。与其他成像方法相比，MSI无需化合物标记，可实现非靶向分析。

研究背景

脂质在生物体生命过程中承担着多种重要作用，多数脂质是由脂肪酸合成而来。成年哺乳动物体内的细胞通常由血液中摄取脂肪酸，而脂肪、肝脏以及癌细胞还可以Acetyl-CoA为底物，从头合成脂肪酸^[2]。Acetyl-CoA经过一系列代谢反应，可以生成含有16个碳的饱和脂肪酸棕榈酸（16:0），之后棕榈酸发生碳链延长或去饱和反应生成不同的饱和、不饱和脂肪酸，从而影响脂质组成。而Acetyl-CoA同样有多种来源，除了葡萄糖经由TCA循环生成的柠檬酸在ACLY作用下生成Acetyl-CoA以外，在缺氧环境下，葡萄糖后续代谢产物丙酮酸会转化为乳酸，从而无法合成Acetyl-CoA、进入脂肪酸合成途径。在此情况下，谷氨酰胺可通过还原羧化反应生成柠檬酸，进而合成Acetyl-CoA ^[3,4] 。另有文献报道，缺氧环境下的癌细胞还可以将乙酸作为脂肪酸合成的前体 ^[5,6] 。而Acetyl-CoA除了作为脂肪酸合成底物以外，对于蛋白翻译后修饰、基因表达等均有重要作用。通过监控脂肪酸合成和Acetyl-CoA代谢间的互动可以帮助我们深入理解癌细胞的生存状态。

分析方法

大气压MALDI成像分析是通过AP-SMALDI5离子源配合Q Exactive plus高分辨质谱仪实现的。激光像素设置为 10×10 μm，激光衰减器角度设置为33°。质谱在负离子模式下采用一级全扫描和全离子碎裂（AIF）扫描模式。AIF模式的隔离范围为 m/z 600-1000，扫描范围为m/z 100-400，分辨率 140k，最大注入时间500 ms，碰撞能量NC 25%。（图1）

图1. 单细胞代谢流质谱成像分析流程

（点击查看大图）

MALDI分析前后，分别应用显微镜检测，确定细胞影像位置及MALDI消融标记位置。通过检测MALDI的消融标记，将其与细胞影像叠加，并通过应用数学公式进行解卷积，从而整合显微镜图像和MALDI图像。实现了应用MALDI成像质谱检测到的单细胞分子轮廓。（图2）

图2. 整合显微镜和MALDI-MS分析结果实现单细胞质谱成像（点击查看大图）

鼠肝癌细胞常氧-低氧模型单细胞成像分析

鼠肝癌细胞在添加25 mM的¹²C-葡萄糖或U-¹³C-葡萄糖后，用含1mM醋酸、2 mM谷氨酰胺和10%透析胎牛血清的无葡萄糖DMEM细胞培养基培养，在37°C、5% CO₂的培养箱中在常氧（20% O₂）或低氧（0.5% O₂）条件下培养72小时。选择72小时的时间点是为了确保棕榈酸的同位素标记已经达到稳态。（图3）在低氧条件下培养的细胞被表达绿色荧光蛋白（GFP）标记。在共培养实验中，常氧和低氧细胞使用胰酶分离，每种条件下混合10000个细胞，在同一张玻璃片上进行培养，并在固定之前允许其附着3小时。

图3. 由稳定同位素标记的¹³C₆-葡萄糖生成细胞质Acetyl-CoA以及后续的脂肪酸和脂质合成途径

（点击查看大图）

通过质谱一级全扫描分析，质谱成像共检测到64种脂质，包括磷脂酸（PA）、磷脂酰肌醇（PI）、磷脂酰乙醇胺（PE）、磷脂酰丝氨酸（PS）等。具体脂质鉴定结果经过了常规LCMS脂质分析确认。在AIF模式下，检测到了11种含量最高的脂肪酸，相应检测结果同样与常规LCMS分析结果相符。

为了验证本方法，研究人员检测了常氧-低氧培养的鼠肝癌细胞混合样本。通过对氨基酸同位素峰的定量分析，发现¹³C标记的棕榈酸（M>0）主要在正常细胞中检出，而缺氧细胞中的棕榈酸以未标记状态（M+0）为主。通过GFP标记结果的对照，证明了本方法可以通过同位素峰分布有效识别不同培养状态的细胞。

图4. 在常氧（GFP阴性）和低氧（GFP阳性）条件下的原代鼠肝癌细胞共培养模型的显微镜和质谱成像结果

（点击查看大图）

图5. 通过GFP标记验证识别不同培养模式细胞的准确性（点击查看大图）

单细胞Acetyl-CoA池标记水平分析

研究人员使用了两种表达不重叠的shRNA序列（ACLYkd oligo1和ACLYkd oligo 2）细胞系以及一个对照组细胞系。通过使用1 μg/mL的四环素处理细胞72小时实现了ACLY沉默。

质谱成像数据是以10 μm的像素大小获得的，每个细胞的平均面积为550μm²，平均每个细胞有12个像素。通过应用二项式模型计算每个细胞的acetyl-CoA池标记程度p值，从而量化细胞质中acetyl-CoA池中从葡萄糖衍生的同位素标记acetyl-CoA的比例。测试结果与预期相符，ACLYkd细胞中的acetyl-CoA池标记水平低于对照组。值得注意的是，两种ACLYkd细胞之间的差异非常明显。ACLYkd oligo1的结果呈双峰分布，p值的差异明显较大，表明该细胞系存在两个亚群体。其中一个模式显示的p值与对照组相近，说明存在一个“沉默失败”的细胞亚群。ACLYkd oligo1第二个模式具有的p值明显则低于ACLYkd oligo 2，表明ACLYkd oligo 1中还存在一个“强沉默”的亚群，在这些细胞中，沉默效率非常高，导致acetyl-CoA同位素标记比例大幅降低。在ACLYkd oligo 2中，acetyl-CoA池的标记程度以及GFP报告基因强度显示出更均一的分布。M+2峰是最能表现出ACLYkd oligo1细胞中“强沉默”群体的低acetyl-CoA标记表型的质谱峰。M+8峰则为对照组细胞的特征标记峰。M+2和M+8之间的差异可以作为显示异质性的指标，用于展示葡萄糖对细胞质中acetyl-CoA的相对贡献。因此，¹³C-SpaceM能够检测ACLY敲降细胞中的异质性，并识别不同的亚群体。这种单细胞和空间异质性无法通过整体分析揭示，显示了¹³C-SpaceM方法的独特优势。

图6. 细胞ACLY敲降后acetyl-CoA的同位素标记程度分析（点击查看大图）

肿瘤组学中氨基酸合成异质性的空间组学分析

研究人员分析了从横向植入表达突变型异柠檬酸脱氢酶（IDH）和红色荧光蛋白（RFP）的GL261胶质瘤细胞的小鼠大脑组织切片。在采集组织前的48小时，小鼠被喂食未标记的或含有U-13C葡萄糖的液体饮食。

首先，研究人员分析了¹²C-葡萄糖饮食的肿瘤携带小鼠大脑切片中的酯化脂肪酸组成。通过比较质谱TIC与显微镜明场和荧光成像，发现整个大脑（包括肿瘤区域）的质谱离子响应很高（图7a）。测试过程中，肿瘤区域与组织切片的其余部分分别采用10μm和50μm激光分辨率进行分析。对不同脂肪酸的空间分析揭示了在非肿瘤携带的脑半球组织中，脂肪酸丰度存在高度的异质性，我们可以仅根据它们的脂肪酸组成来识别的某些结构，如胼胝体和前连合部，这两个区域都富含油酸（18:1）且棕榈酸（16:0）、硬脂酸（18:0）和花生四烯酸（20:4）的含量低。有趣的是，尽管棕榈酸、油酸、硬脂酸和花生四烯酸在肿瘤和周围的大脑组织中的含量相似，肉豆蔻酸（14:0）和棕榈酸（16:1）在肿瘤组织中则明显增加。与大脑其它部分相比，肿瘤中必需脂肪酸亚麻油酸（18:2）和α/γ亚麻酸（18:3）也明显增高。

之后，研究人员分析了喂食含有U-¹³C葡萄糖饮食的小鼠肿瘤组织，从肿瘤组织中选择性分离出的5种主要从头合成的脂肪酸的同位素分布（图7c）。三种饱和脂肪酸肉豆蔻酸（14:0）、棕榈酸（16:0）和硬脂酸（18:0）的13C摄入丰度较高，同位素分布最大分别可至M+10，M+12和M+14。其中，肉豆蔻酸M+0的强度极低，几乎完全源自脂肪酸从头合成。由于肉豆蔻酸对一些重要信号蛋白的翻译后修饰很重要，这一发现表明胶质瘤可能选择性地上调肉豆蔻酸的合成以促进自身生长。相比之下，两种单不饱和脂肪酸，棕榈酸（16:1）和油酸（18:1）的M+0同位素的相对丰度较高。硬脂酸和油酸的M+2同位素丰度明显增加，表明它们是由未标记的前体（即棕榈酸和棕榈酸）延长形成的。研究人员进一步利用棕榈酸的同位素分布计算acetyl-CoA池中源自葡萄糖的比例，发现肿瘤组织内的该比例同样具有显著的空间异质性（图7d）。

图7. 小鼠脑胶质瘤组织内部脂肪酸代谢空间异质性分析（点击查看大图）