朱國淵（1989-），副研究員，云南省科技特派員，云南省天然橡膠產(chǎn)業(yè) 技術(shù)體系病蟲害崗位核心成員，云南省天然橡膠可持續(xù)利用研究重點實 驗室（籌）成員，主要從事熱帶作物害蟲防控技術(shù)及天敵昆蟲資源發(fā)掘利 用研究工作。主持省部級項目2項，作為主要成員參與國家重點研發(fā)計 劃項目、省級科研項目10余項，在《Mitochondrial DNAPart B：Resources》 《Plant Disease》《南方農(nóng)業(yè)學(xué)報》等國內(nèi)外核心期刊上發(fā)表學(xué)術(shù)論文20余 篇；作為主要發(fā)明人獲授權(quán)發(fā)明專利5項、實用新型專利2項；作為主要起 草人制定云南省地方標(biāo)準(zhǔn)2項、西雙版納州地方標(biāo)準(zhǔn)3項；作為完成人獲 得云南省科學(xué)技術(shù)進(jìn)步獎三等獎1項，西雙版納州科學(xué)技術(shù)進(jìn)步獎三 等獎1項；參編技術(shù)手冊1本。

摘要：【目的】分析天敵昆蟲叉角厲蝽[Eocanthecona fircellata（Wolff）]不同發(fā)育階段腸道細(xì)菌群落組成及多樣 性，并預(yù)測其功能，為深入了解腸道細(xì)菌對叉角厲蝽生長發(fā)育和營養(yǎng)獲取方式的影響及開發(fā)新的害蟲生物防治策略 提供理論依據(jù)?！痉椒ā窟x取叉角厲蝽1齡、2齡、5齡若蟲和7日齡雌、雄成蟲，分別提取其腸道細(xì)菌基因組DNA。利 用 單分子實時測序技術(shù)對腸道細(xì)菌16S rRNA序列全長進(jìn)行測序，統(tǒng)計操作分類單元（OTU） 數(shù)量，分析腸道細(xì)菌群落組 成結(jié)構(gòu)與多樣性，并采用PICRUSt 預(yù)測腸道細(xì)菌功能?！窘Y(jié)果】從叉角厲蝽15個腸道樣品中共得到202401條有效序 列，聚類為222個OTUs，其中共有OTUs為5個。共注釋到19門32綱62目103科154屬188種細(xì)菌，優(yōu)勢門、綱、目、科、屬、種分別為變形菌門（Proteobacteria，72.22% ）、y-變形菌綱（Gammaproteobacteria，72.21%）、腸桿菌目（Enterobacterales， 72.11%）、腸桿菌科（Enterobacteriaceae，49.39%）、西地西菌屬（Cedecea，49.24%） 和奈氏西地西菌（C.neteri，49.24%）。多樣性和差異分析結(jié)果顯示，腸道細(xì)菌菌群組成在叉角厲蝽不同發(fā)育階段存在顯著差異（Plt;0.05， 下同）。線性判 別分析結(jié)果顯示，1齡、2齡若蟲和7日齡雌、雄成蟲腸道具有顯著優(yōu)勢的菌群，而5齡若蟲腸道不存在顯著優(yōu)勢菌群（Pgt;0.05）， 且不同發(fā)育階段腸道的優(yōu)勢菌群在科、屬、種3級分類階元水平上亦不同。KEGG代謝途徑分析發(fā)現(xiàn)，叉角 厲蝽不同發(fā)育階段腸道菌群的功能相似，主要涉及新陳代謝（74.1 7%～75.59%）和環(huán)境信息處理（10.53%～11.03%），二級代謝功能通路主要為全局和概述圖譜、糖類代謝和氨基酸代謝等。【結(jié)論】叉角厲蝽不同發(fā)育階段腸道細(xì)菌群落組 成和多樣性存在顯著差異，優(yōu)勢菌分布呈動態(tài)變化。

關(guān)鍵詞：叉角厲蝽；腸道微生物；細(xì)菌多樣性；單分子實時測序

中圖分類號：S476.2;Q965 文獻(xiàn)標(biāo)志碼：A 文章編號：2095-1191（2024）01-0013-11

Composition，diversity and function prediction of the gut bacterial community in different developmental stages of Eocanthecona furcellata（Wolff）

ZHU Guo-yuan'2，LIU Bing-shun3，ZHANG Man' ，SHI Chun-lan'，TANG Yong-yu1，LAN Ming-xian'，WU Guo-xing'，Ci-na-you-ruo'，GAO Xil°

（'Plant Protection College，Yunnan Agricultural University，Kunming，Yunnan 650201，China;2Yunnan Institute of Tropical Crops，Jinghong，Yunnan 666100，China;'Lufeng Agricultural Technical Extension Center， Chuxiong，Yunnan 651200，China）

Abstract ：[Objective]The purpose of the study was to analyze the composition and diversity of gut bacterial commu- nities in different developmental stages of the natural enemy insect Eocanthecona furcellata（Wolff），and to predict their functions，in order to provide a theoretical basis for a deeper understanding of the influenceof gut bacteria on the growth， development and nutrient acquisition patterns of E.furcellata，as well as for the development of a novel biocontrol strate- gies for the pests.【Method]The gut bacterial genomic DNA was extracted from the first instar nymphs，second instar nymphs，fifth instar nymphs，7-day-old female and male adults of E.furcellata，respectively.The full-length of 16S rRNA sequence of E.fucellata gut bacteria was sequenced using the single-molecule real-time sequencing technique.The number of operational taxonomic unit（OTU）was counted，the composition and diversity of gut bacterial communities were analyzed，and the functions of gut bacteria were predicted by PICRUSt.【Result]Atotal of 202401 effective se- quences were obtained from 15 gut samples ofE.fiurcellata，and clustered into 222 OTUs，of which 5OTUs were shared by them.A total of 19 phyla，32 classes，62 orders，103 families，154 genera and 188 species bacteria were annotated.The dominant phylum，class，order，family，genus and species of the bacteria wereProteobacteria（72.22%），Gammaproteo- bacteria（72.21%），Enterobacterales（72.11%），Enterobacteriaceae（49.39%），Cedecea（49.24%），and C.neteri（49.24%）， respectively.The results of the diversity and variance analysis revealed significant differences in the composition of the gut bacterial flora in different developmental stages of E.fuacellata（Plt;0.05，thesame below）.The linear discriminant analysis results showed that significantly dominant flora existed in the guts of first instar nymphs，second instar nymphs and 7-day-old female and male adults，while no significantly dominant flora existed in the gut of fifth instar nymphs（P gt;0.05），and the dominant flora in the gut ofdifferent developmental stages varied at the level of 3 taxonomic catego- ries，namely，family，genus and species.KEGG metabolic pathway analysis revealed that the functions of the gut flora of E.fiurcellata were similar in different developmental stages，mainly involving metabolism（74. 17%-75.59%）and environ- mental information processing（10.53%-11.03%）.The secondary metabolic pathways were mainly global and overview mapping，sugar metabolism，and amino acid metabolism，etc.【Conclusion]Significant differences in the composition and diversity of gut bacterial community are observed in different developmental stages of E.fiucellata，and the distribu- tion of dominant bacteria shows dynamic changes.

Key words：Eocanthecona furcellata（Wolff）;gut microorganism;bacterial diversity;single-molecule real-time se-quencing

Foundation items：Yunnan Agricultural Basic Research Joint Special Project（202101BD070001-050，202301BD 070001-137）;Reserve Talents Project for Young and Middle-aged Academic and Technical Leaders of Yunnan（202205AC160077）;Science and Technology Plan Project of Yunnan Tobacco Company（2023530000241012）

0 引言

【研究意義】叉角厲蝽[Eocanthecona furcellata（Wolff）] 隸屬于半翅 目（Hemiptera） 蝽 科（Pentatomi- dae） 益蝽亞科（Asopinae） 厲蝽屬（Eocanthecona）， 可捕食草地貪夜蛾（Spodoptera frugiperda）、斜紋夜蛾（S.litura） 和柑桔尺蛾（Hyposidra talaca）等40多種 農(nóng)林害蟲的幼蟲，是一種在生物防治方面具有重要作用和應(yīng)用潛力的天敵昆蟲（范悅莉等，2019;宮靖 垚等，2019;Keerthi et al.，2020;Sarkar et al.，2021）。叉角厲蝽1齡若蟲僅刺吸嫩枝梢、嫩葉汁液，2～5齡 若蟲及成蟲均能捕食害蟲，若蟲齡期越大捕食能力 越強（姚明勇，2019;張曼等，2022）。昆蟲腸道作為 一種特殊生境，棲息于此的大量微生物對昆蟲的物 種形成、營養(yǎng)代謝、生殖代謝及免疫等具有重要意義（王四寶和曲爽，2017;張振宇等，2017;Song et al.，2022）。同時，昆蟲腸道微生物物種多樣性和特異性 受宿主的系統(tǒng)發(fā)育、不同發(fā)育時期、食物和營養(yǎng)等的 影響（Colman et al.，2012;Yun et al.，2014）。鑒 于 此，利用高通量測序技術(shù)研究昆蟲腸道微生物的結(jié) 構(gòu)組成和多樣性，已成為研究熱點（Fernandes et al.， 2014）。然而截至目前，關(guān)于叉角厲蝽腸道細(xì)菌與其 不同發(fā)育階段間的關(guān)聯(lián)鮮有報道。因此，研究叉角 厲蝽不同發(fā)育階段腸道細(xì)菌的組成、多樣性 及顯著 差異菌群，并分析腸道細(xì)菌對其生長發(fā)育的影響，對 該天敵昆蟲人工飼料中微生物添加劑的開發(fā)與應(yīng)用 以及大規(guī)模人工擴繁和利用具有重要意義?！厩叭搜?究進(jìn)展】目前關(guān)于昆蟲腸道微生物的研究主要涉 及鱗翅目（Lepidoptera）、雙翅目（Diptera）、膜翅目（Hymenoptera）、蜚蛾目（Blattaria）、鞘翅目（Coleop- tera）、直翅目（Orthoptera） 和半翅目等昆蟲（鄭林宇 等，2022），其中對半翅目蝽類的大多數(shù)研究集中于 蝽次目（Infraorder Pentatomomorpha）的植食性 蝽及 與異翅亞目（Heteropteran） 昆蟲譜系相關(guān)的腸道微 生物區(qū)系（Shan et al.，2021;Lietal.，2022），如斑腹刺益蝽（Podisus maculiventris）、中黑盲蝽（Adelphoco-ris suturalis）、荔枝蝽（Tessaratoma papillosa）和綠盲 蝽（Apolygus lucorum）等（Ahsaei et al.，2020;Hu et al.，2020;Xue et al.，2021;Liu et al.，2022;Xue et al.， 2023） 。Xue 等（2021）基于Illumina HiSeq高通量測 序平臺分析了中黑盲蝽1～5齡若蟲和1、6、9日齡雌、雄成蟲腸道細(xì)菌組成及多樣性，發(fā)現(xiàn)上述發(fā)育階段 中黑盲蝽腸道內(nèi)的葡萄球菌屬（Staphylococcus）、不 動桿菌屬（Acinetobacter）、假單胞菌屬（Pseudomo- nas）和棒桿菌屬（Corynebacterium）等微生物的相對 豐度呈顯著動態(tài)變化。安哲等（2022）通過相似方法 研究了棉田綠盲蝽若蟲和雌、雄成蟲體內(nèi)共生菌的 群落組成與多樣性，結(jié)果顯示，綠盲蝽不同發(fā)育時期 體內(nèi)共生菌的相對豐度和多樣性指數(shù)不同，且成蟲 期的物種豐富度與群落多樣性均大于若蟲期。韓順 財?shù)龋?022）亦采用該方法分析了采自云南、福建、廣 東、廣西、海南和臺灣的荔枝蝽成蟲和若蟲體內(nèi)細(xì)菌 組成及多樣性，發(fā)現(xiàn)6?。ㄗ灾螀^(qū)）荔枝蝽成蟲間細(xì) 菌多樣性差異不顯著，但1齡若蟲間細(xì)菌組成及其 占比差異顯著;成蟲體內(nèi)優(yōu)勢菌為腸桿菌，而多數(shù) 1齡若蟲體內(nèi)優(yōu)勢菌為鏈霉菌。Kuang 等（2022）利用高通量測序技術(shù)分析3種飼料飼養(yǎng)的叉角厲蝽雌 成蟲中腸腸道細(xì)菌 16S rRNA 序列V3～V4區(qū)域，結(jié) 果發(fā)現(xiàn)變形菌門（Proteobacteria） 和厚壁菌門（Fir-micutes）為優(yōu)勢菌門，沙雷氏菌屬（Serratia）（人工飼料處理的相對豐度較低）、腸球菌屬（Enterococcus）和腸桿菌科（Enterobacteriaceae）的一個未培養(yǎng)細(xì)菌 屬為優(yōu)勢菌屬，且不同飼料飼養(yǎng)的叉角厲蝽雌成蟲 中腸腸道微生物區(qū)系組成存在顯著差異。Liu 等 （2022）利用依賴培養(yǎng)（Culture-dependent）和第三代測 序技術(shù)評估了荔枝蝽腸道細(xì)菌群落多樣性和動態(tài)變 化，結(jié)果表明，7個發(fā)育階段中存在5個細(xì)菌門類，變 形菌門為優(yōu)勢菌門，泛菌屬（Pantoea）為優(yōu)勢菌屬； 卵階段的細(xì)菌群落最復(fù)雜；細(xì)菌群落在卵期、1～3齡 若蟲期、4～5齡若蟲期和成蟲期發(fā)揮著不同的作用?！颈狙芯壳腥朦c】目前，有關(guān)叉角厲蝽的研究報道主 要集中在生物學(xué)特性、捕食功能反應(yīng)、控害能力、觸 角轉(zhuǎn)錄組及嗅覺相關(guān)基因分析、唾液腺轉(zhuǎn)錄組及差 異分析、實時熒光定量PCR內(nèi)參基因篩選與驗證、不同食物對雌成蟲中腸細(xì)菌群落及多樣性影響等 方面（高平等，2021;張曼等，2022;趙航等，2022; Kuang et al.，2022;Pan et al.，2022），尚未涉及其不同 齡期若蟲和成蟲腸道細(xì)菌群落組成及多樣性?！緮M解 決的關(guān)鍵問題】基于PacBio測序平臺，使用單分子實 時（Single-molecule real-time，SMRT）測序技術(shù)對叉 角厲蝽1齡、2齡、5齡若蟲和7日齡雌、雄成蟲腸道 細(xì)菌16S rRNA序列全長進(jìn)行測序，并采用PICRUSt 預(yù)測其腸道細(xì)菌功能，旨在探明該蟲不同發(fā)育階段 腸道細(xì)菌群落組成及多樣性差異，為深入了解腸道 細(xì)菌對叉角厲蝽生長發(fā)育和營養(yǎng)獲取方式的影響及 開發(fā)新的害蟲生物防治策略提供理論依據(jù)。

1 材料與方法

1.1 供試?yán)ハx

叉角厲蝽采自云南省玉溪市元江縣濱江大道小 燕村。在云南省昆明市北郊云南農(nóng)業(yè)大學(xué)植物保護(hù) 學(xué)院養(yǎng)蟲室內(nèi)用黃粉蟲蛹飼養(yǎng)，建立穩(wěn)定種群。養(yǎng) 蟲室內(nèi)溫度（26±1）℃、相對濕度（65±10）%、光周期 14L：10D。

1.2 叉角厲蝽腸道解剖

根據(jù)食性、捕食量及產(chǎn)卵前期，分別挑選健康、生長發(fā)育一致且已饑餓處理24 h 的叉角厲蝽1齡、2齡、5齡若蟲和7日齡雌、雄成蟲，解剖腸道。每組 3次重復(fù)，每重復(fù)20頭。解剖前將試蟲于-20℃冰 箱 中冷凍3～5 min， 使其麻痹后置于超凈工作臺中，先 用75%酒精進(jìn)行蟲體體表消毒20 s， 再用無菌水清 洗3～5次，后轉(zhuǎn)至盛有無菌磷酸鹽緩沖液（PBS， pH 7.2～7.5）的培養(yǎng)皿中，置于體視顯微鏡下解剖。

為便于解剖叉角厲蝽腸道，可先剪去試蟲胸足 和成蟲雙翅。解剖時，將試蟲正放，先用普通鑷子夾 住其頭、胸部，后用尖頭鑷子輕輕刺入背部，向后慢慢拉開表皮層，露出內(nèi)容物，再用鑷子夾住嗦囊處， 向后拉出食道，剔除唾液腺、馬氏管、肛門等組織，取 出完整腸道。用無菌水沖洗后，置于1.5 mL離心管 中，蓋上蓋子并用封口膜包封，將其置于液氮中速凍 3 h后取出，轉(zhuǎn)至-80℃冰箱保存?zhèn)溆谩?/p>

1.3 叉角厲蝽腸道細(xì)菌基因組DNA提取

采用磁珠法并按基因組DNA提取試劑盒 [DP812，天根生化科技（北京）有限公司]說明書操 作步驟提取腸道細(xì)菌基因組DNA。用NanoDrop 2000型超微量分光光度計（賽默飛世爾科技有限公 司）檢測基因組DNA濃度和純度。

1.416S rRNA序列PCR擴增選取27F（5'-AGRGTTTGATYNTGGCTCAG-3'） 和1492R（5'-TASGGHTACCTTGTTASGACTT-3'）引 物，用Veriti 96-Well型梯度基因擴增儀（Applied Biosystems公司）擴增16SrRNA 序列全長，PCR反 應(yīng)體系30.0 μL：KOD OnePCR Master Mix 15.0μL， DNA模板1.5 μL，正、反向引物各1.5 μL， 無核酸酶 水10.5μL。擴增程序：95℃預(yù)變性2 min;98℃10s， 55℃30s，72℃1 min 30s，進(jìn)行25個循環(huán)；72℃延 伸2 min。

經(jīng)Qubit 4.0 Q33216型熒光定量儀（賽默飛世爾 科技有限公司）和1.8%瓊脂糖凝膠電泳分別對擴增 產(chǎn)物的濃度和條帶進(jìn)行檢測后，送至北京百邁客 生物科技有限公司使用三代測序平臺PacBic Sequel Ⅱ進(jìn)行測序。

1.5 統(tǒng)計分析

對 PacBio原始下機 subreads 進(jìn)行校正，得到 cCS（Circular consensus sequencing）序列；利用lima v1.7.0， 通過barcode對CCS序列進(jìn)行識別，得到原始 序列數(shù)（Raw CCS）;利用cutadapt v2.7識別正、反向 引物，丟棄不包含引物的CCS序列，對CCS長度進(jìn) 行過濾，丟棄不滿足長度閾值（1200～1650 bp） 的 序列，得到去除引物后序列（Clean CCS）;利用 UCHIME v8.1去除嵌合體序列，得到最終有效序 列數(shù)（Effective CCS）。使用USEARCH v10.0（http：// drive5.com/uparse/）（Edgar，2013） 在相似性97%的 水平上對序列進(jìn)行聚類并得到操作分類單元（Opera- tional taxonomic unit，OTU）。

使用Ribosomal Database Project Classifier的樸 素貝葉斯算法對97%相似水平的OTU 代表序列進(jìn) 行分類學(xué)分析（Wanget al.，2007），與Silva v138數(shù)據(jù) 庫（https：//www.arb-silva.de/）進(jìn)行比對，在各分類階 元上統(tǒng)計物種組成。使用QIIME 2分析菌群豐富度（ACE 、Chao指數(shù)）和多樣性（Simpson、Shannon指數(shù) ）（Bolyen et al.，2019）。利 用Metastats（White et al.，2009） 分析組間豐度差異顯著的菌群?；贐ray- Curtis距離算法，用非度量多維標(biāo)度（Non-metric multidimensional scaling，NMDS）分析比較不同組 間群落組成的相似性，并用置換多元方差分析（Permutational multivariate analysis of variance， PERMANOVA）檢驗差異顯著性。

通過PICRUSt2 v2.3.0（Langille et al.，2013）預(yù) 測叉角厲蝽不同發(fā)育階段腸道菌群KEGG代謝通 路，使用STAMP v2.0.8（Parks et al.，2014）分析樣品 間腸道菌群功能的差異。利用R 語言vegan v2.3 作圖。

2 結(jié)果與分析

2.1 叉角厲蝽不同發(fā)育階段腸道細(xì)菌測序結(jié)果

基于PacBio 測序平臺，采用SMRT測序技術(shù)分 析叉角厲蝽1齡、2齡、5齡若蟲和7日齡雌、雄成蟲 腸道細(xì)菌16S rRNA序列全長，結(jié)果從15個樣品中 共獲得203678條Raw CCS，平均每樣品13579條 CCS; 經(jīng)質(zhì)量過濾、嵌合體去除，得到202401條 Effective CCS，平均長度為1469.2 bp（表1）。基于 Effective CCS進(jìn)行聚類分析，15個樣品共得到222個 OTUs， 其中共有OTUs為5個，說明各樣品細(xì)菌群 落組成存在差異。Shannon指數(shù)曲線顯示，隨著測 序數(shù)量的增加，曲線趨于平緩（圖1），說明本研究測序數(shù)據(jù)量充足。此外，所有樣品的覆蓋度指數(shù)均在 0.9970以上，表明具備足夠的測序深度和可靠的數(shù) 據(jù)質(zhì)量。

2.2 叉角厲蝽腸道主要細(xì)菌群落組成

對特征序列進(jìn)行分類學(xué)注釋，共得到19門32綱 62目103科154屬188種細(xì)菌。在門分類水平上，主 要菌群有變形菌門（72.22%）、厚壁菌門（27.52%）、擬桿菌門（Bacteroidetes，0.13%）、疣微菌門（Verruco- microbiota，0.04%） 和放線菌門（Actinobacteriota， 0.02%）;在綱分類水平上，主要有γ-變形菌綱（Gam- maproteobacteria，72.21%）、芽孢桿菌綱（Bacilli， 27.41%）、擬桿菌綱（Bacteroidia，0.13%）、梭菌綱（Clostridia，0.11%）、疣微菌綱（Verrucomicrobiae， 0.04%）、阿爾法變形菌綱（Alphaproteobacteria，0.02%） 和放線桿菌綱（Actinobacteria，0.02%）等；在目分類 水平上，主要有腸桿菌目（Enterobacterales，72.11%）、乳桿菌目（Lactobacillales，27.32%）、擬桿菌目（Bacte- roidales，0.11%）、伯克氏菌目（Burkholderiales，0.06%）、梭菌目（Clostridiales，0.05%）、蟲原體目（Entomo- plasmatales，0.05%）、假單胞菌目（Pseudomonadales， 0.04%）和疣微菌目（Verrucomicrobiales，0.04%）等； 在科分類水平上，主要有腸桿菌科（49.39%）、腸球菌 科（Enterococcaceae，26.60%）、耶爾森菌科（Yersi- niaceae，19.20%） 、摩根菌科（Morganellaceae，3.49%）、乳桿菌科（Lactobacillaceae，0.40%）、鏈球菌科（Strep- tococcaceae，0.31%）、鼠桿狀菌科（Muribaculaceae， 0.06%）、螺原體科（Spiroplasmataceae，0.05%）、梭菌 科（Clostridiaceae，0.05%）和阿克曼氏菌科（Akker- mansiaceae，0.03%） 等；在屬分類水平上，主要有西地 西菌屬（Cedecea，49.24%） 、腸球菌屬（26.60%）、沙雷氏菌屬（19.20%）、變形桿菌屬（Proteus，2.24%）、摩根氏菌屬（Morganella，1.24%）、乳球菌屬（Lacto-coccus，0.19%）、慢生乳桿菌屬（Lentilactobacillus， 0.17%）、埃希氏一志賀氏菌屬（Escherichia-Shigella， 0.15%）、鏈球菌屬（Streptococcus，0.12%）和乳植桿菌屬（Lactiplantibacillus，0.10%）等；在種分類水平上，主要有奈氏西地西菌（Cedecea neteri，49.24%）、糞腸球菌（Enterococcus faecalis，25.03%） 、Serratia endosymbiont of Nilaparvata lugens（19.20%） 、奇 異 變形桿菌（Proteus mirabilis，2.24%）、蒼黃腸球菌（Enterococcus gilvus，1.56%）、摩氏摩根氏菌（Mor-ganella morganii，1.24%）、乳酸乳球菌（Lactococcus lactis，0. 18%）、布赫納氏慢生乳桿菌（Lentilactobacil- lus buchneri，0.17%）、大腸埃希氏菌（Escherichia coli，0.15%）和植物乳植桿菌（Lactiplantibacillus plantarum，0.10%）等（圖2）。

2.3 叉角厲蝽不同發(fā)育階段腸道細(xì)菌多樣性分析結(jié)果

2.3.1 腸道細(xì)菌Alpha 多樣性分析 由圖3可知，叉 角厲蝽從低齡若蟲到高齡若蟲及成蟲的發(fā)育過程 中，其腸道細(xì)菌豐富度指數(shù)（ACE指數(shù)和Chaol指 數(shù)）減小，多樣性指數(shù)（Simpson指數(shù)和 Shannon指 數(shù)）逐漸增大。ACE 指數(shù)在1齡若蟲與雌成蟲之間 存在顯著差異（Plt;0.05， 下同）。Chaol 指數(shù)在1齡 若蟲和其他發(fā)育階段之間均存在顯著差異，即1齡 若蟲的Chaol 指數(shù)顯著高于2齡和成蟲，極顯著高于5齡若蟲（Plt;0.01， 下同）。1齡若蟲的Simpson 指 數(shù) 顯著低于2齡和雌成蟲，極顯著低于5齡若蟲；5齡 若蟲的Simpson指數(shù)顯著高于雄成蟲；雌、雄成蟲間 Simpson 指數(shù)差異顯著。雌成蟲的Shannon 指數(shù)顯著高于1齡和2齡若蟲，極顯著高于雄成蟲；5齡若蟲的 Shannon指數(shù)顯著高于2齡若蟲和雄成蟲。

2.3.2 腸道細(xì)菌Beta多樣性分析 NMDS分析結(jié) 果（圖4）顯示，不同樣品按照不同發(fā)育階段聚集，且 1齡若蟲與2齡、5齡若蟲及雌、雄成蟲距離存在明顯 分離，Stress=0.0448（lt;0.05），說明NMDS分析具有 很好的代表性。PERMANOVA分析結(jié)果（圖5）顯 示，叉角厲蝽不同發(fā)育階段腸道細(xì)菌菌群結(jié)構(gòu)差異 極 顯 著（F=16.627，R2=0.869，P=0.001）。

2.4 叉角厲蝽不同發(fā)育階段腸道顯著差異菌群 分析結(jié)果

采 用LEFSe（Linear discriminant analysis effect size）分析叉角厲蝽不同發(fā)育階段腸道菌群差異（LDAgt;4.0），結(jié)果（圖6）顯示，1齡、2齡若蟲和雌、雄 成蟲腸道存在顯著優(yōu)勢的菌群，而5齡若蟲腸道不 存在顯著優(yōu)勢的菌群（Pgt;0.05，下同），且不同發(fā)育 階段腸道的優(yōu)勢菌群在科、屬、種3級分類水平上亦 不 同 。在 科 分 類 水 平 ， 腸桿 菌 科 顯 著 存 在 于 1 齡若蟲階段，耶爾森菌科顯著存在于2齡若蟲階段，摩 根菌科顯著存在于雌成蟲階段。在屬分類水平，西 地西菌屬顯著存在于1齡若蟲階段，沙雷氏菌屬顯 著存在于2齡若蟲階段，摩根菌屬顯著存在于雄成 蟲階段，變形桿菌屬顯著存在于雌成蟲階段。在種 分類水平，糞腸球菌、奈氏西地西菌、奇異變形桿 菌 、蒼黃腸球菌、Serratia endosymbiont of Nilaparvata lugens 5個菌為所有發(fā)育階段腸道共有菌；奈氏西地 西菌顯著存在于1齡若蟲階段，糞腸球菌和Serratia endosymbiont of Nilaparvatalugens顯著存在2齡若蟲階段，蒼黃腸球菌和奇異變形桿菌顯著存在于雌 成蟲階段，摩氏摩根氏菌顯著存在于雄成蟲階段。

2.5 叉角厲蝽不同發(fā)育階段腸道菌群功能基因 預(yù)測結(jié)果

KEGG代謝途徑分析結(jié)果顯示，叉角厲蝽不同 發(fā)育階段腸道菌群的功能相似，主要涉及新陳代謝（74.17%～75.59%）、環(huán)境信息處理（10.53%～11.03%）、遺傳信息處理（5.71%～6.11%）、細(xì)胞過程（3.78%～4. 3 2 % ） 、人 類 疾 ?。? . 0 7 % ～ 3 . 3 5 % ） 和 有 機 體 系 統(tǒng)（1.32%～1.39%）（圖7-A） 。進(jìn) 一 步利用單因素方差 分析發(fā)現(xiàn)，不同發(fā)育階段腸道新陳代謝通路在功能 豐度上存在顯著差異，即2齡若蟲的代謝功能顯著 高于5齡若蟲和雌成蟲，極顯著高于雄成蟲和1齡若 蟲（圖7-B）。

二級水平上功能通路預(yù)測結(jié)果有43種，其中二""""" 級代謝功能通路12種，依次為全局和概述圖譜（38.87%～39.29%）、糖類代謝（9.82%～10.86%）、氨基 酸代謝（6.29%～6.76%）、能量代謝（3.57%～3.86%）、輔助因子和維生素代謝（3.69%～3.81%）、核苷酸代 謝（3.13%～3.57%）、脂質(zhì)代謝（1.85%～2.37%）、其他氨基酸代謝（1.71%～1.88%）、外源化學(xué)物質(zhì)生物 降解與代謝（1.50%～1.63%）、聚糖生物合成與代謝（1.26%～1.34%）、萜類和聚酮類化合物代謝（0.87%～ 1 . 0 2 % ） 及 其 他 次 生 代 謝 產(chǎn) 物 生 物 合 成（0.70%～ 0.74%）。

3 討論

叉角厲蝽具有繁殖力強、捕食范圍廣和捕食能 力強等優(yōu)勢，已成為生物防治領(lǐng)域應(yīng)用前景廣闊的 一種天敵昆蟲。本研究基于PacBio測序平臺，利用 SMRT 測序技術(shù)分析了叉角厲蝽不同發(fā)育階段的腸 道細(xì)菌多樣性及其組成，并探索了該蟲與其腸道細(xì) 菌的協(xié)同進(jìn)化關(guān)系，結(jié)果共獲得222個用于物種分 類的OTUs，涵蓋了19門32綱62目103科154屬188種，其中優(yōu)勢菌門為變形菌門和厚壁菌門，優(yōu)勢菌屬 為西地西菌屬、腸球菌屬和沙雷氏菌屬。本研究結(jié) 果與Kuang等（2022）對叉角厲蝽雌成蟲中腸細(xì)菌組 成的研究結(jié)論相似，但優(yōu)勢菌群的相對豐度和優(yōu)勢 種有一定差異，可能與飼養(yǎng)代數(shù)、雌成蟲日齡、取樣 部位以及取樣方式等有關(guān)。Xue 等（2023）報道綠盲 蝽6和9日齡雌成蟲的優(yōu)勢屬為乳球菌屬，而雌成蟲 的變形菌門相對豐度從1日齡開始顯著下降且9日 齡最低。本研究中叉角厲蝽雌成蟲腸道細(xì)菌的樣本 取自7日齡已交配且饑餓處理24 h的整個腸道，而 交配及饑餓處理對其腸道細(xì)菌組成及豐富度的影響 機制有待進(jìn)一步研究。

Xue等（2021）報道中黑盲蝽1齡和2齡若蟲腸 道細(xì)菌群落多樣性指數(shù)最高，而從3齡若蟲開始下降。綠盲蝽1齡和2齡若蟲體內(nèi)細(xì)菌群落多樣性和 豐富度亦最高，且顯著高于其他發(fā)育階段（Xue et al.，2023）; 但安哲等（2022）發(fā)現(xiàn)綠盲蝽若蟲期體內(nèi) 微生物的物種豐富度和群落多樣性小于成蟲期。本 研究發(fā)現(xiàn)叉角厲蝽從1齡若蟲到成蟲性成熟的生長 發(fā)育過程中，其腸道細(xì)菌豐富度指數(shù)逐漸降低，而其 多樣性指數(shù)逐漸升高，與前人研究結(jié)果（Luo et al.， 2021;Li et al.，2022）存在一定差異，可能與天敵昆蟲 體內(nèi)微生物群落受其取食習(xí)性和食物類型的影響有 關(guān)。中黑盲蝽和綠盲蝽具有相近的習(xí)性和生態(tài)功 能，不僅具有植食性，還兼具肉食性（Liet al.，2020）， 但相比于捕食獵物，它們更偏愛取食植物（李文靜等，2015），而叉角厲蝽除1齡若蟲僅刺吸植物汁液， 2齡若蟲處于植食和捕食的過渡期，其他齡期的若蟲 和成蟲均為捕食性。據(jù)報道，腸道微生物菌群對食 物變化引起的腸道環(huán)境變化需要代謝適應(yīng)過程（Tilocca et al.，2017）。因此，本研究中叉角厲蝽腸 道菌群的變化亦可能是為了逐步適應(yīng)新的腸道 環(huán)境。

本研究發(fā)現(xiàn)叉角厲蝽不同發(fā)育階段腸道優(yōu)勢菌 群不同，而某些菌僅存在于特定的發(fā)育時期，如奈氏 西地西菌主要存在于1齡若蟲階段、摩氏摩根氏菌 存在于5齡若蟲和雌、雄成蟲中。這些細(xì)菌的存在 和顯著變化可能與叉角厲蝽的生理代謝有關(guān)，從而 幫助不同發(fā)育階段的叉角厲蝽適應(yīng)不同的食物源和 生態(tài)環(huán)境（Engel and Moran，2013）。已有研究表明， 變形菌門細(xì)菌與宿主的捕食、生長發(fā)育和適合度等 有關(guān)（Nishiwaki et al.，2007;Kuang et al.，2022）。沙 雷氏菌屬能為宿主提供維生素和氨基酸，幫助其降 解纖維素、單萜或雙萜等物質(zhì)（Boone et al.，2013;胡霞等，2018），摩氏摩根氏菌在一定程度上能促進(jìn)獵 物死亡（Nishiwaki et al.，2007），變形桿菌屬有助于 吸引雌成蟲產(chǎn)卵（Chaudhury et al.，2010），而西地西 菌屬可分泌有機酸、降低環(huán)境pH 。在其他昆蟲的研 究中也有類似發(fā)現(xiàn)，如雜食或肉食中黑盲蝽的腸道 菌群中沙雷氏菌 屬、西地西菌屬和泛菌屬的豐度 較高，而植食中黑盲蝽中其他屬的豐度則較高（成 艷霞，2022）。與 植食性蝽相比，獵蝽科有一個特 有的腸道微生物區(qū)系，主要由腸球菌屬和不同種 類的變形桿菌組成（Li et al.，2022）。中華通草蛉（Chrysoperla sinica）卵、初孵幼蟲和成蟲都有相似 的微生物，其中假單胞菌屬在1齡幼蟲（初孵和1日 齡）、2齡幼蟲、3齡幼蟲和蛹期出現(xiàn)，在成蟲期未出 現(xiàn)；檸檬酸桿菌屬（Citrobacter） 、醋酸桿菌屬（Aceto-bacter）、朝井桿菌（Asaia sp.）和Rosenbergiella sp.僅見于雄蟲，而Cosenzaea myxofaciens僅見于雌蟲（趙 輝，2018;Zhao et al.，2019）。

腸道細(xì)菌對宿主代謝和健康具有重要的調(diào)節(jié)作 用（李青晏等，2020）。本研究通過對叉角厲蝽不同 發(fā)育階段腸道菌群功能注釋，發(fā)現(xiàn)新陳代謝為主要 功能且在不同發(fā)育階段存在顯著差異，其中2齡 若蟲的代謝功能豐度最高，5齡若蟲次之，可能歸因 于叉角厲蝽2齡若蟲處于食性轉(zhuǎn)換期且2齡和5齡 若蟲均處于不斷進(jìn)食的生長期，新陳代謝功能發(fā)揮 主要作用，則有助于該蟲利用食物營養(yǎng)并促進(jìn)其生 長發(fā)育。然而，基因功能預(yù)測并不能完全反映叉角 厲蝽腸道菌群的實際功能，預(yù)測到的這些功能有待 進(jìn)一步解析和驗證。

本研究雖然明確了叉角厲蝽不同發(fā)育階段腸道 細(xì)菌組成和多樣性，但不同細(xì)菌譜系在宿主間的傳 播方式和對宿主食物差異的響應(yīng)機制及其對宿主適 合度的影響等方面均有待深入研究。

4 結(jié)論

叉角厲蝽腸道菌群組成和多樣性在其不同發(fā)育 階段均存在顯著差異，優(yōu)勢菌分布呈動態(tài)變化。腸 道菌群中新陳代謝功能的豐度最高，二級代謝功能 通路主要為全局和概述圖譜、糖類代謝和氨基酸代 謝等。

參考文獻(xiàn)：

安哲，薛匯，高雪珂，朱香鎮(zhèn)，雒 瑜，崔金杰，馬德英.2022. 棉田綠盲蝽微生物多樣性研究[J]. 環(huán)境昆蟲學(xué)報，44（5）：1063-1070.[AnZ，Xue H，Gao XK，Zhu XZ，Luo J Y，Cui JJ，Ma DY.2022.Microbial diversity of Apolygus lucorumin cotton field[J].Journal of Environmental Ento- mology，44（5）：1063-1070.]doi：10.3969/j.issn.1674-0858 2022.05.2.

成艷霞.2022.食物差異導(dǎo)致中黑盲蝽腸道菌群紊亂并致死 宿主昆蟲的機制研究[D]. 武漢：湖北大學(xué).[Chen YX 2022.Variation of gut microbiota caused by an imbalance diet is detrimental to bugs survival[D].Wuhan：Hubei University.]

范悅莉，谷星慧，洗繼東，陳雪梅，陸永躍，張立猛，周文兵. 2019.叉角厲蝽對草地貪夜蛾的捕食功能反應(yīng)[J].環(huán)境 昆蟲學(xué)報，41（6）：1175-1180.[Fan Y L，Gu XH，Xian J"D，ChenXM，Lu YY，Zhang LM，Zhou WB.2019.Func-tional response of Eocanthecona furcellata（Hemiptera： Pentatomidae）to Spodoptera frugiperda（Lepidoptera："Noctuidae）[J].Journal of Environmental Entomology，41（6）：1175-1180.]doi：10.3969/j.issn.1674-0858.2019.06.06.

高平，廖賢斌，李麗芳，蘭明先，趙航，陳斌，吳國星，高熹. 2021.叉角厲蝽唾液腺轉(zhuǎn)錄組及差異分析[J]. 云南農(nóng)業(yè)大學(xué)學(xué)報（自然科學(xué)），36（1）：29-38.[GaoP，Liao XB， LiLF，Lan MX，Zhao H，Chen B，Wu GX，Gao X.2021. Transcriptome and difference analysis of Eocanthecona"furcellata（Wolff）[J].Journal of Yunnan Agricultural Uni- versity（Natural Science），36（1）：29-38.]doi：10.12101/j.issn.1004-390X（n）.201911011.

宮靖垚，陳科偉，溫健，劉杰，朱永基.2019.人工飼料飼養(yǎng)的 叉角厲蝽對昆蟲的捕食能力評價[J]. 環(huán)境昆蟲學(xué)報，41（3）：47 1-478.[Gong JY，Chen KW，Wen J，Liu J，Zhu Y"J.2019.Predatory capacity of Eocanthecona fucellate"（Wolff）reared with artificial diets[J].Journal of Environ mental Entomology，41（3）：471-478.]doi：10.3969/j.issn. 1674-0858.2019.03.6.

韓順財，蔣敬濤，程林，劉家生，李彤，彭凌飛.2022.中國六省 區(qū)荔枝蝽體內(nèi)細(xì)菌組成及多樣性分析[J]. 植物保護(hù)學(xué) 報，49（3）：749-757.[Han SC，Jiang JT，Cheng L，Liu J S，Li T，Peng LF.2022.Bacterial composition and diver-sity in the litchi stink bug，Tessaratoma papillosa（Hemip-tera：Tessaratomidae）from six provincesof China[J].Jour- nal of Plant Protection，49（3）：749-757.]doi：10.13802j. cnki.zwbhxb.2022.2020233.

胡霞，傅慧靜，李俊楠，林中平，張飛萍.2018.松墨天牛幼蟲 腸道纖維素降解細(xì)菌的分離與鑒定[J]. 福建農(nóng)林大學(xué)學(xué) 報（自然科學(xué)版），47（3）：322-328.[Hu X，F(xiàn)uHJ，LiJN， Lin Z P，Zhang F P.2018.lsolation and identification of cellulolytic bacteria associated with the gut of Monocha- mus alternatus larvae[J].Journal of Fujian Agriculture and Forestry University（Natural Science Edition），47（3）：322- 328.]doi：10.13323/ j.cnki.j.fafu（nat.sci.）.2018.03.009.

李青晏，唐運林，蔣睿軒，張永紅，朱峰，白興榮，顧偌鉞，吳燕 燕，吳玉嬌，陳潔，李田，李春峰，韋俊宏，潘國慶，周澤 揚.2020.云南地區(qū)草地貪夜蛾腸道細(xì)菌的分離及鑒定 [J].西南大學(xué)學(xué)報（自然科學(xué)版），42（1）：1-8.[LiQY，

Tan YL，Jiang RX，Zhang YH，Zhu F，BaiXR，Gu RC WuYY，Wu YJ，Chen J，LiT，LiCF，Wei JH，Pan GQ， Zhou ZY.2020.Isolation and identification of gut bacteria"of Spodoptera frugiperda feeding on maize in Yunnan"China[J].Journal of Southwest University（Natural Scien- ce），42（1 ）： 1-8.]doi：10.13718/j.cnki.xdzk.2020.01.00 1.

李文靜，苑鶴，陸宴輝，李強，吳孔明.2015.中黑盲蝽對四季 豆豆莢和棉鈴蟲卵的取食選擇性[J]. 植物保護(hù)，41（4）：29-34.[LiWJ，Yuan H，Lu Y H，LiQ，Wu K M.2015.

Diet selection of Adelphocoris suturalis on Phaseolus vul- garis pods and Helicoverpa armigera eggs[J].Plant Pro-tection，41（4）：29-34.]doi：10.3969/j.issn.0529-1542.2015. 04.005.

王四寶，曲爽.2017.昆蟲共生 菌及其在病蟲害防控中的應(yīng)用 前 景[J]. 中 國科學(xué) 院院刊，32（8）：863-872.[Wang SB， Qu S.2017.Insect symbionts and their potential applica- tion in pest and vector-borne disease control[J].Bulletin ofChinese Academy of Sciences，32（8）：863-872.]doi： 10.16418/j.issn.1000-3045.2017.08.008.

姚明 勇.2019.叉角厲蝽生物學(xué)及捕食作用研究[D]. 貴陽： 貴州大學(xué).[Yao MY.2019.Studies on biology and preda cious function of Eocanthecona fiurcellata（Wolff）[D]. Guiyang：Guizhou University.]doi：10.27047/d.cnki.ggu- du.2019.000066.

張曼，高平，趙航，周辰彥，梁晨，湯永玉，邢孔政，吳國星，高 熹.2022.捕食性天敵叉角厲蝽生長發(fā)育、繁殖及各蟲態(tài) 形態(tài)特征觀察[J]. 南方農(nóng)業(yè)學(xué)報，53（4）：1078-1087. [Zhang M，Gao P，Zhao H，Zhou CY，Liang C，TangYY， Xing K Z，Wu G X，Gao X.2022.Development，fecundity and morphological characteristics of the predatory Eocan-thecona fiucellata（Wolff）[J].Journal of Southern Agri-culture，53（4）：1078-1087.]doi：10.3969 .issn.2095-1191.2022.04.021.

張振宇，圣平，黃勝威，趙永順，張宏宇.2017.昆蟲腸道微生 物的多樣性、功能及應(yīng)用[J]. 生物資源，39（4）：231-239. [Zhang ZY，Sheng P，HuangS W，Zhao YS，Zhang HY. 2017.Diversity，function and application of insect gut microbiota[J].Biotic Resources，39（4）：231-239.]doi：10. 14188/j.ajsh.2017.04.001.

趙輝.2018.中華通草蛉內(nèi)生菌多樣性研究[D]. 石河子：石 河子大學(xué).[Zhao H.2018.The studyon endophytic bacte-ria diversity of the Chrysoperla sinica（Tjeder）[D].Shi-hezi：Shihezi University.]

趙航，吳國星，湯永玉，張曼，梁晨，邢孔政，蘭明先，高熹. 2022. 叉角厲蝽觸角轉(zhuǎn)錄組及嗅覺相關(guān)基因分析[J]. 環(huán) 境昆蟲學(xué)報，44（5）：1205-1217.[Zhao H，Wu GX，Tang Y Y，Zhang M，Liang C，Xing K Z，Lan M X，Gao X. 2022.Analysis of the antennal transcriptome and olfaction-related genes of Eocanthecona fucellata[J].Journal of"Environmental Entomology，44（5）：1205-1217.]doi：10. 3969/j.issn.1674-0858.2021.05.15.

鄭林宇，倫才智，柳麗君，李志紅.2022.昆蟲共生菌調(diào)控宿主 生長發(fā)育和生殖的研究進(jìn)展[J]. 植物保護(hù)學(xué)報，49（1）： 207-219.[Zheng L Y，Lun C Z，Liu L J，Li Z H.2022. Influencesof insect symbionts on host growth，develop- ment andreproduction：A review[J].Journal of Plant Pro- tection，49（1）：207-219.]doi：10.13802/j.cnki.zwbhxb.2022."2022824.

Ahsaei S M，Hosseininaveh V，Talaei-Hassanloui R，Zamani M.2020.Contribution of bacterial gut symbionts to diges- tion and development in Podisus maculivenris（Hemip- tera：Pentatomidae）[J].The Proceedings of the National Academy of Sciences，India Section B：Biological Scien- ces，90：959-967.doi：10.1007/s400 11-020-01165-3.

Bolyen E，Rideout JR，Dillon MR，Bokulich NA，Abnet CC， Al-Ghalith G A，Alexander H，Alm E J，Arumugam M， Asnicar F，Bai Y，Bisanz JE，Bittinger K，Brejnrod A，Bri- slawn CJ，BrownCT，Callahan B J，Caraballo-Rodriguez AM，Chase J，Cope EK，Silva R D，Diener C，Dorrestein PC，Douglas G M，Caporaso J G.2019.Reproducible， interactive，scalableand extensible microbiome data science using QIIME 2[J].Nature Bioechnology，37：852-857. doi：10.1038/s41587-019-0209-9.

Boone CK，Keefover-Ring K，Mapes A C，Adams A S，Bohl- mann J，Raffa K F.2013.Bacteria associated with a tree- killing insect reduce concentrations of plant defense com- pounds[J].Journal of Chemical Ecology，39：1003-1006.doi： 10.1007/s10886-013-0313-0.

Chaudhury M F，Skoda S R，Sagel A，Welch J B.2010.Vola- tiles emitted from eight wound-isolated bacteria differen- tially attract gravid screwworms（Diptera：Calliphoridae） to oviposit[J].Journal of Medical Entomology，47（3）：349- 354.doi：10.1603/me09235.

Colman DR，Toolson EC，Takacs-Vesbach C D.2012.Do diet and taxonomy influence insect gut bacterial communities？ [J].Molecular Ecology，21（20）：5124-5137.doi：10.1111/ j.1365-294X.2012.05752.x.

Edgar R C.2013.UPARSE：Highly accurate OTU sequences from microbial amplicon reads[J].Nature Methods，10：996- 998.doi：10.1038/nmeth.2604.

Engel P，Moran N A.2013.The gut microbiota of insects- diversity in structure and function[J].FEMS Microbio- logy Reviews，37（5）：699-735.doi：10.1111/1574-6976.12 025.

Fernandes A D，Reid JN，Macklaim JM，McMurrough T A， Edgell D R，Gloor G B.2014.Unifying the analysis of high-throughput sequencing datasets：Characterizing RNA- seq，16S rRNA gene sequencing and selective growth ex- periments by compositional data analysis[J].Microbiome， 2（1）：15.doi：10.1186/2049-2618-2-15.

Hu Y，Xie H G，Gao MZ，Huang P，Zhou H L，Ma Y B，Zhou MY，LiangJY，Yang J，Lü ZY.2020.Dynamic of compo- sition and diversity of gut microbiotain Triatoma rubrofas- ciata in different developmental stages and environmenta conditions[J].Frontiers inCellular and Infection Microbio- logy，10：587708.doi：10.3389/fcimb.2020.587708.

Keerthi MC，Sravika A，Mahesha H S，Gupta A，Bhargavi HA， Ahmed S.2020.Performance of the native predatory bug"Eocanthecona furcellata（Wolff）（Hemiptera：Pentatomi-dae ），on the fall armyworm，Spodoptera frugiperda（J.E. Smith）（Lepidoptera：Noctuidae），and its limitation under field condition[J].Egyptian Journal of Biological Pest Con- trol，30：69.doi：10.1 186/s41938-020-00272-7.

Kuang ZL，Wen J，Zhu YJ，He XF，Chen KW.2022.Dietary association with midgut microbiota components of Eocan-thecona furcellata（Wolff）[J].Diversity，14（12）：1130.doi"10.3390/d14121130.

Langille M G I， Zaneveld J，Caporaso J G，MeDonald D， Knights D，Reyes J A，Clemente J C，Burkepile D E， Thurber R LV，Knight R，Beiko R G，Huttenhower C. 2013.Predictive functional profiling of microbial commu nities using 16S rRNA marker gene sequences[J].Nature Biotechnology，31（9）：814-82 1.doi：10.1038/nbt.2676.

LiGN，Sun JJ，Meng YJ，Yang CF，Chen Z，WuYF，Tian L， Song F，CaiWZ，Zhang X，Li H.2022.The impact of environmental habitats and diets on the gut microbiota diversity of true bugs（Hemiptera：Heteroptera）[J].Biology， 11（7）：1039.doi：10.3390 /biology11071039

LiWJ，Wang L L，Jaworsk CC，Yang F，Liu B，Jiang YY，Lu YH，WuKM，Desneux N.2020.Theoutbreaksof nontar- get mirid bugs promote arthropod pest suppression in Bt cotton agroecosystems[J].Plant Biotechnology Journal，18（2）：322-324.doi：10.1111/pbi.13233.

Liu ZH，YangZW，Zhang J，LuoJY，Men Y，Wang YH，Xie """"" Q.2022.Stage correlationof symbiotic bacterial commu- "" nity and function in the development of litchi bugs " （Hemiptera：Tessaratomidae ）[J].Antonie van Leeuwen- """" hoek，115（1）：125-139.doi： 10.1007/s10482-021-01685-6. Luo J，Chen Y X，Guo L B，Wang A，Lu M，XuL T.2021. """" Variation of gut microbiota caused by an imbalance diet is "" detrimentalto bugs'survival[J].Science of the Total Envi- """""" ronment，771：144880.doi：10.1016/j.scitotenv.2020.144880. Nishiwaki H，Ito K，Shimomura M，Nakashima K，Matsuda K.

2007.Insecticidal bacteria isolated from predatory larvae of the antlion species Myrmeleon bore（Neuroptera：Myr- meleontidae）[J].Journal of Invertebrate Pathology，96（1）： 80-88.doi：10.1016/j.jip.2007.02.007.

Pan YN，Zhao R N，F(xiàn)u D，Yu C，Pan C N，ZhouW，Chen WL. 2022.Assessment of suitable reference genes for qRT-PCR"normalization in Eocanthecona fiurcellata（Wolff）[J].In-"sects， 13（9）：773.doi：10.3390/insects13090773.

Parks D H，Tyson G W，Hugenholtz P，Beiko R G.2014. STAMP：Statistical analysis of taxonomic and functional profiles [J].Bioinformatics，30（21）：3123-3124.doi：10. 1093/bioinformatics/btu494.

Sarkar S，Babu A，Chakraborty K，Deka B，Roy S.2021.Eocan- theconafurcellata（Wolf）（Hemiptera：Pentatomidae），a"potential biocontrol agent of theblack inchworm，Hypasidra talaca Walker（Lepidoptera：Geometridae）infesting tea [J].Phytoparasitica，49（2）：363-376.doi：10.1007/s12600-021-00888-x.

Shan HW，Wu W，Sun ZT，Chen J P，LiHJ.2021.The gut microbiota of the insect infraorder Pentatomomorpha（Hemiptera：Heteroptera）for the light of ecology and evo lution[J].Microorganisms，9（2）：464.doi： 10.3390/mic- roorganisms9020464.

Song XM，Zhong ZW，Gao L，Weiss B L，Wang JW.2022. Metabolic interactions between disease-transmitting vec- tors and theirmicrobiota[J].Trends in Parasitology，38（8）： 697-708.doi：10.1016/ j.pt.2022.05.002

Tilocca B，Burbach K，HeyerC M E，Hoelzle L E，Mosenthin R，Stefanski V，Camarinha-Silva A，Seifert J.2017.

Dietary changes in nutritional studies shape the structural and functional composition of the pigs'fecal microbiome— From days to weeks[J].Microbiome，5：144.doi：10.1186/ s40168-017-0362-7.

Wang Q，Garrity G M，Tiedje J M，Cole J R.2007.Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy[J].Applied and Environ- mental Microbiology，73（16）：5264-5267.doi：10.1128/ AEM.00062-07.

White J R，Nagarajan N，Pop M.2009.Statistical methods for detecting differentially abundant features in clinical meta- genomic samples[J].PLoS Computational Biology，5（4）： e1000352.doi： 10.1371/journal.pcbi.1000352.

Xue H，Zhu XZ，Wang L，Zhang K X，LiDY，Ji JC，NiuL， Gao XK，Luo JY，Cui J J.2023.Dynamics and diversity of symbiotic bacteria in Apolygus lucorum at different developmental stages[J].Journal of Cotton Research，6：5. doi：10.1186/s42397-023-00142-1.

Xue H，Zhu XZ，Wang L，Zhang KX，Li D Y，JiJC，Niu L， Wu CC，Gao XK，Luo JY，Cui JJ.2021.Gut bacterial diversity in different life cycle stages of Adelphocoris suturalis（Hemiptera：Miridae）[J].Frontiers in Microbio- logy，12：670383.doi：10.3389/fmicb.202 1.670383.

Yun JH，Roh SW，Whon TW，Jung M J，Kim M S，Park D S， Yoon C，Nam YD，KimYJ，ChoiJH，KimJY，Shin NR， Kim S H，Lee W Bae JW.2014.Insects gut bacterial diversity deermined by host environmental habitat，diet， developmental stage and phylogeny of host [J].Applied and Environmental Microbiology，80（17）：5254-5264.doi："10.1128/AEM.01226-14.

Zhao CC，Zhao H，Zhang S，LuoJY，Zhu XZ，Wang L，Zhao P，Hua HX，CuiJJ.2019.The developmental stage sym- bionts of the pea aphid-feeding Chrysoperla sinica（Tjeder） [J].Frontiers in Microbiology，10：2454.doi：10.3389/ fmicb.2019.02454.

（責(zé)任編輯 麻小燕）