田余波， 任艷*， 劉海波， 毛羽， 吳發(fā)明

(1.西南民族大學(xué)藥學(xué)院，民族藥物研究所，成都610041； 2. 遵義醫(yī)科大學(xué)藥學(xué)院，貴州遵義563003)

芫菁科Meloidae昆蟲分布廣泛，是豆類作物常見害蟲(李曉婷，2010；潘昭等，2011)，同時(shí)該科昆蟲多具有重要藥用價(jià)值，應(yīng)用歷史悠久(劉海波等，2019)，如大斑溝芫菁Hycleusphaleratus(Pallas, 1781)和眼斑溝芫菁Hycleuscichorii(Linnaeus, 1758)被《中華人民共和國(guó)藥典》收錄(任艷等，2018；國(guó)家藥典委員會(huì)，2020)，其他有相關(guān)藥用記載還包括綠芫菁Lyttacaraganae(Pallas, 1781)、鋸角豆芫菁Epicautagorhami(Marseul, 1873)等(王雪梅，2007；王騰蛟等，2018)。芫菁科昆蟲在新疆、甘肅、寧夏、河北、內(nèi)蒙古等多地廣泛分布(張潤(rùn)楊，2018；楊益春，2019)，而哈思山區(qū)是甘肅芫菁科昆蟲聚集的典型代表區(qū)域(賀春貴等，2005)。該類昆蟲分泌的斑蝥素及其衍生物具有抗腫瘤、增強(qiáng)免疫、殺蟲等作用，醫(yī)療價(jià)值極高(Yan，1993；曾瑤波等，2016；顏道宇等，2017；李濟(jì)森等，2019)。作為重要的藥用昆蟲資源，當(dāng)前市場(chǎng)需求量日益增長(zhǎng)，在人工養(yǎng)殖方面存在瓶頸，無(wú)法規(guī)?；?，臨床使用仍以野外收集為主，過度不合理采集使其野生數(shù)量日漸減少(陳仕江，2004；王雪梅等，2007；莫讓瑜，2014)，而目前有關(guān)芫菁科昆蟲分布特征的研究尚無(wú)詳盡報(bào)道。

昆蟲種群空間分布是昆蟲種群的重要生態(tài)學(xué)屬性之一(楊現(xiàn)明等，2020)，不僅對(duì)了解昆蟲種群的爆發(fā)、擴(kuò)散具有重要的意義，而且在抽樣技術(shù)及監(jiān)測(cè)、預(yù)報(bào)中也發(fā)揮了一定的指導(dǎo)作用(魯武鋒等，2018)。對(duì)芫菁科藥用昆蟲聚集典型代表區(qū)域進(jìn)行系統(tǒng)調(diào)查，探索其空間分布規(guī)律和生活習(xí)性，對(duì)該類藥用昆蟲資源的保護(hù)和綜合開發(fā)利用具有重要的作用和實(shí)際意義。因此，本文采用隨機(jī)取樣法對(duì)芫菁科藥用昆蟲的分布進(jìn)行調(diào)查，利用聚集度指標(biāo)法、Taylor冪函數(shù)法及Iwao回歸分析法對(duì)昆蟲的空間分布類型及抽樣技術(shù)進(jìn)行研究，以期為芫菁科藥用昆蟲資源的監(jiān)測(cè)、保護(hù)和綜合開發(fā)奠定基礎(chǔ)，也為防治該科昆蟲因數(shù)量過度造成農(nóng)業(yè)危害提供參考。

1 材料與方法

1.1 研究地概況

哈思山區(qū)是以哈思山自然保護(hù)區(qū)為中心的農(nóng)林牧交錯(cuò)帶，地處甘肅省白銀市靖遠(yuǎn)縣境內(nèi)，屬隴中黃土高原農(nóng)牧林區(qū)，地理坐標(biāo)104°13′40″～104°35′00″E，36°58′20″～37°2′40″N，海拔2 400～3 017 m。年均氣溫6～7 ℃，年均降水量350～410 mm，降水極不均勻，多集中于7—9月(楊萬(wàn)軍，2003；張啟立等，2019)。全縣現(xiàn)有耕地7.58萬(wàn)hm2，草原面積37.33萬(wàn)hm2，其中哈思山區(qū)草原主要為高寒草甸類、溫性草原類和溫性荒漠草原類(董四娟，2012)。該區(qū)域植被覆蓋較好，但森林植被類型少，同時(shí)由于地處旱地農(nóng)業(yè)向荒漠牧地過渡邊界，荒漠化較嚴(yán)重，主要有石質(zhì)荒漠、沙質(zhì)荒漠和鹽化荒漠(舒海元，1997；李連惠，2003)。

1.2 調(diào)查方法

本研究于2019年7月，當(dāng)?shù)剀据伎评ハx頻繁活動(dòng)時(shí)期，根據(jù)哈思山區(qū)概況及地形地貌特征，選擇農(nóng)田型、草地型、森林型和荒漠型4種具有代表性的生態(tài)類型樣地各3塊，每塊樣地設(shè)置5 m2(n=7)正方形樣方，樣方間距約為500 m。部分采樣生境如圖1 所示，記錄樣方內(nèi)芫菁科主要藥用昆蟲及其所分布植物的種類和數(shù)量。采用實(shí)驗(yàn)室自制昆蟲捕捉裝置進(jìn)行樣方內(nèi)目標(biāo)昆蟲的捕捉，進(jìn)行分類學(xué)鑒定。

芫菁科主要藥用昆蟲空間分布型測(cè)定取樣方法：選取14個(gè)樣方(農(nóng)田型和草地型樣地各7個(gè)樣方)，采用隨機(jī)抽樣法取樣，每個(gè)樣方隨機(jī)選取20株植株，記錄植株上芫菁科昆蟲數(shù)量，并進(jìn)行空間分布型測(cè)定。

1.3 分析方法

1.3.1 空間分布型采用聚集度指標(biāo)法(王智，2007；賈彥霞等，2019)、Taylor冪函數(shù)法(Taylor，1961；楊紫涵等，2020)和Iwao回歸分析法(Iwao，1968；劉延超等，2018)對(duì)芫菁科主要藥用昆蟲的空間分布型進(jìn)行測(cè)定。

聚集度指標(biāo)法：計(jì)算得到各樣方昆蟲的平均蟲口密度(m)、樣本方差(S2)、平均擁擠度m*。再按照聚集度指標(biāo)的計(jì)算公式分別求出擴(kuò)散系數(shù)(C)：C=S2/m、聚塊性指數(shù)m*/m、負(fù)二項(xiàng)分布k值：k=m2/(S2-m)、叢生指數(shù)(I)：I=(S2-m)/m、聚集度指數(shù)(Ca)：Ca=1/k=(S2-m)/m2和La指標(biāo)：La=m-m/S2+1。

Taylor冪函數(shù)法：根據(jù)Taylor(1961)提出的m與S2的回歸關(guān)系公式：S2=lga+blgm，式中，a和b分別表示抽樣因素和聚集特征。當(dāng)lga=0、b=1時(shí)，種群為隨機(jī)分布；當(dāng)lga>0、b=1時(shí)，種群在一切密度下呈聚集分布，且不依賴于密度；當(dāng)lga>0、b>1時(shí)，種群在一切密度下呈聚集分布，且依賴于密度聚集；當(dāng)lga<0、b<1時(shí)，種群為均勻分布，且密度越高，分布越均勻。

Iwao回歸分析法：m*-m回歸模型為：m*=α+βm，m*和m表示平均擁擠度與平均密度。式中，α為分布基本成分的分布性質(zhì)(平均擁擠程度)，當(dāng)α=0時(shí)，分布的基本成分為單個(gè)個(gè)體；當(dāng)α>0時(shí)，個(gè)體間相互吸引，分布的基本成分為個(gè)體群；當(dāng)α<0時(shí)，個(gè)體間相互排斥。β為基本成分的空間分布型，當(dāng)β=1時(shí)，種群為隨機(jī)分布；當(dāng)β>1時(shí)，種群為聚集分布；當(dāng)β

1.3.2 聚集原因分析應(yīng)用Blackith(1958)提出的聚集均數(shù)λ對(duì)芫菁科主要藥用昆蟲的聚集因素進(jìn)行分析。聚集均數(shù)公式為：λ=mγ/2k，式中，k為各樣方昆蟲負(fù)二項(xiàng)分布指數(shù)，m為平均值，γ為df=2k、P=0.05時(shí)的χ2值。當(dāng)λ≥2時(shí)，聚集主要是昆蟲的本身習(xí)性和環(huán)境條件綜合影響的結(jié)果；當(dāng)λ<2時(shí)，聚集是由環(huán)境因素引起(Blackith，1958；閻雄飛等，2019)。

1.3.3 抽樣技術(shù)根據(jù)Iwao(1968)提出的理論抽樣數(shù)計(jì)算公式：N=t2/D2[(α+1)/m+β-1]，計(jì)算芫菁科主要藥用昆蟲在不同蟲口密度下的最適抽樣數(shù)。式中，N為最適抽樣數(shù)，D為允許誤差，取t=1.96(95%可靠度)，m為平均蟲口密度(張治科等，2018)。

1.4 數(shù)據(jù)統(tǒng)計(jì)分析

運(yùn)用SPSS 20.0進(jìn)行數(shù)據(jù)統(tǒng)計(jì)分析，運(yùn)用Heat Mapper進(jìn)行聚類熱點(diǎn)圖分析。

2 結(jié)果與分析

2.1 哈思山區(qū)芫菁科主要藥用昆蟲的分布

2.1.1 芫菁科主要藥用昆蟲種類及分布情況對(duì)所有樣地中的芫菁科藥用昆蟲進(jìn)行鑒定，常見有西北斑芫菁Mylabrissibirica(Fischer von Waldheim, 1823)、麗斑芫菁M.speciosa(Pallas, 1781)、圓點(diǎn)斑芫菁M.aulica(Ménétriés, 1832)、暗頭豆芫菁Epicautaobscurocephala(Reitter, 1873)、扁角豆芫菁E.impressicornis(Pic, 1913)和綠芫菁Lyttacaraganae(Pallas, 1781)。調(diào)查發(fā)現(xiàn)，在不同類型樣地中芫菁科主要藥用昆蟲數(shù)量及種類均存在一定差異(表1)。該區(qū)域芫菁科昆蟲分布于農(nóng)田型樣地和草地型樣地，且農(nóng)田型樣地分布數(shù)量遠(yuǎn)多于草地型樣地。西北斑芫菁、麗斑芫菁、圓點(diǎn)斑芫菁、暗頭豆芫菁、扁角豆芫菁及綠芫菁在2種類型樣地中均有分布，其中，西北斑芫菁在2種樣地中最常見、數(shù)量最多，圓點(diǎn)斑芫菁次之。森林型和荒漠型樣地中均未發(fā)現(xiàn)芫菁科昆蟲分布。有研究報(bào)道，芫菁科昆蟲啃食沙棘Hippophaerhamnoides、沙棗Elaeagnusangustifolia、錦雞兒Caraganamicrophylla等的葉片(王文和，2015)，但本研究的沙棘林及其他森林型樣地中均未發(fā)現(xiàn)芫菁科昆蟲分布。

表1 樣地類型與昆蟲分布Table 1 Plot types and insect distribution

2.1.2 芫菁科主要藥用昆蟲在植物上的分布調(diào)查發(fā)現(xiàn),芫菁科主要藥用昆蟲分布的植物種類較多(表2)，共10科17種植株，集中于菊科Asteraceae、豆科Leguminosae和唇形科Lamiaceae植物。其中，西北斑芫菁分布的植物種類最豐富，共計(jì)11種，是區(qū)域內(nèi)最常見的芫菁類昆蟲。根據(jù)昆蟲在不同植物上的分布情況，對(duì)6種芫菁類昆蟲進(jìn)行聚類分析，從分布熱點(diǎn)圖(圖2)可以看出，芫菁類昆蟲可分為2類，西北斑芫菁和圓點(diǎn)斑芫菁為一類，綠芫菁、扁角豆芫菁、麗斑芫菁和暗頭豆芫菁為另一類。

表2 植物類型與芫菁科昆蟲分布Table 2 Plant types and insect distribution of Meloidae

這是由于西北斑芫菁和圓點(diǎn)斑芫菁對(duì)植物的選擇相似，2種昆蟲多分布在同一樣地，均主要取食密花香薷Elsholtziadensa、紫苜蓿Medicagosativa、乳苣Mulgediumtataricum等植物。而其他4種芫菁分布相對(duì)零散，其中綠芫菁和扁角豆芫菁分布的植物種類相對(duì)較少，綠芫菁主要取食野豌豆Viciasepium，該地區(qū)綠芫菁活動(dòng)時(shí)間集中在豆類植物開花期間，是危害豆類作物生產(chǎn)最主要害蟲；扁角豆芫菁主要分布在黃芪Astragalusmembranaceus植株上，同時(shí)該地區(qū)大量種植紫苜蓿作為牧草，其在紫苜蓿地中也分布；麗斑芫菁主要分布在串鈴草Phlomismongolica、冰草Agropyroncristatum等上；暗頭豆芫菁分布范圍較廣，在農(nóng)田型和草地型樣地中均有分布，主要分布在河北石頭花Gypsophilatschiliensis和白花枝子花Dracocephalumheterophyllum的植株上，但總體數(shù)量較少。

2.2 芫菁科主要藥用昆蟲的空間分布型

2.2.1 聚集度指標(biāo)測(cè)定結(jié)果對(duì)各樣方的調(diào)查數(shù)據(jù)進(jìn)行聚集度指標(biāo)計(jì)算(表3)，各樣方的m*>m、C>1、I>0、Ca>0、m*/m>1、La/m>1、k>0，各項(xiàng)指標(biāo)表明芫菁科主要藥用昆蟲的空間分布型為聚集分布。

表3 芫菁科主要藥用昆蟲聚集度指標(biāo)Table 3 Aggregation indexes of the main medicinal insects of Meloidae

2.2.2 Taylor冪函數(shù)法則分析將m和S2進(jìn)行對(duì)數(shù)轉(zhuǎn)換，分析二者之間線性關(guān)系(圖3：a)，芫菁科主要藥用昆蟲Taylor冪函數(shù)方程為：lgS2=0.204 2+1.103 8lgm(r=0.886)，式中l(wèi)ga=0.204 2>0、b=1.103 8>1，芫菁科主要藥用昆蟲種群在一切密度下呈聚集分布，且具密度依賴性，聚集強(qiáng)度隨種群密度的升高而增加，與應(yīng)用聚集度指標(biāo)法測(cè)定的結(jié)果一致。

2.2.3 Iwao回歸分析將m和m*進(jìn)行線性回歸(圖3：b)，m*-m模型回歸方程為：m*=0.642 5+1.069 7m(r=0.916)，由α=0.642 5>0、β=1.0697>1可知，芫菁科主要藥用昆蟲的空間分布型為聚集分布，種群分布的基本成分是個(gè)體群，個(gè)體間相互吸引，與采用聚集度指標(biāo)法和Taylor冪法則分析法測(cè)定的種群空間分布型的結(jié)果一致。

2.2.4 聚集原因分析根據(jù)昆蟲的聚集原因分析(表4)，各樣方中芫菁科主要藥用昆蟲的λ>2，表明其聚集原因可能是由其本身的取食、產(chǎn)卵等習(xí)性及溫度、日照等周圍環(huán)境條件綜合影響或其中之一引起。

表4 芫菁科主要藥用昆蟲聚集因素分析Table 4 Aggregation factor analysis of the main medicinal insects of Meloidae

2.3 抽樣技術(shù)

2.3.1 最適抽樣數(shù)利用理論抽樣數(shù)計(jì)算公式，取置信水平t=1.96(概率為0.95)，芫菁科主要藥用昆蟲理論抽樣模型為：N=(6.310/m+0.268)/D2。當(dāng)D為0.1、0.2和0.3時(shí)，m在1、2、4、6、8、10、15……50的條件下，計(jì)算得出不同密度條件下的最適抽樣數(shù)(表5)。在允許誤差一致的情況下，隨著平均蟲口密度的增加，芫菁科主要藥用昆蟲的最適抽樣數(shù)逐漸減少；在平均蟲口密度相同的情況下，隨著允許誤差的增大，芫菁科主要藥用昆蟲的最適抽樣數(shù)逐漸減少。

表5 芫菁科主要藥用昆蟲不同密度下最適抽樣數(shù)Table 5 Optimum sampling number under different densities of the main medicinal insects of Meloidae

表6 不同防治指標(biāo)下芫菁科主要藥用昆蟲序貫抽樣模型Table 6 Sequential sampling model of the main medicinal insects of Meloidae under different control indexes

表7 芫菁科主要藥用昆蟲序貫抽樣結(jié)果Table 7 Sequential sampling results of the main medicinal insects of Meloidae

3 討論

在甘肅哈思山區(qū)不同類型樣地中，調(diào)查發(fā)現(xiàn)的芫菁科藥用昆蟲共3屬6種，分布于農(nóng)田型和草地型樣地中，西北斑芫菁、麗斑芫菁、圓點(diǎn)斑芫菁、暗頭豆芫菁、扁角豆芫菁及綠芫菁在2種類型樣地中均有分布。在調(diào)查過程中發(fā)現(xiàn)，該區(qū)域草地型樣地中適宜芫菁科昆蟲取食的植物較農(nóng)田型樣地少，多分布在禁牧區(qū)和保護(hù)區(qū)草地，整體分布數(shù)量遠(yuǎn)少于農(nóng)田型樣地。此外，有研究發(fā)現(xiàn)芫菁科昆蟲與蝗蟲、蜜蜂存在著明顯的相關(guān)性(胡子淵，2020)，其幼蟲主要取食蝗蟲卵或蜜蜂幼蟲，而農(nóng)田中密集分布的作物及其他植物通常吸引了大量的蝗蟲和蜜蜂。這可能是造成該區(qū)域芫菁科主要藥用昆蟲在農(nóng)田型和草地型樣地中分布數(shù)量差異的原因。

昆蟲的空間分布型是由昆蟲種群的生物學(xué)、生態(tài)學(xué)特性和所處生境條件所決定的空間分布方式(王小武等，2017；魯武鋒等，2018)。本研究通過聚集度指標(biāo)法、Taylor冪函數(shù)法及Iwao回歸法分析哈思山區(qū)芫菁科藥用昆蟲的分布特征，結(jié)果均表明芫菁科藥用昆蟲的空間分布型為聚集分布，分布的基本成分是個(gè)體群，與李亞林(2010)研究中描述芫菁科昆蟲具有群聚性相符合。芫菁科昆蟲為復(fù)變態(tài)昆蟲，雌蟲產(chǎn)卵量大，平均產(chǎn)卵達(dá)500余粒(劉洋洋，2016)。幼蟲為肉食性，主要取食蝗蟲卵，部分取食蜂類幼蟲。成蟲為植食性，喜食植物的花和嫩葉，具有集群危害的特性，調(diào)查過程中甚至發(fā)現(xiàn)不同種類的芫菁科昆蟲會(huì)聚集在一起取食，如西北斑芫菁與圓點(diǎn)斑芫菁多分布在同一樣地，這可能屬于營(yíng)養(yǎng)群集或“婚期”群集(楊玉霞，2007)。芫菁科藥用昆蟲的聚集原因究竟是本身的取食、產(chǎn)卵、交配等習(xí)性及溫度、日照等周圍環(huán)境條件綜合影響，還是其中之一所引起，仍需進(jìn)一步研究。在空間分布型研究的基礎(chǔ)上，本研究對(duì)芫菁科昆蟲的抽樣技術(shù)進(jìn)行分析；在同一允許誤差內(nèi)，隨著平均蟲口密度的增加，最適抽樣數(shù)逐漸減少。不同蟲口密度下的最適理論抽樣數(shù)隨成蟲數(shù)量的增大而減小。因而在調(diào)查時(shí)，通過預(yù)估植株上芫菁科藥用昆蟲的蟲口密度，確定抽樣的數(shù)量，當(dāng)蟲口密度較低時(shí)，應(yīng)適當(dāng)增加長(zhǎng)植株的抽樣株數(shù)，反之，則可適當(dāng)減少。同時(shí)還確定了芫菁科昆蟲的理論抽樣模型，建立了不同防治指標(biāo)下芫菁科昆蟲的序貫抽樣模型及一定允許誤差和不同蟲口密度下的序貫抽樣結(jié)果表。

昆蟲是動(dòng)物藥的重要來源之一，是中藥的重要組成，占有獨(dú)特地位。然而，與藥用植物相比，有關(guān)藥用昆蟲資源方面的研究相對(duì)薄弱(韋桂寧等，2019；張星賢等，2019)。由于其動(dòng)態(tài)的生活史，在時(shí)間和空間上的不確定導(dǎo)致研究的不易性，嚴(yán)重制約藥用昆蟲資源的開發(fā)利用。芫菁科藥用昆蟲種類眾多，形成的中藥或民間藥種類豐富，如中藥品種斑蝥、青娘子、葛上亭長(zhǎng)、地膽等(黃文華，2001)，但是現(xiàn)代藥學(xué)研究集中在大斑溝芫菁和眼斑溝芫菁分泌的斑蝥素及其衍生物方面(劉亞楠，2013)。一方面可能因過度利用造成2種芫菁野生資源一定程度的破壞，另一方面卻使得同類資源大量浪費(fèi)。因此對(duì)于藥用昆蟲而言，針對(duì)其同科同屬或者近緣種屬大類資源進(jìn)行系統(tǒng)研究，既能夠?yàn)獒t(yī)藥產(chǎn)業(yè)提供豐富的資源，也不易造成資源的過度使用和破壞。

本研究建立了芫菁科主要藥用昆蟲農(nóng)田型樣地理論抽樣模型及序貫抽樣模型，分析了該地區(qū)芫菁科主要藥用昆蟲分布特征和生活習(xí)性，研究結(jié)果為該類藥用昆蟲資源的監(jiān)測(cè)提供參考，一方面有利于其可持續(xù)應(yīng)用，一方面可為農(nóng)業(yè)生產(chǎn)中該類昆蟲數(shù)量過度造成的危害提供防治理論依據(jù)。