張夢蝶 錢路兵 澤桑梓 楊 斌 李宗波
(1.西南林業(yè)大學(xué)生物多樣性保護(hù)學(xué)院 云南省森林災(zāi)害預(yù)警與控制重點實驗室 昆明 650224; 2.云南省林業(yè)和草原有害生物防治檢疫局 昆明 650224)
昆蟲表皮碳?xì)浠衔?Cuticular hydrocarbons, CHCs)是指廣泛存在于昆蟲上表皮中一類由長鏈烷烴和烯烴所組成的混合物,其性質(zhì)穩(wěn)定,加之成分和含量在種間種內(nèi)的差異性,即使是復(fù)合種、近緣種、隱形種、同種個體的不同發(fā)育階段或社會品級等也有明顯差異,對于昆蟲化學(xué)分類、物種演化、信息通訊與調(diào)控、胚后發(fā)育、農(nóng)林害蟲治理及其益蟲利用等研究均具有重要的啟示作用(Pageetal., 1997; Blomquistetal., 2010; Yoonetal., 2012; Derstineetal., 2018; Bersonetal., 2019; 李群臣等, 2019)。云南切梢小蠹(Tomicusyunnanensis)隸屬于鞘翅目象蟲總科(Curculionoidea)小蠹科(Scolytinae)切梢小蠹屬(Tomicus),又稱為縱坑切梢小蠹或云南松縱坑切梢小蠹,是一種為害云南松的蛀干害蟲,自20世紀(jì)80年代暴發(fā)危害以來,累積致死的活立木蓄積量超過20萬hm2,嚴(yán)重影響了云南松林的可持續(xù)發(fā)展和生態(tài)安全(李麗莎 等, 1997; Kirkendalletal., 2008; 沈紹偉等, 2018)。當(dāng)前,雖然已有一套關(guān)于云南切梢小蠹的綜合治理措施,例如噴灑殺蟲劑和生物制劑、釋放天敵昆蟲、信息素和寄主氣味或餌木誘殺、非寄主氣味干擾、營造混交林等(李浩然等, 2000; 王海林等, 2002; 李雙成等, 2003; 葉輝等, 2006; 岳鋒等, 2011; 王大偉等, 2015; Wuetal., 2019),但仍無法控制該蟲擴(kuò)散蔓延的趨勢,這說明在云南切梢小蠹治理體系中存在著缺失環(huán)節(jié),特別是關(guān)于云南切梢小蠹的發(fā)育生物學(xué)特征。但對于蛀干、潛葉、土壤中等營隱蔽式生活和危害的昆蟲,其發(fā)育生物學(xué)難以直接觀察,特別是幼蟲的齡數(shù)不易確定; 而通過人工室內(nèi)飼養(yǎng)導(dǎo)致幼蟲生活環(huán)境的改變,常常引起齡數(shù)出現(xiàn)較大變化,從而無法精確了解寄主幼蟲發(fā)育齡數(shù)及發(fā)育進(jìn)度。由于這種基礎(chǔ)研究的缺失,會極大地降低害蟲科學(xué)預(yù)測預(yù)報的準(zhǔn)確性,隨之降低有害生物的綜合治理效能(陳輝, 2002; Castexetal., 2018),如遷飛高峰措置(Peak shift displacement)使得噴灑的真菌制劑殺蟲效果不佳,或者釋放的天敵昆蟲無法及時找到合適寄主(王海林等, 2002; 葉輝等, 2006)。同時,隨著全球化進(jìn)程的加快,國際交往日益頻繁,并考慮到全球氣候變化的大背景(Castexetal., 2018; Liuetal., 2019),該蟲的擴(kuò)散及危害潛能會被放大,產(chǎn)生的生態(tài)破壞性將難以估量。本研究旨在通過林間大量抽樣調(diào)查,運(yùn)用形態(tài)數(shù)據(jù)測定自然條件下云南切梢小蠹幼蟲的齡數(shù),結(jié)合不同齡級幼蟲的表皮碳?xì)浠衔镞M(jìn)行定性定量分析,為深入了解云南切梢小蠹生物學(xué)特征及其科學(xué)治理該蟲危害提供依據(jù)。
云南切梢小蠹幼蟲采集于曲靖市沾益區(qū)九龍山林場(103°7′15″ E,25°0′35″ N),為該區(qū)域內(nèi)林木害蟲的優(yōu)勢種群,主要在樹干中上部蛀干繁殖(沈紹偉等, 2018)。于2017年11月至2019年5月,分批次帶回新蛀入或已在樹干內(nèi)發(fā)育的云南切梢小蠹幼蟲。為了更準(zhǔn)確地區(qū)分齡數(shù),將位于坑道兩側(cè)的卵(60粒)置于培養(yǎng)皿中的濕潤濾紙上,室溫條件下保濕培養(yǎng),待其孵化后作為1齡幼蟲(圖1A),而老齡幼蟲以蛹室形成且有木屑封堵為判斷依據(jù)(圖1B)。將每次采集的云南切梢小蠹幼蟲分別放入75%的乙醇溶液中保存,共收集到各齡幼蟲1 426頭。參考象蟲總科幼蟲齡數(shù)的劃分標(biāo)準(zhǔn)(Mizelletal., 1979; Castaeda-Vildzolaetal., 2016),用蔡氏體式鏡(V20,German)對各個幼蟲進(jìn)行拍照,用Image J(Version 1.62,https:∥imagej.nih.gov)測量其頭殼寬和體長,并用掃描電鏡(Zeiss sigma 300, German)觀察不同齡數(shù)幼蟲的微觀形態(tài),測其上顎長與寬。
圖1 云南切梢小蠹初孵幼蟲和末齡幼蟲的判定Fig. 1 Determination of larval instars in Tomicus yunnanensis using the freshly hatched larvae from eggs and last larvae blocked in the pupal cellsA1: 坑道兩側(cè)的卵 Eggs on the both sides of gallery ; A2: 初孵幼蟲 The freshly hatched larvae from eggs; B1: 封堵在蛹室(虛線)內(nèi)的老齡幼蟲 The last instar larvae are blocked by the wood chips (the arrow) in the pupal cell (the dashed line); B2: 末齡幼蟲 Morphology of the last instar larvae.
采用浸提法對云南切梢小蠹不同齡數(shù)幼蟲的表皮碳?xì)浠衔镞M(jìn)行提取(Pageetal., 1997; Blomquistetal., 2010)。首先將從木段中解剖出來的幼蟲立即置入-20 ℃冰箱內(nèi)殺死,超聲波清洗后(冰浴),依據(jù)頭殼寬確定齡數(shù),每齡重復(fù)15次,每次測定60只,用400 μL正己烷洗滌2次,第1次10 min(靜置5 min,震蕩1 min,再靜置3.5 min,震蕩30 s),第2次靜置1 min,然后將混合液放置在2 mL棕色試劑瓶內(nèi),-20 ℃保存待用。在GC-MS(Agilent 7890A-5975C,USA)分析以前,先將樣品用氮氣吹干,滴入少量正己烷,接著加入內(nèi)標(biāo)化合物正十四烷(n-C14)和正十六烷(n-C16)各2 μL(2 ng·μL-1),定容至50 μL。接著取1 μL定溶液從進(jìn)樣口注入,質(zhì)譜柱型號為極性HP-5 MS(30 m × 0.25 mm × 0.25 μm)。升溫程序: 起始溫度為40 ℃,以10 ℃·min-1升溫至150 ℃,以5 ℃·min-1升溫至300 ℃,保持10 min。高純度氮氣作為運(yùn)載體,流速為1 mL·min-1,進(jìn)樣口溫度250 ℃,柱前壓100 kPa,進(jìn)樣采用不分流模式。MS條件: EI源,70 eV,傳輸線溫度250 ℃,離子源溫度230 ℃,質(zhì)量掃描范圍全掃描。采用NIST 08譜庫檢索法、峰面積標(biāo)準(zhǔn)化法(Aitchison, 1986)、科瓦特指數(shù)法(Marquesetal., 2000)及其與C8—C40烷烴混合標(biāo)準(zhǔn)品(500 μg·mL-1,AccuStandard,USA)比對進(jìn)行鑒定(計算公式如下),化合物含量根據(jù)內(nèi)標(biāo)的量計算。
Zip= ln[Aip/g(Ap)]。
式中:Zip為轉(zhuǎn)換后待鑒定化合物的峰面積;Aip為待鑒定化合物的峰面積;g(Ap)為各組分色譜峰面積的幾何均值。
KI=100n+100 ×[tR(x)-tR(n)]/[tR(n+1)-tR(n)]。
式中:KI為科瓦特指數(shù);n為待鑒定化合物前面正構(gòu)烷烴的碳原子數(shù);x為待測化合物碳原子數(shù);tR(x)為保留時間;tR(n)和tR(n+1)分別為待鑒定化合物色譜峰前一個和后一個正構(gòu)烷烴的保留時間。
所有的數(shù)據(jù)統(tǒng)計分析均在開源軟件R(Version 3.5.1,https:∥cran.r-project.org/)中進(jìn)行,制圖用GraphPad Prism(Version 8.0.2)完成。對測定的頭殼寬度,首先進(jìn)行頻數(shù)分布分析,間隔最小設(shè)置為0.02 mm,接著再用體長、上顎長和寬等進(jìn)行驗證,以此確定云南切梢小蠹幼蟲齡數(shù)。齡數(shù)劃分后,進(jìn)一步根據(jù)測定指標(biāo)的均值、標(biāo)準(zhǔn)誤差、變幅、變異系數(shù)、Brooks指數(shù)、Crosby指數(shù)等確定最佳齡數(shù),若變異系數(shù)小于20%、Crosby指數(shù)小于10%,則表明分齡指標(biāo)合理(Castaeda-Vildzolaetal., 2016),接著用回歸分析擬合齡數(shù)與體長、頭殼寬間的關(guān)系,用Tukey氏多重比較法對比各齡幼蟲指標(biāo)的差異性。對于表皮碳?xì)浠衔?,利用FactoMineR軟件包進(jìn)行主成分分析(Lêetal., 2008),選用極大方差法旋轉(zhuǎn)因子,并基于相關(guān)系數(shù)矩陣計算主成分得分系數(shù),提取特征值大于1的成分,最后再用RandomForest軟件包進(jìn)行化合物貢獻(xiàn)量分析(Mccormicketal., 2014),提取基尼指數(shù)(MeanDecreaseGini)大于2的成分; 用Kruskal-Wallis test進(jìn)行各齡數(shù)幼蟲化合物總量和關(guān)鍵化合物量的多重比較。
幼蟲頭殼寬度頻次分布具有3個界線明顯的峰區(qū),而每個峰區(qū)中又有明顯的峰值,分別位于0.50、0.825和1.04 mm處(圖2),因幼蟲發(fā)育過程中需解除外骨骼限制而產(chǎn)生蛻皮,每次蛻皮幼蟲即會突長一次,進(jìn)入下一齡,所以從頻次圖上可以推測云南切梢小蠹幼蟲共有3齡。進(jìn)一步來看,幼蟲頭殼寬度與齡數(shù)間的變異幅度是比較小的(7%~14.6%),且Crosby指數(shù)小于0.1,并經(jīng)統(tǒng)計檢驗和回歸分析均達(dá)到了極顯著水平(P<0.000 1; 表1),表明將云南切梢小蠹幼蟲齡數(shù)劃分為3齡是合理準(zhǔn)確的。掃描電鏡觀察顯示云南切梢小蠹幼蟲共有12體節(jié),其表面微觀形態(tài),包括前額心形盾及其表面剛毛數(shù)量和排列形式、額須上的味覺感器、上顎及其齒數(shù)、內(nèi)唇等均沒有不同,僅長度與寬度有所增加(圖3; 表1)。當(dāng)用體長、上顎長和寬進(jìn)行齡數(shù)擬合時,盡管這些指標(biāo)均可將云南切梢小蠹幼蟲劃分為3齡,統(tǒng)計檢驗和回歸分析也顯示有顯著差異(體長:F=3 383.6,P<0.000 1; 上顎長:F=740.82,P<0.000 1; 上顎寬:F=91.395,P<0.000 1; df=1),但幼蟲體長的變異系數(shù)超過20%,2齡和3齡幼蟲的上顎寬度間差異不顯著(表1)。
圖2 云南切梢小蠹幼蟲頭殼寬度頻次分布Fig. 2 Frequency distribution of the head-capsule widths of the larval T. yunnanensis
表1 云南切梢小蠹幼蟲各項測量指標(biāo)的分齡與統(tǒng)計①Tab.1 Measurement and statistical analysis of the four various indexes for the larval instars of T. yunnanensis
圖3 云南切梢小蠹不同齡數(shù)幼蟲的微觀形態(tài)Fig. 3 Micrograph of larval instars of T. yunnanensisA—C: 1~3齡幼蟲 The instar numbers of the larvae from the first instar to the third instar; 字母后的數(shù)字1~3分別示整體、頭殼背面和腹面 The numerals behind the capital letter indicate the ventral view of the larval body, the dorsal and ventral view of larval head capsule, respectively.
共從云南切梢小蠹幼蟲中鑒定出19種表皮碳?xì)浠衔?,主要為正鏈烷烴、支鏈烷烴和烯烴3類烴類物質(zhì),含量分別為每蟲4.41 ng± 0.67 ng、0.26 ng ±0.04 ng和1.22 ng± 0.18 ng,各烴類物質(zhì)在不同齡級幼蟲中分布相差較大(表2)。其中,正鏈烷烴以正二十五烷和正二十七烷為代表,占比分別為55.02% ± 8.20%和21.37% ± 3.19%; 支鏈烷烴中以13-甲基-二十九烷含量最高,占整個支鏈烷烴的94.04%±14.02%; 烯烴類以1, 21-二十二烷二烯為主,占比為63.03%±0.39%。
表2 云南切梢小蠹不同齡數(shù)幼蟲的表皮碳?xì)浠衔锓N類、含量及其主成分分析與隨機(jī)森林判定①Tab.2 Mean relative content of cuticular hydrocarbon and their results of principal component scores and random forest analysis on the larval instars of T. yunnanensis
云南切梢小蠹初孵幼蟲(1齡幼蟲)的表皮碳?xì)浠衔镉?個,碳數(shù)介于10~19之間,C10—C14和C15—C19各占近50%; 2齡幼蟲有8個化合物,碳數(shù)為10~29,其中C20—C24和C25—C29為主要成分,分別占19.98%和71.11%; 而老熟幼蟲(3齡幼蟲)則有17個,以C25—C29和C30—C35為主,分別占61.94%和12.74%(表2)。隨著齡級增加,幼蟲表皮碳?xì)浠衔锏暮糠謩e為每蟲2.32 ng±0.59 ng、44.40 ng±11.49 ng和280.75 ng±72.48 ng,齡級間有著顯著差異(χ2=39.13,P<0.000 1)(圖4)。主成分分析顯示第一主成分和第二主成分聯(lián)合貢獻(xiàn)率為81.7%,其中正十二烷(7.46%)、正十三烷(7.44%)、2,6,11-三甲基-十二烷(3.48%)、正十八烷(7.67%)等12種化合物超過第一主成分貢獻(xiàn)量,正二十三烷(13.32%)和正二十五烷(12.78)是第二主成分的重要貢獻(xiàn)者,不同齡數(shù)幼蟲的化合物分離程度高,齡數(shù)與齡數(shù)間無重疊區(qū)域(表2、 圖5)。在幼蟲表皮碳?xì)浠衔飿?gòu)成的重要性方面,隨機(jī)森林分析顯示2,6,11-三甲基-十二烷、正十五烷、正十九烷、正二十三烷和正二十五烷5種化合物具有較大的重要性(表2)。其中,2,6,11-三甲基-十二烷僅存在于1齡幼蟲,而正十五烷在3齡幼蟲中相對含量僅為0.1%,為一種微量化合物。在不同齡級幼蟲上,這5種化合物的含量上有著顯著差異(P<0.000 1),2,6,11-三甲基-十二烷為1齡幼蟲的標(biāo)識性物質(zhì),而正二十三烷和正二十五烷可作為2齡和3齡幼蟲的識別物(表2、圖6)。
圖5 云南切梢小蠹不同齡數(shù)幼蟲表皮碳?xì)浠衔锏闹鞒煞值梅諪ig. 5 Scores plot of principal components extracted from cuticular hydrocarbons of the larval instar of T. yunnanensis
圖6 云南切梢小蠹不同齡數(shù)幼蟲表皮中5種重要碳?xì)浠衔锏暮縁ig. 6 Mean contents of five important compounds of cuticular hydrocarbons on the larval instar of T. yunnanensis
昆蟲表皮碳?xì)浠衔锸抢ハx化學(xué)分類中采用時間最久的生物學(xué)指標(biāo)(Bioindicators)(Jackson, 1970; Blomquistetal., 2010; 李群臣等, 2019),如早在20世紀(jì)70年代,Jackson (1970)就發(fā)現(xiàn)碳數(shù)為23~26的正二十三烷、順-9-二十三烯、3-甲基-二十三烷、11-甲基-二十三烷、13-甲基-二十五烷等5種長鏈烷烴可輕松將3種大蠊(Periplaneta)區(qū)分開。隨著化學(xué)分析技術(shù)及其多元統(tǒng)計分析方法的大量應(yīng)用,發(fā)現(xiàn)表皮碳?xì)浠衔锍蓞^(qū)分蜚蠊目外,還可用來區(qū)分鞘翅目(Pageetal., 1997)、鱗翅目(Hay-Roeetal., 2007)、雙翅目(Barbosaetal., 2017)、半翅目(Yoonetal., 2012)、等翅目(Blomquistetal., 2010)、同翅目(Raboudietal., 2005)、膜翅目(Soaresetal., 2017)、直翅目(Brozaetal., 2000)等眾多昆蟲,并且一些表皮碳?xì)浠衔锓从沉似渑吆蟀l(fā)育狀態(tài)(Zhuetal., 2006; Blomquistetal., 2010; Kuoetal., 2012; Yoonetal., 2012)。例如,漸變態(tài)的點蜂緣蝽(Riptortuspedestris)卵的標(biāo)識物是3-甲基-二十七烷,若蟲是正二十九烷,成蟲是正二十七烷、正三十一烷、13,17-/15,19-二甲基-三十三烷(Yoonetal., 2012); 而全變態(tài)的紅顏金蠅(Chrysomyarufifacies)幼蟲期主要是正烷烴類,隨齡級增加短鏈碳?xì)浠衔锖砍手笖?shù)級減少,代表性物質(zhì)是正二十五烷和正二十九烷(Zhuetal., 2006)。不僅如此,形態(tài)極為相似的近似種鱗翅目夜蛾科的棉鈴蟲與煙青蟲2種幼蟲表皮碳?xì)浠衔锝M分也出現(xiàn)區(qū)別,并且材料為5~6齡幼蟲碳數(shù)分布24~32與云南切梢小蠹幼蟲3齡幼蟲主要碳數(shù)分布25~35極為相似(高明媛等, 1999)。從云南切梢小蠹幼蟲表皮碳?xì)浠衔锏慕Y(jié)果來看,不同齡級幼蟲的化合物組成差別較大,特別是對分類有重要意義的2,6,11-三甲基-十二烷、正十五烷、正十九烷、正二十三烷和正二十五烷這5種化合物,良好的分離效果可能反映了其在齡級中的動態(tài)變化,也可能與總含量的指數(shù)級增加有關(guān)(y=e0.086 5x-0.629,R=0.911)(圖4)。盡管一些化合物含量極低,但其在化合物整體組成中不可或缺,顯示了一些新的成分分析方法,如隨機(jī)森林、非度量多維尺度分析(non-metric multidimensional scaling,NMDS)、典型辨別分析(canonical discriminant analysis,CDA)等,對于探索多元數(shù)據(jù)背后所隱含生物學(xué)含義具有重要的作用(Mccormicketal., 2014)。另一方面,云南切梢小蠹幼蟲在樹干內(nèi)發(fā)育,且每一個幼蟲均有各自獨立的坑道進(jìn)行取食,考慮到樹皮提供的保護(hù)及相對高含水量的環(huán)境,防止體內(nèi)水分流失應(yīng)不是營隱蔽生活昆蟲所面臨的問題,更多應(yīng)該與其生理活動和代謝水平有關(guān)(Blomquistetal., 2010),正如1齡幼蟲的烷烴類化合物僅有8個,且碳數(shù)均小于20,而隨著齡級的增加,老熟幼蟲有17種,而碳數(shù)也達(dá)到最大的35個,其發(fā)育進(jìn)度中化合物種類與碳數(shù)增加的結(jié)果與棲境相似昆蟲在發(fā)育進(jìn)度的變化相符(Zhuetal., 2006; Kuoetal., 2012; Yoonetal., 2012),但與裸露生活的點蜂緣蝽、葉甲、胡蜂、齒小蠹成蟲等差別較大(Pageetal., 1997; Yoonetal., 2012; Zhangetal., 2014; Soaresetal., 2017)。另外,昆蟲表皮碳?xì)浠衔镦滈L每增加一個仲碳原子(-CH2-),其熔點會提升2 ℃,進(jìn)而會增加幼蟲體壁的堅硬程度,對溫度的容忍能力亦會提升(Menzeletal., 2017),這也很好地解釋了隨著昆蟲齡數(shù)的增加,包括云南切梢小蠹在內(nèi)的眾多昆蟲的表皮碳?xì)浠衔镦滈L為何增加的現(xiàn)象。
用幼蟲表皮碳?xì)浠衔锟奢p松地將云南切梢小蠹3個齡級幼蟲區(qū)分開,比用頭殼寬度、體長、上顎長與寬的分離效果更好(表1、圖5)。昆蟲的生長周期表現(xiàn)出明顯的周期性,特別是鱗翅目(Godinetal., 2002)、鞘翅目(Mizelletal., 1979; Castaeda-Vildzolaetal., 2016)、雙翅目(Richardietal., 2013)、膜翅目(Vilhelmsen, 2011)等幼蟲頭殼骨化程度高且相對穩(wěn)定的昆蟲,身體各部的長寬會按指數(shù)級增長,局部的生長狀況可以指示整個蟲體的生長狀況。根據(jù)云南切梢小蠹幼蟲的4個測定指標(biāo),雖然均可將其劃分為3齡,但因各齡幼蟲體長和上顎寬重疊度較大,加之顎的長與寬測定相對困難,以它們作為齡數(shù)劃分依據(jù),存在諸多的限制條件,因此,出頭殼寬度仍不失是一個良好的形態(tài)指標(biāo)。但是,林間大量抽樣過程中仍會出現(xiàn)明顯錯誤,特別是頭殼寬為0.60~0.675 mm和0.9~0.94 mm這些介于齡數(shù)之間的蟲體樣本,而用幼蟲表皮碳?xì)浠衔锬軌蛎黠@克服這些缺點,特別是用基尼指數(shù)大于2的5種化合物或者對齡期有重要標(biāo)識的2,6,11-三甲基-十二烷、正二十三烷和正二十五烷。未來,將幼蟲表皮碳?xì)浠衔锏难芯拷Y(jié)果與便攜式氣質(zhì)聯(lián)用儀、電子鼻、表面增強(qiáng)拉曼光譜等快速分析鑒定方法相結(jié)合(Castexetal., 2018; Turlingsetal., 2018),會實現(xiàn)云南切梢小蠹發(fā)育進(jìn)度的實時監(jiān)測,為提升預(yù)測預(yù)報精度和綜合治理該蟲危害奠定基礎(chǔ)。
云南切梢小蠹幼蟲表皮碳?xì)浠衔镉烧溚闊N、支鏈烷烴和烯烴3類19種烴類物質(zhì)組成,從1齡幼蟲至3齡幼蟲,化合物數(shù)量從8個增加到17個,碳數(shù)從10增加到35,含量從每蟲0.139 ng±0.007 ng增至16.845 ng±1.366 ng,呈指數(shù)級遞增。幼蟲表皮碳?xì)浠衔锟擅黠@將不同齡數(shù)的幼蟲區(qū)分開,齡數(shù)與齡數(shù)間無重疊,劃分標(biāo)準(zhǔn)顯著優(yōu)于頭殼寬、體長、上顎長與寬等形態(tài)指標(biāo)法。2,6,11-三甲基-十二烷、正十五烷、正十九烷、正二十三烷、正二十五烷具有較大的重要性,含量上亦有顯著差異,其中2,6,11-三甲基-十二烷可用于1齡幼蟲的標(biāo)識性物質(zhì),正二十三烷和正二十五烷可識別2齡和3齡幼蟲,這將為研究云南切梢小蠹幼蟲發(fā)育進(jìn)度及監(jiān)測預(yù)報提供科學(xué)依據(jù)。