孫健，肖鵬飛，劉毅華，尹曉輝，朱國(guó)念

1. 浙江農(nóng)林大學(xué)，杭州 311300 2. 浙江農(nóng)林大學(xué)暨陽(yáng)學(xué)院，諸暨 311800 3. 中國(guó)林科院亞熱帶林業(yè)研究所，富陽(yáng) 311400 4. 浙江大學(xué)農(nóng)藥與環(huán)境毒理研究所，杭州 310029

三唑酮(triadimefon)是一種甾醇合成抑制劑，因其廣譜性的殺菌作用，可用于防治水稻銹病[1]。在我國(guó)，三唑酮在農(nóng)田中被廣泛使用。施用方式的原因?qū)е氯蛲寐什桓?，最終在環(huán)境中存在殘留。根據(jù)中華人民共和國(guó)環(huán)境保護(hù)部2008年發(fā)布的《雜環(huán)類(lèi)農(nóng)藥工業(yè)水污染物排放標(biāo)準(zhǔn)》，三唑酮原藥生產(chǎn)企業(yè)廢水處理設(shè)施總排放口的排放標(biāo)準(zhǔn)為5 000 μg·L-1。如果按照國(guó)家標(biāo)準(zhǔn)的臨界值濃度排放含有三唑酮的工業(yè)廢水，三唑酮有可能在環(huán)境中表現(xiàn)出短期高濃度暴露[2]。付巖[3]在農(nóng)藥使用期間檢測(cè)到稻田溝渠水中三唑酮的殘留濃度為12 μg·L-1，稻田附近池塘檢測(cè)到三唑酮濃度達(dá)到0.4 μg·L-1。三唑酮的用途廣泛，在國(guó)外，三唑酮可用于草坪防治灰斑病，Vincelli[4]模擬了高爾夫球場(chǎng)草地施用三唑酮后，檢測(cè)周?chē)乇硭修r(nóng)藥的殘留情況，結(jié)果表明，水中三唑酮含量峰值可高達(dá)300 μg·L-1。除了三唑酮在施用地點(diǎn)的環(huán)境濃度有報(bào)道外，在江河湖口也有三唑酮檢出。游明華[5]在九龍江江口豐水期時(shí)檢測(cè)到三唑酮含量達(dá)0.173 μg·L-1，劉娜等[6]在2014—2017年于太湖流域飲用水源地檢測(cè)出三唑酮，濃度為0.002 μg·L-1。

三唑酮可以通過(guò)抑制細(xì)胞色素P-450酶系，以抑制麥角甾醇合成，使新生真菌細(xì)胞壁畸形，達(dá)到抑菌效果[7]。P-450酶系參與生物體多種內(nèi)分泌生物通路，因此，三唑酮可能會(huì)對(duì)水生動(dòng)物造成影響，這個(gè)影響可能是由于甾醇合成受抑制而導(dǎo)致生物內(nèi)分泌紊亂[8]，進(jìn)而對(duì)脊椎動(dòng)物和無(wú)脊椎動(dòng)物產(chǎn)生致畸作用。劉少穎[9]發(fā)現(xiàn)在實(shí)驗(yàn)室條件下三唑酮濃度達(dá)到2 000 μg·L-1和4 000 μg·L-1時(shí)，斑馬魚(yú)胚胎出現(xiàn)畸形。三唑酮進(jìn)入水體后可在非洲爪蟾體內(nèi)累積，干預(yù)神經(jīng)嵴細(xì)胞的遷移，致使爪蟾胚胎鰓弓畸變[10]。Kenneke等[11]的21 d慢性毒性研究表明，當(dāng)三唑酮濃度>50 μg·L-1時(shí)可導(dǎo)致大型溞幼溞蛻皮頻率降低。有研究表明，三唑酮可以推遲大型溞的蛻皮時(shí)間，并對(duì)子代眼部造成畸形[12]。三唑酮可以抑制生物發(fā)育，進(jìn)而導(dǎo)致生物種群增長(zhǎng)受抑制。胡芳華等[13]的研究表明，三唑酮的濃度達(dá)到160 μg·L-1時(shí)，三唑酮可以影響母溞的產(chǎn)溞數(shù)量，并抑制幼溞的體長(zhǎng)。

水環(huán)境中往往含有多種農(nóng)藥，三唑酮可以和其他殺蟲(chóng)劑產(chǎn)生協(xié)同增效毒性作用。Cedergreen和Streibig[14]的研究表明，很多殺蟲(chóng)劑和殺菌劑的組合存在明顯的協(xié)同增效作用，其中，大約96%屬于唑類(lèi)殺菌劑(包括咪唑、吡唑和三唑類(lèi))和神經(jīng)毒劑殺蟲(chóng)劑的組合。Bjergager等[15]通過(guò)聯(lián)合使用咪酰胺和高效氯氰菊酯，發(fā)現(xiàn)該協(xié)同增效作用在浮游動(dòng)物群落水平上同樣存在。淡水浮游動(dòng)物在水生態(tài)系統(tǒng)中起著承上啟下的作用，受到影響后會(huì)影響其他水生動(dòng)物。

因此，研究三唑酮在水環(huán)境中的生態(tài)毒性十分必要。目前，三唑酮的研究主要是其對(duì)動(dòng)物個(gè)體的影響，對(duì)動(dòng)物的群落效應(yīng)研究不多。本文利用標(biāo)準(zhǔn)微宇宙進(jìn)行試驗(yàn)，評(píng)定三唑酮對(duì)淡水浮游動(dòng)物種群的生態(tài)風(fēng)險(xiǎn)。

1 材料與方法(Materials and methods)

1.1 微宇宙試驗(yàn)系統(tǒng)構(gòu)建

室內(nèi)微宇宙系統(tǒng)在投藥前14 d(記為-14 d)開(kāi)始裝配，依次加入4種淡水藻類(lèi)(購(gòu)自中國(guó)科學(xué)院水生生物研究所)和8種浮游動(dòng)物(由浙江大學(xué)農(nóng)藥與生態(tài)毒理研究所提供)，生物種類(lèi)及數(shù)量如表1所示。微宇宙系統(tǒng)靜置孵育4 d后投放農(nóng)藥，投藥當(dāng)日記為0 d。微宇宙系統(tǒng)構(gòu)建工作流程圖如圖1所示。整個(gè)實(shí)驗(yàn)周期，微宇宙系統(tǒng)安置于大型水浴裝置上(圖3(a))。自-14～0 d，為微宇宙系統(tǒng)構(gòu)建平衡期，保持溫度24～26 ℃，光強(qiáng)度3 000 lux，光周期L∶D = 14 h∶10 h。自0 d起，每隔7 d向系統(tǒng)補(bǔ)充藻類(lèi)(1×103cells·mL-1，1 mL)和浮游動(dòng)物(輪蟲(chóng)儲(chǔ)備液1 mL，甲殼綱生物5個(gè)成體)。

圖1 微宇宙系統(tǒng)構(gòu)建工作流程圖注：(a)水生微宇宙系統(tǒng)；(b)大型控溫水浴設(shè)備。Fig. 1 Aquatic microcosm system internal structure and establishment flowsheetNote: (a) Aquatic microcosm system; (b) thermostat water bath equipment.

1.2 受試化學(xué)品、濃度設(shè)置及其檢測(cè)方式

受試化學(xué)品為三唑酮原藥(純度98.5%，由江蘇劍牌股份有限公司提供)。試驗(yàn)設(shè)置6個(gè)濃度，依次為0 μg·L-1(ds0)、156 μg·L-1(ds1)、315 μg·L-1(ds2)、625 μg·L-1(ds3)、1 250 μg·L-1(ds4)和2 500 μg·L-1(ds5)，每個(gè)暴露濃度重復(fù)3次。在投藥后2 h采集所有處理組的水樣，在投藥后1、3、5、7、10、14、21、28和35 d采集最高濃度處理組的水樣，檢測(cè)三唑酮在水中的殘留量。

所有水樣通過(guò)有機(jī)溶劑萃取，采用氣相色譜法檢測(cè)其濃度，色譜條件為：Agilent 6 890氣相色譜儀，HP-5色譜柱，ECD檢測(cè)器，進(jìn)樣口溫度240 ℃，檢測(cè)器溫度250 ℃，柱溫初始溫度70 ℃，進(jìn)樣體積1 μL，載氣為高純氮(99.999%)，恒壓180 kPa，梯度升溫程序?yàn)橹鶞?0 ℃保持2 min，15 ℃·min-1升至190 ℃保持1 min，5 ℃·min-1升至240 ℃。進(jìn)樣模式為不分流進(jìn)樣。在此檢測(cè)條件下，三唑酮的保留時(shí)間為9.4 min(圖2)。三唑酮的回收率為93.6% (SD=0.1)。

圖2 三唑酮檢測(cè)的氣相色譜圖Fig. 2 Gas chromatography of triadimefon

1.3 水質(zhì)理化參數(shù)測(cè)定

水質(zhì)理化指標(biāo)包括酸堿度(pH)、電導(dǎo)率(EC)、溶解氧濃度(DO)和可溶性有機(jī)碳濃度(DOC)。酸堿度和電導(dǎo)率的測(cè)定均在每次生物樣品采集前進(jìn)行。為排除采集水樣對(duì)溶解氧濃度的干擾，溶解氧濃度測(cè)定在水樣采集前一日完成(即試驗(yàn)進(jìn)行的第1、6、13、20、27和34天)，測(cè)定2次，分別在光周期開(kāi)始前(在8:00—9:00內(nèi)完成)，定義為系統(tǒng)中溶解氧濃度最低值(DOd)；暗周期開(kāi)始前(在22:00—23:00內(nèi)完成)，定義為系統(tǒng)中溶解氧濃度最高值(DOL)。計(jì)算得到每日溶氧消耗量(DOδ= DOL-DOd)。溶氧消耗量可以間接表征微宇宙系統(tǒng)整體的代謝情況。可溶性有機(jī)碳樣品采集與生物樣品采集同時(shí)進(jìn)行，采集50 mL水樣經(jīng)過(guò)0.4 μm水性濾膜過(guò)濾，使用CN200TOC-VCPH FA有機(jī)碳測(cè)定儀檢測(cè)水樣中有機(jī)碳濃度。

表1 構(gòu)建標(biāo)準(zhǔn)化微宇宙的浮游動(dòng)物及淡水藻類(lèi)種類(lèi)來(lái)源Table 1 Freshwater algae source and zooplankton for standardized microcosm

1.4 生物樣品采集與檢測(cè)

在0 d投藥前以及投藥后7、14、21、28和35 d采集上覆水，采用0.5 L深水采樣器(圖3(c))采集上覆水樣品。分3次采集水樣，共采集1.5 L上覆水，混勻后保留0.5 L分樣品用于浮游生物鏡檢。向樣品中滴加2 mL魯哥氏液(Lugol’s solution)固定染色，靜置48 h后用虹吸管(管頭包有孔徑為25 μm絹紗)排去上清液，濃縮至20 mL左右，用10%甲醛溶液定容至30 mL。濃縮樣品裝入棕色塑料樣品瓶，4 ℃冷藏待鏡檢。體視鏡下采用5 mL和0.1 mL生物技術(shù)框(圖3(b))鏡檢計(jì)數(shù)浮游動(dòng)物。

圖3 自制的微宇宙系統(tǒng)設(shè)備注：(a)大型控溫水浴設(shè)備；(b)浮游生物計(jì)數(shù)板；(c)深水采樣器。Fig. 3 Microcosmic test equipmentNote: (a) thermostat water bath equipment; (b) plankton counting chamber; (c) water sampler.

1.5 數(shù)據(jù)分析

使用SPSS 22.0對(duì)水質(zhì)參數(shù)、浮游動(dòng)物數(shù)據(jù)進(jìn)行描述性統(tǒng)計(jì)分析。對(duì)浮游動(dòng)物群落結(jié)構(gòu)變化采用Canoco?5.0軟件分析，分析內(nèi)容包括(1)對(duì)微宇宙系統(tǒng)內(nèi)的浮游動(dòng)物組成進(jìn)行主成分分析PCA[16-18]，獲得浮游動(dòng)物在微宇宙內(nèi)的第一排序軸取值(CaseRZ.)；(2)對(duì)0～35 d內(nèi)浮游動(dòng)物豐度的變化進(jìn)行主效應(yīng)曲線分析PRC[19-22]，并進(jìn)行蒙特卡洛置換檢驗(yàn)(Monte Carlo permutation test)，隨機(jī)置換599次，置信限為95%，以驗(yàn)證標(biāo)準(zhǔn)化微宇宙中浮游動(dòng)物隨時(shí)間變化的顯著性。浮游動(dòng)物種群的最大無(wú)作用劑量(NOECspecies)的計(jì)算方法為通過(guò)方差分析得到與對(duì)照組差異不顯著的最大濃度；浮游動(dòng)物群落的最大無(wú)作用劑量(NOECcommunity)的計(jì)算方法為：按照PCA分析得到2個(gè)獨(dú)立的統(tǒng)計(jì)變量CaseR1.和CaseR2.，將三唑酮濃度和取樣時(shí)間設(shè)為自變量，各處理組及不同采樣日樣品的CaseR1.和CaseR2.設(shè)為因變量，取每個(gè)取樣樣品與對(duì)照組在95%置信限水平下無(wú)顯著性差異的最高濃度。

2 結(jié)果(Results)

2.1 三唑酮在水體中的濃度變化

投藥后2 h水中三唑酮濃度的分析結(jié)果顯示，三唑酮實(shí)際暴露濃度為名義濃度的65.3%～86.1%(表2)。最高濃度處理組的三唑酮消解動(dòng)態(tài)結(jié)果如圖4所示。三唑酮的消解動(dòng)力學(xué)方程為cDT= 2705.5e-0.12t，cDT為時(shí)間為t時(shí)的三唑酮濃度，t為時(shí)間，R2= 0.9782，半衰期DT1/2為8.28 d。

表2 三唑酮處理組的理論和實(shí)際濃度Table 2 Nominal and actual concentrations of triadimefon treatment

圖4 三唑酮在標(biāo)準(zhǔn)化微宇宙上覆水中的消解動(dòng)力學(xué)曲線注：cDT為時(shí)間為t時(shí)的三唑酮濃度，t為時(shí)間。Fig. 4 Kinetics curve of triadimefon in standardized aquatic microcosm overlying waterNote: When the time is t, the concentration of triadimefon is cDT.

2.2 三唑酮對(duì)水質(zhì)參數(shù)的影響

在試驗(yàn)周期內(nèi)測(cè)定了所有微宇宙系統(tǒng)的光周期和暗周期的溶解氧濃度，并在根據(jù)溶氧檢測(cè)數(shù)據(jù)計(jì)算微宇宙系統(tǒng)中暗周期溶解氧消耗量(DOc)。三唑酮處理組的DOc變化趨勢(shì)與對(duì)照組一致(P>0.05，Duncan test)。各個(gè)微宇宙系統(tǒng)的pH值始終在7.08～7.65，保持中性。方差分析顯示，各系統(tǒng)酸堿度在處理組和時(shí)間上沒(méi)有顯著性差異(P>0.05，Duncan test)。各個(gè)微宇宙系統(tǒng)的電導(dǎo)率(EC)在694～794 μS·cm-1范圍內(nèi)無(wú)顯著性差異(P>0.05，Duncan test)。各微宇宙系統(tǒng)的可溶性有機(jī)碳濃度(DOC)范圍在39.40～44.09 mg·L-1，方差分析顯示無(wú)顯著性差異(P>0.05，Duncan test)。由此判斷三唑酮在(119.06 ± 11.16)～(2 078.88 ± 153.18) μg·L-1濃度下對(duì)標(biāo)準(zhǔn)微宇宙系統(tǒng)水質(zhì)指標(biāo)沒(méi)有影響。

2.3 三唑酮對(duì)單物種豐度的影響

試驗(yàn)開(kāi)始(0 d)時(shí)，各微宇宙系統(tǒng)間8種浮游動(dòng)物數(shù)量無(wú)顯著性差異，平均豐度分別為：無(wú)節(jié)幼蟲(chóng)(327 ± 126.4) ind·L-1、鋸緣真劍水蚤(107 ± 41.4) ind·L-1、中華薄殼介(5 ± 2.9) ind·L-1、隆線溞(10 ± 7.2) ind·L-1、低額溞(7 ± 3.2) ind·L-1、平突船卵溞(10 ± 6.0) ind·L-1、盤(pán)腸溞(9 ± 5.7) ind·L-1、萼花臂尾輪蟲(chóng)(45 ± 32.2) ind·L-1和單趾輪蟲(chóng)(52 ± 25.5) ind·L-1。隨著暴露時(shí)間的延長(zhǎng)，三唑酮對(duì)枝角類(lèi)生物的影響比較明顯，三唑酮處理組中的枝角類(lèi)生物豐度大多低于同期對(duì)照組(表3)，當(dāng)三唑酮濃度>104.97 μg·L-1(實(shí)際暴露濃度，下同)，在7～21 d時(shí)，處理組枝角類(lèi)生物豐度與對(duì)照組差異顯著。

表3 三唑酮對(duì)枝角類(lèi)豐度的影響Table 3 Effects of triadimefon on Cladocerans secies abundance

2.4 三唑酮對(duì)浮游動(dòng)物種群的影響

對(duì)照組和各三唑酮處理組中浮游動(dòng)物組成變化經(jīng)PRC分析的結(jié)果如圖5所示。置換檢驗(yàn)結(jié)果表明，三唑酮處理組中浮游動(dòng)物組成隨時(shí)間變化達(dá)到顯著水平(Pseudo-F = 12.9，P= 0.002)。圖5中Resp.1反映的是不同物種對(duì)三唑酮的敏感性差異。本研究中浮游生物對(duì)三唑酮的敏感性依次為：盤(pán)腸溞>平突船卵溞>隆線溞>萼花臂尾輪蟲(chóng)>低額溞≥鋸緣真劍水蚤≥無(wú)節(jié)幼蟲(chóng)≥單趾輪蟲(chóng)≥中華薄殼介。三唑酮處理對(duì)單趾輪蟲(chóng)以及介蟲(chóng)沒(méi)有明顯的作用(Resp.1接近0)；鋸緣真劍水蚤與其無(wú)節(jié)幼蟲(chóng)受到三唑酮輕微的影響；三唑酮對(duì)枝角類(lèi)生物具有顯著的負(fù)效應(yīng)(Resp.1為1～4.8)。

圖5 三唑酮對(duì)浮游動(dòng)物組成結(jié)構(gòu)的影響注：圖例中的濃度值為三唑酮理論濃度/實(shí)際濃度。Fig. 5 Effects of triadimefon on zooplankton assemblageNote: The concentration is theoretical concentration/actual concentration in legend.

筆者研究發(fā)現(xiàn)，三唑酮對(duì)室內(nèi)標(biāo)準(zhǔn)微宇宙試驗(yàn)的浮游動(dòng)物群落沒(méi)有出現(xiàn)崩潰式的影響[23]。三唑酮暴露下的浮游動(dòng)物群落演替過(guò)程可分為：演替前期(表現(xiàn)農(nóng)藥直接效應(yīng))、演替后期(表現(xiàn)農(nóng)藥間接效應(yīng))和平衡期[24]。在試驗(yàn)初期(即0～7 d)三唑酮對(duì)枝角類(lèi)生物產(chǎn)生明顯的負(fù)效應(yīng)，對(duì)橈足類(lèi)生物則產(chǎn)生輕微的負(fù)效應(yīng)。三唑酮對(duì)浮游動(dòng)物群落的影響逐漸減小。試驗(yàn)開(kāi)始14 d后，三唑酮對(duì)浮游動(dòng)物群落的影響不明顯。在試驗(yàn)21 d后除了最高濃度處理組，其余的處理組生物豐度有顯著的回升。

三唑酮對(duì)整個(gè)室內(nèi)微宇宙的顯著效應(yīng)主要在試驗(yàn)初期，在該時(shí)間段三唑酮對(duì)隆線溞的NOEC為537.98 μg·L-1(實(shí)際濃度，下同)；對(duì)低額溞、平突船卵溞和盤(pán)腸溞的NOEC均為1 042.97 μg·L-1。其他階段，三唑酮對(duì)其他生物的NOEC均>2 078.88 μg·L-1。經(jīng)過(guò)PCA分析得到統(tǒng)計(jì)變量CaseR1.和CaseR2.。CaseR1.體現(xiàn)了輪蟲(chóng)物種豐度信息，CaseR2.主要體現(xiàn)了三唑酮對(duì)其他甲殼類(lèi)浮游生物的影響。兩者所代表的生物變量組成接近，對(duì)變量總體的解釋能力分布為30.87%和20.02%。各物種對(duì)CaseR1.和CaseR2.排序軸的貢獻(xiàn)度如表4所示。對(duì)CaseR1.和CaseR2.進(jìn)行方差分析，結(jié)果表明，在浮游動(dòng)物水平上三唑酮僅在試驗(yàn)早期7 d內(nèi)有明顯作用，其N(xiāo)OECcommunity為537.98 μg·L-1，對(duì)三唑酮處理組浮游動(dòng)物主成分分析如圖6所示。三唑酮整個(gè)試驗(yàn)周期NOECcommunity>2 078.88 μg·L-1(表5)

圖6 三唑酮處理組中浮游動(dòng)物組成結(jié)構(gòu)主成分分析注：(a)和(b)象限內(nèi)物種種群數(shù)量與三唑酮濃度呈正相關(guān)；(c)和(d)象限內(nèi)的呈負(fù)相關(guān)。Fig. 6 Principal component analysis of triadimefon treated zooplankton communityNote: Species in (a) and (b) quadrants are positive correlation with concentration of triadimefon, and those in (c) and (d) quadrants are negative correlation.

3 討論(Discussion)

三唑酮對(duì)微宇宙系統(tǒng)的直接效應(yīng)主要表現(xiàn)在施藥后的0～7 d。從PRC圖分析可以看出，浮游動(dòng)物對(duì)三唑酮的敏感程度排序主要是枝角類(lèi)>橈足類(lèi)≥輪蟲(chóng)類(lèi)≥介形蟲(chóng)。其中，枝角類(lèi)為最敏感的動(dòng)物，這與Yin等[25]對(duì)咯菌腈暴露的室外微宇宙研究的結(jié)果類(lèi)似。三唑酮對(duì)整個(gè)室內(nèi)微宇宙的顯著效應(yīng)主要在試驗(yàn)初期(即0～7 d)，在該時(shí)間段三唑酮對(duì)隆線溞的NOEC為537.98 μg·L-1(實(shí)際濃度，下同)；對(duì)低額溞、平突船卵溞和盤(pán)腸溞的NOEC均為1 042.97 μg·L-1。其他階段三唑酮對(duì)其他生物的NOEC均>2 078.88 μg·L-1。方差分析表明，三唑酮僅在試驗(yàn)早期7 d內(nèi)對(duì)微宇宙系統(tǒng)的生物群落有明顯作用，其N(xiāo)OECcommunity為537.98 μg·L-1，這個(gè)值與三唑酮對(duì)大型溞的NOEC值相同，表明三唑酮對(duì)微宇宙系統(tǒng)的NOECcommunity與系統(tǒng)內(nèi)最敏感浮游動(dòng)物的NOECspecies有關(guān)。其余時(shí)間段三唑酮對(duì)微宇宙試驗(yàn)系統(tǒng)的NOECcommunity均>2 078.88 μg·L-1。根據(jù)文獻(xiàn)報(bào)道，在三唑酮使用期間，稻田附近溝渠水三唑酮濃度為12 μg·L-1，這個(gè)值低于本實(shí)驗(yàn)三唑酮暴露初期的NOECcommunity(即537.98 μg·L-1)。因此，三唑酮對(duì)浮游動(dòng)物群落產(chǎn)生的生態(tài)效應(yīng)很小。

表4 三唑酮處理組浮游動(dòng)物對(duì)主成分排序軸的貢獻(xiàn)度及相關(guān)性Table 4 Species contribution and correlation for PCA axis in triadimefon treatment

表5 三唑酮對(duì)浮游動(dòng)物物種及其群落的最大無(wú)作用濃度(NOEC)Table 5 The no observed effect concentration (NOEC) of triadimefon to zooplankton species and its community (μg·L-1)

◆