響應(yīng)槲皮素誘導的美國白蛾P(guān)450基因及CYP6ZB2功能分析

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

響應(yīng)槲皮素誘導的美國白蛾P(guān)450基因及功能分析

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

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

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

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

0 引言

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

1 材料與方法

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

1.1 槲皮素人工飼料的配制

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

1.2 供試昆蟲

美國白蛾幼蟲為實驗室人工飼料飼養(yǎng)的種群，人工飼料的配制參照曹利軍等[16]。飼養(yǎng)條件設(shè)定為溫度（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（諾唯贊，中國南京）。

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

1.4 響應(yīng)槲皮素誘導的P450基因鑒定與系統(tǒng)發(fā)育分析

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

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

1.5 P450基因響應(yīng)槲皮素的表達分析

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

1.6 CYP6ZB2的功能驗證

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

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

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

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

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

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

1.7 數(shù)據(jù)分析

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

2 結(jié)果

2.1 響應(yīng)槲皮素誘導的P450基因的篩選與鑒定

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

表2 響應(yīng)槲皮素誘導的P450基因的BLASTX結(jié)果

Down：表達下調(diào)down-regulated expression；Up：表達上調(diào)up-regulated expression

2.2 P450基因的系統(tǒng)發(fā)育分析

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

Bmor：家蠶B. mori；Harm：棉鈴蟲 H. armigera；Sfru：草地貪夜蛾 S. frugiperda；Slit：斜紋夜蛾S. litura；Tni：粉紋夜蛾 T. ni。不同顏色代表不同P450基因進化枝；紅色字體基因代表本研究鑒定的5個P450基因；分支上的數(shù)字表示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基因的劑量效應(yīng)

RT-qPCR結(jié)果顯示，、和的表達均被不同濃度槲皮素抑制，而較低濃度（0.5%—2.0%）的槲皮素誘導的表達顯著上調(diào)。、和的表達隨著槲皮素濃度升高，抑制作用逐漸減弱，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結(jié)果）與轉(zhuǎn)錄組測序（RNA-seq）結(jié)果進行比較發(fā)現(xiàn)，在0.5%槲皮素誘導水平下，5個基因在RT-qPCR與RNA-seq中的表達水平基本一致（圖2-F），驗證了RT-qPCR檢測結(jié)果的準確性。

2.4 CYP6ZB2沉默效率

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

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

3 討論

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

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

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

3.2 美國白蛾P(guān)450基因的表達具有顯著的劑量效應(yīng)

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

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

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

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

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

4 結(jié)論

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

[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] 賴城玲, 張珺, 申屠旭萍, 郝培應(yīng), 龐琨, 俞曉平. 植物次生代謝物對植食性昆蟲防御作用的研究進展. 應(yīng)用昆蟲學報, 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種次生代謝物對美國白蛾幼蟲生長發(fā)育及解毒酶活性的影響[D]. 北京: 北京林業(yè)大學, 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] 曹利軍, 楊帆, 唐思瑩, 陳敏. 適合三種鱗翅目昆蟲的一種人工飼料配方. 應(yīng)用昆蟲學報, 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] 陶蓉, 李慧, 孫宇航, 于曉航, 朱晗, 郝德君. 美國白蛾內(nèi)參基因的鑒定及篩選. 林業(yè)科學, 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基因在氯蟲苯甲酰胺不同抗性品系小菜蛾中的表達及功能分析. 中國農(nóng)業(yè)科學, 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

國家重點研發(fā)計劃（2021YFD1400300）

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

（責任編輯 岳梅）