尿液DNA甲基化和羟甲基化动态作为职业化学暴露候选生物标志物——一项探索性初步研究

《International Journal of Molecular Sciences》：Urinary DNA Methylation and Hydroxymethylation Dynamics as Candidate Biomarkers of Occupational Chemical Exposure—An Exploratory Pilot Study

　　

DNA甲基化（5-mdC）及其氧化衍生物反映了环境暴露与基因组调控之间的相互作用，但其对职业因素的短期响应性尚未得到充分表征。这项探索性初步研究评估了标准工作周是否与不同职业环境中的尿液表观遗传和氧化DNA修饰谱的可测量变化相关。研究人员采用配对暴露前/后设计，考察了代表三种暴露场景的队列——基于氰基丙烯酸酯的美容服务、甲苯二异氰酸酯（TDI）暴露和聚氨酯（PUR）泡沫加工。定量了尿液中的5-甲基-2′-脱氧胞苷（5-mdC）、5-羟甲基胞嘧啶（5-hmCyt）、5-羟甲基-2′-脱氧胞苷（5-hmdC）、8-氧代-7,8-二氢-2′-脱氧鸟苷（8-oxodG）、5-羟甲基-2′-脱氧尿苷（5-hmdU）和5-羟甲基尿嘧啶（5-hmUra），并推导出表观遗传周转率（5-hmdC/5-mdC）。尽管观察到的变化在错误发现率（FDR）校正后均未达到正式统计学显著性（p_FDR > 0.05），但各队列中均观察到暴露特异性、方向一致的非显著性趋势。异氰酸酯暴露组表现出5-mdC和5-hmdC的协调性降低，同时伴有5-hmdU升高。PUR加工组表现出最显著的反应，包括5-hmdU显著升高和周转率增加。相比之下，美容行业的中位数变化最小，但个体间变异性较大。这些初步发现表明，尿液表观遗传和氧化DNA标志物可能反映短期暴露相关的生物学趋势，值得进一步评估。

**研究背景、问题与目的**：脱氧核糖核酸（DNA）的化学与结构稳定性是细胞稳态的基础。现代分子生物学认为，从生理健康向病理状态（包括癌变）的转变，是遗传突变、表观遗传失调和氧化损伤协同作用的结果。DNA甲基化（5-mdC）是表观遗传调控的基石，而5-羟甲基-2′-脱氧胞苷（5-hmdC）作为“第六碱基”与转录活性及染色质可塑性相关。环境与职业应激源常触发活性氧（ROS）过量产生，导致鸟嘌呤氧化形成8-氧代-7,8-二氢-2′-脱氧鸟苷（8-oxodG）以及嘧啶氧化产物如5-羟甲基-2′-脱氧尿苷（5-hmdU）和5-羟甲基尿嘧啶（5-hmUra）。目前，短期职业暴露对尿液表观遗传和氧化DNA修饰动态的影响尚未充分阐明。国际劳工组织（ILO）和世界卫生组织（WHO）估计，全球超过10亿工人暴露于有害物质，每年约100万人死于职业相关癌症。因此，需要灵敏、非侵入性的“早期预警”生物标志物，以在临床症状出现前检测生物学压力。本研究旨在评估标准工作周后，三个不同行业工人尿液表观遗传标志物（5-mdC、5-hmCyt、5-hmdC）和氧化损伤产物（8-oxodG、5-hmdU、5-hmUra）的动态变化，探讨这些分子特征是否可作为职业遗传毒性应激的早期预警指标。该研究发表在《International Journal of Molecular Sciences》。

**主要关键技术与方法**：本研究招募了三个职业队列：美容行业（n=13，接触氰基丙烯酸酯）、异氰酸酯暴露组（n=5，有效配对4人）和聚氨酯（PUR）泡沫加工组（n=4，接触聚合物粉尘和三聚氰胺）。采用配对暴露前/后设计，在标准工作周前后收集晨尿。主要技术方法包括：二维超高效液相色谱-串联质谱（2D UPLC-MS/MS）定量除5-hmUra外的修饰核苷；5-hmUra采用液相色谱/气相色谱-同位素稀释质谱（LC/GC-MS）定量。所有数据经肌酐校正。统计方法包括Wilcoxon符号秩检验、Benjamini–Hochberg错误发现率（FDR）校正、效应量（r）计算和主成分分析（PCA）。

**研究结果**：
2.1 不同职业组尿液DNA修饰的定量分析：通过定量分析，异氰酸酯组（组2）表现出最显著的表观遗传改变，5-mdC和5-hmdC的中位数和四分位距（IQR）在暴露后均显著降低。美容行业（组1）以个体间异质性为特征，尤其是基线5-hmCyt水平分布广泛，暴露后分布变窄且中位数小幅下降。PUR泡沫加工组（组3）的5-hmCyt呈相反趋势，暴露后中位数上升。8-oxodG在各组中相对稳定，而嘧啶来源的氧化标志物对特定工业应激源敏感：组3的5-hmdU中位数和IQR显著升高，组2也呈平行上升趋势；组3的5-hmUra水平一致升高，组1则保持较高浓度但方差减小。
2.2 工作周诱导的变化：通过计算中位数百分比变化（%Δ），确认了行业特异性分子轨迹。组1（美容）的表观遗传标志物中位数变化微小（5-mdC +0.7%），但个体间变异大，同时5-hmCyt（-23.2%）和5-hmUra（-21.8%）中位数下降。组2（异氰酸酯）表现出表观遗传标志物的同步消耗：5-mdC下降17.2%，5-hmdC下降29.1%，5-hmCyt下降22.8%，而5-hmdU中位数升高103.3%。组3（PUR加工）呈现最显著反应：5-hmdU中位数升高221.5%，5-hmCyt升高108.7%，5-hmdC升高14.1%，但5-mdC仅轻度升高（+3.75%），8-oxodG下降14.9%。
2.3 表观遗传周转率分析：通过计算5-hmdC与5-mdC的摩尔比（[5-hmdC]/[5-mdC] × 10^-3）作为表观遗传周转率，发现组2（异氰酸酯）的周转率下降16.3%，表明TET介导的氧化通路线受抑制，导致表观遗传停滞。组3（PUR加工）的周转率升高133.3%，提示TET通路过度激活，表观遗传高速周转。组1（美容）的周转率保持稳定（+6.5%），表明全球甲基化与去甲基化平衡维持。FDR校正后，所有变化均未达到统计学显著性（p_FDR > 0.05），但效应量分析（r值）显示组2的5-mdC（r=0.913）和组3的5-hmUra（r=0.913）具有非常大的效应量。PCA分析显示，第一主成分（PC1）区分了体尺寸与表观遗传标志物变化，第二主成分（PC2）反映了年龄和生活方式对氧化修饰的影响，但组3的个体年龄较大且吸烟比例高，可能混淆了尿嘧啶衍生物的变化。

**讨论总结与结论翻译**：讨论部分指出，PUR加工组中5-hmdU的221.5%激增是最显著发现，但需谨慎解释，因该组年龄较大且吸烟比例高，无法直接归因于职业暴露。异氰酸酯组中5-mdC和5-hmdC的同步消耗及周转率下降提示表观遗传停滞，可能源于异氰酸酯对DNMT和TET酶的抑制。美容行业的高个体间变异性表明需关注个体增量变化而非群体均值。研究还讨论了5-hmdC和5-mdC的尿液特征在临床转化中的意义，但强调短期波动不能直接外推至慢性疾病风险。

结论翻译：本研究揭示了不同职业群体在单个标准工作周内尿液DNA修饰谱的快速变化。这些探索性发现提示了两种遗传毒性压力的分子特征，即氧化嘧啶不稳定性和表观遗传停滞。尽管泡沫切割组中尿液5-hmdU的明显积累是一个值得注意的发现，但由于该组中年龄和吸烟的混杂效应，其直接归因于颗粒物暴露仍不明确。相反，异氰酸酯组中5-mdC和5-hmdC的同步消耗可能提示基因组甲基化维持的破坏。通过2D UPLC-MS/MS精确定量的这些表观遗传和氧化标志物的趋势，表明它们作为非侵入性早期指标的潜在实用性。在更大、匹配更完善的队列中进行进一步验证以考虑人口统计学和生活方式因素后，这些分子特征可能最终为临床病理发生前的系统性生物学压力早期监测提供窗口。将对这些基因组“第六碱基”的评估纳入职业健康研究，可能最终有助于风险评估的精准化和长期疾病预防策略。