王秀琴，歐后丁，于曉飛，芶劍渝，劉健鋒，申修賢，劉明宏，楊茂發(fā),*

（1.貴州大學(xué)昆蟲(chóng)研究所，貴州山地農(nóng)業(yè)病蟲(chóng)害重點(diǎn)實(shí)驗(yàn)室，貴陽(yáng) 550025；2.貴州大學(xué)煙草學(xué)院，貴陽(yáng) 550025；3.貴州省煙草公司遵義市公司，貴州 遵義 563000）

煙蚜Myzus persicae（Sulzer）又名桃蚜，屬半翅目Hemiptera 蚜科Aphididae 瘤蚜屬M(fèi)yzus，為世界性害蟲(chóng)，在我國(guó)各煙區(qū)均有分布。煙蚜是為害煙草的一種主要害蟲(chóng)，具有繁殖速度快、世代周期短、分布廣[1]和遷飛、擴(kuò)散能力強(qiáng)等特點(diǎn)[2]，為典型的r-型害蟲(chóng)[3]。通過(guò)刺吸煙株韌皮部汁液影響煙株生長(zhǎng)發(fā)育，受害嚴(yán)重的煙葉不僅產(chǎn)量降低，品質(zhì)也受到嚴(yán)重影響[3]。此外，煙蚜是黃瓜花葉病毒（Cucumber mosaic virus，CMV）、馬鈴薯Y 病毒（Potato virus Y，PVY）、煙草線(xiàn)條病毒（Tobacco streak virus，TSV）等100 多種病毒的傳播媒介[2,4]，排泄的蜜露還能誘發(fā)煤污病，使煙葉減產(chǎn)乃至絕收。利用生物防治控制害蟲(chóng)種群是被大力提倡的環(huán)境友好型措施。煙蚜的天敵昆蟲(chóng)主要有寄生性天敵和捕食性天敵兩類(lèi)，如煙蚜繭蜂[1-2]，七星瓢蟲(chóng)[5]、異色瓢蟲(chóng)[5-6]、南方小花蝽[7]、食蚜癭蚊[8-10]等。近年來(lái)，利用多種天敵協(xié)同防控害蟲(chóng)的策略已成為研究熱點(diǎn)。如在煙蚜低密度條件下，煙蚜繭蜂能快速建立種群，降低初始蚜量；煙蚜種群增大時(shí)，引入異色瓢蟲(chóng)能與煙蚜繭蜂協(xié)同降低煙蚜種群[6]。因此，捕食性天敵與寄生性天敵的協(xié)同增效作用是提高害蟲(chóng)防控效果的有效途徑之一，其中將食蚜癭蚊開(kāi)發(fā)作為煙蚜生物防治的有效補(bǔ)充具有重要意義。

食蚜癭蚊Aphidoletes aphidimyza(Rondani)隸屬雙翅目Diptera 癭蚊科Cecidomyiidae，具有分布廣，適應(yīng)性、繁殖力、搜尋能力和分散能力較強(qiáng)等特點(diǎn)[11-12]。捕食蚜蟲(chóng)種類(lèi)多，能捕食85 種蚜蟲(chóng)[13]，因而在生物防治中具有廣闊的開(kāi)發(fā)利用前景。蚜蟲(chóng)密度較大時(shí)，1 頭食蚜癭蚊幼蟲(chóng)可以捕食幾十頭蚜蟲(chóng)[14-15]。食蚜癭蚊不僅直接取食蚜蟲(chóng)體液，還具有獵殺作用，且獵殺的蚜蟲(chóng)量常高于捕食的蚜蟲(chóng)量[16]。因此，食蚜癭蚊是一種重要的捕食性天敵昆蟲(chóng)[13]，常用于對(duì)青椒[17]、甜椒[18-20]、蘿卜[21]等溫室作物上煙蚜的防控，也有食蚜癭蚊對(duì)煙田煙蚜防控的報(bào)道[9]。

目前國(guó)內(nèi)外對(duì)食蚜癭蚊的生物學(xué)特性[22-24]、繁殖與儲(chǔ)存[25-29]等內(nèi)容進(jìn)行了大量的研究，同時(shí)對(duì)食蚜癭蚊功能反應(yīng)的報(bào)道也逐漸增多[15,24,30-32]，但是關(guān)于食蚜癭蚊不同齡期幼蟲(chóng)對(duì)不同齡期煙蚜的捕食功能反應(yīng)還尚未報(bào)道。因此，本試驗(yàn)系統(tǒng)測(cè)定了食蚜癭蚊幼蟲(chóng)對(duì)煙蚜的捕食作用和種內(nèi)干擾效應(yīng)，從而明確不同齡期食蚜癭蚊幼蟲(chóng)對(duì)不同齡期煙蚜的控害潛能，以期為田間利用食蚜癭蚊防控?zé)熝撂峁├碚撘罁?jù)。

1 材料與方法

1.1 供試材料

煙蚜：2017 年5 月采集于貴州省畢節(jié)市林泉鎮(zhèn)煙草基地，在人工氣候室以煙草（K326）連續(xù)飼養(yǎng)多代后供試。

食蚜癭蚊：2017 年5 月采于貴州省遵義市樂(lè)山鎮(zhèn)煙草基地，以煙蚜連續(xù)飼養(yǎng)多代后供試。

1.2 飼養(yǎng)條件

于寧波江南儀器廠(chǎng)RXZ-380A-LED 人工氣候箱內(nèi)進(jìn)行飼養(yǎng)。溫度（25±1）℃，相對(duì)濕度（70±5）%，光周期為14 h 光照、10 h 黑暗。

1.3 試驗(yàn)方法

1.3.1 食蚜癭蚊對(duì)煙蚜的捕食功能反應(yīng) 1～2 齡、3齡和4 齡煙蚜均設(shè)10、20、30、40、50 和60 頭6個(gè)密度，將不同密度的煙蚜分別與1 頭饑餓24 h 處理的食蚜癭蚊2 齡或3 齡幼蟲(chóng)組合，放入直徑為35 mm 的培養(yǎng)皿中，并用帕拉膜封口，以昆蟲(chóng)針在帕拉膜表面扎孔保持通風(fēng)透氣，放入人工氣候箱內(nèi)繼續(xù)飼養(yǎng)，24 h 后于體視鏡下檢查，記錄被捕食的煙蚜。每個(gè)處理重復(fù)10 次。

1.3.2 食蚜癭蚊3 齡幼蟲(chóng)捕食煙蚜的種內(nèi)干擾反應(yīng)3 種煙蚜蟲(chóng)齡（1～2 齡、3 齡和4 齡）與饑餓24 h的食蚜癭蚊3 齡幼蟲(chóng)分別以20∶1、40∶2、60∶3、80∶4 和100∶5 比例處理放入直徑35 mm 的培養(yǎng)皿中（處理方法同1.3.1），24 h 后在體視鏡下檢查，記錄各處理的被捕食量。共15 個(gè)處理，每處理重復(fù)8 次。

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

1.4.1 捕食功能反應(yīng) 捕食功能反應(yīng)的數(shù)據(jù)分析通過(guò)兩個(gè)步驟實(shí)現(xiàn)。采用Logistc 回歸法對(duì)功能反應(yīng)的原始數(shù)據(jù)進(jìn)行分析并確定類(lèi)型[33]，利用多項(xiàng)式函數(shù)（1）對(duì)獵物被捕食的比例與獵物的初始密度進(jìn)行回歸分析：

式中：Na為被捕食的獵物，N0為獵物的初始密度，參數(shù)P0、P1、P2和P3采用最大似然法估計(jì)（SAS v9.4）。根據(jù)參數(shù)正負(fù)判定功能反應(yīng)曲線(xiàn)類(lèi)型：P1=0，為Ⅰ型功能反應(yīng)；P1＜0，為Ⅱ型功能反應(yīng)；P1＞0，P2＜0 則為Ⅲ型功能反應(yīng)[33-34]。由于本試驗(yàn)中沒(méi)有補(bǔ)充獵物，即獵物密度隨著時(shí)間推移而減小。因此，本試驗(yàn)采用（2）式的隨機(jī)捕食模型[35]進(jìn)行擬合：

式中：Na為被捕食的獵物，N0為獵物的初始密度，a為瞬時(shí)攻擊率（搜尋效率），Th為捕食者處理1 頭獵物所用的時(shí)間，T為獵物暴露于捕食者的時(shí)間（24 h）。參數(shù)估計(jì)采用PROCNLIN程序（SAS v9.4）。

1.4.2 種內(nèi)干擾反應(yīng) 根據(jù)（3）Hassell 模型[36-37]模擬食蚜癭蚊受到的種內(nèi)干擾系數(shù)：

式中E為平均捕食率，Q 為搜尋常數(shù)，P為捕食者密度，m為干擾系數(shù)，Na 為捕食的獵物總數(shù)，N為獵物密度。方程的擬合度用卡方（χ2）檢驗(yàn)判別。

對(duì)不同密度下食蚜癭蚊的日均捕食量進(jìn)行正態(tài)分布檢驗(yàn)，服從正態(tài)分布后進(jìn)行單因素方差（One-Way ANOVA）分析，方差齊性的數(shù)據(jù)用LSD法做多重比較，方差非齊性的數(shù)據(jù)用Tamhane 法做多重比較（SPSS 20.0）。

數(shù)據(jù)采用Excel 2016 進(jìn)行整理，SPSS 20.0 及SAS v9.4 軟件分析，SigmaPlot 14.0 軟件作圖。

2 結(jié)果

2.1 食蚜癭蚊對(duì)煙蚜的捕食功能反應(yīng)

利用Logistic 模型進(jìn)行回歸的試驗(yàn)結(jié)果見(jiàn)表1和表2。當(dāng)2 齡食蚜癭蚊幼蟲(chóng)捕食1～2 齡煙蚜?xí)r，無(wú)論是用二次還是三次回歸模型擬合的P1值均小于0，但不顯著（表1）。因此，需結(jié)合這兩個(gè)模型中的似然比進(jìn)一步判斷功能反應(yīng)類(lèi)型，因二次回歸模型的=60.23，p=0.36；三次回歸模型的=59.18，p=0.36，故判斷2 齡食蚜癭蚊捕食1～2齡煙蚜的功能反應(yīng)為類(lèi)型Ⅱ。當(dāng)捕食3 齡或4 齡煙蚜?xí)r，二次回歸模型擬合的P1值均小于0 且顯著。因此，判斷該齡期的食蚜癭蚊對(duì)3 齡或4 齡煙蚜的功能反應(yīng)類(lèi)型為Ⅱ型；而三次回歸模型擬合的P1值雖然均小于0 且系數(shù)不顯著，但因三次回歸模型對(duì)3 齡或4 齡煙蚜擬合的似然比分別是=14.27，p=1.00；=20.07，p=1.00，所以2 齡食蚜癭蚊對(duì)3 齡或4 齡煙蚜的功能反應(yīng)類(lèi)型也為Ⅱ型。同理，3齡食蚜癭蚊對(duì)1～2 齡、3 齡和4 齡煙蚜的功能反應(yīng)類(lèi)型為Ⅱ型（表2）。

研究結(jié)果表明，食蚜癭蚊2 齡和3 齡幼蟲(chóng)對(duì)各齡期煙蚜的捕食量均隨煙蚜密度的增加而增大。煙蚜密度增加到一定程度時(shí)，食蚜癭蚊捕食量的增加趨于緩慢，即食蚜癭蚊對(duì)煙蚜的捕食功能反應(yīng)呈負(fù)加速曲線(xiàn)，捕食量和獵物密度之間的關(guān)系為逆密度制約，這再一次證明了食蚜癭蚊2 齡或3 齡幼蟲(chóng)的捕食功能反應(yīng)類(lèi)型為Ⅱ型（圖1 和圖2）。

表1 2 齡食蚜癭蚊幼蟲(chóng)對(duì)煙蚜的捕食率與初始密度的logistic 回歸結(jié)果Table 1 Maximum likelihood estimates via logistic regression of proportion of aphids consumed by second instar larvae Aphidoletes aphidimyza and initial aphid densities

表2 3 齡食蚜癭蚊幼蟲(chóng)對(duì)煙蚜的捕食率與初始密度的logistic 回歸結(jié)果Table 2 Maximum likelihood estimates via logistic regression of proportion of aphids consumed by 3rd instar larvae Aphidoletes aphidimyza and initial aphid densities

圖1 食蚜癭蚊2 齡幼蟲(chóng)對(duì)不同齡期煙蚜的捕食功能反應(yīng)Fig.1 Functional response of second instar larvae Aphidoletes aphidimyza under different aphid densities (Mean±SE)

圖2 食蚜癭蚊3 齡幼蟲(chóng)對(duì)不同齡期煙蚜的捕食功能反應(yīng)Fig.2 Functional response of third instar larvae Aphidoletes aphidimyza under different aphid densities (Mean±SE)

Rogers’s 隨機(jī)捕食方程對(duì)2 齡食蚜癭蚊捕食不同齡期煙蚜的參數(shù)估計(jì)見(jiàn)表3。2 齡食蚜癭蚊對(duì)煙蚜的平均瞬時(shí)攻擊率（a）為0.019 9±0.003 7，對(duì)1～2 齡煙蚜最大（0.027 0），是4 齡煙蚜（0.014 1）的1.91 倍。2 齡食蚜癭蚊對(duì)煙蚜的平均處理時(shí)間為(7.415 5±2.690 9) h，1～2 齡煙蚜的是2.052 5 h，是3 齡煙蚜的0.211 倍，是4 齡煙蚜的0.196 倍。2 齡食蚜癭蚊對(duì)1～2 齡煙蚜的理論最大捕食量為11.69頭，而對(duì)4 齡煙蚜的理論最大捕食量?jī)H為2.29 頭。因此，2 齡食蚜癭蚊對(duì)1～2 齡煙蚜的處理時(shí)間短，理論最 大捕食量最大，由此說(shuō)明其對(duì)低齡煙蚜的控制作用最佳。

Rogers’s 隨機(jī)捕食方程對(duì)3 齡食蚜癭蚊捕食不同齡期煙蚜的參數(shù)估計(jì)見(jiàn)表4。3 齡食蚜癭蚊對(duì)不同齡期煙蚜的瞬時(shí)攻擊率（a）為1～2 齡>3 齡>4 齡，處理時(shí)間Th為3 齡>4 齡>1～2 齡，理論最大捕食量T/Th為1～2 齡>3 齡>4 齡。因此，3 齡食蚜癭蚊對(duì)低齡煙蚜（1～2 齡）的控制作用也最強(qiáng)。

2.2 食蚜癭蚊3 齡幼蟲(chóng)捕食煙蚜的種內(nèi)干擾反應(yīng)

由表5 可見(jiàn)，食蚜癭蚊3 齡幼蟲(chóng)對(duì)不同齡期煙蚜的日均捕食量隨著食蚜癭蚊密度的增加而顯著增大（1～2 齡：F4,35=72.89，p<0.001；3 齡：F4,35=334.96，p<0.001；4 齡：F4,35=145.07，p<0.001），但平均捕食率隨著食蚜癭蚊密度的增加而逐漸減小。食蚜癭蚊捕食1～2齡煙蚜的干擾反應(yīng)方程：E=0.355 8P-0.763，χ2=1.262，df=4，p=0.868＞0.05；捕食3 齡煙蚜的干擾反應(yīng)方程：E=0.2184P-0.5736，χ2=0.807，df=4，p=0.089＞0.05；捕食4 齡煙蚜的干擾反應(yīng)方程：E=0.195 7P-0.6415，χ2=5.80，df=4，p=0.215＞0.05；經(jīng)卡方（χ2）適合性檢驗(yàn)，3 個(gè)方程的理論值與觀(guān)測(cè)值均無(wú)顯著差異，表明3 個(gè)方程均可描述種內(nèi)干擾對(duì)食蚜癭蚊捕食煙蚜的影響。因此，食蚜癭蚊3 齡幼蟲(chóng)捕食煙蚜的種內(nèi)干擾系數(shù)平均值為（0.659 4±0.055 4），且捕食不同齡期煙蚜的種內(nèi)干擾作用不同，順序?yàn)?～2 齡（0.76）＞4 齡（0.64）＞3 齡（0.57）。

表3 2 齡食蚜癭蚊幼蟲(chóng)對(duì)煙蚜的捕食功能反應(yīng)參數(shù)Table 3 Estimates of instantaneous attack rate(a) and handling time (Th) of second instar larvae Aphidoletes aphidimyza

表4 3 齡食蚜癭蚊幼蟲(chóng)對(duì)煙蚜的捕食功能反應(yīng)參數(shù)Table 4 Estimates of instantaneous attack rate(a) and handling time (Th) of third instar larvae Aphidoletes aphidimyza

3 討 論

本研究結(jié)果表明，食蚜癭蚊對(duì)煙蚜各齡期的捕食量均隨著獵物密度的增加而增大，2 齡幼蟲(chóng)對(duì)低齡（1～2 齡）煙蚜的捕食量最大，3 齡幼蟲(chóng)是控制高齡（3～4 齡）煙蚜的最佳蟲(chóng)期，且捕食量明顯多于2 齡。這與林清彩等[24]報(bào)道的食蚜癭蚊幼蟲(chóng)隨著齡期的增長(zhǎng)其捕食量也明顯增加的結(jié)果基本一致。因此，在煙蚜發(fā)生初期，釋放2 齡幼蟲(chóng)不僅能控制煙蚜的初始密度，而且利于自身發(fā)育至3 齡，最終將煙蚜控制在一個(gè)低水平的范圍內(nèi)。而食蚜癭蚊在田間的釋放方式多為蛹[9]。因此，提前釋放食蚜癭蚊蛹可以更好地控制煙蚜。

當(dāng)2 齡食蚜癭蚊幼蟲(chóng)捕食3 齡或4 齡煙蚜?xí)r，二次Logistic回歸模型擬合的P1值均小于0且顯著，而三次Logistic 回歸模型擬合的P1值雖然均小于0，但不顯著。由此推測(cè)在使用Logistic 回歸模型擬合功能反應(yīng)曲線(xiàn)類(lèi)型參數(shù)時(shí)，若使用三次Logistic 回歸模型擬合的參數(shù)不顯著，可降元為二次Logistic回歸模型，再對(duì)參數(shù)進(jìn)行重新擬合；若參數(shù)仍不顯著，需再結(jié)合模型中的似然比最終綜合判斷功能反應(yīng)的類(lèi)型[33-34]。結(jié)果表明，食蚜癭蚊各齡幼蟲(chóng)（2齡或3 齡）對(duì)不同齡期煙蚜的捕食功能反應(yīng)均屬于類(lèi)型Ⅱ。這與食蚜癭蚊對(duì)其他蚜蟲(chóng)如甘藍(lán)蚜Brevicoryne brassicae、薔薇長(zhǎng)管蚜Macrosiphum rosae[30]、玉米蚜Rhopalosiphum maidis[15,24]、棉蚜Aphis gossypii[32]的捕食功能反應(yīng)類(lèi)型一致，表明功能反應(yīng)Ⅱ型是分析食蚜癭蚊捕食蚜蟲(chóng)功能反應(yīng)的最佳方法，可用類(lèi)型Ⅱ擬合該天敵對(duì)蚜蟲(chóng)的功能反應(yīng)。每頭食蚜癭蚊幼蟲(chóng)一生取食40～60 頭甘藍(lán)蚜Bravicoryne brassicae，28 頭蘋(píng)果黃蚜（又名繡線(xiàn)菊蚜）Aphis citricola，25 頭蠶豆蚜Aphis fabae，28頭桃蚜Myzus persicae或13 頭豌豆蚜Acyrthosiphon pisum[16]，而本試驗(yàn)結(jié)果發(fā)現(xiàn)3 齡幼蟲(chóng)對(duì)1～2 齡、3齡和4 齡煙蚜的最大捕食量分別為12、9 和10 頭，原因可能與蚜蟲(chóng)種類(lèi)、大小和密度相關(guān)，甚至與食蚜癭蚊的生理狀態(tài)有關(guān)。

在一定的空間范圍內(nèi)，天敵昆蟲(chóng)常常對(duì)鄰近的同種個(gè)體存在著明顯的干擾反應(yīng)，使其平均捕食率下降[38]。食蚜癭蚊3 齡幼蟲(chóng)對(duì)不同齡期煙蚜的捕食也受到較強(qiáng)的種內(nèi)干擾，在煙蚜密度和捕食空間一定的情況下，食蚜癭蚊3 齡幼蟲(chóng)的平均捕食率隨其自身密度增加而降低。這與異色瓢蟲(chóng)捕食枸杞木虱會(huì)受到自身密度的干擾[39]和七星瓢蟲(chóng)捕食草地貪夜蛾時(shí)表現(xiàn)出種內(nèi)干擾使捕食率呈現(xiàn)下降的趨勢(shì)[40]一致。同樣在寄生性天敵中也存在種內(nèi)干擾作用，如沙地木虱跳小蜂[41]、枸杞木虱嚙小蜂[42]寄生枸杞木虱時(shí)均受到種內(nèi)干擾作用。因此，根據(jù)益害比合理地釋放食蚜癭蚊，有助于對(duì)煙蚜的最佳防控。

由于本研究在恒溫恒濕的條件下進(jìn)行，且只對(duì)食蚜癭蚊不同齡期的幼蟲(chóng)捕食量進(jìn)行了相關(guān)試驗(yàn)，而當(dāng)食蚜癭蚊處于自然環(huán)境條件時(shí)，對(duì)煙蚜的捕食量受到日照、風(fēng)、溫度、雨水等因素的影響。因而在田間釋放食蚜癭蚊，需綜合考慮自身的捕食效率、田間環(huán)境條件、獵物密度等因素，最終確定適當(dāng)?shù)尼尫帕?、釋放次?shù)和釋放方式。

4 結(jié) 論

研究結(jié)果表明，食蚜癭蚊不同齡期的幼蟲(chóng)對(duì)煙蚜各齡期的捕食量均隨著煙蚜密度的增加而增大，其捕食功能反應(yīng)均屬于Ⅱ型，3 齡食蚜癭蚊幼蟲(chóng)對(duì)煙蚜的控制效能最高，試驗(yàn)建立的模型能較好地反應(yīng)不同齡期食蚜癭蚊對(duì)各齡煙蚜的潛在捕食能力，對(duì)實(shí)踐上利用食蚜癭蚊防治煙蚜具有一定的指導(dǎo)意義。