《Frontiers in Immunology》：Integrative multi-omics dissection identifies ACO2, KLF5, and IMP4 as central regulators of the mitochondrial–immune axis in ulcerative colitis
编辑推荐：
溃疡性结肠炎（UC）是一种慢性复发性炎症性肠病，其免疫代谢失调机制尚不明确。本研究针对此问题，整合了孟德尔随机化（MR）、多队列转录组学、机器学习与单细胞RNA测序（scRNA-seq）等技术，系统解析UC的线粒体-免疫调控网络。研究发现并验证了ACO2、KLF5、IMP4（及AGPS）是连接线粒体代谢与免疫调控的核心枢纽基因，为理解UC发病机制提供了新的免疫代谢视角，并提示了潜在的治疗靶点和生物标志物。
想象一下，肠道内壁正经历着一场无声的战争。一边是维持秩序、修复损伤的“守军”细胞，另一边是不断攻击、引发持续炎症的“敌军”免疫细胞。溃疡性结肠炎（Ulcerative colitis, UC）就是这场慢性、复发性战争的写照，其特点是结肠黏膜的持续性炎症和上皮屏障的破坏。尽管近年来生物疗法取得进展，但仍有相当一部分患者应答不佳或不稳定，这凸显了疾病背后机制的复杂性与异质性。越来越多的证据表明，细胞代谢重编程是连接遗传易感性、免疫失调和上皮屏障功能障碍的关键纽带。线粒体功能异常、氧化磷酸化改变、脂质重塑等代谢失调，如何精准调控免疫反应和上皮稳态，从而推动UC的发生发展？其中的核心“调控开关”——那些连接代谢与免疫的关键基因及其细胞机制，仍是亟待解析的黑箱。
为系统回答这些问题，一项发表于《Frontiers in Immunology》的研究采用了创新的整合多组学策略。研究人员没有局限于单一视角，而是构建了一个从基因到细胞、从人群到模型的立体分析框架。他们首先利用大规模全基因组关联研究（Genome-wide association study, GWAS）和表达数量性状位点（Expression quantitative trait loci, eQTL）数据，通过孟德尔随机化（Mendelian randomization, MR）和基于汇总数据的孟德尔随机化（Summary-data–based Mendelian randomization, SMR）来寻找与UC风险存在潜在因果关系的基因。接着，他们整合了多个独立的批量转录组数据集进行交叉验证，并运用机器学习算法筛选出具有诊断潜力的核心枢纽基因。为了揭示这些基因在具体细胞类型中的功能，研究进一步结合了免疫浸润分析、人源UC结肠组织的单细胞RNA测序（Single-cell RNA sequencing, scRNA-seq）数据，以及葡聚糖硫酸钠（Dextran sulfate sodium, DSS）诱导的小鼠结肠炎模型进行验证。基因集富集和网络分析则用于探索这些基因可能参与的代谢和免疫通路。
这项宏大的研究得出了哪些具体结果呢？
3.1 通过整合MR、SMR和转录组学分析鉴定UC相关基因
研究人员首先通过MR分析，初步筛选出889个与UC存在全基因组显著因果关联的基因。进一步的SMR分析结合共定位检验，从中精确定位出242个高可信度的候选基因。与此同时，对三个独立UC转录组数据集的荟萃分析识别出5564个差异表达基因（Differentially expressed genes, DEGs）。将两组基因取交集，最终锁定47个在遗传和转录水平上均与UC强相关的枢纽基因。共表达网络分析显示这些基因聚集在特定的功能模块中，与免疫激活、干扰素信号通路等相关。
3.2 基于机器学习的UC诊断特征识别与验证
为了评估这47个基因的诊断价值并进一步筛选，研究团队构建了包含11种机器学习算法的分类器。结果发现，弹性网络正则化广义线性模型（glmnet）性能最优。通过分析模型的特征重要性，四个基因——ACO2、KLF5、IMP4和AGPS—— consistently 排名靠前，被确定为预测UC的Top特征基因。SHAP（SHapley Additive exPlanations）值分析直观展示了这4个基因对模型预测的贡献。
3.3 因果与表达分析确认ACO2、KLF5、IMP4和AGPS在UC中的功能相关性
进一步的MR因果推断显示，AGPS和KLF5对UC风险有正向因果效应，而ACO2和IMP4呈负向效应。在表达层面，荟萃分析表明这四个基因在UC组织中的平均表达水平均呈下降趋势。SMR位点图显示，这四个基因的GWAS信号与eQTL信号存在清晰的共定位，支持共享的因果变异同时影响基因表达和UC易感性。
3.4 四个枢纽基因的免疫细胞关联与浸润分析
通过多种免疫解卷积算法（如xCell, EPIC）量化结肠黏膜中的免疫细胞组成，研究发现UC黏膜的免疫微环境发生了重塑，例如内皮细胞、髓样树突状细胞和CD8+T细胞富集，而中性粒细胞、调节性T细胞和B细胞减少。相关性分析揭示了四个枢纽基因表达与特定免疫细胞丰度之间的密切联系。例如，ACO2表达与NK细胞、CD8+T细胞、单核/巨噬细胞呈正相关；KLF5与Th2细胞、巨噬细胞相关；IMP4与B细胞、活化的髓样树突状细胞相关；AGPS则与内皮细胞正相关。
3.5 单细胞转录组学分析揭示关键基因的免疫细胞特异性表达与关联
利用人源UC结肠免疫细胞单细胞数据集的分析，t-SNE（t-distributed Stochastic Neighbor Embedding）图谱清晰展示了健康与UC样本中免疫细胞簇的差异。特征图显示，ACO2和KLF5在巨噬细胞和部分T细胞亚群中相对富集。相关性网络分析进一步表明，IMP4、ACO2和AGPS的表达与循环T细胞呈正相关，而与巨噬细胞呈负相关，提示了线粒体基因表达与免疫细胞组成重塑之间的潜在联系。
3.6 差异表达与通路富集分析揭示四个基因的免疫细胞特异性功能改变
对比不同免疫细胞簇中基因的表达水平发现，在UC中，ACO2和IMP4在循环单核细胞、巨噬细胞和浆细胞中表达显著降低，而KLF5和AGPS在活化的CD4+FOS高表达和CD8+固有层T细胞亚群中表达升高。通路富集分析显示，与这些基因显著相关的基因集合富集在“抗原加工与提呈”、“氧化磷酸化”和“糖酵解”等通路上，明确了其在免疫和代谢交叉网络中的核心地位。
3.7 小鼠DSS诱导结肠炎模型中四个基因动态表达与功能关联的单细胞验证
在小鼠DSS结肠炎模型的单细胞数据中，UMAP（Uniform Manifold Approximation and Projection）可视化展示了正常、急性期和慢性期结肠中上皮与免疫细胞的动态变化。条形图分析（基于所有免疫细胞合并）显示，Aco2和Klf5的表达在急性期显著降低，并在慢性期进一步下降；Imp4在急性期下降但在慢性期显著上调；Agps则无显著变化。热图展示了这些基因在不同细胞类型和疾病阶段的特异性表达模式。网络分析和KEGG（Kyoto Encyclopedia of Genes and Genomes）通路富集表明，这些基因主要与AMPK信号通路、白细胞跨内皮迁移、产热等免疫-代谢通路相关。
3.8 DSS诱导小鼠结肠炎模型中枢纽基因表达的验证
通过建立DSS诱导的小鼠结肠炎模型并进行qPCR（Quantitative PCR）验证，实验结果表明，与对照组相比，模型小鼠结肠组织中Aco2、Klf5和Imp4的mRNA水平显著降低，而Agps无显著变化。这与人源数据中的发现一致，进一步在实验动物层面证实了这些基因在结肠炎发生发展中的表达变化。
综上所述，这项整合多组学的研究成功绘制了UC中线粒体代谢与免疫调控交互作用的基因蓝图。研究的主要结论在于，通过系统性的因果推断、转录组验证和细胞定位，鉴定出ACO2、KLF5和IMP4是调控UC“线粒体-免疫轴”的核心枢纽基因。ACO2作为三羧酸循环（Tricarboxylic acid cycle, TCA cycle）关键酶，其下调提示UC中存在线粒体能量代谢障碍；KLF5作为上皮修复和免疫调节的关键转录因子，其表达变化关联上皮完整性与Th17免疫应答；IMP4参与核糖体生物合成，其动态变化可能连接细胞增殖/翻译调控与炎症反应。尽管AGPS在转录水平变化不显著，但其在机器学习模型中的重要性提示其在脂质代谢重塑中可能扮演功能性角色。
讨论部分深入阐释了这些发现的意义。该研究不仅将已知与UC相关的遗传位点与具体的功能基因（如ACO2, KLF5）联系起来，还挖掘了IMP4等此前在UC中未被充分重视的基因，拓宽了对疾病机制的理解。更重要的是，它超越了传统的“免疫”或“代谢”单维度视角，通过多组学数据整合，实证了在UC发病过程中，线粒体代谢通路（如氧化磷酸化）的失调与特定免疫细胞（如巨噬细胞、T细胞）的浸润和激活存在基因层面的紧密耦合。这为“免疫代谢”理论在UC中的具体实施提供了新的分子证据。这些枢纽基因作为潜在的生物标志物，有望用于疾病分层或疗效预测；同时，它们也可能成为干预免疫代谢失衡、恢复黏膜稳态的新颖治疗靶点。当然，研究也指出了其局限性，如主要基于计算和关联分析，未来需要更多的功能实验和更接近人体疾病的模型（如类器官）来验证这些基因的因果机制和治疗潜力。总体而言，这项研究为深入理解UC的复杂发病机制开辟了新途径，并为开发更精准的诊断和治疗策略奠定了重要的科学基础。
