轉(zhuǎn)雙價(jià)基因棉SGK321不同秸稈還田量對(duì)土壤線(xiàn)蟲(chóng)群落的影響

李 靜，蔡連賀，劉瑞華，風(fēng) 春，李 剛，楊殿林，趙建寧*

（1.東北農(nóng)業(yè)大學(xué)資源與環(huán)境學(xué)院，哈爾濱 150030；2.農(nóng)業(yè)農(nóng)村部環(huán)境保護(hù)科研監(jiān)測(cè)所，農(nóng)業(yè)農(nóng)村部產(chǎn)地環(huán)境污染防控重點(diǎn)實(shí)驗(yàn)室/天津市農(nóng)業(yè)環(huán)境與農(nóng)產(chǎn)品安全重點(diǎn)實(shí)驗(yàn)室，天津 300191；3.中國(guó)農(nóng)業(yè)科學(xué)院棉花研究所，河南 安陽(yáng) 455000；4.呼和浩特第九中學(xué)，呼和浩特 010000）

全球轉(zhuǎn)基因作物種植面積不斷擴(kuò)大，已經(jīng)從1996年的170萬(wàn)hm2上升至2017年的1.898億hm2[1]。棉花是全球主要經(jīng)濟(jì)、纖維作物之一，2018年，我國(guó)轉(zhuǎn)基因棉花的種植面積約300萬(wàn)hm2，占全國(guó)棉花播種總面積的80%以上，是我國(guó)種植面積最大的轉(zhuǎn)基因作物。但轉(zhuǎn)基因作物在為社會(huì)帶來(lái)巨大經(jīng)濟(jì)效益的同時(shí)，其潛在的生態(tài)環(huán)境安全風(fēng)險(xiǎn)也被越來(lái)越多的人們所關(guān)注[2-3]。

土壤線(xiàn)蟲(chóng)是地球上數(shù)量豐富的后生動(dòng)物之一，在維持土壤生態(tài)系統(tǒng)穩(wěn)定、促進(jìn)物質(zhì)循環(huán)和能量流動(dòng)等方面發(fā)揮著重要作用[4]，常被作為監(jiān)測(cè)農(nóng)業(yè)管理措施干擾和評(píng)價(jià)土壤環(huán)境質(zhì)量的重要指示生物[5-7]。轉(zhuǎn)Bt作物的目的基因可以通過(guò)花粉、秸稈還田、根系分泌等方式進(jìn)入土壤[8-9]，會(huì)大幅增加土壤中Bt蛋白的含量[10-12]，進(jìn)而可能會(huì)對(duì)土壤動(dòng)物如土壤線(xiàn)蟲(chóng)產(chǎn)生影響。段杉[13]研究發(fā)現(xiàn)，SGK321與石遠(yuǎn)321相比在蕾期和吐絮期的線(xiàn)蟲(chóng)數(shù)量存在顯著差異，但其他一些研究并未發(fā)現(xiàn)轉(zhuǎn)Bt棉花對(duì)土壤線(xiàn)蟲(chóng)群落有顯著影響[14-16]。對(duì)轉(zhuǎn)Bt基因玉米MON810的研究表明，其土壤線(xiàn)蟲(chóng)密度顯著高于常規(guī)玉米[17]。Gu等[18]研究發(fā)現(xiàn)秸稈還田對(duì)總豐度無(wú)顯著影響，但隨著秸稈還田量的增加，食細(xì)菌性線(xiàn)蟲(chóng)豐度逐漸增高，雜食性線(xiàn)蟲(chóng)豐度逐漸減少。有研究發(fā)現(xiàn)，秸稈還田后土壤食細(xì)菌性線(xiàn)蟲(chóng)、食真菌性線(xiàn)蟲(chóng)和有較高c-p值的雜食/捕食性線(xiàn)蟲(chóng)數(shù)量增加，而植食性線(xiàn)蟲(chóng)數(shù)量降低[19-20]。牟文雅等[21]研究發(fā)現(xiàn)秸稈對(duì)線(xiàn)蟲(chóng)數(shù)量及物種數(shù)有一定影響，可改變線(xiàn)蟲(chóng)群落組成。國(guó)內(nèi)外已有很多文獻(xiàn)報(bào)道了轉(zhuǎn)基因作物對(duì)土壤線(xiàn)蟲(chóng)群落結(jié)構(gòu)的影響和傳統(tǒng)作物秸稈還田對(duì)土壤線(xiàn)蟲(chóng)群落結(jié)構(gòu)的影響。秸稈還田是現(xiàn)代農(nóng)業(yè)生產(chǎn)過(guò)程中的一項(xiàng)重要技術(shù)，還田后的秸稈腐解對(duì)于促進(jìn)作物生長(zhǎng)、改善和提高土壤質(zhì)量具有重要作用。因此研究轉(zhuǎn)基因棉花秸稈還田腐解后對(duì)土壤線(xiàn)蟲(chóng)群落結(jié)構(gòu)的影響具有十分重要的意義。

本研究以轉(zhuǎn)雙價(jià)（Bt+CpTI）基因棉SGK321及其常規(guī)棉石遠(yuǎn)321為對(duì)象，比較不同棉花秸稈還田量下、不同生育期中兩種棉花土壤線(xiàn)蟲(chóng)密度及其群落的動(dòng)態(tài)變化情況，探討轉(zhuǎn)雙價(jià)（Bt+CpTI）基因棉秸稈還田對(duì)土壤線(xiàn)蟲(chóng)群落結(jié)構(gòu)的影響，以期深入了解轉(zhuǎn)基因棉對(duì)土壤生態(tài)環(huán)境的影響和構(gòu)建轉(zhuǎn)基因棉生態(tài)安全監(jiān)測(cè)指標(biāo)，為轉(zhuǎn)基因棉種植的土壤生態(tài)安全性評(píng)價(jià)提供科學(xué)依據(jù)。

1 材料與方法

1.1 試驗(yàn)材料

供試棉花品種為轉(zhuǎn)雙價(jià)（Bt+CpTI）基因棉SGK321及其常規(guī)棉石遠(yuǎn)321（簡(jiǎn)稱(chēng)SGK321和石遠(yuǎn)321），均由中國(guó)農(nóng)業(yè)科學(xué)院植物保護(hù)研究所提供。

供試花盆上口直徑20 cm、底部直徑16 cm、高23 cm。供試土壤為潮土，取自天津市武清區(qū)常規(guī)農(nóng)田土壤，其基本理化性質(zhì)：全磷含量0.79 g·kg-1，全氮含量 0.63 g·kg-1，有機(jī)質(zhì)含量 18 g·kg-1，pH 8.24。每盆裝過(guò)2 mm篩后充分混勻的干基土5.6 kg，為模擬正常土壤溫度，將花盆埋入土壤中，埋深15 cm。

1.2 試驗(yàn)設(shè)計(jì)

試驗(yàn)在農(nóng)業(yè)農(nóng)村部環(huán)境保護(hù)科研監(jiān)測(cè)所（天津市南開(kāi)區(qū)，39°5′8.28″N，117°8′39.12″E）網(wǎng)室內(nèi)進(jìn)行。試驗(yàn)區(qū)屬溫帶大陸性季風(fēng)氣候，年均氣溫和年均降水量分別為15℃和600 mm。該試驗(yàn)是不完全隨機(jī)區(qū)組實(shí)驗(yàn)，設(shè)置5個(gè)還田梯度，按照土層厚度換算而來(lái)。每公頃0～20 cm耕層土壤約有2.25×106kg，還田梯度分別為 0、2250、4500、6750 kg·hm-2和 9000 kg·hm-2，分別代表不還田、半量還田、全量還田、1.5倍還田量和2倍還田量。換算成盆栽試驗(yàn)的還田梯度為0（對(duì)照）、1.0、2.0、3.0 g·kg-1和 4.0 g·kg-1。試驗(yàn)所用秸稈為SGK321和石遠(yuǎn)321的秸稈，用粉碎機(jī)粉碎后過(guò)20目篩，按照試驗(yàn)設(shè)置分別與土壤均勻混合，用于轉(zhuǎn)基因棉花SGK321及其常規(guī)棉石遠(yuǎn)321棉花的種植。每個(gè)還田處理設(shè)有4個(gè)重復(fù)，每個(gè)重復(fù)5盆，共200盆盆栽。于2015年4月24日播種，每盆放置6粒棉花種子，待出苗后定苗至每盆2株，日常管理過(guò)程中不噴灑農(nóng)藥。分別于播種后60 d（蕾期）和90 d（花鈴期）兩個(gè)時(shí)期進(jìn)行土壤樣品采集。

1.3 土壤樣品的采集

去除花盆表面雜草和落葉，用土鉆取距離棉花主莖2 cm的土壤，深度15 cm，每盆取2鉆，5盆混為一個(gè)樣品。樣品混勻后除去植物根系、殘?jiān)?，放入自封袋帶回?shí)驗(yàn)室，于4℃保存，并在兩周內(nèi)完成線(xiàn)蟲(chóng)分離。

1.4 線(xiàn)蟲(chóng)的分離與鑒定

土壤線(xiàn)蟲(chóng)分離采用Baermann淺盤(pán)法[22]：從采集的混勻土壤中取出50 g鮮土，在25℃的環(huán)境下靜止分離48 h，分離后的土壤浸出液用500目篩過(guò)濾，將分離出的線(xiàn)蟲(chóng)經(jīng)60℃水浴殺死后，保存于4%的福爾馬林溶液中[23]。根據(jù)土壤的含水量，折合成每100 g干土土壤線(xiàn)蟲(chóng)密度，然后每個(gè)樣品根據(jù)鏡檢法[24]隨機(jī)選取100條線(xiàn)蟲(chóng)，使用Nikon公司生產(chǎn)的TS100熒光側(cè)置顯微鏡進(jìn)行鑒定。如果樣品中線(xiàn)蟲(chóng)不足100條，則所有樣品都要鑒定。根據(jù)土壤線(xiàn)蟲(chóng)的頭部形態(tài)學(xué)特征，將線(xiàn)蟲(chóng)劃分為植食性線(xiàn)蟲(chóng)（Plant-parasites）、食細(xì)菌性線(xiàn)蟲(chóng)（Bacterivores）、食真菌性線(xiàn)蟲(chóng)（Fungivores）和雜食/捕食性線(xiàn)蟲(chóng)（Omnivores-predators）等4個(gè)營(yíng)養(yǎng)類(lèi)群；并根據(jù)線(xiàn)蟲(chóng)的食性和生活策略（r-策略和k-策略）將不同線(xiàn)蟲(chóng)類(lèi)群賦予1～5的c-p值[25]。

1.5 線(xiàn)蟲(chóng)生態(tài)指標(biāo)分析

研究采用生態(tài)學(xué)評(píng)價(jià)指數(shù)：

（1）Shannon-Wiener多樣性指數(shù)：H′=-∑PilnPi，Pi=ni/N。式中Pi為樣品中屬于第i種的個(gè)體的比例；ni為第i類(lèi)群的個(gè)體數(shù)；N為所有類(lèi)群的個(gè)體總數(shù)[26]。

（2）Pielou 均勻度指數(shù)：J=H′/lnS。式中H′為Shannon-Wiener多樣性指數(shù)；S為鑒定分類(lèi)單元的數(shù)目[27]。

（3）Margalef豐富度指數(shù)：SR=（S-1）/lnN。式中S為鑒定分類(lèi)單元的數(shù)目；N為鑒定的線(xiàn)蟲(chóng)個(gè)體數(shù)目[28]。

（4）線(xiàn)蟲(chóng)通路指數(shù)：NCR=NB/（NB+NF）。式中NB指食細(xì)菌性線(xiàn)蟲(chóng)的土壤線(xiàn)蟲(chóng)密度；NF指食真菌性線(xiàn)蟲(chóng)的土壤線(xiàn)蟲(chóng)密度[29]。

（5）瓦斯樂(lè)斯卡指數(shù)：WI=（NF+NB）/NPP。式中NF和NB同上；NPP為植食性線(xiàn)蟲(chóng)密度[30]。

（6）自由生活線(xiàn)蟲(chóng)成熟度指數(shù)：MI=∑（vi×fi）。式中vi為自由生活線(xiàn)蟲(chóng)的c-p值；fi為第i種線(xiàn)蟲(chóng)的個(gè)體數(shù)占自由生活線(xiàn)蟲(chóng)的比例[31]。

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

數(shù)據(jù)分析采用SPSS軟件（16.0版，SPSS Inc．），采用T檢驗(yàn)分析轉(zhuǎn)基因棉與常規(guī)棉之間的差異，不同還田處理與不同生育期生態(tài)指數(shù)差異采用單因子方差分析法（One-way ANOVA），品種、秸稈還田量和生育期之間的交互作用采用多因子方差分析法（Threeway ANOVA），處理組間顯著檢驗(yàn)采用Duncan法。

2 結(jié)果與分析

2.1 土壤線(xiàn)蟲(chóng)密度

SGK321土壤線(xiàn)蟲(chóng)密度僅在蕾期3.0 g·kg-1還田處理時(shí)顯著高于石遠(yuǎn)321（P＜0.05），其他處理差異不顯著（圖1）。不還田時(shí)，兩個(gè)時(shí)期SGK321土壤線(xiàn)蟲(chóng)密度均低于石遠(yuǎn)321，但無(wú)顯著差異。

如圖1所示，蕾期，隨著秸稈還田量的增加，石遠(yuǎn)321土壤線(xiàn)蟲(chóng)密度逐漸增加，在4.0 g·kg-1還田處理時(shí)密度達(dá)到最大值。SGK321土壤線(xiàn)蟲(chóng)密度在3.0 g·kg-1還田處理時(shí)達(dá)到最大值，且秸稈還田后土壤線(xiàn)蟲(chóng)密度均高于不還田處理?；ㄢ徠?，隨著秸稈還田量的增加，石遠(yuǎn)321和SGK321土壤線(xiàn)蟲(chóng)密度均逐漸增加，在4.0 g·kg-1還田處理時(shí)密度達(dá)到最大值，同時(shí)SGK321的線(xiàn)蟲(chóng)密度均高于石遠(yuǎn)321，但無(wú)顯著差異。

花鈴期土壤線(xiàn)蟲(chóng)密度相比于蕾期，石遠(yuǎn)321土壤線(xiàn)蟲(chóng)密度在0、2.0 g·kg-1處理中減少，其他處理中增加；SGK321土壤線(xiàn)蟲(chóng)密度在3.0 g·kg-1還田處理中減少，其他處理中增加。

2.2 土壤線(xiàn)蟲(chóng)群落組成

如表1和表2所示，SGK321兩個(gè)時(shí)期共鑒定出土壤線(xiàn)蟲(chóng)19科34屬，其中食細(xì)菌性線(xiàn)蟲(chóng)13屬，食真菌性線(xiàn)蟲(chóng)3屬，植食性線(xiàn)蟲(chóng)11屬，雜食/捕食性線(xiàn)蟲(chóng)7屬。石遠(yuǎn)321兩個(gè)時(shí)期共鑒定出土壤線(xiàn)蟲(chóng)19科37屬，其中食細(xì)菌性線(xiàn)蟲(chóng)13屬，食真菌性線(xiàn)蟲(chóng)3屬，植食性線(xiàn)蟲(chóng)12屬，雜食/捕食性線(xiàn)蟲(chóng)9屬。SGK321和石遠(yuǎn)321的優(yōu)勢(shì)屬均是螺旋屬Helicotylenchus，其中SGK321土壤中螺旋屬在蕾期2.0 g·kg-1與1.0、3.0、4.0 g·kg-1還田處理之間存在顯著差異（P＜0.05）。石遠(yuǎn)321土壤中螺旋屬在花鈴期4.0 g·kg-1與1.0、2.0 g·kg-1還田處理之間存在顯著差異（P＜0.05）。

蕾期，擬麗突屬Acrobeloides、原桿屬Protorhabditis等6屬在兩個(gè)棉花品種各還田處理之間均存在顯著差異（P＜0.05）；SGK321土壤線(xiàn)蟲(chóng)中莖屬Ditylenchus、墊刃屬Tylenchus、異皮屬Heterodera、銼齒屬M(fèi)ylonchulus等10屬在各還田處理之間存在顯著差異（P＜0.05），石遠(yuǎn)321土壤線(xiàn)蟲(chóng)中孔咽屬Aporcelaimus、真矛線(xiàn)屬Eudorylaimus等7屬在各還田處理之間存在顯著差異（P＜0.05）?；ㄢ徠?，鹿角唇屬Cervidellus、莖屬等5屬在兩個(gè)棉花品種各還田處理之間均存在顯著差異（P＜0.05）；SGK321土壤線(xiàn)蟲(chóng)中絲尾墊刃屬Filenchus和原桿屬在各還田處理之間存在顯著差異（P＜0.05）；石遠(yuǎn)321土壤線(xiàn)蟲(chóng)中棱咽屬Prismatolaimus、矮化屬Tylenchorhynchus等10屬在各還田處理之間存在顯著差異（P＜0.05）。石遠(yuǎn)321土壤線(xiàn)蟲(chóng)中擬麗突屬、棱咽屬和矮化屬等5屬在兩個(gè)生育期各還田處理中均存在顯著差異（P＜0.05）；SGK321土壤線(xiàn)蟲(chóng)中絲尾墊刃屬、原桿屬、鹿角唇屬、墊刃屬和莖屬在兩個(gè)生育期各還田處理中均存在顯著差異（P＜0.05）。

與不還田相比，三等齒屬Pucephalobus、麗突屬Acrobeles等11屬只出現(xiàn)在秸稈還田處理中。其中劍屬Xiphinema、三等齒屬Pucephalobus、微矛線(xiàn)屬M(fèi)icrodorylaimus、前矛線(xiàn)屬Prodorylaimus僅出現(xiàn)在石遠(yuǎn)321的還田處理中（表1），麗突屬僅出現(xiàn)在SGK321的4.0 g·kg-1還田處理中（表2），其他屬分布在兩個(gè)時(shí)期的不同還田處理中。

2.3 土壤線(xiàn)蟲(chóng)的營(yíng)養(yǎng)類(lèi)群結(jié)構(gòu)

SGK321和石遠(yuǎn)321土壤線(xiàn)蟲(chóng)在采樣的兩個(gè)生育期均以食細(xì)菌性線(xiàn)蟲(chóng)（蕾期：27.07%～39.93%；花鈴期：32.73%～52.7%）和植食性線(xiàn)蟲(chóng)（蕾期：25.63%～55.03%；花鈴期：14.98%～42.89%）為主要營(yíng)養(yǎng)類(lèi)群（圖2）。不還田時(shí)，蕾期SGK321中食細(xì)菌性和食真菌性線(xiàn)蟲(chóng)豐度高于石遠(yuǎn)321，雜食/捕食性和植食性線(xiàn)蟲(chóng)豐度低于石遠(yuǎn)321；花鈴期SGK321中雜食/捕食性線(xiàn)蟲(chóng)豐度顯著高于石遠(yuǎn)321（P＜0.05）。

蕾期，在2.0 g·kg-1還田處理時(shí)，除石遠(yuǎn)321土壤中食細(xì)菌線(xiàn)蟲(chóng)豐度顯著高于SGK321（P＜0.05），食真菌線(xiàn)蟲(chóng)豐度高于SGK321，但無(wú)顯著差異外，其他處理SGK321土壤中食細(xì)菌性和食真菌性線(xiàn)蟲(chóng)豐度均高于石遠(yuǎn)321，但無(wú)顯著差異；兩種棉花土壤中食細(xì)菌性和食真菌性線(xiàn)蟲(chóng)在不同秸稈還田處理間均無(wú)顯著差異。SGK321的植食性線(xiàn)蟲(chóng)豐度在2.0 g·kg-1還田處理時(shí)顯著高于石遠(yuǎn)321（P＜0.05）；石遠(yuǎn)321土壤中植食性線(xiàn)蟲(chóng)豐度在不還田時(shí)顯著高于3.0 g·kg-1和4.0 g·kg-1還田處理（P＜0.05）；SGK321土壤中植食性線(xiàn)蟲(chóng)豐度在秸稈還田后各處理間均存在顯著差異（P＜0.05）。石遠(yuǎn)321中雜食/捕食性線(xiàn)蟲(chóng)豐度均高于SGK321，但無(wú)顯著差異，兩種棉花不同秸稈還田處理間也均無(wú)顯著差異。

花鈴期，除2.0 g·kg-1還田處理外，石遠(yuǎn)321的食細(xì)菌性線(xiàn)蟲(chóng)豐度均高于SGK321，且在4.0 g·kg-1還田處理時(shí)存在顯著差異（P＜0.05）；石遠(yuǎn)321土壤中食細(xì)菌性線(xiàn)蟲(chóng)豐度在4.0 g·kg-1還田處理時(shí)顯著高于2.0 g·kg-1和不還田處理（P＜0.05）；SGK321土壤中食細(xì)菌性線(xiàn)蟲(chóng)豐度在3.0 g·kg-1還田處理時(shí)顯著高于1.0 g·kg-1和不還田處理（P＜0.05）。石遠(yuǎn)321土壤中食真菌性線(xiàn)蟲(chóng)豐度均高于SGK321，在3.0 g·kg-1還田處理時(shí)存在顯著差異（P＜0.05），且在3.0 g·kg-1還田處理時(shí)顯著高于其他還田處理（P＜0.05）；SGK321土壤中食真菌性線(xiàn)蟲(chóng)豐度在各處理間無(wú)顯著差異。除2.0 g·kg-1還田處理，SGK321土壤中植食性線(xiàn)蟲(chóng)豐度均高于石遠(yuǎn)321，且在3.0、4.0 g·kg-1還田處理時(shí)存在顯著差異（P＜0.05）；石遠(yuǎn)321土壤中植食性線(xiàn)蟲(chóng)豐度在2.0 g·kg-1還田處理時(shí)顯著高于4.0 g·kg-1還田處理（P＜0.05）；SGK321土壤中植食性線(xiàn)蟲(chóng)豐度在1.0 g·kg-1還田處理時(shí)顯著高于其他處理（P＜0.05）。石遠(yuǎn)321土

壤中雜食/捕食性線(xiàn)蟲(chóng)豐度在1.0、3.0 g·kg-1還田處理時(shí)低于SGK321，但差異不顯著；石遠(yuǎn)321土壤中雜食/捕食性線(xiàn)蟲(chóng)豐度在2.0 g·kg-1還田處理時(shí)顯著高于3.0 g·kg-1還田處理（P＜0.05）；SGK321土壤中雜食/捕食性線(xiàn)蟲(chóng)豐度在不還田處理時(shí)顯著高于3.0、4.0 g·kg-1還田處理（P＜0.05）。

表1 蕾期不同還田處理土壤線(xiàn)蟲(chóng)的優(yōu)勢(shì)度Table 1 The dominance of soil nematodes under different treatments at budding stage

表2 花鈴期不同還田處理土壤線(xiàn)蟲(chóng)的優(yōu)勢(shì)度Table 2 The dominance of soil nematodes under different treatments at flowering and boll-forming stage

圖2 蕾期和花鈴期不同秸稈還田處理下各營(yíng)養(yǎng)類(lèi)群豐度Figure 2 Abundance of different nutrient groups in different returning amounts at budding stage and flowering and boll-forming stage

2.4 土壤線(xiàn)蟲(chóng)群落生態(tài)指數(shù)

SGK321的SR指數(shù)在蕾期不還田處理時(shí)顯著高于石遠(yuǎn)321（P＜0.05）（表3）。石遠(yuǎn)321的J指數(shù)在蕾期3.0 g·kg-1還田處理和花鈴期1.0 g·kg-1還田處理時(shí)顯著高于SGK321（P＜0.05）。H′指數(shù)兩品種間不存在顯著性差異。石遠(yuǎn)321不同還田處理中的SR、H′指數(shù)在蕾期4.0 g·kg-1還田處理時(shí)顯著高于1.0 g·kg-1還田處理（P＜0.05），J指數(shù)在各處理間均無(wú)顯著差異。SGK321不同還田處理中的H′指數(shù)在蕾期3.0 g·kg-1還田處理時(shí)顯著高于2.0 g·kg-1還田處理（P＜0.05），SR、J指數(shù)在各處理間均無(wú)顯著差異。

MI指數(shù)[25]是根據(jù)線(xiàn)蟲(chóng)的群落結(jié)構(gòu)特征評(píng)價(jià)人類(lèi)活動(dòng)對(duì)線(xiàn)蟲(chóng)群落結(jié)構(gòu)的影響，進(jìn)而可以表明土壤環(huán)境所受的干擾程度。石遠(yuǎn)321的MI指數(shù)在蕾期2.0 g·kg-1還田處理時(shí)顯著高于SGK321，在花鈴期1.0、4.0 g·kg-1還田處理時(shí)顯著高于 SGK321（P＜0.05）。蕾期不同還田處理中SGK321土壤線(xiàn)蟲(chóng)的MI指數(shù)在1.0、3.0 g·kg-1還田處理時(shí)顯著高于不還田、2.0 g·kg-1還田處理（P＜0.05）；花鈴期石遠(yuǎn)321不同還田處理的MI指數(shù)在4.0 g·kg-1還田處理時(shí)顯著高于不還田、1.0、2.0 g·kg-1還田處理（P＜0.05）。線(xiàn)蟲(chóng)的營(yíng)養(yǎng)類(lèi)群結(jié)構(gòu)與土壤碎屑食物網(wǎng)密切相關(guān)，常用食細(xì)菌和食真菌線(xiàn)蟲(chóng)的比值（即NCR）探測(cè)土壤有機(jī)質(zhì)的分解途徑[32]。SGK321的NCR指數(shù)在花鈴期3.0 g·kg-1還田處理時(shí)顯著高于石遠(yuǎn)321（P＜0.05）。不同還田處理中石遠(yuǎn)321的NCR指數(shù)在蕾期2.0 g·kg-1還田處理時(shí)顯著高于不還田處理（P＜0.05）；在花鈴期2.0 g·kg-1還田處理時(shí)顯著高于3.0 g·kg-1還田處理（P＜0.05）。不同還田處理中SGK321土壤線(xiàn)蟲(chóng)的NCR指數(shù)在蕾期和花鈴期各秸稈還田處理間均無(wú)顯著差異。WI指數(shù)用于反映土壤線(xiàn)蟲(chóng)群落組成和土壤的健康程度[21]，石遠(yuǎn)321的WI指數(shù)在蕾期2.0 g·kg-1還田處理和花鈴期1.0、4.0 g·kg-1還田處理時(shí)顯著高于 SGK321（P＜0.05）。不同還田處理中石遠(yuǎn)321的WI指數(shù)在花鈴期4.0 g·kg-1還田處理時(shí)顯著高于0、1.0、2.0 g·kg-1還田處理（P＜0.05）。SGK321的WI指數(shù)在蕾期3.0 g·kg-1還田處理時(shí)顯著高于不還田和2.0 g·kg-1還田處理（P＜0.05）。

花鈴期相比于蕾期，石遠(yuǎn)321中除SR、NCR指數(shù)外，其他生態(tài)指數(shù)均升高；SGK321中除H′、J指數(shù)外，其他生態(tài)指數(shù)均降低（表3）。

2.5 棉花品種、秸稈還田量和生育期對(duì)線(xiàn)蟲(chóng)數(shù)量及生態(tài)指標(biāo)的交互效應(yīng)

由表4可知，棉花生育期對(duì)植食性線(xiàn)蟲(chóng)和雜食/捕食性線(xiàn)蟲(chóng)數(shù)量有極顯著影響；棉花品種對(duì)線(xiàn)蟲(chóng)總量和植食性線(xiàn)蟲(chóng)數(shù)量有極顯著影響；棉花秸稈還田量對(duì)線(xiàn)蟲(chóng)總量、食細(xì)菌性線(xiàn)蟲(chóng)及食真菌性線(xiàn)蟲(chóng)數(shù)量有極顯著影響，對(duì)雜食/捕食性線(xiàn)蟲(chóng)數(shù)量有顯著影響。多因素方差分析表明：棉花品種和生育期對(duì)食真菌性和雜食/捕食性線(xiàn)蟲(chóng)數(shù)量有顯著交互效應(yīng)；棉花品種和秸稈還田量對(duì)線(xiàn)蟲(chóng)總量和植食性線(xiàn)蟲(chóng)數(shù)量有極顯著交互效應(yīng)，對(duì)食細(xì)菌性線(xiàn)蟲(chóng)有顯著交互效應(yīng)。

由表5可知，除NCR外，棉花生育期對(duì)SR、H′、J、MI和WI指數(shù)均有極顯著影響；棉花品種對(duì)J指數(shù)有極顯著影響，對(duì)WI指數(shù)有顯著影響；棉花秸稈還田量對(duì)MI、WI指數(shù)有極顯著影響，對(duì)SR、NCR指數(shù)有顯著影響。多因素方差分析表明：棉花品種和生育期對(duì)SR、MI指數(shù)有極顯著交互效應(yīng)，對(duì)H′、NCR、WI指數(shù)有顯著交互效應(yīng)；棉花品種和秸稈還田量對(duì)WI指數(shù)有顯著交互效應(yīng)；棉花品種、生育期和秸稈還田量對(duì)WI指數(shù)有顯著交互效應(yīng)。

3 討論

土壤線(xiàn)蟲(chóng)因其獨(dú)特的生物特性及其在土壤群落結(jié)構(gòu)中的重要作用，而成為指示土壤生態(tài)系統(tǒng)結(jié)構(gòu)和功能變化的有力工具[33]。陳群英等[34]連續(xù)兩年的轉(zhuǎn)Bt水稻田間種植未對(duì)土壤線(xiàn)蟲(chóng)總密度及各營(yíng)養(yǎng)類(lèi)群的線(xiàn)蟲(chóng)密度產(chǎn)生顯著影響。本研究表明，轉(zhuǎn)雙價(jià)基因棉和常規(guī)棉土壤線(xiàn)蟲(chóng)密度在不還田時(shí)無(wú)顯著差異，這與Yang等[16]的結(jié)果一致，說(shuō)明轉(zhuǎn)雙價(jià)基因棉花的種植對(duì)土壤線(xiàn)蟲(chóng)密度沒(méi)有顯著影響。秸稈還田后，兩種棉花的土壤線(xiàn)蟲(chóng)密度隨著秸稈還田量的增加而增加，這說(shuō)明秸稈的施入，為線(xiàn)蟲(chóng)提供了更加充足的營(yíng)養(yǎng)物質(zhì)和能量來(lái)源，有利于土壤線(xiàn)蟲(chóng)的生長(zhǎng)繁殖[21]。

表3 蕾期和花鈴期不同還田處理土壤線(xiàn)蟲(chóng)的生態(tài)指數(shù)Table 3 Ecological indicators of the soil nematode of the different treatments at budding stage and flowering and boll-forming stage

表4 棉花品種、秸稈還田量和生育期對(duì)線(xiàn)蟲(chóng)數(shù)量的交互效應(yīng)Table 4 Interactive effects of cotton varieties，straw returning amount and growth period on the number of nematode

表5 棉花品種、秸稈還田量和生育期對(duì)線(xiàn)蟲(chóng)生態(tài)指標(biāo)的交互效應(yīng)Table 5 Interactive effects of cotton varieties，straw returning amount and growth period on nematode ecological indicators

本試驗(yàn)中兩個(gè)棉花品種，5個(gè)還田處理共鑒定出線(xiàn)蟲(chóng)19科38屬。秸稈還田對(duì)螺旋屬和擬麗突屬等優(yōu)勢(shì)屬影響不大，但增加線(xiàn)蟲(chóng)稀有屬的數(shù)量。說(shuō)明秸稈還田主要是通過(guò)影響稀有類(lèi)群的數(shù)量進(jìn)而影響線(xiàn)蟲(chóng)群落結(jié)構(gòu)[35]。與不還田處理相比，秸稈還田處理中出現(xiàn)了三等齒屬、麗突屬等屬，表明秸稈還田增加了土壤線(xiàn)蟲(chóng)的豐富度[21]，提高了線(xiàn)蟲(chóng)多樣性[36]。

土壤線(xiàn)蟲(chóng)群落的營(yíng)養(yǎng)類(lèi)群可以有效地反映土壤線(xiàn)蟲(chóng)食物網(wǎng)營(yíng)養(yǎng)級(jí)關(guān)系及能量流動(dòng)途徑[29]。戚琳等[37]發(fā)現(xiàn)轉(zhuǎn)Bt水稻的田間種植增加植食性線(xiàn)蟲(chóng)豐度。但對(duì)轉(zhuǎn)雙價(jià)基因棉的研究未發(fā)現(xiàn)其改變線(xiàn)蟲(chóng)的營(yíng)養(yǎng)類(lèi)群[16]。本研究表明，轉(zhuǎn)雙價(jià)基因棉的種植也沒(méi)有改變土壤線(xiàn)蟲(chóng)的營(yíng)養(yǎng)類(lèi)群結(jié)構(gòu)，各處理間線(xiàn)蟲(chóng)群落組成不存在顯著差異，植食性和食細(xì)菌性線(xiàn)蟲(chóng)是優(yōu)勢(shì)營(yíng)養(yǎng)類(lèi)群。食真菌性線(xiàn)蟲(chóng)在兩個(gè)棉花品種所有處理中均屬于最低的營(yíng)養(yǎng)類(lèi)群，且秸稈還田后豐度低于不還田，這可能是因?yàn)槊藁ń斩捘举|(zhì)化很強(qiáng)，細(xì)菌比較難以利用，真菌在分解棉花秸稈中的作用更加突出。秸稈還田后，常規(guī)棉雜食/捕食性線(xiàn)蟲(chóng)豐度增加，這說(shuō)明秸稈的施用為雜食/捕食性線(xiàn)蟲(chóng)提供食物來(lái)源，增加了土壤中原生動(dòng)物等的密度[38]。

利用線(xiàn)蟲(chóng)的生態(tài)指數(shù)可以很好地反映不同秸稈還田處理間土壤線(xiàn)蟲(chóng)多樣性和群落結(jié)構(gòu)的差異程度[21]。本研究表明，蕾期石遠(yuǎn)321的SR、H′指數(shù)在4.0 g·kg-1與1.0、2.0 g·kg-1還田處理之間存在顯著差異，SGK321的H′指數(shù)在2.0 g·kg-1與3.0 g·kg-1還田處理之間存在顯著差異，這與Zhang等[39]的研究結(jié)果相同，說(shuō)明秸稈還田后，多樣性有一定的增加趨勢(shì)，在一定還田范圍內(nèi)能夠提高土壤線(xiàn)蟲(chóng)豐富度。但花鈴期兩種棉花品種之間土壤線(xiàn)蟲(chóng)H′、SR指數(shù)均沒(méi)有顯著差異，表明作物生育期對(duì)土壤線(xiàn)蟲(chóng)的影響大于轉(zhuǎn)雙價(jià)基因棉[34]。MI指數(shù)越低，表明土壤環(huán)境受到干擾程度越大[21]。本研究中，蕾期SGK321的MI指數(shù)在3.0 g·kg-1還田處理與不還田、2.0 g·kg-1還田處理之間均存在顯著差異，在花鈴期兩種棉花的MI指數(shù)在3.0、4.0 g·kg-1還田處理顯著高于不還田處理，說(shuō)明較高的秸稈還田量對(duì)土壤環(huán)境的影響相對(duì)更小。蕾期轉(zhuǎn)雙價(jià)基因棉土壤線(xiàn)蟲(chóng)的MI指數(shù)普遍高于常規(guī)棉，說(shuō)明轉(zhuǎn)雙價(jià)基因棉種植對(duì)土壤線(xiàn)蟲(chóng)的干擾較小，這與楊志國(guó)等[40]的研究結(jié)果相同。NCR指數(shù)值為0，代表土壤有機(jī)質(zhì)分解是完全依靠真菌分解途徑；若NCR值為1，則表示完全依靠細(xì)菌分解途徑[29]。WI指數(shù)值越大，表明土壤食微線(xiàn)蟲(chóng)豐度越高，土壤健康程度越高。本研究中，NCR值在0.68～0.90間，均大于0.5，表明各秸稈還田處理的盆栽土壤主要依靠細(xì)菌分解途徑分解有機(jī)物，且NCR在秸稈還田條件下顯著高于不還田處理，這表明秸稈還田能促進(jìn)食細(xì)菌性線(xiàn)蟲(chóng)繁殖并抑制植食性線(xiàn)蟲(chóng)的生長(zhǎng)。除蕾期0、2.0 g·kg-1還田處理和花鈴期1.0 g·kg-1還田處理外，WI值均大于1，表明盆栽所用土壤健康狀況普遍良好。

多因素方差分析表明，生育期對(duì)植食性線(xiàn)蟲(chóng)和雜食/捕食性線(xiàn)蟲(chóng)數(shù)量有極顯著影響，說(shuō)明生育期會(huì)引起土壤線(xiàn)蟲(chóng)密度的變化，這與李修強(qiáng)等[41]結(jié)論一致，對(duì)線(xiàn)蟲(chóng)SR、H′、J、MI和WI指數(shù)有極顯著影響，說(shuō)明生育期對(duì)線(xiàn)蟲(chóng)群落結(jié)構(gòu)組成有影響。品種對(duì)線(xiàn)蟲(chóng)總量和植食性線(xiàn)蟲(chóng)數(shù)量有極顯著影響，對(duì)J、WI指數(shù)有影響，可能是由于轉(zhuǎn)Bt作物外源基因的插入導(dǎo)致根系分泌物成分發(fā)生變化，進(jìn)而影響線(xiàn)蟲(chóng)群落的數(shù)量和組成[8]。秸稈還田量對(duì)線(xiàn)蟲(chóng)數(shù)量和生態(tài)指標(biāo)均有一定影響，說(shuō)明秸稈還田增加線(xiàn)蟲(chóng)豐度，豐富線(xiàn)蟲(chóng)種類(lèi)，有利于土壤生態(tài)穩(wěn)定。生育期與秸稈還田量和生育期、品種與秸稈還田量對(duì)線(xiàn)蟲(chóng)群落沒(méi)有交互作用，說(shuō)明對(duì)線(xiàn)蟲(chóng)數(shù)量的影響是相互抵消的。

轉(zhuǎn)基因作物秸稈還田引起的土壤生態(tài)系統(tǒng)的變化是一個(gè)緩慢的、長(zhǎng)期的過(guò)程，一年的盆栽試驗(yàn)無(wú)法得到全面的反映，需要開(kāi)展長(zhǎng)期的試驗(yàn)研究，以期更全面地評(píng)價(jià)轉(zhuǎn)基因作物秸稈還田對(duì)土壤生態(tài)環(huán)境的影響。

4 結(jié)論

（1）轉(zhuǎn)雙價(jià)基因棉的種植并未改變土壤線(xiàn)蟲(chóng)的營(yíng)養(yǎng)類(lèi)群組成，兩種棉花不同還田處理下的多樣性指數(shù)（H′）、線(xiàn)蟲(chóng)通路指數(shù)（NCR）均無(wú)顯著差異。

（2）秸稈不還田時(shí)，轉(zhuǎn)雙價(jià)基因棉和常規(guī)棉土壤線(xiàn)蟲(chóng)密度無(wú)顯著差異。秸稈還田后隨著秸稈還田量的增加，蕾期轉(zhuǎn)雙價(jià)基因棉土壤線(xiàn)蟲(chóng)數(shù)量先升高后降低，常規(guī)棉土壤線(xiàn)蟲(chóng)數(shù)量逐漸增加?；ㄢ徠趦煞N棉花土壤線(xiàn)蟲(chóng)數(shù)量逐漸升高。

（3）多因素方差分析表明，生育期、品種和秸稈還田量對(duì)線(xiàn)蟲(chóng)數(shù)量和生態(tài)指數(shù)均有一定的影響，生育期與品種和品種與秸稈還田量對(duì)線(xiàn)蟲(chóng)數(shù)量和生態(tài)指數(shù)有一定的交互效應(yīng)，而生育期與秸稈還田量和生育期、品種及秸稈還田量三者之間沒(méi)有交互作用。