基于Monte Carlo模擬法對(duì)液態(tài)乳中黃曲霉毒素M1的風(fēng)險(xiǎn)評(píng)估

耿夢(mèng)夢(mèng)，徐明芳,*，王 陽，黎 明，陳耕南

（1.暨南大學(xué)生命科學(xué)技術(shù)學(xué)院，廣東 廣州 510632；2.暨南大學(xué)應(yīng)急管理研究中心，廣東 廣州 510632）

中國是液態(tài)乳的巨大消費(fèi)國，液態(tài)乳和乳制品對(duì)各個(gè)年齡段的人適用，是因?yàn)橐簯B(tài)乳是人體內(nèi)的主要營養(yǎng)來源之一。根據(jù)聯(lián)合國糧食及農(nóng)業(yè)組織預(yù)測(cè)，2015年世界鮮奶產(chǎn)量達(dá)到8.01億 t，增長1.5%，中國原料乳產(chǎn)量大致為3 700萬 t左右，全國液態(tài)乳產(chǎn)量為2 074.63萬 t，同比增長4.32%[1]。生鮮液態(tài)乳中有害物質(zhì)的來源主要有霉菌毒素污染、農(nóng)藥殘留、獸藥殘留、重金屬污染、微生物污染和人為添加劑。而黃曲霉毒素作為霉菌毒素之一就是其中最為重要的一種[2]。

黃曲霉毒素是由一類黃曲霉（Aspergillus flavus）和寄生曲霉（Aspergillus parasiticus）等真菌產(chǎn)生的次級(jí)代謝產(chǎn)物。迄今為止，至少有300多種不同的黃曲霉毒素被發(fā)現(xiàn)，包括黃曲霉毒素B1（aflatoxin B1，AFB1）、B2（aflatoxin B2，AFB2）、G1（aflatoxin G1，AFG1）、G2（aflatoxin G2，AFG2）和M1（aflatoxin M1，AFM1）等。其中AFM1是由AFB1在肝臟中經(jīng)微粒體細(xì)胞色素P450羥基化代謝的產(chǎn)物，當(dāng)奶牛食用了被AFB1污染的飼料后，其生產(chǎn)的液態(tài)乳及乳制品中就會(huì)污染。AFM1的基本結(jié)構(gòu)是由二氫呋喃環(huán)和氧雜萘鄰?fù)獦?gòu)成，進(jìn)入機(jī)體后能誘發(fā)p53基因特定密碼子的特定堿基突變，從而誘發(fā)癌變，是人類原發(fā)性肝癌的主要致病因素之一。除此之外，AFM1還具有致畸、致突變和免疫抑制等非致癌軀體有害作用。世界衛(wèi)生組織（World Health Organization，WHO）的國際癌癥研究機(jī)構(gòu)根據(jù)其細(xì)胞毒性及致癌效應(yīng)，將其定義為Ⅱ類致癌物[4]。對(duì)此，國際組織及相關(guān)國家制定了食品中AFM1的限量標(biāo)準(zhǔn)。歐盟規(guī)定，液態(tài)乳和乳制品中AFM1的含量不能超過0.05 μg/kg。美國與我國對(duì)乳及乳制品中AFM1的限量標(biāo)準(zhǔn)為0.5 μg/kg，嬰兒乳品中不得檢出[5]。

暴露評(píng)估是風(fēng)險(xiǎn)評(píng)估的核心步驟，分為點(diǎn)評(píng)估和概率評(píng)估兩種。點(diǎn)評(píng)估模型采用食物高消費(fèi)量和污染物高殘留量進(jìn)行計(jì)算，一般按照y=計(jì)算，其中y為某種危害物的人群暴露量（攝入量）；Xi,0.95為第i類食物消費(fèi)量分布的95%分位數(shù)；Ci,max為第k類食物中某一危害物的最大殘留量；p為消費(fèi)食物種類數(shù)目；為被評(píng)估人群的平均體質(zhì)量。

點(diǎn)評(píng)估法一般選取第95百分位數(shù)或第97.5百分位數(shù)用以反映高暴露量人群的情況。Michlig等[6]通過點(diǎn)評(píng)估法評(píng)估了在大型貯罐液態(tài)乳中AFM1的污染水平，并進(jìn)一步確定了阿根廷原料乳暴露評(píng)估與AFM1污染有關(guān)的風(fēng)險(xiǎn)因子，認(rèn)為商業(yè)飼料污染是最主要的風(fēng)險(xiǎn)因子。目前，對(duì)黃曲霉毒素的暴露評(píng)估主要采用的是點(diǎn)評(píng)估法，該法簡(jiǎn)單易行，評(píng)估成本低，但忽略了個(gè)體差異，結(jié)構(gòu)較為保守，無法量化個(gè)體水平消費(fèi)量和食品中化學(xué)物水平的變異，也無法對(duì)參數(shù)估計(jì)的不確定性做出說明[7]。近年來，由美國Palisade公司開發(fā)的基于Monte Carlo模擬技術(shù)的風(fēng)險(xiǎn)分析軟件@Risk已被應(yīng)用于安全定量風(fēng)險(xiǎn)評(píng)估工作中[8]。Monte Carlo模擬是應(yīng)用于人體暴露評(píng)估的最為廣泛的概率評(píng)估法，此方法一般先擬合分布，然后進(jìn)行抽樣，可以認(rèn)為是從總體中抽樣，能夠有效地量化變異性。美國環(huán)境保護(hù)署在風(fēng)險(xiǎn)分析政策中將其定為風(fēng)險(xiǎn)分析基本方法[9]，概率暴露評(píng)估是使用概率模型表現(xiàn)群體中不同風(fēng)險(xiǎn)等級(jí)可能性或表征暴露評(píng)估中的不確定性，應(yīng)用隨機(jī)（概率）模型建立和分析風(fēng)險(xiǎn)的不同情形[10-12]。由于概率評(píng)估的結(jié)果能反映風(fēng)險(xiǎn)分布更為真實(shí)的情況，因此Monte Carlo方法越來越多地應(yīng)用到食品安全風(fēng)險(xiǎn)分析中[13-14]。

健康風(fēng)險(xiǎn)評(píng)估過程中將有毒物質(zhì)分為基因毒物質(zhì)和軀體毒物質(zhì)，其中基因毒性物質(zhì)包括放射性污染物和致癌性化學(xué)物，軀體毒物質(zhì)為非致癌性有毒化學(xué)物質(zhì)，根據(jù)有毒物質(zhì)的類別，評(píng)價(jià)模型分為化學(xué)致癌物風(fēng)險(xiǎn)模型、非化學(xué)致癌物風(fēng)險(xiǎn)模型和放射性物質(zhì)風(fēng)險(xiǎn)模型[15-17]。AFM1作為一種遺傳毒性致癌物，在風(fēng)險(xiǎn)評(píng)估中被視為無閾值化合物，即不存在不產(chǎn)生致癌效應(yīng)的最小安全劑量（no observed adverse effect level，NOAEL）[18]。2005年，食品添加劑聯(lián)合專家委員會(huì)（Joint FAO/WHO Expert Committee on Food Additives，JECFA）在第64次會(huì)議上首次提出了將暴露邊界比（margin of exposure，MOE）應(yīng)用于遺傳毒性致癌物風(fēng)險(xiǎn)描述上[19]，即用實(shí)驗(yàn)動(dòng)物致癌劑量（或人群致癌劑量）和人群攝入量的比值來描述風(fēng)險(xiǎn)。AFM1除能促進(jìn)腫瘤的發(fā)生外，還具有非致癌毒性，目前國際組織對(duì)黃曲霉毒素進(jìn)行的風(fēng)險(xiǎn)評(píng)估還沒有建立其健康指導(dǎo)值，但有文獻(xiàn)報(bào)道了同樣具有致肝癌毒性的微囊藻毒素的每日耐受攝入量（tolerable daily intake，TDI）[20-24]，評(píng)估微囊藻毒素的非致癌健康風(fēng)險(xiǎn)，為AFM1的非致癌健康風(fēng)險(xiǎn)評(píng)估提供了參考。

本研究基于液態(tài)乳中AFM1污染調(diào)查得到的數(shù)據(jù)和中國居民營養(yǎng)與健康調(diào)查數(shù)據(jù)庫中食物消費(fèi)量和人口學(xué)數(shù)據(jù)，構(gòu)建暴露評(píng)估模型，對(duì)消費(fèi)液態(tài)乳對(duì)AFM1的暴露水平進(jìn)行概率評(píng)估，并結(jié)合歐盟食品安全局（European Food Safety Authority，EFSA）基于人群流行病學(xué)數(shù)據(jù)建立的劑量-反應(yīng)關(guān)系所提出的BMDL10（對(duì)照組中出現(xiàn)10%肝癌發(fā)生率的95%統(tǒng)計(jì)學(xué)可信區(qū)間下限值）[25]對(duì)其產(chǎn)生的致癌風(fēng)險(xiǎn)進(jìn)行描述，并參考微囊藻毒素的TDI值，設(shè)定AFM1的TDI值范圍，評(píng)估非致癌健康風(fēng)險(xiǎn)，以期為相關(guān)部門制定風(fēng)險(xiǎn)管理政策提供科學(xué)依據(jù)和數(shù)據(jù)支撐。

1 材料與方法

1.1 材料與試劑

液態(tài)乳樣品的采集：綜合考慮居民日常飲乳習(xí)慣及市場(chǎng)消費(fèi)情況等因素，本研究以不同品種奶源液態(tài)乳為研究對(duì)象，于2015年12月至2016年5月隨機(jī)采樣，總共采379 個(gè)樣品，其中，樣品A 91 個(gè)、樣品B 72 個(gè)、樣品C 72 個(gè)、樣品D 72 個(gè)、樣品E 72 個(gè)（A、B、C、D、E分別代表液態(tài)乳的不同品種，相對(duì)應(yīng)的數(shù)值表示該品種液態(tài)乳的采樣次數(shù)）。

AFM1標(biāo)準(zhǔn)品（純度＞99.8%） 以色列Fermentek公司；乙腈、甲醇（色譜純） 天津市科密歐化學(xué)試劑有限公司。

1.2 儀器與設(shè)備

1100高效液相色譜（high performance liquid chromatography，HPLC）儀（包括數(shù)據(jù)分析平臺(tái)、熒光檢測(cè)器） 美國Agilent公司；AFM1免疫親和柱（柱容量100 ng） 北京華安麥科有限公司；超純水儀 美國Thermo公司。

1.3 方法

1.3.1 樣品制備

根據(jù)GB 5413.37—2010《食品安全國家標(biāo)準(zhǔn) 乳和乳制品中黃曲霉毒素M1的測(cè)定》測(cè)定樣品乳中的AFM1含量，精確量取50 mL液態(tài)乳樣品，置于離心管中，于35 ℃恒溫水浴鍋中溫浴15 min。10 ℃、6 000 r/min離心10 min，恢復(fù)至室溫。量取40 mL上清液分次轉(zhuǎn)移到氣控操作架上20 mL注射器中，調(diào)節(jié)壓力，控制試樣以2～3 mL/min穩(wěn)定的流速過柱，直至2～3 mL空氣通過柱體。向上述注射器內(nèi)加入10 mL水，以穩(wěn)定的流速洗柱，棄去全部流出液，并使2～3 mL空氣通過柱體。準(zhǔn)確加入2.0 mL乙腈至上述注射器，以1～2 mL/min穩(wěn)定流速將親和柱上的AFM1洗脫下來，收集全部洗脫液于錐形管中。然后用氮?dú)庠?0 ℃下將洗脫液蒸發(fā)至體積50 μL，再用水定容到1 mL，過0.22 μm微孔過濾器后轉(zhuǎn)移至樣品瓶，供HPLC分析。

1.3.2 HPLC檢測(cè)條件

色譜柱：C18柱（250 mm×4.6 mm，5 μm）；流動(dòng)相：25%乙腈水溶液；流速1.0 mL/min；熒光檢測(cè)器：發(fā)射波長360 nm、激發(fā)波長440 nm；柱溫30 ℃。

1.3.3 風(fēng)險(xiǎn)評(píng)估方法

1.3.3.1 暴露評(píng)估模型的構(gòu)建

本研究基于Monte Carlo模擬和Bootstrap抽樣技術(shù)，采用概率評(píng)估方法構(gòu)建暴露評(píng)估模型。參照美國環(huán)境保護(hù)署（United States Environmental Protection Agency，US EPA）發(fā)布的化學(xué)污染物健康風(fēng)險(xiǎn)評(píng)估模型，采用日暴露量（chronic daily intake oral，CDIoral）模擬評(píng)估居民飲用液態(tài)乳暴露AFM1的水平[26]。CDIoral根據(jù)公式（1）計(jì)算。

式中：CDIoral為居民消費(fèi)液態(tài)乳對(duì)AFM1的日吸收量/（μg/（kg?d））；CW為液態(tài)乳中AFM1的含量/（μg/kg）；IR為日均攝入量；EF為暴露頻率（365 d/年）；ED為暴露持續(xù)時(shí)間（70 年）；AT為平均時(shí)間（70×365 d）；BW為人群平均體質(zhì)量/kg[27]。

1.3.3.2 風(fēng)險(xiǎn)描述模型的構(gòu)建

基因毒性致癌物的風(fēng)險(xiǎn)描述模型：AFM1能誘發(fā)p53基因特定密碼子的堿基突變，從而誘發(fā)癌變，為基因毒性致癌物[28]?；蚨拘灾掳┪餅闊o閾值化合物，即不存在NOAEL。因此，本研究將MOE應(yīng)用于基因毒性致癌物的風(fēng)險(xiǎn)描述，根據(jù)CDIoral模擬評(píng)估，應(yīng)用MOE對(duì)液態(tài)乳中AFM1進(jìn)行風(fēng)險(xiǎn)描述[18]，MOE根據(jù)公式（2）計(jì)算。

式中：P O D為致癌效應(yīng)分離點(diǎn)（p o i n t o f departure），基于人群流行病學(xué)數(shù)據(jù)確定的數(shù)據(jù)中，參考EFSA建立的劑量-反應(yīng)關(guān)系，BMDL10（對(duì)照組中出現(xiàn)10%肝癌發(fā)生率的95%統(tǒng)計(jì)學(xué)可信區(qū)間下限值）為0.87（μg/（kg?d））[25]。MOE≥100，表示液態(tài)乳中AFM1的致癌風(fēng)險(xiǎn)是可以接受的；MOE＜100，表示液態(tài)乳中AFM1的致癌風(fēng)險(xiǎn)超過了可以接受的限度，應(yīng)當(dāng)采取適當(dāng)?shù)娘L(fēng)險(xiǎn)管理措施。

軀體毒性致癌物的風(fēng)險(xiǎn)描述模型：AFM1除了具有促進(jìn)腫瘤發(fā)生的毒性外，還具有皮膚致敏、免疫抑制等有害作用[29]，本研究以危害指數(shù)（hazard index，HI）來表征液態(tài)乳中AFM1對(duì)居民的非致癌風(fēng)險(xiǎn)。世界衛(wèi)生組織（World Health Organization，WHO）對(duì)AFM1的ADI值并未作出規(guī)定，根據(jù)USEPA規(guī)定，在未規(guī)定ADI值情況下，其值可用TDI代替。本研究以同為致肝癌毒素的微囊藻的TDI（0.04 μg/（kg?d））為參考[20]，設(shè)定AFM1的TDI值分別為0.005、0.010、0.020、0.040、0.060 μg/（kg?d），根據(jù)公式（3）分別計(jì)算居民的HI[18]。

HI＞1，表示液態(tài)乳中AFM1的非致癌風(fēng)險(xiǎn)是可以接受的；HI＜1，表示液態(tài)乳中AFM1的非致癌風(fēng)險(xiǎn)超過了可以接受的限度，應(yīng)當(dāng)采取適當(dāng)?shù)娘L(fēng)險(xiǎn)管理措施。

1.4 數(shù)據(jù)統(tǒng)計(jì)分析

1.4.1 風(fēng)險(xiǎn)評(píng)估軟件

以上模擬抽樣過程在美國Palisade公司的基于Monte Carlo模擬技術(shù)的@Risk7.0專業(yè)風(fēng)險(xiǎn)分析軟件環(huán)境下運(yùn)行，且評(píng)估中使用的各種參數(shù)對(duì)應(yīng)的概率分布采用@Risk7.0提供的標(biāo)準(zhǔn)分布函數(shù)來顯示。

1.4.2 風(fēng)險(xiǎn)評(píng)估相關(guān)數(shù)據(jù)處理

低于檢測(cè)限的AFM1污染數(shù)據(jù)按照WHO和US EPA建議的數(shù)據(jù)處理方法，即1/2檢出限替換處理[7]。采用SPSS 13.0統(tǒng)計(jì)軟件的單因素方差分析比較差異顯著性。每個(gè)液態(tài)乳樣品平行測(cè)定3 次，用Excel 2013計(jì)算結(jié)果。

變異系數(shù)又稱“標(biāo)準(zhǔn)差率”，是衡量不同液態(tài)乳中各觀測(cè)值變異程度的一個(gè)統(tǒng)計(jì)量，采用標(biāo)準(zhǔn)差與平均值的比值（相對(duì)值）來表示。

2 結(jié)果與分析

2.1 HPLC檢測(cè)液態(tài)乳中AFM1

圖1 AFM1標(biāo)準(zhǔn)品（a）與液態(tài)乳樣品（b）的HPLC圖Fig.1 HPLC chromatogram of AFM1 standard (a) and milk samples (b)

AFM1標(biāo)準(zhǔn)品及不同乳源液態(tài)乳樣品通過HPLC檢測(cè)的色譜圖如圖1，AFM1標(biāo)準(zhǔn)品的出峰時(shí)間為5.7 min，以此做定性分析。

2.2 液態(tài)乳中AFM1污染的時(shí)間分布

圖2 液態(tài)乳中AFM1污染的時(shí)間分布Fig.2 Temporal distribution of AFM1 concentrations in contaminated milk

由圖2可知，從2015年12月開始，液態(tài)乳樣品中AFM1含量逐漸上升，到次年3月份升到最高，隨后略有下降。這是因?yàn)?月份為南方的回南天，氣候溫暖、濕度大，適合黃曲霉菌的大量繁殖，到4～5月份后，溫度逐漸升高，天氣炎熱干燥，AFM1的污染程度略微下降。對(duì)于不同奶源液態(tài)乳來說，樣品D的污染水平較高，A和C相當(dāng)，B的污染水平最低。

表1 液態(tài)乳中AFM1的污染水平比較Table1 Comparison of AFM1 levels in milk samples

由表1可知，樣品C中AFM1污染水平范圍明顯寬于其他品牌的液態(tài)乳，最大污染水平達(dá)到0.62 μg/kg；樣品B、D、A污染水平范圍次之；E液態(tài)乳污染水平范圍最窄（污染水平較集中），僅為0.07～0.31 μg/kg。不同奶源液態(tài)乳中AFM1的污染平均值經(jīng)單因素方差分析表明，樣品D與B的污染水平差異顯著。

2.3 液態(tài)乳中AFM1的暴露評(píng)估

2.3.1 液態(tài)乳中AFM1污染水平的分布擬合

利用基于Monte Carlo模擬技術(shù)的@Risk7.0軟件對(duì)不同液態(tài)乳中AFM1污染水平監(jiān)測(cè)數(shù)據(jù)進(jìn)行分布擬合。結(jié)果顯示，AFM1污染數(shù)據(jù)比較符合Pareto、Expon、Laplace等分布類型。運(yùn)用@Risk7.0軟件自帶的赤池信息量準(zhǔn)則（Akaike’s information criterion，AIC）、貝葉斯信息準(zhǔn)則（Bayesian information criterions，BIC）、Chi-Squared、Anderson-Darling和Kolmogorov-Smirnov統(tǒng)計(jì)檢驗(yàn)方法對(duì)上述分布的擬合度進(jìn)行檢驗(yàn)。

表2 5 個(gè)樣品液態(tài)乳中AFM1最佳擬合分布及相關(guān)參數(shù)Table2 Optimum fitting distribution of AFM1 concentrations in five milk samples and related parameters

表3 兩種理論分布對(duì)樣品分布擬合的相關(guān)參數(shù)Table3 Parameters of two theoretical distribution models for milk sample

圖3 不同液態(tài)乳中AFM1含量數(shù)據(jù)分布擬合的概率密度曲線Fig.3 Probability density curve for distribution fitting of AFM1 concentration in different milk samples

綜合考慮擬合結(jié)果（表2），確定樣品中AFM1污染數(shù)據(jù)最佳概率擬合分布類型有Pareto和Expon分布?；跇悠稡、D的污染水平達(dá)顯著差異，以此為例，擬合后各參數(shù)值比較結(jié)果見表3，Pareto比Expon更適合液態(tài)乳樣品的分布。因此，5 個(gè)樣品分布分別記為Risk Pareto（2.661 5，0.029）、Risk Pareto（3.765 2，0.029）、Risk Pareto（3.370 8，0.029）、Risk Pareto（2.854 6，0.029）和Risk Pareto（2.423 5，0.029）（圖3）。

2.3.2 居民消費(fèi)不同液態(tài)乳AFM1的日均暴露量計(jì)算

采用1.3.3.1節(jié)中確定的液態(tài)乳中AFM1污染水平的暴露模型，結(jié)合相關(guān)的暴露參數(shù)，利用@Risk7.0風(fēng)險(xiǎn)分析軟件，分別隨機(jī)從不同樣品液態(tài)乳中AFM1的污染水平分布中抽取數(shù)值計(jì)算AFM1膳食暴露量概率分布，每次模擬過程迭代10 000 次，暴露結(jié)果如圖4所示。

圖4 飲用液態(tài)乳途徑的AFM1日均暴露量概率分布Fig.4 Probabilistic distribution of average daily intake of AFM1 through milk consumption

百分位數(shù)是個(gè)類似平均值的統(tǒng)計(jì)學(xué)概念，是把所有數(shù)值從小到大排列后，按照百分?jǐn)?shù)進(jìn)行劃分（表4），百分位數(shù)10%表示整個(gè)暴露量區(qū)間從低到高至數(shù)據(jù)10%時(shí)的暴露量，依次類推，100%即最高暴露量的位點(diǎn)。暴露評(píng)估的結(jié)果顯示，消費(fèi)不同液態(tài)乳途徑的AFM1日暴露量的平均值依次為0.000 23、0.000 20、0.000 21、0.000 22、0.000 25 μg/（kg·d），中位數(shù)依次為0.000 19、0.000 17、0.000 18、0.000 18、0.000 19 μg/（kg·d），第99百分位數(shù)的AFM1暴露量分別為0.000 82、0.000 49、0.000 57、0.000 73、0.000 97 μg/（kg·d），其他百分位數(shù)概率下的暴露量如表4所示。

表4 居民日飲用液態(tài)乳途徑AFM1暴露量百分位數(shù)值Table4 Average daily intake of AFM1 from milk products

2.4 樣品E中AFM1的風(fēng)險(xiǎn)描述

2.4.1 樣品E中AFM1的致癌性風(fēng)險(xiǎn)描述

由2.3.2節(jié)中AFM1分別在不同樣品液態(tài)乳中的日均暴露量來看，飲用樣品E液態(tài)乳的日均暴露量在各個(gè)百分位數(shù)均高于其他液態(tài)乳樣品。因此本研究以樣品E為例進(jìn)行風(fēng)險(xiǎn)描述。采用1.3.3.2節(jié)中致癌性風(fēng)險(xiǎn)的計(jì)算公式，結(jié)合暴露參數(shù)，利用@Risk軟件，分別隨機(jī)從樣品E液態(tài)乳中AFM1污染水平分布中抽取數(shù)值計(jì)算不同人群消費(fèi)樣品E液態(tài)乳暴露AFM1的致癌性風(fēng)險(xiǎn)概率，每次模擬循環(huán)10 000次，風(fēng)險(xiǎn)評(píng)估結(jié)果如表5所示。

表5 不同人群通過樣品E液態(tài)乳暴露AFM1致癌性風(fēng)險(xiǎn)Table5 Carcinogenic risk associated with dietary exposure to AFM1 from milk sample E for different populations

根據(jù)表5結(jié)果，不同人群消費(fèi)樣品E液態(tài)乳暴露AFM1的致癌性風(fēng)險(xiǎn)值MOE均大于100，說明在抽樣檢測(cè)的時(shí)間段內(nèi)，不同人群通過飲用液態(tài)乳暴露AFM1對(duì)健康的致癌性風(fēng)險(xiǎn)很小。對(duì)于樣品E液態(tài)乳，不同人群間致癌性風(fēng)險(xiǎn)評(píng)估經(jīng)t檢驗(yàn)表明，結(jié)果存在極顯著差異（t＝0.000 1＜0.01）。不同年齡組人群中，幼兒（2～3、4～6 歲）和兒童（7～13 歲）3 個(gè)年齡組人群由于體質(zhì)量較輕，而液態(tài)乳攝入量相對(duì)較高，致癌性風(fēng)險(xiǎn)明顯高于其他年齡組人群，需重點(diǎn)關(guān)注；青少年（14～17 歲）和60歲以上的老年人的致癌性風(fēng)險(xiǎn)相當(dāng)；18～59 歲的成年人的致癌性風(fēng)險(xiǎn)最小。對(duì)于不同地域的人群，城市居民通過液態(tài)乳攝入AFM1的致癌性風(fēng)險(xiǎn)值均顯著高于農(nóng)村人群，是由于城市居民的飲奶量遠(yuǎn)高于農(nóng)村人群。而對(duì)于不同性別的人群，由于男性在體質(zhì)量和液態(tài)乳攝入量方面均高于女性，經(jīng)配對(duì)樣本t檢驗(yàn)得兩個(gè)人群的風(fēng)險(xiǎn)不存在顯著差異（t＝0.276＞0.05）。

2.4.2 樣品E中AFM1的非致癌性風(fēng)險(xiǎn)描述

采用1.3.3.2節(jié)中AFM1的非致癌性風(fēng)險(xiǎn)模型，計(jì)算不同人群消費(fèi)樣品E液態(tài)乳暴露AFM1的致癌性風(fēng)險(xiǎn)概率，每次模擬循環(huán)10 000 次，風(fēng)險(xiǎn)評(píng)估結(jié)果如表6所示。

表6 不同人群通過樣品E液態(tài)乳對(duì)AFM1非致癌性風(fēng)險(xiǎn)Table6 Non-carcinogenic health risk associated with dietary exposure to AFM1 from milk sample E for different populations

由表6可知，不同人群飲用樣品E液態(tài)乳攝入AFM1的非致癌性風(fēng)險(xiǎn)HI值均大于1，（與致癌性風(fēng)險(xiǎn)是相同的），說明在抽樣檢測(cè)的時(shí)間段內(nèi)，通過飲用液態(tài)乳攝入AFM1對(duì)健康的非致癌性風(fēng)險(xiǎn)很小。在設(shè)定的5 個(gè)不同水平的TDI值下，非致癌風(fēng)險(xiǎn)最高的農(nóng)村女幼兒（2～3 歲）HI值為15.60，其他不同人群的HI值均高于顯著高于15.60，說明非致癌風(fēng)險(xiǎn)很小。但在特定區(qū)域特定狀況下對(duì)于幼兒和兒童的風(fēng)險(xiǎn)，仍需引起相關(guān)部門的高度關(guān)注。

3 結(jié) 論

本次調(diào)查的379 個(gè)液態(tài)乳樣品中，AFM1的檢出率為19.79%，平均含量為0.19 μg/kg，最高污染水平為0.62 μg/kg。在檢測(cè)的不同液態(tài)乳樣品中，樣品E的檢出率最高，為23.61%，平均含量為0.18 μg/kg，而樣品D的平均污染水平最高，為0.22 μg/kg。樣品B的變異系數(shù)最大，為98.78%，分布較廣泛，樣品E的分布最集中，變異系數(shù)為38.80%。

在構(gòu)建液態(tài)乳中AFM1暴露評(píng)估模型的基礎(chǔ)上，結(jié)合相關(guān)暴露參數(shù)，采用基于Monte Carlo模擬技術(shù)的@Risk7.0風(fēng)險(xiǎn)評(píng)估軟件，分別對(duì)液態(tài)乳中AFM1的污染水平進(jìn)行分布擬合，并通過5 種擬合度檢驗(yàn)，選出擬合最優(yōu)分布，依次為Risk Pareto（2.661 5，0.029）、Risk Pareto（3.765 2，0.029）、Risk Pareto（3.370 8，0.029）、Risk Pareto（2.854 6，0.029）和Risk Pareto（2.423 5，0.029）。結(jié)合相關(guān)暴露參數(shù)，計(jì)算了通過飲用液態(tài)乳攝入AFM1的日均暴露量。

運(yùn)用@Risk軟件，分別以MOE和HI值表征飲用液態(tài)乳途徑的AFM1的致癌性和非致癌性風(fēng)險(xiǎn)。結(jié)果表明，不同人群通過飲用不同奶源的液態(tài)乳攝入AFM1的致癌性風(fēng)險(xiǎn)值均顯著大于100，非致癌風(fēng)險(xiǎn)值在TDI取值0.005 μg/（kg?d）時(shí)仍均大于1，在抽樣檢測(cè)的時(shí)間段內(nèi)，僅通過飲用液態(tài)乳攝入AFM1對(duì)健康的風(fēng)險(xiǎn)很小。

[1] 李勝利, 姚琨, 曹志軍, 等. 國家奶業(yè)體系: 2015年奶牛產(chǎn)業(yè)技術(shù)發(fā)展報(bào)告[J]. 中國畜牧雜志, 2016, 52(2): 34-39.

[2] CAMPAGNOLLO F B, GANEV K C, KHANEGHAH A M, et al. The occurrence and effect of unit operations for dairy products processing on the fate of aflatoxin M1: a review[J]. Food Control, 2016, 68:310-329. DOI:10.1016/j.foodcont.2016.04.007.

[3] SCAGLIONI P T, BADIALE-FURLONG E. Rice husk as an adsorbent: a new analytical approach to determine af l atoxins in milk[J].Talanta, 2016, 152: 423-431. DOI:10.1016/j.talanta.2016.02.042.

[4] FALLAH A A, FAZLOLLAHI R, EMAMI A. Seasonal study of aflatoxin M1contamination in milk of four dairy species in Yazd, Iran[J]. Food Control, 2016, 68: 77-82. DOI:10.1016/j.foodcont.2016.03.018.

[5] 衛(wèi)生部. 食品中真菌毒素限量: GB 2761-2011[S]. 北京: 中國標(biāo)準(zhǔn)出版社, 2011: 2-4.

[6] MICHLIG N, SIGNORINI M, GAGGIOTTI M, et al. Risk factors associated with the presence of aflatoxin M1in raw bulk milk from Argentina[J]. Food Control, 2016, 64: 151-156. DOI:10.1016/j.foodcont.2015.12.025.

[7] 葉孟亮, 聶繼云, 徐國鋒, 等. 蘋果中4 種常用農(nóng)藥污染水平及其膳食暴露評(píng)估[J]. 中國農(nóng)業(yè)科學(xué), 2016, 49(7): 1289-1302.DOI:10.3864/j.issn.0578-1752.2016.07.007.

[8] 白新明. 蔬菜農(nóng)藥殘留對(duì)人體健康急性風(fēng)險(xiǎn)概率評(píng)估研究[J]. 食品科學(xué), 2014, 35(5): 208-212. DOI:10.7506/spkx1002-6630-201405041.

[9] 劉瀟威. 農(nóng)產(chǎn)品中重金屬風(fēng)險(xiǎn)評(píng)估的研究與進(jìn)展[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào), 2007, 26(1): 15-18.

[10] 田文, 裴斐, 文良, 等. 食用油中有機(jī)污染物的篩查與風(fēng)險(xiǎn)評(píng)估[J].食品科學(xué), 2016, 37(24): 203-208. DOI:10.7506/spkx1002-6630-201624032.

[11] 董慶利. 蒸煮米飯中蠟樣芽孢桿菌的定量暴露評(píng)估[J]. 食品科學(xué),2013, 34(21): 306-310. DOI:10.7506/spkx1002-6630-201321061.

[12] 楊岳平, 胡長鷹, 李克亞, 等. 毒理學(xué)關(guān)注閾值方法在食品接觸材料風(fēng)險(xiǎn)評(píng)估中的應(yīng)用[J]. 食品科學(xué), 2015, 36(23): 334-339.DOI:10.7506/spkx1002-6630-201523060.

[13] 錢永忠, 李耘, 周德慶, 等. 農(nóng)產(chǎn)品質(zhì)量安全風(fēng)險(xiǎn)評(píng)估-原理、方法和應(yīng)用[M]. 北京: 中國標(biāo)準(zhǔn)出版社, 2007: 113-139.

[14] 王永杰, 賈東紅. 健康風(fēng)險(xiǎn)評(píng)價(jià)中的不確定性分析[J]. 環(huán)境工程,2003, 21(6): 66-69. DOI:10.13205/j.hjgc.2003.06.022.

[15] US Environmental Protection Agency. Guideline for Exposure Assessment[R]. Washington D.C.: US EPA, 1992.

[16] 陳勝軍, 李來好, 楊賢慶, 等. 我國水產(chǎn)品安全風(fēng)險(xiǎn)來源與風(fēng)險(xiǎn)評(píng)估研究進(jìn)展[J]. 食品科學(xué), 2015, 36(17): 300-304. DOI:10.7506/spkx1002-6630-201517055.

[17] 吳黎明, 李耘, 余林生. 蜂產(chǎn)品中有害物暴露評(píng)估方法[J]. 食品科學(xué),2011, 32(11): 324-328.

[18] 羅祎. 食品安全風(fēng)險(xiǎn)分析化學(xué)危害評(píng)估[M]. 北京: 中國質(zhì)檢出版社,2012: 83-90.

[19] World Health Organization. Evaluation of certain food contaminants:sixty-fourth report of the Joint FAO/WHO Expert Committee on Food Additives[J]. WHO Technical Report Series, 2006, 930: 8-25.

[20] 王偉琴. 飲用水源水中微囊藻毒素的遺傳毒性與健康風(fēng)險(xiǎn)評(píng)價(jià)[D].杭州: 浙江大學(xué), 2010: 6-17.

[21] 王超, 彭濤, 呂怡兵, 等. 江南某城市飲用水及其水源水中微囊藻毒素調(diào)查及初步健康風(fēng)險(xiǎn)評(píng)價(jià)[J]. 環(huán)境化學(xué), 2014, 33(7): 1237-1238.DOI:10.7524/j.issn.0254-6108.2014.07.027.

[22] 詹曉靜, 向壘, 李彥文, 等. 農(nóng)田土壤中微囊藻毒素污染特征及風(fēng)險(xiǎn)評(píng)價(jià)[J]. 中國環(huán)境科學(xué), 2015, 35(7): 2129-2136.

[23] 楊曉紅, 蒲朝文, 張仁平, 等. 水體微囊藻毒素污染對(duì)人群的非致癌健康風(fēng)險(xiǎn)[J]. 中國環(huán)境科學(xué), 2013, 33(1): 181-185.

[24] 孟元華, 朱鵬飛, 龔燕, 等. 太湖水產(chǎn)品中微囊藻毒素-RR的污染狀況及初步健康風(fēng)險(xiǎn)評(píng)估[J]. 食品安全質(zhì)量檢測(cè)學(xué)報(bào), 2016, 7(5):1794-1797.

[25] 白藝珍, 丁小霞, 李培武, 等. 應(yīng)用暴露限值法評(píng)估中國花生黃曲霉毒素風(fēng)險(xiǎn)[J]. 中國油料作物學(xué)報(bào), 2013, 35(2): 211-216. DOI:10.7505/j.issn.1007-9084.2013.02.017.

[26] 張建英. 飲用水氯化消毒副產(chǎn)物污染控制技術(shù)及健康風(fēng)險(xiǎn)評(píng)價(jià)的研究[D]. 杭州: 浙江大學(xué), 2006: 52-53.

[27] 金水高. 中國居民營養(yǎng)與健康狀況調(diào)查報(bào)告之十: 2002營養(yǎng)與健康狀況數(shù)據(jù)集[M]. 北京: 人民衛(wèi)生出版社, 2008: 347-378.

[28] WILLIAMS G M. DNA reactive and epigenetic carcinogens[J].Experimental and Toxicologic Pathology, 1992, 44: 457-463.

[29] ZHANG J, ZHENG N, LIU J, et al. Aflatoxin B1and aflatoxin M1induced cytotoxicity and DNA damage in differentiated and undifferentiated Caco-2 cells[J]. Food and Chemical Toxicology, 2015,83: 54-60. DOI:10.1016/j.fct.2015.05.020.