响应面优化QuEChERS 结合UHPLC-MS/MS 法测定乌鳢肌肉中四环素类抗生素

李丽春，刘书贵，林嘉薇，尹 怡,*，郑光明，赵 城，马丽莎，单 奇，戴晓欣，魏琳婷，周 豪

（1.中国水产科学研究院珠江水产研究所，广东广州 510380；2.农业农村部外来入侵水生生物防控重点实验室，广东广州 510380；3.广东省水产动物免疫技术重点实验室，广东广州 510380）

四环素类抗生素是一类对常见的革兰氏阳性菌、阴性菌均具有较强抗菌作用的广谱性抗生素，被广泛用于动物性疾病的预防治疗[1]。近年来随着水产养殖业的快速发展，集约化高密度养殖越来越普遍，导致养殖环境水体恶化，各种疾病频发。养殖者常超剂量长时间使用抗菌药物进行病害防治。四环素类抗菌药作为水产养殖中的限用药物，在淡水养殖中广泛使用。水产养殖中不合理使用的四环素类抗生素，会通过水生动物的吸收和代谢转移至水体环境，给环境安全带来风险[2-3]。通过水生动物的生物富集，药物在其组织中残留，随着人们食用水产品而进入人体，给人体健康带来危害[4-5]。因此，美国、欧盟、日本等制定了四环素类抗生素在不同食品中的最大残留限量（MRL），四环素类药物在可食动物中的MRL 值分别为美国：2.0～12 mg/kg[6]，欧盟：100～600 μg/kg[7]，日本：0.2～1.2 mg/kg（总和）[8]，我国规定四环素类总残留量不超过100 μg/kg[9]。

目前动物性食品中四环素类抗菌药的残留检测方法主要有薄层色谱法[10]、化学发光法[11]、免疫分析法[12-13]、高效液相色谱法[8,14-15]和液相色谱串联质谱法[16-17]。其前处理大多采用固相萃取（SPE）[18-21]、固相分散萃取（SPME）[19,22]、超声辅助提取（UAE）[15]、加速溶剂萃取（ASE）等[17,23]较为繁琐复杂的方法。Quick、Easy、Cheap、Effective、Rugged、Safe（QuEChERS）前处理方法是近年来迅速发展的一种快速前处理方法，具有快速、简单、低成本、耐用、可靠、安全和多类型化合物同时检测的特点，在药物残留检测中广泛应用[18-19,24-25]，动物源食品中抗菌药物残留测定方面也有报道[23-24]，采用QuEChERS 前处理方法可以大大提高检测效率，获得更便捷和更准确的分析方法。目前四环素类药物残留检测方法大多采用固相萃取方法，操作繁琐且费时，本研究拟采用响应面法对影响方法准确性和精密度的前处理条件进行优化，有望建立一种简单快速且准确性和精密度较高的乌鳢肌肉中土霉素、金霉素、四环素和强力霉素的QuEChERS 结合UPLC-MS/MS 检测方法。并采用建立的方法测定广东主要养殖池塘的乌鳢样品中四环素残留量，分析其膳食风险，为乌鳢养殖中四环素类抗菌药物的风险评估及监管提供有效依据和技术支撑。

1 材料与方法

1.1 材料与仪器

空白乌鳢样品 来自于中国水产科学研究院珠江水产研究所养殖基地，养殖过程未使用过四环素类药物；实际检测的乌鳢样品 均采集于广东主要养殖区域的养殖池塘；盐酸土霉素（OTC，纯度96.5%）、盐酸金霉素（CTC，纯度93%）、盐酸多西环素（DC，纯度98.7%）、盐酸四环素（TC，纯度96.5%） 标准品，德国Dr.Ehrenstorfer 公司；氘代四环素（TC-D6，纯度>80%） 加拿大Toronto Research Chemicals（TRC）公司；二水合乙二胺四乙酸二钠（Na2EDTA·2H2O）、氯化钠（NaCl）、十二水合磷酸氢二钠（Na2HPO4·12H2O）、一水合柠檬酸（C6H8O7·H2O）、盐酸（HCl） 分析纯，广州化学试剂厂；甲酸、N-丙基乙二胺（PSA，40-63 μm）、中性氧化铝（Alumina-N，100-300 MESH）、弗罗里硅土（Florisil，150-250 μm）

德国CNW 公司；乙酸铵 色谱纯，阿拉丁试剂（上海）有限公司；十八烷基键合硅胶吸附剂（C18，40-60 μm）、石墨化碳黑（GCB，120-400 MESH） 天津Agela 公司；乙腈、甲醇 色谱纯，美国Burdick & Jackson 公司；实验用水为超纯水。

1290-6470 超高效液相色谱串联三重四极杆质谱仪，配有电喷雾离子源（ESI） 美国安捷伦公司；MS3 旋涡混匀 德国IKA 公司；Direct-Q3 超纯水机 美国Millipore 公司；TDL-5-A 离心机 上海安亭公司；N-EVAP-12 氮吹仪 美国Organomation 公司；3k15 高速离心机 美国Sigma 公司；KQ5200 型超声仪 昆山市超声仪器有限公司。

1.2 实验方法

1.2.1 标准溶液及试剂配制 标准溶液配制：准确称取适量的TC、OTC、CTC、DC 和内标TC-D65 种标准品，分别用甲醇溶解并定容，配制成质量浓度为100 μg/mL 的单标储备液。准确移取四种单标标准储备液各1 mL 用甲醇稀释并定容至10 mL，摇匀，配制成10 μg/mL 的混合标准中间溶液，将100 μg/mL的TC-D6内标溶液用甲醇稀释成10 μg/mL 标准中间液，于-18 ℃下避光保存备用。

0.1 mol/L Na2EDTA-McIlvaine（pH4.5）缓冲溶液配制[25]：分别准确称取37.2 g Na2EDTA·2H2O、10.9 g Na2HPO4·12H2O 和12.9 g C6H8O7·H2O，用1.00 L 超纯水充分溶解，加入适量盐酸调节pH 为4.5，于4 ℃冰箱保存备用，使用前放至室温。

1.2.2 前处理方法 提取：准确称取2.00 g 匀浆后的肌肉样品于50 mL 具塞离心管中，加入2 mL 0.1 mol/L Na2EDTA-McIlvaine（pH4.5）缓冲溶液，涡旋混匀后加入10 mL 1%甲酸酸化乙腈和2.0 g 氯化钠，涡旋振荡2 min，5000 r/min 离心5 min，上清液转移至50 mL 离心管中待净化。

净化：提取液中加入500 mg 无水硫酸钠、150 mg PSA 和200 mg C18，涡旋1 min，5000 r/min离心5 min，将上清液转移至50 mL 玻璃管中于45 °C 下氮气吹干，然后用1 mL 15%甲醇水充分溶解后，10000 r/min 离心5 min，过0.22 μm 微孔滤膜，待UHPLC-MS/MS 检测。

1.2.3 色谱条件 色谱柱：Agilent Zorbax SB-C18（2.1 mm×50 mm，18 µm）；流动相A 为0.5%甲酸水，B 为甲醇；流速：0.40 mL/ min；柱温：30 ℃；进样体积：5 µL，梯度洗脱程序列于表1。

表1 梯度洗脱程序Table 1 Program of gradient elution

1.2.4 质谱条件 离子源为电喷雾离子源正离子模式（ESI+）；毛细管电压3000 V；雾化器温度325 ℃；鞘气温度350 ℃；扫描模式为多反应监测（MRM），其他MRM 参数见表2。

表2 质谱参数Table 2 Mass spectrum parameter

1.2.5 单因素实验 为了建立高效、准确的乌鳢肌肉中四环素类化合物残留分析方法，针对影响四种四环素类化合物回收率的因素，通过空白基质加标法，优化提取剂类型，分别采用10 mL 甲醇、乙腈、酸化乙腈、酸化甲醇和酸化乙腈+2 mL Na2EDTA-McIlvanie 缓冲液对加标的样品进行提取，比较其回收率；确定提取剂类型后对提取剂体积进行优化选择，分别采用5、10、15 和20 mL 提取剂提取乌鳢肌肉中目标化合物，以回收率为100%确定提取剂体积；分别采用100 mg 的PSA、C18、GCB、中性氧化铝和弗罗里硅土对提取液进行净化，选择回收率较高且无杂质峰干扰的净化剂；对选择的净化剂用量进行优化选择，分别采用50、100、150、200、250 和300 mg 净化剂净化，比较其回收率和净化效果，每个处理设3 次平行。

1.2.6 响应面优化 通过单因素实验获得对乌鳢肌肉中四环素类化合物回收率影响较大的因素进行响应面优化，以期获得最优的前处理方法。分别以PSA 用量（X1）、C18用量（X2）和Na2EDTA-McIlvanie缓冲液体积（X3）为自变量，四种四环素的各自回收率为应变量，采用Box-Behnken Design（BBD）法设计3 因素3 水平优化试验[26]，试验设计因素与水平如表3 所示。

表3 Box-Behnken 试验设计的编码水平和相应的自变量实际水平Table 3 Coded levels and corresponding actual levels ofindependent variables used for Box-Behnken design

1.2.7 方法学验证 采用QuEChERS 结合UHPLCMS/MS 法测定水产品中4 种四环素类药物的残留量，用空白样品加标法评估方法检出限（LOD）和定量限（LOQ），在空白基质样品中添加4 种化合物，分别以信噪比（S/N）大于或等于3 和10 的最低浓度为方法LOD 和LOQ[27]。在空白基质中加入一定量的标准溶液后，采用建立的方法进行测定，通过加标样品的浓度与实际添加浓度的比值计算加标回收率，分别以LOQ、3 倍LOQ 和5 倍LOQ 为加标水平，每个水平重复6 次，通过计算加标回收率和相对标准偏差验证方法的准确性和精密度；通过分析3 d 内6 个添加样品在5.0 μg/kg 浓度下的回收率，计算相对标准偏差评价日间精密度。

采用提取后加标法，通过比较基质提取物和溶剂在0.5、1.0、2.0、5.0、10、20、50 和100 ng/mL 浓度下的校准曲线斜率评价基质效应（ME）。基质效应按照以下公式计算[28]：

式中，Sm为基质匹配标准溶液所得曲线的斜率，Ss为纯溶剂标准溶液所得曲线的斜率。当ME 的值在-20%到+20%之间时，认为其基质效应是可接受的；如果它们介于-50%～-20%或+20%～+50%之间，则被认为是中等基质效应；如果这些值低于-50%或高于+50%则认为具有较强的基质效应[28]。

1.2.8 实际样品检测及风险评估 在广东乌鳢主要养殖区域的养殖池塘采集乌鳢样品，测定其肌肉中4 种四环素类药物残留量。采用风险商（HQ，hazard quotient）对膳食中四环素类抗生素暴露量进行风险评估的表征。以人体重量为60 kg 的每日允许受摄入量（ADI）为指标，对乌鳢中四环素类化合物的残留风险进行评估。为了评估膳食风险，根据Diao 等[29]和Ramaswamy 等[30]的方法计算了膳食暴露量（EDI）和风险商（HQ）值。通过将人均日消费量与公式（2）中所示的药物残留量相乘，计算出膳食暴露量（EDI）。通过计算公式（3）中所示的风险商（HQ）进行膳食风险评估[31-32]。

式中，EDI 为膳食暴露量（μg/（kg·bw·d）），C 为样品中药物残留浓度（µg/kg），M 为人均消费量（g），BW 为成年平均体重60 kg；ADI 为每日允许摄入量（μg/（kg·bw·d））。当HQ 值大于1 时，风险被认为是不可接受的，而HQ 值小于1 则表示对人类的潜在风险较低[31]。

1.3 数据处理

实验数据采用Origin 软件2021 版本进行统计分析和图形输出，并使用Design-Expert 软件12.0 版本（State Ease, Inc.）进行响应面方差分析（ANOVA）和响应面图形输出，以确定线性项、二次项和交互项的回归系数，所有实验均进行3 次重复。

2 结果与分析

2.1 色谱条件的优化

四环素类化合物分子中的双键和含氧或氮原子的基团可与其它物质相互作用[33]，不仅可与金属离子形成二价、三价化合物，与蛋白质紧密结合，还可以与C18色谱柱上的硅醇基强烈作用，导致分离样品时出现拖尾峰[34]。因此需要选择合适的流动相，从而避免螯合物的形成及色谱柱的吸附，同时可以消除残留的金属离子[35-36]，获得较好的分离度和理想的色谱峰。ESI+离子模式下在流动相中加入一定量的甲酸，可以增加目标物的离子化状态，提高其离子化效率，达到改善峰形，提高分离度的效果。为此对常用的酸化试剂甲酸、乙酸铵和乙酸进行了优化比较，结果表明在甲醇-0.5%甲酸水中四种四环素类化合物可获得较高的响应且峰形尖锐对称。同时对质谱的电离源、锥孔电压和碰撞能等参数进行优化，最终优化值见表2。通过对色谱条件的优化获得的四环素类化合物的色谱图峰形较好（图1）。

图1 4 种四环素类药物在加标浓度为5 μg/kg 乌鳢肌肉样品中的MRM 谱图Fig.1 MRM chromatograms of the 4 tetracycline drugs in muscle samples of Channa argus at 5 μg/kg supplemented concentration

2.2 前处理条件优化

2.2.1 提取溶剂选择 为了建立简单、通用和高效的样品前处理方法，获得较为洁净的提取物，对分析物具有较低的离子抑制，本研究对提取溶剂进行了优化选择。QuEChERS 前处理方法常用的提取剂主要有甲醇、乙腈和酸化乙腈，因此采用加标回收法比较了甲醇、乙腈、酸化乙腈和酸化甲醇对水产品中四环素类化合物的提取效率（提取剂体积均为10 mL）。

不同提取剂对四种目标物的提取回收率见图2，结果表明，采用酸化乙腈和Na2EDTA-McIlvanie 缓冲液作为提取剂，四种四环素类化合物的提取效率较高，其加标回收率在80%～120%之间，可满足药物残留检测的要求。由于四环素类化合物可与基质中的二价金属离子（Ca2+、Mg2+、Zn2+和Cu2+）结合形成强配合物[37]，需要在提取时阻止二价金属离子和蛋白质在提取过程中与四环素类化合物作用[8,38]。为了使目标物从金属配合物中释放出来，须在基质中加入EDTA或柠檬酸等螯合剂。在样品提取时加入Na2EDTAMcIlvanie（pH4.5）缓冲液作为金属离子螯合剂，可提高基质中四环素类化合物的提取效率。因此选择Na2EDTA-McIlvanie 缓冲液和1%甲酸乙腈作为四环素类药物的提取剂。

图2 提取剂对四环素类化合物回收率的影响Fig.2 Effect of extraction agent on recovery of tetracyclines

同时比较了Na2EDTA-McIlvanie 缓冲液体积分别为1、2、3、4 和5 mL 时四种化合物的回收率，结果如图3A 所示，当缓冲液体积为2 mL 时可获得较高的回收率。此外比较了5、10、15 和20 mL 提取剂对四种目标物回收率的影响（图3B），结果表明，5 mL 1%甲酸乙腈难以充分提取回收率较低，当酸化乙腈体积大于10 mL 时回收率在80%～110%之间，因此选择1%甲酸乙腈的体积为10 mL。

图3 Na2EDTA-McIlvanie 和提取剂体积对四环素类化合物回收率的影响Fig.3 Effect of Na2EDTA-McIlvanie and extraction volume on recovery of tetracyclines

2.2.2 净化剂优化 QuEChERS 方法常用的净化剂有PSA、C18、GCB、中性氧化铝（Alumina-N）和弗罗里硅土（Florisil），PSA 通常用于去除提取液中的脂类和糖类物质，C18具有良好的除脂能力，GCB 可去除提取液中的色素成分，但对含苯环官能团的化合物有较强的吸附作用。本研究优化比较了PSA、C18、GCB、中性氧化铝（Alumina-N）和弗罗里硅土（Florisil）5 种吸附剂对4 种四环素类药物的净化效果（图4A）。在空白乌鳢肌肉样品中加入同一浓度水平（5.0 µg/kg）的四环素类混合标准溶液，分别用100 mg 的上述5 种吸附剂净化，比较其回收率。结果表明，GCB、中性氧化铝（Alumina-N）和弗罗里硅土对四种四环素类抗菌药具有极强的吸附性能，三种吸附剂中四环素类药物的回收率低于20%，而在吸附剂PSA 和C18净化下可获得较好回收率（80%～110%），因此选择PSA 和C18作为前处理方法的吸附剂进行优化。

图4 吸附剂种类和用量对加标回收率的影响Fig.4 The influence of type and dosage of cleaning agents on recoveries

分别对吸附剂PSA 和C18的用量进行优化比较，四环素类药物在不同剂量吸附剂下的回收率结果见图4B 和图4C。结果表明，吸附剂用量对回收率有较大影响。从图4B 可知，吸附剂C18的用量对四种四环素类药物的回收率有较大影响，C18用量从50 mg 增加至200 mg 时回收率逐渐增加，当C18用量为200 mg 时回收率较为理想（90.3%～97.8%），但C18用量大于250 mg 时，回收率明显降低。PSA 用量同样也影响回收率，当PSA 用量为150 mg 时，四种药物回收率在94.3%～96.8%之间，但是PSA 用量在100 mg 以下或200 mg 以上时回收率偏低。因此选择C18用量分别为150、200 和250 mg，PSA 用量为100、150 和200 mg 进行响应面优化。

2.3 响应面优化

在单因素实验结果基础上，选择PSA 用量（X1）、C18用量（X2）和Na2EDTA-McIlvanie 缓冲液体积（X3）进行三因素三水平的Box-Benken 响应面优化试验（BBD），响应值为四环素（Y1）、土霉素（Y2）、金霉素（Y3）和多西环素（Y4）的回收率，结果见表4。

表4 Box-Benken 响应面优化试验方案与结果Table 4 Experimental scheme and results of Box-Behnken design

2.3.1 模型建立及方差分析 对BBD 实验设计的结果进行了模型拟合和统计分析，以预测最佳QuEChERS 条件。采用Design-Expert 12.0 软件将实验数据拟合成二阶多项式模型，通过多元回归分析各自变量与响应之间的关系，获得4 种四环素类药物的拟合回归方程：

四环素（ TC） ： Y1=-146.72+0.55X1+1.82X2+31.54X3-0.0004X1X2+0.082X1X3+0.030X2X3-0.0024 X12-0.0046X22-12.42X32

土霉素（OTC）：Y2=-127.14+0.966X1+1.24X2+29.99X3-0.0019X1X2+0.0498X1X3+0.0121X2X3-0.0033X12-0.0031X22-10.08X32

金霉素（CTC）：Y3=-166.93+1.63X1+1.08X2+47.98X3+0.0003X1X2+0.037X1X3-0.0246X2X3-0.0059 X12-0.0030X22-12.53X32

多西环素（DC）：Y4=-158.37+1.30X1+1.37X2+32.36X3-0.0006X1X2+0.0585X1X3+0.0355X2X3-0.0044 X12-0.0035X22-12.30X32

为了检验模型的拟合程度，对模型进行方差分析和显著性检验，方差分析由FisherF检验进行推导[39]，结果见表5。一般而言F值越高，P值越低，回归模型越显著[40]，四种药物回归模型F值在14.59～46.43 之间（自由度df=9），P值均远低于0.01，回归模型差异极显著（表5）。四种药物回归方程的回归系数R2分别为0.9574（TC）、0.9835（OTC）、0.9494（CTC）和0.9575（DC），调整回归系数R2adj 分别为0.9027（TC）、0.9623（OTC）、0.8843（CTC）和0.9028（DC），表明模型拟合良好。且四种药物模型失拟值不显著，F值较低（0.0738～0.7510），P值均大于0.05（0.5762～0.9709），表明回归模型在本研究中可行[39]，具有较好的准确性和可靠性，可用来分析和预测各因素对4 种目标物回收率的影响。

表5 四环素类药物回收率响应面模型及预测结果方差分析Table 5 ANOVA of response surface model and predicted results for recovery of tetracyclines

2.3.2 响应面分析 为了最大限度的提高方法回收率，根据拟合方程进行响应面分析，绘制各变量之间的交互作用三维响应面曲线，四种四环素类药物的回收率受PSA 用量（X1）、C18用量（X2）和Na2EDTAMcIlvanie 缓冲液体积（X3）的影响如图5（TC:A1～A3，OTC:B1～B3，CTC:C1～C3，DC:D1～D3）所示。响应面的斜率越大，表明该因素对响应的影响更大[41]；响应面曲面可直观反映各因素间的交互作用对目标物回收率的影响，图中因素的曲面越陡表明该因素对目标物回收率影响越大。如由图5 可知，对于四环素和土霉素PSA 用量（X1）和Na2EDTA-McIlvanie缓冲液体积（X3）的响应面曲面比较陡，表明两者间交互作用显著；而PSA 用量（X2）和C18用量（X2）的响应面曲线比较缓，表明两者之间交互作用不明显。

图5 各因素交互作用对回收率影响的响应面图Fig.5 Response surface estimated showing the interaction of two different parameters on the recovery

通过Design-Expert 12.0 软件分析，以四环素类化合物的回收率均为100%，确定最佳QuEChERS条件为PSA：146.4 mg、C18：185.5 mg 和Na2EDTAMcIlvanie 缓冲液体积：1.95 mL，该最优条件下的四环素、土霉素、金霉素和多西环素的回收率预测值分别为：100%、102%、99.4%、99.4%。调整实际条件最佳QuEChERS 条件为PSA：150 mg、C18：200 mg和Na2EDTA-McIlvanie 缓冲液体积：2.0 mL，采用优化后的参数进行验证实验，经6 次平行实验后测得平均回收率分别为：101±1.7%（TC）、97.7±3.8%（OTC）、101±3.1%（CTC）和95.7±2.8%（DC），四种化合物实际回收率与模型预测值差异不明显，说明拟合模型优化出的QuEChERS 条件较为准确。

2.4 方法学验证

2.4.1 方法的线性、检出限、定量限 在上述优化的仪器条件下，分析一系列不同浓度的基质样品，以各化合物的质量浓度为横坐标，定量离子峰面积为纵坐标，绘制标准曲线。结果表明，四种四环素类药物在线性范围内均呈良好的线性关系（r>0.999）。添加不同浓度混合标准溶液至空白样品进行测试，分别以目标峰S/N=3 和10 确定检出限（LOD）和定量限（LOQ），结果如表6 所示。在0.5～100 ng/mL 浓度范围内，4 种四环素类化合物的线性关系良好，相关系数均大于0.999；检出限为0.25～0.75 μg/kg，定量限为0.84～2.00 μg/kg，满足国内外药物残留检测要求。

表6 四环素类化合物线性方程、相关系数、检出限、定量限和基质效应Table 6 Linear equation, correlation coefficient, LOD, LOQ and matrix effects for tetracyclines

2.4.2 回收率和精密度 采用基质加标法对方法回收率和精密度进行了考察，分设置LOQ、3 倍LOQ和5 倍LOQ 三种添加浓度水平，按照本方法确定的前处理方法和仪器条件进行测定，每个浓度水平设置6 次重复，结果见表7。四种四环素类化合物日内和日间回收率在89.5%～118%之间，精密度小于10%，与现有文献报道的方法相比（表8），本方法回收率和精密度较好，可满足国内外药物残留检测要求。

表7 四环素类抗生素在乌鳢肌肉中的加标回收率和精密度（n=6）Table 7 Recovery and precision of tetracyclines in the muscle of Channa argus (n=6)

表8 本方法与其他报道的水产品中四环素类药物残留检测方法比较Table 8 Comparison of the proposed method with other reported methods for the determination of TCs in fish samples

2.4.3 基质效应 采用空白基质提取后加标法，用空白乌鳢肌肉样品按照优化的前处理方法进行前处理后，分别采用定溶液（15%甲醇水）和空白基质溶液配制0.5～100 ng/mL 的标准溶液，按照优化好的仪器条件进行测定，计算乌鳢基质中四环素类化合物的基质效应，结果见表6。四种四环素类化合物的基质效应在-36.0%～5.30%之间，其中四环素存在基质增强效应，其他三种物质存在基质抑制效应，为减少基质效应对检测结果的影响，在实际样品检测中应采用基质标准曲线进行定量校准。

2.5 实际样品测定及风险分析

从广东省乌鳢主要养殖区域采集乌鳢肌肉样品33 个，按照本研究优化的方法检测样品中四环素类化合物的残留量，并通过计算风险商（hazard quoti-ent，HQ）对4 种四环素类化合物的残留风险进行风险评估。乌鳢样品中四环素类化合物残留结果见表9。

表9 乌鳢肌肉中四环素类抗生素残留浓度和风险商Table 9 The residual concentrations and HQ of tetracyclines in the muscle of Channa argus

以我国制定的抗生素每日允许摄入量（Acceptable Daily Intake，ADI μg/（kg·bw·d））为标准进行评价，四环素、金霉素和土霉素ADI 均为30 μg/（kg·bw·d），多西环素ADI 为3.0 μg/（kg·bw·d）。根据国家统计局数据，2020 年我国人均水产品消费量为38.08 g/d，通过计算可知，土霉素HQ 值在1.07×10-4～1.78×10-4之间，四环素HQ 值在7.51×10-5～1.36×10-4，金霉素HQ 值在5.33×10-5～7.31×10-4，多西环素HQ 值为2.28×10-4～3.07×10-3。所有四环素类化合物HQ 值均远低于1.0，表明通过食用广东养殖乌鳢摄入四环素类药物对人体的健康没有影响，风险较低，膳食安全性较高。

3 结论

本研究采用响应面法优化了四环素类抗生素在乌鳢肌肉中的QuEChERS 结合超高效液相色谱-串联质谱（UHPLC-MS/MS）检测方法，样品用10 mL 1%甲酸酸化乙腈和2 mL 0.1 mol/L Na2EDTA-McIlvaine（pH4.5）缓冲溶液提取，经150 mg PSA 和200 mg C18净化后，采用UHPLC-MS/MS 测定，内标法定量。结果表明，该方法简便快捷、溶剂用量少，且准确性和精密度较高。在0.5～100 ng/mL 浓度范围内线性良好，相关系数均大于0.999；检出限为0.25～0.75 μg/kg，定量限为0.84～2.00 μg/kg。在三个添加浓度水平下四种四环素类化合物的日内和日间回收率在89.5%～118%之间，精密度小于10%，满足国内外药物残留检测要求，可用于乌鳢肌肉中四环素类药物残留检测。本方法的定量限较已有报道的方法低，回收率、精密度和基质效应与现有方法基本一致。采用优化的快速检测方法检测了广东省乌鳢主要养殖区域的33 个乌鳢样品，通过计算风险商（HQ）进行风险评估，结果表明广东乌鳢中四环素的膳食风险较低，其残留对人体健康无影响。