PAH 暴露的代谢组学指纹——在动态培养 3D 细胞球中识别毒理学生物标志物
📋 基本信息
- 英文标题: Metabolomic fingerprints of PAH exposure - identifying toxicological biomarkers in dynamically cultured 3D cell spheroids
- 作者: Anja Bosnjakovic¹, Thomas O. Eichmann², Alja Stern³,⁴, Morten Arendt Rasmussen⁵,⁶, Mario Lovric¹,⁵,⁷,⁸, a *, Bojana Zegura³,⁴, a *
- 机构:
- ¹ 克罗地亚应用生物人类学中心,萨格勒布
- ² 格拉茨医科大学核心设施质谱中心
- ³ 斯洛文尼亚国家生物研究所遗传毒理学与癌症生物学系
- ⁴ 卢布尔雅那大学生物技术学院
- ⁵ 哥本哈根大学食品科学系
- ⁶ 哥本哈根大学儿童哮喘前瞻性研究
- ⁷ The Lisbon Council,布鲁塞尔
- ⁸ 奥西耶克约西普·尤拉伊·斯特罗斯梅尔大学食品技术学院
- 期刊: bioRxiv preprint
- 发表日期: 2025-07-16
- DOI: 10.1101/2025.07.10.663939
- 许可证: CC BY 4.0
- 资助: Horizon Europe (EDIAQI #101057499; CutCancer #10050740), 克罗地亚科学基金会 (MOBODL-2023-12-5502, MOBDOK-2023-1077), 斯洛文尼亚研究署 (P1-0245)
🔗 原文链接
- PDF 原文: SeaFile 下载
- bioRxiv: https://www.biorxiv.org/content/10.1101/2025.07.10.663939v2
- 补充数据: https://doi.org/10.5281/zenodo.15851269
📖 摘要
环境暴露于多环芳烃(PAHs)会导致代谢功能障碍,但仍需要可靠的生物标志物来评估其对人类健康的影响。本研究使用成熟 21 天的人类 HepG2 球体(一种具有代谢能力的三维 3D 肝脏模型),评估了对分级、非细胞毒性浓度的苯并a芘(BaP)和苯并b荧蒽(BBF)在 24 小时和 96 小时暴露后的代谢反应。非靶向液相色谱-质谱(LC–MS)代谢组学结合多变量和网络分析,鉴定了化合物和时间特异性的代谢特征。在 24 小时时,没有代谢物显示显著变化。相比之下,在 96 小时时,两种 PAHs 一致地改变了 7 种稳健的代谢物,这些代谢物与多胺代谢、膜动力学、线粒体能量和 DNA 修复途径相关。网络分析显示 BBF 引起的变化比 BaP 更广泛和更连贯,表明不同的毒代动力学。这些发现强调了延长暴露在揭示代谢紊乱中的重要性,并支持一组候选生物标志物用于未来的低剂量研究和改善空气传播毒物的风险评估。
🧪 研究设计
3D 球体培养
使用 CelVivo ClinoStar 2™ 旋转生物反应器系统培养 HepG2 3D 球体:

培养流程:
- 使用 AggreWell™ 400 μm 24 孔板预形成细胞聚集体(200,000 细胞/孔)
- 离心 100g × 5 分钟,过夜培养
- 转移至 ClinoReactor™,在 ClinoStar 2™ 系统中动态培养
- 初始转速 10 rpm,随球体生长逐步调整
- 每 48/48/72 小时更换 90% 培养基
- 第 10 天转移至新的 ClinoReactor™
- 培养 21 天后进行 PAH 暴露实验
暴露条件
| 化合物 | 浓度(24h) | 浓度(96h) | 溶剂对照 |
|---|---|---|---|
| BaP | 20 μM / 2 μM | 10 μM / 1 μM | 0.2% / 0.1% DMSO |
| BBF | 20 μM / 2 μM | 10 μM / 1 μM | 0.05% / 0.025% DMSO |
每个实验组设 3 个独立重复池(Exp. 1.1–4.3)。
分析方法
- 细胞活力: CellTiter-Glo® 2.0 ATP 检测
- 代谢组学: 非靶向 LC-MS(QExactive Focus,正负离子 DDA 模式)
- 统计分析: ANOVA + Benjamini-Hochberg FDR 校正
- 网络分析: Spearman 相关 + NetworkX 社区检测
- PCA: z-score 标准化后降维
- 通路富集: MetaboAnalyst 6.0 ORA 模块
📊 核心发现
1. 细胞活力无显著毒性

- 24 小时:所有条件细胞活力 > 90%(除 BBF 20 μM ≈ 80%)
- 96 小时:无细胞毒性,活力与对照相当或略高
- 所有条件均满足 ≥ 75% 活力阈值,适合后续代谢组学分析
2. PCA 分析揭示化合物特异性代谢响应

高浓度 96 小时处理样本的主成分分析:
- PC1 区分代谢响应强度(值越高 = 响应越强)
- PC2 区分响应类型(BBF → 正值,BaP → 负值)
- BBF 高浓度(红色向量):聚类在右上象限,响应最强且一致
- BaP 高浓度(蓝色向量):聚类在右下象限,响应强但模式不同
- 对照组(灰色向量):位于左侧,代谢变化最小
- 低浓度组(绿色/橙色):位于左侧,效应较弱或更多变
25 个关键代谢物通过余弦相似性筛选(≥ 0.95):
- 12 个与 BBF_H profile 最相似
- 13 个与 BaP_H profile 最相似
3. 火山图显示时间依赖性代谢变化

| 处理组 | 显著变化代谢物数量 |
|---|---|
| BBF 24h | 0 |
| BBF 96h | 20 (FDR p < 0.05) |
| BaP 24h | 2 |
| BaP 96h | 24 |
关键发现:24 小时暴露几乎无显著代谢变化,96 小时才出现明显的代谢重编程。
4. 七种稳健代谢物生物标志物
| 代谢物 | PAH | 浓度 | 时间 | 变化方向 | 生物学意义 |
|---|---|---|---|---|---|
| 2-氨基己二酸 | BBF | 10 μM, 1 μM | 96h | ↓ 下调 | 赖氨酸分解代谢,TCA 循环连接 |
| N-乙酰亚精胺 | BBF, BaP | 10 μM, 1 μM | 96h | ↑ 上调 | 多胺代谢,细胞生长/死亡调控 |
| 甘油磷酸肌醇 | BBF | 10 μM, 1 μM | 96h | ↓ 下调 | 膜脂质重塑,GPI 锚定蛋白 |
| 尿苷 | BBF | 10 μM, 1 μM | 96h | ↓ 下调 | 嘧啶核苷酸合成 |
| 戊酰肉碱 (C5) | BBF | 10 μM, 1 μM | 96h | ↓ 下调 | 线粒体 β-氧化,脂肪酸转运 |
| S-腺苷甲硫氨酸 | BBF | 10 μM, 1 μM | 96h | ↓ 下调 | 甲基供体,表观遗传调控 |
| β-柠檬酰谷氨酸 | BBF | 10 μM, 1 μM | 96h | ↑ 上调 | 铁载体,线粒体保护 |
5. 代谢物相关网络

网络特征:
- 节点 = 代谢物,边 = 高置信度相关(|r| > 0.6)
- 深橙色 = 稳健响应代谢物(PCA 对齐 + 统计显著)
- 浅橙色 = PCA 对齐代谢物
- 中心节点提示潜在的调控或整合作用
6. 通路富集分析

BBF 暴露组(蓝色):
- 甲基组氨酸代谢(最强信号)
- 磷脂酰胆碱生物合成
- 超长链脂肪酸 β-氧化
- 亚精胺和精胺生物合成
BaP 暴露组(红色):
- 肉碱合成
- 嘌呤代谢
- 信号较弱,通路更特异
7. BBF 通路网络

BBF 组 96 小时的通路网络显示:
- 氨基酸代谢和甲基化相关通路形成连接簇
- 涉及甲硫氨酸代谢、甜菜碱代谢、组氨酸代谢、甘氨酸和丝氨酸代谢
- 暗示协调的代谢响应,而非孤立的通路命中
- S-腺苷甲硫氨酸 (SAM) 在甲基化和多胺生物合成中起核心作用
🔬 技术创新
CelVivo ClinoStar 系统的优势
- 动态培养:模拟体内血流,提供更好的营养和氧气分布
- 长期培养:支持 21 天成熟期,确保代谢能力完全发育
- 均匀球体:通过"园艺"步骤去除不规则球体,提高培养一致性
- 生理相关性:3D 结构更接近体内肝脏微环境
- 代谢能力增强:相比 2D 培养,异生物质代谢酶活性显著提升
代谢组学策略
- 非靶向 LC-MS:无偏倚检测所有代谢物变化
- 多变量分析:识别整体代谢模式
- 网络分析:揭示代谢物之间的关联和通路扰动
- PCA + 余弦相似性:筛选与处理模式最匹配的代谢物
💡 核心创新点
1. 延长暴露时间的重要性
首次证明 96 小时暴露才能检测到显著代谢变化,24 小时可能不足以揭示 PAHs 的代谢毒性。这对传统体外毒性试验设计(通常 24-48 小时)提出挑战。
2. BBF 的毒性被低估
尽管 BaP 是 IARC 1 类致癌物,是监管框架的参考化合物,但本研究发现 BBF 引起的代谢扰动更广泛、更一致,提示 BBF 的毒性可能被低估。
3. 3D 球体模型优势
使用成熟的 3D 肝脏模型评估环境毒物,比传统 2D 模型更具生理相关性,特别是对于需要代谢激活的化合物。
4. 候选生物标志物
鉴定了 7 种稳健的代谢物作为 PAHs 暴露的潜在生物标志物,可用于:
- 低剂量 PAHs 暴露研究
- 改善空气传播毒物的风险评估
- 整合到 ILCR(终生癌症风险)评估中
🎯 代谢物详细解读
N-乙酰亚精胺(↑ 上调)
- 多胺代谢中间体,调控细胞生长、分化和死亡
- 在 BaP 和 BBF 暴露后均上调,提示多胺代谢紊乱是 PAHs 暴露的共同响应
- 可能反映多胺周转受损,指示细胞应激和稳态失衡
甘油磷酸肌醇(↓ 下调)
- 肌醇脂质的降解产物,参与膜生物合成和细胞信号传导
- BBF 暴露后下调,提示膜脂质重塑或 GPI 锚定蛋白降解增加
戊酰肉碱 C5(↓ 下调)
- 短链酰基肉碱,参与支链氨基酸分解代谢和线粒体 β-氧化
- 下调提示脂肪酸转运和线粒体能量产生受损
β-柠檬酰谷氨酸(↑ 上调)
- 内源性铁载体,可向线粒体乌头酸酶递送 Fe(II)
- 上调可能是细胞适应性响应,维持线粒体功能并减轻氧化损伤
S-腺苷甲硫氨酸 SAM(↓ 下调)
- 主要甲基供体,参与 DNA/RNA/蛋白质/脂质甲基化
- 下调提示甲基化能力受损,可能影响表观遗传调控
尿苷(↓ 下调)
- 嘧啶核苷酸合成中间体
- 下调反映线粒体功能障碍和氧化应激抑制嘧啶生物合成
2-氨基己二酸(↓ 下调)
- 赖氨酸分解代谢产物,连接 TCA 循环
- 下调提示赖氨酸分解代谢抑制或线粒体功能障碍
⚠️ 局限性
- preprint 状态:尚未经过同行评审
- 单一细胞系:仅使用 HepG2(肝癌细胞系),缺乏原代肝细胞的复杂性
- 体外模型:无法完全模拟体内肝脏的代谢和解毒过程
- 仅测试两种 PAHs:环境中存在多种 PAHs 混合物
- 固定浓度和时间:24h/96h,可能不反映真实世界的慢性低剂量暴露
- FDR 校正后通路不显著:原始 p 值有生物学意义,但统计显著性不足
- 相关分析排除部分代谢物:可能遗漏生物学相关信息
🎯 学术启示
对毒理学研究的影响
- 暴露时间设计:未来体外毒性研究应考虑更长的暴露时间(≥96 小时)
- 3D 模型推广:3D 球体模型可能成为环境毒物评估的标准工具
- 代谢组学应用:非靶向代谢组学是识别环境暴露生物标志物的有效策略
对风险评估的意义
- BBF 毒性重估:BBF 的代谢扰动比 BaP 更广泛,可能需要重新评估其风险等级
- 低剂量研究:鉴定的候选生物标志物可用于低剂量 PAHs 暴露研究
- ILCR 整合:代谢生物标志物可能整合到终生癌症风险评估中
对 Celvivo 产品的相关性
- ClinoStar 系统验证:本研究成功使用 ClinoStar 系统进行 21 天长周期培养,证明其在毒理学研究中的实用性
- 肝脏毒性评估:可作为肝脏毒性评估的体外模型验证案例
- 代谢组学整合:展示了 3D 培养系统与代谢组学分析的结合应用
- 差异化优势:动态培养相比静态培养能更好地维持细胞代谢功能
📚 关键词
代谢组学 · 多环芳烃 · 3D 细胞球 · 肝脏模型 · 毒理学生物标志物 · 苯并a芘 · 苯并b荧蒽 · LC-MS · 动态培养 · 风险评估 · Celvivo · ClinoStar
联系人|李东伟
📞 18045130660
🌐 主页:https://www.biokingdom.top