楊蓉，李娜，饒凱鋒，馬梅,*，王子健

1.北京市南水北調(diào)水質(zhì)監(jiān)測(cè)中心，北京100195

2.中國(guó)科學(xué)院生態(tài)環(huán)境研究中心中國(guó)科學(xué)院飲用水科學(xué)與技術(shù)重點(diǎn)實(shí)驗(yàn)室，北京100085

環(huán)境混合物的聯(lián)合毒性研究方法

楊蓉1,2，李娜2，饒凱鋒2，馬梅2,*，王子健2

1.北京市南水北調(diào)水質(zhì)監(jiān)測(cè)中心，北京100195

2.中國(guó)科學(xué)院生態(tài)環(huán)境研究中心中國(guó)科學(xué)院飲用水科學(xué)與技術(shù)重點(diǎn)實(shí)驗(yàn)室，北京100085

人類和其他生物通過(guò)多重途徑暴露于大量的環(huán)境污染物中。單一化合物由于其在環(huán)境中非常低的濃度使得人們忽視其可能的健康危害,然而,多種化合物在環(huán)境中同時(shí)存在,就意味著更加復(fù)雜的聯(lián)合作用的可能性。Something from nothing現(xiàn)象引起了科學(xué)界對(duì)于環(huán)境混合物加和作用和潛在協(xié)同作用的關(guān)注,致使許多研究聯(lián)合毒性的數(shù)學(xué)和生物模型的提出和應(yīng)用。本文綜述了這些研究方法,以及未來(lái)化學(xué)混合物研究和管理面臨的挑戰(zhàn)。

混合物；聯(lián)合作用；作用方式；模型預(yù)測(cè)

人類和其他生物終生暴露在不同來(lái)源、不同種類的多種環(huán)境污染物中,單一污染物環(huán)境暴露在實(shí)際環(huán)境中并不存在。檢測(cè)各污染物的環(huán)境濃度水平、并評(píng)價(jià)其毒性和健康風(fēng)險(xiǎn)一直是科學(xué)界和大眾的關(guān)注焦點(diǎn)。傳統(tǒng)方法中通常是通過(guò)對(duì)目標(biāo)化合物的分析定量測(cè)定其環(huán)境濃度,并根據(jù)其毒性數(shù)據(jù)判斷其潛在的健康影響。雖然儀器和方法學(xué)的發(fā)展帶來(lái)了檢測(cè)限的降低和更多痕量物質(zhì)的發(fā)現(xiàn),但儀器分析檢測(cè)到的目標(biāo)化合物種類有限,可能只占種類繁多的環(huán)境污染物中很少一部分,會(huì)導(dǎo)致我們忽略了未知污染物的存在和毒性效應(yīng)。例如,Escher等[1]研究了包括回用水、飲用水、雨水在內(nèi)的9個(gè)環(huán)境水樣造成的氧化應(yīng)激反應(yīng),指出已知化合物的貢獻(xiàn)小于0.1%,而超過(guò)99.9%的氧化應(yīng)激都由未知化合物造成。另一方面,環(huán)境污染及其健康影響的調(diào)查與評(píng)價(jià),其最終目的是保護(hù)生態(tài)系統(tǒng)和人類健康。從這個(gè)角度考慮,污染物濃度本身并沒(méi)有實(shí)質(zhì)意義,但是其所產(chǎn)生的生物毒性或?qū)】档母蓴_效應(yīng)才是應(yīng)該受到關(guān)注的重點(diǎn)。因此在生態(tài)和健康風(fēng)險(xiǎn)評(píng)估中,暴露評(píng)價(jià)和危害評(píng)價(jià)均是整體研究中的不可或缺的部分。

毒理學(xué)研究中的毒性測(cè)試則關(guān)注整體生物效應(yīng),通常是測(cè)定某個(gè)具體污染物或混合物樣品針對(duì)具體受試生物和測(cè)試終點(diǎn)的毒性表現(xiàn)。在20世紀(jì),絕大多數(shù)的毒理學(xué)研究和風(fēng)險(xiǎn)評(píng)價(jià)都停留在單一化合物的水平上[2],因?yàn)獒槍?duì)單一物質(zhì)的毒性研究花費(fèi)低廉且簡(jiǎn)單易行。然而,已有大量證據(jù)證明混合物的聯(lián)合毒性并不能簡(jiǎn)單地通過(guò)組分濃度或毒性當(dāng)量加合來(lái)描述,而是表現(xiàn)出相互影響、相互作用,會(huì)產(chǎn)生諸如加和、協(xié)同或拮抗等聯(lián)合毒作用。此外,聯(lián)合毒性還會(huì)隨著組分濃度水平、濃度比、暴露時(shí)間、路徑和順序、以及測(cè)試生物的狀態(tài)發(fā)生變化。由此可見(jiàn),混合物的聯(lián)合毒性研究是一個(gè)挑戰(zhàn)性與必要性并存的課題?；趩我晃镔|(zhì)的毒性研究和管理策略很難用來(lái)解決復(fù)雜的污染問(wèn)題,而研究聯(lián)合效應(yīng)對(duì)進(jìn)行風(fēng)險(xiǎn)評(píng)價(jià)、改進(jìn)環(huán)境標(biāo)準(zhǔn)和環(huán)境容量、開(kāi)展預(yù)防評(píng)價(jià)對(duì)策都有很強(qiáng)的推動(dòng)作用和現(xiàn)實(shí)意義。

對(duì)聯(lián)合作用的描述及其評(píng)價(jià)方法的探討已日益引起研究人員和環(huán)境管理人員的重視,提出了許多基于生理學(xué)、毒理學(xué)和生物統(tǒng)計(jì)學(xué)的方法模型并加以應(yīng)用優(yōu)化。近年來(lái),美國(guó)環(huán)保署(United States Environmental Protection Agency,USEPA)、毒物與疾病登記署(Agency for Toxic Substances and Disease Registry,ATSDR)、經(jīng)濟(jì)合作與發(fā)展組織(Organization for Economic Cooperation and Development, OECD)和世界衛(wèi)生組織(World Health Organization, WHO)等[3-7]一系列具有影響力的環(huán)境機(jī)構(gòu)和組織也開(kāi)始發(fā)布混合物毒性研究的準(zhǔn)則和標(biāo)準(zhǔn)。本文綜述了這一領(lǐng)域的進(jìn)展和主要成就,包括用來(lái)預(yù)測(cè)聯(lián)合毒性的數(shù)學(xué)和生物學(xué)模型,以及未來(lái)化學(xué)混合物研究和管理面臨的挑戰(zhàn)。

1 聯(lián)合作用的分類(Classification of combined effects)

Bliss[8]于1939年發(fā)表的論文標(biāo)志著混合污染物聯(lián)合毒性研究的開(kāi)始,其中首次提出了聯(lián)合作用的分類方法:獨(dú)立作用(independent joint action)、相似作用(similar joint action)和協(xié)同作用(synergistic action)。Bliss強(qiáng)調(diào)了聯(lián)合作用分類的關(guān)鍵之處,即一種組分的存在是否干擾混合物中其他組分的效果。由此,聯(lián)合作用可被分為無(wú)相互作用(non-interative,包括獨(dú)立作用和相似作用)和相互作用(interactive,此處的“協(xié)同作用”包含了傳統(tǒng)意義上的協(xié)同和拮抗)2種。在無(wú)相互作用的進(jìn)一步分類中,也考慮了混合物組分的作用位點(diǎn)和路徑是否相同。

表1 Bliss的聯(lián)合作用分類法Table 1 Classification of combined effects by Bliss

表1說(shuō)明Bliss僅將協(xié)同作用定義為對(duì)加和的偏離,但并未對(duì)此開(kāi)展進(jìn)一步討論。在此基礎(chǔ)上, Plackett和Hewlett[9-10]根據(jù)作用方式(mode of action, MoA)的不同,對(duì)相互作用做了進(jìn)一步分類。復(fù)雜相似作用(complex similar action)和關(guān)聯(lián)作用(dependent action)兩個(gè)概念被分別用來(lái)描述機(jī)理相似和機(jī)理不同的混合物組分的相互作用(見(jiàn)表2)。然而,Plackett和Hewlett的分類法仍存在不足之處,它只適用二元混合物,不能用于多組分混合物、尤其是復(fù)雜環(huán)境樣品的分類。因此,Ashford[11]在總結(jié)了以上分類法的特性和優(yōu)缺點(diǎn)后,將各種聯(lián)合作用基于亞系統(tǒng)、作用位點(diǎn)和混合物組分重新分類,提出了將組分相似性、關(guān)聯(lián)性、一致性等信息同時(shí)納入框架的新方法(見(jiàn)表3)。與前述分類法不同,新引入的部分相似作用(partially similar)和部分關(guān)聯(lián)作用(partially dependent)2個(gè)概念更加符合環(huán)境樣品的真實(shí)情況,也為聯(lián)合作用提供了更加細(xì)致的描述。

表2 Plackett和Hewlett的聯(lián)合作用分類法Table 2 Classification of combined effects by Plackett and Hewlett

在Ashford分類法提出的同一時(shí)期,WHO和USEPA也發(fā)布了各自推薦的聯(lián)合作用分類,主要參考了半世紀(jì)前Bliss分類法的思路。WHO[12]將聯(lián)合作用分成獨(dú)立作用(independent joint action)、加和作用(additive effects)、協(xié)同作用(synergistic joint action)和拮抗作用(antagonistic joint action)4類；USEPA[5]則將最后兩者合并,分為劑量加和(dose addition)、效應(yīng)加和(response addition)和相互作用(interaction)3類。30多年來(lái),WHO和USEPA的巨大影響力使其推薦的分類法被普遍采用。但這也帶來(lái)了以下問(wèn)題:作為自然界中廣泛存在的現(xiàn)象,相互作用的研究不應(yīng)當(dāng)僅僅滿足于諸如“偏離加和”、“協(xié)同(大于加和)”和“拮抗(小于加和)”的現(xiàn)象描述,而應(yīng)該更加深入到作用方式和機(jī)理的探討。因此,雖然WHO和EPA對(duì)本領(lǐng)域中關(guān)鍵概念的推廣和標(biāo)準(zhǔn)化作出了很大貢獻(xiàn),但他們的觀點(diǎn)并不能代表混合物毒性研究的最新進(jìn)展。

圖1 基于數(shù)據(jù)質(zhì)量的混合物風(fēng)險(xiǎn)評(píng)價(jià)Fig.1 Flow chart for evaluating mixtures based on data quality

表3 Ashford的聯(lián)合作用分類法Table 3 Classification of combined effects by Ashford

2 無(wú)相互作用模型(Non-interactive combined effects)

USEPA于1986年發(fā)布的化學(xué)混合物風(fēng)險(xiǎn)評(píng)價(jià)指南是混合物毒性研究領(lǐng)域的重要文檔。其針對(duì)混合物毒性研究提出了3個(gè)層次的問(wèn)題,以指導(dǎo)不同數(shù)據(jù)質(zhì)量的混合物進(jìn)行研究方法的選擇[5-6]:(1)是否是高質(zhì)量數(shù)據(jù)?(2)是否有完全混合物數(shù)據(jù)?(3)混合物組分信息是否已知?因此,可以由數(shù)據(jù)質(zhì)量的評(píng)價(jià)開(kāi)始,評(píng)估已有數(shù)據(jù)并選擇合適的風(fēng)險(xiǎn)評(píng)價(jià)方法,如Teuschler[13]的流程圖所示(圖1)。完全混合物及充分相似混合物是把混合物看作一個(gè)整體加以評(píng)價(jià),在混合物組分信息已知的情況下,可以通過(guò)對(duì)其組分的分析來(lái)評(píng)價(jià)混合物,這更利于聯(lián)合毒性的分析和預(yù)測(cè)、以及作用模式機(jī)理和相互作用的探討。在基于組分的預(yù)測(cè)方法中,大多數(shù)研究聯(lián)合效應(yīng)的模型都與2個(gè)基本概念相關(guān),是由Loewe和Bliss分別提出的濃度加和模型(concentration addition,CA)和效應(yīng)加和模型(response addition,RA)。在相互作用不存在或組分濃度低到可以忽略相互作用的情況下,濃度加和和效應(yīng)加和都被USEPA推薦為風(fēng)險(xiǎn)評(píng)價(jià)的通用方案,方法簡(jiǎn)便可靠,應(yīng)用廣泛[8,14]。這兩者的區(qū)別在于混合物每一組分是否具有毒理學(xué)相似性。混合物各組分可能具有相同或不同的作用方式,結(jié)合相同或不同的作用位點(diǎn)(target of action, ToA),這都會(huì)影響對(duì)評(píng)價(jià)方法和預(yù)測(cè)模型的選擇。

濃度加和也被稱為劑量加和[5]和簡(jiǎn)單相似作用[8]。它基于混合物組分共享一個(gè)作用方式的假設(shè),認(rèn)為在混合物各組分結(jié)構(gòu)類似、生理功能和毒理行為相似、只在毒性效力上有所區(qū)別、某組分與其他組分無(wú)相互作用。因此,在調(diào)整各組分的毒性強(qiáng)度后,我們可以將各單一化合物的劑量(或濃度)相加,由此計(jì)算聯(lián)合效應(yīng)[15]。嚴(yán)格地說(shuō),濃度加和只可用于所有組分競(jìng)爭(zhēng)結(jié)合同一分子位點(diǎn)的情況[16],但這一苛刻的條件可以適當(dāng)放寬至各組分可以引起相同生物響應(yīng)[17]。濃度加和說(shuō)明即使各化合物濃度低于其各自無(wú)觀察效應(yīng)濃度(no-observed effect concentration,NOEC),也可能會(huì)產(chǎn)生顯著的聯(lián)合作用。其可表示為:

其中,pi表示i化合物濃度在混合物總濃度中的比例,ECx,i表示i化合物產(chǎn)生x%效應(yīng)時(shí)的濃度,比如各化合物的半數(shù)效應(yīng)濃度(half-effective concentration,EC50),ECx,mix表示混合物產(chǎn)生x%效應(yīng)時(shí)的濃度[18]。顯然,不論單個(gè)化合物的濃度是否低于無(wú)觀察效應(yīng)濃度,混合物中每個(gè)化合物對(duì)總體的毒效應(yīng)均有貢獻(xiàn)。

濃度加和是混合物風(fēng)險(xiǎn)評(píng)價(jià)中使用最多的方法,它具有簡(jiǎn)單的數(shù)學(xué)表達(dá)式,方便易用,傾向在某種程度上高估了聯(lián)合作用,這有利于潛在風(fēng)險(xiǎn)的識(shí)別、符合預(yù)防原則。并且當(dāng)復(fù)雜混合物的毒性信息不完全已知時(shí),濃度加和可以被用作通用解決方法,因?yàn)槟Ｐ皖A(yù)測(cè)和實(shí)驗(yàn)觀察到的結(jié)果總體而言相差不大[19-20]。由這個(gè)方法衍生出來(lái)的方法眾多,包括危害指數(shù)法(hazard index,HI)[21-22]、毒性單位法(toxic unit,TU)[23-24]、相對(duì)強(qiáng)度因子法(relative potency factor,RPF)[6]、毒性當(dāng)量因子法(toxicity equivalence factor,TEF)[25]、總體負(fù)荷法(total body residue,TBR)[26]等,需要根據(jù)毒性機(jī)理和毒理學(xué)相似性方面的相關(guān)信息進(jìn)行選擇[15]。

濃度加和模型及其衍生方法在環(huán)境研究和管理中得到廣泛的應(yīng)用。有許多研究者利用濃度加和研究混合物毒性,受試生物包括細(xì)菌、酵母、細(xì)胞、嚙齒類和魚類,測(cè)試終點(diǎn)涵蓋生物發(fā)光[27]、酶活性[28-29]、細(xì)胞增殖[30-31]、卵黃蛋白原水平[32-33]、mRNA 表達(dá)[34]、子宮增重[35]等。另外,作為濃度加和的衍生物,標(biāo)準(zhǔn)化合物當(dāng)量法也被成功應(yīng)用于描述聯(lián)合作用。例如,用二噁英的TEF表示混合體系相對(duì)于2, 3,7,8-TCDD的芳烴受體結(jié)合效應(yīng)[25,36],用雌激素的EEF表示混合體系相對(duì)于雌二醇的雌激素受體結(jié)合效應(yīng)[37],用RPF表示有機(jī)磷農(nóng)藥混合物相對(duì)于甲胺磷的乙酰膽堿酯酶抑制效應(yīng)[38],而用BaP-TEF表示化合物的遺傳毒物相對(duì)于苯并[α]芘的致癌性[39]。

另一方面,效應(yīng)加和是每個(gè)混合物組分效應(yīng)的加和,也稱獨(dú)立作用[8,40]或簡(jiǎn)單非相似作用[15]。此時(shí),混合物中的化合物被認(rèn)為彼此獨(dú)立地與不同的系統(tǒng)或靶位點(diǎn)發(fā)生作用、不影響其他物質(zhì)的毒性。效應(yīng)加和通常被描述為:

其中,E(cmix)表示混合物的聯(lián)合毒效應(yīng),E(ci)表示i化合物單獨(dú)作用時(shí)的毒效應(yīng)。與濃度加和概念中的“不論濃度大小,所有組分均對(duì)總毒性有所貢獻(xiàn)”不同,效應(yīng)加和模型認(rèn)為各化合物的濃度必須在其環(huán)境安全水平以上才能對(duì)混合物的總毒性有所貢獻(xiàn)[41]。當(dāng)混合物中各化合物具有不同的作用方式、且每個(gè)組分的存在不會(huì)影響其他化合物毒性時(shí),效應(yīng)加和對(duì)聯(lián)合毒性有更好的預(yù)測(cè)能力[42]。

事實(shí)上,環(huán)境中化合物作用完全相同或完全不同的情況很少發(fā)生,濃度加和和效應(yīng)加和用于實(shí)際樣品都存在局限性。因此Junghans[43]在更接近實(shí)際混合物的探索中提出了整合濃度加和/效應(yīng)加和于一體的兩步預(yù)測(cè)(two-step prediction,TSP)模型,其表達(dá)式如下:

顯然,上式為濃度加和和效應(yīng)加和表達(dá)式的結(jié)合。作用方式的化合物先各自組成數(shù)個(gè)集合,分別用濃度加和預(yù)測(cè)集合內(nèi)效應(yīng)；之后,集合間聯(lián)合作用的再使用效應(yīng)加和預(yù)測(cè)。兩步預(yù)測(cè)模型已被應(yīng)用于污水處理廠出水的評(píng)價(jià),被證明可用于評(píng)價(jià)復(fù)雜混合物的效應(yīng)[44]。

無(wú)相互作用模型在實(shí)踐中得到廣泛應(yīng)用,很長(zhǎng)一段時(shí)間內(nèi)幾乎被用于混合物毒性描述和預(yù)測(cè)的通用解決方案。表4總結(jié)了使用濃度加和和效應(yīng)加和進(jìn)行的部分聯(lián)合毒性研究,顯然,混合物組分的作用方式和最優(yōu)模型選擇具有一定相關(guān)性。

表4 濃度加和和效應(yīng)加和模型在混合物毒理研究中的應(yīng)用Table 4 Application of CA and RA in mixture toxicity researches

然而,研究者發(fā)現(xiàn)現(xiàn)實(shí)世界里無(wú)相互作用的普適性并不令人滿意。例如,在針對(duì)雌激素混合物的研究中,采用 MCF-7細(xì)胞增殖[53-55]、海龜性別決定[56-57]和小鼠子宮增重[58]作為測(cè)試終點(diǎn),也都發(fā)現(xiàn)了協(xié)同或拮抗效應(yīng)。Vijver等[59]綜述了2010年之前發(fā)表的文獻(xiàn),發(fā)現(xiàn)模型預(yù)測(cè)是無(wú)相互作用的只占全部數(shù)據(jù)的13%。對(duì)有機(jī)有毒物質(zhì),加和模型有較好的預(yù)測(cè)效果,但不同金屬離子之間主要是協(xié)同和拮抗作用。由此可知,實(shí)際環(huán)境中存在的相互作用更加普遍。這就意味著預(yù)測(cè)相互作用需要有新的方法學(xué)突破,應(yīng)該重點(diǎn)考慮不同的靶點(diǎn)和不同的作用方式之間的交互影響[60]。理論上,這類聯(lián)合毒性的研究更需要數(shù)學(xué)模型、生理學(xué)、系統(tǒng)毒理學(xué)、計(jì)算機(jī)預(yù)測(cè)等多方面的支持與合作,目前處于發(fā)展階段。

3 常見(jiàn)的聯(lián)合毒性實(shí)驗(yàn)設(shè)計(jì)(Common experimental design for combined effects)

濃度加和效應(yīng)加和模型使用雖廣,但其最大的缺點(diǎn)在于它們忽視了組分之間可能的相互作用。相互作用是對(duì)加和性(additivity)的偏離的經(jīng)驗(yàn)性描述,會(huì)影響混合物的毒理學(xué)應(yīng)答和整體毒性,一般分為協(xié)同和拮抗效應(yīng)。許多研究者研究了混合物組分之發(fā)生的相互作用頻率及強(qiáng)度,證明在真實(shí)世界里相互作用是普遍存在的[59,61]。當(dāng)組分的濃度和濃度比發(fā)生變化時(shí),可能會(huì)有更微妙和復(fù)雜的相互作用發(fā)生[62]。因此,研究者提出了諸多用于研究聯(lián)合毒性、分辨加和、協(xié)同和拮抗的實(shí)驗(yàn)設(shè)計(jì)方法,在此予以綜述。

3.1 等效線圖

在特定效應(yīng)水平對(duì)二元混合物進(jìn)行分析時(shí),等效線圖是區(qū)分加和效應(yīng)和相互作用的經(jīng)典統(tǒng)計(jì)方法[14]。在此圖上,x軸和y軸分別是組分A和B的毒性單位,用混合物中該組分濃度與其單獨(dú)作用時(shí)特定效應(yīng)水平下濃度(一般是EC50)的比值表示(Ci/ EC50i)。當(dāng)混合物兩組分之間沒(méi)有相互作用,實(shí)驗(yàn)點(diǎn)會(huì)形成一條筆直的“加和線”,如圖2所示。如果實(shí)驗(yàn)點(diǎn)出現(xiàn)在加和線的上方和下方,即認(rèn)為有拮抗和協(xié)同作用發(fā)生。

圖2 等效線圖(I:加和效應(yīng)；II:協(xié)同效應(yīng)；III:拮抗效應(yīng))Fig.2 Isobologram(I:additivity;II:synergism;III:antagonism)

盡管等效線圖可以直觀地表示聯(lián)合效應(yīng),它也存在很多方面的局限:只可用于二元混合物在某一特定效應(yīng)水平上的評(píng)估,因此需要開(kāi)展大量工作來(lái)精確地測(cè)定每一組分的等效應(yīng)濃度[63]。然而,在試驗(yàn)方法學(xué)的發(fā)展下,某一效應(yīng)水平上的毒性數(shù)據(jù)——如EC50已經(jīng)不是評(píng)估和比較混合物毒性的唯一指標(biāo)。環(huán)境中化合物濃度隨空間和時(shí)間變化巨大,只研究單一濃度比的傳統(tǒng)方法難以反映環(huán)境現(xiàn)實(shí),而增加組分和濃度組的數(shù)量又會(huì)大大提高工作量。此外,實(shí)驗(yàn)設(shè)計(jì)方式和實(shí)驗(yàn)產(chǎn)生的誤差會(huì)很大程度地影響基于等效線設(shè)計(jì)而得出結(jié)論的置信度水平,需要在統(tǒng)計(jì)學(xué)方法的幫助下得到清晰、可控、數(shù)據(jù)使用充分的實(shí)驗(yàn)設(shè)計(jì)[64]。

3.2 混合物毒性指數(shù)法

混合物毒性指數(shù)法是通過(guò)公式計(jì)算定量判斷模型預(yù)測(cè)值和實(shí)際觀測(cè)值之間關(guān)系的方法。針對(duì)二元混合物,大多數(shù)指數(shù)法是由等效線圖推導(dǎo)得出的,與其代數(shù)等價(jià)。例如Finney[65]的K值法,定義K為實(shí)驗(yàn)得到的EC50和濃度加和預(yù)測(cè)到的EC50的比值。從理論角度講,K=1說(shuō)明毒物具有加和作用,K＜1或K＞1說(shuō)明它們有拮抗或協(xié)同作用。但由于實(shí)驗(yàn)誤差和實(shí)驗(yàn)動(dòng)物狀態(tài)的波動(dòng),在0.4到2.5之間的K值一般被認(rèn)為表示加和。之后,也有其他研究者推薦了不同的K值加和區(qū)間,或者提出了更完善的K值評(píng)價(jià)系統(tǒng)[66-67]。

其他可行并常用的指數(shù)法包括毒性單位、加和指數(shù)(additive index,AI)、混合毒性指數(shù)(mixture toxicity index,MTI)、相似參數(shù)(similarity parameter,λ)等,它們?cè)诮耆员挥糜诼?lián)合毒性評(píng)價(jià)[68-69],其定義式、特性、應(yīng)用等詳細(xì)闡述可參見(jiàn)曾鳴等[70]的綜述。與等效線圖相比,指數(shù)法可用于多元化合物的聯(lián)合毒性預(yù)測(cè),并使大量毒性數(shù)據(jù)最終以一個(gè)清晰、簡(jiǎn)明的數(shù)字呈現(xiàn),從而更受歡迎[71]。但是,指數(shù)法具有與等效線圖相似的缺點(diǎn),其僅適用于具有單調(diào)劑量-效應(yīng)曲線、EC50單指標(biāo)即可用來(lái)確定整個(gè)作用區(qū)域的化合物。若混合物由數(shù)種劑量-效應(yīng)曲線為U型或倒U型的組分構(gòu)成,其聯(lián)合毒性將不能用指數(shù)法計(jì)算。

3.3 析因分析

為避免傳統(tǒng)實(shí)驗(yàn)設(shè)計(jì)上的種種局限,研究者開(kāi)始采用統(tǒng)計(jì)學(xué)方法研究化學(xué)混合物的聯(lián)合毒性。其中,析因分析是一種判斷多組分體系是否存在相互作用的定性方法,它可以同時(shí)觀察不同因子的作用,分析因子之間的相互作用,提供各化合物各濃度間的數(shù)學(xué)關(guān)系。與等效線圖和指數(shù)法不同,析因設(shè)計(jì)能夠分析獲得的所有數(shù)據(jù)點(diǎn),并不要求提供各化合物精確的劑量-效應(yīng)方程,因此可以提高試驗(yàn)效率,更全面地利用實(shí)驗(yàn)數(shù)據(jù),減少工作量和實(shí)驗(yàn)動(dòng)物用量,得到高質(zhì)量的結(jié)果[72-73]。采用析因分析研究聯(lián)合毒性的思路為:(1)獲得各組分的單一劑量-效應(yīng)曲線；(2)選擇合適的濃度進(jìn)行混合物析因設(shè)計(jì)(2水平2因子是最簡(jiǎn)單的)；(3)通過(guò)雙因素方差分析并結(jié)合均數(shù)輪廓圖判斷混合物組分間有無(wú)交互作用。若估計(jì)邊緣均數(shù)輪廓圖中曲線互相平行,認(rèn)為混合物各組分間表現(xiàn)為相加,兩條曲線互相遠(yuǎn)離表示協(xié)同,曲線靠近或交叉為拮抗。例如,孫立偉等[33]利用析因分析研究了雌二醇、壬基酚、雙酚A所組成的二元混合物的雌激素作用,發(fā)現(xiàn)估計(jì)邊緣均數(shù)輪廓圖的所有線條基本平行,說(shuō)明組分間不存在相互作用。這一方法也可以擴(kuò)展到三元混合物的情況[74-75],以大幅度減少多組分、多組合、多濃度水平下混合物設(shè)計(jì)和毒性測(cè)試的工作量。

4 生物學(xué)和生理學(xué)模型(Biological and physiological models)

4.1 效應(yīng)終點(diǎn)與混合物作用靶點(diǎn)和作用方式

在毒理學(xué)研究中,決定用什么生物進(jìn)行暴露、暴露后檢測(cè)哪些指標(biāo),這些方面的考慮和研究何種類化合物或其混合物的毒性同樣重要。當(dāng)生物體暴露在環(huán)境樣品中,各種化合物通過(guò)飲食、呼吸和皮膚進(jìn)入機(jī)體和組織,通過(guò)血液、淋巴等組織液運(yùn)輸達(dá)到靶位點(diǎn)。通過(guò)與相應(yīng)的生物分子結(jié)合干擾正常生理過(guò)程。它們依次導(dǎo)致可逆的生理生化改變和不可逆的病理改變,并隨之產(chǎn)生病變、疾病和死亡。因此,作為環(huán)境混合物的作用對(duì)象,生物體是一個(gè)過(guò)于復(fù)雜的系統(tǒng)。化合物在不同機(jī)體層面上都可引起效應(yīng),有的還能在不同靶位點(diǎn)引起不同的效應(yīng),所以評(píng)價(jià)一個(gè)化合物毒性會(huì)有很多個(gè)測(cè)試終點(diǎn)可供選擇。

對(duì)同樣的受試物和同樣的測(cè)試生物,不同的毒理學(xué)終點(diǎn)具有不同的敏感性。例如,針對(duì)銅對(duì)布氏雙尾藻的生長(zhǎng)抑制作用,Rijstenbil等[76]發(fā)現(xiàn)當(dāng)所選擇的效應(yīng)終點(diǎn)不同時(shí)會(huì)獲得不同的劑量-效應(yīng)關(guān)系。也就是說(shuō),在給定的銅濃度下,親脂配體、有性繁殖、胞溶現(xiàn)象等效應(yīng)終點(diǎn)會(huì)表現(xiàn)出強(qiáng)度不同的響應(yīng)。而金小偉等[77]在評(píng)價(jià)壬基酚的生態(tài)風(fēng)險(xiǎn)時(shí),選擇存活、生長(zhǎng)、生化指標(biāo)以及繁殖作為測(cè)試終點(diǎn)會(huì)使得實(shí)驗(yàn)的靈敏性逐漸升高,相應(yīng)使預(yù)測(cè)無(wú)效應(yīng)濃度(predicted no effect concentration,PNEC)值逐漸降低。因此在毒性研究中測(cè)試終點(diǎn)的選擇是非常重要的,生物測(cè)試方法的開(kāi)發(fā)也正在囊括越來(lái)越多的毒理學(xué)效應(yīng)終點(diǎn)。

然而,污染物具有不同的致毒機(jī)制和作用方式。單一化合物進(jìn)入生物體內(nèi)都會(huì)引起不同位點(diǎn)、不同水平的反應(yīng),混合物、尤其是作用機(jī)制復(fù)雜的混合物導(dǎo)致的聯(lián)合效應(yīng)將由于ToA和MoA的不同而更難以分析和預(yù)測(cè)。例如在一個(gè)復(fù)雜混合物中,有機(jī)污染物和重金屬離子組分對(duì)于生物體的作用靶點(diǎn)可能相差甚大,這大大增加了研究混合物毒性的難度。傳統(tǒng)意義上對(duì)于化合物急性毒性的分類相對(duì)簡(jiǎn)單,主要是炎癥性損傷和反應(yīng)性損傷；而當(dāng)涉及慢性毒性或健康效應(yīng)會(huì)讓聯(lián)合毒性研究舉步維艱。例如,多氯聯(lián)苯和對(duì)硫磷都是典型有機(jī)污染物,但它們的致毒機(jī)理分別是抑制免疫系統(tǒng)和抑制膽堿酯酶；而作為重金屬,鉛的作用位點(diǎn)是神經(jīng)海馬細(xì)胞,鎘的主要作用位點(diǎn)是骨骼細(xì)胞,硒和鉈的主要作用位點(diǎn)則分別是膠原蛋白代謝和Na/K-ATP酶。在很多環(huán)境樣品毒性評(píng)價(jià)中,實(shí)驗(yàn)人員把測(cè)試生物暴露在復(fù)雜混合物中檢測(cè)急性毒性。但由于存在眾多的致死途徑,導(dǎo)致結(jié)果復(fù)雜和無(wú)法解釋。例如,一項(xiàng)4種抗生素聯(lián)合暴露魚腥藻和綠藻的研究發(fā)現(xiàn)聯(lián)合毒性呈現(xiàn)出協(xié)同效應(yīng)[78],但這一結(jié)果只針對(duì)此研究所用的混合物,這種結(jié)論對(duì)預(yù)測(cè)模型的開(kāi)發(fā)和外推至人類的風(fēng)險(xiǎn)評(píng)估幫助有限。因此,為了解釋復(fù)雜混合物的效應(yīng),研究者已經(jīng)把注意力轉(zhuǎn)到了基于作用機(jī)制的研究思路上。此時(shí)關(guān)注的重點(diǎn)是受體、作用方式和機(jī)理,而不是一種或一類化合物對(duì)整體動(dòng)物或細(xì)胞的毒性。從這一角度,“改變攝入和排泄過(guò)程的化合物”、“抑制某特定酶的化合物”等概念會(huì)比“酚類”、“重金屬”更易于在聯(lián)合毒性領(lǐng)域進(jìn)行分類和開(kāi)展研究。在Groten等[15]的綜述里,“簡(jiǎn)單混合物”被定義為10個(gè)或以下化合物的集合；但在新的框架下,它也許會(huì)表示包含了相同作用方式物質(zhì)的混合物,因?yàn)樗慕M分之間不存在相互作用,討論起來(lái)相對(duì)簡(jiǎn)單。隨著角度的變化,很多概念和思考方式都會(huì)不同。這個(gè)基本概念預(yù)示著未來(lái)的聯(lián)合毒性研究將更多側(cè)重于與毒理學(xué)研究同步[79]。在本實(shí)驗(yàn)室最近發(fā)表的2篇文章中,也強(qiáng)調(diào)了基于相同作用方式和作用靶點(diǎn)研究聯(lián)合毒性的重要性[48,80]。

4.2 基于過(guò)程的研究框架:對(duì)危害指數(shù)的修正

近年來(lái),系統(tǒng)毒理學(xué)方法越來(lái)越多地應(yīng)用于相互作用的判斷和預(yù)測(cè)中。從生物學(xué)過(guò)程的角度看,相互作用會(huì)在3個(gè)層次上發(fā)生[62]:一是生物體外的相互作用,排放到環(huán)境中的大多數(shù)污染物都會(huì)與某些環(huán)境成分或互相發(fā)生作用,以改變暴露水平和最終被生物體吸收的量。例如土壤、沉積物顆粒對(duì)可溶性腐殖質(zhì)的吸附和解吸,藥物和酸性物質(zhì)影響金屬的生物可利用性等[81-82]。二是化合物在體內(nèi)吸收、分配、代謝、排出(aborption-distribution-metabolism-excretion,ADME)過(guò)程中的相互作用(毒物代謝動(dòng)力學(xué))。當(dāng)組分B改變組分A的攝入和/或消除,解毒或激活作用會(huì)發(fā)生,引起組分A在靶位點(diǎn)的濃度改變。如丙磺舒可通過(guò)降低青霉素的腎臟排泄加強(qiáng)其藥效,和鋅可通過(guò)降低鉛的吸收抑制鉛中毒[83]。三是靶位點(diǎn)的相互作用(毒物作用動(dòng)力學(xué))。如果混合物各組分和相同位點(diǎn)競(jìng)爭(zhēng)或非競(jìng)爭(zhēng)性結(jié)合,組分B的存在將會(huì)影響組分A和靶位點(diǎn)的結(jié)合,影響其毒性。例如,阿托品可通過(guò)結(jié)合乙酰膽堿受體位點(diǎn)拮抗有機(jī)磷的毒性,他莫昔芬由于可與雌激素競(jìng)爭(zhēng)結(jié)合受體被用于雌激素陽(yáng)性乳腺癌的治療。后兩個(gè)層次在聯(lián)合效應(yīng)研究中更受重視,一些毒物代謝動(dòng)力學(xué)和毒物作用動(dòng)力學(xué)的修飾也被用來(lái)反映化學(xué)組份在體內(nèi)的濃度變化和相互作用。因此,有研究者使用生理藥代動(dòng)力學(xué)模型(physiologically-based pharmacokinetic model,PBPK)和證據(jù)加權(quán)法(weight of evidence,WOE)對(duì)危害指數(shù)HI進(jìn)行修正,使基于濃度加和的HI成為基于相互作用的IbHI(Interaction-based HI),這也是結(jié)合相互作用描述聯(lián)合毒性的方法之一。兩者分別對(duì)機(jī)體內(nèi)靶位點(diǎn)的暴露濃度和相互作用存在時(shí)相應(yīng)的區(qū)別進(jìn)行修正,從而使危害指數(shù)更接近實(shí)際情況。

在基于組分信息研究聯(lián)合毒性時(shí),組分濃度、作用方式和作用強(qiáng)度等信息對(duì)開(kāi)展評(píng)價(jià)都很關(guān)鍵。在離體毒性測(cè)試中,測(cè)試生物是細(xì)菌、真菌和細(xì)胞,因此一般默認(rèn)暴露濃度即為靶位點(diǎn)濃度。但對(duì)于比較高等的生物體,化合物進(jìn)入體內(nèi)后要經(jīng)歷復(fù)雜的ADME過(guò)程,直接用暴露濃度計(jì)算聯(lián)合效應(yīng)會(huì)造成較大誤差,很不利于毒性評(píng)價(jià)和離體試驗(yàn)結(jié)果的外推。由于靶位點(diǎn)濃度很難獲得,需要用暴露濃度或給藥量進(jìn)行模型推導(dǎo)。PBPK模型就是建立在生理學(xué)和解剖學(xué)上的一種整體模型,研究化合物在體內(nèi)如何分配,屬于毒物代謝動(dòng)力學(xué)范疇[84]。與經(jīng)典的房室模型不同,在PBPK模型中,生物體被看作一系列與組織和器官對(duì)應(yīng)的隔室,以血液、淋巴液、膽汁等為傳輸工具連接[85]。每個(gè)隔室里的化合物濃度都被看作時(shí)間和暴露量的函數(shù),從而可以計(jì)算[86]。它提供了化合物在血液和各個(gè)隔室里的濃度-時(shí)間關(guān)系,反映了物理參數(shù)的變化如何影響ADME過(guò)程,也表現(xiàn)出高濃度到低濃度、急性毒性到慢性毒性的可外推性[87],在化合物風(fēng)險(xiǎn)評(píng)估中有越來(lái)越多的應(yīng)用[88]。有研究報(bào)道在高濃度混合物引起中央神經(jīng)系統(tǒng)毒性時(shí),經(jīng)PBPK修正的危害指數(shù)比基于簡(jiǎn)單濃度加和的危害指數(shù)更大[89]。這說(shuō)明了在混合物風(fēng)險(xiǎn)評(píng)價(jià)中考慮多組分的藥代動(dòng)力學(xué)相互作用的必要性。PBPK模型是一個(gè)分室模型,其對(duì)系統(tǒng)預(yù)測(cè)的置信度水平與分室的復(fù)雜程度有關(guān)。例如,即使分室到細(xì)胞水平,仍然沒(méi)有解決預(yù)測(cè)不同污染物進(jìn)入細(xì)胞后與不同生物大分子之間的相互作用問(wèn)題。

WOE是USEPA推薦的相互作用存下HI的修正方法[6]。它是一種基于經(jīng)驗(yàn)觀察和機(jī)理研究的定性判斷,可以在給定暴露條件下對(duì)相互作用進(jìn)行分類和整理[40,90]。WOE把多組分化合物拆解成二元組合,通過(guò)判斷相互作用方向、毒理學(xué)作用機(jī)制、毒理學(xué)意義的重要性和其他修正因子將每個(gè)組合的相互作用信息轉(zhuǎn)化為用數(shù)字表示的權(quán)重,然后對(duì)HI進(jìn)行修正。在一篇有關(guān)WOE應(yīng)用的綜述文章里, 380個(gè)來(lái)自ATSDR文檔的二元混合物被按照上述3個(gè)方面分組并賦予相應(yīng)權(quán)重[91]。例如,鋅和鉛的聯(lián)合作用表現(xiàn)為拮抗,其作用機(jī)制清晰,毒理學(xué)重要性強(qiáng),此合物被歸入＜IA組。另一個(gè)例子是TCDD和p, p’-DDE的相互作用,p,p’-DDE的存在對(duì)TCDD并無(wú)影響,對(duì)其作用機(jī)制的了解不足,毒理學(xué)重要性較弱,因此屬于=IIIC組。然而,WOE的應(yīng)用經(jīng)常顯示出我們對(duì)二元混合物機(jī)理的了解遠(yuǎn)遠(yuǎn)不足。未來(lái),實(shí)驗(yàn)數(shù)據(jù)、專家判斷和計(jì)算機(jī)工具的結(jié)合將為WOE在混合物風(fēng)險(xiǎn)評(píng)價(jià)里的應(yīng)用提供更好的支持。WOE方法是一個(gè)管理學(xué)方法,并不需要科學(xué)理論的支持。

4.3 動(dòng)態(tài)能量平衡理論

如前所述,化合物進(jìn)入體內(nèi)會(huì)與多重位點(diǎn)發(fā)生作用,導(dǎo)致多種作用方式的影響。致死、生長(zhǎng)抑制、干擾繁殖、致畸致癌,每一類效應(yīng)都包含諸多測(cè)試終點(diǎn)可供進(jìn)行生物測(cè)試。因此,測(cè)試終點(diǎn)的選擇始終是混合物毒性評(píng)價(jià)的重要問(wèn)題,諸多測(cè)試終點(diǎn)的結(jié)果無(wú)法統(tǒng)一、難以外推也一直是研究者的思考重點(diǎn)。在此情況下,動(dòng)態(tài)能量平衡(dynamic energy budget, DEB)理論的提出在一定程度上指出了聯(lián)合毒性研究的新方向。DEB是使用能量作為效應(yīng)端點(diǎn)的研究方法。生物體從外界物質(zhì)中獲得能量,之后這些能量被用于生長(zhǎng)、成熟、繁殖等一系列生物過(guò)程。特定的生物對(duì)能量如何獲得和分配會(huì)保持一個(gè)特定的模式,但當(dāng)暴露在有毒化合物中時(shí),此分配過(guò)程就會(huì)受到干擾、甚至被破壞[92]。DEB把生命周期里的測(cè)試終點(diǎn)都綜合到能量上,以解釋為何混合物毒性能夠隨時(shí)間變化,并在不同測(cè)試終點(diǎn)有不同表現(xiàn)[93]。致死、半致死、生長(zhǎng)繁殖等任何作用終點(diǎn)都與能量流相關(guān),而任何毒物單獨(dú)或聯(lián)合作用都會(huì)擾亂能量的分配過(guò)程。因此,通過(guò)結(jié)合與時(shí)間相關(guān)的實(shí)驗(yàn)結(jié)果、各種相互作用參數(shù)和不同測(cè)試終點(diǎn)的效應(yīng),我們就能解釋混合物組分的特定相互作用。在這個(gè)模型中,能量化身為國(guó)際通用貨幣,影響被認(rèn)為是一個(gè)消費(fèi)過(guò)程。

前文中2個(gè)無(wú)相互作用模型——濃度加和和效應(yīng)加和都是數(shù)學(xué)模型,均將生物體看作黑箱,不關(guān)注作用方式和作用位點(diǎn),針對(duì)不同的測(cè)試終點(diǎn)結(jié)果通常被分開(kāi)解釋；PBPK和WOE具有很強(qiáng)的物種特異性,無(wú)法克服不同種群和個(gè)體間物理、化學(xué)、生化和生理參數(shù)的差異,導(dǎo)致外推困難。而作為毒物代謝動(dòng)力學(xué)和毒物作用動(dòng)力學(xué)的結(jié)合,DEB理論包含了各種與新陳代謝機(jī)制有關(guān)的參數(shù),這樣的綜合分析方法在分析毒理學(xué)內(nèi)在機(jī)理上是不可缺少的。在一項(xiàng)熒蒽和芘對(duì)大型溞不同測(cè)試終點(diǎn)的聯(lián)合毒性研究中利用DEB理論描述了4個(gè)維度(熒蒽濃度,芘濃度,時(shí)間,響應(yīng))和3種測(cè)試終點(diǎn)(存活率,體長(zhǎng),繁殖率)的信息,最終證明2種多環(huán)芳烴作用于相同靶位點(diǎn),其混合物的聯(lián)合毒性表現(xiàn)為加和效應(yīng)[94]。這是目前可以預(yù)測(cè)未經(jīng)測(cè)試的化合物和暴露情境下聯(lián)合毒性的有限選擇。

5 展望(Perspective)

聯(lián)合毒性研究是一個(gè)極具挑戰(zhàn)性但又有很強(qiáng)的現(xiàn)實(shí)意義的課題。從20世紀(jì)上半葉至今,科學(xué)界在無(wú)相互作用模型的應(yīng)用上收獲了許多有價(jià)值的數(shù)據(jù)；而近年來(lái)隨著對(duì)生物體作用位點(diǎn)和作用途徑的研究,又發(fā)展了一系列基于生物學(xué)和生理學(xué)的模型來(lái)將復(fù)雜問(wèn)題簡(jiǎn)單化。這些研究可以看作是一個(gè)良好的開(kāi)端,但是尚沒(méi)有看到能夠解決實(shí)際問(wèn)題的希望。

混合物研究的難點(diǎn)在于當(dāng)環(huán)境復(fù)合污染體系的組分眾多,各組分濃度變化幅度巨大、且作用機(jī)理各不相同。同時(shí),由于單一污染物可以作用于不同靶點(diǎn)、或通過(guò)靶點(diǎn)間是生物學(xué)交叉作用(cross-talk)來(lái)影響最終的效應(yīng)。發(fā)展預(yù)測(cè)能力和外推能力更好的新方法和新模型是未來(lái)的研究方向。如前所述,目前的研究能夠證明在單一靶點(diǎn)或單一作用方式上的聯(lián)合毒性基本符合濃度加和原理,并已經(jīng)廣泛應(yīng)用于諸如二噁英類污染物的環(huán)境管理。因此未來(lái)的研究工作不僅需要研究復(fù)雜體系中各組分的相互作用,更需要了解這種相互作用的毒理學(xué)機(jī)制,認(rèn)識(shí)混合體系中單個(gè)組分的作用方式,及其發(fā)揮聯(lián)合作用的過(guò)程機(jī)理。已經(jīng)有很多研究者利用定量活性結(jié)構(gòu)關(guān)系(quantitative structure-activity relationship,QSAR)模型預(yù)測(cè)未經(jīng)過(guò)毒性測(cè)試的新型化合物的毒性,以此評(píng)價(jià)其在混合物中的聯(lián)合效應(yīng)類型[95],這也成為聯(lián)合毒性研究極有潛力的分支之一?？梢灶A(yù)期,環(huán)境組學(xué)、計(jì)算毒理學(xué)方法與靶向毒性測(cè)試技術(shù)的聯(lián)合應(yīng)用、多學(xué)科知識(shí)的融合才可能會(huì)導(dǎo)致新的突破。

[1]Escher B I,van Daele C,Dutt M,et al.Most oxidative stress response in water samples comes from unknown chemicals:The need for effect-based water quality trigger values[J].Environmental Science&Technology,2013,47 (13):7002-7011

[2]Yang R S H.Introduction to the toxicology of chemical mixtures[M]//Toxicology of Chemical Mixtures.New York:Academic Press,1994:1-10

[3]Agency for Toxic Substances and Disease Registry(ATSDR).Guidance manual for the assessment of joint toxic action of chemical mixtures[R].Atlanta,GA:ASTDR and USDHHS,2004

[4]Organization for Economic Cooperation and Development (OECD).Harmonised integrated classification system for human health and environmental hazards of chemical substances and mixtures[R].Paris:OECD,2001

[5]United States Environmental Protection Agency(USEPA).Guidelines for the health risk assessment of chemical mixtures[R].Washington DC:Risk Assessment Forum, 1986

[6]United States Environmental Protection Agency(USEPA). Supplementary guidance for conducting health risk assessment of chemical mixtures[R].Washington DC:Risk Assessment Forum,2000

[7]World Health Organization(WHO).Assessment of combined exposures to multiple chemicals:Report of a WHO/ IPCS international workshop on aggregate/cumulative risk assessment[R].Washington DC:WHO,2009

[8]Bliss C I.The toxicity of poisons applied jointly[J].Annals of Applied Biology,1939,26(3):585-615

[9]Hewlett P S,Plackett R L.A unified theory for quantal responses to mixtures of drugs:Non-interactive action[J]. Biometrics,1959,15(4):591-610

[10]Plackett R L,Hewlett P S.A comparison of two approaches to the construction of models for quantal responses to mixtures of drugs[J].Biometrics,1967,23(1):27-44

[11]Ashford J R.General models for the joint action of mixtures of drugs[J].Biometrics,1981,37(3):457-474

[12]World Health Organization(WHO).Health effects of combined exposures in the work environment[R].Geneva:WHO,1981

[13]Teuschler L.Deciding which chemical mixtures risk assessment methods work best for what mixtures[J].Toxicology and Applied Pharmacology,2007,223(2):139-147

[14]Loewe S.Die quantitativen probleme der pharmarkologic [J].Ergebnisse der Physiologie,1928,27(1):47-187

[15]Groten J P.Mixtures and interactions[J].Food and Chemical Toxicology,2000,38(S1):S65-S71

[16]P?ch G.Combined Effects of Drugs and Toxic Agents: Modern Evaluation in Theory and Practice[M].Beilin: Springer,1993

[17]Berenbaum M C.What is synergy?[J].Pharmacological Reviews,1989,41(2):93-141

[18]Zhang Y H,Liu S S,Song X Q,et al.Prediction for the mixture toxicity of six organophosphorus pesticides to the luminescent bacterium Q67[J].Ecotoxicology and Environmental Safety,2008,71(3):880-888

[19]Ven Leeuwen C J.Ecotoxicological Effects,in Risk Assessment of Chemicals:An Introduction[M].Dordrecht: Kluwer Academic Publishers,1995:175–237

[20]Berenbaum M C.The expected effect of a combination of agents:The general solution[J].Journal of Theoretical Biology,1985,114(3):413-431

[21]Teuschler L K,Hertzberg R C.Current and future risk assessment guidelines,policy,and methods development for chemical mixtures[J].Toxicology,1995,105(2– 3): 137-144

[22]Evans R M,Scholze M,Kortenkamp A.Examining the feasibility of mixture risk assessment:A case study using a tiered approach with data of 67 pesticides from the Joint FAO/WHO Meeting on Pesticide Residues(JMPR)[J]. Food and Chemical Toxicology,2015,84:260-269

[23]Sprague J B,Ramsay B A.Lethal levels of mixed copper–zinc solutions for juvenile salmon[J].Journal of the Fisheries Research Board of Canada,1965,22(2):425-432

[24]Marking L L,Dawson V K.Method for assessment of toxicity or efficacy of mixtures of chemicals[R].Washington DC:US Fish and Wildlife Service,1975

[25]Van den Berg M,Birnbaum L,Bosveld A T C,et al.Toxic equivalency factors(TEFs)for PCBs,PCDDs,PCDFs for humans and wildlife[J].Environmental Health Perspectives,1998,106(12):775-792

[26]Van Loon W M G M,Verwoerd M E,Wijnker F G,et al. Estimating total body residues and baseline toxicity of complex organic mixtures in effluents and surface waters [J].Environmental Toxicology and Chemistry,1997,16 (7):1358-1365

[27]Fulladosa E,Murat J C,Villaescusa I.Study on the toxicity of binary equitoxic mixtures of metals using the luminescent bacteria as a biological target[J].Chemosphere, 2005,58(5):551-557

[28]Payne J,Rajapakse N,Wilkins M,et al.Prediction and assessment of the effects of mixtures of four xenoestrogens [J].Environmental Health Perspectives,2000,108(10): 983-987

[29]Silva E,Rajapakse N,Kortenkamp A.Something from" nothing"-eight weak estrogenic chemicals combined at concentrations below NOECs produce significant mixture effects[J].Environmental Science&Technology,2002, 36(8):1751-1756

[30]Arcaro K F,Vakharia D D,Yang Y,et al.Lack of synergy by mixtures of weakly estrogenic hydroxylated polychlorinated biphenyls and pesticides[J].Environmental Health Perspectives,1998,106(S4):1041-1046

[31]Payne J,Scholze M,Kortenkamp A.Mixtures of four organochlorines enhance human breast cancer cell proliferation[J].Environmental Health Perspectives,2001,109(4): 391-397

[32]Brian J V,Harris C A,Scholze M,et al.Accurate prediction of the response of freshwater fsh to a mixture of estrogenic chemicals[J].Environmental Health Perspectives,2005,113(6):721-728

[33]Sun L W,Zha J M,Wang Z J.Interactions between estrogenic chemicals in binary mixtures investigated using vitellogenin induction and factorial analysis[J].Chemosphere,2009,75(3):410-415

[34]Ribeiro C,Urbatzka R,Castro L F C,et al.In vitro exposure of Nile tilapia(Oreochromis niloticus)testis to estrogenic endocrine disrupting chemicals:mRNA expression of genes encoding steroidogenic enzymes[J].Toxicology Mechanisms and Methods,2012,22(1):47-53

[35]Ramamoorthy K,Vyhlidal C,Wang F,et al.Additive estrogenic activities of a binary mixture of 2',4',6'-trichloro-and 2',3',4',5'-tetrachloro-4-biphenylol[J].Toxicology and Applied Pharmacology,1997,147(1):93-100

[36]Pohjanvirta R,Tuomisto J.The toxic equivalency principle and its application in dioxin risk assessment[M]//The AH Receptor in Biology and Toxicology.Hoboken:John Wiley&Sons,Inc,2011:317-330

[37]Jaro?ová B,Bláha L,Giesy J P,et al.What level of estrogenic activity determined by in vitro assays in municipal waste waters can be considered as safe?[J].Environment International,2014,64:98-109

[38]Vuc?inic'S,Antonijevic'B,Brkic'D.Occupational and environmental aspects of organophosphorus compounds [M]//Basic and Clinical Toxicology of Organophosphorus Compounds.London:Springer,2014:213-244

[39]Jung K H,Yan B,Chillrud S N,et al.Assessment of benzo(a)pyrene-equivalent carcinogenicity and mutagenicity of residential indoor versus outdoor polycyclic aromatic hydrocarbons exposing young children in New York City[J].International Journal of Environmental Research and Public Health,2010,7(5):1889-1900

[40]Mumtaz M M,DeRosa C T.Approaches and challenges in risk assessment of chemical mixtures[M]//Toxicology of Chemical Mixtures.San Diego:Academic Press,1994: 565-598

[41]Faust M,Altenburger R,Backhaus T,et al.Joint algal toxicity of 16 dissimilarly acting chemicals is predictable by the concept of independent action[J].Aquatic Toxicology,2003,63(1):43-63

[42]Backhaus T,Arrhenius A,Blanck H.Toxicity of a mixture of dissimilarly acting substances to natural algal communities:Predictive power and limitations of independent action and concentration addition[J].Environmental Science&Technology,2004,38(23):6363-6370

[43]Junghans M.Studies on combination effects of environmentally relevant toxicants[D].Bremen:University of Bremen,2004:78-106

[44]Ra J S,Lee B C,Chang N I,et al.Estimating the combined toxicity by two-step prediction model on the complicated chemical mixtures from wastewater treatment plant effluents[J].Environmental Toxicology and Chemistry,2006,25(8):2107-2113

[45]Altenburger R,Walter H,Grote M.What contributes to the combined effect of a complex mixture?[J].Environmental Science&Technology,2004,38(23):6353-6362

[46]Chen C Y,Lu C L.An analysis of the combined effects of organic toxicants[J].Science of the Total Environment, 2002,289(1–3):123-132

[47]Fulladosa E,Murat J-C,Villaescusa I.Study on the toxicity of binary equitoxic mixtures of metals using the luminescent bacteriaVibrio fischerias a biological target[J]. Chemosphere,2005,58(5):551-557

[48]Yang R,Li N,Ma M,et al.Combined effects of estrogenic chemicals with the same mode of action using an estrogen receptor binding bioassay[J].Environmental Toxicology and Pharmacology,2014,38(3):829-837

[49]Ermler S,Scholze M,Kortenkamp A.The suitability of concentration addition for predicting the effects of multicomponent mixtures of up to 17 anti-androgens with varied structural features in an in vitro AR antagonist assay [J].Toxicology and Applied Pharmacology,2011,257(2): 189-197

[50]Olmstead A W,LeBlanc G A.Joint action of polycyclic aromatic hydrocarbons:Predictive modeling of sublethal toxicity[J].Aquatic Toxicology,2005,75(3):253-262

[51]Hodges G,Roberts D W,Marshall S J,et al.Defining the toxic mode of action of ester sulphonates using the joint toxicity of mixtures[J].Chemosphere,2006,64(1):17-25

[52]Baylay A J,Spurgeon D J,Svendsen C,et al.A metabolomics based test of independent action and concentration addition using the earthwormLumbricus rubellus[J].Ecotoxicology,2012,21(5):1436-1447

[53]Suzuki T,Ide K,Ishida M.Response of MCF-7 human breast cancer cells to some binary mixtures of oestrogenic compounds in-vitro[J].Journal of Pharmacy and Pharmacology,2001,53(11):1549-1554

[54]Rajapakse N,Silva E,Scholze M,et al.Deviation from additivity with estrogenic mixtures containing 4-nonylphenol and 4-tert-octylphenol detected in the E-SCREEN assay[J].Environmental Science&Technology,2004,38 (23):6343-6352

[55]Silva E,Rajapakse N,Scholze M,et al.Joint effects of heterogeneous estrogenic chemicals in the E-Screen-exploring the applicability of concentration addition[J]. Toxicological Sciences,2011,122(2):383-394

[56]Arnold S F,Bergeron J M,Tran D Q,et al.Synergistic responses of steroidal estrogens in vitro(yeast)and in vivo(turtles)[J].Biochemical and Biophysical Research Communications,1997,235(2):336-342

[57]Bergeron J M,Crews D,McLachlan J A.PCBs as environmental estrogens:Turtle sex determination as a biomarker of environmental contamination[J].Environmental Health Perspectives,1994,102(9):780-781

[58]Tinwell H,Ashby J.Sensitivity of the immature rat uterotrophic assay to mixtures of estrogens[J].Environmental Health Perspectives,2004,112(5):575-582

[59]Vijver M G,Peijnenburg W J G M,Snoo G R D.Toxicological mixture models are based on inadequate assumptions[J].Environmental Science&Technology,2010(44): 4841-4842

[60]Backhaus T,Faust M.Predictive environmental risk assessment of chemical mixtures:A conceptual framework [J].Environmental Science&Technology,2012,46(5): 2564-2573

[61]Cedergreen N,Christensen A M,Kamper A,et al.A review of independent action compared to concentration addition as reference models for mixtures of compounds with different molecular target sites[J].Environmental Toxicology and Chemistry,2008,27(7):1621-1632

[62]Spurgeon D J,Jones O A H,Dorne J-L C M,et al.Systems toxicology approaches for understanding the joint effects of environmental chemical mixtures[J].Science of the Total Environment,2010,408(18):3725-3734

[63]Cassee F R,S u¨hnel J,Groten J P,et al.Toxicology of chemical mixtures[M]//General and Applied Toxicology. London:Macmillan Reference Ltd,1999:303-319

[64]Groten J P,Tajima O,Feron V J,et al.Toxicology of chemical mixtures:International perspective[J].Environmental Health Perspectives,1998,106(S6):1281-1289

[65]Finney D J.Probit Analysis[M].London:Cambridge U-niversity Press,1952

[66]Keplinger M L,Deichmann W B.Acute toxicity of combinations of pesticides[J].Toxicology and Applied Pharmacology,1967,10(3):586-595

[67]Smyth Jr H F,Weil C S,West J S,et al.An exploration of joint toxic action:Twenty-seven industrial chemicals intubated in rats in all possible pairs[J].Toxicology and Applied Pharmacology,1969,14(2):340-347

[68]Li X F,Zhou Q X,Luo Y,et al.Joint action and lethal levels of toluene,ethylbenzene,and xylene on midge (Chironomus plumosus)larvae[J].Environmental Science and Pollution Research,2013,20(2):957-966

[69]Mori I C,Arias-Barreiro C R,Koutsaftis A,et al.Toxicity of tetramethylammonium hydroxide to aquatic organisms and its synergistic action with potassium iodide[J]. Chemosphere,2015,120:299-304

[70]曾鳴,林志芬,尹大強(qiáng),等.混合污染物聯(lián)合毒性研究進(jìn)展[J].環(huán)境科學(xué)與技術(shù),2009(2):80-86 Zeng M,Lin Z F,Yin D Q,et al.Progress on joint effect of mixture pollutants[J].Environmental Science&Technology,2009(2):80-86(in Chinese)

[71]Altenburger R,Nendza M,Schüürmann G.Mixture toxicity and its modeling by quantitative structure-activity relationships[J].Environmental Toxicology and Chemistry, 2003,22(8):1900-1915

[72]Groten J P,Tajima O,Feron V J,et al.Statistically designed experiments to screen chemical mixtures for possible interactions[J].Environmental Health Perspectives, 1998,106(S6):1361-1365

[73]Greco W R,Bravo G,Parsons J C.The search for synergy:A critical review from a response surface perspective [J].Pharmacological Reviews,1995,47(2):331-385

[74]Charles G D,Gennings C,Zacharewski T R,et al.An approach for assessing estrogen receptor-mediated interactions in mixtures of three chemicals:A pilot study[J]. Toxicological Sciences,2002,68(2):349-360

[75]李小朋,張琛,李娟,等.基于析因設(shè)計(jì)的多種重金屬對(duì)發(fā)光菌聯(lián)合毒性研究[J].環(huán)境科學(xué)與技術(shù),2012,35 (12):169-174 Li X P,Zhang C,Li J,et al.Joint toxicity of multi-heavymetal to photobacterium phosphoreum based on factorial design[J].Environmental Science&Technology,2012,35 (12):169-174(in Chinese)

[76]Rijstenbil J W,Gerringa L J A.Interactions of algal ligands,metal complexation and availability,and cell responses of the diatomDitylum brightwelliiwith a gradual increase in copper[J].Aquatic Toxicology,2002,56(2): 115-131

[77]Jin X W,Wang Y Y,Jin W,et al.Ecological risk of nonylphenol in China surface waters based on reproductive fitness[J].Environmental Science&Technology,2014,48 (2):1256-1262

[78]González-Pleiter M,Gonzalo S,Rodea-Palomares I,et al. Toxicity of five antibiotics and their mixtures towards photosynthetic aquatic organisms:Implications for environmental risk assessment[J].Water Research,2013,47 (6):2050-2064

[79]Altenburger R,Scholz S,Schmitt-Jansen M,et al.Mixture toxicity revisited from a toxicogenomic perspective[J]. Environmental Science&Technology,2012,46(5):2508-2522

[80]Yang R,Li N,Rao K F,et al.Combined action of estrogen receptor agonists and antagonists in two-hybrid recombinant yeast in vitro[J].Ecotoxicology and Environmental Safety,2015,111(1):228-235

[81]章強(qiáng),辛琦,強(qiáng)麗媛,等銅及其與四環(huán)素的聯(lián)合暴露對(duì)斑馬魚胚胎的毒性效應(yīng)[J].生態(tài)毒理學(xué)報(bào),2015,10(5): 35-46 Zhang Q,Xin Q,Qiang L Y,et al.Developmental toxicity of Cu and tetracycline on zebrafish embryos(Danio Rerio)[J].Asian Journal of Ecotoxicology,2015,10(5):35-46(in Chinese)

[82]吳寅,吳永貴,馬嵐,等.水體酸化條件下Cu(II)對(duì)斑馬魚胚胎的毒性效應(yīng)[J].生態(tài)毒理學(xué)報(bào),2013,8(3):389-394 Wu Y,Wu Y G,Ma L,et al.Toxic effect of Cu(II)on zebrafish embryos under water acidification condition[J]. Asian Journal of Ecotoxicology,2013,8(3):389-394(in Chinese)

[83]United States Department of Health and Human Services (USDHHS).Guidance manual for the assessment of joint toxic action of chemical mixtures(DRAFT)[R].Washington DC:USDHHS,2001

[84]Teorell T.Kinetics of distribution of substances administered to the body.I:The extravascular modes of administration[J].Arch Int Pharmacodyn,1937,57:205-225

[85]Leahy D E.Progress in simulation modelling for pharmacokinetics[J].Current Topics in Medicinal Chemistry, 2003,3(11):1257-1268

[86]Parrott N,Jones H,Paquereau N,et al.Application of full physiological models for pharmaceutical drug candidate selection and extrapolation of pharmacokinetics to man [J].Basic and Clinical Pharmacology and Toxicology, 2005,96(3):193-199

[87]Roy A.Mechanistic modeling of transport and metabolism in physiological systems[D].Washington DC:Rutgers University,1997:1-7

[88]Mumtaz M,Fisher J,Blount B,et al.Application of physiologically based pharmacokinetic models in chemical risk assessment[J].Journal of Toxicology,2012,DOI:10. 1155/2012/904603.

[89]Haddad S,Béliveau M,Tardif R,et al.A PBPK modeling-based approach to account for interactions in the health risk assessment of chemical mixtures[J].Toxicological Sciences,2001,63(1):125-131

[90]Hertzberg R C,Rice G,Teuchler L K.Methods for health risk assessment of combustion mixtures[M]//Hazardous Waste Incineration:Evaluating the Human Health and Environmental Risks.Boca Raton,FL:Lewis Publishers, 1999:105-148

[91]Pohl H R,Mumtaz M M,Scinicariello F,et al.Binary weight-of-evidence evaluations of chemical interactions—15 years of experience[J].Regulatory Toxicology and Pharmacology,2009,54(3):264-271

[92]Nisbet R M,Muller E B,Lika K,et al.From molecules to ecosystems through dynamic energy budget models[J]. Journal of Animal Ecology,2000,69(6):913-926

[93]Jager T,Vandenbrouck T,Baas J,et al.A biology-based approach for mixture toxicity of multiple endpoints over the life cycle[J].Ecotoxicology(London,England),2010, 19(2):351-361

[94]Baas J,Jager T,Kooijman B.A review of DEB theory in assessing toxic effects of mixtures[J].Science of the Total Environment,2010,408(18):3740-3745

[95]Kim J,Kim S.State of the art in the application of QSAR techniques for predicting mixture toxicity in environmental risk assessment[J].SAR and QSAR in Environmental Research,2015,26(1):41-59

Review on Methodology for Environmental Mixture Toxicity Study

Yang Rong1,2,Li Na2,Rao Kaifeng2,Ma Mei2,*,Wang Zijian2
1.Beijing Water Quality Monitoring Center for South-to-North Water Diversion,Beijing 100195,China
2.Key Laboratory of Drinking Water Science and Technology,Research Center for Eco-Environmental Sciences,Chinese Academy of Sciences,Beijing 100085,China

25 November 2015 accepted 8 January 2016

Both human and wild animals are simultaneously exposed to numbers of environmental pollutants,and via multiple exposure pathways.The adverse effect of single chemical could sometimes be ignored due to its low environmental concentration.However the effects of co-existing multiple chemicals could not be negligible because of the so-called“something from nothing”phenomenon.Therefore the joint toxicity of a mixture arouse to be a hot topic,especially on the additive and potential synergistic actions of environmental mixtures.A numbers of mathematical or biological models have been developed and applied to investigate the combined effects of environmental pollutants.In this review paper,widely used and newly emerging models were evaluated on their suitability to the circumstances,and the future challenges in toxicology of chemical mixture research and their applications in environmental management were discussed.

environmental mixture;combined effects;mode of action;model prediction

2015-11-25 錄用日期:2016-01-08

1673-5897(2016)1-001-13

X171.5

A

10.7524/AJE.1673-5897.20151125002

楊蓉,李娜,饒凱鋒,等.環(huán)境混合物的聯(lián)合毒性研究方法[J].生態(tài)毒理學(xué)報(bào),2016,11(1):1-13

Yang R,Li N,Rao K F,et al.Review on methodology for environmental mixture toxicity study[J].Asian Journal of Ecotoxicology,2016,11(1):1-13(in Chinese)

自然科學(xué)基金重點(diǎn)課題(No.51290283)；自然科學(xué)基金青年基金(No.21307131)；國(guó)家高技術(shù)研究發(fā)展計(jì)劃項(xiàng)目(No.2014AA06A506)

楊蓉(1987-),女,博士,研究方向?yàn)樗鷳B(tài)毒理學(xué),E-mail:ygrg1987@gmail.com；

),E-mail:mamei@rcees.ac.cn

簡(jiǎn)介:馬梅(1967—)，女，博士，研究員，博士生導(dǎo)師，主要從事水生態(tài)毒理學(xué)研究。