转录组指导的代谢建模揭示轻度认知障碍与阿尔茨海默病进展中星形胶质细胞特异性脆弱性

《Frontiers in Bioinformatics》：Transcriptome-informed metabolic modeling reveals astrocyte-specific vulnerabilities in mild cognitive impairment and Alzheimer’s disease progression

摘要

引言：星形胶质细胞对于维持神经元稳态至关重要，但其在轻度认知障碍（MCI）和阿尔茨海默病（AD）不同疾病阶段中的特异性贡献仍然认识不足。本研究旨在采用整合转录组学与基因组尺度代谢建模方法，研究从对照到MCI再到AD连续谱中的星形胶质细胞相关转录与代谢改变。

方法：研究人员分析了海马CA1组织（GSE28146）在4种临床状态下的转录组谱：对照、早期MCI、晚期MCI和AD。采用无监督去卷积算法CDSeq推断星形胶质细胞相关表达程序，并通过经典标志物富集、与外部参考特征的相关性分析以及与独立单核RNA测序（RNA-seq）星形胶质细胞伪总体数据集的比较进行验证。随后将推断得到的表达谱整合到经人工校正的人类星形胶质细胞基因组尺度代谢模型中，构建条件特异性模型，并使用通量平衡分析（FBA）和通量变异分析（FVA）进行分析。

结果：分析结果支持在疾病进展过程中，星形胶质细胞相关转录程序和代谢程序发生阶段依赖性重塑。早期MCI主要与信号传导和应激适应性改变相关，而晚期MCI和AD则表现出更广泛的突触支持、氧化还原稳态以及炎症相关程序紊乱。模型预测显示，从对照到晚期MCI，源于生物量的维持代理指标逐步下降，而在AD阶段随后出现部分回升，这提示其更可能反映向反应性样星形胶质细胞状态的代偿性转变，而非功能的完全恢复。通量变异分析显示，疾病各阶段的代谢灵活性均有所降低，尤其体现在谷氨酸-谷氨酰胺循环、谷胱甘肽/氧化还原代谢、糖酵解-丙酮酸代谢、胆固醇处理以及一碳/叶酸代谢等方面。

讨论：这些发现支持如下观点：在从健康衰老向AD转变的过程中，星形胶质细胞会发生进行性且具有阶段特异性的代谢重编程。氧化还原调控、神经递质循环和线粒体功能方面的早期改变，可能促成神经元早期脆弱性。本研究强调了以星形胶质细胞为中心的通路可作为未来实验验证和治疗探索的潜在靶点。

该文发表于《Frontiers in Bioinformatics》，围绕阿尔茨海默病（AD）连续病程中星形胶质细胞（astrocyte，承担中枢神经系统代谢支持、神经递质回收与氧化还原稳态维持的重要胶质细胞）的阶段特异性功能失衡展开。研究背景在于，既往神经退行性疾病研究长期偏重神经元损伤与死亡，但越来越多证据表明，轻度认知障碍（MCI）向AD进展并非单一神经元病变，而是涉及神经元、星形胶质细胞、小胶质细胞、少突胶质细胞及血管相关细胞等多细胞网络的动态重塑。其中，星形胶质细胞在谷氨酸-谷氨酰胺循环、能量代谢、胆固醇转运、抗氧化防御和突触支持中居于枢纽位置，因此其早期异常很可能参与疾病起始与放大。然而，目前关于星形胶质细胞在对照、早期MCI、晚期MCI及AD不同阶段中究竟经历何种转录重塑和代谢重编程，仍缺乏系统且细胞类型分辨的机制性解析。传统转录组研究虽揭示了炎症、线粒体功能、脂质代谢及突触通路异常，但多停留于描述层面，难以直接映射到细胞代谢能力和功能状态改变，因此有必要引入基因组尺度代谢模型（GEM）与约束性建模框架，对转录变化进行功能翻译。

为解决上述问题，研究人员构建了一个“转录组到代谢”的整合计算流程，以海马CA1区样本为核心对象，试图在疾病连续谱上提炼星形胶质细胞特异性脆弱环节。研究结论显示，星形胶质细胞并非在AD晚期才发生异常，而是在早期MCI阶段即已出现信号转导、应激适应、氧化还原调节和代谢灵活性下降等变化；随着病程推进，这些异常进一步扩展到突触支持、谷氨酸代谢、谷胱甘肽（GSH）相关抗氧化网络、糖酵解-丙酮酸轴、胆固醇代谢和一碳/叶酸代谢等核心功能模块。整体上，研究支持星形胶质细胞在MCI至AD转变过程中经历进行性、阶段依赖性的代谢重编程，这一过程可能与神经元脆弱性增加密切相关，并为未来靶向星形胶质细胞代谢干预提供了候选通路。

研究人员使用的主要技术方法可概括为以下几类：首先，对GEO数据库中的转录组数据集进行比较筛选，最终选用海马CA1组织微阵列数据集GSE28146，覆盖对照、早期MCI、晚期MCI和AD样本；其次，采用主成分分析（PCA）和线性判别分析（LDA）评估分组可分性，并开展差异表达与功能富集分析；再次，利用无监督去卷积算法CDSeq从混合组织转录信号中推断星形胶质细胞相关表达谱，并结合经典标志物、Allen Institute参考谱和独立单核RNA测序伪总体数据进行多层验证；最后，将去卷积得到的星形胶质细胞表达谱整合到人工校正的人类星形胶质细胞GEM中，借助通量平衡分析（FBA）和通量变异分析（FVA）构建并比较各阶段条件特异性代谢模型。

在研究结果部分，论文首先以“3.1 Transcriptomic differentiation and functional enrichment across disease stages”为题，说明海马CA1区转录组在疾病阶段间呈现连续性重排。PCA显示对照与AD样本分居投影空间两端，早期MCI与晚期MCI位于中间区域，提示MCI并非离散类别，而是具有明确生物学过渡意义的状态。差异基因功能富集进一步揭示分阶段分子特征：早期MCI主要富集于酶联受体信号、VEGFA–VEGFR2信号、膜组织和轴突导向，反映早期应答与适应性改变；晚期MCI转向突触传递、跨突触信号和囊泡介导通讯等通路，提示神经网络支持功能受损；AD阶段则出现突触囊泡循环、胞吞作用、丝裂原活化蛋白激酶（MAPK）信号、氧化应激反应和神经炎症通路的广泛紊乱，显示疾病已进入更深层次的系统失衡。

在“3.2 Astrocyte-specific transcriptomic remodeling”部分，研究人员利用CDSeq从总体转录组中解析出6个潜在细胞成分，包括星形胶质细胞、神经元、少突胶质样细胞、小胶质样细胞、内皮样细胞及一个混合胶质群体。星形胶质细胞成分通过胶质纤维酸性蛋白（GFAP）、水通道蛋白4（AQP4）和乙醛脱氢酶1家族成员L1（ALDH1L1）等经典标志物得到确认。结果显示，正常对照中的星形胶质细胞表达程序偏向突触维持、能量代谢、囊泡信号和抗氧化防御；而在早期与晚期MCI中，转向炎症、细胞骨架重塑和应激反应程序；到AD阶段，神经递质循环和氧化还原稳态相关基因明显下调，提示稳态支持型和神经保护型功能进一步衰退。

在“3.3 External validation of astrocyte-specific deconvolution”部分，研究人员使用独立的内嗅皮层单核RNA测序星形胶质细胞数据进行外部验证。尽管该外部数据的分期体系与GSE28146基于简易精神状态检查量表（MMSE）的分组不完全一致，但相关性分析显示，对应疾病阶段之间存在较高的一致性趋势；同时，在通路层面观察到共同的功能主题：早期阶段涉及细胞外重塑、缺氧反应和早期应激适应，中间阶段突出蛋白稳态与蛋白折叠，晚期阶段则呈现信号、激素相关过程及钙介导通路异常。作者据此认为，尽管基因层面的逐一对应有限，但功能模块层面的收敛支持去卷积所得星形胶质细胞程序具有生物学可信度。

在“3.4 Context-specific astrocyte metabolic models”部分，研究人员将星形胶质细胞特异性转录谱整合进人工校正的人类星形胶质细胞GEM，构建了对照、早期MCI、晚期MCI和AD四个条件特异性模型。模型规模覆盖6520种代谢物和10586个反应，并使用类似E-Flux的表达约束策略调节反应通量边界。FBA结果显示，以生物量为基础的维持代理指标从对照到晚期MCI持续下降，在AD阶段出现部分回升。论文明确指出，这种回升不应简单解释为功能恢复，更可能代表病理条件下代谢通量重新分配后的代偿状态，与反应性样星形胶质细胞表型相一致。FVA则显示，随着疾病进展，允许的通量范围持续收窄，提示代谢灵活性逐级下降。

在“3.5 Comparative validation of GEM-derived flux changes against APOE-dependent astrocyte metabolic profiling”部分，作者将模型预测与APOE4相关人诱导多能干细胞（iPSC）来源星形胶质细胞同位素代谢实验结果进行对照，以评估模型的生物学相关性。比较结果表明，在谷氨酸代谢模块中，模型预测谷氨酰胺酶通量明显下降；在线粒体谷氨酸代谢层面，NAD依赖和NADP依赖谷氨酸脱氢酶反应均出现通量方向反转，提示氧化还原耦联和线粒体代谢重组；三羧酸循环（TCA）相关反应则整体表现为通量幅度下降，尤其是与琥珀酸相关的代谢。这些结果与实验观察在方向和通路层面呈部分一致，支持谷氨酸代谢、线粒体功能和中心碳代谢重编程是星形胶质细胞异常的重要轴线。

讨论部分系统整合了上述发现。论文指出，MCI至AD病程可被理解为从早期信号失衡与适应性应答，逐步发展为突触传递障碍、氧化还原失衡、脂质处理紊乱和代谢刚性增强的连续过程。特别值得注意的是，代谢灵活性的丧失在早期MCI已可检测，说明星形胶质细胞代谢约束的增加可能先于明显神经变性而出现。谷氨酸-谷氨酰胺循环受损提示神经递质回收支持功能下降，可能增加兴奋性毒性风险；谷胱甘肽与氧化还原相关途径异常提示抗氧化缓冲能力减弱；胆固醇代谢异常则可能影响膜稳态、突触形成与神经元胆固醇供给；一碳/叶酸代谢下降则涉及甲基化能力、核苷酸合成和烟酰胺腺嘌呤二核苷酸磷酸（NADPH）供应不足。作者同时强调，表达约束型FBA本质上预测的是满足约束条件的可行通量分布，而非真实酶活性或细胞内代谢物浓度；此外，去卷积源于混合组织并依赖模型假设，微阵列伪计数转换及转录-蛋白-代谢层间偏差均构成研究局限。因此，本文结果更适合作为机制假设生成和优先靶点筛选依据，而非最终因果定论。

论文结论部分可译为：本研究整合无监督CDSeq与基因组尺度代谢建模，基于海马CA1数据构建了具有疾病阶段分辨率的星形胶质细胞特异性代谢模型，包括对照、早期MCI、晚期MCI和AD，从而实现了将基因表达变化机制性转化为预测代谢后果。总体而言，结果支持这样一种预测：沿MCI–AD连续谱，星形胶质细胞会发生早期且阶段依赖性的代谢重编程，这与疾病相关代谢异常相一致，其特征包括代谢灵活性下降，以及涉及氧化还原稳态（谷胱甘肽）、谷氨酸-谷氨酰胺循环、能量代谢（糖酵解/丙酮酸）、脂质/胆固醇处理和一碳/叶酸代谢等核心神经支持通路的破坏。AD阶段生物量来源维持代理指标表面上的部分回升，更可能反映与反应性样星形胶质细胞状态一致的代偿性转变，而非真实的功能恢复。尽管这些预测仍需进一步实验验证，并有待与更多组学层面整合，但所提出框架为识别星形胶质细胞早期代谢脆弱性、并为未来旨在理解或潜在调控痴呆进展的研究优先排序相关通路和靶点，提供了可复现的策略。