抗旱和耐鹽堿轉(zhuǎn)基因棉花對(duì)土壤線蟲(chóng)群落的影響

郭佳惠， 李 剛， 趙建寧， 楊殿林， 閆鳳鳴， 修偉明①

(1.農(nóng)業(yè)部環(huán)境保護(hù)科研監(jiān)測(cè)所， 天津 300191； 2.河南農(nóng)業(yè)大學(xué)植物保護(hù)學(xué)院， 河南 鄭州 450002)

土壤是農(nóng)業(yè)生態(tài)系統(tǒng)中物質(zhì)與能量交換的樞紐,土壤生物在土壤物質(zhì)能量遷移轉(zhuǎn)化過(guò)程中具有特殊的功能和作用。土壤線蟲(chóng)在土壤中分布廣泛,以土壤中的細(xì)菌、真菌、其他線蟲(chóng)及植物組織為食,有分解有機(jī)物、養(yǎng)分,維持土壤生態(tài)平衡等重要作用,對(duì)土壤生態(tài)環(huán)境的變化極為敏感。土壤線蟲(chóng)作為土壤生態(tài)系統(tǒng)的重要成員,已成為指示土壤生態(tài)環(huán)境變化的指示生物[1]。

轉(zhuǎn)基因育種是一種新型育種方式,具有周期短、效率高等特點(diǎn),已成為作物育種界最重要的技術(shù)。自從1996年轉(zhuǎn)基因作物首次商品化以來(lái),截至2015年,全球轉(zhuǎn)基因作物的種植面積已達(dá)1.797億hm2,比1996年增加100多倍[2]。隨著人口數(shù)量的增多,耕地面積的不斷減少,為提高土地利用率,應(yīng)對(duì)工業(yè)生產(chǎn)對(duì)棉花質(zhì)量的要求、棉花種植技術(shù)的變革和種植區(qū)域的轉(zhuǎn)移問(wèn)題,通過(guò)生物技術(shù)培育抗旱、耐鹽堿、抗病、抗除草劑、新型抗蟲(chóng)、優(yōu)質(zhì)高產(chǎn)等新性狀棉花品種已成為新需求。自2008年以來(lái),在轉(zhuǎn)基因新品種保護(hù)項(xiàng)目的大力支持下,國(guó)內(nèi)涌現(xiàn)出一批新型轉(zhuǎn)基因棉花新材料[3]。轉(zhuǎn)基因作物的大面積種植可能會(huì)對(duì)自然和農(nóng)業(yè)生態(tài)系統(tǒng)產(chǎn)生負(fù)面影響[4]。轉(zhuǎn)基因作物對(duì)土壤非靶標(biāo)生物的負(fù)面影響便是其中一個(gè)重要的潛在風(fēng)險(xiǎn)。有關(guān)轉(zhuǎn)基因作物對(duì)非靶標(biāo)生物土壤線蟲(chóng)影響的研究已有大量報(bào)道,但結(jié)果卻不盡相同。多數(shù)研究結(jié)果表明轉(zhuǎn)基因作物對(duì)土壤線蟲(chóng)多樣性無(wú)影響[5-8]。但也有研究指出,種植轉(zhuǎn)基因作物會(huì)對(duì)土壤線蟲(chóng)個(gè)別指標(biāo)及屬種產(chǎn)生影響[9-11]。以往有關(guān)轉(zhuǎn)基因作物對(duì)土壤線蟲(chóng)群落影響的研究多集中于轉(zhuǎn)Bt基因抗蟲(chóng)作物[12-13]。由于干旱、鹽堿等性狀屬于多基因控制的原因,所以提高作物抗逆水平一直是轉(zhuǎn)基因研究領(lǐng)域的重點(diǎn)和難點(diǎn),抗旱、耐鹽堿等新性狀轉(zhuǎn)基因作物研發(fā)較晚且多尚處于研發(fā)階段,繼而對(duì)土壤線蟲(chóng)群落影響的研究也相對(duì)較少。因此,為評(píng)價(jià)抗逆轉(zhuǎn)基因棉對(duì)棉田土壤生態(tài)環(huán)境的安全性,筆者通過(guò)試驗(yàn)研究了抗旱、耐鹽堿2種性狀抗逆轉(zhuǎn)基因棉花對(duì)棉田土壤線蟲(chóng)群落的影響,以期為抗逆轉(zhuǎn)基因棉花的推廣應(yīng)用提供依據(jù)。

1 材料與方法

1.1 試驗(yàn)地概況與樣地設(shè)置

試驗(yàn)地位于天津市武清區(qū)的中國(guó)農(nóng)業(yè)科學(xué)院武清轉(zhuǎn)基因生物農(nóng)田生態(tài)系統(tǒng)影響野外科學(xué)觀測(cè)試驗(yàn)站(39°21′50″ N,117°12′53″ E)。該區(qū)域海拔為6.3 m,地處華北平原東北部,屬暖溫帶濕潤(rùn)氣候區(qū)。土壤類型屬潮土,土壤全氮含量為0.56 g·kg-1,全磷含量為0.66 g·kg-1,有機(jī)質(zhì)含量為11.42 g·kg-1。

1.2 供試材料與樣地設(shè)置

抗旱轉(zhuǎn)基因棉花品系013011及其受體對(duì)照棉花品系TH2、耐鹽堿轉(zhuǎn)基因棉花品系013018及其受體對(duì)照棉花品系Q,均由中國(guó)農(nóng)業(yè)科學(xué)院棉花研究所提供。試驗(yàn)采用隨機(jī)區(qū)組設(shè)計(jì),每個(gè)棉花品種設(shè)3個(gè)重復(fù),共12個(gè)小區(qū),每個(gè)小區(qū)規(guī)格均為10 m×10 m,小區(qū)間均設(shè)10 m寬普通玉米隔離帶。棉花全生育期采用常規(guī)耕作管理,不噴施殺蟲(chóng)劑。試驗(yàn)地周圍設(shè)100 m寬普通玉米隔離帶。

1.3 樣品采集

分別于2014和2015年棉花苗期、蕾期、花鈴期、吐絮期和衰老期采樣。采樣時(shí),去除地表雜草和枯枝落葉,分別在各小區(qū)隨機(jī)選取5株棉花,用土壤采樣器緊貼棉花主莖采集0～10 cm表土,混勻后置于冰盒中帶回實(shí)驗(yàn)室,剔除土樣中根系及有機(jī)殘?bào)w后分離土壤線蟲(chóng)。

1.4 土壤線蟲(chóng)的分離與鑒定

采用Baermann淺盤(pán)法分離線蟲(chóng)[11]:稱取50.0 g鮮土置于Baermann淺盤(pán)中的濾紙上,23 ℃條件下分離48 h,然后過(guò)30 μm孔徑鋼篩收集線蟲(chóng),并于60 ℃條件下熱殺死后保存在40 g·L-1福爾馬林溶液中。

用顯微鏡對(duì)線蟲(chóng)進(jìn)行計(jì)數(shù),將獲得數(shù)據(jù)根據(jù)質(zhì)量比例換算為100 g干土中線蟲(chóng)數(shù)量,質(zhì)量比例為總線蟲(chóng)數(shù)與烘干后土壤質(zhì)量比值的100倍。隨機(jī)選擇100條線蟲(chóng),根據(jù)線蟲(chóng)形態(tài)分類鑒定法進(jìn)行分類鑒定[14],如果線蟲(chóng)總數(shù)不足100條,則需全部鑒定。

1.5 數(shù)據(jù)處理與分析

1.5.1線蟲(chóng)群落生態(tài)指數(shù)的計(jì)算

將所鑒定的線蟲(chóng)按食性分為食細(xì)菌線蟲(chóng)(Ba)、食真菌線蟲(chóng)(Fu)、捕食雜食性線蟲(chóng)(OP)和植食性線蟲(chóng)(PP)[15],并計(jì)算其生態(tài)指數(shù)。成熟度指數(shù)(IM)計(jì)算公式:IM=∑vi×fi。其中,vi為自由生活的線蟲(chóng)屬i的c-p值;fi為該線蟲(chóng)屬i的比例[16]。植食性線蟲(chóng)指數(shù)(IPP)計(jì)算公式[16]:IPP=∑vj×fj。其中,vj為植食性線蟲(chóng)屬j的c-p值;fj為該線蟲(chóng)屬j的比例[16]。線蟲(chóng)通路指數(shù)(RNC)計(jì)算公式:RNC=B/(B+F)。其中,B為食細(xì)菌線蟲(chóng)數(shù)量;F為食真菌線蟲(chóng)數(shù)量[16]。Wasilewska指數(shù)(IW)計(jì)算公式:IW=(B+F)/P。其中,P為植食性線蟲(chóng)數(shù)量[17]。

1.5.2統(tǒng)計(jì)分析

所有線蟲(chóng)數(shù)量均轉(zhuǎn)化為100 g干土中的數(shù)量,并計(jì)算每個(gè)線蟲(chóng)屬的相對(duì)豐度和優(yōu)勢(shì)度。所有數(shù)據(jù)采用Excel 2010和SPSS 16.0軟件處理,用單因素分析方法對(duì)各指標(biāo)進(jìn)行方差分析。

2 結(jié)果與分析

2.1 不同處理棉田土壤線蟲(chóng)營(yíng)養(yǎng)類群的豐度變化

在2014—2015年所采集土壤樣品中,共鑒定出土壤線蟲(chóng)34屬,其中,食細(xì)菌線蟲(chóng)13屬,食真菌線蟲(chóng)3屬,捕食性線蟲(chóng)9屬,植食性線蟲(chóng)9屬(表1)。

表1抗旱和耐鹽堿轉(zhuǎn)基因棉及其受體棉棉田土壤線蟲(chóng)組成和相對(duì)豐度

Table1Compositionandrelativeabundanceofsoilnematodesinthesoilsundertransgeniccottonsandtheircorrespondingnon-transgenicone

屬營(yíng)養(yǎng)類群c-p013011TH2013018Q相對(duì)豐度/%分級(jí)相對(duì)豐度/%分級(jí)相對(duì)豐度/%分級(jí)相對(duì)豐度/%分級(jí) 擬麗突屬(Acrobeloides)Ba23.0±0.6++2.8±0.4++4.5±0.8++5.2±0.8++ 無(wú)咽屬(Alaimus)Ba43.2±0.7++3.5±0.6++3.4±1.0++3.1±0.7++ 麗突屬(Acrobeles)Ba20.7±0.2+0.4±0.2+1.4±0.4++0.7±0.2+ 頭葉屬(Cephalobus)Ba223.3±2.6+++24.8±3.2+++21.8±3.0+++22.7±2.4+++ 鹿角唇屬(Cervidellus)Ba21.5±0.3++1.2±0.4++3.1±0.7++1.6±0.4++ 板唇屬(Chiloplacus)Ba20—0—0.2±0.1+0.2±0.1+ 真頭葉屬(Eucephalobus)Ba220.2±2.3+++16.3±1.9+++17.3±2.2+++17.2±2.6+++ 中桿屬(Mesorhabditis)Ba11.0±0.3a+0.4±0.1b+2.0±0.7++1.5±0.6++ 單宮屬(Monhystera)Ba10—0—0—0— 盆咽屬(Panagrolaimus)Ba17.4±1.3++4.7±0.7++6.2±0.7++6.0±0.7++ 棱咽屬(Prismatolainus)Ba30.1±0.1+0.1±0.1+0.5±0.2+0.3±0.1+ 原桿屬(Protorhabditis)Ba14.9±1.4++3.4±0.9++5.9±1.5++6.4±1.2++ 小桿屬(Rhabditis)Ba20—0.1±0.1+0.8±0.5+0.2±0.1+ 真滑刃屬(Aphelenchus)Fu21.4±0.5++1.5±0.5++0.7±0.2+0.7±0.2+ 短矛屬(Doryllium)Fu40.3±0.1+0.2±0.1+0.5±0.3+0.1±0.1+ 莖屬(Ditylenchus)Fu22.8±0.8++2.3±0.6++2.4±0.7++2.0±0.4++ 孔咽屬(Aporcelaimus)OP50.2±0.1+1.0±0.3+0.7±0.3+1.2±0.6++ 真矛線屬(Eudorylaimus)OP40.7±0.3+2.3±1.2++0.7±0.2+1.8±0.6++ 中矛線屬(Mesodorylaimus)OP50—0—0—0— 微矛線(Microdorylaimus)OP40—0—0.5±0.3+0.1±0.1+ 擬矛線(Dorylaimoides)OP40—0.2±0.2+0.2±0.1+0.3±0.1+ 銼齒屬(Mylonchulus)OP40.1±0.1+0.2±0.1+0.1±0.1+0.1±0.1+ 擬桑尼屬(Thorneella)OP40.1±0.1+0—0—0— 盤(pán)咽屬(Discolaimus)OP50.5±0.2+0.5±0.3+0.4±0.2+0.7±0.3+ 桑尼屬(Thornia)OP40.1±0.1A+0.2±0.1B+0.3±0.2+0.3±0.1+ 絲尾墊刃屬(Filenchus)PP20—0—0.1±0.1+0— 螺旋屬(Helicotylenchus)PP324.5±3.7+++28.7±4.2+++20.7±2.7+++22.8±3.1+++ 短體屬(Pratylenchus)PP30.2±0.1+0—0.2±0.1+0.2±0.1+ 矮化屬(Tylenchorhynchus)PP30.7±0.2+0.4±0.2+0.9±0.4+1.0±0.3+ 墊刃屬(Tylenchus)PP30—0.3±0.2+0.3±0.1+0.7±0.3+ 散香屬(Beleodorus)PP21.7±0.5++1.5±0.5++2.8±0.7++1.8±0.4++ 裸矛屬(Psilenchus)PP30.2±0.1+0.2±0.1+0.8±0.4+0.9±0.3+ 偽墊刃屬(Nothotylenchus)PP20—0—0—0— 滑刃屬(Aphelenchoides)PP22.0±0.7++3.0±0.7++1.8±0.4++1.7±0.6++

數(shù)據(jù)為平均值±標(biāo)準(zhǔn)誤。Ba為食細(xì)菌線蟲(chóng),Fu為食真菌線蟲(chóng),OP為捕食雜食性線蟲(chóng),PP為植食性線蟲(chóng)。+++指相對(duì)豐度>10%,優(yōu)勢(shì)屬; ++指相對(duì)豐度介于1%～10%之間,普通屬; +指相對(duì)豐度<1%,稀有屬?！啊北硎緹o(wú)數(shù)據(jù)。同一行英文小寫(xiě)字母不同表示在α=0.05水平上差異顯著,英文大寫(xiě)字母不同表示在α=0.01水平上差異顯著(F檢驗(yàn))。

抗旱轉(zhuǎn)基因棉013011棉田土壤樣品中鑒定出線蟲(chóng)25屬,對(duì)照棉TH2棉田土壤樣品中鑒定出線蟲(chóng)26屬。耐鹽堿轉(zhuǎn)基因棉013018棉田土壤樣品中鑒定出線蟲(chóng)30屬,對(duì)照棉Q棉田土壤樣品中鑒定出線蟲(chóng)29屬。在所有處理棉田中均為優(yōu)勢(shì)屬的有:頭葉屬(Cephalobus)相對(duì)豐度(100 g干土中某屬線蟲(chóng)總數(shù)占總線蟲(chóng)數(shù)的比例)分別為23.3%、24.8%、21.8%和22.7%,真頭葉屬(Eucephalobus)相對(duì)豐度分別為20.2%、16.3%、17.3%和17.2%,螺旋屬(Helicotylenchus)相對(duì)豐度分別為24.5%、28.7%、20.7%和22.8%?？购缔D(zhuǎn)基因棉013011棉田土壤中桿屬(Mesorhabditis)線蟲(chóng)相對(duì)豐度顯著高于其對(duì)照受體棉TH2(P<0.05),抗旱轉(zhuǎn)基因棉田013011棉田土壤桑尼屬(Thornia)相對(duì)豐度顯著低于其對(duì)照受體棉TH2(P<0.01)(表1)。

由表2可知,在2014—2015年間,抗旱轉(zhuǎn)基因棉013011與對(duì)照棉TH2棉田土壤線蟲(chóng)總豐度(豐度為100 g干土中線蟲(chóng)數(shù)量),食真菌線蟲(chóng)、捕食雜食性線蟲(chóng)、植食性線蟲(chóng)豐度均無(wú)顯著差異,但抗旱轉(zhuǎn)基因棉013011棉田土壤食細(xì)菌線蟲(chóng)豐度顯著高于其受體對(duì)照棉(P<0.05)。耐鹽堿轉(zhuǎn)基因棉013018與對(duì)照棉Q棉田土壤線蟲(chóng)總數(shù)及食細(xì)菌線蟲(chóng)、食真菌線蟲(chóng)、捕食雜食性線蟲(chóng)和植食性線蟲(chóng)的豐度無(wú)顯著差異。

方差分析結(jié)果表明,抗旱轉(zhuǎn)基因棉013011和對(duì)照棉TH2棉田土壤線蟲(chóng)總數(shù)和各營(yíng)養(yǎng)類型線蟲(chóng)不受采樣時(shí)間的影響,采樣時(shí)間和性狀間也沒(méi)有交互作用(表3)。但耐鹽堿轉(zhuǎn)基因棉013018與對(duì)照棉Q棉田土壤線蟲(chóng)總數(shù)及食細(xì)菌線蟲(chóng)、捕食雜食性線蟲(chóng)豐度在不同采樣時(shí)間之間均存在顯著或極顯著差異(P<0.05或P<0.01),而采樣時(shí)間和性狀間沒(méi)有交互作用(表3)。

表2轉(zhuǎn)基因棉與對(duì)照棉棉田土壤不同食性線蟲(chóng)的豐度

Table2Abundancesofdifferenttrophicgroupsofsoilnematodesinsoilsundertransgeniccottonandtheircorrespondingnon-transgenicone

轉(zhuǎn)基因棉與其對(duì)照相比,同列英文小寫(xiě)字母相同表示兩者間某指標(biāo)在α=0.05水平上無(wú)顯著差異(F檢驗(yàn)),不同表示存在顯著差異。

表3性狀及采樣時(shí)間對(duì)棉田土壤線蟲(chóng)豐度的方差分析

Table3Varianceanalysisofabundancesofsoilnematodesincottonfieldsrelativetocottontraitandsamplingtime

*和**分別表示在α=0.05和0.01水平上存在顯著差異。

2.2 不同處理棉田土壤線蟲(chóng)群落的生態(tài)指數(shù)變化

方差分析表明,抗旱轉(zhuǎn)基因棉013011與對(duì)照棉TH2處理棉田土壤線蟲(chóng)群落在多樣性指數(shù)(H)、優(yōu)勢(shì)度指數(shù)(C)和均勻度指數(shù)(J)上均無(wú)顯著差異(表4)。但抗旱轉(zhuǎn)基因棉013011棉田土壤線蟲(chóng)屬的物種豐富度(S)顯著低于對(duì)照棉TH2(P<0.05),2個(gè)處理間線蟲(chóng)群落相似度達(dá)0.74±0.02。耐鹽堿轉(zhuǎn)基因棉013018與對(duì)照棉Q棉田土壤線蟲(chóng)的S、H、C和J均無(wú)顯著差異,2個(gè)處理間線蟲(chóng)群落相似度達(dá)0.66±0.02(表4)。4個(gè)處理棉田土壤線蟲(chóng)多樣性指數(shù)(S、H、C和J)在不同采樣時(shí)間之間均無(wú)顯著差異(表5)。

2014—2015年,抗旱、抗鹽堿轉(zhuǎn)基因棉分別與其對(duì)照棉之間棉田土壤線蟲(chóng)群落生態(tài)指數(shù)均無(wú)顯著差異(表4)?？购缔D(zhuǎn)基因棉013011及其對(duì)照棉TH2棉田土壤線蟲(chóng)群落生態(tài)指數(shù)RNC、IM和IPP在不同采樣時(shí)間之間存在極顯著差異(P<0.01,表5)。而耐鹽堿轉(zhuǎn)基因棉田013018及其對(duì)照棉田土壤線蟲(chóng)群落生態(tài)指數(shù)在不同采樣時(shí)間無(wú)顯著差異(表5)。

表4轉(zhuǎn)基因棉與對(duì)照棉田土壤線蟲(chóng)的群落生態(tài)指數(shù)

Table4Soilnematodecommunityecologicalindicesofthesoilsundertransgenicandnon-transgeniccottons

性狀棉花品系SHCJRNCIWIMIPP 抗旱01301111.15±0.54a1.87±0.05a0.22±0.01a0.42±0.02a0.92±0.02a5.18±1.79a1.39±0.07a0.84±0.12a TH212.55±0.34b1.82±0.05a0.24±0.01a0.41±0.01a0.92±0.02a5.73±3.10a1.40±0.11a0.98±0.13a 耐鹽堿01301814.10±1.07a2.04±0.08a0.19±0.02a0.46±0.02a0.94±0.01a3.94±0.76a1.46±0.07a0.77±0.09a Q13.10±0.81a2.01±0.08a0.20±0.02a0.44±0.02a0.95±0.01a4.19±1.05a1.45±0.07a0.83±0.10a

S為物種豐富度指數(shù),H為多樣性指數(shù),C為優(yōu)勢(shì)度指數(shù),J為Pielou均勻度指數(shù),RNC為線蟲(chóng)通路指數(shù),IW為Wasilewska指數(shù),IM為成熟度指數(shù),IPP為植食性線蟲(chóng)指數(shù)。轉(zhuǎn)基因棉與其對(duì)照相比,同列英文小寫(xiě)字母相同表示兩者間某指標(biāo)在α=0.05水平上無(wú)顯著差異(F檢驗(yàn)),不同表示存在顯著差異。

表5性狀和采樣時(shí)間對(duì)棉田土壤線蟲(chóng)群落生態(tài)指數(shù)的影響

Table5Soilnematodecommunityecologicalindicesofthesoilsincottoneldsrelativetocottontraitandsamplingtime

性狀因素SHCJRNCIWIMIPP 抗旱性狀4.61*0.421.630.560.000.020.010.67時(shí)間0.090.670.011.268.07**3.009.98**6.70**性狀?時(shí)間0.021.010.060.250.140.230.730.12 耐鹽堿性狀0.540.110.280.790.100.040.000.14時(shí)間0.440.010.220.260.470.172.281.16性狀?時(shí)間0.660.020.050.221.270.350.670.34

S為物種豐富度指數(shù),H為多樣性指數(shù),C為優(yōu)勢(shì)度指數(shù),J為Pielou均勻度指數(shù),RNC為線蟲(chóng)通路指數(shù),IW為Wasilewska指數(shù),IM為成熟度指數(shù),IPP為植食性線蟲(chóng)指數(shù)。*和**分別表示在α=0.05和0.01水平上存在顯著差異。

3 討論與結(jié)論

通過(guò)2 a的田間調(diào)查研究,2種抗逆轉(zhuǎn)基因棉田土壤總的線蟲(chóng)群落結(jié)構(gòu)及多樣性指數(shù)與其相應(yīng)受體對(duì)照田無(wú)顯著差異,但與對(duì)照田相比,抗旱轉(zhuǎn)基因棉田土壤食細(xì)菌線蟲(chóng)豐度發(fā)生顯著變化,尤其是中桿屬(Mesorhabditis)線蟲(chóng)豐度顯著增加(P<0.05)。由于此次調(diào)查結(jié)果中抗旱轉(zhuǎn)基因棉與其對(duì)照田存在顯著差異的僅為個(gè)別食性及種屬〔表2,食細(xì)菌線蟲(chóng)群落P<0.05;表1,中桿屬(Mesorhabditis)線蟲(chóng)P<0.05,桑尼屬(Thornia)線蟲(chóng)P<0.01〕,推測(cè)這可能是抗旱轉(zhuǎn)基因棉對(duì)棉田線蟲(chóng)的影響存在食性間、種屬間差異,但仍需進(jìn)一步調(diào)查研究。另外,MANACHINI等[21]也曾指出Bt 油菜(Cry1Ac)與非Bt處理組相比,食真菌線蟲(chóng)比例升高,植食性線蟲(chóng)比例降低。戚琳等[13]發(fā)現(xiàn)種植Bt水稻提高了土壤線蟲(chóng)中植食性線蟲(chóng)和捕食雜食性線蟲(chóng)的比例。吳剛等[22]發(fā)現(xiàn)轉(zhuǎn)Cry1Ab/c基因水稻華恢1號(hào)顯著降低中桿屬線蟲(chóng)種群數(shù)量,但顯著增加鉤唇屬(Diploseapter)線蟲(chóng)種群數(shù)量。與常規(guī)對(duì)照親本玉米比較發(fā)現(xiàn),轉(zhuǎn)植酸酶基因玉米處理總的線蟲(chóng)群落多樣性及豐度無(wú)明顯變化,但食細(xì)菌類線蟲(chóng)豐度及各多樣性指數(shù)增加,植食類線蟲(chóng)豐度及各多樣性指數(shù)降低[10]。由此可見(jiàn),有些轉(zhuǎn)基因植物確實(shí)可以對(duì)不同食性線蟲(chóng)豐度產(chǎn)生影響。因此,對(duì)轉(zhuǎn)基因植物生態(tài)安全性進(jìn)行評(píng)價(jià)時(shí),選取敏感食性、種屬的線蟲(chóng)作為評(píng)價(jià)指標(biāo)是值得進(jìn)一步深入研究的。

[1] HAEGEMAN A,MANTELIN S,JONES J T,etal.Functional Roles of Effectors of Plant-Parasitic Nematodes[J].Gene,2012,492(1):19-31.

[2] JAMES C.20th Anniversary (1996 to 2015) of the Global Commercialization of Biotech Crops and Biotech Crop Highlights in 2015[R].Ithaca,NY:ISAAA,2015.

[3] 雒珺瑜,張帥,呂麗敏,等.轉(zhuǎn)ACO2 基因優(yōu)質(zhì)棉生長(zhǎng)特性及其對(duì)田間昆蟲(chóng)群落的影響[J].生態(tài)學(xué)報(bào),2016,36(19):6113-6121.[LUO Jun-yu,ZHANG Shuai,Lü Li-min,etal.Effects of Transgenic ACO2 Cotton on Growth Characteristics and Insect Community[J].Acta Ecologica Sinica,2016,36(19):6113-6121.]

[4] 李建平,肖琴,王吉鵬.轉(zhuǎn)基因作物的風(fēng)險(xiǎn)分析及我國(guó)的應(yīng)對(duì)策略[J].中國(guó)食物與營(yíng)養(yǎng),2013,19(1):5-9.[LI Jian-ping,XIAO Qin,WANG Ji-peng.Risk Analysis for Genetically Modified Crops and China′s Countermeasures[J].Food and Nutrition in China,2013,19(1):5-9.]

[5] SAXENA D,STOTZKY G.Bacillusthuringiensis(Bt) Toxin Released From Root Exudates and Biomass of Bt Corn Has no Apparent Effect on Earthworms,Nematodes,Protozoa,Bacteria,and Fungi in Soil[J].Soil Biology and Biochemistry,2001,33(9):1225-1230.

[6] GRIFFITHS B S,HECKMANN L H,CAUL S,etal.Varietal Effects of Eight Paired Lines of Transgenic Bt Maize and Near-Isogenic Non-Bt Maize on Soil Microbial and Nematode Community Structure[J].Plant Biotechnology Journal,2007,5(1):60-68.

[7] LI X G,LIU B.A 2-Year Field Study Shows Little Evidence that the Long-Term Planting of Transgenic Insect-Resistant Cotton Affects the Community Structure of Soil Nematodes[J].PLoS One,2013,8(4):e61670.

[8] 肖能文,劉勇波,陳群英,等.轉(zhuǎn)基因油菜與野芥菜混作對(duì)土壤線蟲(chóng)群落的影響[J].生態(tài)學(xué)報(bào),2015,35(18):6189-6198.[XIAO Neng-wen,LIU Yong-bo,CHEN Qun-ying,etal.Effect of Interspecific Mixed Cropping Between Transgenic Oilseed Rape and Wild Mustard on Soil Nematode Communities[J].Acta Ecologica Sinica,2015,35(18):6189-6198.]

[9] 楊志國(guó),趙建寧,李剛,等.耐草甘膦轉(zhuǎn)基因大豆對(duì)土壤線蟲(chóng)多樣性的影響[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào),2013,329(11):2199-2205.[YANG Zhi-guo,ZHAO Jian-ning,LI Gang,etal.Effect of Glyphosate-Tolerant Soybean on Diversity of Soil Nematodes[J].Journal of Agro-Environment Science,2013,32(11):2199-2205.

[10] 郭維維,趙宗潮,蘇營(yíng),等.轉(zhuǎn)植酸酶基因玉米種植對(duì)土壤線蟲(chóng)群落的影響[J].應(yīng)用生態(tài)學(xué)報(bào),2014,25(4):1107-1114.[GUO Wei-wei,ZHAO Zong-chao,SU Ying,etal.Effects of Phytase Transgenic Corn Planting on Soil Nematode Community[J].Chinese Journal of Applied Ecology,2014,25(4):1107-1114.]

[11] 風(fēng)春,李剛,修偉明,等.復(fù)合性狀轉(zhuǎn)基因棉花對(duì)土壤線蟲(chóng)多樣性的影響[J].棉花學(xué)報(bào),2015,27(3):268-274.[FENG Chun,LI Gang,XIU Wei-ming,etal.Effects of Genetically Modified Stacked Cottons on Diversity of Soil Nematodes[J].Cotton Science,2015,27(3):268-274.]

[12] YANG B,CHEN H,LIU X,etal.Bt Cotton Planting Does Not Affect the Community Characteristics of Rhizosphere Soil Nematodes[J].Applied Soil Ecology,2014,73:156-164.

[13] 戚琳,陳法軍,劉滿強(qiáng),等.三種轉(zhuǎn)Bt水稻短期種植對(duì)土壤微生物生物量和線蟲(chóng)群落的影響[J].生態(tài)學(xué)雜志,2013,32(4):975-980.[QI Lin,CHEN Fa-jun,LIU Man-qiang,etal.Effects of Short-Term Planting of Three Transgenic Bt Rice Lines on Soil Microbial Biomass and Soil Nematode Assemblage[J].Chinese Journal of Ecology,2013,32(4):975-980.]

[14] 尹文英,胡圣豪,沈韞芬,等.中國(guó)土壤動(dòng)物檢索圖鑒[M].北京:科學(xué)出版社,1988:51-89.[YIN Wen-ying,HU Sheng-hao,SHEN Yun-fen,etal.Pictorial Keys to Soil Animals of China[M].Beijing:Science Press,1988:51-89.]

[15] YEATES G W.Nematodes as Soil Indicators:Functional and Biodiversity Aspects[J].Biology and Fertility of Soils,2003,37(4):199-210.

[16] BONGERS T.The Maturity Index:An Ecological Measure of Environmental Disturbance Based on Nematode Species Composition[J].Oecologia,1990,83(1):14-19.

[17] WASILEWSKA L.The Relationship Between the Diversity of Soil Nematode Communities and the Plant Species Richness of Meadows[J].Ekologia Polska,1997,45(3/4):719-732.

[18] SHANNON C E.A Mathematical Theory of Communication[J].Bell System Technical Journal,1948,27(4):379-423.

[19] YEATES G W,BONGERS T.Nematode Diversity in Agroecosystems[J].Agriculture,Ecosystems & Environment,1999,74(1/2/3):113-135.

[20] PIELOU E C.Mathematical Ecology[M].New York,USA:Wiley,1977：385.

[21] MANACHINI B,LANDI S,FIORE M C,etal.First Investigations on the Effects of Bt-TransgenicBrassicanapusL. on the Trophic Structure of the Nematofauna[J].IOBC/WPRS Bulletin,2004,27:103-108.

[22] 吳剛,李俊生,肖能文,等.轉(zhuǎn)Bt基因水稻對(duì)土壤跳蟲(chóng)、線蟲(chóng)和螨類種群數(shù)量的影響[J].湖北植保,2012(5):6-9.[WU Gang,LI Jun-sheng,XIAO Neng-wen,etal.Effects of the Transgenic Bt Rice Lines on Soil Springtails,Nematodes and Mites Populations[J].Hubei Plant Protection,2012(5):6-9.]