響應槲皮素誘導的美國白蛾P450基因及CYP6ZB2功能分析

邵紅揚，孟香，張濤，陳敏

響應槲皮素誘導的美國白蛾P450基因及功能分析

邵紅揚，孟香，張濤，陳敏

北京林業大學林學院/林木有害生物防治北京市重點實驗室，北京 100083

【】明確美國白蛾（）響應槲皮素誘導的差異表達P450基因及系統發育關系，探究美國白蛾P450基因對槲皮素誘導的劑量效應，以及P450基因在美國白蛾對槲皮素的解毒適應中的作用。【】基于槲皮素誘導的美國白蛾中腸轉錄組數據，從中篩選顯著差異表達的P450基因，應用MEGA X軟件對5個P450 基因及其他近緣物種的P450基因進行系統進化關系分析；通過RT-qPCR技術檢測5個P450基因在不同槲皮素濃度（0.5%、1.0%、2.0%和4.0%）處理下的表達模式；克隆全長cDNA序列，利用RNA干擾（RNAi）技術探明基因沉默對美國白蛾體重及存活率的影響。【】從轉錄組的差異表達基因中篩選到5個顯著響應槲皮素誘導的P450基因，經P450命名委員會分別被命名為、、、和（GenBank登錄號：ON032363—ON032367）。系統發育分析表明，這5個基因被聚類到、、和線粒體CYP進化枝。RT-qPCR結果表明，P450基因響應槲皮素脅迫具有顯著的劑量效應，、、和的表達均被不同濃度槲皮素抑制，的表達在0.5%—2.0%的槲皮素誘導下顯著上調，而在4.0%的槲皮素濃度下被顯著抑制。克隆的全長cDNA序列，具有完整的開放閱讀框（ORF），編碼500個氨基酸。RNAi結果表明，注射兩種劑量（500和1 000 ng）ds均可顯著抑制的表達，注射1 000 ng ds48 h后干擾效率達到最高（96.9%）。被干擾的美國白蛾幼蟲經0.5%槲皮素處理后體重增長量低于對照組，存活率（53.33%）分別為未處理組的61.5%和注射雙蒸水組的72.7%（＜0.05）。【】篩選鑒定到美國白蛾響應槲皮素誘導的5個P450基因，槲皮素對5個基因的誘導表達具有明顯的劑量效應，說明不同的P450基因對槲皮素的響應模式存在差異。RNAi能高效沉默，0.5%槲皮素對被干擾的美國白蛾幼蟲的生長發育具有一定抑制作用，表明在美國白蛾對槲皮素的解毒適應中發揮著重要作用。

美國白蛾；槲皮素；P450基因；表達分析；RNA干擾；基因功能

0 引言

【研究意義】美國白蛾（）屬鱗翅目（Lepidoptera）目夜蛾科（Erebidae）（原燈蛾科Arctiidae），是典型的多食性食葉昆蟲，其寄主植物可達600多種[1]。由于其具有環境適應能力強、寄主種類多、危害性大等特點，被農業農村部和國家林草局列為重點管理的外來入侵物種（農業農村部公告第567號）。目前，美國白蛾在我國已擴展到14個省（直轄市）的611 個縣（區）（國家林業和草原局2022年第5號公告），且在我國南方地區擴散加速，這不僅增加了防控難度，還給我國農林業造成了巨大經濟損失[2]。槲皮素（quercetin）是黃酮類植物次生代謝物，在美國白蛾寄主植物如桑（）[3]、楊樹（）[4]等植物中廣泛存在。研究表明槲皮素對多種昆蟲具有抗蟲活性，且其抗蟲活性與濃度密切相關[5]。植食性昆蟲具有完善的解毒酶系以催化分解寄主食物中的有毒次生代謝物從而減少對自身的傷害[6]，其中細胞色素P450單加氧酶（P450）在植食性昆蟲對寄主植物次生代謝物的解毒適應中發揮了主要作用[7]。因此，研究美國白蛾響應槲皮素誘導表達的P450基因及功能，對于明確美國白蛾對寄主植物中槲皮素的適應機制具有重要意義。【前人研究進展】細胞色素P450單加氧酶（P450）是一個重要的解毒酶基因家族，主要通過羥基化、脫烷基化等方式參與昆蟲對內源（如激素等）和外源性化合物（如植物次生代謝物、殺蟲劑等）的解毒代謝[8]。在植食性昆蟲中，P450基因共分為4個進化枝：、、和線粒體CYP。進化枝由大量昆蟲P450基因組成，這些P450基因被細分為11個亞家族：、、、、、、、、、和，這些成員參與昆蟲對各種植物次生代謝物以及殺蟲劑抗藥性相關的解毒作用[9]。目前，多種農林業害蟲的P450基因被成功鑒定，如小菜蛾（）[10]、斜紋夜蛾（）[11]、草地貪夜蛾（）[11]和甜菜夜蛾（）[12]等。前期研究發現，0.5%槲皮素導致美國白蛾幼蟲的發育歷期顯著延長，幼蟲成活率、化蛹率、羽化率以及產卵量等各項指標顯著低于對照組，同時，一定濃度的槲皮素能夠提高美國白蛾體內多種解毒酶的活性，且對美國白蛾P450酶的誘導作用具有顯著的劑量效應[13]。【本研究切入點】不同昆蟲響應各種植物次生代謝物質的解毒酶基因存在差異。美國白蛾體內響應槲皮素誘導的P450基因的種類以及這些基因對槲皮素的響應模式和機制尚不清楚。【擬解決的關鍵問題】基于前期構建的0.5%槲皮素誘導的美國白蛾中腸轉錄組，從中篩選差異表達的P450基因，通過生物信息學分析，結合RT-qPCR技術，研究這些基因響應槲皮素誘導的劑量效應，進一步克隆響應槲皮素誘導的關鍵基因，并運用RNA干擾（RNA interference，RNAi）技術驗證其在槲皮素的解毒適應中的功能，為闡明美國白蛾對槲皮素的適應性機制提供科學依據。

1 材料與方法

試驗于2020年9月至2022年4月在北京林業大學林木有害生物防治北京市重點實驗室完成。

1.1 槲皮素人工飼料的配制

依據槲皮素對美國白蛾幼蟲的毒力測定[13]，以及美國白蛾喜食寄主（如桑葉）的槲皮素含量[14-15]，將槲皮素的誘導濃度范圍設置為0.5%—4.0%。槲皮素人工飼料的配制參照潘忠玉[13]的方法，以添加等量DMSO溶液的人工飼料作為對照飼料（control）。

1.2 供試昆蟲

美國白蛾幼蟲為實驗室人工飼料飼養的種群，人工飼料的配制參照曹利軍等[16]。飼養條件設定為溫度（25±1）℃，相對濕度55%—70%，光周期16L﹕8 D。

1.3 主要試劑和儀器

試劑：槲皮素（純度≥98%；源葉生物，中國上海），二甲基亞砜（DMSO）、50×TAE buffer、DL 2000 DNA Ladder、GelRed核酸染料、6×loading buffer、2×Taq PCR Master MIX（百瑞極，中國北京），瓊脂糖（貝晶生物，中國上海），RNeasy?Plus Mini Kit（Qiagen，德國），PrimeScriptTMRT reagent Kit with gDNA Eraser、TB Green?TMⅡ、PrimeSTAR?HS DNA Polymerase（TaKaRa，中國大連），TIANgel Midi Purification Kit、TIANprep Rapid Mini Plasmid Kit（天根，中國北京），?-Blunt3 Cloning Kit（全式金，中國北京），T7 RNAi Transcription Kit（諾唯贊，中國南京）。

儀器：恒溫培養箱（寧波江南儀器廠，中國），電泳儀（六一生物，中國北京），體視顯微鏡（徠卡，德國）PCR儀、CFX96型實時熒光定量PCR檢測系統（Bio-Rad，美國），凝膠成像系統（Aplegen，美國），超微量紫外分光光度計NanoDrop 8000（Thermo，美國），超凈臺（蘇州凈化，中國），臺式離心機、搖床、超低溫冰箱（Thermo，美國），普通冰箱（海爾，中國青島），微量注射器（哈美頓博納圖斯，瑞士），恒溫水浴鍋（上海精宏，中國），電子天平（Sartorius，德國）。

1.4 響應槲皮素誘導的P450基因鑒定與系統發育分析

基于課題組前期構建的槲皮素（0.5%）誘導的美國白蛾中腸轉錄組（GenBank 登錄號：PRJNA798672）以及正常飼養的美國白蛾的中腸轉錄組，從差異表達基因中篩選美國白蛾顯著響應槲皮素誘導的P450基因。通過NCBI的ORF Finder（https://www. ncbi.nlm.nih.gov/orffinder/）預測候選P450基因的開放閱讀框（open reading frame，ORF），并翻譯成氨基酸序列用于后續構建系統發育樹。通過與NCBI的NR數據庫（https://blast.ncbi.nlm.nih.gov/Blast.cgi）進行BLASTX比對，鑒定候選基因是否為P450基因。將鑒定結果為P450基因的序列信息提交至P450命名委員會（Standardized Cytochrome P450 Nomenclature Committee）進行命名。

使用MEGA X軟件將鑒定的P450基因的氨基酸序列與5種鱗翅目昆蟲（家蠶、棉鈴蟲、草地貪夜蛾、斜紋夜蛾、粉紋夜蛾）的氨基酸序列進行多序列比對。使用PhyloSuite軟件[17]中的IQ-TREE程序，采用最大似然法構建系統發育樹，各分支置信度經由bootstrap 法1 000 次循環檢驗[18-19]。使用在線工具 iTOL（https://itol.embl.de/）對發育樹進行美化。

1.5 P450基因響應槲皮素的表達分析

取剛蛻皮的美國白蛾4齡幼蟲，饑餓處理12 h后，分別轉入不同濃度的槲皮素人工飼料（0.5%、1.0%、2.0%和4.0%）和對照飼料（control），飼養24 h后收集試蟲，解剖中腸，按照RNeasy?Plus Mini Kit試劑盒的說明書提取各樣品總RNA，10頭幼蟲為一個重復，設置3個生物學重復。總RNA的完整性采用1%瓊脂糖凝膠電泳檢測，濃度和純度使用超微量紫外分光光度計NanoDrop 8000檢測。將完整性、純度和濃度均合格的總RNA按照PrimeScriptTMRT reagent Kit with gDNA Eraser試劑盒的說明書合成cDNA。使用Primer Premier 5.0軟件設計P450基因的RT-qPCR引物，和作為內參基因[20]，引物均由生工生物工程（上海）股份有限公司合成。通過10倍梯度稀釋cDNA模板的方法，計算各引物的擴增效率，RT-qPCR試驗所用引物序列見表1。根據TB Green?TMⅡ試劑盒配制RT-qPCR反應體系（25 μL）：TB GreenⅡ 12.5 μL，正反引物（10 μmol·L-1）各1 μL，cDNA模板2 μL，滅菌水8.5 μL。熱循環條件設定為95℃預變性30 s；95℃ 5 s，60℃ 30 s，循環40次；最后將RT-qPCR產物升溫至95℃制作熔解曲線，檢測是否有引物二聚體。空白對照以無菌水代替cDNA模板，用于檢測體系是否出現污染。RT-qPCR反應在CFX96型熒光定量PCR儀中運行。采用2-ΔΔCT法[21]計算基因的相對表達量。每個處理設置3個生物學重復，每個生物學重復設置3個技術重復。

1.6 CYP6ZB2的功能驗證

5個P450基因中僅的表達水平在不同濃度槲皮素下顯著升高，因此進一步對該基因進行功能分析。

1.6.1合成 使用Primer Premier 5.0軟件設計的全長引物，在全長引物的5′端加上T7啟動子序列用于合成dsRNA模板，引物均由生工生物工程（上海）股份有限公司合成，引物序列見表1。以常規人工飼料飼養的美國白蛾4齡幼蟲cDNA為模板，使用的全長引物進行PCR擴增。50mL PCR反應體系：5×PrimeSTAR Buffer 10mL，dNTP Mixture 4mL，正反向全長引物各1mL，cDNA 2mL，PrimeSTAR HS DNA Polymerase 0.5mL，ddH2O 31.5mL。PCR運行程序設置為98℃ 10 s，55℃ 5 s，72℃ 90 s，30個循環。PCR產物經1%的瓊脂糖凝膠電泳驗證后，使用TIANgel Midi Purification Kit進行回收純化，將純化產物與?-Blunt3 Cloning Kit克隆載體連接，然后轉化到-T1 Phage Resistant化學感受細胞，經過藍白斑篩選后，將陽性單克隆菌液送至北京擎科生物科技有限公司進行測序，將測序正確的單克隆菌液擴大培養，按照TIANprep Rapid Mini Plasmid Kit說明書抽提菌液質粒。

以上一步獲得的質粒為模板，使用添加T7啟動子序列的全長引物進行PCR擴增，反應體系和程序設置同上。PCR產物經1%的瓊脂糖凝膠電泳檢測后，進行回收純化，純化產物即為合成dsRNA的模板，按照T7 RNAi Transcription Kit的說明書合成dsRNA，將初步得到dsRNA采用磁珠法進行純化，純化的dsRNA經電泳驗證和濃度測定，稀釋至500和1 000 ng·mL-1保存于-20℃冰箱用于后續注射試驗。

表1 RT-qPCR、基因克隆和dsRNA合成所用引物

下劃線部分為T7啟動子序列The underlined letters represent T7 promoter sequence

1.6.2 RNAi驗證的功能 為了驗證是否參與美國白蛾對槲皮素的解毒代謝過程，根據預試驗結果，選擇500和1 000 ng作為dsRNA的注射劑量。挑選剛蛻皮1 d的美國白蛾4齡幼蟲，使用微量注射器吸取500和1 000 ng的dsRNA分別注射到幼蟲腹部倒數第二個節間膜位置，以注射等體積的雙蒸水（ddH2O）作為對照（control）。每種劑量注射30頭幼蟲為一個重復，設置3個重復。將注射后的試蟲在不含槲皮素的人工飼料中飼養48和96 h后，通過RT-qPCR檢測的沉默效率。RT-qPCR方法同1.5。將注射1 000 ng ds2的美國白蛾轉入含有0.5%槲皮素人工飼料中飼養，統計試蟲10 d內的體重增長量和存活率，以未處理（CK）和注射等體積雙蒸水（ddH2O）的幼蟲作為對照，每30頭為一個重復，設置3個生物學重復。

1.7 數據分析

試驗數據均以平均值±標準誤表示，采用Excel統計分析數據，使用SPSS 24.0中的單因素方差分析不同濃度槲皮素誘導P450基因的相對表達量和沉默效率的差異顯著性，以及RNAi后美國白蛾幼蟲存活率的差異顯著性，并使用Duncan氏新復極差檢驗法進行多重比較（＜0.05）。使用獨立樣本檢驗分析轉錄組中P450基因的差異表達驗證結果，使用重復檢測方差分析RNAi后美國白蛾幼蟲的平均體重增長量和差異顯著性（＜0.05）。采用GraphPad Prism 8.0和Origin 2017作圖。

2 結果

2.1 響應槲皮素誘導的P450基因的篩選與鑒定

從0.5%槲皮素誘導的美國白蛾中腸轉錄組數據中，共篩選出5個顯著響應槲皮素誘導的差異表達P450基因，經BLASTX鑒定，除（GenBank登錄號：BBD13397.1）基因外，其他4個基因均為美國白蛾物種中新鑒定的P450基因，鑒定結果見表2。轉錄組的差異表達分析結果顯示，在0.5%的槲皮素誘導下，、、和的表達水平顯著降低；而的表達水平升至對照的2.69倍。5個P450基因中僅具有完整的ORF，編碼500個氨基酸，其余基因均缺失3′端，編碼蛋白范圍為146—507個氨基酸。P450命名委員會將5個基因分別命名為、、、和。5個P450基因的序列均已提交至NCBI數據庫，GenBank登錄號為ON032363—ON032367。

表2 響應槲皮素誘導的P450基因的BLASTX結果

Down：表達下調down-regulated expression；Up：表達上調up-regulated expression

2.2 P450基因的系統發育分析

美國白蛾響應槲皮素誘導的5個P450基因的系統發育樹見圖1。結果表明，5個基因分別聚類到P450家族的4個進化枝，其中，和聚類到進化枝；聚類到線粒體CYP進化枝；聚類到進化枝；聚類到進化枝的亞家族。從系統發育樹可知，美國白蛾與棉鈴蟲和斜紋夜蛾序列同源性較高，美國白蛾與棉鈴蟲聚為一支，序列同源性較高。

Bmor：家蠶B. mori；Harm：棉鈴蟲 H. armigera；Sfru：草地貪夜蛾 S. frugiperda；Slit：斜紋夜蛾S. litura；Tni：粉紋夜蛾 T. ni。不同顏色代表不同P450基因進化枝；紅色字體基因代表本研究鑒定的5個P450基因；分支上的數字表示1000次重復自舉檢驗值Different colors represent different cytochrome P450 gene clades. The red font genes represent the five cytochrome P450 genes identified in this study. The number at tree node represents bootstrap value from 1000 replicates

2.3 槲皮素誘導5個P450基因的劑量效應

RT-qPCR結果顯示，、和的表達均被不同濃度槲皮素抑制，而較低濃度（0.5%—2.0%）的槲皮素誘導的表達顯著上調。、和的表達隨著槲皮素濃度升高，抑制作用逐漸減弱，3個基因均在0.5%槲皮素濃度下表達水平降至最低，分別為對照組的26%（＜0.0001，圖2-A）、32%（＜0.001，圖2-B）和10%（＜0.001，圖2-C）。的表達隨著槲皮素濃度的增高，抑制作用逐漸加強但差異并不顯著，在4.0%槲皮素濃度下抑制作用最強，表達水平為對照組的5%（＜0.0001，圖2-D）。的表達水平在0.5%—2.0%的槲皮素濃度下顯著升高，分別為對照組的2.33、2.04和2.31倍，然而在4.0%的槲皮素濃度下被顯著抑制，為對照組的11%（＜0.01，圖2-E）。

將以上5個P450基因在0.5%槲皮素誘導下的表達水平（RT-qPCR結果）與轉錄組測序（RNA-seq）結果進行比較發現，在0.5%槲皮素誘導水平下，5個基因在RT-qPCR與RNA-seq中的表達水平基本一致（圖2-F），驗證了RT-qPCR檢測結果的準確性。

2.4 CYP6ZB2沉默效率

美國白蛾4齡幼蟲的RNA干擾效率檢測結果如圖3所示，結果表明，合成的ds片段可有效抑制靶基因的表達。注射500和1 000 ng ds48 h后，美國白蛾中腸的相對表達量分別比對照組降低了76.3%和96.9%（＜0.001，圖3-A）；注射96 h后，美國白蛾中腸的相對表達量分別比對照組顯著降低了93.9%和95.9%（＜0.0001，圖3-B）。由此可知，注射48 h后，注射劑量1 000 ng比500 ng的沉默效率稍高，但隨著作用時間的延長（至96 h），沉默效率相當。

圖3 注射dsCYP6ZB2 48 h（A）和96 h（B）時CYP6ZB2的相對表達量

2.5 0.5%槲皮素脅迫對CYP6ZB2沉默的美國白蛾幼蟲體重和存活率的影響

圖4為美國白蛾4齡幼蟲注射ds后，用0.5%槲皮素人工飼料飼養10 d的單頭平均體重變化情況。未處理（CK）組和注射ddH2O組取食含0.5%槲皮素人工飼料第6天體重開始增加，而注射ds的美國白蛾幼蟲直到第8天體重才開始增加。沉默的美國白蛾體重增長量第6—10天顯著低于未處理（CK）組（＜0.05）；沉默的美國白蛾在第2、6、8天體重增長量顯著低于注射ddH2O組（＜0.05），表明的表達被干擾后，美國白蛾幼蟲對0.5%槲皮素的解毒適應能力降低，從而導致生長發育受到一定影響，而隨著沉默時間延長，第10天ds注射組體重增長量雖小于ddH2O注射組，但二者沒有顯著性差異，可能是RNAi效果減弱引起的。

取食0.5%槲皮素人工飼料10 d后，注射ds組存活率（53.33%）顯著降低，分別為未處理（CK）組和注射ddH2O組的61.5%和72.7%（＜0.05，圖5），可見，沉默的幼蟲在0.5%槲皮素脅迫下死亡率顯著升高，表明在美國白蛾對槲皮素的適應中發揮了重要作用。

3 討論

3.1 P450在昆蟲解毒代謝植物次生代謝物中起重要作用

P450是一類廣泛分布于生命有機體的代謝酶系，在植食性昆蟲的生長發育、激素生物合成、解毒代謝等過程中發揮著多種生理功能[22-23]。本研究從美國白蛾的中腸轉錄組中篩選出了響應槲皮素誘導而顯著差異表達的5個P450基因，除外，其余均為美國白蛾物種中新鑒定的解毒酶基因。本研究中唯一受槲皮素誘導表達顯著上調的，與同樣屬于亞家族的棉鈴蟲親緣關系較近，不同的是棉鈴蟲的表達會被2-十三烷酮誘導，而經槲皮素誘導后無顯著變化[24]。與本研究結果相似，Hafeez等[25]的研究表明0.2%槲皮素可誘導甜菜夜蛾中腸中的表達，最高誘導水平為對照組的6.38倍，且被干擾的甜菜夜蛾在槲皮素處理后死亡率顯著增加。因而，可以推測美國白蛾的可能同甜菜夜蛾的一樣在槲皮素的解毒代謝中發揮著重要作用。

柱上不同字母表示差異顯著

3.2 美國白蛾P450基因的表達具有顯著的劑量效應

昆蟲P450基因的過表達是對植物次生代謝物和殺蟲劑產生抗性的最常見機制之一[23]。本研究中5個P450基因在不同濃度槲皮素誘導后mRNA表達水平顯著變化，RT-qPCR結果表示槲皮素對美國白蛾的P450基因的表達具有劑量效應，說明這些差異基因可能參與美國白蛾幼蟲對槲皮素的解毒代謝。其中，表現為較低濃度（0.5%—2.0%）被誘導，較高濃度（4.0%）被抑制。其余4個基因在不同槲皮素劑量處理下表達量均下調，表明槲皮素對不同P450基因的誘導作用不同。Li等[26]也發現槲皮素對P450基因的抑制作用比誘導作用更為常見，在其研究的6個家蠶P450基因中，除表達水平保持不變外，其余基因均受到槲皮素的抑制。目前關于槲皮素抑制昆蟲P450基因表達的機制還未見明確的研究報道，由于昆蟲具有復雜的解毒代謝系統，可能這些表達下調的基因并非美國白蛾應對槲皮素解毒代謝的主要基因，而需要其他相關基因的上調表達以應對槲皮素的解毒代謝。

3.3 沉默CYP6ZB2能有效降低美國白蛾幼蟲對槲皮素的適應能力

RNAi是一種通過將雙鏈RNA（dsRNA）引入昆蟲來降低目的基因表達的技術，目前已經成功用于多種昆蟲P450基因功能的研究[27]。dsRNA引入昆蟲體內的方式有口服或注射兩種，其中注射法具有注射效率高，容易控制dsRNA的劑量等優勢，廣泛用于實驗室小規模水平的干擾[28-29]。本研究采用注射的方法顯著降低了美國白蛾幼蟲mRNA的表達水平，且注射劑量越大，沉默效率越高，可見，RNAi可成功用于美國白蛾P450基因功能分析和驗證。

槲皮素能抑制昆蟲生長、減少化蛹和羽化率，嚴重時甚至會導致死亡[30]。通常，桑葉中槲皮素含量會隨著葉片的不同發育階段發生變化，最高濃度可達0.06%[15]，在這一濃度下美國白蛾能正常生長，完成全部的生命周期，本研究中所使用的槲皮素濃度為自然條件下桑葉中槲皮素最高濃度的8.3倍，在這一濃度下注射ds的美國白蛾幼蟲體重增加量和存活率顯著低于未處理組和注射ddH2O組，表達降低使美國白蛾對槲皮素的敏感性增加，導致適應性降低。體重增長量結果還顯示沉默的美國白蛾幼蟲從體重減少恢復到體重增加的時間較對照組時間久，表明沉默會導致幼蟲生長發育遲緩。這與課題組前期利用槲皮素對美國白蛾幼蟲進行生物測定的結果一致，0.5%槲皮素處理使美國白蛾幼蟲發育歷期比對照組延長20.75%（＜0.05），顯著抑制了美國白蛾的生長[13]。這可能是因為較高的槲皮素濃度會使美國白蛾消耗更多的能量用于解毒，相應的也會減少生長發育方面的能量分配，導致幼蟲出現體重增加量減少以及死亡率顯著升高的現象。但由于本研究設置的槲皮素濃度梯度有限，美國白蛾幼蟲對槲皮素耐受的濃度范圍還不明確，在將來的研究中可以增加槲皮素飼養試驗的濃度梯度設置，對參與槲皮素解毒代謝的多基因進行沉默，以便更好地探究美國白蛾對槲皮素代謝機制和途徑。

本試驗中基因沉默的美國白蛾幼蟲在第6—8天體重增加量顯著低于未處理組和注射ddH2O組，但在第10天與注射ddH2O組體重增加量無顯著差異，這可能是RNAi效果持續性減弱引起的，研究表明RNAi對基因的干擾方式不是敲除，而是敲低，干擾只能持續一段時間，隨后基因抑制的效果會減弱甚至消失[31]。Liu等[32]運用RNAi技術注射dsRNA-經檢測在第5天沉默效率最高，為25.3%，隨后減弱，至第8天RNAi對mRNA水平的表達無任何影響。本研究中檢測了48和96 h的基因沉默效率，但是更長時間的干擾效率未知，為了提高RNAi效果，還可以采用多次注射的措施提高RNAi效率，延長其干擾時間[32]。由于上述因素限制，為了多角度探究對美國白蛾幼蟲生長發育以及解毒適應槲皮素的作用，可進一步應用CRISPR/Cas9基因編輯技術探究的功能，以尋求美國白蛾P450基因響應槲皮素的調控元件及其對應的轉錄因子。

4 結論

在槲皮素誘導的美國白蛾中腸轉錄組中篩選出的5個P450基因分別屬于、、和線粒體CYP進化枝。槲皮素對美國白蛾不同P450基因的誘導模式存在差異。的表達被干擾后美國白蛾幼蟲對0.5%槲皮素的適應能力降低，導致幼蟲體重增長緩慢，死亡率顯著升高，表明該基因在美國白蛾對槲皮素的解毒適應過程中發揮重要作用。

[1] MOJARAB-MAHBOUBKAR M, SENDI J J, MAHMOODI N. The sweet wormwood essential oil and its two major constituents are promising for a safe control measure against fall webworm. Pesticide Biochemistry and Physiology, 2022, 184: 105124. doi: 10.1016/j. pestbp.2022.105124.

[2] NING J, LU P F, FAN J T, REN L L, ZHAO L L. American fall webworm in China: A new case of global biological invasions. The Innovation, 2022, 3(1): 100201. doi: 10.1016/j.xinn.2021.100201.

[3] 蔣立娣, 宣貴達, 吳好好, 李麗萍. 桑葉提取物中槲皮素和山萘酚的含量測定. 浙江大學學報(理學版), 2009, 36(6): 705-707, 713. doi: 10.3785/j.issn.1008-9497.2009.06.023.

JIANG L D, XUAN G D, WU H H, LI L P. Determination of quercetin and kaempferol in folium mori extract after hydrolysis by hydrochloric acid. Journal of Zhejiang University (Science Edition)2009, 36(6): 705-707, 713. doi: 10.3785/j.issn.1008-9497.2009.06.023. (in Chinese)

[4] WANG Z Y, NUR F A, MA J Y, WANG J G, CAO C W. Effects of poplar secondary metabolites on performance and detoxification enzyme activity of. Comparative Biochemistry and Physiology, Part C, 2019, 225: 108587. doi: 10.1016/j.cbpc.2019.108587.

[5] 賴城玲, 張珺, 申屠旭萍, 郝培應, 龐琨, 俞曉平. 植物次生代謝物對植食性昆蟲防御作用的研究進展. 應用昆蟲學報, 2022, 59(5): 969-978. doi: 10.7679/j.issn.2095-1353.2022.100.

LAI C L, ZHANG J, SHENTU X P, HAO P Y, PANG K, YU X P. Review of how secondary metabolites defend plants against herbivorous insects. Chinese Journal of Applied Entomology, 2022, 59(5): 969-978. doi: 10.7679/j.issn.2095-1353.2022.100. (in Chinese)

[6] HAFEEZ M, LIU S S, YOUSAF H K, JAN S, WANG R L, FERNáNDEZ-GRANDON G M, LI X W, GULZAR A, ALI B, REHMAN M, ALI S, FAHAD M, LU Y B, WANG M. RNA interference-mediated knockdown of a cytochrome P450 gene enhanced the toxicity of-cypermethrin in xanthotoxin-fed larvae of(Hübner). Pesticide Biochemistry and Physiology, 2020, 162: 6-14. doi: 10.1016/j.pestbp.2019.07.003.

[7] YU H Z, WEN D F, WANG W L, GENG L, ZHANG Y, XU J P. Identification of genes putatively involved in chitin metabolism and insecticide detoxification in the rice leaf folder () larvae through transcriptomic analysis. International Journal of Molecular Sciences, 2015, 16(9): 21873-21896. doi: 10.3390/ijms160921873.

[8] 陳高滿, 陳展博, 葛輝, 楊雪清, 王小奇. 蘋果蠹蛾細胞色素P450基因和的克隆及表達分析. 昆蟲學報, 2020, 63(8): 941-951. doi: 10.16380/j.kcxb.2020.08.004.

CHEN G M, CHEN Z B, GE H, YANG X Q, WANG X Q. Cloning and expression analysis of cytochrome P450 genesand(Lepidoptera: Tortricidae). Acta Entomologica Sinica, 2020, 63(8): 941-951. doi: 10.16380/j.kcxb.2020.08.004. (in Chinese)

[9] FEYEREISEN R. Evolution of insect P450. Biochemical Society transactions, 2006, 34(6): 1252-1255. doi: 10.1042/BST0341252.

[10] YU L Y, TANG W Q, HE W Y, MA X L, VASSEUR L, BAXTER S W, YANG G, HUANG S G, SONG F Q, YOU M S. Characterization and expression of the cytochrome P450 gene family in diamondback moth,(L.). Scientific Reports, 2015, 5: 8952. doi: 10.1038/srep08952.

[11] AMEZIAN D, NAUEN R, GOFF G L. Comparative analysis of the detoxification gene inventory of four majorpest species in response to xenobiotics. Insect Biochemistry and Molecular Biology, 2021, 138: 103646. doi: 10.1016/j.ibmb.2021.103646.

[12] WANG R L, LIU S W, BAERSON S R, QIN Z, MA Z H, SU Y J, ZHANG J E. Identification and functional analysis of a novel cytochrome P450 geneassociated with pyrethroid detoxification inHübner. International Journal of Molecular Sciences, 2018, 19(3): 737. doi: 10.3390/ijms19030737.

[13] 潘忠玉. 3種次生代謝物對美國白蛾幼蟲生長發育及解毒酶活性的影響[D]. 北京: 北京林業大學, 2020.

PAN Z Y. Effects of three secondary metabolites on the growth and development and detoxification enzyme activities in(Lepidoptera: Arctiidae)[D]. Beijing: Beijing Forestry University, 2020. (in Chinese)

[14] 王海燕, 楊金龍, 谷山林, 王小燕, 王介平, 周嬋, 呂金鳳, 曾秀. 桑葉槲皮素提取物抗氧化活性研究. 絲綢, 2018, 55(3): 15-20. doi: 10.3969/j.issn.1001-7003.2018.03.003.

WANG H Y, YANG J L, GU S L, WANG X Y, WANG J P, ZHOU C, Lü J F, ZENG X. Study on antioxidant activity of quercetin extract from mulberry leaves. Journal of Silk, 2018, 55(3): 15-20. doi: 10.3969/j.issn.1001-7003.2018.03.003. (in Chinese)

[15] 孫蓮, 孟磊, 陳堅, 馬季, 胡瑞, 賈殿增. 毛細管電泳法測定桑葉中的黃酮類成分——蘆丁和槲皮素. 色譜, 2001, 19(5): 395-397.

SUN L, MENG L, CHEN J, MA J, HU R, JIA D Z. Determination of rutin and quercetin in mulberry leaves by high performance capillary electrophoresis. Chinese Journal of Chromatography, 2001, 19(5): 395-397. (in Chinese)

[16] 曹利軍, 楊帆, 唐思瑩, 陳敏. 適合三種鱗翅目昆蟲的一種人工飼料配方. 應用昆蟲學報, 2014, 51(5): 1376-1386. doi: 10.7679/j.issn. 2095-1353.2014.165.

CAO L J, YANG F, TANG S Y, CHEN M. Development of an artificial diet for three lepidopteran insects. Chinese Journal of Applied Entomology, 2014, 51(5): 1376-1386. doi: 10.7679/j.issn. 2095-1353.2014.165. (in Chinese)

[17] ZHANG D, GAO F L, JAKOVLI? I, ZOU H, ZHANG J, LI W X, WANG G T. PhyloSuite: An integrated and scalable desktop platform for streamlined molecular sequence data management and evolutionary phylogenetics studies. Molecular ecology resources, 2020, 20(1): 348-355. doi: 10.1111/1755-0998.13096.

[18] GUINDON S, DUFAYARD J F, LEFORT V, ANISIMOVA M, HORDIJK W, GASCUEL O. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Systematic Biology, 2010, 59(3): 307-321. doi: 10.1093/ sysbio/syq010.

[19] MINH B Q, NGUYEN M A, von HAESELER A. Ultrafast approximation for phylogenetic bootstrap. Molecular Biology and Evolution, 2013, 30(5): 1188-1195. doi: 10.1093/molbev/mst024.

[20] 陶蓉, 李慧, 孫宇航, 于曉航, 朱晗, 郝德君. 美國白蛾內參基因的鑒定及篩選. 林業科學, 2019, 55(9): 111-120. doi: 10.11707/j. 1001-7488.20190912.

TAO R, LI H, SUN Y H, YU X H, ZHU H, HAO D J. Identification and screening of internal reference genes of(Lepidoptera: Arctiidae). Scientia Silvae Sinicae, 2019, 55(9): 111-120. doi: 10.11707/j.1001-7488.20190912. (in Chinese)

[21] LIVAK K J, SCHMITTGEN T D. Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCTmethod. Methods, 2001, 25(4): 402-408. doi: 10.1006/meth.2001.1262.

[22] 尹飛, 李振宇, SAMINA S, 林慶勝. P450基因在氯蟲苯甲酰胺不同抗性品系小菜蛾中的表達及功能分析. 中國農業科學, 2022, 55(13): 2562-2571. doi: 10.3864/j.issn.0578-1752.2022.13.007.

YIN F, LI Z Y, SAMINA S, LIN Q S. Expression and function analysis of cytochrome P450 genes inwith differentchlorantraniliprole resistance. Scientia Agricultura Sinica, 2022, 55(13): 2562-2571. doi: 10.3864/j.issn.0578-1752.2022.13.007. (in Chinese)

[23] WANG R L, XIA Q Q, BAERSON S R, REN Y, WANG J, SU Y J, ZHENG S C, ZENG R S. A novel cytochrome P450gene in(Lepidoptera: Noctuidae) and its potential role in plant allelochemical detoxification. Journal of Insect Physiology, 2015, 75: 54-62. doi: 10.1016/j.jinsphys.2015.02.013.

[24] LIU X N, LIANG P, GAO X W, SHI X Y. Induction of the cytochrome P450 activity by plant allelochemicals in the cotton bollworm,(Hübner). Pesticide Biochemistry and Physiology, 2006, 84(2): 127-134. doi:10.1016/j.pestbp.2005.06.002.

[25] HAFEEZ M, QASIM M, ALI S, YOUSAF H K, WAQAS M, ALI E, AHMAD M A, JAN S, BASHIR M A, NOMAN A, WANG M, GHARMH H A, KHAN K A. Expression and functional analysis of P450 gene induced tolerance/resistance to lambda-cyhalothrin in quercetin fed larvae of beet armyworm(Hübner). Saudi Journal of Biological Sciences, 2020, 27(1): 77-87. doi: 10.1016/j.sjbs.2019.05.005.

[26] LI B, ZHANG H, NI M, WANG B B, LI F C, XU K Z, SHEN W D, XIA Q Y, ZHAO P. Identification and characterization of six cytochrome P450 genes belonging toandgene families in the silkworm,. Molecular Biology Reports, 2014, 41(8): 5135-5146. doi: 10.1007/s11033-014-3379-z.

[27] KATSAVOU E, RIGA M, IOANNIDIS P, KING R, ZIMMER C T, VONTAS J. Functionally characterized arthropod pest and pollinator cytochrome P450s associated with xenobiotic metabolism. Pesticide Biochemistry and Physiology, 2021, 181: 105005. doi: 10.1016/j. pestbp.2021.105005.

[28] JIN M, LIAO C, FU X, HOLDBROOK R, WU K, XIAO Y. Adaptive regulation of detoxification enzymes into different host plants. Insect Molecular Biology, 2019, 28(5): 628-636. doi: 10.1111/imb.12578.

[29] LI H, JIANG W H, ZHANG Z, XING Y R, LI F. Transcriptome analysis and screening for potential target genes for RNAi-mediated pest control of the beet armyworm,. PloS one, 2013, 8(6): e65931. doi: 10.1371/journal.pone.0065931.

[30] CHEN C Y, HAN P, YAN W Y, WANG S Y, SHI X Y, ZHOU X G, DESNEUX N, GAO X W. Uptake of quercetin reduces larval sensitivity to lambda-cyhalothrin in. Journal of Pest Science, 2018, 91(2): 919-926. doi: 10.1007/s10340-017-0933-1.

[31] HUVENNE H, SMAGGHE G. Mechanisms of dsRNA uptake in insects and potential of RNAi for pest control: A review. Journal of Insect Physiology, 2010, 56(3): 227-235. doi: 10.1016/j.jinsphys. 2009.10.004.

[32] LIU S H, DING Z P, ZHANG C W, YANG B J, LIU Z W. Gene knockdown by intro-thoracic injection of double-stranded RNA in the brown planthopper,. Insect Biochemistry and Molecular Biology, 2010, 40(9): 666-671. doi: 10.1016/j.ibmb.2010.06.007.

Analysis of Cytochrome P450 Genes in Response to Quercetin and Function of

SHAO HongYang, MENG Xiang, ZHANG Tao, CHEN Min

College of Forestry, Beijing Forestry University/Beijing Key Laboratory for Forest Pest Control, Beijing 100083

【】The objective of this study is to clarify the differentially expressed cytochrome P450 genes and phylogenetic relationships in response to quercetin in, and to explore the dosage effect of cytochrome P450 genes in response to quercetin, and the role of cytochrome P450 genes in detoxification adaptation to quercetin in.【】Based on the transcriptome data of quercetin-induced midgut of, the significantly differentially expressed cytochrome P450 genes were screened. MEGA X software was used to analyze the phylogenetic relationships of the target genes with cytochrome P450 genes of other related species. RT-qPCR was used to detect the expression of the five cytochrome P450 genes in different quercetin concentrations (0.5%, 1.0%, 2.0% and 4.0%). The full-length cDNA sequence ofwas cloned, and the effect ofgene-silenced on the body weight and survival rate ofwas investigated by RNA interference (RNAi) technique.【】 Five cytochrome P450 genes were identified from the differentially expressed genes in the transcriptome, and were named,,,andby Standardized Cytochrome P450 Nomenclature Committee (GenBank accession numbers: ON032363-ON032367). Phylogenetic tree analysis showed that these five genes were divided into four clades (,,and the mitochondrial CYP clade). RT-qPCR results showed that cytochrome P450 genes had a significant dosage effect in response to quercetin, and the expressions of,,andwere inhibited by different concentrations of quercetin. The expression ofwas significantly increased at 0.5%-2.0% concentration of quercetin, but it was significantly inhibited at 4.0% quercetin. The full-length cDNA sequence ofwas cloned.had a complete open reading frame (ORF) encoding 500 amino acids. RNAi results showed that the expression level ofwas significantly inhibited by injection of two doses(500 and 1 000 ng) ds, andthe maximum reduction inexpression level was 96.9% after injection of 1 000 ng dsfor 48 h. The weight increase ofgene-silencedlarvae fed on 0.5% quercetin was lower than that of the control group. Furthermore, the survival rate (53.33%) was 61.5% of that of the untreated group and 72.7% of that of the ddH2O group (＜0.05).【】Five cytochrome P450 genes were identified in the transcriptome data of quercetin-induced midgut of, and quercetin had a significant dosage effect on the five genes, indicating that the response patterns of different cytochrome P450 genes to quercetin were different. RNAi can significantly inhibit the expression level of, and quercetin has a significant negative impact on the growth of, suggesting thatplays an important role in quercetin detoxification in.

; quercetin; cytochrome P450; expression analysis; RNA interference (RNAi); gene function

2022-12-20；

2023-01-31

國家重點研發計劃（2021YFD1400300）

邵紅揚，E-mail：hongyangshao2021@163.com。孟香，E-mail：1210987080@qq.com。邵紅揚和孟香為同等貢獻作者。通信作者陳敏，E-mail：minch@bjfu.edu.cn

（責任編輯 岳梅）