克罗恩病（Crohn′s disease, CD）是一种慢性复发性胃肠道炎症性疾病，以透壁性炎症、肉芽肿形成为特征，并可导致狭窄、瘘管和穿孔等严重并发症。作为一种典型的复杂疾病，CD的发病涉及遗传易感性、环境触发因素、免疫失调和肠道微生物组改变等多重因素的复杂相互作用。尽管过去二十年通过全基因组关联研究（GWAS）在解析CD的遗传结构方面取得了显著进展，已发现超过200个易感位点，涉及先天免疫、自噬、适应性免疫应答和上皮屏障功能等相关基因，但遗传因素仅能解释疾病遗传力的一小部分。与此同时，代谢重编程被认为是CD肠道炎症的核心特征，短链脂肪酸（SCFAs）代谢异常、乳酸积累等变化在实验模型和患者队列中均有报道。肠道微生物组的失调，包括微生物多样性降低、有益共生菌减少和机会致病菌增多，也是CD发病机制的关键组成部分。然而，大多数现有研究仅针对遗传、代谢或微生物特征进行孤立考察，缺乏将这些分子层面联系起来的整合分析，这阻碍了对CD发病机制的系统理解以及精准治疗策略的开发。
为此，研究人员在《HUMAN MUTATION》上发表了一项综合性多组学研究，旨在通过整合GWAS荟萃分析、配对肠组织非靶向代谢组学分析、基于RNA-seq的组织微生物组表征、机器学习特征筛选、通路水平网络分析以及独立队列的16S rRNA基因测序验证，系统阐明CD的病理生理学机制。
为开展此项研究，作者主要运用了以下几项关键技术方法：首先，对已发表的炎症性肠病（IBD）GWAS数据进行系统性整理和荟萃分析。其次，前瞻性招募了6名活动性CD患者，收集其病变组织和相邻正常外观黏膜的配对活检样本，进行非靶向代谢组学分析。第三，利用公开的RNA-seq数据集（GSE193677），对376个肠组织样本（包括115例严重CD和261例健康对照）进行基于未比对序列的组织微生物组分析。最后，通过机器学习整合多组学数据，并利用独立队列的16S rRNA基因测序和代谢物定量进行验证。
研究结果如下：
3.1. GWAS荟萃分析揭示CD的遗传结构
综合分析共鉴定出24个与CD易感性显著相关的遗传变异，分布于13个染色体区域。其中，NOD2基因座（位于16号染色体）和IL23R基因座（位于1号染色体）显示出最强的关联信号。效应大小分析显示，NOD2变异（如rs2066845）的疾病风险最高（OR=5.2），而IL23Rrs11209026则显示出保护效应（OR=0.26）。功能通路分类表明，适应性免疫通路相关的变异占比最高（27.8%），其次是先天免疫通路（16.7%）。
3.2. 组织代谢组学分析揭示CD中SCFA代谢改变
对6名CD患者的配对肠组织样本进行非靶向代谢组学分析发现，病变组织与对照组织的代谢谱存在明显分离。差异代谢物分析显示，乳酸在CD病变组织中显著升高（1.73倍，p=0.028），其次是琥珀酸和丙酸。代谢通路富集分析显示，糖酵解相关通路（包括瓦博格效应、丙酮酸代谢和糖异生）显著富集，提示炎症CD黏膜中存在增强的有氧糖酵解和代谢重编程。
3.3. 组织微生物组失调以微生物多样性降低为特征
对GSE193677数据集的分析显示，CD组织的微生物α多样性（香农指数）显著低于健康对照。差异丰度分析在物种水平鉴定出24个丰度发生显著改变的微生物类群。在CD组织中，有益共生菌（如Prevotella copri和Faecalibacterium prausnitzii）丰度降低，而机会致病菌（如Escherichia coli和Bacillus cereus）丰度增加。
3.4. 多组学整合揭示遗传、代谢与微生物组间的跨域相互作用
整合遗传变异、代谢物谱和微生物组丰度数据后，主成分分析显示CD患者与健康对照在多组学图谱上出现部分分离。跨组学相关性分析鉴定了49个显著的组学间关系，包括33个正相关和16个负相关，表明宿主遗传背景影响肠道微生物定植模式。
3.5. 机器学习鉴定出具有高判别力的关键多组学生物标志物
随机森林机器学习分析显示，遗传特征对疾病分类的总重要性最高，其次是微生物组特征。构建的多组学随机森林模型在疾病分类中表现出色，曲线下面积（AUC）达0.977。特征重要性排名显示，多基因风险评分（PRS）是最具判别力的单个特征，其次是研究入选年龄、CARD9基因的rs10870077变异。微生物类群中，Limosilactobacillus reuteri、Klebsiella variicola等较为重要。代谢物中，丙酸是唯一进入前15位的重要特征。
3.6. 通路水平分析揭示多组学信号汇聚于免疫和代谢通路
网络分析显示微生物类群间存在协调互作。通路富集分析表明，多组学信号显著富集于先天免疫、适应性免疫、自噬和屏障功能等IBD相关生物过程。疾病特征分析显示，NOD2rs2066845和IL23Rrs11209026等关键遗传变异在CD中风险等位基因频率或表达升高，乳酸水平在CD中显著升高。由PRS、L. reuteri和丙酸构成的多组学生物标志物组合中，遗传特征贡献了最大的判别力。
3.7. 16S rRNA测序验证微生物组-代谢物轴
在独立队列中通过16S rRNA基因测序验证了关键微生物类群的改变，并分析了其与代谢物的关系。结果证实E. coli和B. cereus等在CD患者中丰度增加。相关分析显示，B. cereus的丰度与乳酸浓度呈强正相关（ρ=0.82, p=0.019）。乳酸比色法检测也证实CD粪便样本中乳酸浓度显著高于对照。这建立了特定微生物类群（特别是B. cereus）与CD中乳酸积累之间的机制联系，表明菌群失调通过增强乳酸生产加剧肠道炎症。
结论与讨论
本项研究通过系统性的多组学整合，全面描绘了CD的遗传架构、黏膜代谢重编程和组织微生物组失调特征，并揭示了这些分子层面之间的相互作用。研究发现，先天和适应性免疫通路基因的遗传变异是CD易感性的核心；病变肠组织表现出明显的糖酵解代谢特征，乳酸积累尤为突出；组织微生物组多样性降低，伴随着共生菌减少和机会致病菌增多。更重要的是，多组学数据整合揭示了宿主遗传背景与微生物定植模式之间的显著关联，并且机器学习模型利用多基因风险评分、特定微生物类群和代谢物构建的生物标志物组合，实现了对CD近乎完美的分类。通路水平分析进一步表明，遗传、代谢和微生物信号共同汇聚于免疫和代谢相关通路。最后，通过独立队列验证，明确了B. cereus等微生物与乳酸积累之间的强相关性，为“微生物-代谢物”轴在CD发病中的作用提供了直接证据。
这项研究的意义在于，它超越了单一组学研究的局限，首次在CD中系统整合了GWAS、配对组织代谢组和黏膜微生物组数据，提供了疾病病理生理学的系统级视图。所鉴定出的关键多组学生物标志物和相互作用的通路，不仅深化了对CD复杂发病机制的理解，也为未来开发基于患者分子特征的精准分型、风险预测和靶向治疗策略奠定了坚实的基础。研究结果突出表明，靶向宿主免疫、能量代谢和微生物稳态的交叉点，可能是未来干预CD的有效途径。
