作者及机构
本项研究由来自Broad Institute of MIT and Harvard、哈佛大学陈曾熙公共卫生学院、纽约市立大学等机构的跨学科团队共同完成,通讯作者为Curtis Huttenhower。研究于2024年4月发表在*Molecular Systems Biology*期刊(Volume 20, Issue 4)。
炎症性肠病(IBD,inflammatory bowel disease)作为一种复杂的肠道微生态相关疾病,其病理机制与肠道微生物代谢密切相关。尽管近年来微生物组研究取得了重要进展,但人类粪便代谢组中仍有超过90%的代谢特征(metabolic features)未被注释,这些”化学暗物质”可能包含关键的免疫调节分子。传统代谢组学分析面临两大挑战:1) 高分辨率质谱检测到的数万个代谢特征中仅有约0.7%能与标准品匹配;2) 即使通过质谱数据库比对获得”推定质量匹配”(putative mass-match)的代谢物,其注释准确性也难以保证。
研究团队由此开发了MACARRON(Metabolome Analysis and Combined Annotation Ranks to Prioritize Novel Bioactives)计算方法,旨在从复杂的微生物群落代谢组中优先筛选具有潜在生物活性的分子。该方法整合了化学特征、生态学特性和流行病学关联三个维度的信息,在人类微生物组计划二期(iHMP/HMP2)的546份粪便代谢组数据中进行验证,重点关注IBD相关代谢物的发现。
研究分析了iHMP项目中106名参与者(80名IBD患者和26名非IBD对照)的546份纵向粪便代谢组数据(平均每人5.2个时间点),共包含81,867个代谢特征。通过三步过滤:
- 去除14,634个可能为加合物或碎片的冗余特征
- 筛选在至少一种表型(克罗恩病CD、溃疡性结肠炎UC或非IBD)中检出率≥70%的37,201个”普遍特征”
- 将特征分为三类:596个标准品(内部600种化合物库匹配)、41,305个推定质量匹配(与HMDB数据库相似质量)和16%完全未知特征
采用创新性的分步聚类策略:
1) 相关性计算:使用双权中相关(biweight midcorrelation)分析不同表型分层下的代谢物丰度相关性
2) 模块识别:基于WGCNA算法生成355个共丰度模块(大小33-1,116个特征),其中124个为”注释模块”(含至少1个标准品)
3) 化学同质性验证:46%的注释模块显示≥75%的化学类别一致性(如同一模块中标准品均属胆汁酸类)
MACARRON通过三级评分系统对代谢特征进行排序:
1) 丰度相对锚定值(AVA):计算未知特征与同模块标准品的丰度比值
2) 表型关联:混合效应模型评估代谢物在CD活动期/UC vs 非IBD的差异丰度(q<0.05)
3) 综合优先级评分:整合AVA、效应量和q值的百分位数生成0-1的meta-rank
选取优先级评分最高的烟酰胺核糖苷(nicotinamide riboside, NR)进行体内验证:
- 急性DSS模型:3%葡聚糖硫酸钠处理5天后,NR治疗组(1000 mg/kg)结肠炎组织学评分显著改善(p=0.014)
- 慢性DSS模型:三周期梯度浓度DSS诱导后,NR组显示出更强的保护效果(p=0.007)
研究发现:
- 12,220个代谢特征在IBD中显著扰动(q<0.05),其中200个为标准品(如初级胆汁酸增多、次级胆汁酸减少)
- 通过共丰度模块将15,481个未知特征与已知代谢物关联,覆盖33个化学类别,占粪便代谢组的23%
- 鉴定出1,000+个潜在生物活性分子,包括2672个与IBD表型强关联的未知特征
除已知的胆汁酸、短链脂肪酸外,系统揭示了:
- 胆红素衍生物:15个与粪胆素(stercobilin)共变且显著耗竭的特征,质荷比差异提示可能存在新型微生物修饰(如氧基/二氧基加成)
- 维生素代谢紊乱:维生素B3(烟酸)和B5(泛酸)途径显著下调,伴随NAD+合成酶基因表达降低
- 微生物特异分子:对羟基苯乙酸(p-hydroxyphenylacetate)等微生物代谢产物在IBD中异常富集
技术层面:MACARRON首次实现了:
医学层面:
研究者指出,未来可将MACARRON与:
1) 代谢反应预测工具(如Biotransformer)结合,提高未知物注释精度
2) 宏基因组数据联用,区分宿主与微生物来源分子
3) 扩展至其他微生态相关疾病(如代谢综合征、肿瘤免疫治疗)的研究
该项工作不仅为IBD机制研究提供了新视角,更建立了从海量代谢组数据中挖掘功能分子的通用框架,推动了微生物组研究从相关性分析向因果机制发现的跨越。