虱螨脲對棉鈴蟲葉酸一碳庫代謝途徑的影響

楊林 劉彩月 韓小強(qiáng) 陳明慧 陳閱新 張艷聰 周科瑾

摘要 為了探究虱螨脲對棉鈴蟲的作用機(jī)理，用虱螨脲（200 mg/L）處理棉鈴蟲3齡幼蟲，通過轉(zhuǎn)錄組測序、葉酸含量測定、葉酸飼喂試驗(yàn)和分子對接技術(shù)分析了虱螨脲對棉鈴蟲葉酸一碳庫代謝途徑的影響。結(jié)果表明，棉鈴蟲幼蟲差異基因數(shù)量隨虱螨脲處理后時(shí)間的延長而增多，主要富集在葉酸一碳庫、氧化還原活性、幾丁質(zhì)結(jié)合、氨基酸合成、嘌呤合成與代謝以及離子跨膜運(yùn)輸?shù)韧緩缴希渲腥~酸一碳庫代謝通路在處理后0.5 h極顯著富集（P=6.35×10-5）。棉鈴蟲幼蟲的葉酸含量隨著虱螨脲處理后時(shí)間的延長而顯著降低，處理后24 h葉酸含量為4.53 μg/100 g，約為對照組的1/2;用含葉酸的飼料飼喂虱螨脲處理的棉鈴蟲幼蟲，能夠明顯降低其中毒癥狀。分子對接結(jié)果表明，胞質(zhì)10-甲酰基四氫葉酸脫氫酶（cytosolic 10-formyltetrahydrofolate dehydrogenase）與虱螨脲共形成12個(gè)氫鍵和6個(gè)鹵鍵，具有最低的結(jié)合能（-10.4 kcal/mol），推測其是虱螨脲的潛在作用靶標(biāo)之一。研究結(jié)果進(jìn)一步豐富了虱螨脲的作用機(jī)理，為以葉酸途徑為靶標(biāo)的殺蟲劑開發(fā)提供了理論指導(dǎo)。

關(guān)鍵詞 虱螨脲;?棉鈴蟲;?葉酸一碳庫;?轉(zhuǎn)錄組;?分子對接

中圖分類號： S 482.37; Q 965.9; TP 391

文獻(xiàn)標(biāo)識碼：?A

DOI：?10.16688/j.zwbh.2021626

Abstract In order to explore the action mechanism of lufenuron on Helicoverpa armigera， the 3rd instar larvae of H.armigera were treated with lufenuron （200 mg/L）. The effects of lufenuron on folate and one-carbon metabolism of H.armigera were analyzed by transcriptome sequencing， folate content determination， folate feeding test and molecular docking technology. The number of differentially expressed genes （DEGs） increased with the extension of time after treatment， and their functions were mainly enriched in one-carbon pool by folate， oxidation-reduction process， chitin binding， amino acid synthesis， purine synthesis and metabolism， and ion transport across membranes. Among them， the one-carbon pool by folate pathway （P=6.35×10-5） was extremely significantly enriched in 0.5 h after treatment. The folate content in H.armigera larvae decreased significantly after treated with lufenuron， and the folate content of H.armigera treated with lufenuron was 4.53 μg/100 g after 24 h， about 1/2 of the control. The poisoning symptoms of H.armigera treated with lufenuron were obviously reduced by feeding H.armigera with folate. Molecular docking results showed that cytosolic 10-formyltetrahydrofolate dehydrogenase formed 12 hydrogen bonds and 6 halogen bonds with lufenuron. They had the lowest binding energy （-10.4 kcal/mol）， indicating that they were the potential targets of lufenuron. Our study helps to understand the action mechanism of lufenuron， and provides a theoretical guidance for the development of insecticides targeting folate pathway.

Key words lufenuron;?Helicoverpa armigera;?one-carbon pool by folate;?transcriptome;?molecular docking

虱螨脲（lufenuron）是一種選擇性好，對環(huán)境和非靶標(biāo)生物安全的苯甲酰脲類昆蟲生長調(diào)節(jié)劑，具有觸殺和胃毒作用。對鱗翅目［1］、雙翅目［2］、鞘翅目［3］以及虱目［4］等害蟲有出色的防效，適于防治對擬除蟲菊酯和有機(jī)磷農(nóng)藥產(chǎn)生抗性的害蟲，廣泛用于棉花、玉米、甜菜、馬鈴薯、葡萄、柑橘等作物。一般認(rèn)為，虱螨脲干擾鱗翅目昆蟲的正常生長發(fā)育，但其具體的作用機(jī)理卻一直沒有定論。

關(guān)于虱螨脲作用機(jī)理的研究較多。如Dean等發(fā)現(xiàn)虱螨脲引起貓櫛頭蚤Ctenocephalides felis表皮細(xì)胞中細(xì)胞質(zhì)減少，線粒體、核糖體、高爾基體等細(xì)胞器溶解，中腸細(xì)胞的分化被抑制［5］。Douris等認(rèn)為虱螨脲直接作用于黑腹果蠅Drosophila melanogaster的幾丁質(zhì)合成酶［6］。此外，虱螨脲能夠抑制卵幾丁質(zhì)的合成，進(jìn)而影響卵的呼吸代謝及胚胎發(fā)育過程中DNA和蛋白質(zhì)代謝，致使卵內(nèi)胚胎缺乏幾丁質(zhì)而不能孵化或孵化后隨即死亡［7-10］。以上研究表明，虱螨脲的作用機(jī)理比較復(fù)雜。

棉鈴蟲Helicoverpa armigera （Hübner）是為害棉、糧、油、蔬菜等多種作物的重要害蟲，在全世界范圍內(nèi)造成了重大的經(jīng)濟(jì)損失［11］。為了進(jìn)一步研究虱螨脲的作用機(jī)理，本研究通過轉(zhuǎn)錄組測序分析虱螨脲對棉鈴蟲3齡幼蟲代謝途徑的影響，在此基礎(chǔ)上利用鼠李糖桿菌法測定虱螨脲處理后不同時(shí)間棉鈴蟲體內(nèi)葉酸含量的變化，并使用葉酸飼喂虱螨脲處理后的棉鈴蟲，觀察棉鈴蟲癥狀表現(xiàn)。進(jìn)一步結(jié)合分子對接技術(shù)探究虱螨脲與棉鈴蟲葉酸一碳庫代謝途徑相關(guān)酶的作用關(guān)系，旨在深入地了解虱螨脲殺蟲作用機(jī)理，為棉鈴蟲的可持續(xù)防控提供指導(dǎo)。

1?材料與方法

1.1?供試材料

供試?yán)ハx：棉鈴蟲H.armigera幼蟲，購于河南省濟(jì)源白云實(shí)業(yè)有限公司。

供試試劑：97%虱螨脲原藥，武漢遠(yuǎn)成共創(chuàng)科技有限公司;N，N-二甲基甲酰胺和Triton X-100，天津市光復(fù)科技發(fā)展有限公司;胰蛋白酶-EDTA消化液、98%葉酸，索萊寶科技有限公司。

供試儀器：昆蟲生長培養(yǎng)箱（Adaptis A1000IN），加拿大康威龍公司;實(shí)時(shí)熒光PCR儀（ABI 7500）、酶標(biāo)儀（Multiskan FC），Thermo Fisher Scientific公司。

1.2?試驗(yàn)方法

1.2.1?轉(zhuǎn)錄組樣本制備

將虱螨脲原藥溶于適量的N，N-二甲基甲酰胺中，配制成10 g/L的母液，用0.01% Triton X-100溶液稀釋成200 mg/L的供試藥液。以0.01% Triton X-100溶液為溶劑對照。采用浸蟲法處理棉鈴蟲［12］。將棉鈴蟲3齡幼蟲分別在200 mg/L虱螨脲藥液和0.01% Triton X-100中浸漬5 s，取出后用濾紙吸干多余藥液，轉(zhuǎn)入昆蟲生長培養(yǎng)箱，在溫度（26±1）℃，相對濕度（70±10）%，光周期L∥D=14 h∥10 h條件下單頭飼養(yǎng)并接入人工飼料。浸蟲時(shí)棉鈴蟲幼蟲口器保持在液面以上，減少胃毒作用干擾。分別取處理后0、0.5、1、3、6、12 h和24 h的棉鈴蟲幼蟲，在冰上放置2 min，用大頭針固定蟲體于蠟盤上，用鑷子、剪刀取下表皮組織，用預(yù)冷的0.9% NaCl溶液清洗干凈，-78°C保存后送至北京諾禾致源科技股份有限公司進(jìn)行轉(zhuǎn)錄組測序。每個(gè)處理重復(fù)3次，每重復(fù)30頭試蟲。

1.2.2?轉(zhuǎn)錄組數(shù)據(jù)分析

測序獲得的所有clean reads與棉鈴蟲參考基因組 （Harm_1.0）［13］比對，使用featureCounts計(jì)算映射到每個(gè)基因的reads數(shù)。根據(jù)基因的長度計(jì)算每個(gè)基因的FPKM （fragments per kilobase of transcript per million mapped fragments）。使用DESeq2 1.20.0進(jìn)行虱螨脲各處理組與CK組兩兩組合之間基因的差異表達(dá)分析，通過clusterProfiler 3.4.4軟件實(shí)現(xiàn)差異表達(dá)基因的GO富集分析和KEGG通路分析。

1.2.3?差異基因的篩選和RT-qPCR驗(yàn)證

由轉(zhuǎn)錄組分析得到差異表達(dá)基因后，將|log2（Foldchange）|>0 且Padj≤0.05 作為差異表達(dá)基因的篩選標(biāo)準(zhǔn)［14］。選取差異表達(dá)基因中與葉酸一碳庫代謝途徑相關(guān)的4個(gè)基因進(jìn)行轉(zhuǎn)錄組數(shù)據(jù)的驗(yàn)證，從NCBI網(wǎng)站下載目的基因序列，使用Primer-BLAST模塊設(shè)計(jì)基因引物（表1）。將棉鈴蟲RNA反轉(zhuǎn)錄為cDNA。qPCR按照ChamQ Universal SYBR qPCR Master Mix熒光定量試劑盒（南京諾唯贊生物科技股份有限公司）步驟進(jìn)行。熒光定量PCR反應(yīng)體系（10 μL）：10 mmol/L上、下游引物各0.5 μL，2×ChamQ Universal SYBR qPCR Master Mix 5 μL，RNase-Free ddH2O 3 μL、cDNA模板 1 μL。反應(yīng)程序?yàn)椋?5℃預(yù)變性10 min；95℃變性15 s，60℃延伸1 min，循環(huán)40次。以棉鈴蟲的β-actin為內(nèi)參基因。每個(gè)樣本設(shè)置3次重復(fù)。2-ΔΔCt法計(jì)算基因的相對表達(dá)量。

1.2.4?棉鈴蟲葉酸含量測定

按照1.2.1的方法處理試蟲，分別取處理后0.5、1、3、6 h和24 h的棉鈴蟲整蟲液氮研磨制備樣品。分別取1 g制備的棉鈴蟲樣品和 20 mg 豬胰酶（豬胰酶P1750，Sigma-Aldrich公司）加入50 mL離心管中，加入30 mL磷酸鹽緩沖液，搖勻后37℃黑暗條件下孵育4 h;95℃水浴30 min，12 500 r/min離心5 min，吸取上清液過0.2 μm 無菌濾膜后備用。葉酸含量測定和葉酸標(biāo)準(zhǔn)曲線的制備方法按照試劑盒（葉酸試劑盒P1001，德國VitaFast公司）說明書進(jìn)行，葉酸含量計(jì)算公式：葉酸（μg/100 g）=吸光值對應(yīng)標(biāo)準(zhǔn)曲線的濃度×稀釋倍數(shù)/樣品重量。

1.2.5?葉酸飼喂試驗(yàn)

選取發(fā)育狀態(tài)良好的棉鈴蟲3齡幼蟲，饑餓處理24 h后，按照1.2.1的方法處理試蟲。處理后分別飼喂葉酸含量為0和0.1%（質(zhì)量百分比）的飼料0.3 g/頭，每濃度重復(fù)3次，每重復(fù)12頭試蟲。待試蟲取食完含葉酸飼料后繼續(xù)飼喂正常飼料，藥后第7天記錄試蟲癥狀。

1.2.6?分子對接

通過轉(zhuǎn)錄組測序結(jié)果，結(jié)合葉酸含量測定和飼喂試驗(yàn)，挑選了4個(gè)在虱螨脲處理后下調(diào)或連續(xù)下調(diào)表達(dá)的基因（表2），從NCBI（https：∥www.ncbi.nlm.nih.gov/）下載目的基因編碼的氨基酸序列，導(dǎo)入在線服務(wù)器Swiss-model （https：∥swissmodel.expasy.org/）進(jìn)行同源建模，再將同源建模結(jié)果上傳到SAVES 6.0 （https：∥saves.mbi.ucla.edu） 在線服務(wù)器，使用Verify3D和拉氏圖（PROCHECK）模塊對蛋白模型進(jìn)行結(jié)構(gòu)評價(jià)。使用AutoDock Vina軟件進(jìn)行分子對接，蛋白質(zhì)模型和虱螨脲加氫、去水、加電荷等操作，對接參數(shù)exhaustiveness=24進(jìn)行精細(xì)對接。用Discovery Studio 2021 client進(jìn)行虱螨脲與蛋白質(zhì)相互作用分析，使用PyMOL 2.3軟件繪制虱螨脲和蛋白的結(jié)合圖［15］。

2?結(jié)果與分析

2.1?轉(zhuǎn)錄組結(jié)果分析

2.1.1?差異表達(dá)基因分析和GO富集分析

相比于對照，虱螨脲處理棉鈴蟲幼蟲后0 h共有1 151個(gè)基因差異表達(dá)，其中上調(diào)表達(dá)基因400個(gè)，下調(diào)表達(dá)基因751個(gè);處理后0.5 h共有1 053個(gè)基因差異表達(dá)，其中上調(diào)表達(dá)基因364個(gè)，下調(diào)表達(dá)基因689個(gè);處理后1 h共有1 573個(gè)基因差異表達(dá)，其中上調(diào)表達(dá)基因691個(gè)，下調(diào)表達(dá)基因882個(gè);處理后 3 h共有1 486個(gè)基因差異表達(dá)，其中上調(diào)表達(dá)基因677個(gè)，下調(diào)表達(dá)基因809個(gè);處理后6 h共有939個(gè)基因差異表達(dá)，其中上調(diào)表達(dá)基因385個(gè)，下調(diào)表達(dá)基因554個(gè);處理后12 h共有1 685個(gè)基因差異表達(dá)，其中上調(diào)表達(dá)基因843個(gè)，下調(diào)表達(dá)基因842個(gè);處理后24 h共有2 174個(gè)基因差異表達(dá)，其中上調(diào)表達(dá)基因1 167個(gè)，下調(diào)表達(dá)基因1 007個(gè)。

GO分析篩選出69個(gè)有顯著差異的分類條目（Padj＜0.05），主要?dú)w屬于生物過程、細(xì)胞組分和分子功能。在生物過程上，氧化還原過程、碳水化合物衍生物代謝、幾丁質(zhì)代謝和藥物代謝條目富集基因較多，分別有314、109、74個(gè)和72個(gè)基因。

藥物代謝過程涉及外來化學(xué)物質(zhì)如藥物在蟲體內(nèi)的化學(xué)反應(yīng)和途徑，氧化還原過程涉及一個(gè)或多個(gè)電子與質(zhì)子的轉(zhuǎn)移。幾丁質(zhì)屬于高分子多糖類物質(zhì)，是昆蟲表皮結(jié)構(gòu)的重要組成成分，幾丁質(zhì)代謝過程涉及氧化還原過程和碳水化合物衍生物代謝。生物過程的富集結(jié)果表明虱螨脲處理后，棉鈴蟲幾丁質(zhì)代謝出現(xiàn)異常，虱螨脲可能干擾了氧化還原過程。

在細(xì)胞組分上，細(xì)胞外區(qū)域條目被富集多次，表明細(xì)胞外區(qū)域的正常生化反應(yīng)受到影響。在分子功能上，結(jié)構(gòu)分子活性、表皮結(jié)構(gòu)成分和輔酶結(jié)合條目富集到基因較多，表明試蟲表皮結(jié)構(gòu)物質(zhì)的代謝受到影響。

2.1.2?差異表達(dá)基因的KEGG富集

KEGG分析發(fā)現(xiàn)葉酸一碳庫途徑在每個(gè)處理組中都有顯著富集，其中0.5 h時(shí)極顯著富集（P=6.35×10-5）。有1 053個(gè)差異基因被富集到93個(gè)代謝通路，共獲得12條顯著性代謝通路（表3），包括葉酸一碳庫、DNA復(fù)制、錯(cuò)配修復(fù)、嘌呤代謝、酪氨酸代謝、苯丙氨酸代謝、氨基酸生物合成、N-多糖生物合成、乙醛酸二羧酸代謝、壽命調(diào)節(jié)通路以及昆蟲激素生物合成，其中葉酸一碳庫代謝通路最顯著。葉酸一碳庫主要與葉酸及其衍生物的一碳單位代謝過程相關(guān)，在DNA復(fù)制、核酸甲基化修飾和一碳單元轉(zhuǎn)運(yùn)等方面發(fā)揮重要作用［16］。將富集到的葉酸一碳庫的差異基因進(jìn)行組間比較時(shí)發(fā)現(xiàn)，在所有處理組中參與葉酸一碳庫途徑的差異表達(dá)基因共有10個(gè)，這些基因在虱螨脲處理后呈現(xiàn)不同的表達(dá)模式（圖1）。三功能嘌呤生物合成蛋白腺苷-3 （trifunctional purine biosynthetic protein adenosine-3）、雙功能嘌呤生物合成蛋白水解酶（bifunctional purine biosynthesis protein PURH）和胞質(zhì)絲氨酸羥甲基轉(zhuǎn)移酶（cytosolic serine hydroxy methyltransferase）基因都連續(xù)上調(diào)表達(dá)，氨甲基轉(zhuǎn)移酶（aminomethyltransferase）、二氫葉酸還原酶（dihydrofolate reductase）、胸苷酸合酶（thymidylate synthase）和C-1-四氫葉酸合酶（C-1-tetrahydrofolate synthase）基因都先上調(diào)表達(dá)，然后下調(diào)表達(dá)，最后又上調(diào)表達(dá);5-甲酰四氫葉酸環(huán)連接酶（5-formyltetrahydrofolate cyclo-ligase）基因的表達(dá)先下調(diào)，然后上調(diào)，最后又下調(diào)表達(dá);雙功能亞甲基四氫葉酸脫氫酶 （bifunctional methylene tetrahydrofolate dehydrogenase）和胞質(zhì)10-甲?；臍淙~酸脫氫酶（cytosolic 10-formyltetrahydrofolate dehydrogenase）基因表達(dá)呈現(xiàn)持續(xù)的下調(diào)表達(dá)模式。

2.2?差異基因的RT-qPCR驗(yàn)證

KEGG通路分析表明，葉酸一碳庫代謝途徑相關(guān)基因在虱螨脲處理的棉鈴蟲轉(zhuǎn)錄本中顯著的差異表達(dá)（圖2）。cytosolic 10-formyltetra-hydrofolate dehydrogenase基因（FDH）是葉酸代謝酶基因，主要功能是調(diào)節(jié)葉酸量和葉酸結(jié)合甲基量以及參與嘌呤從頭合成等［17］。5-formyltetrahydrofolate cyclo-ligase基因（ALDH1L1）是唯一的將5-甲酰四氫葉酸催化成5，10-甲酰四氫葉酸的酶基因［18］。bifunctional methylenetetrahydrofolate dehydrogenase基因（FTHFDH）是一類依賴于NADP，主要催化四氫葉酸衍生物之間相互轉(zhuǎn)化的基因［19］。C-1-tetrahydrofolate synthase基因（FTHFS）是葉酸一碳庫途徑中的中心酶基因，將5，10-甲酰四氫葉酸催化成5，10-次甲基四氫葉酸［20］。虱螨脲處理后，棉鈴蟲3齡幼蟲的這4個(gè)酶的基因大致呈現(xiàn)先上調(diào)表達(dá)后下調(diào)表達(dá)的趨勢，其中HaFDH（圖2a）的FPKM值變化趨勢和RT-qPCR一致，在處理后24 h基因表達(dá)量顯著下降;HaALDH1L1 （圖2b）也具有類似的規(guī)律。HaFTHFS （圖2c）基因FPKM值和RT-qPCR試驗(yàn)所測值在處理后6 h和12 h處表現(xiàn)不一致，但整體趨勢兩者相似。 HaFTHFDH （圖2d）在處理后0～12 h FPKM值變化趨勢和RT-qPCR試驗(yàn)所測值相同，但在處理后24 h基因FPKM值增大而RT-qPCR試驗(yàn)所測值變化不大。

2.3?虱螨脲對棉鈴蟲葉酸含量的影響

為了進(jìn)一步驗(yàn)證虱螨脲對棉鈴蟲葉酸代謝的影響，采用鼠李糖桿菌法測定了虱螨脲處理后棉鈴蟲幼蟲體內(nèi)的葉酸含量（葉酸的標(biāo)準(zhǔn)曲線方程為y=0.547 2x-0.000 2，R2=0.998 7）。虱螨脲處理棉鈴蟲后，蟲體內(nèi)的葉酸含量隨著處理后時(shí)間的延長而下降，處理后6 h和24 h葉酸含量分別為4.27 μg/100 g和4.53 μg/100 g，約是CK組葉酸含量的一半;除處理后0.5 h的葉酸含量高于CK外，其余時(shí)間點(diǎn)都顯著低于CK（圖3）。

虱螨脲處理后棉鈴蟲表現(xiàn)出典型的“雙頭囊”癥狀，部分蟲體出現(xiàn)蛻皮困難，新表皮和舊表皮同時(shí)存在，表皮顏色變黑，毛片突出（圖4a）。虱螨脲處理過的棉鈴蟲再飼喂含葉酸飼料后，未出現(xiàn)雙頭囊癥狀，中毒癥狀明顯減輕，表明葉酸對虱螨脲處理的棉鈴蟲具有治療效果（圖4b）。

2.4?分子對接

以上研究結(jié)果表明，虱螨脲可干擾棉鈴蟲葉酸一碳庫途徑基因的表達(dá)，并顯著影響棉鈴蟲體內(nèi)的葉酸含量，為了進(jìn)一步探究虱螨脲影響棉鈴蟲葉酸一碳庫代謝途徑的作用機(jī)理，使用分子對接技術(shù)探究虱螨脲處理后葉酸一碳庫途徑的4個(gè)基因所編碼蛋白與虱螨脲的關(guān)系。結(jié)果表明，虱螨脲對5-formyltetrahydrofolate cyclo-ligase、cytosolic 10-formyltetrahydrofolate dehydrogenase、C-1-tetrahydrofolate synthase、bifunctional methylenetetrahydrofolate dehydrogenase蛋白的結(jié)合能都低于-8 kcal/mol，其中對cytosolic 10-formyltetrahydrofolate dehydrogenase結(jié)合能最低，為-10.4 kcal/mol。結(jié)合能越小表示虱螨脲與蛋白質(zhì)結(jié)合程度越強(qiáng)，故推測cytosolic 10-formyltetrahydrofolate dehydrogenase和虱螨脲的結(jié)合程度強(qiáng)，可能是虱螨脲的潛在作用靶標(biāo)之一。

結(jié)合虱螨脲與蛋白模型組成復(fù)合體時(shí)的相互作用關(guān)系分析，cytosolic 10-formyltetrahydrofolate dehydrogenase和虱螨脲之間的氫鍵和鹵鍵數(shù)最多，共形成12個(gè)氫鍵和6個(gè)鹵鍵（圖5a）。而bifunctional methylenetetrahydrofolate dehydrogenase與虱螨脲只形成6個(gè)氫鍵，不存在鹵鍵（圖5b）。5-formyltetrahydrofolate cyclo-ligase和虱螨脲之間存在6個(gè)氫鍵和6個(gè)鹵鍵（圖5c）。C-1-tetrahydro-folate synthase與虱螨脲形成3個(gè)氫鍵和4個(gè)鹵鍵（圖5d）。由于虱螨脲特殊結(jié)構(gòu)以及氮原子和氯原子所處的特殊位置關(guān)系賦予虱螨脲影響棉鈴蟲葉酸一碳庫途徑中基因表達(dá)和葉酸含量的生物活性。

3?結(jié)論與討論

虱螨脲處理后，棉鈴蟲幼蟲的差異基因GO富集主要集中在結(jié)合蛋白活性、催化活性和代謝過程，這與虱螨脲作用于甜菜夜蛾Spodoptera exigua Hübner的結(jié)果類似［21］。虱螨脲處理貓櫛首蚤后，其表皮細(xì)胞發(fā)生降解，細(xì)胞器數(shù)量顯著少于對照組［5］，表明虱螨脲對DNA復(fù)制過程也產(chǎn)生了抑制作用，阻止表皮細(xì)胞正常增殖。轉(zhuǎn)錄組分析表明，虱螨脲處理棉鈴蟲后幾丁質(zhì)代謝、表皮細(xì)胞組分相關(guān)基因顯著富集，藥物代謝途徑、嘌呤合成途徑、嘌呤核苷酸代謝過程、蛋白質(zhì)合成以及氧化還原途徑也有明顯的富集。以上結(jié)果表明虱螨脲通過多種途徑影響棉鈴蟲的生長發(fā)育。

葉酸一碳庫代謝為昆蟲的生命活動提供一碳單位，發(fā)揮一碳供體和受體的作用［22］，參與嘌呤和胸腺嘧啶的合成以及DNA、RNA和蛋白質(zhì)的甲基化過程［23］。本研究表明，虱螨脲顯著影響棉鈴蟲葉酸一碳庫代謝途徑相關(guān)基因的表達(dá)，虱螨脲處理后棉鈴蟲幼蟲的葉酸含量顯著下降，飼喂含有葉酸的飼料可明顯緩解棉鈴蟲的中毒癥狀。進(jìn)一步的分子對接結(jié)果表明，虱螨脲與5-formyltetrahydrofolate cyclo-ligase、cytosolic 10-formyltetrahydrofolate dehydrogenase、C-1-tetrahydrofolate synthase和bifunctional methylenetetrahydrofolate dehydrogenase等4種酶都有較好的結(jié)合作用。研究表明，通過干擾鱗翅目害蟲的葉酸一碳庫代謝途徑可達(dá)到殺蟲的目的，如葉酸拮抗劑2，4-二氨基嘧啶類化合物對多種鱗翅目害蟲都有較好的防治效果，對美洲棉鈴蟲H.zea （Boddie） LC50為0.5 mg/ L，對煙芽夜蛾Heliothis virescens的LC50為0.2 mg/ L［24］。氨甲蝶呤可通過抑制二氫葉酸還原酶（DHFR）的活性影響果蠅的生殖和后代的生長發(fā)育過程，并造成果蠅的生長發(fā)育異常以及畸形等癥狀［25］。

綜上所述，本研究通過轉(zhuǎn)錄組測序首次分析出葉酸一碳庫代謝途徑為虱螨脲作用于棉鈴蟲的潛在途徑之一，葉酸含量測定和飼喂葉酸試驗(yàn)進(jìn)一步證實(shí)了該機(jī)理，并通過同源模建和分子對接技術(shù)預(yù)測了虱螨脲與葉酸一碳庫代謝途徑相關(guān)蛋白的結(jié)合模式，研究結(jié)果可為研究虱螨脲的作用機(jī)理提供線索，為以葉酸途徑為靶標(biāo)的殺蟲劑開發(fā)提供理論指導(dǎo)。

參考文獻(xiàn)

［1］?CORREIA A A， WANDERLEY-TEIXEIRA V， TEIXEIRA A A C， et al. Microscopic analysis of Spodoptera frugiperda （Lepidoptera： Noctuidae） embryonic development before and after treatment with azadirachtin， lufenuron， and deltamethrin ［J］. Journal of Economic Entomology， 2013， 106（2）： 747-755.

［2］?SAMPSON B J， MARSHALL D A， SMITH B J， et al. Erythritol and lufenuron detrimentally alter age structure of wild Drosophila suzukii （Diptera： Drosophilidae） populations in blueberry and blackberry ［J］. Journal of Economic Entomology， 2017， 110（2）： 530-534.

［3］?COSTA H N， CUNHA F M D， CRUZ G S， et al. Lufenuron impact upon Anthonomus grandis Boheman （Coleoptera： Curculionidae） midgut and its reflection in gametogenesis ［J］. Pesticide Biochemistry and Physiology， 2017， 137： 71-80.

［4］?DAVIS R M， CLEUGH E， SMITH R T， et al. Use of a chitin synthesis inhibitor to control fleas on wild rodents important in the maintenance of plague， Yersinia pestis， in California ［J］. Journal of Vector Ecology， 2008， 33（2）： 278-284.

［5］?DEAN S R， MEOLA R W， MEOLA S M， et al. Mode of action of lufenuron in adult Ctenocephalides felis （Siphonaptera： Pulicidae） ［J］. Journal of Medical Entomology， 1999， 36（4）： 486-492.

［6］?DOURIS V， STEINBACH D， PANTELERI R， et al. Resistance mutation conserved between insects and mites unravels the benzoylurea insecticide mode of action on chitin biosynthesis ［J］. Proceedings of the National Academy of Sciences of the United States of America， 2016， 113（51）： 14692-14697.

［7］?MANSUR J F， FIGUEIRA-MANSUR J， SANTOS A S， et al. The effect of lufenuron， a chitin synthesis inhibitor， on oogenesis of Rhodnius prolixus ［J］. Pesticide Biochemistry and Physiology， 2010， 98（1）： 59-67.

［8］?賈變桃，盧晶晶，袁嘉偉，等. 虱螨脲對小菜蛾成蟲繁殖力及子代種群發(fā)展的影響［J］. 植物保護(hù)學(xué)報(bào)， 2019， 46（3）： 582-588.

［9］?MATSUMURA F. Studies on the action mechanism of benzoylurea insecticides to inhibit the process of chitin synthesis in insects： A review on the status of research activities in the past， the present and the future prospects ［J］. Pesticide Biochemistry and Physiology， 2010， 97（2）： 133-139.

［10］GANGISHETTI U， BREITENBACH S， ZANDER M， et al. Effects of benzoylphenylurea on chitin synthesis and orientation in the cuticle of the Drosophila larva ［J］. European Journal of Cell Biology， 2009， 88（3）： 167-180.

［11］LU Yanhui， WU Kongming， JIANG Yuying， et al. Widespread adoption of Bt cotton and insecticide decrease promotes biocontrol services ［J］. Nature， 2012， 487（7407）： 362-365.

［12］MEYER F， FLOTENMEYER M， MOUSSIAN B. The sulfonylurea receptor sur is dispensable for chitin synthesis in Drosophila melanogaster embryos ［J］. Pest Management Science， 2013， 69（10）： 1136-1140.

［13］PEARCE S L， CLARKE D F， EAST P D， et al. Genomic innovations， transcriptional plasticity and gene loss underlying the evolution and divergence of two highly polyphagous and invasive Helicoverpa pest species ［J/OL］. BMC Biology， 2017， 15（1）： 63.DOI：10.1186/s12915-017-0402-6.

［14］YOUNG M D， WAKEFIELD M J， SMYTH G K， et al. Gene ontology analysis for RNA-seq： accounting for selection bias ［J/OL］. Genome Biology， 2010， 11（2）： R14. DOI：10.1186/gb-2010-11-2-r14.

［15］TROTT O， OLSON A J. Autodock vina： improving the speed and accuracy of docking with a new scoring function， efficient optimization， and multithreading ［J］. Journal of Computational Chemistry， 2009， 31（2）： 455-461.

［16］DUCKER G S， RABINOWITZ J D. One-carbon metabolism in health and disease ［J］. Cell Metabolism， 2017， 25（1）： 27-42.

［17］KRUPENKO N I， SHARMA J， PEDIADITAKIS P， et al. Cytosolic 10-formyltetrahydrofolate dehydrogenase regulates glycine metabolism in mouse liver ［J/OL］. Scientific Reports， 2019， 9（1）： 14937. DOI： 10.1038/s41598-019-51397-1.

［18］MEIER C， CARTER L G， WINTER G， et al. Structure of 5-formyltetrahydrofolate cyclo-ligase from Bacillus anthracis （BA4489） ［J］. Acta Crystallographica Section F Structural Biology and Crystallization Communications， 2007， 63（3）： 168-172.

［19］LEE W H， SUNG M W， KIM J H， et al. Crystal structure of bifunctional 5，10-methylenetetrahydrofolate dehydrogenase/cyclohydrolase from Thermoplasma acidophilum ［J］. Biochemical and Biophysical Research Communications， 2011， 406（3）： 459-463.

［20］CHEN Meijun， ZHAI Jingjie， LIU Yu， et al. Molecular cloning and characterization of C1 tetrahydrofolate （C1-THF） synthase in Bombyx mori， silkworm ［J］. Gene， 2018， 663： 25-33.

［21］陳吉祥， 任相亮， 姜偉麗， 等. 甲氧蟲酰肼與3種藥劑對甜菜夜蛾的聯(lián)合毒力及亞致死效應(yīng)［J］. 中國棉花， 2019， 46（2）： 4-8.

［22］FU T F， MARAS B， BARRA D， et al. A noncatalytic tetrahydrofolate tight binding site is on the small domain of 10-formyltetrahydrofolate dehydrogenase ［J］. Archives of Biochemistry and Biophysics， 1999， 367（2）： 161-166.

［23］KIM D W， HUANG T， SCHIRCH D， et al. Properties of tetrahydropteroylpentaglutamate bound to 10-formyltetrahydrofolate dehydrogenase ［J］. Biochemistry， 1996， 35（49）： 15772-15783.

［24］WIERENGA J M， CULLEN T G， DYBAS J A， et al. Insecticidal activity and mode of action of 2，4-diaminopyrimidines ［J］. Pest Management Science， 2000， 3（56）： 233-236.

［25］AFFLECK J G， WALKER V K. Transgenic rescue of methotrexate-induced teratogenicity in Drosophila melanogaster ［J］. Toxicological Sciences， 2007， 99（2）： 522-531.

（責(zé)任編輯：楊明麗）