编辑推荐：
本研究通过整合生物信息学分析，揭示2型糖尿病（T2D）和慢性乙型肝炎病毒（HBV）感染通过一个由156个基因组成的核心“代谢炎症”转录组特征，协同促进胆管癌（CCA）发生发展的分子机制。该研究鉴定了IL6、TNF、AKT1、STAT3和PPARG等核心枢纽基因，构建了具有独立预后价值的“代谢炎症评分”，为理解高危人群的致癌机制、开发新的预后工具及靶向干预策略提供了分子框架。
1. 背景
胆管癌（Cholangiocarcinoma, CCA）是一种起源于胆管上皮的恶性肿瘤，诊断晚、治疗抵抗性强，5年生存率通常低于20%，是肿瘤学领域面临的严峻临床挑战。其病因复杂，除原发性硬化性胆管炎和肝吸虫感染等已知风险因素外，全球发病率的上升与现代流行病学转变密切相关，特别是与2型糖尿病（Type 2 Diabetes, T2D）和慢性乙型肝炎病毒（Hepatitis B Virus, HBV）感染这两种系统性代谢功能障碍和慢性病毒性肝炎有关。流行病学分析表明，T2D使CCA（尤其是肝内胆管癌iCCA）的风险增加1.5至2.0倍，而HBV感染也已被证实与CCA风险升高相关。这两种疾病在肝脏的病理生理过程中交汇，共同促成了“代谢炎症”状态，即由代谢功能障碍驱动的持续性、低度炎症。尽管流行病学证据确凿，但T2D和HBV与CCA关联背后精确且保守的分子机制仍未明确。本研究旨在通过整合多组学数据，解码这一疾病三角共有的致病网络，检验T2D和慢性HBV感染通过一组核心失调的代谢-炎症驱动基因网络促进胆管癌发生的核心假说。
2. 材料与方法
本研究采用整合生物信息学方法，结合多个公共数据库（TCGA, GEO）的数据，探究T2D、HBV和CCA之间的分子关联。分析流程包括数据获取与预处理、差异表达分析、整合跨条件分析、功能富集、蛋白-蛋白互作网络构建、生存分析和验证。研究使用了TCGA-CHOL数据集（45个样本）和GSE107943数据集（163个样本）作为CCA数据，GSE23343（20个样本）作为T2D肝脏代谢特征数据，GSE58208（102个样本）作为HBV肝脏炎症特征数据。数据经过平台特异性预处理和质量控制，差异表达基因（Differentially Expressed Genes, DEGs）的识别采用DESeq2（用于RNA-Seq）和limma（用于微阵列）方法。通过整合四种关键疾病-对照比较的DEG列表，定义了一个高置信度的核心基因集。随后对核心基因集进行了功能富集分析、蛋白-蛋白互作网络分析、生存分析和基于人类蛋白质图谱的蛋白质水平验证。
3. 结果
3.1. 条件特异性转录组景观揭示病因线索
差异表达分析为每种疾病状态建立了独特但重叠的转录谱。在TCGA-CHOL队列中，鉴定出2347个显著失调基因，上调基因包括MMP7、CEACAM6等已确定的CCA标志物。独立的GSE107943数据集分析结果与TCGA数据集高度一致。T2D肝脏特征表现为IL6、TNF等关键炎症调节因子上调。HBV感染的肝组织显示出强烈的干扰素和抗病毒反应特征，同时伴有IL6、TNF等促炎细胞因子上调。
3.2. 鉴定出保守的核心代谢炎症基因集
整合分析揭示了CCA、T2D和HBV三种病理状态之间共有的156个基因（92个上调，64个下调）核心特征，定义为核心代谢炎症特征。该核心集在疾病状态与正常组织间能可靠区分，并揭示了CCA与T2D、HBV炎症成分之间的部分转录重叠。
3.3. 富集通路凸显代谢-炎症交互作用
核心基因集的功能富集分析揭示了在代谢和炎症界面显著富集的通路。最显著富集的通路是PPAR信号通路，其次是细胞因子-细胞因子受体相互作用、PI3K-Akt信号和TNF信号通路。基因本体分析也证实了这些发现，显示“炎症反应”、“脂质代谢过程”等生物过程以及“细胞因子活性”、“转录因子结合”等分子功能显著富集。
3.4. 网络分析识别中枢枢纽基因与功能模块
对核心基因的蛋白-蛋白互作网络分析揭示了一个包含142个节点和458个相互作用的高置信度网络。多指标中心性分析确定了十个高置信度枢纽基因，其中促炎细胞因子IL6和TNF成为拓扑中心性最高的节点，其次是关键信号转导因子AKT1和STAT3，以及代谢核受体PPARG。通过MCODE算法在网络内识别出三个紧密互连的功能簇：1）炎症信号模块，以IL6、TNF、IL1B和NFKB1为中心；2）代谢调节模块，以PPARG、SREBF1和LEP为中心；3）细胞生长与存活模块，以AKT1、STAT3和MTOR为中心。JUN和NFKB1等连接蛋白被确定位于模块界面，提示了炎症、代谢和增殖信号之间串扰的分子机制。
3.5. 枢纽基因具有预后价值并与临床侵袭性相关
在TCGA-CHOL队列中的生存分析显示，枢纽基因表达水平与患者预后显著相关。促炎枢纽基因IL6和TNF以及增殖枢纽基因STAT3的高表达与较差的总生存期相关，而代谢调节因子PPARG的高表达则与良好的预后相关。此外，IL6和TNF的高表达与晚期肿瘤分期和淋巴结转移显著相关，PPARG的表达与淋巴结转移呈显著负相关。
3.6. 衍生代谢炎症评分是稳健的预后生物标志物
基于五个关键枢纽基因（IL6、TNF、AKT1、STAT3、PPARG）的表达，构建了一个定量代谢炎症评分。在TCGA-CHOL队列中，该评分能有效将患者分层为低、中、高风险组，高风险组患者的总生存期显著更短。在多变量Cox回归分析中，该评分仍是总生存期的独立预测因子。该评分的预后有效性在独立的GSE107943队列中得到了进一步确认。
3.7. 蛋白质水平验证确认转录组失调
使用人类蛋白质图谱进行的免疫组化验证证实了关键枢纽蛋白在胆管癌组织中的失调。与正常胆管上皮相比，CCA样本中IL6和TNF蛋白在肿瘤细胞质和肿瘤相关间质中强表达。AKT1在肿瘤细胞中呈现强烈的胞质和膜染色。PPARG在CCA中表现出显著的核染色丢失，与转录组下调一致，支持了其保护功能丧失的假说。
4. 讨论
本研究首次在代谢炎症概念框架内，通过整合多数据库分析，阐明了T2D、HBV感染和CCA之间的分子互连性。主要发现是鉴定出一个在三种疾病状态下持续失调的、由156个基因组成的保守转录特征。该特征的功能模块化揭示了核心生物程序的协调扰动，这些程序被组织成以炎症信号级联、代谢调节和细胞生长通路为中心的功能模块。
IL6、TNF、AKT1、STAT3和PPARG被确定为顶级枢纽基因，这为理解CCA发病机制提供了关键的机制性见解。这些分子在炎症和代谢子网络中的网络中心性表明，它们充当了将代谢压力转化为致癌信号的分子整合器。观察到的相互失调，特别是促炎介质（IL6、TNF）的激活与关键代谢调节因子（PPARG）的抑制同时发生，为癌症相关恶病质的临床表型和晚期恶性肿瘤中观察到的代谢重编程提供了分子关联。
代谢炎症特征的预后意义证实了其转化相关性。该特征与患者生存期的显著关联，独立于常规临床病理因素，表明对代谢-炎症轴的分子谱分析为标准分期系统增加了预后价值。这一发现在临床中具有潜在效用，特别是对于改善早期CCA患者的风险分层，可能识别出能从强化监测或辅助治疗干预中获益更大的亚组患者。
本研究的发现将T2D和HBV合并症背景下的CCA定性为一种代谢炎症性恶性肿瘤。PPAR（代谢）和细胞因子（炎症）通路在同一基因集中的同时富集，暗示了一个恶性循环：炎症通过下调PPARG等方式抑制代谢稳态，而代谢功能障碍（如脂毒性、胰岛素抵抗）则使炎症信号持续存在。由不同但相互关联的炎症、代谢和生长相关模块组成的模块化网络架构，为这种串扰提供了结构蓝图。IL6和AKT1等枢纽基因位于这些模块的界面，充当分子整合器。在功能模块边界识别的JUN和NFKB1等连接蛋白，指出了维持代谢炎症状态的炎症、代谢和增殖信号之间串扰的关键分子机制。该模型阐释了全身性疾病如何建立一个促进胆道转化的许可性肝脏微环境。
5. 结论
总之，这项系统生物学方法将代谢炎症定义为连接T2D、HBV和CCA的关键机制环节。我们鉴定了一个保守的转录特征及其核心调节因子（IL6、TNF、AKT1、STAT3和PPARG），它们共同协调着一个驱动肿瘤发生的代谢-炎症串扰网络。尽管源自胆管和全肝数据集的组合，但该特征为理解系统性代谢和病毒性疾病如何在原位创造有利于胆道癌变的微环境提供了一个组织水平的框架。此外，衍生的代谢炎症评分代表了一个稳健且独立的预后生物标志物。总而言之，这些发现增进了我们对CCA病因学的理解，为在风险人群中制定新的预防策略奠定了基础，并揭示了包括药物再利用机会在内的可操作治疗靶点。通过严格的生物标志物验证和靶向治疗试验，将这些见解转化为临床实践，有望改善罹患这种毁灭性恶性肿瘤患者的预后。
