《Frontiers in Systems Biology》：Integrative genome-scale metabolic model of GABAergic neurons reveals metabolic signatures across the mild cognitive impairment – Alzheimer disease continuum
编辑推荐：
摘要
引言：轻度认知障碍（mild cognitive impairment, MCI）代表阿尔茨海默病（Alzheimer’s disease, AD）的前驱阶段，但早期神经元功能障碍背后的代谢机制仍未完全阐明。GABA能神经元（GABAergic neu
摘要
引言：轻度认知障碍（mild cognitive impairment, MCI）代表阿尔茨海默病（Alzheimer’s disease, AD）的前驱阶段，但早期神经元功能障碍背后的代谢机制仍未完全阐明。GABA能神经元（GABAergic neurons）维持兴奋–抑制平衡和网络稳定性，在神经退行性变过程中表现出早期易损性，尽管与其功能障碍相关的代谢改变仍未得到充分表征。
方法：研究人员利用去卷积的海马转录组数据，构建了覆盖MCI–AD连续谱的人GABA能神经元情境特异性基因组规模代谢模型（genome-scale metabolic model, GEM）。通过将转录组去卷积与约束性建模（包括通量平衡分析（flux balance analysis, FBA）和通量变异分析（flux variability analysis, FVA））相结合，研究人员在对照（Control）、早期MCI（early MCI, E-MCI）、进展期MCI（advanced MCI, A-MCI）和AD条件下推断出与疾病分期相关的代谢改变。
结果：分析结果表明，能量代谢、谷氨酸–谷氨酰胺–GABA循环、氧化还原稳态、脂质代谢以及神经元–星形胶质细胞代谢相互作用发生进行性重塑。FVA识别出可行通量范围的反应特异性变化，表明可行代谢解空间发生重塑，而非各通路出现一致性收缩。这些预测的代谢改变伴随着GABA能标志物的转录变化，并与独立的代谢组学观察结果定性一致，支持其生物学合理性。
讨论：总体而言，这项工作提供了一个系统层面的计算框架，将转录组改变与GABA能神经元中预测的代谢重塑联系起来，并生成了关于从MCI向AD进展过程中代谢功能障碍的可实验验证假说。
论文解读：整合基因组规模代谢模型揭示GABA能神经元在MCI–AD连续谱中的代谢特征
一、研究背景与科学问题
轻度认知障碍（mild cognitive impairment, MCI）被认为是阿尔茨海默病（Alzheimer’s disease, AD）的前驱阶段，每年约有12%的MCI患者进展为痴呆，因此该阶段被视为早期诊断和干预的关键窗口。然而，尽管脑脊液tau蛋白和淀粉样蛋白-β等生物标志物研究不断推进，现有方法仍难以捕捉疾病进展前系统层面的代谢改变。近年来，神经环路层面的早期功能障碍受到关注，其中兴奋–抑制（excitatory–inhibitory, E/I）平衡的破坏被认为是早期致病事件。GABA能神经元作为抑制性神经元的主要群体，在维持环路稳定和认知功能中发挥核心作用，并在神经退行性变早期表现出选择性易损性。此外，神经元功能与星形胶质细胞支持密切相关，尤其是通过谷氨酸–谷氨酰胺–GABA循环和代谢底物交换实现。在MCI等病理条件下，星形胶质细胞发生功能重编程，损害谷氨酸摄取并改变GABA代谢，进而破坏E/I平衡。然而，由于实验方法难以在体内直接测量特定神经元亚群的细胞内代谢通量，早期GABA能神经元功能障碍的代谢机制仍不清楚。基因组规模代谢模型（genome-scale metabolic model, GEM）提供了一种整合组学数据并推断条件特异性代谢状态的系统生物学框架，但既往模型多为泛组织或泛神经元模型，缺乏细胞类型特异性。
二、研究内容与结论
研究人员利用去卷积的海马转录组数据，构建了人类GABA能神经元的参考GEM，并沿MCI–AD连续谱生成了对照、早期MCI（E-MCI）、进展期MCI（A-MCI）和AD四种情境特异性模型。通过整合转录组去卷积与约束性建模，包括通量平衡分析（FBA）和通量变异分析（FVA），研究发现疾病进展伴随着能量代谢、谷氨酸–谷氨酰胺–GABA循环、氧化还原稳态、脂质代谢及神经元–星形胶质细胞代谢相互作用的进行性重塑。FBA显示，生物量相关通量从对照至AD逐渐降低，且该排序在生物量函数扰动和全局反应边界敏感性分析中保持不变。FVA则揭示反应特异性和分期依赖的可行通量范围改变，表明代谢解空间发生选择性重塑而非均匀收缩。这些预测与独立的代谢组学观察结果定性一致，支持其生物学合理性。该研究提供了将转录组改变与GABA能神经元代谢表型联系起来的系统级计算框架，并生成了可实验验证的假说。
三、主要关键技术方法
研究人员使用公开可得的RNA-seq数据集GSE115565，通过基因–蛋白质–反应（GPR）关联将表达值映射至Recon3D反应，并采用iMAT算法重建参考GABA能神经元GEM。随后利用GSE28146数据集（包含海马CA1区对照、E-MCI、A-MCI和AD样本，共30例）进行CDSeq无监督去卷积，依据GAD1、GAD2和SLC32A1标志物识别GABA能富集转录组分，并以30个独立抑制性神经元标志物进行验证。通过exp2flux算法将分期特异性转录谱整合至参考模型，生成四种情境模型。采用FBA和FVA进行约束性模拟，使用神经元特异性生物量维持反应作为目标函数，营养约束参考DMEM/F12组成。独立验证采用SEA-AD单核RNA测序数据及Lysaker等人的人类iPSC衍生神经元代谢组学数据进行定性比较。
四、研究结果
（一）GABA能神经元特异性代谢模型的结构与功能表征
重建的参考GEM包含1,548种代谢物、2,203个反应和3,290个独特基因。反应分布显示胞质是代谢整合的中心枢纽，谷氨酰胺、谷氨酸和GABA代谢在胞质汇聚。转运反应占40.4%，主要涉及氨基酸交换、单羧酸转运和线粒体载体系统；交换和需求反应占25.69%，表明网络高度依赖细胞外代谢物交换。
（二）转录组去卷积及GABA能神经元组分的定性评估
CDSeq去卷积选择6个潜在组分，通过独立抑制性神经元标志物验证确认所选组分的GABA能身份。与SEA-AD独立snRNA-seq数据的比较显示，相匹配疾病分期之间的相关性较高（Control–Control r=0.42，E-MCI–E-MCI r=0.52，A-MCI–A-MCI r=0.54，AD–AD r=0.45），非匹配分期相关性较低。
（三）疾病进展中的转录动态与功能改变
GABA能标志物分析显示分期依赖性转录重组，GABBR1、GABBR2、MEF2C等基因表达呈负单调趋势，而NOS1、CALB1、BCL11B等呈正单调趋势。功能富集分析显示，GABA能突触传递、神经递质转运、突触调控等相关生物学过程在AD阶段活性显著降低。
（四）情境化对神经元代谢通量的影响
FBA显示最优生物量通量从Control（0.100 mMgDW
?1 h
?1）降至E-MCI（0.086）、A-MCI（0.043）和AD（0.039）。ATP合成酶、己糖激酶、苹果酸脱氢酶、乳酸脱氢酶等能量代谢相关反应通量降低；谷氨酰胺、谷氨酸和GABA代谢相关反应也出现改变。敏感性分析表明相对排序（Control > E-MCI > A-MCI > AD）在±20%扰动下保持不变。
（五）GABA能神经元中的进行性代谢破坏
代谢通路分析显示糖酵解、神经递质代谢、氧化还原稳态和转运过程呈现反应特异性和分期依赖性重塑。GABBR1/GABBR2下调与谷氨酸/GABA代谢降低一致，MEF2C和NR2F2表达下降与能量通路受损平行，而NOS1和NLGN2上调可能反映代偿性转录反应。
（六）通量变异与代谢灵活性
FVA结果不再支持各通路通量范围的一致性收缩，而是显示反应特异性和分期依赖的代谢灵活性改变。葡萄糖代谢、谷氨酰胺交换、GABA转运、脂质代谢和氧化还原相关过程表现出最明显的通量变异改变，而谷氨酸交换和甘油磷酸胆碱代谢相对稳定。
（七）基于独立代谢组学数据的定性合理性评估
模型预测的谷氨酰胺/谷氨酸代谢、α-酮戊二酸、琥珀酸、天冬氨酸、苹果酸–天冬氨酸穿梭及NAD相关代谢改变与Lysaker等人的实验观察定性一致。谷氨酰胺交换、乳酸转运、铵代谢、氯离子转运、GABA转运、谷氨酰胺转运和丙酮酸转运的FVA幅度在AD中较对照降低65%–88%，显示通路水平的一致性。
五、讨论总结与研究结论
讨论部分指出，该研究的主要概念性进展在于将去卷积转录组谱与约束性代谢建模整合，以推断GABA能神经元中条件特异性代谢表型。模型结构以转运和交换反应为主，与GABA能神经元高代谢需求及对星形胶质细胞代谢支持的高度依赖一致。代谢重塑涉及生物能量缺陷、神经递质循环受损、氧化还原失衡和脂质代谢改变，这些过程相互关联并构成系统层面的进行性破坏。研究还发现代谢功能障碍在MCI早期阶段已可检测，提示代谢改变出现在前驱期并随疾病进展逐渐加重。模型验证采用多层策略，包括独立标志物验证、SEA-AD转录组比较和代谢组学定性对照，支持预测的生物学合理性。局限性包括依赖转录组推断、样本量较小、Recon3D非神经元特异性、培养环境营养约束与体内差异、去卷积平台偏差以及稳态假设等。未来需整合空间转录组学、单细胞代谢组学、蛋白质组学和同位素示踪实验，并在患者来源iPSC模型和纵向队列中验证。
研究结论：该研究通过整合细胞类型特异性转录组信息与基因组规模代谢建模，提供了一个系统级框架，用于理解代谢功能障碍如何促进GABA能神经元在从轻度认知障碍向阿尔茨海默病进展过程中的选择性易损性。预测的进展性反应特异性代谢灵活性改变伴随能量代谢、神经递质循环和神经元–星形胶质细胞代谢相互作用的变化，为候选代谢生物标志物和潜在治疗相关代谢通路提供了假设生成基础，但这些通路需经未来实验验证。
