路獻(xiàn)勇， 張 帥， 呂麗敏， 雒珺瑜， 王春義， 崔金杰

中國(guó)農(nóng)業(yè)科學(xué)院棉花研究所，棉花生物學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室，河南 安陽(yáng) 455000

自1990年國(guó)外首次用農(nóng)桿菌介導(dǎo)法獲得Cry1Ab、Cry1Ac轉(zhuǎn)基因棉花植株以來(lái)(Perlaketal.,1990)，編碼Cry1Ac蛋白的基因已成功應(yīng)用于植物基因工程，并培育出了轉(zhuǎn)基因抗蟲(chóng)棉用于大規(guī)模生產(chǎn)。中國(guó)研究人員也于1993年將Cry1Ac基因?qū)氲街袊?guó)主產(chǎn)棉區(qū)的數(shù)個(gè)栽培品種中，獲得了第1代單價(jià)抗蟲(chóng)棉(郭三堆等，1999)。目前中國(guó)主要以種植轉(zhuǎn)單價(jià)Bt基因的抗蟲(chóng)棉花為主，由于轉(zhuǎn)入的Bt基因類(lèi)型單一，長(zhǎng)期種植容易導(dǎo)致棉鈴蟲(chóng)HelicoverpaarmigeraHübner對(duì)Cry1Ac基因的抗性增加，因此，為了延緩靶標(biāo)害蟲(chóng)對(duì)Cry1Ac蛋白的抗性，培育疊加基因的抗蟲(chóng)棉花品種受到廣泛重視。目前，轉(zhuǎn)Cry1Ac+Cry2Ab基因棉已在美國(guó)、澳大利亞、墨西哥等國(guó)家商業(yè)化種植，我國(guó)也培育出了轉(zhuǎn)雙價(jià)基因Cry1Ac+Cry2Ab的抗蟲(chóng)棉(雒珺瑜等，2011)。為了揭示不同Bt蛋白的抗蟲(chóng)效果，本文研究了Cry1Ac、Cry1Ab、Cry2Ab、Cry1Ah、Cry2Aa和Cry1Ca 6種Bt蛋白對(duì)棉鈴蟲(chóng)幼蟲(chóng)室內(nèi)種群存活率、體重和生長(zhǎng)發(fā)育的影響，期望為新一代轉(zhuǎn)基因抗蟲(chóng)棉花的培育提供依據(jù)，為新型轉(zhuǎn)基因棉花材料的抗蟲(chóng)性鑒定和安全評(píng)價(jià)提供參考。

1 材料與方法

1.1 供試棉鈴蟲(chóng)

室內(nèi)敏感品系：2005年6月采自河南省安陽(yáng)縣白壁鎮(zhèn)中國(guó)農(nóng)業(yè)科學(xué)院棉花研究所棉田第2代，在室內(nèi)不用任何藥劑處理，連續(xù)用人工飼料飼養(yǎng)70代。

1.2 供試Bt蛋白

Cry1Ac、Cry1Ab、Cry2Ab、Cry1Ah、Cry2Aa和Cry1Ca 6種Bt蛋白均由中國(guó)農(nóng)業(yè)科學(xué)院植物保護(hù)研究所中?？妻r(nóng)公司(北京)提供。

1.3 棉鈴蟲(chóng)人工飼料

采用卓樂(lè)姒等(1981)的棉鈴蟲(chóng)人工飼料配方經(jīng)稍微改進(jìn)。人工飼料配方如下：黃豆粉200 g，玉米粉300 g，小麥粉300 g，抗壞血酸10 g，酵母粉80 g，甲醛10 mL，36%乙酸60 mL，苯甲酸鈉9 g，瓊脂35 g，水2000 mL。成蟲(chóng)期飼喂10%的蜂蜜水。養(yǎng)蟲(chóng)室每月用紫外燈消毒1次，棉鈴蟲(chóng)蛹用甲醛消毒15 min，以防止棉鈴蟲(chóng)遭受核角體病毒的侵染。

在配制飼料過(guò)程中，當(dāng)飼料內(nèi)的溫度降到55～60 ℃時(shí)，取不同濃度的Bt蛋白分別加到不同處理的人工飼料中，每個(gè)處理的飼料量為100 g；用攪棒迅速攪勻，待飼料凝固后放入4 ℃冰箱備用。使用前切成1 g左右的塊狀，分別放入24孔培養(yǎng)板中用于飼養(yǎng)棉鈴蟲(chóng)幼蟲(chóng)，培養(yǎng)板表面覆蓋白布以防止幼蟲(chóng)逃逸。

1.4 生物測(cè)定方法

1.4.1 不同Bt蛋白對(duì)棉鈴蟲(chóng)初孵幼蟲(chóng)的毒力測(cè)定 測(cè)定包括預(yù)備試驗(yàn)和正式試驗(yàn)2個(gè)部分。預(yù)備試驗(yàn)共設(shè)置5個(gè)質(zhì)量濃度梯度，將供試幼蟲(chóng)的死亡率控制在10%～90%范圍內(nèi)；其中Cry1Ab、Cry1Ac和Cry2Ab的濃度梯度為1.0、0.5、0.1、0.01和0.001 μg·g-1；Cry1Ah、Cry2Aa和Cry1Ca的質(zhì)量濃度梯度為20.0、5.0、2.0、1.0和0.5 μg·g-1。正式試驗(yàn)共設(shè)置7個(gè)濃度梯度，其中Cry1Ac、Cry1Ab和Cry2Ab的濃度梯度為1.0、0.5、0.2、0.1、0.05、0.01和0.001 μg·g-1；Cry1Ah、Cry2Aa和Cry1Ca的濃度梯度為20.0、10.0、5.0、2.5、1.25、0.625和0.3125 μg·g-1。均用5 mL的Na2CO3溶液作為對(duì)照。

毒力測(cè)定前先用pH 9.5的Na2CO3溶液對(duì)毒蛋白進(jìn)行稀釋?zhuān)缓髮⑾♂尯蟮腂t毒蛋白混入飼料中(Beegleetal.,1990)。測(cè)定時(shí)在飼料中接入孵化后1日齡幼蟲(chóng)(24 h±2 h)，在溫度(27±1) ℃、濕度(75±5)%、光周期(光照/黑暗)14/10以及合適的光強(qiáng)條件下飼養(yǎng)。每處理24頭幼蟲(chóng)，重復(fù)6次。飼養(yǎng)7 d后記錄幼蟲(chóng)存活數(shù)，計(jì)算校正死亡率。用SPSS 19.0軟件進(jìn)行Probit回歸分析，計(jì)算LC50值(賈春生，2006)。

1.4.2 不同Bt蛋白對(duì)棉鈴蟲(chóng)幼蟲(chóng)生長(zhǎng)發(fā)育的影響 配制飼料前用pH 9.5的Na2CO3溶液分別將6種蛋白稀釋成100 μg·mL-1，配制飼料時(shí)取稀釋好的Bt蛋白溶液5 mL混入100 g飼料，使其在飼料中的濃度為1.0 μg·g-1，對(duì)照僅混入5 mL的Na2CO3。

測(cè)定設(shè)初孵幼蟲(chóng)、1齡幼蟲(chóng)和2齡幼蟲(chóng)3組試驗(yàn)。每組試驗(yàn)設(shè)6種Bt蛋白和對(duì)照共7個(gè)處理，每處理32頭幼蟲(chóng)，重復(fù)3次。測(cè)定時(shí)在飼料中接入幼蟲(chóng)，放置在與1.4.1相同的條件下飼養(yǎng)。初孵幼蟲(chóng)和1齡幼蟲(chóng)飼養(yǎng)7 d后，分別記錄幼蟲(chóng)死亡率，并用電子天平(精確到0.00001 g)稱(chēng)取所有活蟲(chóng)體重，計(jì)算體重抑制率。

體重抑制率(%)=(對(duì)照組平均體重-處理后平均體重)/對(duì)照組平均體重×100。

2齡幼蟲(chóng)飼養(yǎng)3、5、7和9 d后分別稱(chēng)取幼蟲(chóng)體重，每天記錄幼蟲(chóng)發(fā)育天數(shù)、化蛹數(shù)、蛹重、蛹期天數(shù)和羽化數(shù)，直至羽化完成。試驗(yàn)數(shù)據(jù)用SAS 8.1軟件分析。

2 結(jié)果與分析

2.1 不同Bt蛋白對(duì)棉鈴蟲(chóng)初孵幼蟲(chóng)的毒力測(cè)定

由表1可知，不同Bt蛋白對(duì)棉鈴蟲(chóng)的毒力效應(yīng)不同，用Bt蛋白連續(xù)飼喂7 d后的致死中濃度(LC50)依次為Cry1Ab>Cry1Ac>Cry2Ab>Cry2Aa>Cry1Ah。

表1 不同Bt蛋白處理的棉鈴蟲(chóng)初孵幼蟲(chóng)的致死中濃度(LC50)和線(xiàn)性回歸方程參數(shù)Table 1 Mean lethal concentrations (LC50)and probit regression line parameters of Bt proteins against the newly hatched larvae of cotton bollworm (H.armigera)

2.2 不同Bt蛋白對(duì)棉鈴蟲(chóng)低齡幼蟲(chóng)校正死亡率和體重抑制率的影響

棉鈴蟲(chóng)初孵幼蟲(chóng)和1齡幼蟲(chóng)取食含不同Bt蛋白飼料7 d后，不同處理間的校正死亡率有極顯著差異(初孵幼蟲(chóng)，df=5，F(xiàn)=341.81，P<0.0001；1齡幼蟲(chóng)，df=5，F(xiàn)=329.02，P<0.0001)。不同處理間的體重抑制率也有極顯著差異(初孵幼蟲(chóng)，df=5，F(xiàn)=337.50，P<0.0001；1齡幼蟲(chóng)，df=5，F(xiàn)=57.84，P<0.0001)。Cry1Ac和Cry1Ab蛋白對(duì)棉鈴蟲(chóng)初孵幼蟲(chóng)的校正死亡率極顯著高于Cry2Ab蛋白，與Cry2Ab蛋白相比，Cry1Ah和Cry2Aa蛋白處理的棉鈴蟲(chóng)初孵幼蟲(chóng)的校正死亡率顯著降低，Cry1Ca處理下，初孵幼蟲(chóng)死亡率極顯著低于另外5種蛋白。而對(duì)于1齡棉鈴蟲(chóng)則略有不同，Cry1Ab蛋白顯著高于其余5種蛋白，Cry1Ac和Cry2Ab蛋白極顯著高于Cry1Ah、Cry2Aa和Cry1Ca，而Cry1Ah蛋白極顯著高于Cry2Aa，Cry1Ca蛋白極顯著低于Cry2Aa。Cry1Ac、Cry1Ab和Cry2Ab處理的初孵幼蟲(chóng)體重抑制率極顯著高于Cry1Ah和Cry2Aa蛋白處理，而Cry1Ca處理下的棉鈴蟲(chóng)體重抑制率則極顯著低于Cry1Ah和Cry2Aa蛋白；6種蛋白處理的1齡幼蟲(chóng)的體重抑制率的多重比較結(jié)果與初孵幼蟲(chóng)相同(圖1)。

2.3 不同Bt蛋白處理棉鈴蟲(chóng)2齡幼蟲(chóng)后的體重變化

棉鈴蟲(chóng)2齡幼蟲(chóng)取食含不同Bt蛋白的飼料后，不同處理間幼蟲(chóng)體重有極顯著性差異(1 d，df=6，F(xiàn)=175.29，P<0.0001；3 d，df=6，F(xiàn)=94.77，P<0.0001；7 d,df=6，F(xiàn)=66.42，P<0.0001；9 d，df=6，F(xiàn)=924.33，P<0.0001)。多重比較結(jié)果表明，對(duì)照和Cry1Ca處理之間棉鈴蟲(chóng)的體重在1、3、7和9 d均沒(méi)有顯著差異；Cry1Ac、Cry1Ab處理下的2齡幼蟲(chóng)僅在9 d后體重極顯著高于Cry2Ab蛋白；而Cry1Ah處理下的幼蟲(chóng)體重僅在7 d后顯著低于Cry2Aa，而在1、3和9 d后均是極顯著低于該蛋白。與Cry2Ab相比，Cry1Ah處理下的幼蟲(chóng)體重僅在1 d后極顯著增加，其他時(shí)間段均沒(méi)有顯著差異(表2)。

2.4 不同Bt蛋白處理棉鈴蟲(chóng)2齡幼蟲(chóng)后的生物學(xué)參數(shù)差異

棉鈴蟲(chóng)2齡幼蟲(chóng)取食含不同Bt蛋白的飼料后，不同處理間除了蛹期顯著差異外，其他生物學(xué)參數(shù)均為極顯著差異(11 d后幼蟲(chóng)死亡率，df=6，F(xiàn)=9.47，P=0.0006<0.001；幼蟲(chóng)發(fā)育歷期，df=4，F(xiàn)=183.54，P<0.0001；化蛹率，df=6，F(xiàn)=97.92，P<0.0001；蛹重，df=6，F(xiàn)=27.92，P<0.001；蛹期，df=4，F(xiàn)=6.41，P=0.0129<0.05；羽化率，df=6，F(xiàn)=107.37，P<0.0001)。Cry1Ac和Cry1Ab蛋白處理下11 d后棉鈴蟲(chóng)的死亡率顯著高于其他4種蛋白，差異不顯著。用Cry1Ac和Cry1Ab蛋白處理的棉鈴蟲(chóng)無(wú)法完成整個(gè)幼蟲(chóng)發(fā)育歷期，而用Cry1Ah、Cry2Aa、Cry2Ab、Cry1Ca和對(duì)照處理的棉鈴蟲(chóng)幼蟲(chóng)發(fā)育歷期均有顯著差異，幼蟲(chóng)發(fā)育歷期的順序?yàn)镃K

圖1 1.0 μg·g-1不同Cry蛋白對(duì)棉鈴蟲(chóng)初孵幼蟲(chóng)和1齡幼蟲(chóng)校正死亡率和體重抑制率的影響Fig.1 Abbott-corrected mortality rates and weight inhibition rates caused by active Cry proteins against the newly hatched and first-instar larvae of cotton bollworm (H.armigera) fed with 1.0 μg of toxin /g food圖中數(shù)據(jù)為平均數(shù)±標(biāo)準(zhǔn)誤；柱子上不同字母表示數(shù)據(jù)經(jīng)Duncan氏法多重比較差異顯著，大寫(xiě)字母表示差異極顯著(P<0.01),小寫(xiě)字母表示差異顯著(P<0.05)。The data in the figure are mean±SE, and bars with the different letters above are significantly different by Duncan′s multiple ranging test. Capital letters show highly extremely significant differences (P<0.01), and small letters, show significant differences (P<0.05).

處理Treatment棉鈴蟲(chóng)幼蟲(chóng)體重 Cotton bollworm larval body weight (mg)1 d3 d7 d9 dCry1Ac1.28±0.22dD1.50±0.51dC3.50±1.26dD8.92±1.61dDCry1Ab0.98±0.03dD1.89±0.14dC8.07±2.01dD15.41±3.63dDCry2Ab1.17±0.18dD3.02±0.96cdC30.45±3.30cC68.89±6.67cCCry1Ah1.86±0.17cC4.52±0.20cC33.23±3.26cBC67.35±11.70cCCry2Aa3.13±0.10bB12.45±1.09bB72.97±8.85bB158.87±0.85bBCry1Ca4.46±0.57aA18.87±2.78aA188.09±31.36aA331.63±9.14aACK3.62±0.19aA16.73±0.78aA175.97±32.38aA336.54±11.65aA

表中不同字母表示數(shù)據(jù)經(jīng)Duncan氏法多重比較差異顯著，大寫(xiě)字母表示差異極顯著(P<0.01),小寫(xiě)字母表示差異顯著(P<0.05)。Different letters represent significant differences among treatments when using a Duncan′s multiple ranging test. Capital letters show highly significant differences (P<0.01) and small letters, significant differences (P<0.05).

表3 1.0 μg ·g-1的6種Bt蛋白分別處理棉鈴蟲(chóng)2齡幼蟲(chóng)后生物學(xué)參數(shù)的顯著差異Table 3 Responses of cotton bollworm fed with 1.0 μg of the six types of Bt proteins /g food

表中不同字母表示數(shù)據(jù)經(jīng)Duncan氏法多重比較差異顯著，大寫(xiě)字母表示差異極顯著(P<0.01),小寫(xiě)字母表示差異顯著(P<0.05)。Different letters represent significant differences among treatments when using a Duncan′s multiple ranging test. Capital letters show highly significant differences (P<0.01) and small letters, significant differences (P<0.05).

3 討論

本文系統(tǒng)測(cè)定了Cry1Ac、Cry1Ab、Cry2Ab、Cry1Ah、Cry2Aa和Cry1Ca 6種蛋白對(duì)棉鈴蟲(chóng)幼蟲(chóng)的毒力，國(guó)內(nèi)外也有一些文獻(xiàn)報(bào)道了Cry蛋白對(duì)棉鈴蟲(chóng)毒力(韓嵐嵐等，2008；Leeetal.,1996; Chakrabartietal.,1998; Kaliaetal.,2005; Li & Bouwer,2012)。

在澳大利亞，已報(bào)道只有Cry1Ab、Cry1Ac、Cry2Aa和Cry2Ab蛋白能夠在可接受的劑量下殺死棉鈴蟲(chóng)，而Cry1Ca被認(rèn)為對(duì)棉鈴蟲(chóng)是無(wú)毒的(Liaoetal.,2002)。在印度，Cry1Ac對(duì)棉鈴蟲(chóng)幼蟲(chóng)的毒力最大，Cry1Ab和Cry1Aa次之(Kaliaetal.,2005)。在南非，Cry1Ac和Cry2Aa對(duì)棉鈴蟲(chóng)的毒力最大，Cry1Ab次之，Cry1Aa毒性較小(Li & Bouwer,2012)。在中國(guó)，Cry1Ac被認(rèn)為是對(duì)棉鈴蟲(chóng)幼蟲(chóng)毒力最強(qiáng)的Cry蛋白，Cry2Ab次之，此外，Cry1Ah對(duì)棉鈴蟲(chóng)幼蟲(chóng)也有一定的殺蟲(chóng)活性，但由于純化時(shí)易形成寡聚體導(dǎo)致其殺蟲(chóng)活性不很穩(wěn)定(韓嵐嵐等，2008; 周子珊等，2011)。而本試驗(yàn)中測(cè)得的LC50顯示，Cry1Ab和Cry1Ac對(duì)棉鈴蟲(chóng)幼蟲(chóng)的毒力最強(qiáng)，Cry2Ab次之，Cry1Ah和Cry2Aa對(duì)棉鈴蟲(chóng)幼蟲(chóng)毒性較低。由此可見(jiàn)，不同國(guó)家和地區(qū)的棉鈴蟲(chóng)種群對(duì)同一Bt毒素的敏感性存在差異，其原因可能是長(zhǎng)時(shí)間的地理隔離。

估算棉鈴蟲(chóng)田間種群抗性基因頻率時(shí)，一般采用Wuetal.(1999)確定的1.0 μg·mL-1飼料為診斷濃度。李海強(qiáng)等(2012)對(duì)新疆北部地區(qū)轉(zhuǎn)Bt基因棉外源蛋白的表達(dá)量進(jìn)行了測(cè)定，發(fā)現(xiàn)4個(gè)品種苗期葉片的Bt(Cry1Ac)蛋白含量平均為1.0 μg·g-1左右時(shí)，對(duì)當(dāng)?shù)孛掴徬x(chóng)有較好的控制作用。本試驗(yàn)以相同質(zhì)量濃度(1.0 μg·g-1)的6種Bt蛋白分別對(duì)棉鈴蟲(chóng)初孵幼蟲(chóng)、1齡幼蟲(chóng)和2齡幼蟲(chóng)進(jìn)行處理，發(fā)現(xiàn)Cry1Ab和Cry1Ac對(duì)3種發(fā)育狀態(tài)下的棉鈴蟲(chóng)幼蟲(chóng)毒性均是最強(qiáng)，用這2種蛋白對(duì)2齡幼蟲(chóng)連續(xù)處理，棉鈴蟲(chóng)羽化率為0，不能完成整個(gè)生長(zhǎng)發(fā)育。Cry2Ab對(duì)2齡前的幼蟲(chóng)控制效果較好，但對(duì)2齡幼蟲(chóng)控制較差；而Cry1Ah和Cry2Aa對(duì)2齡前的棉鈴蟲(chóng)毒力接近，LC50幾乎相近，但兩者對(duì)2齡幼蟲(chóng)控制效果卻相差很大，可能是由于Cry1Ah蛋白純化過(guò)程中形成了寡聚體，殺蟲(chóng)活性不穩(wěn)定引起的(周子珊等，2011)。

為了阻止或推遲靶標(biāo)害蟲(chóng)對(duì)Cry1Ac蛋白抗性的產(chǎn)生，人們提出了多基因策略，即將多個(gè)Bt基因同時(shí)轉(zhuǎn)入植物中用來(lái)控制農(nóng)田害蟲(chóng)，Cry1Ac和Cry2Aa、Cry1Ac和Cry2Ab可被同時(shí)轉(zhuǎn)入植物中控制棉鈴蟲(chóng)，因?yàn)檫@2個(gè)蛋白可使用不同的受體結(jié)合位點(diǎn)(Cacciaetal.,2010；Gahanetal.,2005)。對(duì)Cry1Ac產(chǎn)生抗性的棉鈴蟲(chóng)品系，對(duì)Cry2Aa和Cry2Ab是敏感的(Akhursttetal.,2003)。而對(duì)Cry2Ab產(chǎn)生抗性的棉鈴蟲(chóng)品系，則被認(rèn)為對(duì)Cry2Aa也有抗性，但對(duì)Cry1Ac敏感(Mahonetal.,2007)。盡管Cry2Ab、Cry2Aa和Cry1Ah對(duì)棉鈴蟲(chóng)幼蟲(chóng)的毒性不如Cry1Ab和Cry1Ac，但是仍然可以作為控制中國(guó)棉鈴蟲(chóng)幼蟲(chóng)的替代策略。

郭三堆, 崔洪志, 夏蘭芹. 1999. 雙價(jià)抗蟲(chóng)轉(zhuǎn)基因棉花研究. 中國(guó)農(nóng)業(yè)科學(xué), 32(3)： 3-9.

韓嵐嵐, 宋福平, 張杰. 2008. 蘇云金芽孢桿菌殺蟲(chóng)晶體蛋白對(duì)棉鈴蟲(chóng)活性分析. 東北農(nóng)業(yè)大學(xué)學(xué)報(bào), 39(8)： 21-24.

賈春生. 2006. 利用SPSS軟件計(jì)算殺蟲(chóng)劑的LC50. 昆蟲(chóng)知識(shí), 43(3)： 414-417.

李海強(qiáng), 李號(hào)賓, 王冬梅. 2012. 轉(zhuǎn)基因抗棉鈴蟲(chóng)抗草甘膦棉花對(duì)棉鈴蟲(chóng)的抗蟲(chóng)性及對(duì)二種地老虎幼蟲(chóng)生長(zhǎng)發(fā)育的影響. 應(yīng)用昆蟲(chóng)學(xué)報(bào), 49(4)： 882-888.

雒珺瑜, 崔金杰, 辛惠江. 2011. 轉(zhuǎn)Cry1Ac+Cry2Ab基因棉花對(duì)天敵生長(zhǎng)發(fā)育的影響. 中國(guó)棉花, 38(12)： 11-12.

周子珊, 束長(zhǎng)龍, 梁革梅, 蒼晶，宋福平，張杰. 2011. Cry1Ah蛋白活性片段在溶液中的單體和寡聚體殺蟲(chóng)活性分析. 植物保護(hù), 37(3)： 48-51.

卓樂(lè)姒, 黃月蘭, 楊家榮. 1981. 棉鈴蟲(chóng)人工飼料的研究. 昆蟲(chóng)學(xué)報(bào), 24(1)： 108-110.

Akhurstt R J, James W and Bird L J. 2003. Resistance to the Cry1Ac delta-endotoxin ofBacillusthuringiensisin the cotton bollworm,Helicoverpaarmigera(Lepidoptera : Noctuidae).JournalofEconomicEntomology, 96： 1290-1299.

Beegle C C， Hickle L A and Fitch W L. 1990. Bioassay methods for quantification ofBacillusthuringiensisdelta-endotoxins∥Hickle L A and Fitch W L.AnalyticalChemistryofBacillusthuringiensis. Washington, DC: American Chemical Society, 14-21.

CacciaA S, Hemandez-rodriguez C S and Mahon R J. 2010. Binding site alteration is responsible for field-isolated resistance toBacillusthuringiensisCry2A insecticidal proteins in twoHelicoverpaspecies.PLoSONE, 5： e9975.

Chakrabarti S K, Mandaokar A and Kumar P A. 1998. Efficacy of lepidopteran specific delta-endotoxins ofBacillusthuringiensisagainstHelicoverpaarmigera.JournalofInvertebratePathology, 72: 336-337.

Gahan L J, Ma Y T and Coble M L M. 2005. Genetic basis of resistance to Cry1Ac and Cry2Aa inHeliothisvirescens(Lepidoptera: Noctuidae).JournalofEconomicEntomology, 98： 1357-1368.

Lee M K, Curtiss A and Alcantara E. 1996. Synergistic effect of theBacillusthuringiensistoxins CryIAa and CryIAc on the gypsy moth,Lymantriadispar.AppliedandEnvironmentalMicrobiology, 62： 583-586.

Liao C Y, Heckel D G and Akhurst R. 2002. Toxicity ofBacillusthuringiensisinsecticidal proteins forHelicoverpaarmigeraandHelicoverpapunctigera(Lepidoptera : Noctuidae), major pests of cotton.JournalofInvertebratePathology, 8： 55-63.

Li H and Bouwer G. 2012. Toxicity ofBacillusthuringiensisCry proteins toHelicoverpaarmigera(Lepidoptera: Noctuidae) in South Africa.JournalofInvertebratePathology, 109：110-116.

Mahon R J, Olsen K M and Downes S. 2007. Frequency of alleles conferring resistance to the Bt toxins Cry1Ac and Cry2Ab in Australian populations ofHelicoverpaarmigera(Lepidoptera : noctuidae).JournalofEconomicEntomology, 100： 1844-1853.

Perlak F J, Deaton R W and Armstrong T A. 1990. Insect resistant cotton plants.Bio-Technology, 8： 939-943.

Wu K, Guo Y and Lv N. 1999. Geographic variation in susceptibility ofHelicoverpaarmigera(Lepidoptera: Noctuidae) toBacillusthuringiensisinsecticidal protein in China.JournalofEconomicEntomology, 92： 273-278.