Monte Carlo模擬對面制品中鋁膳食暴露風(fēng)險的概率評估

劉瓊瑜,李浩洋,容裕棠,黃碧嘉,連海天,何秀琴,李 蓉

(中山海關(guān)(原中山出入境檢驗檢疫局),廣東中山 528403)

鋁是地殼中含量最多的金屬元素,不屬于機體的必需微量元素。大量的流行病學(xué)調(diào)查和動物實驗表明,長期過量鋁暴露會對人體產(chǎn)生慢性毒作用,包括神經(jīng)毒性[1-2]、生殖毒性[3-4]和肝臟毒性[5]等。食物中的鋁主要來源于食品加工生產(chǎn)過程中使用的含鋁食品添加劑以及食物加工和存儲過程與含鋁容器的接觸等[6]。世界糧農(nóng)組織(FAO)/世界衛(wèi)生組織食品添加劑聯(lián)合專家委員會(JECFA)已于1989年正式將鋁確定為食品污染物加以管理,將鋁的暫定每周耐受攝入量(provisional tolerable weekly intake,PTWI)定為7 mg/(kg·bw),2011年將其調(diào)整為2 mg/(kg·bw)[6]。面制品在我國居民的膳食結(jié)構(gòu)中占有較大比重,而油條等面制品屢屢被檢出鋁殘留量超標(biāo)[7-8],因此有必要對面制品來源的鋁膳食暴露風(fēng)險進(jìn)行評估。

目前,關(guān)于面制品中鋁的膳食暴露評估大部分采用簡單的點評估方法[9-10],即通過一個確定的食品消費量(平均值或某一百分位消費量數(shù)據(jù))乘以殘留物含量(一般是平均值),最后得到一個單一的膳食暴露量值[11]。點評估假定了化學(xué)物濃度恒定、居民對膳食消費量水平也相同,而在現(xiàn)實生活中往往不是這樣。近年來,基于Monte Carlo模擬的非參數(shù)概率評估模型開始應(yīng)用于重金屬、農(nóng)藥殘留、毒素的風(fēng)險評估中,能夠反映評估結(jié)果的變異度和不確定性[12-14],評估結(jié)果更加科學(xué)。本研究首次將概率評估應(yīng)用于面制品中鋁的膳食暴露評估,在對常見市售面制品中的鋁殘留量進(jìn)行檢測分析后,結(jié)合居民面制品消費量數(shù)據(jù)和人口學(xué)數(shù)據(jù),基于Monte Carlo模擬和Bootstrap抽樣技術(shù),采用風(fēng)險分析軟件Crystal Ball構(gòu)建了一個非參數(shù)概率評估模型,對居民面制品來源的鋁膳食暴露量進(jìn)行評估,并對其產(chǎn)生的暴露風(fēng)險進(jìn)行描述,以期為相關(guān)食品的安全監(jiān)管提供參考。

1 材料與方法

1.1 材料與儀器

本次調(diào)查的樣品來源于中山市的超市、農(nóng)貿(mào)市場、個體攤檔、餐館、單位食堂等生產(chǎn)流通領(lǐng)域，為了方便統(tǒng)計評價,將樣品根據(jù)其制作和食用方法分為4類,其中,膨化食品包括薯條、薯片、玉米片等,油炸類面制品包括油炸方便面、油條、炸饅頭、麻花等,西式糕點、餅干類焙烤食品包括面包、糕點、餅干等,蒸煮類面制品包括面條、中式包子及點心等；鋁單元素標(biāo)準(zhǔn)溶液 100 μg/mL,中國計量科學(xué)研究院;硝酸 優(yōu)級純,賽默飛世爾科技(中國)有限公司;過氧化氫 優(yōu)級純,廣州化學(xué)試劑廠;所有用水均為超純水。

ETHOS1型微波消解系統(tǒng) 意大利Milestone公司;Optima 5300DV型電感耦合等離子體發(fā)射光譜儀(ICP-AES) 美國PerkinElmer公司;EHD36型加熱消解裝置 北京萊伯泰科公司;超純水系統(tǒng) 德國默克公司。

1.2 實驗方法

1.2.1 鋁殘留量的測定 按照國家標(biāo)準(zhǔn)GB 5009.182-2017《食品安全國家標(biāo)準(zhǔn) 食品中鋁的測定》[15]中微波消解-電感耦合等離子體發(fā)射光譜法測定面制品中鋁的含量,方法定量限為10 mg/kg,對未檢出的樣品,根據(jù)世界衛(wèi)生組織和美國環(huán)境保護(hù)署推薦的數(shù)據(jù)處理方法,按定量限的一半(5 mg/kg)列入計算[16]。

1.2.1.1 樣品前處理 樣品粉碎后,準(zhǔn)確稱取0.500 g于微波消解內(nèi)罐中,加入5 mL硝酸,浸泡10 min,再加入3 mL過氧化氫,旋緊罐蓋后放入微波消解系統(tǒng)內(nèi),按表1所示消解條件進(jìn)行消解。消解結(jié)束,冷卻后取出,用少量水沖洗內(nèi)蓋,放在加熱消解裝置上于120 ℃加熱趕酸,最后用水定容至50 mL,同時做試劑空白。

表1 微波消解系統(tǒng)工作條件

1.2.1.2 電感耦合等離子體發(fā)射光譜儀操作條件 射頻功率1300 W,等離子體氣流量15.0 L/min,輔助氣流量0.20 L/min,蠕動泵轉(zhuǎn)速1.50 mL/min,霧化器流量0.80 L/min,工作頻率40.68 MHz,積分時間5 s,特征波長396.153 nm。鋁標(biāo)準(zhǔn)工作曲線系列濃度為0.1、0.2、0.4、0.6、0.8、1.0 μg/mL,線性相關(guān)系數(shù)在0.999以上。

1.2.2 暴露評估模型的構(gòu)建 參照J(rèn)ECFA制定的暫定每周耐受攝入量(PTWI)2 mg/(kg·bw),采用日暴露量(EXP)來評估面制品來源的鋁膳食暴露水平。

式中:EXP為日暴露量(μg/(kg·d));C為面制品中鋁的含量(mg/kg);DI為居民平均每標(biāo)準(zhǔn)人日面及其制品的攝入量(kg/d);BW為體重(kg)。

上述公式中的面及其制品攝入量和體重等人口學(xué)數(shù)據(jù)來源于《中國居民營養(yǎng)與健康狀況監(jiān)測2010～2013年綜合報告》[17],2010～2012年我國居民平均每標(biāo)準(zhǔn)人日面及其制品的攝入量為142.8 g/d,選取青少年、青年、中年、老年4類人群進(jìn)行暴露評估,具體體重信息如表2,本研究假設(shè)居民的體重數(shù)據(jù)滿足正態(tài)分布,記為N(μ,σ),μ為均值,σ為標(biāo)準(zhǔn)差。

表2 評估人群的體重信息

暴露評估模型的構(gòu)建使用一款基于Monte Carlo模擬技術(shù)的風(fēng)險分軟件——Crystal Ball。上述公式中的每個變量在Crystal Ball軟件中都需定義符合一種分布,該分布可以通過擬合度檢驗得到。每次Monte Carlo模擬都會從定義的變量概率分布中隨機選擇數(shù)值,然后計算模型結(jié)果,所有的結(jié)果構(gòu)成一個新的分布,Monte Carlo模擬給出該分布的均值、百分位數(shù)等統(tǒng)計量。本研究定義每次模擬過程循環(huán)10000次,并以檢測得到的910個面制品樣本中鋁的殘留量數(shù)據(jù)作為經(jīng)驗樣本,進(jìn)行1000次Bootstrap抽樣,每一次抽樣得到的數(shù)據(jù)再執(zhí)行10000次Monte Carlo模擬,以分析暴露評估結(jié)果的不確定性。

1.2.3 風(fēng)險描述 目前,食品中污染物的風(fēng)險評估經(jīng)常用到風(fēng)險商(HQ,hazard quotient)這個參數(shù)[18-19],本研究基于軟件模擬得到的鋁膳食日暴露量,以JECFA制定的暫定每周耐受攝入量(PTWI)為標(biāo)準(zhǔn),換算為每日可耐受攝入量ADI(Acceptable daily intake),為285.7 μg/(kg·d),通過計算面制品來源的鋁膳食日暴露風(fēng)險商(HQ),表征食用面制品中鋁可能導(dǎo)致的鋁暴露風(fēng)險。風(fēng)險商的計算公式如下:

式中,HQ為面制品來源的鋁膳食暴露風(fēng)險指數(shù),當(dāng)HQ≤1時,說明風(fēng)險在可接受范圍內(nèi);HQ>1,說明存在暴露風(fēng)險。

1.2.4 敏感度分析 Crystal Ball軟件的敏感度分析功能用于分析每個假設(shè)變量對一個特定預(yù)測值的影響程度。為了解面制品中的鋁殘留量、鋁及其制品的攝入量及體重對鋁日暴露量的影響程度,在設(shè)置Monte Carlo模擬參數(shù)時,勾選敏感度分析選項,模擬運行時,敏感度分析會評估每個變量對鋁膳食日暴露量的影響,并將影響程度表示為百分比值的形式,百分比值越大,影響力越大。

1.3 數(shù)據(jù)處理

采用SPSS 19.0對不同種類面制品的鋁殘留量數(shù)據(jù)進(jìn)行統(tǒng)計分析,檢驗水準(zhǔn)α=0.05(雙側(cè)),若p<0.05,則認(rèn)為具有統(tǒng)計學(xué)意義。

2 結(jié)果與分析

2.1 不同類別面制品中鋁殘留量檢測結(jié)果

對膨化食品,油炸面制品,西式糕點、餅干類焙烤食品,蒸煮類面制品4類面制品共910份樣品中鋁的殘留量進(jìn)行檢測,所有數(shù)據(jù)經(jīng)過SPSS非參數(shù)檢驗,結(jié)果顯示,全部樣品的鋁殘留量數(shù)據(jù)呈現(xiàn)偏態(tài)分布(p<0.05),鋁殘留量范圍在ND～1437.7 mg/kg之間,平均值為(15.5±79.4) mg/kg,所有樣品的平均檢出率為12.9%,平均超標(biāo)率為8.5%。不同類別面制品中鋁的殘留量差異具有統(tǒng)計學(xué)意義(χ2=9.585,p=0.022),其中,油炸類面制品中鋁的殘留量均值為(22.7±101.9) mg/kg,是4類面制品中最高的,最高檢測值達(dá)到了714.0 mg/kg,相當(dāng)于鋁殘留限量的71.4倍,與楊志瓊等[20]、焦海濤等[21]的調(diào)查結(jié)果一致;其次是蒸煮類面制品,鋁殘留量均值為(21.3±111.8) mg/kg(見表3)。按(JECFA)規(guī)定的鋁暫定每周耐受攝入量(PTWI)2 mg/(kg·bw)[6],即每日可耐受攝入量(ADI)285.7 μg/(kg·bw),中國成年人體重按60 kg計算,得到每人每日允許攝入的鋁量為17.14 mg,由2010～2012年我國居民平均每標(biāo)準(zhǔn)人日面及其制品的攝入量142.8 g/d[17],計算得到每人每日允許攝入的面制食品鋁殘留量最大為120 mg/kg,略高于國標(biāo)規(guī)定的鋁殘留最大限量100 mg/kg,但這只是針對面制食品,考慮到飲用水、鋁制炊具、含鋁藥物、天然食物中的鋁,以及居民接觸到鋁殘留超標(biāo)的面制食品的可能性,居民存在攝入鋁量超標(biāo)的可能性,關(guān)鍵在于杜絕面制食品中的鋁殘留超標(biāo)。

表3 各類面制品中鋁殘留量檢測結(jié)果

研究發(fā)現(xiàn),4類面制品,油炸類面制品的鋁殘留量均值(22.7±101.9) mg/kg和檢出率(20.0%)都最高,表明在售油條制作過程中普遍使用的膨松劑可能還是明礬,其主要成分為硫酸鋁鉀或硫酸鋁銨,含鋁量高,價格低廉,效果明顯,但用量多靠個人經(jīng)驗,很容易造成鋁超標(biāo)。據(jù)文獻(xiàn)報道,2013年四川省南充市面制品中鋁含量的均值為66.13 mg/kg,總超標(biāo)率為22.0%[8];2006～2011年廈門市面制品中鋁的平均含量為162.89 mg/kg[9];2009～2013年廣西市面制品中鋁的平均含量為157.81 mg/kg,鋁超標(biāo)率為35.29%[10]。相比較,本研究全部面制品的鋁殘留量均值(15.5 mg/kg)、平均超標(biāo)率(8.5%)較低,可能得益于兩方面原因:一是無鋁復(fù)合膨松劑如酵母等的推廣使用;二是食品監(jiān)管部門對含鋁食品添加劑的有力監(jiān)管,大大減少了含鋁食品添加劑的濫用。

運用Crystal Ball軟件對910個鋁殘留量數(shù)據(jù)樣本進(jìn)行擬合優(yōu)度檢驗,統(tǒng)計檢驗方法由軟件自動選擇,結(jié)合SPSS非參數(shù)檢驗的Difference圖、P-P圖、Q-Q圖確定最佳擬合分布為Lognormal分布(圖1),記為Lognormal(7.91,6.02)。該概率分布的具體參數(shù)如表4。由圖1可知,全部樣品的鋁殘留量數(shù)據(jù)大部分落在方法定量限水平附近或以下,這與SPSS非參數(shù)檢驗的偏態(tài)分布結(jié)果一致。

圖1 面制品中鋁殘留量的擬合優(yōu)度檢驗結(jié)果

表4 面制品中鋁殘留量最佳擬合分布的相關(guān)參數(shù)

2.2 暴露評估

本研究采用概率評估方法,結(jié)合《中國居民營養(yǎng)與健康狀況監(jiān)測2010～2013年綜合報告》中不同人群的體重數(shù)據(jù),利用Monte Carlo模擬技術(shù)從面制品的鋁殘留量和體重數(shù)據(jù)中隨機抽樣,模擬計算得到不同人群面制品來源的鋁膳食暴露量的概率分布及不同百分位數(shù),反應(yīng)了鋁殘留量和體重差異造成的暴露評估結(jié)果的變異性;通過1000次Bootstrap抽樣,量化分析了暴露評估結(jié)果的不確定性,使評估結(jié)果更科學(xué)可靠。

由表5可知,(6～17歲)青少年女性面制品來源的鋁膳食日暴露量最高,一般消費量人群(均值)和高消費量人群(日暴露量的P99百分位)對面制品來源的鋁膳食暴露量分別28.16(27.53～28.85)和116.47(103.96～129.23)。(6～17歲)青少年中,男性和女性鋁膳食日暴露量的平均值為26.25(25.59～26.83)和28.16(27.53～28.85),女性較男性要高,女性面臨的鋁暴露風(fēng)險高于男性,因為在同一年齡段,女性的體重一般小于男性,其他3組評估人群也有類似的特點。比較各組不同消費量人群的日暴露量,發(fā)現(xiàn)隨著消費量水平的提高,日暴露量水平明顯提高,不同評估人群鋁膳食日暴露量的P99百分位數(shù)分別為青少年(6～17歲):男性103.89/女性116.47;青年(18～44歲):男性66.36/女性79.36;中年(45～59歲):男性68.23/女性72.41;老年(≥60歲):男性73.80/女性78.89 μg/(kg·d),均低于每日可耐受攝入量(ADI)285.7 μg/(kg·d),說明居民面制品來源的鋁膳食暴露風(fēng)險在可接受范圍內(nèi)。

表5 不同人群面制品來源的鋁膳食暴露評估結(jié)果

2.3 風(fēng)險評估

基于前面模擬得到的鋁膳食日暴露量,以上8類人群面制品來源的鋁膳食暴露風(fēng)險指數(shù)如表6所示。一般消費量人群的日暴露量以均值表示,由表6可知,居民面制品來源的鋁膳食暴露風(fēng)險指數(shù)范圍為0.059～0.099,高消費量人群(日暴露量的P99百分位)的風(fēng)險指數(shù)范圍為0.23～0.41,所有HQ值均小于1,表明鋁膳食暴露風(fēng)險在可接受范圍內(nèi),普通居民通過面制品攝入鋁的健康風(fēng)險較低,由于(6～17歲)青少年的體重比成年人輕,受到的鋁膳食暴露風(fēng)險較大,可加強合理膳食引導(dǎo),避免過多攝入高鋁食品。

表6 不同人群面制品來源的鋁膳食暴露風(fēng)險指數(shù)(HQ)

2.4 敏感度分析

本研究對面制品來源的鋁膳食暴露量進(jìn)行了敏感性分析,結(jié)果見表7。從表7可以看出,面制品中鋁殘留量的敏感性最強,敏感值在所有評估人群中達(dá)到了95%以上,可見決定面制品來源的鋁膳食日暴露量大小的主要因素為面制品中的鋁殘留量,因此,要降低面制品來源的鋁膳食暴露量,應(yīng)該避免過多食用鋁殘留量高的面制品。體重的敏感值為負(fù)數(shù),表明體重對鋁膳食日暴露量起負(fù)面的影響,即體重越大,通過面制品攝入鋁的日暴露量越低。

表7 面制品來源的鋁膳食日暴露量的敏感度分析

3 結(jié)論

檢測的910份面制品樣品中,鋁殘留量范圍在ND～1437.7 mg/kg之間,平均值為15.5 mg/kg,平均檢出率和超標(biāo)率分別為12.9%和8.5%。不同類別面制品的鋁殘留量差異具有統(tǒng)計學(xué)意義,其中,以油炸類面制品的鋁殘留量均值和檢出率最高。

即使是在不同評估人群的高消費量水平(日暴露量的P99百分位),風(fēng)險指數(shù)最大值為0.41,小于1。鋁膳食暴露風(fēng)險在可接受范圍內(nèi),居民通過面制品攝入鋁的健康風(fēng)險較低。青少年(6～17歲)由于體重較輕,有一定潛在風(fēng)險,可對這一敏感群體加以關(guān)注,加強膳食合理引導(dǎo)。

面制品中鋁殘留量的敏感性最強(95%以上),是決定面制品來源的鋁膳食日暴露量大小的主要因素。體重的敏感值為負(fù)數(shù),對鋁膳食日暴露量起負(fù)面的影響,即體重越大,通過面制品攝入鋁的日暴露量越低。