盧柏亦　蔣哲　鄢柳慧　鐘小惠　李晶瑩　黃福鋼　李發(fā)活　莫德原　邱永福

摘要：【目的】篩選出同時響應(yīng)褐飛虱和稻癭蚊取食的基因，揭示水稻幼苗應(yīng)對2種害蟲時基因表達層面的差異，為后續(xù)挖掘廣譜抗蟲基因和解析水稻抗蟲機制打下理論基礎(chǔ)?！痉椒ā恳运酒贩N9311為供試植物，分別對15日齡的水稻幼苗進行褐飛虱和稻癭蚊接蟲處理，觀測幼苗表型，并通過轉(zhuǎn)錄組測序技術(shù)分別對褐飛虱接入后24 h和稻癭蚊接入后7 d的水稻幼苗進行轉(zhuǎn)錄組測序。以水稻日本晴基因組作為參考基因組進行對比，利用FPKM法計算基因表達量，設(shè)定參數(shù)（|log₂FC|>1且P<0.05）篩選差異表達基因。結(jié)合基因差異表達分析和功能富集分析，研究水稻響應(yīng)2種害蟲的機制異同點?！窘Y(jié)果】鑒定出響應(yīng)褐飛虱取食的差異表達基因3963個，響應(yīng)稻癭蚊取食的差異表達基因1206個，有251個差異表達基因同時響應(yīng)2種害蟲，其中108個具有相同表達模式，143個具有相反表達模式。GO功能注釋分析表明，褐飛虱取食顯著影響植物體內(nèi)細胞壁合成和缺水響應(yīng)，而稻癭蚊取食則對植物光合作用的影響更顯著，具有相同表達模式的差異表達基因主要富集在光合作用、水楊酸合成、茉莉酸信號轉(zhuǎn)導(dǎo)和傷害響應(yīng)等生物學(xué)功能，而具有相反表達模式的差異表達基因僅富集在乙醛酸循環(huán)。KEGG信號通路富集分析表明，褐飛虱和稻癭蚊取食均顯著影響植物體內(nèi)次生代謝物的合成，相同表達模式的差異表達基因富集到MAPK信號通路、植物激素信號轉(zhuǎn)導(dǎo)通路及光合作用等6個通路，而相反表達模式的差異表達基因沒有富集的通路?！窘Y(jié)論】水稻應(yīng)對褐飛虱和稻癭蚊的基因表達調(diào)控存在相同點和不同之處，共響應(yīng)2種害蟲的調(diào)控通路可能在水稻抗蟲過程中發(fā)揮重要作用。

關(guān)鍵詞： 水稻;褐飛虱;稻癭蚊;轉(zhuǎn)錄組

中圖分類號： S435.112? ? ? ? ? ? ? ? ? ? ? ? 文獻標志碼： A 文章編號：2095-1191（2022）03-0628-13

Transcriptome analysis of rice seedlings responding to brown planthopper infestation and to gall midge infestation

LU Bai-yi， JIANG Zhe， YAN Liu-hui， ZHONG Xiao-hui， LI Jing-ying，

HUANG Fu-gang， LI Fa-huo， MO De-yuan， QIU Yong-fu

（College of Agriculture， Guangxi University/State Key Laboratory for Conservation and Utilization of Subtropical

Agro-bioresources， Nanning， Guangxi? 530004， China）

Abstract：【Objective】To screen out the functional genes that respond to both brown planthopper （BPH） infestation and rice gall midge （GM） infestation and reveal the difference of gene expression level in rice seedlings responding to the two pests， so as to provide theoretical basis for new broad-spectrum resistance and rice anti-insect mechanism. 【Method】Rice 9311 was the tested rice variety. Fifteen-day-old seedlings were treated with BPH and GM respectively to observe seedling phenotypes and to sample for transcriptome which sequence was performed on rice seedlings that were collected 24 hours after BPH infestation and 7 days after GM infestation， respectively. Taking the genome of Nipponbare （MSU7） as a reference to align the sequences， the differentially expressed genes （DEGs） were screened by calculating the gene expression level in Fragments Per Kilobase per Million （FPKM）， with （|log₂FC|≥1 and P≤0.05） as the parameter. Through gene differential expression analysis and functional enrichment analysis， the similarities and differences of rice response mechanisms to the two pests were identified. 【Result】The results showed that 3963 genes were differentially expressed in BPH infestation and 1206 genes were differentially expressed in GM infestation. And 251 genes were found responsive to both pests， of which 108 genes with the same expression pattern and 143 genes with opposite expression patterns. GO enrichment analysis showed that：BPH infestation significantly influenced the cell wall biogenesis and the response to water deprivation in rice seedlings; seedlings of GM infestation enriched GO terms were mostly chloroplast or photosynthesis-related; the DEGs with the same expression patterns of two pest stresses were enriched to some biological processes such as photosynthesis， salicylic acid biosynthetic process， regulation of jasmonic acid mediated signaling pathway and response to wounding; whereas the DEGs with the opposite expression patterns of two pest stresses were enriched to glyoxy-late cycle only. KEGG enrichment analysis indicated that：Both BPH infestation and GM infestationsignificantly affected the biosynthesis of secondary metabolites in seedlings; the DEGs with the same expression patterns in both treatments were enriched into six pathways mainly involved in MAPK signaling pathway， plant hormone signal transduction pathway， photosynthesis pathway; while no significant KEGG pathway was enriched in the DEGs with opposite expression patterns. 【Conclusion】These results suggest that the common pathways in two pest infestations probably play an important role in rice resistance while similarities and differences exist in the regulations of gene expression between BPH infestation and GM infestation.C3EF4567-7C12-48BB-84A9-C29F2CF52217

Key words： rice; brown planthopper; gall midge; transcriptome

Foundation items： Guangxi Natural Science Foundation （2020GXNSFDA297008）

0 引言

【研究意義】水稻是我國主要糧食作物，在第一次綠色革命后，水稻由高稈向矮稈轉(zhuǎn)變，對水稻危害嚴重的害蟲種類也發(fā)生了明顯變化，且害蟲種類不斷增多，蟲災(zāi)日益頻繁（嚴火其，2021）。褐飛虱（Brown planthopper，BPH）和稻癭蚊（Gall midge，GM）是我國南方災(zāi)害性害蟲，因危害方式（褐飛虱通過刺吸，稻癭蚊通過咀嚼）和危害部位的差異（Li et al.，2020;Shi et al.，2021），水稻對這2種害蟲的防御機制也有所不同，只有極少數(shù)水稻品種可兼顧對2種害蟲的抗性。因此，做好有效的防治，實現(xiàn)對褐飛虱和稻癭蚊災(zāi)害發(fā)生風(fēng)險的有力控制具有重要意義，而防治2種害蟲最科學(xué)、安全的手段是挖掘和利用內(nèi)源性廣譜抗蟲基因（楊利艷等，2020;嚴火其，2021），并以此培育抗蟲品種，從而減少災(zāi)害損失和農(nóng)藥的施用?！厩叭搜芯窟M展】在褐飛虱和稻癭蚊雙抗材料研究方面，已報道并進行遺傳分析的雙抗水稻品種有：江潮、RP^{1976-18-6-4-2}和ARC5833（黃鳳寬等，2006;吳碧球等，2014;秦學(xué)毅等，2016;鄢柳慧等，2020）。RP^{1976-18-6-4-2}對褐飛虱和稻癭蚊的抗性評價為抗，由1對顯性基因和1對隱性基因控制對褐飛虱的抗性，由1對顯性基因控制對稻癭蚊的抗性（黃鳳寬等，2006）。江潮對2種害蟲均表現(xiàn)為高抗，是由1對顯性基因控制對稻癭蚊的抗性、1對隱性基因控制對褐飛虱的抗性，2個基因間獨立遺傳（秦學(xué)毅等，2016）。ARC5833與江潮的抗性表現(xiàn)一致，對2種害蟲表現(xiàn)為高抗，但在抗性基因的遺傳上比江潮復(fù)雜。ARC5833對褐飛虱的抗性由單一顯性基因控制，對稻癭蚊中國II型由1對顯性基因和1對隱性基因控制，而對稻癭蚊中國IV型則由1對顯性基因控制（吳碧球等，2014;鄢柳慧等，2020）。已發(fā)現(xiàn)的3個雙抗褐飛虱和稻癭蚊水稻品種對2種害蟲的控制基因為分開獨立遺傳，增加了后期基因鑒定、抗蟲機制解析及利用的復(fù)雜程度。在水稻抗蟲機制研究方面，組學(xué)研究（主要包括轉(zhuǎn)錄組、蛋白組和代謝組等）作為了解植物對昆蟲取食反應(yīng)的重要技術(shù)途徑之一，已有大量的組學(xué)研究揭示水稻幼苗響應(yīng)褐飛虱取食的機制，卻很少研究水稻對稻癭蚊的響應(yīng)機制。代謝組學(xué)相關(guān)的抗蟲研究表明，GABA支路和莽草酸介導(dǎo)的次生代謝途徑對水稻抗褐飛虱具有重要作用（Liu et al.，2010），而脂肪酸及其衍生物、氨基酸及次生代謝物則在水稻抗稻癭蚊中發(fā)揮重要作用（陳燦等，2021）。利用抗蟲激素組學(xué)對Bph6和Bph9的9311背景的近等基因系（NIL-Bph6和NIL-Bph9，以9311為背景材料）以及對照9311為材料進行檢測，發(fā)現(xiàn)在褐飛虱取食過程中2個近等基因系的SA、JA和JA-Ile上調(diào)顯著且響應(yīng)速度快，而感性品種9311則激素水平變化小且反應(yīng)時間滯后（Zhao et al.，2016;杜波等，2018;Guo et al.，2018）。通過轉(zhuǎn)錄組技術(shù)測定褐飛虱取食48 h后的水稻幼苗基因表達變化，揭示了褐飛虱取食能激發(fā)水稻的應(yīng)激反應(yīng)和保護機制，同時降低營養(yǎng)合成的過程，解釋了褐飛虱為害造成水稻減產(chǎn)的原因（曹征鴻等，2020）。在水稻抗蟲基因挖掘方面，已有研究利用抑制消減雜交技術(shù)（SSH）通過篩選差異表達基因，鑒定出OsERF3、OsWRKY70和OsHI-LOX等對水稻褐飛虱和二化螟抗性均具有顯著影響的基因，這3個基因均能增強水稻對二化螟的抗性，而降低水稻對褐飛虱的抗性，在植株內(nèi)表現(xiàn)出相反的調(diào)節(jié)模式（Zhou et al.，2009;Lu et al.，2011;Li et al.，2015）。OsERF3、OsWRKY70和OsHI-LOX等基因的研究結(jié)果說明，通過組學(xué)研究可篩選出具有調(diào)控水稻對刺吸式昆蟲和咀嚼式昆蟲抗性的基因，實現(xiàn)單一基因調(diào)控水稻對不同種類昆蟲的抗性具有可行性。【本研究切入點】目前，關(guān)于單一基因控制水稻共調(diào)控褐飛虱和稻癭蚊的抗性機理研究尚無文獻報道，系統(tǒng)揭示水稻響應(yīng)刺吸式昆蟲和咀嚼式昆蟲取食間的組學(xué)差異研究也鮮見報道?！緮M解決的關(guān)鍵問題】通過研究水稻在褐飛虱和稻癭蚊脅迫下的轉(zhuǎn)錄水平變化，比較2種脅迫下差異表達基因的功能異同，獲得水稻對2種害蟲取食共響應(yīng)的基因，明確篩選出雙抗基因的可行性，為后續(xù)挖掘廣譜抗蟲基因和解析水稻抗蟲機制打下理論基礎(chǔ)。

1 材料與方法

1. 1 試驗材料

1. 1. 1 供試植物 供試水稻品種為9311，于溫室培育。溫室條件：溫度（26±2）℃，光周期 L∶D=16 h∶8 h，相對濕度（75±5）%。

1. 1. 2 供試害蟲 供試褐飛虱采自廣西南寧市廣西大學(xué)，供試稻癭蚊收集自廣西崇左市大新縣，2種害蟲均在溫室的9311植株上繁殖多代。

1. 2 樣品處理

在褐飛虱取食處理中，每杯（口徑10 cm，深度16 cm）種植10株水稻幼苗，15日苗齡接入褐飛虱2～3齡若蟲，每株苗8頭蟲，即每杯接入80頭蟲。在稻癭蚊取食處理中，每杯種植8株稻苗，15日苗齡接入稻癭蚊成蟲，每杯苗接入4頭蟲。

1. 3 表型獲取

加入褐飛虱若蟲后第3 d，用佳能70D相機拍攝褐飛虱為害后的水稻植株表型;加入稻癭蚊成蟲后第16 d，拍攝稻癭蚊為害后的植株表型。用小米11微距相機拍攝取食中的褐飛虱表型，用體式顯微鏡拍攝取食中的稻癭蚊表型。

1. 4 測序樣品采集

接入褐飛虱2～3齡若蟲后昆蟲立即取食，24 h可觀測到苗尖輕微萎蔫癥狀，72 h可觀測到稻苗完全干枯，因此選取24 h作為褐飛虱處理的取樣時間點。接入稻癭蚊成蟲后前3 d為產(chǎn)卵孵化期，7 d為蟲體快速增長期，因此選擇7 d作為稻癭蚊處理的取樣時間點。褐飛虱取食24 h后，取出水稻，剪取莖稈組織即褐飛虱取食部位作為轉(zhuǎn)錄組測序樣品存于液氮中。接入稻癭蚊后第7 d，取出水稻，剪取靠近根部的莖稈組織即稻癭蚊取食部位存于液氮中。褐飛虱和稻癭蚊處理各設(shè)1個對照組，對照組不做任何處理，各處理設(shè)3個生物學(xué)重復(fù)，總計12個。C3EF4567-7C12-48BB-84A9-C29F2CF52217

1. 5 RNA處理與轉(zhuǎn)錄組測序

采用TRIzol法提取水稻樣品總RNA，經(jīng)瓊脂糖凝膠電泳和微量分光光度計（NanoDrop 2000）檢測樣品合格后用于cDNA文庫構(gòu)建。利用PicoGreen和熒光分光光度計方法定量文庫[雙鏈DNA（dsDNA）定量檢測試劑盒Quant-iT PicoGreen dsDNA Assay Kit購自Invitrogen公司;熒光分光光度計Quantifluor-ST fluorometer購自Promega公司]。使用Agilent 2100生物分析儀（購自Agilent公司，并搭配使用與其配套的DNA試劑盒Agilent High Sensitivity DNA Kit運行DNA樣品分析）對PCR富集片段進行質(zhì)量控制，驗證DNA文庫的片段大小及分布。將合格樣品分別稀釋至10 nmol/L后，逐步稀釋定量至4～5 nmol/L進行上機測序，采用Illumina HiSeqTM 2000平臺完成轉(zhuǎn)錄組測序。建庫和測序均由上海派森諾生物科技有限公司完成。

1. 6 差異基因表達分析

測序原始數(shù)據(jù)序列經(jīng)質(zhì)量分析去除低質(zhì)量序列和接頭序列后，得到可供進行后續(xù)分析的序列。因9311基因組序列組裝質(zhì)量不如日本晴，因此選用水稻基因組注釋數(shù)據(jù)庫（http：//rice.uga.edu/，日本晴基因組）作為參照基因組，使用HTSeq統(tǒng)計比對到每個基因上的Read count值，作為基因原始表達量，然后采用FPKM對表達量進行標準化。采用DESeq對基因表達量進行差異分析，差異表達基因（Differentially expressed genes，DEGs）篩選條件為：|log2 FC|>1且P<0.05。

1. 7 GO功能注釋和KEGG信號通路富集分析

使用基迪奧云平臺對差異基因進行數(shù)據(jù)分析（https：//www.omicshare.com/）。差異表達基因的GO功能注釋分析和KEGG信號通路富集分析分別通過云平臺的動態(tài)GO功能注釋分析云工具和動態(tài)KEGG信號通路富集分析云工具完成。

1. 8 轉(zhuǎn)錄組數(shù)據(jù)可視化

使用基迪奧云平臺對差異基因進行數(shù)據(jù)可視化（https：//www.omicshare.com/）：差異表達基因數(shù)量的可視化由柱狀圖和維恩圖云工具完成;GO功能注釋分析和KEGG信號通路富集分析可視化由動態(tài)富集差異氣泡圖和動態(tài)富集條形圖云工具完成;差異表達基因表達量可視化由動態(tài)熱圖云工具完成。

2 結(jié)果與分析

2. 1 褐飛虱和稻癭蚊取食對9311幼苗的影響

2種害蟲取食對9311水稻品種幼苗產(chǎn)生的影響不同。接入褐飛虱3 d后，幼苗表現(xiàn)枯死癥狀（圖1-A）;圖1-B和圖1-C展示在水稻莖稈外部吸食汁液的褐飛虱若蟲。稻癭蚊對稻苗的為害特點為生長點基部膨大，15 d左右可明顯觀察到稻癭蚊取食后水稻頂葉轉(zhuǎn)變?yōu)槭[管（圖1-D）;圖1-E和圖1-F分別展示接入稻癭蚊成蟲16和7 d后正在水稻蔥管內(nèi)取食的稻癭蚊若蟲。

2. 2 差異表達基因分析結(jié)果

轉(zhuǎn)錄組測序共設(shè)12個樣本，包括褐飛虱處理組、褐飛虱處理的對照組、稻癭蚊處理組和稻癭蚊處理的對照組，各3個樣本。從整體上看，每個樣本最后得到4000萬～4700萬條Reads，Q30為93%～96%。過濾后得到的高質(zhì)量數(shù)據(jù)比對參考基因組，樣本中有83%～87%的測序序列可定位到參考序列上。表明整體測序數(shù)據(jù)質(zhì)量良好，數(shù)據(jù)可用于進一步的轉(zhuǎn)錄組分析。

經(jīng)過褐飛虱處理，共有3963個基因表達發(fā)生明顯變化，其中有1518個基因上調(diào)，2445個基因下調(diào)。而稻癭蚊處理下僅發(fā)現(xiàn)1206個差異表達基因，其中461個基因上調(diào)，745個基因下調(diào)（圖2-A）。共有251個差異表達基因在2種害蟲取食后表達量均發(fā)生明顯變化，在251個基因中有55個基因在2種處理中表達量均上調(diào);53個基因在2種處理中均下調(diào);91個基因在褐飛虱處理中表達上調(diào)、在稻癭蚊處理中表達下調(diào);52個基因在褐飛虱處理中表達下調(diào)、在稻癭蚊處理中表達上調(diào)（圖2-B）。即本研究中對褐飛虱和稻癭蚊具有相同的響應(yīng)模式基因108個，具有相反響應(yīng)模式基因143個。

2. 3 差異表達基因GO功能注釋分析結(jié)果

將差異表達基因進行GO功能注釋分析，發(fā)現(xiàn)響應(yīng)褐飛虱脅迫的基因和響應(yīng)稻癭蚊脅迫的基因在GO功能條目上具有較大差異。響應(yīng)褐飛虱脅迫的基因參與的生物學(xué)過程主要注釋在細胞壁合成、缺水響應(yīng)、碳水化合物代謝、脂肪酸代謝、酸性化學(xué)物質(zhì)響應(yīng)和激素響應(yīng)，基因的表達產(chǎn)物在細胞中的定位注釋在細胞外周和細胞膜，基因的分子功能注釋于水解酶活性和催化活性（圖3-A）。響應(yīng)稻癭蚊脅迫的差異表達基因在光合作用和氧化還原等生物過程富集，在類囊體相關(guān)、光合膜、光合系統(tǒng)、葉綠體、質(zhì)體和DNA聚合酶等細胞組分上富集，在氧化還原活性相關(guān)、四吡咯結(jié)合和輔因子結(jié)合等分子功能富集（圖3-B）。

將2種害蟲的251個共響應(yīng)差異表達基因分為2類進行GO功能注釋分析，一類是在共響應(yīng)中具有相同響應(yīng)模式的108個差異表達基因，另一類是在共響應(yīng)中具有相反響應(yīng)模式的143個差異表達基因。具有相同響應(yīng)模式的差異表達基因在功能條目上主要包括：光合作用、水楊酸合成與代謝、茉莉酸介導(dǎo)的信號通路、細胞過程調(diào)節(jié)和傷害響應(yīng)等生物過程;光合系統(tǒng)、光合膜、類囊體、葉綠體和質(zhì)體等細胞組分（圖4-A）。這些功能條目在整體上與響應(yīng)稻癭蚊脅迫的基因較相似，但在富集程度上也存在一定差異，如響應(yīng)植物激素的基因在相同響應(yīng)模式的基因中比在響應(yīng)稻癭蚊脅迫的基因中更富集（圖4-A）。而在2種害蟲脅迫過程中具有相反響應(yīng)模式的差異表達基因最主要集中在乙醛酸循環(huán)相關(guān)的生物過程，只有乙醛酸循環(huán)、乙醛酸代謝過程和乙醛酸循環(huán)體3個生物過程的GO富集P<0.05（圖4-B）。說明較多的生物學(xué)功能在水稻應(yīng)對2種害蟲時響應(yīng)模式相似，而響應(yīng)模式相反的生物學(xué)功能較少。C3EF4567-7C12-48BB-84A9-C29F2CF52217

2. 4 差異表達基因的KEGG信號通路富集分析結(jié)果

對差異表達基因前20個KEGG信號通路分析發(fā)現(xiàn)，在褐飛虱和稻癭蚊脅迫富集基因最多的前3個通路均為次生代謝物的生物合成、代謝途徑和苯丙素的生物合成。除該3條代謝通路外，在褐飛虱脅迫下基因在以下通路中富集：代謝通路中的脂肪酸合成與代謝、黃酮與類黃酮的生物合成、糖分代謝、酸類代謝及角質(zhì)、木栓質(zhì)和蠟質(zhì)的生物合成等;環(huán)境信息處理通路中的植物激素信號傳導(dǎo);有機系統(tǒng)通路上的植物與病原菌互作（圖5-A）。與褐飛虱的響應(yīng)機制不同，稻癭蚊處理后基因主要參與以下途徑：代謝通路中的光合作用、玉米素合成、氨基酸代謝、亞油酸代謝、各種次生代謝物的生物合成;遺傳信息處理通路中的DNA復(fù)制（圖5-B）。此外，卟啉和葉綠素代謝通路及甘油磷脂代謝通路在2種害蟲取食中的KEGG信號通路富集均為前20，但這2個通路的響應(yīng)程度在褐飛虱與稻癭蚊間存在明顯差異，卟啉和葉綠素代謝通路對稻癭蚊的響應(yīng)程度高于褐飛虱，而甘油磷脂通路對褐飛虱的響應(yīng)程度高于稻癭蚊，且甘油磷脂通路在稻癭蚊響應(yīng)中不屬于顯著富集，富集程度位于閾值P=0.05之下（圖5-A和圖5-B）。

對共響應(yīng)基因進行KEGG信號通路富集分析，發(fā)現(xiàn)在相同響應(yīng)模式的基因中僅有6個通路表現(xiàn)為富集：MAPK信號通路、光合作用—天線蛋白、甜菜紅堿生物合成、植物激素信號轉(zhuǎn)導(dǎo)、氨基糖和核苷酸糖的代謝及光合作用（圖6-A）。分析發(fā)現(xiàn)在相反響應(yīng)模式的基因中沒有可以被顯著富集到的KEGG通路（圖6-B）。綜上分析表明，在統(tǒng)計學(xué)角度水稻內(nèi)存在6個對2種害蟲均具有一致響應(yīng)模式的通路，無響應(yīng)模式相反的通路。

2. 5 對2種害蟲具有相同響應(yīng)模式的基因分析結(jié)果

篩選出轉(zhuǎn)錄組中對2種害蟲有相同響應(yīng)模式的基因，分為共上調(diào)基因和共下調(diào)基因進行表達量可視化及功能注釋。基因的功能注釋由水稻基因組數(shù)據(jù)庫網(wǎng)站（http：//rice.uga.edu/）獲得。利用log2 FC的值來表示表達量變化倍數(shù)的高低，上調(diào)為正值，下調(diào)為負值。由圖5可明顯看出，在共上調(diào)基因中，LOC_Os03g64260和LOC_Os11g47600表達量增加倍數(shù)顯著高于其他基因，前者被注釋為含AP2結(jié)構(gòu)域的蛋白（AP2 domain containing protein），對褐飛虱的響應(yīng)程度更高，而后者被注釋為糖苷水解酶（Glycoside hydrolase），對稻癭蚊的響應(yīng)程度更高。共下調(diào)基因中的LOC_Os10g40520.1被注釋為脂質(zhì)轉(zhuǎn)運家族蛋白前體145（LTPL145），在稻癭蚊取食后表達量顯著下調(diào)，下調(diào)倍數(shù)明顯高于其他共下調(diào)基因。被注釋為表達蛋白（Expressed protein）的基因無論是在共上調(diào)還是共下調(diào)基因中均數(shù)量最多。除了表達蛋白，在共上調(diào)基因中數(shù)量較多的還有被注釋為以下蛋白質(zhì)的基因：含有AP2結(jié)構(gòu)域的蛋白（AP2 domain containing protein）、鈣調(diào)素結(jié)合蛋白（Calmodulin binding protein）、糖苷水解酶（Glycoside hydrolase）和含有ZIM結(jié)構(gòu)域的蛋白質(zhì)（ZIM domain containing protein）（圖7-A）;在下調(diào)基因中注釋為植物脂質(zhì)轉(zhuǎn)運家族蛋白前體的基因數(shù)量最多（LTPL145、LTPL20、LTPL52和LTPL66）（圖7-B）。

3 討論

通過分析褐飛虱和稻癭蚊取食后苗期水稻的轉(zhuǎn)錄組情況，篩選得到3963個褐飛虱取食誘導(dǎo)的差異表達基因和1206個稻癭蚊取食誘導(dǎo)的差異表達基因。響應(yīng)褐飛虱的差異表達基因數(shù)較稻癭蚊的多1倍以上，可能與稻苗的處理時長和傷害水平有關(guān)。褐飛虱處理的幼苗取樣時間短，僅24 h，而稻癭蚊處理為7 d后取樣，植物擁有更長的緩沖平衡時間。不同的取樣時間點，對差異基因數(shù)量的獲得具有較大影響，如水稻品種Taichuang Native 1在褐飛虱取食6 h后引起1104個基因表達差異，取食48 h后則有2523個差異表達基因（曹征鴻等，2020;Liu et al.，2021）。除此之外，褐飛虱對水稻的傷害比稻癭蚊嚴重，同時接入大量的褐飛虱若蟲致使幼苗3 d內(nèi)完全枯死，而稻癭蚊若蟲僅啃食咬小部分組織。因此，褐飛虱取食于水稻而言是突發(fā)的致死性傷害，導(dǎo)致應(yīng)激的基因多于稻癭蚊處理。

GO功能注釋分析表明響應(yīng)褐飛虱取食后差異基因的生物學(xué)功能最顯著富集在細胞壁合成和缺水響應(yīng)，而稻癭蚊取食后差異基因的生物學(xué)功能最顯著富集在光合作用上。差異基因富集的細胞定位與生物學(xué)功能具有密切聯(lián)系，如褐飛虱的差異基因富集于細胞壁和細胞外周，而響應(yīng)稻癭蚊的差異基因富集在光合作用相關(guān)的細胞器，如類囊體、葉綠體和質(zhì)體等。究其原因與水稻應(yīng)對這2種脅迫的特征有關(guān)，褐飛虱是刺吸式昆蟲，吸食水稻汁液導(dǎo)致稻苗缺水枯死，因此引起大量細胞壁合成相關(guān)基因響應(yīng)和稻苗的缺水反應(yīng)（Zhang et al.，2004）;稻癭蚊的取食則促使水稻發(fā)育形成標蔥，影響植株形態(tài)發(fā)育，進而影響光合作用。已知的抗褐飛虱基因bph6和bph30被證實為通過增加細胞壁硬度和厚壁組織的厚度阻礙褐飛虱取食，前者通過胞泌蛋白調(diào)控細胞壁的形成，后者則通過上調(diào)薄壁組織纖維素與半纖維素相關(guān)基因增加厚壁組織細胞壁中纖維素和半纖維素的積累（Guo et al.，2018;Shi et al.，2021）。而水稻對咀嚼式昆蟲一直缺乏有效的抗性資源，因此涉及的抗蟲機理存在大量空白。但已有大量研究表明致癭昆蟲普遍削減寄主植物的光合作用，Huang等（2011，2015）通過測量葉綠素?zé)晒獍l(fā)現(xiàn)紅楠和長葉木姜子癭蚊誘導(dǎo)的蟲癭部位光合效率顯著低于健康部位。然而也有部分致癭昆蟲對植物沒有顯著影響或增強光合作用，如橄欖星室木虱和癭蜂（Fay et al.，1993;Malenovsky et al.，2015）。此外，致癭昆蟲在蟲癭內(nèi)完成整個生活史，對植物產(chǎn)生的影響豐富多樣，研究表明致癭昆蟲還具有使植物細胞壁變薄、胞間連絲增大、細胞器變多等植物細胞重編程的功能（Florentine et al.，2005;Giron et al.，2016;楊孟可等，2021）。本研究中的GO富集分析結(jié)果與前人研究結(jié)論一致，在感性品種9311中，與細胞壁合成和光合作用相關(guān)的基因更多地呈現(xiàn)下調(diào)趨勢，因此研究細胞壁的形成與水稻抗褐飛虱取食之間的關(guān)系有利于揭示水稻抗褐飛虱機理，而稻癭蚊還需進一步研究探索增強植物光合作用以及增強植物細胞壁強度是否有助于抗蟲功能。稻癭蚊引起的水稻標蔥是蟲癭的一種（Agarrwal et al.，2014），大部分關(guān)于致癭昆蟲與植物間互作關(guān)系的研究集中于林木，稻癭蚊引起的水稻蟲癭未有深入研究，利用水稻作為模式植物進一步探明水稻標蔥的形成與稻癭蚊取食之間的關(guān)系是未來的研究方向，研究結(jié)果可為植物響應(yīng)生物脅迫和致癭昆蟲與植物協(xié)同進化提供理論支撐。C3EF4567-7C12-48BB-84A9-C29F2CF52217

KEGG信號通路富集分析表明2種害蟲取食均引起有關(guān)代謝通路、次生代謝物生物合成和苯丙素生物合成的顯著變化。次生代謝物是植物長期進化中應(yīng)對環(huán)境和脅迫所產(chǎn)生的結(jié)果，在植物對昆蟲防御方面發(fā)揮著重要作用（張文輝和劉光杰，2003）。其中揮發(fā)性的次生代謝物可誘集植食性昆蟲的寄生蜂和天敵，間接地進行植物防御（Howe and Jander，2008）。水稻氫過氧化物裂解酶OsHPL3通過調(diào)節(jié)羥磷脂途徑，吸引稻虱纓小蜂寄生從而提高水稻對褐飛虱的抗性（Tong et al.，2012）。但最新研究表明植物次生代謝物不僅能被植物用于間接防御蟲害，反之亦可被昆蟲用于抵御天敵，如褐飛虱偏愛與二化螟共享寄主是因為褐飛虱可借助二化螟誘導(dǎo)水稻揮發(fā)物而規(guī)避稻虱纓小蜂對其卵的寄生（Hu et al.，2020）;該團隊隨后又發(fā)現(xiàn)褐飛虱也能幫助二化螟實現(xiàn)對水稻的反防御，褐飛虱取食不僅可下調(diào)二化螟單獨取食誘導(dǎo)的防御基因，而且能降低水稻利用信息揮發(fā)物吸引稻螟赤眼蜂，甚至做到排斥，極大降低二化螟卵的寄生率（Liu et al.，2021）。褐飛虱與二化螟作為水稻最主要的兩大害蟲，二者間可形成種間互利關(guān)系，攜手對抗水稻。由此可見，植物的次生代謝物合成與昆蟲間的關(guān)系非常復(fù)雜，褐飛虱與稻癭蚊之間是否也具有互利的關(guān)系值得深入研究，為將來制定害蟲生態(tài)防控打下理論基礎(chǔ)。

對差異基因進行維恩分析，發(fā)現(xiàn)水稻對褐飛虱和稻癭蚊的共響應(yīng)基因有251個，其中108個基因表達模式相似，143個表達模式相反。對共響應(yīng)基因進行富集分析，發(fā)現(xiàn)具有相同響應(yīng)模式的差異表達基因具有多個顯著富集的GO功能條目和6個KEGG信號通路，具有相反模式的差異表達基因僅有3個顯著富集的GO功能條目。因此，水稻內(nèi)存在較多的對2種害蟲共響應(yīng)的生物學(xué)功能和信號通路，尋找到一個能共同調(diào)控褐飛虱和稻癭蚊抗性的通路具有可行性。具有相同表達模式的差異表達基因在生物學(xué)功能上主要與光合作用、激素信號轉(zhuǎn)導(dǎo)和MAPK信號通路有關(guān)，其中激素調(diào)節(jié)主要涉及水楊酸的生物合成和茉莉酸信號途徑。已有大量研究表明茉莉酸和水楊酸參與植物響應(yīng)生物脅迫，調(diào)控褐飛虱的抗性（李毅等，2018;Xu et al.，2021）。如褐飛虱取食激活了水楊酸依賴的激素信號通路，并通過水楊酸信號通路抑制胼胝質(zhì)降解，從而增強對褐飛虱的抗性（Ling and Weilin，2016）。另外，研究表明，MAPK參與調(diào)控植株體內(nèi)的水楊酸與茉莉酸信號，在水稻信號傳導(dǎo)中發(fā)揮重要作用（Yuan et al.，2005;Wang et al.，2008;Hu et al.，2011）。Li等（2019）研究表明，OsMAPK20-5基因負調(diào)控水稻對褐飛虱的抗性，該基因的表達量在褐飛虱懷孕的雌蟲和成蟲取食過程中上調(diào)迅速，將該基因沉默后顯著提高植株乙烯和一氧化氮的積累，從而提高植株對褐飛虱成蟲的抗性并降低褐飛虱的產(chǎn)卵水平。關(guān)于激素信號通路和MAPK信號通路在褐飛虱抗性研究中較為廣泛，因這2個通路在2種害蟲中具有相同的響應(yīng)模式，因此后續(xù)在稻癭蚊的抗蟲機理解析中可借鑒褐飛虱的研究進展進行綜合考量。具有相反表達模式的基因的GO功能注釋，3個條目均與乙醛酸循環(huán)相關(guān)，乙醛酸循環(huán)可實現(xiàn)脂肪到糖的轉(zhuǎn)變，褐飛虱取食增強感性材料的乙醛酸循環(huán)（Peng et al.，2016），而稻癭蚊取食則降低乙醛酸循環(huán)相關(guān)基因的表達。此外，代謝組檢測稻癭蚊取食48 h后糖類的含量，其中果糖、葡萄糖和半乳糖均顯著降低，抗性材料的下降程度比感性材料更劇烈（陳燦等，2021）。推測出現(xiàn)這種情形的原因與2種害蟲對植物的傷害相關(guān)，褐飛虱取食后改變了植物體的滲透壓，需要更多的糖分來維持平衡，而稻癭蚊取食后，植物快速生長形成蟲癭，消耗大量糖分。因此，探清乙醛酸循環(huán)在植物抗性機制內(nèi)的功能尚需更深入的研究。

4 結(jié)論

通過水稻幼苗的轉(zhuǎn)錄組分析，篩選出3963個響應(yīng)褐飛虱取食的差異表達基因，1206個響應(yīng)稻癭蚊取食的差異表達基因，其中有251個基因?qū)?種害蟲的取食共響應(yīng)，這些共響應(yīng)基因主要富集在光合作用、激素信號傳導(dǎo)、次生代謝物合成和乙醛酸循環(huán)等通路上。表明水稻幼苗響應(yīng)這2種害蟲具有相同點和不同之處，推測這些共響應(yīng)通路在水稻廣譜抗蟲上可能發(fā)揮著重要作用。

參考文獻：

曹征鴻，賀康，徐樂，湯沈楊，王亞琴，李飛. 2020. 苗期水稻響應(yīng)褐飛虱取食的基因差異表達分析[J]. 應(yīng)用昆蟲學(xué)報，57（4）：911-920. [Cao Z H，He K，Xu L，Tang S Y，Wang Y Q，Li F. 2020. Change in the gene expression of seedling stage rice in response to feeding by the brown planthopper Nilaparvata lugens （St?l）（Hemiptera： Delphaci-dae）[J]. Chinese Journal of Applied Entomology 57（4）：911-920.] doi： 10.7679/j.issn.2095-1353. 2020.093.

陳燦，張宗瓊，夏秀忠，楊行海，農(nóng)保選，張曉麗，徐志健，李丹婷，郭輝，馮銳. 2021. 代謝組學(xué)揭示水稻—稻癭蚊互作的潛在生化標識物[J]. 南方農(nóng)業(yè)學(xué)報，52（10）：2662-2670. [Chen C，Zhang Z Q，Xia X Z，Yang X H，Nong B X，Zhang X L，Xu Z J，Li D T，Guo H，F(xiàn)eng R. 2021. Metabolomics reveals potential biomarkers of rice gall midge（Orseolia oryzae） interactions[J]. Journal of Sou-thern Agriculture，52（10）：2662-2670.] doi：10.3969/j.issn.2095-1191.2021.10.005.C3EF4567-7C12-48BB-84A9-C29F2CF52217

杜波，陳榮智，何光存. 2018. 抗褐飛虱基因的發(fā)掘、鑒定與利用[J]. 生命科學(xué)，30（10）：1072-1082. [Du B，Chen R Z，He G C. 2018. The progress of functional genomics research of rice resistance to insect[J]. Chinese Bulletin of Life Sciences，30（10）：1072-1082.] doi： 10.13376/j.cbls/2016158.

黃鳳寬，韋素美，梁廣文，黃所生，蔣顯斌，羅善昱，李青. 2006. 水稻品種RP1976-18-6-4-2對褐飛虱和稻癭蚊的抗性評價及其遺傳分析[J]. 中國水稻科學(xué)，20（1）：113-115. [Huang F K，Wei S M，Liang G W，Huang S S，Jiang X B，Luo S Y，Li Q. 2006. Genetic analysis on the resistance to the brown planthopper and rice gall midge in rice variety RP1976-18-6-4-2[J]. Chinese Journal of Rice Science，20（1）：113-115.] doi： 10.16819/j.1001-7216.2006.01.021.

李毅，張嘉嬌，杜波，何光存，李家儒. 2018. 水稻與褐飛虱化學(xué)關(guān)系的研究進展[J]. 植物生理學(xué)報，54（4）：528-538. [Li Y，Zhang J J，Du B，He G C，Li J R. 2018. Research progress of chemical interactions between rice and brown planthopper[J]. Plant Physiology Journal，54（4）：528-538.] doi： 10.13592/j.cnki.ppj.2017.0530.

秦學(xué)毅，鄧國富，馮銳，郭輝，劉百龍. 2016. 稻褐飛虱和稻癭蚊雙抗種質(zhì)的抗性基因精細定位及創(chuàng)新利用[Z]. 國家科技成果. [Qin X Y，Deng G F，F(xiàn)eng R，Guo H，Liu B L. 2016. Fine mapping and innovative utilization of resistance genes of rice brown planthopper and rice gall midge double resistant germplasm[Z]. National Scientific and Technological Achievements.]

吳碧球，黃鳳寬，黃所生，韋素美. 2014. 水稻品種ARC5833對稻癭蚊的抗性評價及其遺傳分析[J]. 西南農(nóng)業(yè)學(xué)報，27（4）：1483-1487. [Wu B Q，Huang F K，Huang S S，Wei S M. 2014. Identification and genetic analysis of gall midge resistance in rice variety ARC5833[J]. Southwest China Journal of Agricultural Sciences，27（4）：1483-1487.] doi： 10.16213/j.cnki.scjas.2014.04.017.

鄢柳慧，黃福鋼，舒宛，張啟明，趙能，莫心怡，王馳，曹凱，邱永福. 2020. 秈稻品種‘ARC5833抗褐飛虱基因的遺傳分析與定位[J]. 分子植物育種，18（18）：6038-6043. [Yan L H，Huang F G，Shu W，Zhang Q M，Zhao N，Mo X Y，Wang C，Cao K，Qiu Y F. 2020. Genetic analysis and mapping of brown planthopper resistance gene in indica rice variety ‘ARC5833[J]. Molecular Plant Breeding，18（18）：6038-6043.] doi： 10.13271/j.mpb.018.006038.

嚴火其. 2021. 農(nóng)業(yè)害蟲危害何以越來越嚴重[J]. 中國農(nóng)史，40（3）：3-15. [Yan H Q. 2021. Why are agricultural pests becoming more and more harmful[J]. Agricultural History of China，40（3）：3-15.]

楊利艷，張玉榮，楊雅舒，王美霞，陳保國，趙麗，張麗光，王創(chuàng)云. 2020. 基于RNA-seq數(shù)據(jù)分析玉米抗蟲響應(yīng)基因的可變剪接事件[J]. 河南農(nóng)業(yè)大學(xué)學(xué)報，54（2）：181-188. [Yang L Y，Zhang Y R，Yang Y S，Wang M X，Chen B G，Zhao L，Zhang L G，Wang C Y. 2020. Analysis of alternative splicing events of insect-resistant response genes based on RNA-seq data in Zea mays[J]. Journal of Henan Agricultural University，54（2）：181-188.] doi：10.16445/j.cnki.1000-2340.20200403.009.

楊孟可，劉賽，喬海莉，郭昆，徐榮，徐常青，陳君. 2021. 致癭昆蟲對寄主植物生理和代謝的影響[J]. 昆蟲學(xué)報，64（4）：536-548. [Yang M K，Liu S，Qiao H L，Guo K，Xu R，Xu C Q，Chen J. 2021. Influences of gall-inducing insects on the physiology and metabolism of host plants[J]. Acta Entomologica Sinica，64（4）：536-548.] doi： 10.16380/j.kcxb.2021.04.012.C3EF4567-7C12-48BB-84A9-C29F2CF52217

張文輝，劉光杰. 2003. 植物抗蟲性次生物質(zhì)的研究概況[J]. 植物學(xué)通報，20（5）：522-530. [Zhang W H，Liu G J. 2003. A review on plant secondary substances in plant resistance to insect pests[J]. Chinese Bulletin of Botany，20（5）：522-530.] doi：10.3969/j.issn.1674-3466.2003.05.002.

Agarrwal R，Bentur J S，Nair S. 2014. Gas chromatography mass spectrometry based metabolic profiling reveals biomarkers involved in rice-gall midge interactions[J]. Journal of Integrative Plant Biology，56（9）：837-848. doi： 10.1111/jipb.12244.

Fay P A，Hartnett D C，Knapp A K. 1993. Increased photosynthesis and water potentials in Silphium integrifolium galled by cynipid wasps[J]. Oecologia，93（1）：114-120. doi：10.1007/bf00321200.

Florentine S K，Raman A，Dhileepan K. 2005. Effects of gall induction by Epiblema strenuana on gas exchange，nutrien-ts，and energetics in Parthenium hysterophorus[J]. Biocontrol，50（5）：787-801. doi：10.1007/s10526-004-5525-3.

Giron D，Huguet E，Stone G N，Body M. 2016. Insect-induced effects on plants and possible effectors used by galling and leaf-mining insects to manipulate their host-plant[J]. Journal of Insect Physiology，84：70-89. doi：10.1016/j.jinsphys.2015.12.009.

Guo J P，Xu C X，Wu D，Zhao Y，Qiu Y F，Wang X X，Ou-yang Y D，Cai B D，Liu X，Jing S L，Shangguan X X，Wang H Y，Ma Y H，Hu L，Wu Y，Shi S J，Wang W L，Zhu L L，Xu X，Chen R Z，F(xiàn)eng，Y Q，Du B，He G C. 2018. Bph6 encodes an exocyst-localized protein and confers broad resistance to planthoppers in rice[J]. Nature Genetics，50（2）：297-306. doi：10.1038/s41588-018-0039-6.

Howe G A，Jander G. 2008. Plant immunity to insect herbivores[J]. Annual Review of Plant Biology，59：41-66. doi：10.1146/annurev.arplant.59.032607.092825.

Hu J，Zhou J B，Peng X X，Xu H H，Liu C X，Du B，Yuan H Y，Zhu L L，He G C. 2011. The Bphi008a gene interacts with the ethylene pathway and transcriptionally regulates MAPK genes in the response of rice to brown planthopper feeding[J]. Plant Physiology，156（2）：856-872. doi：10.1104/pp.111.174334.

Hu X Y，Su S L，Liu Q S，Jiao Y Y，Peng Y F，Li Y H，Turlings T C. 2020. Caterpillar-induced rice volatiles provide enemy-free space for the offspring of the brown planthopper[J]. Elife，9：e55421. doi：10.7554/eLife.55421.

Huang M Y，Huang W D，Chou H M，Chen C C，Chen P J，Chang Y T，Yang C M. 2015. Structural，biochemical，and physiological characterization of photosynthesis in leaf-derived cup-shaped galls on Litsea acuminata[J]. BMC Plant Biology，15（1）：61. doi：10.1186/s12870-015- 0446-0.

Huang M Y，Lin K H，Yang M M，Chou H M，Yang C M，Chang Y T. 2011. Chlorophyll fluorescence，spectral proper-ties，and pigment composition of galls on leaves of Machilus thunbergii[J]. International Journal of Plant Sciences，172（3）：323-329. doi：10.1086/658157.C3EF4567-7C12-48BB-84A9-C29F2CF52217

Li J C，Liu X L，Wang Q，Huangfu J Y，Schuman M C，Lou，Y G. 2019. A group D MAPK protects plants from autotoxicity by suppressing herbivore-induced defense signaling[J]. Plant Physiology，179（4）：1386-1401. doi：10.1104/pp. 18.01411.

Li R，Zhang J，Li J C，Zhou G X，Wang Q，Bian W B，Erb M，Lou Y G. 2015. Prioritizing plant defence over growth through WRKY regulation facilitates infestation by non-target herbivores[J]. Elife，4：e04805. doi：10.7554/elife. 04805.020.

Li Y，Mo Y，Li Z H，Yang M，Tang L H，Cheng L，Qiu Y F. 2020. Characterization and application of a gall midge resistance gene （Gm6） from Oryza sativa ‘Kangwenqingzhan[J]. Theoretical and Applied Genetics，133（2）：579-591. doi：10.1007/s00122-019-03488-w.

Ling Y，Weilin Z. 2016. Genetic and biochemical mechanisms of rice resistance to planthopper[J]. Plant Cell Reports，35（8）：1559-1572. doi：10.1007/s00299-016-1962-6.

Liu C X，Hao F H，Hu J，Zhang W L，Wan L L，Zhu L L，Tang H R，He G C. 2010. Revealing different systems responses to brown planthopper infestation for pest susceptible and resistant rice plants with the combined metabonomic and gene-expression analysis[J]. Journal of Proteome Research，9（12）：6774-6785. doi：10.1021/pr100 970q.

Liu Q S，Hu X Y，Su S L，Ning Y S，Peng Y F，Ye G Y，Lou Y G，Turlings T C J，Li Y H. 2021. Cooperative herbivory between two important pests of rice[J]. Nature Communications，12（1）：6772. doi：10.1038/s41467-021-270 21-0.

Lu J，Ju H P，Zhou G X，Zhu C S，Erb M，Wang X P，Wang P，Lou Y G. 2011. An EAR-motif-containing ERF transcription factor affects herbivore-induced signaling，defense and resistance in rice[J]. The Plant Journal，68（4）：583-596. doi：10.1111/j.1365-313x.2011.04709.x.

Malenovsky I，Burckhardt D，Queiroz D L，Isaias R M S，Oliveira D C. 2015. Descriptions of two new Pseudophacopteron species（Hemiptera：Psylloidea：Phacopteronidae）inducing galls on Aspidosperma（Apocynaceae）in Brazil[J]. Acta Entomologica Musei Nationalis Pragae，55（2）：513-538.

Peng L，Zhao Y，Wang H Y Zhang J J，Song C P，Shangguan X X，Zhu L L，He G C. 2016. Comparative metabolomics of the interaction between rice and the brown planthopper[J]. Metabolomics，12（8）：1-15. doi：10.1007/s11306-016-1077-7.

Shi S J，Wang H Y，Nie L Y，Tan D，Zhou C，Zhang Q，Li Y，Du B，Guo J P，Huang J，Wu D，Zheng X H，Guan W，Shan J H，Zhu L L，Chen R Z，Xue L J，Walling L L，He G C. 2021. Bph30 confers resistance to brown planthopper by fortifying sclerenchyma in rice leaf sheaths[J]. Molecular Plant，14（10）：1714-1732. doi：10.1016/j.molp. 2021.07.004.C3EF4567-7C12-48BB-84A9-C29F2CF52217

Tong X H，Qi J F，Zhu X D，Mao B Z，Zeng L J，Wang B H，Li Q，Zhou G X，Xu X J，Lou Y G，He Z H. 2012. The rice hydroperoxide lyase OsHPL3 functions in defense responses by modulating the oxylipin pathway[J]. The Plant Journal，71（5）：763-775. doi：10.1111/j.1365-313x. 2012.05027.x.

Wang Y Y，Wang X L，Yuan H Y，Chen R Z，Zhu L L，He R F，He G C. 2008. Responses of two contrasting genotypes of rice to brown planthopper[J]. Molecular Plant-Microbe Interactions，21（1）：122-132. doi：10.1094/mpmi-21-1-0122.

Xu J，Wang X J，Zu H Y，Zeng X，Baldwin I T，Lou Y G，Li R. 2021. Molecular dissection of rice phytohormone signaling involved in resistance to a piercing-sucking herbivore[J]. New Phytologist，230（4）：1639-1652. doi：10. 1111/nph.17251.

Yuan H Y，Chen X P，Zhu L L，He G C. 2005. Identification of genes responsive to brown planthopper Nilaparvata lugens St?l （Homoptera： Delphacidae） feeding in rice[J]. Planta，221（1）：105-112. doi：10.1007/s00425-004-1422-3.

Zhang F，Zhu L，He G C. 2004. Differential gene expression in response to brown planthopper feeding in rice[J]. Journal of Plant Physiology，161（1）：53-62. doi：10.1078/0176-1617-01179.

Zhao Y，Huang J，Wang Z Z，Jing S L，Wang Y，Ouyang Y D，Cai B D，Xin X F，Liu X，Zhang C X，Pan Y F，Ma R，Li Q F，Jiang W H，Zeng Y，Shangguan X X，Wang H Y，Du B，Zhu L L，Xu X，F(xiàn)eng Y Q，He S Y，Chen R Z，Zhang Q F，He G C. 2016. Allelic diversity in an NLR gene BPH9 enables rice to combat planthopper variation[J]. Proceedings of the National Academy of Sciences，113（45）：12850-12855. doi：10.1073/pnas.1614862113.

Zhou G X，Qi J F，Ren N，Cheng J A，Erb M，Mao B Z，Lou Y G. 2009. Silencing OsHI-LOX makes rice more susceptible to chewing herbivores，but enhances resistance to a phloem feeder[J]. The Plant Journal，60（4）：638-648. doi：10.1111/j.1365-313x.2009.03988.x.

（責(zé)任編輯 麻小燕）C3EF4567-7C12-48BB-84A9-C29F2CF52217