用于研究环境与生物介质中汞甲基化及甲基汞脱甲基化的四同位素加标(Quadruple Isotope Spiking)方法开发

《Analytica Chimica Acta》：DEVELOPMENT OF A QUADRUPLE ISOTOPE SPIKING METHODOLOGY TO STUDY MERCURY METHYLATION AND METHYLMERCURY DEMETHYLATION IN ENVIRONMENTAL AND BIOLOGICAL MEDIA.

　　

摘要：研究人员开发了一种四同位素加标(quadruple isotope spiking)方法并应用于培养实验中Hg物种特异性转化的研究。该方法首次实现了将培养过程中发生的物种间相互转化反应与样品前处理及分析过程中发生的转化反应区分开来。研究中建立了两种基于同位素模式解卷积(isotope pattern deconvolution, IPD)的四同位素加标数学模型，并在模型培养实验中追踪无机Hg(II)的甲基化及甲基汞(methylmercury, MeHg)的脱甲基化。结果表明，两种数学模型得到的Hg化合物浓度及转化反应速率一致，实现了计算程序的内部验证。该方法还被用于测定生物膜(biofilms)、沉积物、淡水及纯浮游植物培养体系中Hg(II)甲基化、MeHg脱甲基化为Hg(II)及MeHg损失(loss)速率，体系Hg浓度范围为pg L-1至μg L-1。在培养初始时间(t0)，应用四同位素加标方法后可完全校正并消除分析过程产生的非预期Hg(II)甲基化和/或MeHg脱甲基化干扰。分析过程中的甲基化可忽略不计(<0.2%)，而所有样品均观察到显著的分析过程MeHg脱甲基化(3%~10%)。天然Hg(II)甲基化仅见于生物膜和沉积物中；所有样品均可区分MeHg脱甲基化与MeHg损失。新生成化合物浓度的整体分析不确定度明显受浓度水平影响，其中实验色谱峰积分不确定度及富集同位素示踪剂同位素丰度不确定度是自然样品中新生成Hg化合物浓度总体分析不确定度的主要来源。本研究建立了一种研究生物及环境介质中Hg化合物转化的新方法，所基于的多重线性回归(mathematical framework based on multiple linear regression)原则上可拓展至其他具至少五种稳定同位素且有商品化物种标准物质(isotopically enriched species standards)的多同位素元素(如锡Sn)，但尚待实验验证。

《Analytica Chimica Acta》论文解读：环境与生物介质中汞甲基化与甲基汞脱甲基化的四同位素加标方法研究

一、研究背景与立项依据

汞(Hg)在环境中的反应性对其在地-气-水圈层的交换及在食物网中的生物累积至关重要，关键过程包括Hg(II)还原、元素汞氧化、Hg(II)甲基化为甲基汞(methylmercury, MeHg)及MeHg脱甲基化(demethylation)。同位素示踪剂是追踪这些耦合转化的有力工具。已有的双同位素加标(double-spiking)方法虽可同时测定Hg(II)甲基化和MeHg脱甲基化速率，但无法区分培养(incubation)过程中发生的物种转化与分析前处理(样品制备、衍生化、提取等)过程中发生的非预期转化(尤其是MeHg在分析过程中的脱甲基化)；此外，传统方法常将MeHg脱甲基化仅按MeHg损失(loss)计算，不能区分脱甲基化生成Hg(II)与还原脱甲基生成Hg(0)或其他损失途径(如生成二甲基汞并挥发)，且要求使用极高富集度(>99%)的昂贵示踪剂。虽有研究尝试用四个示踪剂校正，但缺乏能同时量化培养与分析过程转化因子、校正内源及外源Hg化合物浓度并提供完整数学框架的研究。因此，本研究旨在开发基于同位素模式解卷积(isotope pattern deconvolution, IPD)的四同位素加标(quadruple isotope spiking)方法，将培养过程转化(M₁、D₁)与分析过程转化(M₂、D₂)定量分离，准确测定各介质中Hg(II)甲基化潜力(M)、MeHg脱甲基化潜力(D)及MeHg损失，并给出严格的数学模型与不确定度评估。

二、主要关键技术方法概述

研究人员采用GC-ICP-MS(gas chromatography coupled to inductively coupled plasma mass spectrometry)测定Hg同位素分布。方法核心为向样品中分别加入两种培养示踪剂——富集同位素Hg(II)(^t1Hg(II)，如²⁰²Hg(II))和富集同位素MeHg(^t2MeHg，如²⁰¹MeHg)——用以追踪培养过程转化，以及在样品前处理前加入两种分析示踪剂——富集同位素Hg(II)(^t3Hg(II)，如¹⁹⁸Hg(II))和富集同位素MeHg(^t4MeHg，如¹⁹⁹MeHg)——用以校正分析过程物种间转化。基于测得的各Hg同位素峰面积比，通过多重线性回归求解自然本底Hg(II)、本底MeHg、四种示踪剂及新生成物种(^t1MeHg、^t2Hg(II))的摩尔分数，建立两种数学模型(Method 1假设本底与培养示踪剂反应活性相同；Method 2仅由培养示踪剂计算M₁、D₁，再对本底化合物作质量平衡校正)。先在无基质模型培养(微波诱导转化模拟M₁、D₁和M₂、D₂)中验证两模型一致性，再应用于野外采集的生物膜(biofilms)、沉积物(sediments)、淡水(freshwaters)及实验室纯培养浮游植物(Cyclotella meneghiniana)孵化实验，并用t₀对照样扣除初始非特异性转化，按公式算M(% h^-1)、D(% h^-1)及MeHg loss(% h^-1)，不确定度按Kragten法评估。

三、研究结果

INTRODUCTION（引言要点归纳）

已有单/双同位素稀释可校正分析过程物种转化并定量内源Hg(II)和MeHg，但无法区分培养与分析过程转化，也无法同时获得培养转化产率与新生成物种准确浓度。本研究填补此空白。

EXPERIMENTAL（实验设计要点）

模型培养实验在无基质条件下加入已知量自然丰度Hg(II)、MeHg及^t1Hg(II)、^t2MeHg，经微波加热诱发转化后加^t3Hg(II)、^t4MeHg再次微波以模拟分析过程转化。真实样品(生物膜、沉积物、淡水、浮游植物细胞培养)加入¹⁹⁹Hg(II)和²⁰¹MeHg为培养示踪剂，培养结束后提取/衍生化前加¹⁹⁸Hg(II)和²⁰²MeHg为分析示踪剂，GC-ICP-MS测定。设t₀对照及试剂空白。

Mathematical models for quadruple spiking experiments（四同位素加标数学模型）

基于各Hg同位素实测丰度是自然界本底、^t1Hg(II)、^t2MeHg、^t3Hg(II)、^t4MeHg五类来源同位素丰度的线性组合，用多重线性回归求各来源摩尔分数。Method 1假设本底与培养示踪剂具相同M₁、D₁，直接解算校正后浓度及M₁、D₁、M₂、D₂；Method 2先用培养示踪剂量平衡单独求M₁、D₁及分析过程转化，再对本底做质量平衡校正。两法在无差异反应性时数学等价；若不一致则仅Method 2可靠，两法比对可作诊断工具。

RESULTS AND DISCUSSION — Simulated incubations in the absence of matrix（无基质模型培养结果）

模型实验中二次微波产生显著分析过程脱甲基化(D₂=28±12% in 1:1混合)，未校正导致本底Hg(II)和MeHg回收率虚高(可达130%)，新生成^t2Hg(II)浓度严重高估；经四同位素加标校正后本底化合物回收率为98.1–103.6%(Hg(II))和97.7–105.0%(MeHg)，培养示踪剂回收率94.0–103.2%，新生成物种浓度显著不同于未校正值。M₁、D₁未校正值因D₂干扰出现偏差，而两数学模型(Method 1与Method 2)所得M₁、D₁及化合物浓度一致，证实数学模型内部验证通过。三种^t2MeHg:^t1Hg(II)比例(1:1、3:1、10:1)下方法均稳健。

RESULTS AND DISCUSSION — Incubations with real samples（真实样品培养结果）

所有基质分析过程甲基化M₂<0.2%，分析过程脱甲基化D₂=3–10%(生物膜/沉积物来自HNO₃微波提取+NaBEt₄衍生化；淡水/浮游植物仅来自NaBEt₄衍生化)，四同位素加标可完全校正。t₀时新生成化合物经校正后为零，排除初始非特异性转化干扰。天然Hg(II)甲基化仅见于生物膜和沉积物；所有样品可区分MeHg脱甲基化为Hg(II)(D)与MeHg loss(含还原脱甲基、二甲基汞挥发等)。经校正与未校正相比，浮游植物和沉积物中D(% h^-1)显著差异，生物膜和淡水因t₀与t_f具相似D₂自动抵消故无差异；MeHg loss无差异因其只涉及^t2MeHg消失量。比较D与MeHg loss可推断各基质中脱甲基生成Hg(II)所占比例(如光照生物膜中MeHg loss全来自D，暗条件<50%；淡水暗约30%、光约60%；浮游植物中占~60%)。

Estimation of method detection limits（检出限）

常规检出限：生物膜/沉积物中Hg(II)和MeHg均为30 pg kg^-1；淡水Hg(II) 15 pg L^-1、MeHg 3 pg L^-1；浮游植物细胞培养Hg(II) 50 pg L^-1、MeHg 30 pg L^-1。通过IPD求得的新生成物种及转化产率之实践检出限：M₁0.003%–0.6%，D₁0.3%–1%；生物膜/沉积物中新生成^t1MeHg ~0.03 μg L^-1级，^t2Hg(II) ~0.3 μg L^-1级，淡水中新生成^t2Hg(II)达pg L^-1级。

Evaluation of sources of uncertainty（不确定度评估）

新生成^t1MeHg浓度不确定度主要受MeHg中²⁰²Hg等同位素峰积分影响(贡献达94%)；新生成^t2Hg(II)浓度不确定度主要受培养示踪剂^t1Hg(II)的同位素丰度(²⁰⁰Hg、²⁰¹Hg、²⁰²Hg)及MeHg中²⁰¹Hg、²⁰²Hg峰积分影响。表明富集示踪剂同位素组成精确表征及GC峰面积分辨率与积分质量是方法关键控制点。

四、讨论与结论（翻译浓缩结论部分）

研究人员建立了两种基于多重线性回归的汞物种浓度及转化因子计算方法，互为内部验证。模型培养实验证明经典双同位素加标在未校正分析过程转化(M₂、D₂)时会导致M₁、D₁及新生成物种浓度偏差，而四同位素加标方法可有效校正之。应用于多种环境与生物介质(pg L^-1–μg L^-1Hg浓度范围)时，t₀时无非预期培养过程转化检出；分析过程Hg(II)甲基化可忽略(<0.2%)，而分析过程MeHg脱甲基化显著(3%–10%)并可被完全校正；该方法还可区分MeHg loss与MeHg脱甲基化为Hg(II)。综上，本四同位素加标方法为实验方案验证及复杂天然介质中Hg化合物转化潜力测定提供了可靠工具。所述基于同位素模式解卷积的双数学模型要求目标元素至少有五种可测同位素以容纳本底加四种富集示踪剂物种，在具备相应富集物种标准物质前提下，原则上可拓展至Sn(10种同位素)、Cd(8种)、Se(6种)等多同位素元素，但尚待实验证实。