李凡等

摘要：采用直接觀察法調(diào)查了大豆田間節(jié)肢動(dòng)物群落的多樣性，分析比較了節(jié)肢動(dòng)物群落、天敵亞群落和害蟲(chóng)亞群落的多樣性指數(shù)、均勻性指數(shù)、優(yōu)勢(shì)集中性指數(shù)，研究了抗除草劑轉(zhuǎn)基因大豆對(duì)田間節(jié)肢動(dòng)物群落多樣性的影響。結(jié)果表明，在大豆生長(zhǎng)期間，抗除草劑轉(zhuǎn)基因大豆與受體大豆間各指數(shù)均無(wú)顯著性差異。因此，轉(zhuǎn)基因抗除草劑大豆對(duì)田間節(jié)肢動(dòng)物群落多樣性無(wú)明顯影響。

關(guān)鍵詞：抗除草劑轉(zhuǎn)基因大豆；節(jié)肢動(dòng)物群落；多樣性

中圖分類(lèi)號(hào)：Q958.12+2.5文獻(xiàn)標(biāo)識(shí)號(hào)：A文章編號(hào)：1001-4942（2013）07-0083-04

大豆是我國(guó)傳統(tǒng)的糧食油料兼用作物，在農(nóng)業(yè)生產(chǎn)中占有重要地位。隨著基因工程技術(shù)的飛速發(fā)展，全球轉(zhuǎn)基因植物的種植面積越來(lái)越大，尤其是抗除草劑轉(zhuǎn)基因大豆的發(fā)展最為迅速^[1]。1996年孟山都公司推出的抗草甘膦轉(zhuǎn)基因大豆首先獲準(zhǔn)在美國(guó)推廣種植^[2]，至2011年全球抗除草劑轉(zhuǎn)基因大豆的種植面達(dá)到了7 540萬(wàn)公頃^[3]。目前我國(guó)尚未批準(zhǔn)商品化種植轉(zhuǎn)基因大豆，但是用于原料加工的進(jìn)口轉(zhuǎn)基因大豆量正逐年增加，至2009年進(jìn)口量超過(guò)國(guó)產(chǎn)大豆總量，突破4 255萬(wàn)噸^[4]。雖然抗除草劑轉(zhuǎn)基因大豆的種植能夠減少除草劑的使用總量，降低環(huán)境污染，節(jié)約勞動(dòng)成本，但是其生態(tài)安全性問(wèn)題仍然備受爭(zhēng)議，其中對(duì)環(huán)境生物多樣性的影響就是安全評(píng)價(jià)的重要內(nèi)容之一。本研究的目的就是通過(guò)分析抗除草劑轉(zhuǎn)基因大豆對(duì)田間節(jié)肢動(dòng)物群落多樣性的影響，為轉(zhuǎn)基因大豆的環(huán)境安全性評(píng)價(jià)提供科學(xué)參考。

1材料與方法

1.1試驗(yàn)材料

供試品種：抗除草劑轉(zhuǎn)基因大豆HRS、受體大豆RS、普通栽培大豆品種中黃13，由山東省農(nóng)業(yè)科學(xué)院植物保護(hù)研究所提供。

1.2試驗(yàn)地點(diǎn)

試驗(yàn)地點(diǎn)設(shè)在農(nóng)業(yè)部轉(zhuǎn)基因植物環(huán)境安全監(jiān)督檢驗(yàn)測(cè)試中心（濟(jì)南）試驗(yàn)基地，位于東經(jīng)117°5′27.42″，北緯 36°42′59.58″，屬于暖溫帶大陸性季風(fēng)氣候，年平均氣溫14.7℃，年平均降水量671.1 mm，年無(wú)霜期235天。試驗(yàn)田地勢(shì)平坦，土壤為中性壤土，肥力中等；前茬作物為蔬菜，試驗(yàn)地四周有2.5 m高的隔離圍墻，周?chē)? 000 m范圍內(nèi)無(wú)其他大豆種植。

1.3試驗(yàn)方法

試驗(yàn)設(shè)3個(gè)處理，處理1：抗除草劑轉(zhuǎn)基因大豆HRS；處理2：受體大豆RS；處理3：普通栽培大豆品種中黃13。隨機(jī)區(qū)組設(shè)計(jì)，重復(fù)4次，小區(qū)面積為150 m2，小區(qū)間設(shè)有2 m寬隔離帶。

2011年5月24日播種，人工撒播，播種深度3～4 cm，播種量10～12 g/m2，行距60 cm，株距5～6 cm。按當(dāng)?shù)爻Ｒ?guī)栽培模式進(jìn)行管理，全生育期田間管理除必要的農(nóng)事操作（如鋤草、澆水）外，不噴施農(nóng)藥，以減少人為影響，保證大豆植株正常生長(zhǎng)。

1.4調(diào)查方法

從定苗到成熟，每10天調(diào)查1次。每小區(qū)采用對(duì)角線5點(diǎn)取樣，每點(diǎn)系統(tǒng)調(diào)查相鄰的20株大豆上各種節(jié)肢動(dòng)物的種類(lèi)和數(shù)量。

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

統(tǒng)計(jì)分析各個(gè)調(diào)查時(shí)期不同群落（節(jié)肢動(dòng)物群落、害蟲(chóng)亞群落和天敵亞群落）的結(jié)構(gòu)與組成，計(jì)算出3個(gè)處理各時(shí)期各群落的多樣性指數(shù)、均勻性指數(shù)和優(yōu)勢(shì)集中性指數(shù)。計(jì)算公式如下：

多樣性指數(shù)H=–∑Si=1PilnPi

均勻性指數(shù)J=H/lnS

優(yōu)勢(shì)集中性指數(shù) C=∑Si=1（Ni/N）2

式中：Pi=Ni/N，Ni為第i個(gè)物種的個(gè)體數(shù)，N為群落中所有物種的總個(gè)體數(shù)，S為群落中物種數(shù)。

利用Excel 2003和DPS統(tǒng)計(jì)軟件對(duì)試驗(yàn)數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析，用Duncans新復(fù)極差法進(jìn)行多重比較。

2結(jié)果與分析

2.1對(duì)大豆田節(jié)肢動(dòng)物群落多樣性的影響分析

在大豆生長(zhǎng)期，3個(gè)大豆品種的田間節(jié)肢動(dòng)物群落的多樣性指數(shù)（圖1）、均勻性指數(shù)（圖2）以及優(yōu)勢(shì)集中性指數(shù)（圖3）的變化趨勢(shì)均基本一致。3個(gè)大豆品種的多樣性指數(shù)和均勻性指數(shù)動(dòng)態(tài)均呈現(xiàn)先升后降的趨勢(shì)，而優(yōu)勢(shì)集中性指數(shù)動(dòng)態(tài)變化趨勢(shì)剛好相反，呈先下降后上升態(tài)勢(shì)。在大豆生長(zhǎng)早期，田間節(jié)肢動(dòng)物種類(lèi)較少，多樣性指數(shù)和均勻性指數(shù)較低，而優(yōu)勢(shì)集中性指數(shù)較高。在6月下旬～7月下旬大豆的生長(zhǎng)旺盛期，田間節(jié)肢動(dòng)物種類(lèi)較多，多樣性指數(shù)較高，均在1.2以上；均勻性指數(shù)也較高，均在0.5以上；而優(yōu)勢(shì)集中性指數(shù)較低，均在0.5以下。進(jìn)入9月份隨著大豆生長(zhǎng)期的結(jié)束，節(jié)肢動(dòng)物種類(lèi)逐漸減少，多樣性指數(shù)和均勻性指數(shù)逐步降低，而優(yōu)勢(shì)集中性指數(shù)升高。

在各個(gè)調(diào)查時(shí)期，分別對(duì)3個(gè)處理田間節(jié)肢動(dòng)物群落的生物多樣性指數(shù)、均勻性指數(shù)和優(yōu)勢(shì)集中性指數(shù)進(jìn)行差異顯著性分析，結(jié)果表明整個(gè)調(diào)查期間，轉(zhuǎn)基因大豆HRS與受體大豆RS、普通栽培大豆中黃13之間各指數(shù)均無(wú)顯著性差異，說(shuō)明抗除草劑轉(zhuǎn)基因大豆HRS對(duì)田間節(jié)肢動(dòng)物群落多樣性沒(méi)有明顯影響。

2.2對(duì)大豆田天敵亞群落多樣性的影響分析

由圖4～圖6可以看出，在大豆整個(gè)生長(zhǎng)期，3個(gè)大豆品種田間天敵亞群落的多樣性指數(shù)、均勻性指數(shù)以及優(yōu)勢(shì)集中性指數(shù)的大體趨勢(shì)是一致的。3個(gè)處理的多樣性指數(shù)在大豆生長(zhǎng)早期均較低，在生長(zhǎng)中期較高（在0.8～1.2之間）且趨于穩(wěn)定，生長(zhǎng)后期又有所下降；均勻性指數(shù)隨時(shí)間的推移呈先升后降的趨勢(shì)，而優(yōu)勢(shì)集中性指數(shù)呈先降后升的態(tài)勢(shì)，這主要是由于大豆田間苗期和成熟期天敵物種數(shù)較少。在各個(gè)調(diào)查時(shí)期，分別分析比較各處理田間天敵亞群落的多樣性指數(shù)、均勻性指數(shù)和優(yōu)勢(shì)集中性指數(shù)，結(jié)果表明轉(zhuǎn)基因大豆HRS與受體大豆RS、普通栽培大豆中黃13之間各指數(shù)均差異不顯著，說(shuō)明抗除草劑轉(zhuǎn)基因大豆HRS對(duì)田間天敵亞群落多樣性沒(méi)有明顯影響。

2.3對(duì)大豆田害蟲(chóng)亞群落多樣性的影響分析

通過(guò)調(diào)查分析大豆田間害蟲(chóng)亞群落發(fā)現(xiàn)，3個(gè)處理的多樣性指數(shù)（圖7）、均勻性指數(shù)（圖8）以及優(yōu)勢(shì)集中性指數(shù)（圖9）在整個(gè)生長(zhǎng)期間的變化趨勢(shì)也基本一致。3個(gè)大豆品種的害蟲(chóng)亞群落的多樣性指數(shù)動(dòng)態(tài)均表現(xiàn)為先升后降，在6月下旬～7月間達(dá)到高峰；均勻性指數(shù)的總體趨勢(shì)是大豆生長(zhǎng)中期比早期和晚期都要高，這主要是由于生長(zhǎng)中期田間的害蟲(chóng)種類(lèi)較多。而優(yōu)勢(shì)集中性指數(shù)呈先降后升的態(tài)勢(shì)，生長(zhǎng)早期和后期調(diào)查的害蟲(chóng)物種數(shù)較少，優(yōu)勢(shì)集中性指數(shù)較高，隨著大豆生長(zhǎng)旺盛期田間害蟲(chóng)種類(lèi)增多優(yōu)勢(shì)集中性指數(shù)降低。

在各個(gè)調(diào)查時(shí)期，分別對(duì)3個(gè)大豆品種田間害蟲(chóng)亞群落的多樣性指數(shù)、均勻性指數(shù)和優(yōu)勢(shì)集中性指數(shù)進(jìn)行分析比較，結(jié)果表明轉(zhuǎn)基因大豆HRS與受體大豆RS、普通栽培大豆中黃13之間各指數(shù)均無(wú)顯著性差異，說(shuō)明抗除草劑轉(zhuǎn)基因大豆HRS對(duì)田間害蟲(chóng)亞群落多樣性也沒(méi)有明顯影響。

3討論

本試驗(yàn)通過(guò)對(duì)大豆田間主要節(jié)肢動(dòng)物群落的定期調(diào)查，發(fā)現(xiàn)豆田主要害蟲(chóng)有蚜蟲(chóng)、粉虱、灰飛虱、葉蟬、薊馬等，主要天敵有瓢蟲(chóng)、花蝽、獵蝽、草蛉、蜘蛛等?？钩輨┺D(zhuǎn)基因大豆田間害蟲(chóng)和天敵的發(fā)生種類(lèi)、數(shù)量與受體大豆之間沒(méi)有明顯差異，這與吳奇等（2007）^[5]對(duì)豆田主要害蟲(chóng)的調(diào)查結(jié)果基本一致。

轉(zhuǎn)基因作物的生態(tài)安全性問(wèn)題是目前國(guó)內(nèi)外學(xué)者研究的熱點(diǎn)。劉志誠(chéng)^[6，7]對(duì)轉(zhuǎn)Bt基因水稻的研究顯示，Bt水稻對(duì)稻田節(jié)肢動(dòng)物群落結(jié)構(gòu)無(wú)明顯的負(fù)面影響。吳奇等（2008）^[8]通過(guò)直接觀察法和吸蟲(chóng)器法對(duì)抗草甘膦轉(zhuǎn)基因大豆AG5601以及張卓等（2011）^[9]采用直接觀察法對(duì)轉(zhuǎn)基因耐草甘膦大豆呼交03-263和06-698的研究都表明轉(zhuǎn)基因大豆對(duì)豆田節(jié)肢動(dòng)物群落無(wú)明顯影響。McPherson等（2003）^[10]的研究發(fā)現(xiàn)轉(zhuǎn)基因耐除草劑大豆對(duì)昆蟲(chóng)群落的影響較小。Jasinski等（2003）^[11]的研究也表明轉(zhuǎn)基因大豆對(duì)非靶標(biāo)節(jié)肢動(dòng)物群落沒(méi)有影響。本文通過(guò)對(duì)大豆田間節(jié)肢動(dòng)物群落、天敵亞群落以及害蟲(chóng)亞群落的主要參數(shù)（多樣性指數(shù)、均勻性指數(shù)、優(yōu)勢(shì)集中性指數(shù)）的分析比較，表明轉(zhuǎn)基因大豆與受體品種間各參數(shù)均無(wú)顯著性差異，這與前人的研究結(jié)果基本一致^[8，9]，說(shuō)明轉(zhuǎn)基因抗除草劑大豆對(duì)田間節(jié)肢動(dòng)物群落多樣性沒(méi)有明顯影響。

隨著國(guó)內(nèi)對(duì)進(jìn)口轉(zhuǎn)基因大豆需求量的不斷擴(kuò)大，中國(guó)作為大豆的原產(chǎn)地和品種多樣性的集中地，更應(yīng)該加強(qiáng)對(duì)大豆品種資源的保護(hù)，增強(qiáng)對(duì)轉(zhuǎn)基因大豆生產(chǎn)的安全監(jiān)控和評(píng)價(jià)，建立完善的風(fēng)險(xiǎn)評(píng)估體系，從而正確認(rèn)識(shí)轉(zhuǎn)基因作物的安全性問(wèn)題，保障轉(zhuǎn)基因作物的健康發(fā)展。

參考文獻(xiàn)：

[1]沈曉峰， 欒鳳俠， 陶波. 抗草甘膦轉(zhuǎn)基因大豆生物與環(huán)境安全性[J]. 東北農(nóng)業(yè)大學(xué)學(xué)報(bào)， 2007， 38（3）： 401-404.

[2]周波， 陶波， 欒鳳俠， 等. 抗草甘膦轉(zhuǎn)基因大豆生物安全性綜述[J]. 作物雜志， 2006， 2： 7-9.

[3]James C. Global status of commercialized biotech/GM crops： 2011 [R]. ISAAA， 2011， 43.

[4]索海翠， 熊瑞權(quán)， 馬啟斌， 等. 轉(zhuǎn)基因大豆的潛在風(fēng)險(xiǎn)及發(fā)展趨勢(shì)研究[J]. 廣東農(nóng)業(yè)科學(xué)， 2010， 10： 244-245.

[5]吳奇， 彭煥， 彭可維， 等. 抗除草劑轉(zhuǎn)基因大豆對(duì)豆田主要害蟲(chóng)發(fā)生動(dòng)態(tài)的影響[J]. 植物保護(hù)， 2007， 33（5）： 50-53.

[6]劉志誠(chéng). 轉(zhuǎn)cry1Ab/cry1Ac基因秈稻對(duì)稻田節(jié)肢動(dòng)物群落的影響[J]. 昆蟲(chóng)學(xué)報(bào)， 2003， 46（4）： 454-465.

[7]劉志誠(chéng). 轉(zhuǎn)cry1Ab基因粳稻對(duì)稻田節(jié)肢動(dòng)物群落的影響[J]. 生物安全學(xué)報(bào)， 2011， 20（1）： 69-76.

[8]吳奇， 彭德良， 彭于發(fā). 抗草甘膦轉(zhuǎn)基因大豆對(duì)非靶標(biāo)節(jié)肢動(dòng)物群落多樣性的影響[J]. 生態(tài)學(xué)報(bào)， 2008， 28（6）： 2622-2628.

[9]張卓， 黃文坤， 劉茂炎， 等. 轉(zhuǎn)基因耐草甘膦大豆對(duì)豆田節(jié)肢動(dòng)物群落多樣性的影響[J]. 植物保護(hù)， 2011， 37（6）： 115-119.

[10]McPherson R， Johnson W， Mullinix B J， et al. Influence of herbicide tolerant soybean production system on insect pest populations and pest-induced crop damage [J]. Journal of Economic Entomology， 2003， 96（3）： 690-698.

[11]Jasinski J R， Eisley J B， Young C E， et al. Select nontarget arthropod abundance in transgenic and nontransgenic field crops in Ohio[J]. Environmental Entomology， 2003， 32（2）： 407-413.