點蜂緣蝽滯育和非滯育對嗅覺相關基因表達的影響

摘要

有研究表明滯育和非滯育點蜂緣蝽Riptortus pedestris對聚集信息素的感受存在顯著差異，本研究擬通過測定滯育和非滯育點蜂緣蝽觸角轉錄組，尋找表達差異的嗅覺相關基因，分析與滯育關聯的信息素識別機制。利用高通量測序平臺BGISEQ500對滯育和非滯育點蜂緣蝽雌、雄成蟲觸角進行轉錄組測序，篩選差異表達基因并進行通路富集分析，同時對差異表達的嗅覺基因進行RTqPCR驗證。點蜂緣蝽觸角轉錄組共獲得57 648條unigenes，其中61.05%在七大公共數據庫中成功注釋。注釋到NR數據庫的基因最多（52.90%），匹配到茶翅蝽Halyomorpha halys的序列最多，為55.83%。通過轉錄組數據分析，與非滯育雌蟲相比，滯育雌蟲觸角中共篩選到2 222個差異表達基因，其中上調表達基因1 288個，下調表達基因934個;與非滯育雄蟲相比，滯育雄蟲觸角中共篩選到2 293個差異表達基因，其中上調表達基因1 410個，下調表達基因883個;雌、雄蟲觸角間差異基因較少，分別為315個（滯育）和584個（非滯育）。進一步對嗅覺相關基因進行分析發現：滯育和非滯育雌蟲觸角間10個化學感受蛋白（chemosensory proteins，CSPs）基因和11個氣味結合蛋白（odorantbinding proteins，OBPs）基因差異表達;滯育和非滯育雄蟲觸角有7個CSPs基因和11個OBPs基因差異表達;滯育雌、雄成蟲僅有1個OBP基因表達出現差異;非滯育雌、雄蟲間有4個CSPs基因和3個OBPs基因差異表達。從差異表達的嗅覺相關基因中選取10個基因（4個OBPs，4個CSPs 2個ORs）進行RTqPCR驗證，結果與轉錄組測序相符。研究結果將為今后點蜂緣蝽識別聚集信息素的分子機制及研發更為高效的化學通訊調控技術提供參考。

關鍵詞

點蜂緣蝽; 觸角轉錄組; 生殖滯育; 嗅覺相關基因; 差異表達基因

中圖分類號：

S 476.2

文獻標識碼： A

DOI： 10.16688/j.zwbh.2024416

Effects of diapause and nondiapause on the expression of olfactoryrelated genes in Riptortus pedestris

XIE Minghui， ZHONG Yongzhi， LIN Lulu， ZHANG Guangling， ZHAO Wei， CHEN Haoliang*

（Institute of Plant Protection and Agroproducts Safety， Anhui Academy of Agricultural Sciences， Hefei 230001， China）

Abstract

Previous studies have shown significant differences in the response to aggregation pheromones between diapause and nondiapause individuals of Riptortus pedestris. To explore the molecular mechanisms underlying pheromone recognition associated with diapause， the antennal transcriptomes of diapause and nondiapause male and female adults were sequenced using the BGISEQ500 highthroughput sequencing platform， screened for differentially expressed genes （DEGs）， performed pathway enrichment analysis， and validated differentially expressed olfactoryrelated genes by RTqPCR. A total of 57 648 unigenes were obtained from the antennal transcriptome， with 61.05% successfully annotated in seven public databases. The highest annotation rate was to the NR database （52.90%）， with most matches corresponding to Halyomorpha halys （55.83%）. Transcriptome analysis identified 2 222 DEGs between diapause and nondiapause females （1 288 upregulated， 934 downregulated） and 2 293 DEGs between diapause and nondiapause males （1 410 upregulated， 883 downregulated）. Comparatively fewer DEGs were found between sexes within the same diapause status： 315 DEGs for diapause and 584 DEGs for nondiapause groups. Analysis of olfactoryrelated genes revealed that ten chemosensory protein （CSP） genes and eleven odorantbinding protein （OBP） genes were differentially expressed between diapause and nondiapause females， while seven CSP genes and eleven OBP genes were differentially expressed between diapause and nondiapause males. Only one OBP gene was differentially expressed between diapause females and males， whereas four CSP genes and three OBP genes differed between nondiapause females and males. Ten olfactoryrelated genes （four OBP genes， four CSP genes） and two OR genes were selected for RTqPCR validation， and the results were consistent with the transcriptome data. These results provide a reference for understanding the molecular mechanisms of pheromone recognition in R.pedestris and support the development of more efficient chemical communication regulation technologies.

Key words

Riptortus pedestris; antennal transcriptome; reproductive diapause; olfactoryrelated genes; differentially expressed genes

點蜂緣蝽Riptortus pedestris屬半翅目Hemiptera蛛緣蝽科Alydidae，是重要的刺吸性害蟲［1］。該蟲每年可發生2～3代，成蟲和若蟲均以刺吸植株汁液為食，主要為害大豆、蠶豆、豌豆等豆科植物，可造成植株貪青、空莢癟粒等“癥青”現象［23］。成蟲具有很強的飛行能力，同時可傳播多種病原菌，嚴重影響豆科作物的產量和品質［2，4］。近年來，點蜂緣蝽對我國大豆的為害加重，在河南、河北、山東、安徽、北京、山西、陜西等地均有暴發，且呈逐年加重趨勢，已成為我國大豆生產上亟待解決的重要問題［1， 57］。利用聚集信息素誘殺點蜂緣蝽是有效的綠色防控手段，具有廣闊的應用前景［8］。

昆蟲通過嗅覺感知環境中的氣味，完成寄主定位、求偶、交配、尋找產卵場所和躲避敵害等行為活動［910］。觸角是昆蟲重要的嗅覺感受器官，通過一系列復雜的生物化學反應來實現對各種氣味分子的精確靈敏識別和特異性信號轉導［1011］。觸角感器中的多種蛋白參與構建昆蟲的外周嗅覺系統，主要有化學感受蛋白（chemosensory proteins，CSPs）、氣味結合蛋白（odorantbinding proteins，OBPs）、 味覺受體（gustatory receptors，GRs）、嗅覺受體（olfactory receptors，ORs）、離子型受體（ionotropic receptors，IRs）和感覺神經元膜蛋白（sensory neuron membrane proteins，SNMPs）等［12］。嗅覺相關基因與昆蟲的滯育息息相關，切除果蠅Drosophila melanogaster觸角會縮短滯育成蟲的壽命，降低滯育解除后的繁殖力，嗅覺受體神經元和溫度感應神經元是滯育成功恢復所必需的［13］。尖音庫蚊Culex pipiens介導滯育的核心鐘基因Cycle能夠調控嗅覺相關基因的表達［14］。

光周期是誘導點蜂緣蝽生殖滯育的主要因子，成蟲在短日照（lt;13.5 h）進入生殖滯育狀態［15］。此時，雌蟲卵巢停止發育;雄蟲的附腺囊空癟，不再合成和釋放聚集信息素［16］。點蜂緣蝽的聚集信息素是由雄性成蟲分泌釋放的，主要包括3個組分：異丁酸十四酯（tetradecyl isobutyrate， 14∶iBu）、反2己烯基順3己烯酸酯［（E）2hexenyl （Z）3hexenoate， E2HZ3H］和反2己烯基反2己烯酸酯［（E）2hexenyl （E）2hexenoate， E2HE2H］［17］。點蜂緣蝽雄蟲在非滯育狀態、寄主植物適宜且充足、營養狀態良好的前提下合成和分泌聚集信息素;而其他個體在食物源缺乏時，對聚集信息素的敏感性會顯著增加［1619］。在田間，秋季來臨后日照時間逐漸縮短，食物的數量和質量都隨之降低。此時，進入生殖滯育的點蜂緣蝽雌蟲雖不與雄蟲交配，但對作為食物源信號的聚集信息素更加敏感［15］。由此推斷，參與聚集信息素識別的嗅覺相關基因在滯育和非滯育點蜂緣蝽雌蟲間的表達可能存在差異。

為了探索點蜂緣蝽嗅覺相關基因在聚集信息素識別中的作用，本研究借助高通量測序對滯育和非滯育點蜂緣蝽成蟲觸角轉錄組進行生物學分析，發掘參與嗅覺感受的重要基因，同時篩選出差異表達的嗅覺相關基因。研究結果將為今后點蜂緣蝽識別聚集信息素的分子機制及研發更為高效的化學通訊調控技術提供參考。

1 材料與方法

1.1 供試蟲源

點蜂緣蝽成蟲采集于宿州市農業科學院基地（117.09°E， 33.65°N），在人工氣候室內的養蟲籠（50 cm×50 cm×50 cm，鋁合金框架，尼龍網罩，紗網網孔為50目）內用玉米和四季豆飼養擴繁建立實驗種群。飼養條件為溫度（24±1）℃，相對濕度為（70±5）%，光周期為L∥D=16 h∥8 h。待點蜂緣蝽在室內飼養5代以上后，采集蟲卵并置于新的養蟲籠中，將光周期分別調整為L∥D=16 h∥8 h和L∥D=10 h∥14 h，其他飼養條件不變。光周期L∥D=16 h∥8 h條件下點蜂緣蝽可以正常繁殖;而光周期L∥D=10 h∥14 h條件下點蜂緣蝽進入生殖滯育狀態［1415］。將2個光周期下羽化5 d的雌雄成蟲觸角剪下，經液氮研磨后放入-80℃超低溫冰箱保存待用，每處理15頭，重復3次。

1.2 觸角總RNA的提取

采用TRIzolReagent（Invitrogen， Carlsbad， CA， USA）分別提取滯育和非滯育點蜂緣蝽雌、雄成蟲觸角總RNA，并用DNase I（TaKaRa， Dalian， China）清除RNA樣品中殘留的基因組DNA。用Nano Drop 2000（Thermo Scientific， Wilmington， DE， USA）檢測RNA的純度，使用Fragment Analyzer（Agilent 5300， CA， USA）生物分析儀檢測樣品RNA的濃度和完整性，用于后續轉錄組測序。

1.3 cDNA文庫的構建和測序

對檢測合格的RNA樣品進行cDNA文庫的構建和轉錄組測序，基本流程為：取一定量的RNA樣品，用帶有Oligo（dT）的磁珠富集有polyA尾巴的mRNA;用打斷緩沖液將mRNA片段化;轉錄后合成雙鏈cDNA;修復雙鏈cDNA末端，并在3′末端加上A堿基，配制接頭連接反應體系，使接頭與cDNA連接，對產物進行擴增和質檢。文庫質檢后，將PCR產物熱變性成單鏈，用一段橋式引物將單鏈DNA環化得到單鏈環狀DNA文庫，消化掉未被環化的線性DNA分子;單鏈環狀DNA分子通過滾環復制，形成一個包含多個拷貝的DNA納米球（DNB）。將得到的DNBs采用高密度DNA 納米芯片技術，加到芯片上的網狀小孔內，通過聯合探針錨定聚合技術（cPAS）進行測序。

1.4 轉錄組數據組裝和功能注釋

測序得到的原始數據（raw data）使用過濾軟件SOAPnuke進行數據過濾，去除包含接頭的reads、未知堿基（N堿基）含量大于5%的reads 和低質量的reads，得到高質量數據（clean data）。使用Trinity對clean reads進行de novo組裝，同時通過BUSCO評估組裝質量。通過軟件Transdecoder識別unigene中候選編碼區域，BLAST比對SwissProt并且用Hmmscan搜索Pfam蛋白同源序列，預測編碼區域。對組裝得到的unigene進行功能數據庫注釋。

1.5 差異表達基因分析

使用DESeq2進行組內差異基因分析，條件為差異倍數（foldchange）≥2和qvalue≤0.05。使用PossionDis進行組間差異分析，條件為foldchange≥2和錯誤發現率（1 discovery rate， FDR）≤0.05。根據GO及KEGG注釋結果以及官方分類，將差異基因進行功能分類，使用R軟件中的phyper進行KEGG富集分析，并使用TermFinder包進行GO富集分析。以qvalue≤0.05為閾值，滿足此條件的定義為在候選基因中顯著富集。

1.6 嗅覺相關基因的分析和RTqPCR驗證

從觸角轉錄組中篩選和鑒定OBPs、CSPs、Ors、GRs、IRs 和 SNMPs 等嗅覺基因，將得到的基因序列在NCBI數據庫中進行Blastx比較。選擇10個嗅覺相關基因（4個OBPs，4個CSPs和2個ORs）進行RTqPCR驗證。使用Primer Premier 5設計特異性引物，選擇βtub和EF1a作為內參基因（表1）。每個樣本取1 μg總RNA作為模板反轉錄cDNA，按照

SYBR Primer Script RTPCR Kit （TaKaRa）構建反應體系，并在CFX96系統（BioRad， Hercules， CA， USA）上進行檢測。反應體系：SYBR Premix 12.5 μL， cDNA 2 μL， 10 μmol/L上、下游引物各1 μL， ddH2O 8.5 μL。

反應程序： 95℃ 30 s; 95℃ 30 s， 58℃ 30 s， 72℃ 30 s， 45個循環。每個處理均設3個生物樣品和2個技術重復。

1.7 數據處理與分析

RTqPCR結果通過內參基因標準化，并使用2-ΔΔCT計算所測基因的相對表達量;轉錄組測序采用FPKM（fragments per kilobase million， FPKM）值對基因表達量進行評估，根據不同樣本間的比較計算出相對表達量。使用Graphpad Prism 5進行圖形繪制。

2 結果與分析

2.1 測序數據質量評估

通過對滯育和非滯育點蜂緣蝽雌、雄蟲觸角樣本進行轉錄組測序，通過數據質控和過濾后，每個樣本均獲得6G以上的數據，Clean reads占raw reads百分比在93.97%～97.77%，Q20堿基百分比在97.52%以上，Q30堿基百分比在93.31%以上，GC含量在33.80%～34.25%（表2）。

使用Trinity軟件對轉錄組測序結果進行組裝，共獲得57 648條unigenes，平均長度為1 555 bp，N50為2 848 bp。在所有觸角樣本中均有表達的unigene數量為49 211個，在非滯育雌蟲、非滯育雄蟲、滯育雌蟲和滯育雄蟲中特異表達的unigene數量分別為371、373、401個和319個（圖1）。

2.2 Unigenes的功能注釋

將獲得的unigene在7個數據庫中進行功能注釋，61.05%（34 194個）的unigenes，注釋成功（表3）。

其中，注釋基因數量最高的數據為NR和KEGG，分別為30 494個和24 209個，占比分別為52.90%和41.99%。NR庫注釋比對發現，注釋到茶翅蝽Halyomorpha halys的序列最多，為55.83%;注釋到先前報道的點蜂緣蝽序列的比例為9.18%;其余依次為溫帶臭蟲Cimex lectularius（4.74%）、錘脅蹺蝽Yemma signatus（1.43%）和褐飛虱Nilaparvata lugens（0.96%）（圖2）。

1） NR：非冗余蛋白序列數據庫;NT：核酸序列數據庫;GO：基因本體數據庫;KEGG：京都基因與基因組百科全書;SwissProt：SwissProt蛋白質序列數據庫;Pfam：蛋白質家族數據庫;COG：直系同源蛋白質簇數據庫。

NR： Nonredundant protein sequence database; NT： Nucleotide sequence database; GO： Gene ontology database; KEGG： Kyoto encyclopedia of genes and genomes; SwissProt： SwissProt protein sequence database; Pfam： The protein family database; COG： Clusters of orthologous groups database.

2.3 差異表達基因的GO功能富集分析

與非滯育雌蟲相比，滯育雌蟲觸角中共篩選到2 222個差異表達基因，包括1 288個上調表達基因和934個下調表達基因。與非滯育雄蟲相比，滯育雄蟲觸角中篩選到的差異表達基因2 293個，其中1 410個上調表達，883個下調表達。滯育雌、雄蟲觸角中差異表達的基因較少，僅有315個;非滯育雌、雄蟲觸角中差異表達的基因為584個。

功能注釋，共富集到生物學過程（biologicalprocesses）、細胞組分（cellular component）和分子功能（molecular function）三大類別的46個功能分類中（圖3）。

滯育和非滯育雌蟲差異表達基因注釋到生物學過程、細胞組分和分子功能的百分比分別為28.40%、42.12%和29.47%。在生物學組分中，滯育和非滯育雌蟲差異表達基因主要富集到細胞過程、代謝過程、生物調節等;在細胞組分中，注釋到20種與細胞和機體構建相關，如細胞膜、細胞膜部位和細胞有關的功能分類;而在分子功能中，注釋到9種功能，主要為結

合和催化活性。滯育和非滯育雄蟲差異表達基因注釋到生物學過程（28.54%）、細胞組分（44.43%）、分子功能（27.03%）共45個功能分類，其中差異表達基因顯著富集的功能注釋與滯育和非滯育雌蟲相同。滯育雌雄蟲之間和非滯育雌雄蟲之間差異表達的基因數量較少，其中滯育雌雄蟲之間僅注釋到31個功能分類，非滯育雌雄蟲之間注釋到36個功能分類。

2.4 差異表達基因的KEGG富集分析

將不同處理觸角比對產生的差異表達基因進行KEGG富集分析（圖4），結果發現：滯育和非滯育雌蟲共有947個差異表達基因比對到KEGG中的

235個代謝通路里，富集的通路主要為蛋白代謝（ko04974）、溶酶體（ko04142）和藥物代謝（ko00983）;

滯育和非滯育雄蟲共有1 014個差異表達基因比對到KEGG中的245個代謝通路里，富集的通路主要為藥物代謝、溶酶體和吞噬體（ko04145）;滯育雌雄蟲之間僅有124個差異表達基因比對到KEGG代謝通路，主要涉及昆蟲激素合成（ko00981）、mRNA調控途徑（ko03015）和mTOR信號通路（ko04150）;非滯育雌雄蟲之間共有229個差異表達基因顯著富集到昆蟲激素合成、藥物代謝和溶酶體。

2.5 嗅覺相關基因篩選分析

根據關鍵詞搜索和注釋基因分析，從點蜂緣蝽

觸角轉錄組共鑒定出12個CSPs、37個OBPs、8個GRs、60個ORs、14個IRs和3個SNMPs基因，進

一步對點蜂緣蝽滯育和非滯育觸角轉錄組間的差異表達基因進行分析，我們發現：滯育和非滯育雌蟲觸角轉錄組中有10個差異表達的CSPs基因，其中1個上調表達，9個下調表達;11個OBPs基因，其中3個上調表達，8個下調表達;3個ORs基因上調表達和1個OR基因下調表達;1個IR基因上調表達。滯育與非滯育雄蟲比較，共得到1個上調表達和6個下調表達的表達的OBPs基因;差異表達的ORs基因有16個，7個上調表達，9個下調表達。同是滯育狀態下的雌、雄成蟲僅有1個OBP基因表達出現差異;4個ORs基因上調表達，2個ORs基因下調表達。非滯育狀態下的雌、雄蟲比較發現，4個CSPs基因差異表達，1個基因上調表達和3個基因下調表達;3個OBPs基因上調表達;18個ORs基因差異表達，其中9個上調表達，9個下調表達（表4）。

2.6 差異表達的嗅覺相關基因的RTqPCR驗證

由表4我們可以看出，滯育和非滯育成蟲觸角中差異表達的嗅覺基因較多，特別是OBPs和CSPs，因此，我們選取10個差異表達的嗅覺基因（包括4個OBPs，4個CSPs和2個ORs）進行RTqPCR驗證。結果表明：2個OBPs （CL1466.Contig3和CL418.Contig5）和1個OR（CL2647.Contig1）在滯育雌雄成蟲觸角的表達量均高于非滯育雌雄成蟲;CSP基因（Unigene9693）在滯育雌蟲觸角中上調表達;其他所選OBPs，CSPs和ORs基因在滯育雌雄蟲觸角中均下調表達;RTqPCR結果與轉錄組測序結果基本一致（圖5）。

3 結論與討論

本研究利用BGISEQ500測序平臺對滯育和非滯育點蜂緣蝽雌、雄成蟲觸角的轉錄組進行了測序，共獲得57 648條unigenes，平均長度為1 555 bp，N50為2 848 bp，高于宋月芹等［20］的測序結果（平均長度618 bp，N50為1 013 bp），經過功能注釋和同源搜索，共有34 194條unigenes注釋成功，占所有unigenes的61.05%。而宋月芹等［20］的測序結果得到21 365條注釋成功的unigenes，僅占總unigenes數量的23.15%。這證明我們得到高質量的點蜂緣蝽觸角轉錄組測序及組裝結果。同時，隨著高通量測序技術的發展，測序的深度和靈敏度都有很大提升，加上各大數據庫儲存基因的數量和種類日趨豐富，成功注釋的unigene數量較之前顯著增加。采用NR蛋白數據庫比對，共有30 494條unigenes成功注釋，占比為52.90%。其中，注釋到茶翅蝽的基因比例最多，為55.83%，其次為先前報道的點蜂緣蝽序列（9.18%）。而宋月芹等［20］的比對結果顯示，注釋基因比例最高的是內華達古白蟻Zootermopsis nevadensis（15.22%），其次為赤擬谷盜Tribolium castaneum（5.2%），這主要歸因于內華達古白蟻和赤擬谷盜具有相對完整的全基因組序列和豐富的基因組注釋信息［2122］，而半翅目昆蟲的基因數據還不夠豐富［23］。隨著半翅目昆蟲研究的不斷深入，本次測序點蜂緣蝽觸角轉錄組中注釋到的基因與親緣關系更近的同為半翅目的茶翅蝽相似序列最多。

我們比較了滯育雌、雄蟲，非滯育雌、雌蟲4個處理觸角的轉錄組信息。其中，點蜂緣蝽雌蟲和雄蟲觸角間差異表達基因數量較少，滯育和非滯育狀態下分別為315個和584個;而滯育和非滯育相比，差異表達的基因數量較多，其中雌蟲為2 222個，雄蟲為2 293個。從差異表達的基因數量上來看，滯育和非滯育因素帶來的差異比性別因素的差異要大。昆蟲滯育或休眠的轉錄組分析多見于整個成蟲，涉及激素調控、脂肪酸代謝等［2426］，關注滯育前后觸角中差異基因的研究較少。GO功能注釋結果顯示，差異表達的基因多集中在結合、催化以及細胞過程。KEGG代謝途徑富集分析顯示，滯育和非滯育成蟲觸角差異表達基因富集在蛋白代謝（ko04974）、溶酶體（ko04142）、藥物代謝（ko00983）和吞噬體（ko04145），多數基因在滯育成蟲中是上調表達的，這可能與昆蟲在滯育期間需要積累能量，上調免疫能力來應對滯育時的不良環境，并為滯育解除后的生長發育做準備［27］。

觸角是昆蟲嗅覺系統的主要器官，昆蟲的嗅覺感受和信號轉導依靠各種嗅覺蛋白的共同參與，構成昆蟲高度靈敏發達的嗅覺感受系統［28］。近年來，隨著高通量測序和全基因組測序的高速發展，眾多非模式昆蟲的嗅覺相關基因得到了深度挖掘和研究。本研究通過BLAST和同源搜索，共鑒定出134個嗅覺相關基因，包括12個CSPs基因、37個OBPs基因、8個GRs基因、60個ORs基因、14個IRs基因和3個SNMPs基因。與宋月芹等［20］從點蜂緣蝽轉錄組中鑒定出嗅覺相關基因相比，OBP基因數量較多，OR基因數量較少，其他基因在合理范圍內。在點蜂緣蝽全基因組序列公布后，學者們從基因組中鑒定出更多的嗅覺相關基因［2931］。本研究的轉錄組測序是在點蜂緣蝽全基因組序列公布之前，且只對昆蟲觸角做了轉錄組分析，因此鑒定出的嗅覺相關基因是在觸角中特異性表達的基因，數量較少，后續還要結合全基因組對鑒定出的嗅覺相關基因做進一步分析。

利用FPKM值分析點蜂緣蝽觸角中嗅覺相關基因的表達量，以及這些基因在滯育、非滯育雌、雄觸角中的差異表達情況。我們選取了10個差異表達的嗅覺基因并利用熒光定量 RTqPCR 進行驗證。 非滯育點蜂緣蝽雌、雄成蟲觸角中差異表達的嗅覺基因包括4個CSPs基因和3個OBPs基因，推測這些基因可能在性信息素識別過程中發揮著重要作用［32］。而滯育點蜂緣蝽成蟲兩性間僅有1個OBP基因表達出現差異，這可能是由于雌、雄成蟲進入滯育后，性腺停止發育，不再合成和釋放性信息素，雄蟲觸角中對性信息素的識別和信號轉導等活動也隨之停止，因而嗅覺相關基因在滯育成蟲兩性間的表達差異性降低［1516］。滯育點蜂緣蝽對聚集信息素的敏感度增加［16， 1819］，推測滯育和非滯育點蜂緣蝽觸角中差異表達的嗅覺相關基因與聚集信息素或食物源氣味物質的識別相關，有待后續深入研究。

本研究通過分析點蜂緣蝽滯育和非滯育雌雄蟲觸角的轉錄組差異，較詳細地解析了轉錄組間基因的表達水平、功能分類及代謝通路間的差異。進一步對這些差異表達的基因開展深入研究，將為今后利用嗅覺靶標防治點蜂緣蝽提供參考，并為研發更為高效的化學通訊調控技術奠定基礎。

參考文獻

［1］ 謝明惠， 鐘永志， 嚴丹侃， 等.安徽阜陽點蜂緣蝽發生動態及不同誘捕器的誘集效果比較［J］.中國植保導刊， 2022， 42（11）：3841.

［2］ 高宇， 史樹森.大豆“莢而不實”型“癥青”與蝽類害蟲為害的相關性及防控策略［J］.大豆科學， 2019， 38（4）：650655.

［3］ WEI Zhongyan， GUO Wenbin， JIANG Shanshan， et al. Transcriptional profiling reveals a critical role for GmFT2a in soybean staygreen syndrome caused by the pest Riptortus pedestris ［J］. New Phytologist， 2023， 237（5）： 18761890.

［4］ TSUNODA T， MORIYA S. Measurement of flight speed and estimation of flight distance of the bean bug， Riptortus pedestris （Fabricius） （Heteroptera： Alydidae） and the rice bug， Leptocorisa chinensis Dallas （Heteroptera： Alydida） with a speed sensor and flight mills ［J］. Applied Entomology and Zoology， 2008， 43（3）： 451456.

［5］ LI Jinbu， YIN Maozhu， YAO W C， et al. Genomewide analysis of odorantbinding proteins and chemosensory proteins in the bean bug Riptortus pedestris ［J/OL］. Frontiers in Physiology， 2022， 13： 949607. DOI： 10.3389/fphys.2022.949607.

［6］ 錢鑫鑫， 潘志金， 徐昊， 等.阜陽地區大豆病蟲害發生情況及防控建議［J］.大豆科技， 2022（6）：2429.

［7］ 孫曉輝， 吳斌， 王守海， 等.點蜂緣蝽在山東小麥、玉米、大豆上的季節性發生規律研究［J］.山東農業科學， 2023， 55（3）：153157.

［8］ ENDO N， WADA T， SASAKI R. Seasonal synchrony between pheromone trap catches of the bean bug， Riptortus pedestris （Heteroptera：Alydidae） and the timing of invasion of soybean fields ［J］. Applied Entomology and Zoology， 2011， 46（4）： 477482.

［9］ 王桂榮， 吳孔明， 郭予元.昆蟲感受氣味物質的分子機制研究進展［J］.農業生物技術學報， 2004， 12（6）：720726.

［10］BRITO N F， MOREIRA M F， MELO A C. A look inside odorantbinding proteins in insect chemoreception ［J］. Journal of Insect Physiology， 2016， 95： 5165.

［11］LEAL W S. Pheromone reception ［J］. Topics in Current Chemistry， 2005， 240： 136.

［12］鞏雪芳， 謝壽安， 楊平， 等.花椒窄吉丁觸角轉錄組及嗅覺相關基因分析［J］.昆蟲學報， 2020， 63（10）：11591170.

［13］EASWARAN S， MONTELL J D. A genomewide association study implicates the olfactory system in Drosophila melanogaster diapauseassociated lifespan extension and fecundity ［J/OL］. bioRxiv， 2025， DOI： 10.1101/2024.03.10.584341.

［14］PRABIN D， WEI Xueyan， MEGAN M， et al. Identification of CYCLE targets that contribute diverse features of circadian rhythms in the mosquito Culex pipiens ［J/OL］. Comparative Biochemistry and Physiology Part D： Genomics and Proteomics， 2023， 48： 101140. DOI： 10.1016/j.cbd.2023.101140.

［15］RAHMAN M M， LIM U T. Females of Riptortus pedestris （Hemiptera： Alydidae） in reproductive diapause are more responsive to synthetic aggregation pheromone ［J］. Journal of Economic Entomology， 2016， 109（5）： 20822089.

［16］MIZUTANI N， YASUDA T， YAMAGUCHI T， et al. Pheromone contents and physiological conditions of adult bean bugs， Riptortus pedestris （Heteroptera： Alydidae）， attracted to conspecific males during nondiapause and diapause periods in fields ［J］. Applied Entomology and Zoology， 2008， 43（3）： 331339.

［17］LEAL W S， HIGUCHI H， MIZUTANI N， et al. Multifunctional communication in Riptortus clavatus （Heteroptera： Alydidae）： Conspecific nymphs and egg parasitoid Ooencyrtus nezarae use the same adult attractant pheromone as chemical cue ［J］. Journal of Chemical Ecology， 1995， 21（7）： 973985.

［18］LIM E， KIM J， CHO K S， et al. Influences of food sources on reproduction and aggregation pheromone production of Riptortus pedestris （Hemiptera： Alydidae） in its adult stage ［J］. Entomological Research， 2015， 45（5）： 241246.

［19］SHIN H S， RYU T H， KWON H R. Aggregation pheromone secretion change of Riptortus pedestris （Hemiptera： Alydidae） depending on diet ［J］. Journal of the Faculty of Agriculture Kyushu University， 2015， 60： 357362.

［20］宋月芹， 董鈞鋒， 陳慶霄， 等.點蜂緣蝽觸角轉錄組及化學感受相關基因的分析［J］.昆蟲學報， 2017， 60（10）：11201128.

［21］RICHARDS S， GIBBS R A， WEINSTOCK G M， et al. The genome of the model beetle and pest Tribolium castaneum ［J］. Nature， 2008， 452 （7190）： 949955.

［22］TERRAPON N， LI Cai， ROBERTSON H， et al. Molecular traces of alternative social organization in a termite genome ［J/OL］. Nature Communications， 2014， 5： 3636. DOI： 10.1038/ncomms4636.

［23］AN Xingkui， SUN Liang， LIU Hangwei， et al. Identification and expression analysis of an olfactory receptor gene family in green plant bug Apolygus lucorum （MeyerDür） ［J/OL］. Scientific Reports， 2016， 6： 37870. DOI： 10.1038/srep37870.

［24］CUI Jin， LIN Kejian， XU Linbo， et al. Transcriptome analysis of beet webworm shows that histone deacetylase may affect diapause by regulating juvenile hormone ［J/OL］. Insects， 2022， 13（9）： 835. DOI： 10.3390/insects13090835.

［25］ZHOU Wenzhen， WU Youfang， YIN Zhiyong， et al. Juvenile hormone is an important factor in regulating Aspongopus chinensis Dallas diapause ［J/OL］. Frontiers in Physiology， 2022， 13： 873580. DOI： 10.3389/fphys.2022.873580.

［26］姚知含， 李玲， 龐保平， 等. 沙蔥螢葉甲對蛻皮激素響應的轉錄組分析［J］.環境昆蟲學報， 2023， 45（1）：155162.

［27］韓海斌， 姜珊， 林克劍， 等. 牧草優勢寄生性天敵傘裙追寄蠅滯育相關蛋白分析［J］. 植物保護學報， 2021， 48（6）：13591369.

［28］ZHAO Yunhe， DING Jinfeng， ZHANG Zhengqun， et al. Sexand tissuespecific expression profiles of odorant binding protein and chemosensory protein genes in Bradysia odoriphaga （Diptera： Sciaridae） ［J/OL］. Frontiers in Physiology， 2018， 9： 107. DOI： 10.3389/fphys.2018.00107.

［29］HUANG Haijian， YE Yuxuan， YE Zhuangxin， et al. Chromosomelevel genome assembly of the bean bug Riptortus pedestris ［J］. Molecular Ecology Resources， 2021， 21（7）： 24232436.

［30］LI Kai， ZHANG Xinxin， GUO Jianqiu， et al. Feeding of Riptortus pedestris on soybean plants， the primary cause of soybean stay green syndrome in the HuangHuaiHai river basin ［J］. The Crop Journal， 2019， 7（3）： 360367.

［31］LIU Panjing， GUO Jianglong， WEI Hongyi， et al. Genomewide identification of candidate chemosensory receptors in the bean bug Riptortus pedestris （Hemiptera： Alydidae） and the functional verification of its odorant receptor coreceptor （Orco） in recognizing aggregation pheromone ［J/OL］. Frontiers in Physiology， 2023， 14： 1224009. DOI： 10.3389/fphys.2023.1224009.

［32］楊爽， 趙慧婷， 徐凱， 等.大蠟螟觸角轉錄組及嗅覺相關基因分析［J］.應用昆蟲學報， 2019， 56（6）：12791291.

（責任編輯：楊明麗）