黃曲霉毒素B1誘導(dǎo)大鼠肝細(xì)胞凋亡率劑量效應(yīng)模型的構(gòu)建

黃奧迪，嚴(yán)佳惠，張昭寰,2，劉海泉,2,3，趙 勇,2,3，歐 杰,2,3,*

（1.上海海洋大學(xué)食品學(xué)院，上海 201306；2.上海水產(chǎn)品加工及貯藏工程技術(shù)研究中心，上海 201306；3.農(nóng)業(yè)農(nóng)村部水產(chǎn)品貯藏保鮮質(zhì)量安全風(fēng)險(xiǎn)評(píng)估實(shí)驗(yàn)室（上海），上海 201306）

黃曲霉毒素是一類(lèi)由黃曲霉和寄生曲霉等菌株產(chǎn)生的高毒性次級(jí)代謝物的總稱(chēng)[1]，廣泛存在于谷物、蔬菜、肉類(lèi)等食物中，嚴(yán)重危害食品安全[2]。其中，黃曲霉毒素B1（aflatoxin B1，AFB1）是黃曲霉毒素中出現(xiàn)頻率最高且毒性和致癌性最強(qiáng)的[3]，AFB1主要毒性作用的靶器官是肝臟[4]。Mughal等[5]研究已經(jīng)證實(shí)AFB1通過(guò)引起肝細(xì)胞凋亡和干擾細(xì)胞酶活性引起肝功能障礙，作用機(jī)制涉及腫瘤抑制基因形成DNA加合物，從而導(dǎo)致基因突變，對(duì)人類(lèi)造成致癌性，影響人體健康。

劑量效應(yīng)模型（dose-response model，DRM）是食品風(fēng)險(xiǎn)評(píng)估框架中危害特征描述的重要組成部分，其本質(zhì)在于定量描述暴露于危害物的劑量與發(fā)生不良反應(yīng)概率之間的關(guān)系[6]。目前國(guó)內(nèi)外已開(kāi)展的DRM研究包括沙門(mén)氏菌[7]、副溶血性弧菌[8]、蠟樣芽孢桿菌[9]等，但是國(guó)內(nèi)關(guān)于真菌毒素構(gòu)建DRM的研究較少，且聚焦于AFB1構(gòu)建DRM也缺乏系統(tǒng)性的研究[10]。AFB1生物性危害主要通過(guò)產(chǎn)生相應(yīng)的毒素或者宿主進(jìn)食具有感染性的病原體從而影響人體健康，在我國(guó)肝癌的高發(fā)地區(qū)廣西由于人們經(jīng)常食用含有AFB1的食物，導(dǎo)致肝臟負(fù)擔(dān)增加，從而引發(fā)肝臟疾病。產(chǎn)生黃曲霉毒素的真菌大多存在于氣候炎熱潮濕的地區(qū)，由于收獲前和收獲后的真菌污染，在食物中發(fā)現(xiàn)了黃曲霉毒素，黃曲霉毒素對(duì)食物的污染速度和程度取決于溫度、濕度、土壤和儲(chǔ)存條件[11]。AFB1慢性毒性不僅會(huì)使體質(zhì)量減輕、肝臟相關(guān)指標(biāo)上升，還會(huì)增加腎臟的壓力，造成尿酸濃度升高；而AFB1急性毒性會(huì)誘發(fā)肝臟損傷、腎衰竭、肝細(xì)胞壞死，甚至引起死亡[12]。歐洲食品安全局（European Food Safety Authority，EFSA）前期已經(jīng)按照國(guó)家和年齡分組開(kāi)展了黃曲霉毒素的暴露評(píng)估，并進(jìn)行了定性的風(fēng)險(xiǎn)評(píng)估分析[13]，在此基礎(chǔ)上進(jìn)行AFB1的DRM研究可為其風(fēng)險(xiǎn)評(píng)估提供定量依據(jù)。

接觸0.36 mg/kgmb及以上的AFB1會(huì)導(dǎo)致動(dòng)物和人類(lèi)急性肝損傷和死亡。已有研究報(bào)道有關(guān)AFB1導(dǎo)致肝細(xì)胞凋亡的定性研究，但關(guān)于AFB1急性暴露于肝臟的定量數(shù)據(jù)和有關(guān)劑量效應(yīng)模型的研究很少。因此，本研究以大鼠為研究對(duì)象，通過(guò)設(shè)計(jì)實(shí)驗(yàn)研究不同劑量AFB1對(duì)大鼠肝細(xì)胞凋亡率的影響，并基于實(shí)驗(yàn)數(shù)據(jù)應(yīng)用指數(shù)模型、對(duì)數(shù)模型、Weibull模型、Hill模型等常用模型擬合可能的劑量效應(yīng)關(guān)系，通過(guò)數(shù)學(xué)檢驗(yàn)進(jìn)行參數(shù)比較，選擇最優(yōu)的DRM為黃曲霉毒素風(fēng)險(xiǎn)評(píng)估提供理論參考。

1 材料與方法

1.1 動(dòng)物、材料與試劑

48 只SPF級(jí)雄性Wister大鼠，體質(zhì)量180～200 g，購(gòu)自上海斯萊克實(shí)驗(yàn)動(dòng)物有限責(zé)任公司。飼養(yǎng)于環(huán)境控制的動(dòng)物實(shí)驗(yàn)室內(nèi)，自由飲水飲食，飼養(yǎng)溫度22～24 ℃，相對(duì)濕度50%～55%，光照條件：亮暗周期各12 h。本研究進(jìn)行的動(dòng)物實(shí)驗(yàn)程序均按照上海中醫(yī)藥大學(xué)《實(shí)驗(yàn)動(dòng)物護(hù)理與使用指南》進(jìn)行，并經(jīng)中國(guó)動(dòng)物中心動(dòng)物倫理委員會(huì)和上海中醫(yī)藥大學(xué)批準(zhǔn)，審查批準(zhǔn)編號(hào)：PZSHUTCM210115014。

AFB1上海吉至科技有限公司；二甲基亞砜（dimethyl sulfoxide，DMSO）、磷酸鹽緩沖液（phosphate buffered saline，PBS）生工生物工程（上海）股份有限公司；水合氯醛 阿達(dá)瑪斯試劑有限公司；多聚甲醛、無(wú)水乙醇 國(guó)藥集團(tuán)化學(xué)試劑有限公司；TUNEL染色試劑盒 武漢塞維爾生物科技有限公司。

1.2 儀器與設(shè)備

DW-88L338超低溫保存箱 青島海爾特種電器有限公司；恒溫培養(yǎng)箱 賽默飛世爾科技（中國(guó)）有限公司；HH-S11-2-S恒溫水浴鍋 上海新苗醫(yī)療器械制造有限公司；ZH5102A電子天平 珠恒電子有限公司；JB-L5凍臺(tái) 武漢俊杰電子有限公司；RM2016病理切片機(jī) 上海徠卡儀器有限公司；KD-P組織攤片機(jī)浙江省金華市科迪儀器設(shè)備有限公司；GFL-230烤箱天津市萊玻瑞儀器設(shè)備有限公司；G6004防脫載玻片、TSY-B脫色搖床、WG1066-1組化筆 武漢塞維爾生物科技有限公司；Giotto染色機(jī) 意大利DIAPATH公司；Eclipse E100正置光學(xué)顯微鏡 日本Nikon公司。

1.3 方法

1.3.1 造模及分組

48 只大鼠經(jīng)過(guò)1 周的適應(yīng)性飼養(yǎng)后，隨機(jī)分為4 組（每組12 只）。空白組（CK組）：腹腔注射與AFB1劑量組相同量的4% DMSO溶液；AFB1低劑量組（LD組）：腹腔注射劑量為1 mg/kgmb的AFB1工作液一次；中劑量組（MD組）：腹腔注射劑量為2 mg/kgmb的AFB1工作液一次；高劑量組（HD組）：腹腔注射劑量為3 mg/kgmb的AFB1工作液一次，AFB1溶于4% DMSO溶液。AFB1的劑量根據(jù)參考文獻(xiàn)[14]和預(yù)實(shí)驗(yàn)結(jié)果初步確定。

1.3.2 實(shí)驗(yàn)動(dòng)物處理

造模后的48 h對(duì)大鼠使用10%水合氯醛（350 mg/kgmb）進(jìn)行麻醉處理，待麻醉之后對(duì)大鼠進(jìn)行解剖，分離大鼠的肝臟組織并使用4%多聚甲醛溶液浸泡固定。

1.3.3 肝細(xì)胞凋亡率的測(cè)定（TUNEL染色）

在4%多聚甲醛中固定24 h后，取出肝臟組織進(jìn)行石蠟包埋和切片（4 μm），然后對(duì)石蠟切片脫蠟復(fù)水，將蛋白酶K工作液滴加到組織切片處理，PBS漂洗，之后進(jìn)行破膜處理和室溫平衡，加入TUNEL反應(yīng)混合液（TdT酶∶dUTP∶Buffer＝1∶5∶50（體積比））在37 ℃恒溫箱中孵育2 h，PBS洗滌后滴加3%牛血清白蛋白覆蓋組織，加入雙抗，之后滴加DAPI染液，室溫避光放置10 min，滴加抗熒光猝滅封片液后使用熒光顯微鏡觀察拍照并進(jìn)行細(xì)胞數(shù)量計(jì)算。肝細(xì)胞凋亡率按式（1）計(jì)算：

1.3.4 劑量效應(yīng)模型構(gòu)建

常見(jiàn)的DRM包括指數(shù)模型、對(duì)數(shù)模型、Weibull模型、Hill模型等。以下模型中d為大鼠的注射劑量（mg/kgmb），P為大鼠肝細(xì)胞凋亡率，A、B、C均表示常數(shù)（表1），應(yīng)用Origin Pro for Windows 8.0軟件中的Levenberg-Marquardt法（非線性最小二乘法）進(jìn)行擬合。

表1 常見(jiàn)劑量效應(yīng)模型及其相關(guān)參數(shù)Table 1 Common dose-response models and related parameters

1.3.5 數(shù)學(xué)檢驗(yàn)與評(píng)價(jià)

參照Kasuya等[15]建模數(shù)學(xué)檢驗(yàn)參數(shù)決定系數(shù)R2和Lin等[16]赤池信息量準(zhǔn)則（Akaike information criterion，AIC）、貝葉斯信息準(zhǔn)則（Bayesian information criterion，BIC）進(jìn)行模型擬合度評(píng)價(jià)和最佳模型選擇，相關(guān)表達(dá)式如下：

式中：SSE為方差分析的殘差平方和；SST為總離差平方和。

式中：k為模型中未知參數(shù)個(gè)數(shù)；n為樣本數(shù)量；L為模型中極大似然函數(shù)值。

其中AIC代表實(shí)際值與預(yù)測(cè)值之間的符合程度，理論上AIC越小越好，然而當(dāng)模型的復(fù)雜度提高（參數(shù)個(gè)數(shù)k增大）時(shí)，似然函數(shù)值也會(huì)增大，從而導(dǎo)致AIC變小。當(dāng)k過(guò)大時(shí)，似然函數(shù)值增速減緩，模型過(guò)于復(fù)雜造成過(guò)度擬合現(xiàn)象。針對(duì)該問(wèn)題，BIC引入了與參數(shù)個(gè)數(shù)相關(guān)，與AIC相比更大的懲罰項(xiàng)，考慮了當(dāng)參數(shù)個(gè)數(shù)變大時(shí)，模型精度提高所造成的模型復(fù)雜度過(guò)高的現(xiàn)象。

2 結(jié)果與分析

2.1 AFB1對(duì)大鼠肝細(xì)胞凋亡的影響

圖1為T(mén)UNEL染色結(jié)果部分切片，藍(lán)色部分代表正常肝細(xì)胞，其核膜結(jié)構(gòu)完整；綠色部分代表凋亡肝細(xì)胞，其核膜結(jié)構(gòu)呈現(xiàn)簇狀或碎片狀。由圖1A～D可知，與CK組相比，LD組、MD組和HD組的凋亡率顯著上升，且隨著劑量的提高細(xì)胞凋亡率也有所上升（P＜0.05）。

圖1 AFB1對(duì)大鼠肝細(xì)胞凋亡的影響Fig.1 Effect of AFB1 on apoptosis in rat hepatocytes

2.2 大鼠肝細(xì)胞凋亡率劑量效應(yīng)關(guān)系的擬合

根據(jù)大鼠腹腔注射AFB1誘導(dǎo)肝細(xì)胞凋亡的定量數(shù)據(jù)（圖1），應(yīng)用不同DRM擬合非線性關(guān)系如圖2所示。由圖2可知，除了在迭代次數(shù)范圍內(nèi)未收斂的Weibull模型1以外，其余DRM都達(dá)到了收斂的統(tǒng)計(jì)要求，并且相同的模型出現(xiàn)了幾乎重合的現(xiàn)象。其中DRM所預(yù)測(cè)的半數(shù)有效量（median effective dose，ED50），即能引發(fā)50%細(xì)胞凋亡的劑量均為1.91 mg/kgmb左右，與實(shí)際實(shí)驗(yàn)所用劑量大致相當(dāng)，也與EFSA[17]所預(yù)測(cè)的隨著AFB1濃度提高，肝細(xì)胞癌（hepatocellular carcinoma，HCC）發(fā)生概率顯著升高的趨勢(shì)相同。對(duì)于AFB1造成肝細(xì)胞凋亡率的ED50可通過(guò)動(dòng)物模型毒理學(xué)實(shí)驗(yàn)進(jìn)行進(jìn)一步確認(rèn)，本實(shí)驗(yàn)結(jié)果僅作為預(yù)測(cè)劑量效應(yīng)關(guān)系使用。

圖2 AFB1誘導(dǎo)大鼠肝細(xì)胞凋亡率不同劑量效應(yīng)模型的擬合曲線Fig.2 Fitting curves of different dose-response models for AFB1-induced apoptosis rate in rat hepatocytes

從擬合結(jié)果來(lái)看，Weibull模型、Gamma模型和劑量效應(yīng)曲線分別對(duì)應(yīng)圖2中的曲線G、H、K、L，其決定系數(shù)R2＞0.9960（表2），在擬合效果上優(yōu)于其他模型。相同的模型在外延至高劑量和最小有效量時(shí)幾乎重合，預(yù)測(cè)細(xì)胞凋亡率結(jié)果與Ruggeberg等[18]實(shí)驗(yàn)測(cè)得大鼠靜脈注射AFB1肝細(xì)胞凋亡率結(jié)果相似，證明本實(shí)驗(yàn)DRM的可信性，但食品等殘留AFB1導(dǎo)致肝細(xì)胞凋亡的相關(guān)數(shù)據(jù)需要深入實(shí)驗(yàn)。

表2 AFB1誘導(dǎo)大鼠肝細(xì)胞凋亡率劑量效應(yīng)模型的參數(shù)估計(jì)及數(shù)學(xué)檢驗(yàn)Table 2 Parameter estimation and mathematical testing of doseresponse models for AFB1-induced apoptosis rate in rat hepatocytes

2.3 劑量效應(yīng)模型的數(shù)學(xué)檢驗(yàn)

對(duì)構(gòu)建的DRM進(jìn)行數(shù)學(xué)檢驗(yàn)（表2），對(duì)比參數(shù)決定系數(shù)R2和AIC可知，劑量效應(yīng)曲線（R2＝0.9970）和指數(shù)模型2、3（AIC＝33.91）優(yōu)于其他模型。

將BIC納入考慮的范疇內(nèi)，Gamma模型的BIC值（BIC＝17.19）僅比劑量效應(yīng)曲線的BIC值（BIC＝15.88）大，且AIC值優(yōu)于大多數(shù)模型，參數(shù)決定系數(shù)也處于上述考慮范圍內(nèi)。Gamma模型（R2＝0.9964，AIC＝37.40，BIC＝17.19）優(yōu)于其他模型綜上所述，Gamma模型可作為AFB1誘導(dǎo)大鼠肝細(xì)胞凋亡率最優(yōu)的DRM用于劑量效應(yīng)模型擬合的理論參考。

3 討論

急性暴露于AFB1不僅會(huì)影響肝臟和腎臟等器官[19]，還會(huì)抑制人類(lèi)的免疫系統(tǒng)，使其受到傳染病的危害[20]。據(jù)報(bào)告，嚙齒動(dòng)物口服AFB1的半數(shù)致死量在1～18 mg/kgmb之間[21]，而其他物種的研究表明AFB1的劑量小于1 mg/kgmb也可能導(dǎo)致死亡[22]，證明AFB1不僅會(huì)對(duì)器官造成不可逆的影響，還會(huì)威脅生命健康。EFSA食品鏈污染物小組根據(jù)各項(xiàng)動(dòng)物研究數(shù)據(jù)推斷認(rèn)為AFB1每日攝入量不得超過(guò)0.4 μg/kgmb，由于使用人類(lèi)志愿者代替動(dòng)物模型研究獲得劑量效應(yīng)數(shù)據(jù)的方法存在許多困難，關(guān)于AFB1的人類(lèi)數(shù)據(jù)遠(yuǎn)不如動(dòng)物數(shù)據(jù)那樣豐富，因此一般采用換算系數(shù)的方法外推有關(guān)人類(lèi)的數(shù)據(jù)[23]。

本研究選用大鼠Wister腹腔注射肝細(xì)胞凋亡率的實(shí)驗(yàn)數(shù)據(jù)建立DRM，由于AFB1檢測(cè)報(bào)告大多屬于定性檢測(cè)，關(guān)于大鼠的定量數(shù)據(jù)較少，缺少一定的參考，因此可能存在些許誤差。Williams[24]對(duì)虹鱒魚(yú)進(jìn)行的劑量效應(yīng)研究為AFB1風(fēng)險(xiǎn)評(píng)估模型提供了支持，未來(lái)定量風(fēng)險(xiǎn)評(píng)估必然是風(fēng)險(xiǎn)評(píng)估的新趨勢(shì)所在[25]。

指數(shù)模型是食品中有害物建立DRM的首要選擇，如沙門(mén)氏菌[26]、諾如病毒[27]等微生物均對(duì)指數(shù)模型表現(xiàn)出了良好的擬合性。然而有關(guān)霉菌毒素的DRM很少有研究，Gilbert-Sandoval等[28]基于生理動(dòng)力學(xué)建模的反向劑量法預(yù)測(cè)AFB1在大鼠和人中的急性肝毒性劑量反應(yīng)數(shù)據(jù)，通過(guò)定量體外劑量評(píng)估AFB1急性肝毒性并預(yù)測(cè)體內(nèi)毒性評(píng)估，為本實(shí)驗(yàn)預(yù)測(cè)AFB1對(duì)于肝細(xì)胞凋亡率趨勢(shì)變化提供了參考。通過(guò)擬合模型推導(dǎo)低劑量效應(yīng)存在不準(zhǔn)確的可能性，會(huì)出現(xiàn)低劑量外推數(shù)據(jù)溢出問(wèn)題[29]，因此在將本研究結(jié)果外推至人體或?qū)嶋H應(yīng)用時(shí)應(yīng)謹(jǐn)慎。

建立劑量效應(yīng)模型后有必要采用志愿者實(shí)驗(yàn)進(jìn)行驗(yàn)證，但是文獻(xiàn)中關(guān)于人類(lèi)吸收AFB1的數(shù)據(jù)僅有Jubert等[30]召集3 位男性志愿者口服低劑量（30 ng）數(shù)據(jù)，在暴露1 h后觀測(cè)到血漿中的最大放射性尚未達(dá)到肝細(xì)胞凋亡的最低有效劑量，對(duì)于本研究不具備參考性。同時(shí)AFB1的感染劑量也會(huì)因?yàn)槭荏w、環(huán)境因素、毒素作用等因素導(dǎo)致數(shù)據(jù)不同[31]，因此本研究采用常用的數(shù)學(xué)檢驗(yàn)系數(shù)R2和AIC以及BIC評(píng)價(jià)和選擇。通過(guò)比較，建議選用Gamma模型作為黃曲霉毒素導(dǎo)致肝細(xì)胞凋亡率最優(yōu)劑量效應(yīng)模型。

隨著風(fēng)險(xiǎn)評(píng)估模型的發(fā)展，在建立劑量效應(yīng)模型的時(shí)候可以引入其他信息，例如引入新的參數(shù)TPI（time post inoculation），用于預(yù)測(cè)時(shí)間和劑量同時(shí)作為自變量引起的反應(yīng)變化，從而建立時(shí)間劑量效應(yīng)模型評(píng)估反映隨時(shí)間和劑量的發(fā)展趨勢(shì)，同時(shí)也能結(jié)合宿主因素(如年齡、性別）或病原體因素(如微生物氣溶膠直徑、感染性）等[32]。本研究結(jié)論僅針對(duì)于模型方法論時(shí)證明其可行性較好，在外推應(yīng)用時(shí)應(yīng)謹(jǐn)慎使用[33]。

4 結(jié)論

本研究基于AFB1誘導(dǎo)大鼠肝細(xì)胞凋亡的定量分析結(jié)果，分別應(yīng)用指數(shù)模型、對(duì)數(shù)模型、Weibull模型、Hill模型、Gamma模型、劑量效應(yīng)曲線6 種模型擬合了AFB1誘導(dǎo)大鼠肝細(xì)胞凋亡率的劑量效應(yīng)模型，其中指數(shù)模型、Gamma模型、劑量效應(yīng)曲線表現(xiàn)出良好的擬合性，根據(jù)數(shù)學(xué)檢驗(yàn)結(jié)果，建議選用Gamma模型作為AFB1誘導(dǎo)大鼠肝細(xì)胞凋亡率最優(yōu)的劑量效應(yīng)模型用于AFB1定量風(fēng)險(xiǎn)評(píng)估框架中的危害特征描述。