榴梿重要害蟲(chóng)

摘 要：本研究報(bào)道了榴梿重要害蟲(chóng)——馬來(lái)異裂木虱[Allocarsidara malayensis （Crawford）]在海南的發(fā)現(xiàn)，詳細(xì)描述了其形態(tài)特征，克隆了其COI 基因，描述了其對(duì)榴梿的危害情況，并預(yù)測(cè)了該害蟲(chóng)在我國(guó)的潛在適生區(qū)。通過(guò)25 個(gè)有效分布點(diǎn)位和19 個(gè)氣候環(huán)境因素，采用最大熵（MaxEnt）模型結(jié)合地理信息系統(tǒng)（ArcGIS），預(yù)測(cè)馬來(lái)異裂木虱的潛在適生區(qū)，并使用受試者工作特征（receiver operating characteristic， ROC）曲線分析法進(jìn)行驗(yàn)證，同時(shí)運(yùn)用刀切法分析影響馬來(lái)異裂木虱分布的主要環(huán)境因素。結(jié)果表明：馬來(lái)異裂木虱的潛在分布區(qū)主要在東南亞、南亞、非洲、南美州和北美洲的赤道附近的熱帶地區(qū)。我國(guó)海南、廣西、廣東、福建、臺(tái)灣、云南、四川、貴州均存在馬來(lái)異裂木虱的適生區(qū)，且海南全省及臺(tái)灣部分沿海地區(qū)為高適生區(qū)。在影響馬來(lái)異裂木虱潛在地理分布的環(huán)境因素中，最濕月份的降水量（Bio13）和最濕季度的降水量（Bio16）為最重要影響變量。由于馬來(lái)異裂木虱的適生區(qū)與榴梿的種植區(qū)較為一致，為防止馬來(lái)異裂木虱的進(jìn)一步傳播和擴(kuò)散，保障我國(guó)榴梿產(chǎn)業(yè)的健康發(fā)展，應(yīng)加強(qiáng)馬來(lái)異裂木虱的監(jiān)測(cè)，同時(shí)開(kāi)展該害蟲(chóng)的生物學(xué)、生態(tài)學(xué)特性及消長(zhǎng)規(guī)律研究，研究結(jié)果為該蟲(chóng)經(jīng)濟(jì)高效防治策略的制定提供理論依據(jù)。

關(guān)鍵詞：榴梿；馬來(lái)異裂木虱；形態(tài)特征；COI 基因；MaxEnt 模型；適生區(qū)

中圖分類(lèi)號(hào)：S436.67 文獻(xiàn)標(biāo)志碼：A

榴梿（Durio zibethinus Rumph. ex Murray）俗名榴蓮，屬錦葵科（Malvaceae）榴梿屬（Durio），原產(chǎn)于馬來(lái)西亞、印度尼西亞、文萊，現(xiàn)主要種植于泰國(guó)、馬來(lái)西亞、印度尼西亞等地，其果肉營(yíng)養(yǎng)豐富、味道濃郁且風(fēng)味獨(dú)特，享有“熱帶果王”的美譽(yù)[1-3]。我國(guó)自1958 年開(kāi)始從馬來(lái)西亞、泰國(guó)、越南等東南亞國(guó)家引進(jìn)榴梿種苗進(jìn)行試種，但少有開(kāi)花結(jié)果的報(bào)道[4]。2019 年在我國(guó)海南保亭成功掛果，表明海南部分地區(qū)存在榴梿大面積推廣種植的可能[5]。2024 年在西雙版納成功掛果，使榴梿成為海南乃至我國(guó)熱區(qū)高度關(guān)注的水果，并在海南、云南、廣東、廣西等地掀起了種植熱潮[6-7]。

2024 年4 月，在開(kāi)展海南榴梿害蟲(chóng)調(diào)查中，筆者在保亭什玲鎮(zhèn)的榴梿園發(fā)現(xiàn)一種嚴(yán)重危害榴梿嫩梢的木虱，通過(guò)形態(tài)特征和雌雄蟲(chóng)外生殖器的鑒定，確認(rèn)其為馬來(lái)異裂木虱[Allocarsidaramalayensis （Crawford）]。馬來(lái)異裂木虱分布于馬來(lái)西亞、泰國(guó)、新加坡、印度、印度尼西亞，為榴梿的重要害蟲(chóng)之一，在泰國(guó)、馬來(lái)西亞和新加坡是榴梿害蟲(chóng)防治的主要對(duì)象，澳大利亞將其列為榴梿檢疫對(duì)象之一，此前的文獻(xiàn)記載中尚無(wú)該蟲(chóng)在中國(guó)的分布記錄[8-13]。為快速、準(zhǔn)確地識(shí)別馬來(lái)異裂木虱，本研究詳細(xì)描述了馬來(lái)異裂木虱的形態(tài)特征、為害特點(diǎn)，并對(duì)其COI 基因序列進(jìn)行了研究。

外來(lái)入侵生物潛在分布預(yù)測(cè)是制定物種防控策略以阻止其蔓延的重要途徑之一[14]。最大熵（MaxEnt）模型可通過(guò)已知物種分布點(diǎn)數(shù)據(jù)和環(huán)境變量，精確地預(yù)測(cè)物種的潛在分布，且具有操作簡(jiǎn)便、所需樣本量少、運(yùn)行結(jié)果穩(wěn)定、與實(shí)際分布更相符等特點(diǎn)，是目前最常用的物種分布預(yù)測(cè)模型之一，廣泛應(yīng)用于動(dòng)植物保護(hù)、病蟲(chóng)害防治指導(dǎo)、入侵物種風(fēng)險(xiǎn)評(píng)估等研究中[15-19]。本研究通過(guò)采集和收集的馬來(lái)異裂木虱分布數(shù)據(jù)，結(jié)合相關(guān)環(huán)境因素，利用MaxEnt 模型和地理信息系統(tǒng)（ArcGIS）預(yù)測(cè)馬來(lái)異裂木虱在我國(guó)的潛在適生區(qū)域，以期為預(yù)防其進(jìn)一步擴(kuò)散和制定合理的控制措施提供科學(xué)依據(jù)。

1 材料與方法

1.1 材料

1.1.1 供試?yán)ハx(chóng)

從保亭縣什玲鎮(zhèn)、保亭縣南茂農(nóng)場(chǎng)、三亞市育才生態(tài)區(qū)優(yōu)旗大嘴鳥(niǎo)榴梿一號(hào)基地采回帶木虱的嫩葉，靠接于養(yǎng)蟲(chóng)籠內(nèi)的抽梢榴梿苗中，露天繼續(xù)培養(yǎng)，獲得成蟲(chóng)用于形態(tài)鑒定和分子檢測(cè)。成蟲(chóng)接入抽梢的榴梿苗中，使其自由交配、產(chǎn)卵，用于卵、若蟲(chóng)形態(tài)觀察。

1.1.2 主要儀器

Olympus SZ61 光學(xué)體式顯微鏡，M205FA 全視角攝像顯微鏡，LABGIC TH-Mini微型手持均質(zhì)儀，數(shù)顯恒溫水浴鍋，湘儀CHT210R高速離心機(jī)，賽默飛NANODROP 2000c 超微量分光光度計(jì)，伯樂(lè)C1000 PCR 擴(kuò)增儀，伯樂(lè)電泳儀，SAGECREATION 凝膠成像系統(tǒng)。

1.1.3 主要試劑

乙酸乙酯，血液/細(xì)胞/組織基因組DNA 提取試劑盒[天根生化科技（北京）有限公司]，2×Es Taq MasterMix（康為世紀(jì)生物科技股份有限公司）；6×DNA 上樣緩沖液（北京索萊寶科技有限公司）； 綠色熒光核酸染料（10 000×）（北京索萊寶科技有限公司）；DNAMarker Ⅲ分子量標(biāo)準(zhǔn)（100～5000 bp）（北京蘭杰柯科技有限公司）；50×TAE 緩沖液（北京索萊寶科技有限公司）；瓊脂糖（Biowest）。

1.1.4 數(shù)據(jù)來(lái)源

馬來(lái)異裂木虱的分布數(shù)據(jù)通過(guò)以下方式獲得：（1）實(shí)地考察。2024 年4—6 月，調(diào)查三亞、保亭、五指山、瓊海等地的榴梿園，發(fā)現(xiàn)馬來(lái)異裂木虱即采用GPS 記錄經(jīng)緯度，為其分布預(yù)測(cè)提供數(shù)據(jù)支持。（2）數(shù)據(jù)庫(kù)查詢(xún)物種分布。查詢(xún)的數(shù)據(jù)庫(kù)包括全球生物多樣性信息網(wǎng)絡(luò)（Global Biodiversity Information Facility， GBIF，http：//www.Gbif.org/）和自然愛(ài)好者社區(qū)數(shù)據(jù)庫(kù)（iNaturalist， http：//www. inaturalist.org/）。將收集到的馬來(lái)異裂木虱分布數(shù)據(jù)，去除無(wú)效數(shù)據(jù)后，最終獲得25 個(gè)馬來(lái)異裂木虱的有效分布點(diǎn)位，并將點(diǎn)位數(shù)據(jù)轉(zhuǎn)換成.csv 格式保存。

1.1.5 環(huán)境變量

環(huán)境數(shù)據(jù)來(lái)自WorldClim 網(wǎng)站（http：//www.world.org/）1970—2000 年的氣候變量，得到19 個(gè)分辨率為5 m 的氣候因子，即年平均溫度（Bio1）、平均溫度日較差（Bio2）、等溫性（Bio3）、氣溫季節(jié)性變動(dòng)系數(shù)（Bio4）、最熱月份的最高溫度（Bio5）、最冷月份的最低溫度（Bio6）、年溫度變化范圍（Bio7）、最濕季度的平均溫度（Bio8）、最干季度的平均溫度（Bio9）、最熱季度的平均溫度（Bio10）、最冷季度的平均溫度（Bio11）、年降水量（Bio12）、最濕月份的降水量（Bio13）、最干月份的降水量（Bio14）、降水量季節(jié)性變化（Bio15）、最濕季度的降水量（Bio16）、最干季度的降水量（Bio17）、最暖季度的降水量（ Bio18 ）、最冷季度的降水量（Bio19）。

1.2 方法

1.2.1 形態(tài)學(xué)鑒定

木虱成蟲(chóng)樣品采用乙酸乙酯毒殺后，粘于三角紙片上保存于標(biāo)本盒中。若蟲(chóng)樣品浸泡于95%酒精中。在處理雌雄蟲(chóng)生殖器、后足和頭部等重要結(jié)構(gòu)時(shí)，先使用10%氫氧化鈉溶液溶解肌肉組織，然后用蒸餾水漂洗干凈，最后置于乙醇和甘油混合液中，便于鏡檢。使用Olympus SZ61 光學(xué)體式顯微鏡觀察成蟲(chóng)和若蟲(chóng)；成蟲(chóng)、若蟲(chóng)使用M205FA 全視角攝像顯微鏡拍照，并用Leica LAS X1 軟件測(cè)量體長(zhǎng)、體寬。物種鑒定和分類(lèi)的依據(jù)主要為形態(tài)學(xué)，在分類(lèi)地位上依據(jù)Daniel Burckhardt 的最新分類(lèi)系統(tǒng)[20]。

1.2.2 分子檢測(cè)

用組織基因組DNA 提取試劑盒提取采集的木虱樣本DNA，然后利用無(wú)脊椎動(dòng)物COI 基因通用引物（ LCO1490 5–GGTCAAATCATAAAGATATTGG-3 ， HCO21985–TAAACTTCAGGGTGACCAAAAAATCA-3）[21]進(jìn)行PCR 擴(kuò)增，PCR 反應(yīng)體系：6 μL TemplateDNA（10 ng/μL），25 μL 2×Es Tag PCR MasterMix（Dye），2 μL L-Primer （10 μmol/L），2 μL R-Primer（10 μmol/L），15 μL ddH2O。PCR 擴(kuò)展程序：94 ℃預(yù)變性5 min；94 ℃變性30 s，50 ℃退火30 s，72 ℃延伸30 s，30 個(gè)循環(huán)；72 ℃延伸 2 min。PCR產(chǎn)物用1.5%瓊脂糖凝膠電泳檢測(cè)，并將PCR 產(chǎn)物送生工生物工程（上海）股份有限公司進(jìn)行測(cè)序。

1.2.3 馬來(lái)異裂木虱的適生區(qū)預(yù)測(cè)

MaxEnt 模型預(yù)測(cè)能力驗(yàn)證：MaxEnt 模型預(yù)測(cè)的準(zhǔn)確度和可靠性采用訓(xùn)練集和測(cè)試集數(shù)據(jù)的受試者工作特征（receiver operating characteristic， ROC）曲線進(jìn)行檢驗(yàn)。ROC 曲線與橫坐標(biāo)圍成的面積為曲線下面積（area under curve， AUC），AUC 值越大表示與隨機(jī)分布相距越遠(yuǎn)，環(huán)境變量與預(yù)測(cè)的物種地理分布模型之間相關(guān)性越大，即模型預(yù)測(cè)效果越好[22]。當(dāng)AUC≥0.9，表示模型預(yù)測(cè)結(jié)果優(yōu)秀；當(dāng)0.7≤AUClt;0.9，表示模型預(yù)測(cè)結(jié)果良好；當(dāng)0.5≤AUClt;0.7，表示模型預(yù)測(cè)結(jié)果可接受；當(dāng)AUClt;0.5，表示模型預(yù)測(cè)無(wú)效[23]。

采用MaxEnt 3.4.3 模型軟件結(jié)合地理信息系統(tǒng)（ArcGIS 10.8），預(yù)測(cè)馬來(lái)異裂木虱的潛在適生區(qū)域。利用ArcGIS 10.8 軟件中的SDMtoolbox 插件Conversion Tools 將生物氣候因子環(huán)境圖層轉(zhuǎn)換成MaxEnt 軟件可用的ASCII 格式。將處理的氣候因子?xùn)鸥駭?shù)據(jù)以及分布點(diǎn)位導(dǎo)入MaxEnt3.4.3 軟件，以bootstrap 模式運(yùn)行10 次，模擬結(jié)果取平均值，得到歷史氣候條件下馬來(lái)異裂木虱在全球范圍內(nèi)的潛在適生區(qū)數(shù)據(jù)。將MaxEnt 運(yùn)行的結(jié)果導(dǎo)入ArcGIS 10.8 軟件，在圖層屬性的符號(hào)系統(tǒng)中，將柵格數(shù)據(jù)進(jìn)行分類(lèi)：0～0.1 為非適生區(qū)、0.1～0.4 為低適生區(qū)、0.4～0.7 為中適生區(qū)、0.7～1.0 為高適生區(qū)，并以不同顏色標(biāo)記，得到歷史氣候條件下全球范圍內(nèi)馬來(lái)異裂木虱的適生區(qū)分布圖。在ArcGIS 10.8 軟件中，將馬來(lái)異裂木虱適生區(qū)數(shù)據(jù)進(jìn)行裁切，并導(dǎo)入中國(guó)標(biāo)準(zhǔn)地圖，得到馬來(lái)異裂木虱在我國(guó)的適生區(qū)分布圖。

2 結(jié)果與分析

2.1 鑒別特征

頭約與胸等寬。頭頂寬略大于長(zhǎng)，近方形。頭兩側(cè)側(cè)單眼與觸角基部間各具一明顯的角狀凸起。觸角窩下頰面膨凸，超過(guò)頭頂，略圓，覆蓋額。唇基長(zhǎng)，端部窄圓（圖1A）。觸角約為頭寬的4 倍，第3 節(jié)明顯粗于其他節(jié)。胸部窄。足細(xì)長(zhǎng)。后足脛節(jié)端距5 枚。

2.1.1 雄蟲(chóng)

體黃綠色至黃褐色，體長(zhǎng)2.79～3.47 mm，體寬0.67～0.77 mm。前翅長(zhǎng)卵形，頂角略尖，具翅痣。第1 徑脈（R1）在中部彎折成一鈍角，其后半段與徑分脈（Rs）前半段平行；一條橫脈（rs-m）連接Rs 脈與中脈（M），與M 脈相連處位于M 脈分支處；緣紋2 條，分別位于m1室和m2 室（圖1B）。下生殖板腹緣基部稍平，端部稍突，背緣膨突（圖1C）；陽(yáng)基側(cè)突側(cè)視端部平截，近端部處向后膨突；載肛突后葉發(fā)達(dá)；陽(yáng)莖端節(jié)稍呈波形，端部彎起似勺狀（圖1D，圖1E）。

2.1.2 雌蟲(chóng)

體黃綠色至黃褐色，體長(zhǎng) 3.40～4.04 mm，體寬0.71～0.79 mm。生殖節(jié)約為腹部長(zhǎng)度的2/3；載肛突末端鈍圓稍上翹，環(huán)肛孔長(zhǎng)度超過(guò)載肛突的1/2；亞生殖板窄錐狀，末端尖細(xì)；腹瓣基突長(zhǎng)錐形，基部鈍而端部尖；產(chǎn)卵瓣基本等長(zhǎng)，但背瓣端鈍，而腹瓣端尖（圖1F，圖2）。

2.1.3 卵

呈橢圓形，表面光滑，頂端尖削，基部較圓頓（圖2）。

2.1.4 若蟲(chóng)

若蟲(chóng)共5 個(gè)齡期（圖2），體型和體色隨齡期的變化而變化。5 齡若蟲(chóng)長(zhǎng)2.02～3.00 mm，全身覆白色蠟粉，頭、胸、腹、翅芽、足近黃褐色透明，側(cè)緣橙黃色，頭至尾部具1 條橙色背中線。觸角10 節(jié)，梗節(jié)和柄節(jié)短而粗壯，鞭節(jié)除末端2 節(jié)為黑褐色，其余各節(jié)從基部黃褐色逐漸加深到淺黑褐色。翅芽末端長(zhǎng)至第2 腹節(jié)，腹部側(cè)緣分泌白色蠟絲，且以靠近腹部末端的第2 對(duì)蠟絲最長(zhǎng)。

2.2 分子檢測(cè)

利用COI 通用引物對(duì)木虱樣品的DNA 進(jìn)行PCR 擴(kuò)增，電泳檢測(cè)顯示清晰的單一條帶（圖3）。比對(duì)3 地樣品的COI 序列， 相似性為99.39%～99.85%，表明采集到的標(biāo)本為同一物種。去除引物與低質(zhì)量序列，獲得645 bp 的COI 序列（GenBank 登錄號(hào)為PQ550000）。

2.3 新擬中文名

最初置于喬木虱屬（ Tenaphalara ）， 為T(mén)enaphalara malayensis，此后的學(xué)者均按此學(xué)名進(jìn)行重新描述[24]。但在1985 年，YANG 等[12]在新發(fā)現(xiàn)了喬木虱屬的4 個(gè)新種后，通過(guò)形態(tài)學(xué)分析，對(duì)其原分類(lèi)地位產(chǎn)生了質(zhì)疑，認(rèn)為其應(yīng)為一個(gè)新屬。1987 年HOLLIS[25]正式將其歸為新建立的Allocarsidara 屬，屬征為頭頂側(cè)緣及后緣為圓形，位于頭頂前方的側(cè)單眼與觸角基之間各具一發(fā)達(dá)的瘤凸。世界共有Allocarsidara 屬6 種，A.malayensis 為該屬在中國(guó)分布的首次發(fā)現(xiàn)，根據(jù)其拉丁詞詞源含義擬中文名為“馬來(lái)異裂木虱”。

2.4 為害特征

馬來(lái)異裂木虱主要在榴梿的抽梢期為害，成蟲(chóng)將卵散產(chǎn)在尚未展開(kāi)的嫩葉葉肉中，葉表面可見(jiàn)黃色光滑小鼓包，卵孵化后形成凹陷的黃褐色斑。成蟲(chóng)和若蟲(chóng)取食榴梿嫩梢汁液，導(dǎo)致葉片出現(xiàn)黃色斑點(diǎn)，葉片發(fā)育不良、偏小，甚至扭曲變形，嚴(yán)重時(shí)葉片干枯，提前脫落，影響榴梿長(zhǎng)勢(shì)。此外，成蟲(chóng)、若蟲(chóng)均會(huì)分泌蜜露，引起煤污病，影響榴梿光合作用（圖4）。

2.5 馬來(lái)異裂木虱潛在地理分布MaxEnt 模型精度評(píng)價(jià)及主要環(huán)境因素分析

利用ROC 曲線AUC 值對(duì)模型預(yù)測(cè)結(jié)果進(jìn)行檢驗(yàn)。AUC 值越接近1 說(shuō)明模擬精確性越高，一般認(rèn)為AUC＞0.9，預(yù)測(cè)效果可信。且AUC 值越大，標(biāo)準(zhǔn)差越小，模型的穩(wěn)定性越高，且隨機(jī)取樣對(duì)模型影響小[26]。對(duì)MaxEnt 模型預(yù)測(cè)馬來(lái)異裂木虱適生區(qū)分布結(jié)果進(jìn)行檢驗(yàn)（圖5），結(jié)果表明，歷史氣候條件下，模型的AUC 平均值為0.991±0.003，表明該模型預(yù)測(cè)結(jié)果的可信度與穩(wěn)定性較高，可用于預(yù)測(cè)馬來(lái)異裂木虱的潛在地理分布。

基于MaxEnt 模型正則化訓(xùn)練增益的刀切法檢驗(yàn)結(jié)果表明（圖6），僅此變量中正則化訓(xùn)練增益最高的4 個(gè)環(huán)境因素是最冷月份的最低溫度（Bio6）、年溫度變化范圍（Bio7）、最濕月份的降水量（Bio13）和最濕季度的降水量（Bio16），表明馬來(lái)異裂木虱對(duì)這些因子具有較高的依賴(lài)性。且最濕季度的降水量（Bio16）和最濕月份的降水量（Bio13）對(duì)模型的貢獻(xiàn)最高，貢獻(xiàn)率分別為25.3%和15.5%。除此變量中正則化訓(xùn)練增益值最短的是最暖季度的降水量（Bio18）和平均溫度日較差（Bio2），表明這2 個(gè)環(huán)境因素對(duì)馬來(lái)異裂木虱分布的影響具有最多其他變量所沒(méi)有的信息，即這2 個(gè)環(huán)境因素在馬來(lái)異裂木虱適生區(qū)的預(yù)測(cè)中具有更多獨(dú)有信息。

2.6 馬來(lái)異裂木虱適生區(qū)預(yù)測(cè)分析

馬來(lái)異裂木虱在全球的適生區(qū)預(yù)測(cè)結(jié)果如圖7 所示，馬來(lái)異裂木虱潛在世界分布區(qū)主要在東南亞、南亞、非洲、南美州和北美洲的赤道附近的熱帶地區(qū)。我國(guó)南方一些省份存在馬來(lái)異裂木虱的適生區(qū)，其中云南、四川、貴州的部分地區(qū)為馬來(lái)異裂木虱的低適生區(qū)；廣西、廣東、福建的大部地區(qū)為馬來(lái)異裂木虱的低適生區(qū)，小部地區(qū)為中適生區(qū)；海南省全境為高適生區(qū)；臺(tái)灣沿海地區(qū)存在高適生區(qū)（圖8）。

3 討論

通過(guò)形態(tài)分類(lèi)鑒定，確定了嚴(yán)重危害海南榴梿嫩葉的木虱為馬來(lái)異裂木虱，是該物種在我國(guó)的首次記錄。MaxEnt 模型預(yù)測(cè)結(jié)果表明，我國(guó)海南、臺(tái)灣、廣東、廣西、福建、云南、貴州和四川均存在馬來(lái)異裂木虱的適生區(qū)。

木虱是一類(lèi)寄主專(zhuān)一性強(qiáng)的昆蟲(chóng)，寄主通常限定在某一植物科甚至屬之內(nèi)[27-28]。已有報(bào)道表明，馬來(lái)異裂木虱分布于馬來(lái)西亞、泰國(guó)、新加坡、印度、印度尼西亞，但尚無(wú)文獻(xiàn)記載該蟲(chóng)在中國(guó)的分布記錄[8， 10-13]。馬來(lái)異裂木虱僅在榴梿嫩梢上取食和產(chǎn)卵，且常將卵產(chǎn)在尚未展開(kāi)的榴梿嫩葉葉肉中，尚未展開(kāi)的榴梿葉葉面緊貼，且葉背布滿青銅色的貼生鱗片[4， 29]，難以發(fā)現(xiàn)葉片上的卵，由此推測(cè)該蟲(chóng)可能是隨引進(jìn)的榴梿種苗進(jìn)入我國(guó)。

生態(tài)位模型是利用物種已知的分布信息和相關(guān)環(huán)境變量，根據(jù)一定的算法運(yùn)算構(gòu)建模型，判斷物種的生態(tài)需求，并將運(yùn)算結(jié)果投射至不同的時(shí)間和空間以預(yù)測(cè)物種的分布[30]。常見(jiàn)的物種生態(tài)模型包括規(guī)則集遺傳算法（GARP）、生態(tài)位因子分析（EnFA）、生物氣候模型（BIOCLIM）、區(qū)域環(huán)境模型（DOMAIN）和最大熵模型（MaxEnt）等[31]。MaxEnt 是一種基于最大熵原理的生態(tài)位評(píng)價(jià)模型，即使在物種分布數(shù)據(jù)較少的情況下，仍有較高的準(zhǔn)確率[32]。本研究采用MaxEnt 模型預(yù)測(cè)馬來(lái)異裂木虱在我國(guó)的適生區(qū)，結(jié)果表明，海南、臺(tái)灣、廣東、廣西、福建、云南、四川均存在馬來(lái)異裂木虱的適生區(qū)，這與我國(guó)榴梿的種植區(qū)分布較為一致，表明馬來(lái)異裂木虱在這些地區(qū)具有較大的分布空間和擴(kuò)散潛力。目前的調(diào)查結(jié)果表明馬來(lái)異裂木虱僅在海南少數(shù)榴梿園內(nèi)發(fā)現(xiàn)，為防止馬來(lái)異裂木虱在我國(guó)的進(jìn)一步擴(kuò)散蔓延，建議在進(jìn)行榴梿種苗引進(jìn)和調(diào)運(yùn)時(shí)需加強(qiáng)對(duì)馬來(lái)異裂木虱的檢測(cè)。

盡管馬來(lái)異裂木虱在泰國(guó)、馬來(lái)西亞等國(guó)家是榴梿的重要害蟲(chóng)，但目前僅查詢(xún)到有關(guān)藥劑篩選、施藥方式和生物學(xué)研究的報(bào)道，且不同地區(qū)的馬來(lái)異裂木虱在發(fā)育歷期、成蟲(chóng)壽命、產(chǎn)卵量等有所差異[2， 11， 13， 29， 33-35]。因此建議開(kāi)展馬來(lái)異裂木虱生物學(xué)及生態(tài)學(xué)特性的研究，挖掘本地天敵資源，掌握其擴(kuò)散機(jī)制和發(fā)生規(guī)律，建立適合我國(guó)榴梿生產(chǎn)實(shí)際的馬來(lái)異裂木虱防控措施。

參考文獻(xiàn)

[1] VERHEIJ E W M， CORONEL R E. Plant resource ofSouth-East Asia No. 2： edible fruits and nuts[M]. Leiden：Backhuys Publishers， 1991： 157-161.

[2] SAICHOL K， APINYA W， YOSSAPOL P， ROBERT E P.The durian： botany， horticulture， and utilization[M]//FERGUSON A R， PAULL R E. Horticultural reviews：Volume 47. Hoboken： John Wiley amp; Sons， Inc.， 2020：125-211.

[3] 梁廣勤， 胡學(xué)難， 趙菊鵬. 榴蓮山竹生產(chǎn)與病蟲(chóng)害防治[M].北京： 中國(guó)農(nóng)業(yè)出版社， 2016： 2-18.

LIANG G Q， HU X N， ZHAO J P. Production and pestscontrol of Durio zibethinus and Garcinia mangostana[M].Beijing： China Agriculture Press， 2016： 2-18. （in Chinese）

[4] 周兆禧， 林興娥. 中國(guó)榴蓮栽培技術(shù)[M]. 北京： 中國(guó)農(nóng)業(yè)科學(xué)技術(shù)出版社， 2022： 3-6.

ZHOU Z X， LIN X E. Durian cultivation techniques inChina[M]. Beijing： China Agricultual Science and TechnologyPress， 2022： 3-6. （in Chinese）

[5] 貓山王、金枕榴蓮在海南保亭試種成功[N]. 南方農(nóng)村報(bào)，2019-07-09（07）.

Musang King and Mon Thong were planted successful inBaoting of Hainan province[N]. Southern Rural Daily，2019-07-09（07）. （in Chinese）

[6] 羅江， 陳凱姿， 王軍鋒. 國(guó)產(chǎn)榴蓮嶄露頭角“榴蓮自由”更近了嗎？[N]. 經(jīng)濟(jì)參考報(bào)， 2024-07-18（004）.

LUO J， CHEN K Z， WAGN J F. Domestic durian is emerging，“Durian freedom” is closer？[N]. Economic InformationDaily， 2024-07-18（004）. （in Chinese）

[7] 竇崇尹. 云南景洪一村民成功試種榴蓮[N]. 云南日?qǐng)?bào)，2024-09-26.

DOU C Y. A villager successfully grows durian in Jinghongof Yunnan province[N]. Yunnan Daily， 2024-09-26. （inChinese）

[8] Global Biodiversity Information Facility[DB/OL]. [2024-08-20].https：//www.gbif.org/zh/.[9] Australian Quarantine and Inspection Service. Final importrisk analysis on the importation of fresh durian fruit （Duriozibethinus Murray） from the kingdom of Thailand[M]. Australia：Canberra， 1999： 18-30.

[10] SUTTIAROM S， WONGNIKONG W， PLODKONBURI W，MANASMUNKONG B. Efficacy of some insecticides forcontrolling durian psyllids （Allocaridara malayensis Crawford）[R]. Annual Research Report， Department of Agriculture，Bangkok， Thailand. 2012： 1036-1044.

[11] PE?A J E， SHARP J L， WYSOKI M. Tropical fruit pestsand pollinators： biology， economic importance， natural enemiesand control[M]. Trowbridge： Cromwell Press， 2002：316.

[12] YANG C K， LI F S. Notes on Tenaphalara and Tyora withdescriptions of five new species from China[J]. Entomotaxonomia， 1985， 7： 205-217.

[13] SUTTIAROM S， SIRISINGH S， NAMRUANGSRI W. Pestmanagement of durian insect[C]. Proceedings of the 2nd nationalplant protection conference. Bangkok （Thailand）：Plant protection for environmental quality， Thailand WeedScience Association， Thailand Entomology and ZoologyAssociation， Thailand Plant Protection Association， ThaiAssociation-Agricultural Chemistry Business， 1995： 116-131.

[14] UDEN D R， ALLEN C R， ANGELER D G， ANGELE D G，CORRAL L， FRICKE K A. Adaptive invasive species distributionmodels： a framework for modeling incipient invasions[J]. Biological Invasions， 2015， 10（17）： 2831-2850.

[15] PHILLIPS S J， ANDERSON R P， SCHAPIRE R E. Maximumentropy modeling of species geographic distributions[J].Ecological Modelling， 2016， 190（3/4）： 231-259.

[16] 孫紅云， 徐亮勝， 馮浩， 王程利， 黃麗麗. 基于MaxEnt 模型預(yù)測(cè)蘋(píng)果樹(shù)腐爛病在中國(guó)的潛在地理分布[J]. 西北農(nóng)業(yè)學(xué)報(bào)， 2020， 29（3）： 461-466.

SUN H Y， XU L S， FENG H， WANG C L， HUANG L L.Prediction for potential geographic distribution of Valsa maliin China based on MaxEnt model[J]. Acta Agriculturae BorealiOccidentalis Sinica， 2020， 29（3）： 461-466. （in Chinese）

[17] 張丹麗， 陳小艷， 李敏， 袁娟娟， 姜坤，卜文俊. 基于MaxEnt 模型預(yù)測(cè)黑頭叉胸花蝽在全球的潛在適生區(qū)[J].昆蟲(chóng)學(xué)報(bào)， 2024， 67（5）： 683-691.

ZHANG D L， CHEN X Y， Ll M， YUAN J J， JANG K， BUW J. Prediction of potential suitable habitats of Amphiareusobscuriceps （Hemiptera： Anthocoridae） in the world basedon MaxEnt model[J]. Acta Entomologica Sinica， 2024， 67（5）：683-691. （in Chinese）

[18] 趙清， 陳超， 范詩(shī)羽， 張利軍， 張虎芳. 基于MaxEnt 模型預(yù)測(cè)又角厲蝽在全球的潛在地理分布[J]. 植物保護(hù)學(xué)報(bào)，2022， 49（4）： 1210-1216.

ZHAO Q， CHEN C， FAN S Y， ZHANG L J， ZHANG H F.Prediction of potential distributionof predatory stink bugEocanthecona furcellata based on MaxEnt model[J]. Journalof Plant Diseases and Protection， 2022， 49（4）： 1210-1216. （inChinese）

[19] 梁莉， 冼曉青， 趙浩翔， 郭建洋， 劉萬(wàn)學(xué). 基于MaxEnt 模型的白緣象甲潛在地理分布區(qū)識(shí)別[J]. 昆蟲(chóng)學(xué)報(bào)， 2022，65（10）： 1334-1342.

LIANG L， XIAN X Q， ZHAO H X， GUO J Y， LIU W X.Estimation of the potential geographical distribution ofNaupactus leucoloma （Coleoptera： Curculionidae） based onthe MaxEnt model[J]. Acta Entomologica Sinica， 2022，65（10）： 1334-1342. （in Chinese）

[20] BURCKHARDT D， OUVRARD D， PERCY D M. An updatedclassification of the jumping plant-lice （Hemiptera：Psylloidea） integrating molecular and morphological evidence[J]. European Journal of Taxonomy， 2021， 736：137-182.

[21] VRIJENHOEK R. DNA primers for amplification of mitochondrialcytochrome c oxidase subunit I from diversemetazoan invertebrates[J]. Molecular Marine Biology andBiotechnology， 1994， 3（5）： 294-299.

[22] 齊國(guó)君， 陳婷， 高燕， 雷妍圓， 呂利華. 基于Maxent 的大洋臀紋粉蚧和南洋臀紋粉蚧在中國(guó)的適生區(qū)分析[J]. 環(huán)境昆蟲(chóng)學(xué)報(bào)， 2015， 37（2）： 219-223.

QI G J， CHEN T， GAO Y， LEI Y Y， LYU L H. Potentialgeographic distribution of Planococcus minor and P. lilacinusin China based on Maxent[J]. Joumal of EnrironmentalEntomology， 2015， 37（2）： 219-223. （in Chinese）

[23] QIN Z， ZHANG J E， DITOMMASO A， WANG R L， WU RS. Predicting invasions of Wedelia trilobata （L.） Hitchc. withMaxent and GARP models[J]. Journal of Plant Research，2015， 128（5）： 763-775.

[24] CRAWFORD D L. A monograph of the jumping plant-liceor Psyllidae of the New World[J]. Bulletin of the UnitedStates National Museum， 1914， 85： 186.

[25] HOLLIS D M. A review of the Malvales-feeding psyllidfamily Carsidaridae （Homoptera）[J]. Bulletin of the BritishMuseum of Natural History， 1987， 56： 87-127.

[26] 陳新美， 雷淵才， 張雄清， 賈宏炎. 樣本量對(duì)MaxEnt 模型預(yù)測(cè)物種分布精度和穩(wěn)定性的影響[J]. 林業(yè)科學(xué)， 2012，48（1）： 53-59.

CHEN X M， LEI Y C， ZHANG X Q， JIA H Y. Effects ofsample sizes on accuracy and stablity of maximum entropymodel in predicting species distribution[J]. Scientia SilvaeSinicae， 2012， 48（1）： 53-59. （in Chinese）

[27] MAUCK K E， GEBIOLA M， PERCY D M. The hidden secretsof Psylloidea： biology， behavior， symbionts， and ecology[J]. Annual Review of Entomology， 2024， 69： 277-302.

[28] OUVRARD D， CHALISE P， PERCY D M. Host-plant leapsversus host-plant shuffle： a global survey reveals contrastingpatterns in an oligophagous insect group （Hemiptera， Psylloidea）[J]. Systematics and Biodiversity， 2015， 13（5）： 434-454.

[29] ANUTRAKUNCHAI S， SUTTHIAPAPONG T， CHANSIRIP， THONGKAMNGAM T. Effect of indigenous Beauveriabassiana on the control of Allocaridala maleyensis （Crawford）in durian plantation areas in eastern region[J]. KingMongkut’s Agricultural Journal， 2023， 41（3）： 213-222.[30] 朱耿平， 劉國(guó)卿， 卜文俊， 高玉葆. 生態(tài)位模型的基本原理及其在生物多樣性保護(hù)中的應(yīng)用[J]. 生物多樣性， 2013，21（1）： 90-98.

ZHU G P， LIU G Q， BU W J， GAO Y B. Ecological niche modeling and its applications in biodiversity conservation[J].Biodiversity Science， 2013， 21（1）： 90-98. （in Chinese）

[31] 李響， 張成福， 賀帥， 王雨晴， 苗林. MaxEnt 模型綜合應(yīng)用研究進(jìn)展分析[J]. 綠色科技， 2020（14）： 14-17.

LI X， ZHANG C F， HE S， WANG Y Q， MIAO L. Researchprogress analysis on the comprehensive application of Max-Ent model[J]. Journal of Green Science and Technology，2020（14）： 14-17. （in Chinese）

[32] 孫雪婷. CMIP6 氣候變化情景下黑色枝小蠹全球潛在適生區(qū)分布預(yù)測(cè)[J]. 昆蟲(chóng)學(xué)報(bào)， 2023， 66（3）： 369-380.

SUN X T. Prediction of global potential suitable area of theblack twig borer， Xylosandrus compactus （Coleoptera： Scolytinae）under future climate change scenarios of CMIP6[J].Acta Entomologica Sinica， 2023， 66（3）： 369-380. （in Chinese）

[33] GADUG S， HUSSEIN M Y. The biology of the durian carsidaridTenaphalara malayensis Crawf. （Homoptera： Carsidaridae）[J]. Pertanika Journal of Science amp; Technology，1987， 10（2）： 243-246.

[34] TIGVATTNANONT S， PRAMUAL C. Economic importanceof the durio psyllid， Tenaphalara malayensis Crawfordin Thailand[J]. Kaen Kaset， 1990， 18： 152-159.

[35] WIANGSAMUT B， THONGKAMNGAM T， MAKHONPASC， AMLOY P， WIANGSAMUT M E L， ANUTRAKINCHAIS. The management practice methods for Allocarsidara malayensis（Crawford） by spraying insecticide with AI droneand long hose pump sprayers in durian orchards[J]. InternationalJournal of Agricultural Technology， 2024， 20（4）：1709-1728.

基金項(xiàng)目 海南省重點(diǎn)研發(fā)項(xiàng)目（No. ZDYF2023XDNY178，No. ZDYF2024XDNY281）。