海南本地抗病小米辣響應象耳豆根結線蟲侵染的miRNA鑒定及其靶基因預測

摘" 要：象耳豆根結線蟲（Meloidogyne enterolobii）是一種對辣椒生產危害較大的病原物，使用抗病品種是防治該病害最為經濟、環保的手段。深入理解根結線蟲抗性調控機理將有利于加速抗病育種進程。本課題組前期研究中對海南本地抗病小米辣種質CF25進行了接種前后的轉錄組數據分析，挖掘了與根結線蟲抗性相關的代謝通路。為進一步探索抗性表達調控機理，本研究對CF25接種前后sRNA文庫進行了高通量測序，共檢測到133個差異表達miRNAs。以P≤0.05和|log₂（FC）|≥3為篩選標準，鑒定出33個顯著差異表達的miRNAs，其中包括19個已知miRNAs和14個新發現miRNAs。對33個顯著差異表達的miRNAs進行靶基因預測，共得到373個靶基因，其中miR5658-z預測靶基因數量最多。GO功能和KEGG代謝通路富集分析顯示靶基因主要富集于植物-病原菌互作、信號轉導等通路，推測miRNA介導的信號轉導及植物防御反應可能是CF25抗象耳豆根結線蟲的重要原因。為初步驗證差異表達miRNA和靶基因間的調控關系，選取6個差異表達miRNAs及其靶基因進行qRT-PCR表達驗證，得出miRNA與靶基因的表達差異趨勢與高通量測序結果一致，基本符合miRNA負向調控靶基因表達的規律。以上研究結果表明miRNA可能在辣椒對根結線蟲的防御反應中起到了重要作用，為后續深入解析miRNA介導的根結線蟲抗性機理奠定基礎。

關鍵詞：灌木辣椒；抗象耳豆根結線蟲；sRNA測序；miRNA；靶基因

中圖分類號：S436.418 """""文獻標志碼：A

Identification and Target Gene Prediction of miRNA Responding to the Infection of Meloidogyne enterolobii in Local Resistant Xiaomila Pepper of Hainan， China

GAO Chang¹， ZHANG Zhiyuan²， LIU Ziji³， CHEN Meihong⁴， ZHU Jie^1*， CAO Zhenmu^3*

<！--[if ！supportLists]-->1. <！--[endif]-->School of Tropical Agriculture and Forestry （School of Agriculture and Rural Affairs / School of Rural Revitalization）， Hainan University / Hainan Key Laboratory for Quality Regulation of Tropical Horticultural Crops， Haikou， Hainan 570228， China; 2. Hainan Research Institute， Zhejiang University， Sanya， Hainan 572025， China; 3. Institute of Tropical Crops Genetic Resources， Chinese Academy of Tropical Agricultural Sciences， Haikou， Hainan 571101， China; 4. Hainan Binfen Horticulture Co.， Ltd.， Haikou， Hainan 571133， China

Abstract： Meloidogyne enterolobii is a pathogen causing significant damages to pepper production. Utilization of resistant cultivars is the most economical and environmentally friendly method to control this disease. A deeper understanding of the regulatory mechanisms of resistance to root-knot nematodes could facilitate the breeding process of disease-resistant cultivars. In previous studies， analysis of transcriptomic data of Hainan local M. enterolobii-resistant pepper germplasm CF25 before and after inoculation were carried out to explore metabolic pathways closely related to nematode resistance. In order to further explore the regulatory mechanisms of resistance， sRNA libraries of CF25 before and after inoculation with M. enterolobii were sequenced to identify a total of 133 differentially expressed miRNAs. Based on the criteria of P≤0.05 and |log₂（FC）|≥3， 33 significantly differentially expressed miRNAs were identified， including 19 known miRNAs and 14 newly discovered miRNAs. Target gene prediction for the 33 miRNAs yielded a total of 373 target genes， with miR5658-z having the highest number of predicted target genes. GO and KEGG analyses showed that the target genes were mainly enriched in pathways related to plant-pathogen interactions and signal transduction. It is speculated that miRNA-mediated signal transduction and plant defense responses might be important reasons for the resistance of CF25 to M. enterolobii. To preliminarily verify the regulatory relationships between differentially expressed miRNA and the target genes， six differentially expressed miRNA and the target genes were selected for qRT-PCR validation. The results showed that the differential expression trends of miRNA and target genes were consistent with the high-throughput sequencing results and generally conformed to the negative regulation pattern between miRNA and the target genes. The results of this study suggest that miRNA may play significant roles in the defense response to root-knot nematodes and lay the foundation for furtherly elucidating the mechanisms of miRNA-mediated nematode resistance in pepper.

Keywords： Capsicum frutescens; Meloidogyne enterolobii resistance; sRNA sequencing; miRNA; target gene

DOI： 10.3969/j.issn.1000-2561.2025.07.005

辣椒（Capsicum spp.）是茄科（Solanaceae）辣椒屬（Capsicum）一年生或多年生草本植物，距今已經有6000多年的栽培歷史^[1]。辣椒種植方法簡單、產值高，果實富含辣椒素、類胡蘿卜素和類黃酮等營養成分，不僅能滿足人們對調味的需求，還在醫藥、軍事、農業方面都有一定應用價值^[2]。由于全球氣候改變和農業文明現代化等多種因素的影響，辣椒面臨多種病蟲害的威脅。其中根結線蟲（Meloidogyne spp.）被認為是一類重要的辣椒土傳病害，可導致其產量和品質嚴重下降^[3]。

根結線蟲是一類透明無脊椎植物病原物，能侵染植物根部使其形成膨大根結，從而破壞根系對水分和養分的吸收，影響植物的生長發育，甚至導致植物提早死亡。根結線蟲種類繁多，其中象耳豆根結線蟲（Meloidogyne enterolobii）廣泛分布于熱帶、亞熱帶地區，對辣椒的致病性和侵染能力較強，且已經克服N、Me1、Me3/Me7等辣椒抗根結線蟲基因^[4-5]，在辣椒生產上造成的經濟損失日益嚴重^[6]，因此需要深入研究辣椒抗病機理，發掘新的抗病基因，為抗病育種工作提供理論依據和基因資源。

microRNA（miRNA）是一類具有調控功能的內源性非編碼小分子RNA，其長度為21～24 nt。miRNA可以通過堿基互補配對識別相關靶基因mRNA，進而降解該mRNA或者抑制其翻譯，達到負調控基因表達的效果^[7]。miRNA廣泛參與各植物生物學過程，在防御應答以及寄主-病原物互作中均起到重要作用。有研究發現馬鈴薯遭大麗輪枝菌（Verlicillium dahliae）侵染后，miR482下調表達，負向調控NBS-LRR基因上調表達，從而提高其抗病性^[8]。水稻感染稻瘟病菌（Magna-porthe oryzae）后，miR398b下調表達，從而調控超氧化物歧化酶基因上調表達，加強了對稻瘟病的抗性^[9]。由miR397靶向的漆酶基因CA07g11100在辣椒疫霉菌（Phytophthora capsici）侵染的前期上調，而后轉為下調，推測該基因可能參與辣椒與疫霉菌的互作^[10]。由此可見，miRNA可以通過調節靶基因的表達來影響植物對病原物的抗性^[11]。然而，目前有關辣椒抗根結線蟲miRNA-靶基因的研究工作相對較少，從miRNA層面探究根結線蟲與辣椒互作的分子機制將加深人們對作物抗根結線蟲機理的認識，為后續抗病育種工作提供較為新穎的切入點。

小米辣原產于云南、海南等地區的野生、半野生灌木辣椒（Capsicum frutescens）。其植株粗大近似鋤柄，根系發達，結果期長，產量較高^[12]。在長期自然選擇的壓力下，小米辣逐步演化出了耐高溫潮濕、抗病、耐瘠薄和耐弱光的特性^[13]，對許多生物及非生物脅迫都具有良好的抗性，能生長在極其惡劣的環境條件下^[14-18]。為了改良我國現有辣椒栽培種，挖掘小米辣在抗病育種中的潛力，本課題組前期從21份海南本地小米辣資源中篩選出7份象耳豆根結線蟲高抗種質。通過酸性品紅染色和石蠟切片等組織病理學方法，發現高抗種質CF25對象耳豆根結線蟲的抗性主要體現在抗侵入和限制取食兩方面。通過對CF25和感病對照CF29進行接種前后的防御酶測定，得出CF25在接種后多數時間點的PAL和POD活性顯著高于感病對照，推測較高的PAL和POD活性可能導致包括木質素在內的更多苯丙烷類次生代謝產物的積累。本課題組前期發現接種第24小時，抗、感材料中線蟲數量分別為（0.67±0.94）條和（0.33±0.47）條，不存在顯著差異（P=0.68）。然而，當接種第4天時，CF25單株根系中僅發現（0.67±0.94）條線蟲，而感病對照CF29中約存在（11±3.27）條線蟲，數量差異達到了顯著水平（P=0.01）。不僅如此，通過觀察接種第1、2、3、4天的線蟲形態，發現接種第4天線蟲略有加粗，推測已開始建立取食點。為兼顧對阻止侵入和抑制取食分子機理的研究，本課題組對高抗材料CF25和感病對照CF29進行接種前后（第0天，第4天）的轉錄組測序分析，得出抗、感材料接種前后差異基因主要富集于苯丙烷生物合成、植物-病原體互作等代謝通路^[19-20]。然而，如何從這些代謝通路中篩選出與抗性緊密關聯的基因是一個亟待解決的問題，且相關基因的表達調控機理尚未得到解析。為此，本研究以高抗象耳豆根結線蟲種質CF25為研究對象，構建接種前后（接種第0天和第4天）sRNA文庫，利用Illumina二代測序技術開展miRNA測序，篩選差異表達miRNA并進行靶基因的預測和初步驗證，為進一步研究辣椒中miRNA介導的根結線蟲抗性機制提供重要依據。

1" 材料與方法

1.1" 材料

選取海南本地收集的象耳豆根結線蟲抗性小米辣資源CF25為試驗對象。將種子在55"℃溫水中浸泡15 min，28"℃恒溫催芽。種子露白后播種到用滅菌土填充的21孔穴盤中，每穴播種1粒種子。待辣椒苗長出6片真葉，對每株種苗接種1000條象耳豆根結線蟲二齡幼蟲，在接種第0天（R0，對照）和第4天（R4）分別選取大小均一的種苗各3株，將根系洗凈，迅速用濾紙吸干水分，投入液氮速凍，送至廣州基迪奧生物科技有限公司進行sRNA文庫構建和高通量測序。

1.2" 方法

1.2.1" miRNA的提取和文庫構建 "用Trizol法從樣本中提取總RNA，聚丙烯酰胺凝膠（PAGE）電泳切膠選擇18～30 nt范圍的條帶，回收sRNA。根據sRNA的結構特點，分別連接3￠和5￠接頭，對sRNA進行反轉錄和PCR擴增，最后使用PAGE膠回收并純化條帶，完成文庫構建。構建好的文庫使用Agilent2100以及ABI StepOnePlus Real-Time PCR System（Life Technologies）進行質量和產量的檢測，合格后上機測序。

1.2.2" miRNA鑒定" 對測序原始數據（raw data）進行清洗，去除不帶有3'接頭的、測序效果不好的、不含插入片段或插入片段長度小于18 nt的、70%以上堿基為poly A的低質量序列，得到用于后續分析的純凈序列（clean reads）。將純凈序列與中國辣椒（Capsicum chinense）參考基因組數據（http：//peppergenome.snu.ac.kr/download.php）進行比對，得出比對基因組情況。將測序數據與miRBase中數據進行比對分析得出已知miRNA情況。利用miREvo和mirdeep2軟件預測樣品中的novel miRNA。

1.2.3" miRNA差異表達分析" 對所有miRNA用TPM進行歸一化處理，以P≤0.05和|log₂（FC）|≥3為標準進行顯著差異表達miRNA的篩選工作。

1.2.4" 靶基因預測及功能分析" 使用Patmatch_ v 1.2軟件進行miRNA和靶基因的互補配對分析，再通過程序篩選預測得到最終結果，同時對miRNA靶基因進行GO分析和KEGG富集分析。將差異表達蛋白向GO數據庫（http：//www.gene-ontology.org/）的各term映射，并計算每個term的蛋白數，從而得到具有某個GO功能的蛋白列表及蛋白數目。GO功能分析：應用超幾何檢驗，找出與整個背景蛋白相比，在差異表達蛋白中顯著富集的GO條目，將計算得到的P值通過校正之后，以corrected-P≤0.05為閾值，滿足此條件的GO term定義為在差異表達蛋白中顯著富集的GO term。KEGG富集分析：應用超幾何檢驗，找出與背景蛋白相比，在差異表達蛋白中顯著性富集的Pathway，然后經過多重檢驗校正后，選擇Q≤0.05的Pathway定義為在差異表達蛋白中顯著富集的Pathway。

1.2.5" miRNA及其靶基因的qRT-PCR驗證" 選取6個顯著差異miRNA及其靶基因進行qRT-PCR驗證，引物序列如表1所示。qRT-PCR反應體系為PCR Master Mix 10"μL，上/下游引物（10"μmol/L）各0.4 μL，ddH₂O 8.6 μL，cDNA模板1 μL。擴增條件為：95"℃ 3 min；95"℃ 5 s，60"℃ 30 s，45個循環。內參為U6和Actin。采用2^–ΔΔCt法計算目的基因的相對表達量。使用Microsoft Excel 2016軟件對相對表達量數據進行統計和作圖。

2" 結果與分析

2.1" sRNA高通量測序數據分析

分別構建抗病小米辣種質CF25接種象耳豆根結線蟲第0天（Mi-R0-1、Mi-R0-2、Mi-R0-3）

和第4天（Mi-R4-1、Mi-R4-2、Mi-R4-3）的sRNA文庫，進行高通量測序。由圖1可知，分別從Mi-R0和Mi-R4中獲得39 810 053和34 003 912條原始序列。去污處理后分別得到33"851"394和29"645"658條純凈序列，分別占序列總數的85%和87%。所有樣品的sRNA總豐度為63 497 052條，其中包括rRNA 10 351 257條（16.3%）、SnRNA 277 536條（0.44%）、snoRNA 87 954條（0.14%）、tRNA 499 442條（0.79%）、exon sense 6 709 635條（10.57%）、已知miRNAs 9 149 130條（14.41%）、新發現miRNAs 384 833條（0.61%）、genome others 19 396 665條（30.55%）、未注釋序列16 640 600條（26.21%）（圖2A）。將用于后續分析的clean reads和中國辣椒（Capsicum chinense）參考基因組數據進行比對，結果如圖2B所示。總比對率為73.6%，Mi-R0的平均比對率為78.12%，Mi-R4的平均比對率為66.41%，6個文庫的比對率均達到65%以上。綜上，本研究測得的sRNA序列質量較好，適合開展后續miRNA研究工作。

2.2 "miRNA鑒定

對sRNA測序數據進行鑒定，共檢測到1023個miRNAs，其中包括已知miRNAs 362個和新發現的miRNAs 661個。對這些miRNA進行長度分析，發現其長度主要集中在20～24 nt之間，其中21 nt和24 nt的miRNAs數量最多，分別為233個和502個（圖3），符合植物miRNA的長度特點。對miRNA的堿基偏好性進行分析，發現在第1、3、14、15、18、19、24位點更偏好于尿嘧啶（U），均達到了50%以上；在第2、4、8、9、13、16、20、21、22、23位點更偏好于胞嘧啶（C），均達到51%以上；在第5、6、11、12位點更偏好于鳥嘌呤（G），均達到65%以上；在第7、10、17位點更偏好于腺嘌呤（A），均達到66%以上（圖4A）。miRNA的首位堿基偏好是尿嘧啶（U），在

長度18～22 nt的miRNA中占有比例均超過80%。在24"nt長度的miRNA中5￠端首位堿基腺嘌呤 （A）含量較高，占比為60.60%，其次為尿嘧啶（U），占比為2.50%，最后是鳥嘌呤（G），占比為36.86%（圖4B）。

2.3 "接種前后差異表達miRNA分析

為檢測處理組（Mi-R4）和對照組（Mi-R0）之間以及組內生物學重復間變異度大小，對miRNA測序數據進行主成分（PCA）分析，得到6個主成分。如圖5A所示，主成分PC1和PC2能解釋miRNA轉錄組數據55%的差異。處理組和對照組樣本表現出明顯的分離趨勢，且距離較遠，說明處理組和對照組miRNA的差異比較大。而組內各生物學重復間的miRNA數據表現出一定的聚集性，說明組內樣品間重復性相對較好。采用Pearson相關性系數對6個樣本的相關性進行分析，相關系數越接近1，代表樣本之間相關性越高。結果顯示，組內生物學重復間的相關性均為1，說明組內樣本均一度較高。而組間樣本的相關性處于0.98～0.99之間，表明存在一定差異（圖5B）。綜上所述，處理組和對照組miRNA相對差異較大，而組內差異較小，可進行后續差異表達miRNA的鑒定工作。

如圖6A所示，處理組與對照組間鑒定出133個差異表達miRNAs，其中上調表達的為60個，下調表達的為73個。由此可見，處理組下調表達miRNA比上調表達miRNA數量更多。以P≤0.05和|log₂（FC）|≥3為篩選標準，鑒定出33個顯著差異表達的miRNAs，其中包括19個已知miRNAs和14個新發現miRNAs（圖6B）。表達差異顯著的19個已知miRNAs中，miR827-x、miR827-y、miR5530-z、miR8010-z、miR10518-z、miR6023-z、miR399-y、miR5261-z、miR399-x、miR399-z、miR5185-y等11個miRNAs為上調表達，而miR5658-z、miR1448-z、miR8028-x、miR477-x、miR169-y、miR7712-y、miR477-y、miR395等8個miRNAs則表現出下調表達趨勢（圖6C）。而在新發現的14個顯著差異表達miRNAs中，上調表達的miRNAs有5個，分別為novel-m0242-5p、novel-m0617-3p、novel-m0667-5p、novel-m0307-5p和novel-m0253-5p；顯著下調表達的miRNAs有9個，分別為novel-m0465-5p、novel-m0280-5p、novel-m0548-5p、novel-m0279-3p、novel-m0238- 5p、novel-m0128-3p、novel-m0610-3p、novel- m0395-5p和novel-m0033-3p（圖6D）。

2.4" 差異表達miRNA靶基因預測

對133個差異表達miRNAs進行靶基因預測，得到1781個靶基因，其數量大約是差異表達miRNA數量的13倍。對33個顯著差異表達的miRNAs進行靶基因預測，得到373個靶基因，其中在19個已知miRNAs中預測到296個靶基因。miR5658-z預測靶基因數量最多，為83個。其次是miR5185-y，預測到38個靶基因。miR8028-x預測到37個靶基因，排名第三。而miR5530-z僅預測到了個靶基因。miR8010-z和miR399-x沒有預測到靶基因。新發現的14個miRNAs中共預測到77個靶基因。novel-m0238-5p預測的靶基因最多，為33個。其次為novel-m0617-3p，預測到17個靶基因。novel-m0253-5p預測到13個靶基因，排名第三。而novel-m0667-5p和novel-m0610- 3p分別只預測到1個靶基因（表2）。novel- m0465-3p、novel-m0280-5p、novel-m0548-5p、novel-m0279-3p、novel-m0395-5p、novel-m0033-3p這6個新發現miRNAs沒有預測到靶基因（表2）。

2.5" 差異表達miRNA靶基因GO富集和KEGG代謝通路分析

進一步對差異顯著miRNA的靶基因進行GO富集和KEGG代謝通路分析。將靶基因的GO注釋按照生物學過程（biological process，BP）、細胞組分（cell component，CC）和分子功能（molecular function，MF）分為3類，并將這3類的條目按照注釋到的差異基因數目由大到小進行排序和繪圖。BP中富集的條目最多，主要集中于細胞過程和代謝過程；其次為MF，主要集中于結合和催化活性；富集條目最少的為CC，主要集中于細胞解剖實體和蛋白質復合物等（圖7）。富集條目最多的類目依次為結合、細胞過程、代謝過程和催化活性。

對差異顯著miRNA的靶基因進行KEGG通路富集分析，發現多個代謝通路產生了較大差異。其中，靶基因主要富集于植物-病原體互作途徑（plant-pathogen interaction），其次為MAPK植物信號轉導途徑（MAPK signaling pathway-plant）、嘌呤代謝（purine metabolism）和單環類抗生素生物合成（monobactam biosynthesis）等通路（圖8）。富集于植物-病原體互作的靶基因最多，達到了27個。靶基因數量占第二位是MAPK植物信號轉導途徑，富集了11個靶基因。對靶基因富集的KEGG通路進行分類，可以得到不同分類層級的通路。這33個顯著差異表達miRNAs通過68個靶基因在58條生命途徑中起作用，其中環境適應（environmental adaptation）和信號轉導（signal transduction）這2條生命途徑中富集的靶基因最多（圖9）。由此可見，顯著差異miRNA的靶基因主要富集于植物防御反應和信號轉導等KEGG通路或層級，推測CF25抗象耳豆根結線蟲的分子機理可能主要基于miRNA介導的植物防御反應及其信號轉導相關基因的表達調控。

2.6 "差異表達miRNA及其靶基因的表達驗證

為驗證差異表達miRNA及其靶基因間的調控關系，利用qRT-PCR對6個差異表達miRNAs（miR1448-z、miR8028-x、miR5658-z、novel- m0128-3p、miR5185-y、novel-m0253-5p）及其對應靶基因（CC.CCv1.2.scaffold436.45、CC.CCv1. 2.scaffold876.5、CC.CCv1.2.scaffold703.8、CC. CCv1.2.scaffold785.9、CC.CCv1.2.scaffold19.56）進行表達驗證。如圖10所示，與對照組相比，miR1448-z、miR8028-x、miR5658-z、novel-m0128- 3p在接種象耳豆根結線蟲后呈現下調表達趨勢，

其靶基因CC.CCv1.2.scaffold436.45（RPS2）、CC. CCv1.2.scaffold876.5（ATAS3）、CC.CCv1.2.scaffold 703.8（CCD7）呈現上調表達趨勢，符合miRNA負調控靶基因的特點。而miR5185-y、novel- m0253-5p則是呈現上調表達趨勢，其靶基因CC. CCv1.2.scaffold785.9（ADK）、CC.CCv1.2. scaffold19.56（FLS2）呈現下調表達趨勢。由此可見，所選擇的6個miRNAs以及5個靶基因的qRT-PCR表達驗證結果與基于RNA-seq的基因表達量分析結果基本一致，說明測序結果準確可靠，且這6個miRNAs可能與相關靶基因存在一定的調控關系。

3" 討論

本研究通過高通量測序對抗病小米辣種質CF25接種象耳豆根結線蟲前后的miRNA表達差異進行了探索，篩選出33個表達差異顯著的miRNAs。其中一些miRNA可能與抗病性有一定聯系。miR5658在接種象耳豆根結線蟲4 d后呈現顯著下調表達趨勢。而miR5658在黃瓜綠斑駁花葉病毒（Cucumber green mottle mosaic virus，

CGMMV）侵染黃瓜1dpi時也出現了下調表達^[21]。在毛果楊中miR1448被認為是與抗病相關的miRNA^[22]，在本研究中miR1448呈現下調表達趨勢。miR395可通過調節硫酸鹽代謝來提高水稻對病原體的免疫能力^[23]，在本研究中miR395呈現下調表達趨勢。miR399在抗性品種四季橘（Citrus microcarpa）感染柑橘潰瘍病菌（Xanthomonas citri subsp.citri）5 d時，表達量顯著增加。在柑橘中過表達miR399可使潰瘍病病斑面積顯著減小，病情指數顯著降低，這表明miR399與柑橘的抗潰瘍病密切相關^[24]，在本研究中miR399在抗性辣椒種質CF25接種根結線蟲4 d時也呈現上調表達。過表達miR827可以提高煙草對辣椒脈斑駁病毒（ChiVMV）感染的抵抗力^[25]，本研究中miR827呈現上調表達。還有一些本研究中差異顯著的miRNA在非生物脅迫研究中有所報道，如miR169能響應干旱、低溫、高鹽和高溫^[26]。miR 10518-z、miR5261-z、miR5530-z、miR7712-y、miR8010-z和新發現的novel-m0242-5p、novel- m0617-3p雖未見相關功能的研究報道，但是在小米辣被象耳豆根結線蟲侵染后出現特異性表達，|log₂（FC）|值分別高達7.73、7.64、8.00、8.71、7.75、9.56、8.19，猜測它們在辣椒抗根結線蟲這一過程中可能發揮了重要的作用，值得進一步開展研究。

為了探究miRNA在辣椒抗象耳豆根結線蟲過程中的機制和作用，本研究對其靶基因進行了預測和功能分析。miR5658是33個顯著差異表達miRNA中注釋靶基因最多的miRNA。其靶向的可能與抗病性有關的基因為MYB20、MYB60、MYB306等。MYB家族是植物中最大的轉錄因子家族之一，其成員通過N端保守的MYB結構域與下游基因啟動子區的特異性位點相結合，進而調控這些基因的轉錄^[27-29]，在植物的生長發育、激素信號傳遞以及逆境脅迫響應中發揮關鍵作用^[30-31]。在潰瘍病病菌侵染時，MYB20在抗性品種四季橘（Citrus microcarpa）中表現出顯著的差異表達，感病后其表達量上調了3倍，參與了柑橘對潰瘍病的抗性調控^[32]。此外，MYB20是木質素生物合成途徑的調節因子，在次生細胞壁形成過程中調控木質素的生物合成^[33]。木質素作為植物細胞壁的關鍵成分之一，能夠在植物抵御病原物入侵的過程中起到屏障作用^[34]。MYB60、MYB306雖然在植物抗病方面未見報道，但在非生物脅迫中報道較多，如MYB60參與了細胞中活性氧的平衡調節過程，并以此影響氣孔開閉，響應木薯^[35]和海島棉^[36]的干旱脅迫。而MYB306在燕子花中對低溫脅迫產生了較大的響應^[37]。新預測的14個顯著差異表達的miRNA中，novel-m0238-5p是預測靶基因最多的miRNA，可靶向ART2基因。ART2是一種C₂H₂型鋅指轉錄因子，能在水稻中通過調控細胞壁相關蛋白Os01g0827300上調表達來增強細胞壁的木質化，從而提高植物鋁耐受性^[38]。miRNA可通過靶向抗病基因，讓植物產生抗病反應。本研究預測的靶基因主要富集于植物-病原菌互作、信號轉導等通路。而相關的靶基因主要為FLS2、RPM1、RPS2等。

基于前期轉錄組測序數據，得出CF25接種第0天和第4天的顯著差異表達基因主要富集于代謝途徑、次生代謝物的生物合成等途徑^[20]。結合本研究miRNA差異表達數據進行聯合分析，初步鑒定出19個顯著差異表達miRNA。依據基因功能，可將這些miRNA關聯的轉錄組靶基因大體分為五大類。第一類為代謝相關基因，包括與色素次生代謝相關的PKS1、氮代謝相關的GGT1、嘌呤代謝相關的ADK、色氨酸代謝相關的ASAT3和脂質代謝相關的AAE2等；第二類為生物合成相關基因，包括淀粉生物合成酶APS1、海藻糖合成酶TPS11和獨角金內酯合成酶CCD7等。第三類為生長發育相關基因，如花發育相關基因AGL27、生物鐘調節基因LHY和DNA復制相關基因RFC3等；第四類為與應激響應相關的基因，如SBT4.14和SIN2等；第五類為植物抗病基因，如TGA1、RPP13、FLS2和RPS2等。由此可見，FLS2和RPS2可能是CF25在miRNA調控下防御象耳豆根結線蟲的重要基因。novel-m0128- 5p在接種后表達量降低了1.37倍，其潛在靶基因是FLS2。鞭毛感知蛋白FLS2（flagellin sensing 2）是屬于XII家族的富亮氨酸重復受體樣激酶，FLS2可通過識別鞭毛蛋白表位flg22，對病菌入侵做出迅速反應^[39]。FLS2是植物免疫系統的重要組成蛋白^[40]，植物體內缺失FLS2蛋白會導致植物對病菌的免疫性下降^[41]。FLS2是蒜芥茄（Solanum sisymbriifolium Lam.）響應黃萎病（eggplant verticillium wilt）脅迫的重要基因^[42]，也是中華獼猴桃抗潰瘍病（Pseudomonas syringae pv. actinidiae， Psa）的候選抗病基因^[43]。在本研究中，基于前期轉錄組數據可以看出，抗病小米辣CF25在接種象耳豆根結線蟲第4天后，FLS2上調表達，推測novel-m0128-5p通過負調控靶基因FLS2增強了CF25的免疫系統，從而提高其抗病能力。目前抗病R蛋白大多數屬于NBS-LRR類^[44]，NB-ARC屬于NBS-LRR類蛋白中的一個分支，可以通過結合ADP或ATP來參與植物抗病蛋白的構象變化，進而激活植物免疫響應與病原菌識別^[45]。RPS2是典型的NB-ARC類抗病基因，在擬南芥抗丁香假單胞菌（Pseudomonas syringae）中發揮重要作用^[46]，此外，RPS2在水稻中過表達后，對真菌病原體稻瘟病菌和細菌病原體水稻白葉枯病菌產生正響應^[47]。在本研究中，RPS2關聯到了miR6023-z、miR1448和miR8028三個miRNAs，是關聯miRNA數量最多的靶基因。其中，miR1448與其靶基因RPS2^[48]都被認為與抗病有關，在本研究中，miR1448可能通過下調表達從而引起RPS2上調表達發揮抗病作用。

綜上所述，本研究通過分析抗象耳豆根結線蟲小米辣資源CF25在接種前后的sRNA測序數據，得出miRNA可能調控了植物抗病及信號轉導方面的基因表達變化，并由此介導了CF25抗病反應的發生。本研究獲得了較多差異表達miRNA，但它們具體如何通過調節靶基因在根結線蟲和寄主的互作中發揮抗性作用，需要進一步通過試驗來驗證。還可在此基礎上深入挖掘在辣椒與根結線蟲互作機制中發揮重要調節作用的miRNA，為后續輔助抗病育種工作奠定基礎。

4" 結論

本研究對抗性小米辣材料CF25接種象耳豆根結線蟲前后的sRNA數據進行分析，鑒定出33個表達差異顯著的miRNA。其中包括一些已知與作物抗病性相關的miRNA，如miR5658、miR1448、miR395、miR399和miR827等。對這33個表達差異顯著的miRNA進行靶基因預測，得到373個靶基因。經GO功能和KEGG代謝通路富集分析，得出靶基因主要富集于植物-病原體互作以及信號轉導等代謝途徑。推測miRNA可能通過調控靶基因參與這些代謝通路發揮抗病功能。本研究為深入了解辣椒中miRNA介導的根結線蟲抗性機制提供有益的參考數據。

參考文獻

<！--[if ！supportLists]-->[1]"""""" <！--[endif]-->鄒學校， 朱凡. 辣椒的起源、進化與栽培歷史[J]. 園藝學報， 2022， 49（6）： 1371-1381.

ZOU X X， ZHU F. Origin， evolution and cultivation history of the pepper[J]. Acta Horticulturae Sinica， 2022， 49（6）： 1371-1381. （in Chinese）

<！--[if ！supportLists]-->[2]"""""" <！--[endif]-->梁文薈， 單慶云， 王瑾， 王兵， 張美玲， 萬彧， 鄒學校， 熊程. 辣椒果實辣椒素類物質生物合成調控及其影響因素研究進展[J]. 中國蔬菜， （2025-01-08）[2025-01-20].

https：// doi.org/10.19928/j.cnki.1000-6346.2025.2005.

LIANG W H， SHAN Q Y， WANG J， WANG B， ZHANG M L， WAN Y， ZOU X X， XIONG C. Research progress on biosynthesis regulation and influencing factors of capsaicinoids in pepper fruit[J]. China Vegetables， （2025-01-08）[2025-01- 20]."https：//doi.org/10.19928/j.cnki.1000-6346.2025.2005."（in

Chinese）

<！--[if ！supportLists]-->[3]"""""" <！--[endif]-->MAGA J A， TODD P H. Critical reviews in food Science amp; Nutrition[J]. Capsicum， 1975， 6（2）： 177-199.

<！--[if ！supportLists]-->[4]"""""" <！--[endif]-->KIEWNICK S， DESSIMOZ M， FRANCK L. Effects of the Mi-1 and the N root-knot nematode resistance gene on infection and reproduction of Meloidogyne enterolobii on tomato and pepper cultivars[J]. Journal of Nematology， 2009， 41： 134-139.

<！--[if ！supportLists]-->[5]"""""" <！--[endif]-->PINHEIRO J B， BOITEUX L S， ALMEIDA M R A， PEREIRA R B， GALHARDO L C S， CARNEIRO R M D G. First report of Meloidogyne enterolobii in capsicum root stocks carrying the Me1 and Me3/Me7 genes in central Brazil[J]. Nematropica， 2015， 45（2）： 184-188.

<！--[if ！supportLists]-->[6]"""""" <！--[endif]-->王劍， 宋志強， 成飛雪， 程菊娥， 張德詠， 劉勇. 湖南省辣椒上首次發現象耳豆根結線蟲[J]. 植物保護， 2015， 41（4）： 180-183.

WANG J， SONG Z Q， CHENG F X， CHENG J E， ZHANG D Y， LIU Y. First report of Meloidogyne enterolobii on pepper in Hunan province[J]. Plant Protection， 2015， 41（4）： 180-183. （in Chinese）

<！--[if ！supportLists]-->[7]"""""" <！--[endif]-->HWANG D G， PARK J H， LIM J Y， KIM D， CHOI Y， KIM S， REEVES G， YEOM S I， LEE J S， PARK M， KIM S， CHOI I Y， CHOI D， SHIN C. The hot pepper （Capsicum annuum） microRNA transcriptome reveals novel and conserved targets： a foundation for understanding MicroRNA functional roles in hot pepper[J]. PLoS One， 2013， 8（5）： e64238.

<！--[if ！supportLists]-->[8]"""""" <！--[endif]-->YANG L， MU X， LIU C， CAI J， SHI K， ZHU W， YANG Q. Overexpression of potato miR482e enhanced plant sensitivity to Verticillium dahliae infection[J]. Journal of Integrative Plant Biology， 2015， 57（12）： 1078-1088.

<！--[if ！supportLists]-->[9]"""""" <！--[endif]-->LI Y， CAO X， ZHU Y， YANG X， ZHANG K， XIAO Z， WANG H， ZHAO J， ZHANG L， LI G. Osa-miR398b boosts H₂O₂ production and rice blast disease-resistance via multiple superoxide dismutases[J]. New Phytologist， 2019， 222（3）： 1507-1522.

<！--[if ！supportLists]-->[10]""" <！--[endif]-->LI Y， WANG N， GUO J， ZHOU X J， BAI X Y， Azeem M， ZHU L Y， CHEN L， CHU M， WANG H， CHENG W. integrative transcriptome analysis of mRNA and miRNA in pepper’s response to Phytophthora capsici infection[J]. Biology， 2024， 13（3）： 186.

<！--[if ！supportLists]-->[11]""" <！--[endif]-->PARK J H， SHIN C. The role of plant small RNAs in NB-LRR regulation[J]. Briefings in Functional Genomics. 2015， 14（4）： 268-274.

<！--[if ！supportLists]-->[12]""" <！--[endif]-->楊德榮. 大有開發前途的灌木辣椒[J]. 北京農業， 1994， 9： 35.

YANG D R. The shrub chili pepper， which holds great potential for development[J]. Beijing Agriculture， 1994， 9： 35. （in Chinese）

<！--[if ！supportLists]-->[13]""" <！--[endif]-->鄧明華， 文錦芬. 云南灌木辣椒資源[J]. 辣椒雜志， 2009， 7（1）： 36-37.

DENG M H， WEN J F. Yunnan shrub chili resources[J]. Journal of China Capsicum， 2009， 7（1）： 36-37.

<！--[if ！supportLists]-->[14]""" <！--[endif]-->LIU S， LI W， WU Y， CHEN C， LEI J， CHRISTIAN S. De novo transcriptome assembly in chili pepper （Capsicum fiutescens） to identify genes involved in the biosynthesis of capsaicinoids[J]. PLoS One， 2013， 8（1）： e48156.

<！--[if ！supportLists]-->[15]""" <！--[endif]-->KIM J S， JEE H J， GWAG J G， KIM C K， SHIM C K. Evaluation on red pepper germplasm lines （Capsicum spp.） for resistance to anthracnose caused by Colletotrichum acutaum[J]. The Plant Pathology Journal， 2010， 26（3）： 273-279.

<！--[if ！supportLists]-->[16]""" <！--[endif]-->方榮， 周坤華， 馬輝剛， 陳學軍， 何烈干， 繆南生. 中國灌木辣椒種質農藝性狀鑒定與疫病抗性評價[J]. 植物遺傳資源學報， 2014， 15（1）： 186-191.

FANG R， ZHOU K H， MA H G， CHEN X J， HE L G， MIU N S. Identification of agronomic traits and evaluation of resistance to phytophthora capsici in C. frutescens in China[J]. Journal of Plant Genetic Resources， 2014， 15（1）： 186-191. （in Chinese）

<！--[if ！supportLists]-->[17]""" <！--[endif]-->LIU S， LI W， WU Y， CHEN C， LEI J. De novo transcriptome assembly in chili pepper （Capsicum frutescens） to identify genes involved in the biosynthesis of capsaicinoids[J]. PLoS One， 2013， 8（1）： e48156.

<！--[if ！supportLists]-->[18]""" <！--[endif]-->莊燦然， 呂金殿， 梁耀琦. 中國干制辣椒[M]. 北京： 中國農業出版社， 1995.

ZHUANG C R， LYU J D， LIANG Y Q. Chinese chili[M]. Beijing： China Agricultural Press， 1995. （in Chinese）

<！--[if ！supportLists]-->[19]""" <！--[endif]-->陳梅紅， 黃銀香， 林瓏， 高暢， 涂容生， 朱婕， 劉子記. 灌木辣椒（Capsicum frutescens）抗象耳豆根結線蟲的組織學及生理生化特征分析[J]. 熱帶作物學報， 2024， 45（9）： 1902-1911.

CHEN M H， HUANG Y X， LIN L， GAO C， TU R S， ZHU J， LIU Z J. Analysis of histological and physiological-biochemical characteristics of Capsicum frutescens on the resistance to Meloidogyne enterolobii[J]. Chinese Journal of Tropical Crops， 2024， 45（9）： 1902-1911. （in Chinese）

<！--[if ！supportLists]-->[20]""" <！--[endif]-->陳梅紅. 抗象耳豆根結線蟲海南本地小米辣資源的篩選及其作用機理研究[D]. 海口： 海南大學，

2024.

CHEN M H. Screening of Hainan endemic Xiaomila resources for resistance to Meloidogyne enterolobii and the study of mechanisms of resistance[D]. Haikou： Hainan University， 2024. （in Chinese）

<！--[if ！supportLists]-->[21]""" <！--[endif]-->LIANG C Q， LIU H W， HAO J J， LI J Q， LUO L X. Expression profiling and regulatory network of cucumber microRNAs and their putative target genes in response to cucumber green mottle mosaic virus infection[J]. Archives of Virology， 2019， 164 （4）： 1121-1134.

<！--[if ！supportLists]-->[22]""" <！--[endif]-->ZHAO J P， DIAO S， ZHANG B Y， NIU B Q， WANG Q L， WAN X C， LUO Y Q. Phylogenetic analysis and molecular evolution patterns in the MIR482-MIR1448 polycistron of Populus L.[J]. PloS One， 2012， 7（10）： e47811.

<！--[if ！supportLists]-->[23]""" <！--[endif]-->YANG Z Y， HUI S G， LV Y， ZHANG M J， CHEN D， TIAN J J， ZHANG H T， LIU H B， CAO J B， XIE W Y， WU C Y， WANG S P， YUAN M. miR395-regulated sulfate metabolism exploits pathogen sensitivity to sulfate to boost immunity in rice[J]. Molecular plant， 2021， 15（4）： 671-688.

<！--[if ！supportLists]-->[24]""" <！--[endif]-->王兆昊， 郭興茹， 張樂歡， 何永睿， 陳善春， 姚利曉. csi-miR399響應柑橘潰瘍病菌侵染的表達模式及其抗病性分析[J]. 中國農業科學， 2023， 56（8）： 1484-1493.

WANG Z H， GUO X R， ZHANG L H， HE Y R， CHEN S C， YAO L X. Expression pattern of csi-miR399 in response to Xanthomonas citri subsp. citri Infection and Its disease resistance analysis[J]. Scientia Agricultura Sinica， 2023， 56（8）： 1484-1493. （in Chinese）

<！--[if ！supportLists]-->[25]""" <！--[endif]-->TANG R X， YANG Y F， JI C L， SU Y S， JIAO B L，YUAN B W， YANG X Y， XI D H. MiR827 positively regulates the resistance to chilli veinal mottle virus by affecting the expression of FBPase in Nicotiana benthamiana[J]. Physiologia Plantarum， 2024， 176（3）： e14375.

<！--[if ！supportLists]-->[26]""" <！--[endif]-->黃渺， 宋澤洋， 蔡駿一， 南紅， 肇瑾. 辣椒miR169家族在非生物脅迫下的表達及啟動子順式作用元件分析[J]. 西部林業科學， 2023， 52（1）： 165-174.

HUANG M， SONG Z Y， CAI J Y， NAN H， ZHAO J. Analysis on expression under abiotic stress and the promoter cis-acting elements of miR169 family in pepper （Capsicum annuum L.）[J]. Journal of West China Forestry Science， 2023， 52（1）： 165-174. （in Chinese）

<！--[if ！supportLists]-->[27]""" <！--[endif]-->AMBAWAT S， SHARMA P， YADAV N R， YADAV R C. MYB transcription factor genes as regulators for plant responses： an overview[J]. Physiology and Molecular Biology of Plants， 2013， 19（3）： 307-321.

<！--[if ！supportLists]-->[28]""" <！--[endif]-->DUBOS C， STRACKE R， GROTEWOLD E， WEISSHAAR B， MARTIN C， LEPINIEC L. MYB transcription factors in Arabidopsis[J]. Trends in Plant Science， 2010， 15（10）： 573- 581.

<！--[if ！supportLists]-->[29]""" <！--[endif]-->BALDONI E， GENGA A， COMINELLI E. Plant MYB Transcription factors： their role in drought response mechanisms[J]. International Journal of Molecular Sciences， 2015， 16（7）： 15811-15851.

<！--[if ！supportLists]-->[30]""" <！--[endif]-->PIREY R E M， BUROW M. Regulation of MYB and bHLH transcription factors： a glance at the protein level[J]. Molecular Plant， 2015， 8： 378-388.

<！--[if ！supportLists]-->[31]""" <！--[endif]-->GOLLDACK D， LI C， MOHAN H， PROBST N. Tolerance to drought and salt stress in plants： unraveling the signaling networks[J]. Frontiers in Plant Science， 2014， 5： 151.

<！--[if ！supportLists]-->[32]""" <！--[endif]-->范海芳， 彭蘊， 張慶雯， 何永睿， 鄒修平， 李強， 彭愛紅， 許蘭珍， 雷天剛， 陳善春， 姚利曉. 柑橘R2R3-MYB類轉錄因子CitMYB20的克隆與表達分析[J]. 分子植物育種， 2019， 17（10）： 3199-3207.

FAN H F， PENG Y， ZHANG Q W， HE Y R， ZOU X P， LI Q， PENG A H， XU L Z， LEI T G， CHEN S C， YAO L X. Cloning and expression of citrus R2R3-MYB transcription factor CitMYB20[J]. Molecular Plant Breeding， 2019， 17（10）： 3199-3207. （in Chinese）

<！--[if ！supportLists]-->[33]""" <！--[endif]-->GENG P， ZHANG S， LIU J Y， ZHAO C H， WU J， CAO Y P， FU C X， HAN X， HE H， ZHAO Q. MYB20， MYB42， MYB43， and MYB85 regulate phenylalanine and lignin biosynthesis during secondary cell wall formation[J]. Plant Physiology， 2020， 182（3）： 1272-1283.

<！--[if ！supportLists]-->[34]""" <！--[endif]-->賓金華， 潘瑞熾. 茉莉酸甲酯誘導煙草幼苗抗病與過氧化物酶活性和木質素含量的關系[J]. 應用與環境生物學報， 1999（2）： 49-53.

BIN J H， PAN R Z. The relationship of the disease resistance of tobacco seedlings induced by methyl jasmonate with peroxidase activity and lignin content[J]. Chinese Journal of Applied and Environmental Biology， 1999（2）： 49-53. （in Chinese）

<！--[if ！supportLists]-->[35]""" <！--[endif]-->徐子寅， 于曉玲， 鄒良平， 趙平娟， 李文彬， 耿夢婷， 阮孟斌. 木薯MYB轉錄因子基因MeMYB60表達特征分析及其互作蛋白篩選[J]. 作物學報， 2023， 49（4）： 955-965.

XU Z Y， YU X L， ZOU L P， ZHAO P J， LI W B， GENG M T， RUAN M B. Expression pattern analysis and interaction protein screening of cassava MYB transcription factor MeMYB60[J]. Acta Agronomica Sinica， 2023， 49（4）： 955- 965. （in Chinese）

<！--[if ！supportLists]-->[36]""" <！--[endif]-->高巍， 劉會利， 田新權， 張慧， 宋潔， 楊勇， 龍璐， 宋純鵬. 海島棉轉錄因子基因GbMYB60的克隆、表達及其抗逆性分析[J]. 作物學報， 2016， 42（9）： 1342-1351.

GAO W， LIU H L， TIAN X Q， ZHANG H， SONG J， YANG Y， LONG L， SONG C P. Cloning， expression， and functional analysis of transcription factor gene GbMYB60 in cotton[J]. Acta Agronomica Sinica， 2016， 42（9）： 1342-1351. （in Chinese）

<！--[if ！supportLists]-->[37]""" <！--[endif]-->楊娟. 燕子花IlMYB306基因克隆及其抗逆功能初探[D]. 哈爾濱： 東北林業大學， 2019.

YANG J. Cloning and primary functional analysis on stress resistance of IlMYB306 gene in Iris laevigata[D]. Harbin： Northeast Forestry University， 2019. （in Chinese）

<！--[if ！supportLists]-->[38]""" <！--[endif]-->CHE J， TSUTSUI T， YOKOSHO K， YAMAJI N， MA J F. Functional characterization of an aluminum （Al）-inducible transcription factor， ART2， revealed a different pathway for Al tolerance in rice[J]. New Phytologist， 2018， 220（1）： 209- 218.

<！--[if ！supportLists]-->[39]""" <！--[endif]-->SUN Y D， LI L， MACHO A P， HAN Z F， HU Z H， ZIPFEL C， ZHOU J M， CHAI J J. Structural basis for flg22-induced activation of the Arabidopsis FLS2-BAK1 immune complex[J]. Science， 2013， 342（6158）： 624-628.

<！--[if ！supportLists]-->[40]""" <！--[endif]-->TRDá L， BOUTROT F， CLAVERIE J， BRULé D， DOREY S， POINSSOT B. Perception of pathogenic or beneficial bacteria and their evasion of host immunity： pattern recognition receptors in the frontline[J]. Frontiers in Plant Science， 2015， 6： 219.

<！--[if ！supportLists]-->[41]""" <！--[endif]-->劉沁松. 植物免疫受體FLS2的胞內運輸機制[J]. 生命科學， 2021， 33（5）： 576-581.

LIU Q S. Intracellular trafficking of the plant immune receptor FLS2[J]. Chinese Bulletin of Life Sciences， 2021， 33（5）： 576-581. （in Chinese）

<！--[if ！supportLists]-->[42]""" <！--[endif]-->孫茂， 吳麗艷， 龔亞菊， 鮑銳， 桂敏， 黎志彬， 杜光輝. 蒜芥茄SsFLS2基因的克隆及其表達分析[J]. 西南農業學報， 2024， 37（8）： 1669-1676.

SUN M， WU L Y， GONG Y J， BAO R， GUI M， LI Z B， DU G H. Cloning and expression analysis of SsFLS2 gene in Solanum sisymbriifolium Lam[J]. Southwest China Journal of Agricultural Sciences， 2024， 37（8）： 1669-1676. （in Chinese）

<！--[if ！supportLists]-->[43]""" <！--[endif]-->匡美美， 李黎， 馬建偉， 劉原， 蔣鴻霏， 雷瑞， 滿玉萍， 王一帆， 黃波， 王彥昌， 劉世彪. 利用中華獼猴桃雜交后代轉錄組測序篩選抗潰瘍病相關基因[J]. 園藝學報， 2024， 51（8）： 1743-1757.

KUANG M M， LI L， MA J W， LIU Y， JIANG H F， LEI R， MAN Y P， WANG Y F， HUANG B， WANG Y C， LIU S B. Screening of kiwifruit canker resistance-related genes based on the transcriptome sequencing analysis of hybrids from actinidia chinensis[J]. Acta Horticulturae Sinica， 2024， 51（8）： 1743-1757. （in Chinese）

<！--[if ！supportLists]-->[44]""" <！--[endif]-->KIM E H， JANE E P. Deciphering plant-pathogen communication： fresh perspectives for molecular resistance breeding[J]. Current Opinion in Biotechnology， 2003， 14（2）： 177- 193.

<！--[if ！supportLists]-->[45]""" <！--[endif]-->LEAH M， TAN X P， PATRICE K， RICHARD W M. Plant NBS-LRR proteins： adaptable guards[J]. Genome Biology， 2006， 7（4）： 212.

<！--[if ！supportLists]-->[46]""" <！--[endif]-->KUNKEL B N， BENT A F， DAHLBECK D， INNES R W， STASKAWISC B J. RPS2， an Arabidopsis disease resistance locus specifying recognition of Pseudomonas syringae strains expressing the avirulence gene avrRpt2[J]. The Plant Cell， 1993， 5（8）： 865-875.

<！--[if ！supportLists]-->[47]""" <！--[endif]-->FLOR H. Current status of the gene for gene concept[J]. Annual Review of Phytopathology， 1971（9）： 275-296.

<！--[if ！supportLists]-->[48]""" <！--[endif]-->GURURANI A M， VENKATESH J， UPADHYAYA P C， NOOKARAJU A， PANDEY S K， PARK S W. Plant disease resistance genes： current status and future directions[J]. Physiological and Molecular Plant Pathology， 2012， 78： 51- 65.