bioRxiv 精读:PAH暴露代谢组学指纹——动态培养3D细胞球
📥 原文PDF下载:点击下载
📋 文献基本信息
| 项目 | 内容 |
|---|---|
| 标题 | Metabolomic fingerprints of PAH exposure – identifying toxicological biomarkers in dynamically cultured 3D cell spheroids |
| 中文 | PAH暴露的代谢组学指纹——在动态培养3D细胞球中识别毒理学生物标志物 |
| 期刊 | bioRxiv preprint(同行评审中) |
| DOI | 10.1101/2025.07.10.663939 |
| 发表日期 | 2025年7月16日(预印本) |
| 作者 | Anja Bosnjakovic, Thomas O. Eichmann, Alja Stern, Morten Arendt Rasmussen, Mario Lovric, Bojana Zegura |
| 机构 | 克罗地亚人类学研究所、奥地利格拉茨医科大学、斯洛文尼亚国家生物研究所、丹麦哥本哈根大学、克罗地亚奥西耶克大学 |
🎯 核心科学问题
多环芳烃(PAH)暴露的生物标志物识别——传统 2D 模型 vs 动态 3D 球体模型
- PAH 是广泛存在的环境污染物(有机物不完全燃烧产生)
- BaP(苯并a芘)是 IARC 1 类致癌物;BBF(苯并b荧蒽)是 IARC 2B 类致癌物
- 传统 2D 培养的肝细胞代谢能力严重不足,无法准确反映体内代谢响应
核心发现:
- 24 小时暴露:没有代谢物显示显著变化
- 96 小时暴露:两种 PAH 均一致改变 7 种稳健代谢物
- BBF 比 BaP 诱导更广泛、更连接的网络变化 → 不同的毒性动力学特征
🔬 技术方案
动态培养 3D 球体
| 参数 | 详情 |
|---|---|
| 设备 | ClinoStar 2™ 旋转生物反应器系统 |
| 细胞 | HepG2 人肝细胞癌细胞系(野生型 TP53) |
| 球体数量 | ~1200 个球体/反应器 |
| 培养体积 | 10 mL |
| 转速 | 10 RPM(随球体生长调整) |
| 培养周期 | 21 天成熟球体(动态培养) |
| 换液方案 | 每 48/48/72 小时更换 90% 培养基 |
PAH 暴露方案
| 化合物 | 高浓度 | 低浓度 | 24h 暴露 | 96h 暴露 |
|---|---|---|---|---|
| BaP(苯并a芘) | 20 μM | 2 μM | ✅ | 10 μM / 1 μM ✅ |
| BBF(苯并b荧蒽) | 20 μM | 2 μM | ✅ | 10 μM / 1 μM ✅ |
检测平台
- 非靶向 LC-MS/MS 代谢组学(QExactive Focus,Thermo Fisher)
- 148 种代谢物被检测
- 统计分析:ANOVA + BH-FDR 校正
- 网络分析:Spearman 相关 |r| > 0.6
- 通路富集:MetaboAnalyst 6.0
📊 核心实验数据
1. 细胞毒性

Figure 2: PAH细胞毒性检测
| 条件 | 细胞活力 vs 对照 | 结论 |
|---|---|---|
| 96h BaP 1μM | > 90% | 无毒性 |
| 96h BaP 10μM | > 90% | 无毒性 |
| 96h BBF 1μM | > 90% | 无毒性 |
| 96h BBF 10μM | > 90% | 无毒性 |
结论:所有用于代谢组学分析的暴露条件细胞活力均 > 75%(预定义阈值),代谢变化非细胞毒性驱动。
2. PCA 分析(96h 高剂量组)

Figure 3: PCA分析 — BBF和BaP代谢响应分离
- PC1 分离代谢响应强度(右侧 = 强响应)
- PC2 区分响应类型(上方 = BBF 型,下方 = BaP 型)
- BBF 高剂量(红色向量)→ 右上象限(PC1+, PC2+),强且一致
- BaP 高剂量(蓝色向量)→ 右下象限(PC1+, PC2-),强但方向不同
3. 7 种稳健代谢物标志物
| 代谢物 | BaP 96h | BBF 96h | 变化方向 | 关联通路 |
|---|---|---|---|---|
| N-乙酰亚精胺 | ↑ | ↑ | 一致上调 | 多胺代谢 |
| 2-氨基己二酸 | - | ↓↓ | BBF 下调 | 氨基酸分解 |
| 甘油磷酸肌醇 | - | ↓↓ | BBF 下调 | 膜动力学 |
| 乳清苷(orotidine) | - | ↓↓ | BBF 下调 | 嘧啶代谢/DNA修复 |
| S-腺苷甲硫氨酸(SAM) | - | ↓↓ | BBF 下调 | 甲基化/表观遗传 |
| 戊酰肉碱(C5) | - | ↓↓ | BBF 下调 | 脂肪酸氧化 |
| β-柠檬酰谷氨酸 | - | ↑↑ | BBF 上调 | 三羧酸循环 |

Figure 4: Volcano图 — 不同PAH暴露组显著变化代谢物
4. 通路富集分析
| 通路 | BaP | BBF | 说明 |
|---|---|---|---|
| 肉碱合成 | ✅ 最高富集 | - | BaP 特异性 |
| 甲基组氨酸代谢 | - | ✅ 最高富集 | BBF 最显著 |
| 磷脂酰胆碱生物合成 | - | ✅ | 膜结构 |
| 超长链脂肪酸 β-氧化 | - | ✅ | 线粒体能量 |
| 亚精胺/精胺生物合成 | - | ✅ | 多胺代谢 |
| 甲硫氨酸代谢 | - | 网络集群 | 甲基化核心 |

Figure 5: 代谢物相关性网络

Figure 6: 通路富集分析
关键发现:BBF 的代谢扰动更广泛、更系统性,涉及氨基酸代谢、甲基化、脂质代谢的多个连通通路。BaP 效应更聚焦(仅 2 条通路),但富集比也很高。
💡 创新点
- 首次将 ClinoStar 2 动态培养系统应用于毒理学代谢组学研究
- 21 天成熟 HepG2 球体 → 代谢能力显著优于 2D 单层培养
- 24h vs 96h 对比:证明延长暴露揭示传统短期检测无法捕捉的代谢干扰
- BBF 毒性机制首次系统刻画:比 BaP 更广泛但方向不同
- 非靶向代谢组学 + 网络分析 + 通路富集 → 无偏发现生物标志物
- 7 种稳健代谢物标志物:为未来低剂量暴露研究提供候选标志物
📈 关键量化数据速查表
| 指标 | 数值 | 意义 |
|---|---|---|
| 球体培养周期 | 21 天 | 成熟球体 |
| 球体数量/反应器 | ~1200 | 高通量 |
| 检测代谢物数 | 148 | 非靶向代谢组学 |
| 24h 显著变化代谢物 | BaP 2 个,BBF 0 个 | 短期无法检测 |
| 96h 显著变化代谢物 | BaP 24 个,BBF 20 个 | 长期暴露揭示效应 |
| 稳健标志物数量 | 7 种 | 双浓度一致变化 |
| 细胞活力(96h) | 均 > 90% | 非毒性驱动 |
| BBF 网络连通度 | 显著高于 BaP | 更广泛扰动 |
📌 局限性(论文自述)
- 仅 HepG2 细胞系:单一肝细胞系,缺乏原代肝细胞验证
- 预印本未经同行评审:仍在评审中
- FDR 校正后通路富集不显著:原始 p 值有意义但未通过多重检验校正
- 无 2D vs 3D 直接对比:虽然声称 3D 代谢能力更强,但无实验数据直接对比
联系人|李东伟
📞 18045130660
🌐 主页:https://www.biokingdom.top