馬凱恒 何佳凝 謝牧洪 任 斌 黃倩雪 楊桂燕*

（1. 西北農(nóng)林科技大學(xué)林學(xué)院陜西省核桃工程技術(shù)研究中心，楊凌 712100；2. 西北農(nóng)林科技大學(xué)林學(xué)院陜西省經(jīng)濟(jì)植物資源開(kāi)發(fā)利用重點(diǎn)實(shí)驗(yàn)室，楊凌 712100）

植物抗病基因在識(shí)別病原菌和抗病信號(hào)傳導(dǎo)的病害防御過(guò)程中起重要作用。近年，利用轉(zhuǎn)座子標(biāo)簽技術(shù)、定位克隆和抗病基因同源序列法等從植物中克隆到大量的抗病基因。這些基因的克隆、鑒定和功能分析對(duì)于深入研究植物抗病機(jī)制、通過(guò)基因工程提高植物抗病性具有十分重要的意義。（nonexpressor of pathogenesis-related genes 1）基因是植物防御信號(hào)轉(zhuǎn)導(dǎo)途徑的重要調(diào)控因子，是植物水楊酸（Salicylic acid，SA）信號(hào)通路的重要組成部分，參與了SA 與茉莉酸（Jasmonic acid，JA）∕乙烯（Ethylene，ET）信號(hào)途徑的相互作用；因此是植物響應(yīng)病原菌侵染的核心因子。典型的基因含有BTB和ANK-r結(jié)構(gòu)域，這與NPR1 同其他轉(zhuǎn)錄因子互作密切相關(guān)。如，與bZIP 轉(zhuǎn)錄因子中的TGA 家族之間的互作，通過(guò)BTB 結(jié)構(gòu)域與TGA2 的抑制區(qū)結(jié)合使TGA2 抑制病程相關(guān)基因（Pathogenesis-related gene，PR）表達(dá)的功能缺失。目前，很多植物中的基因被陸續(xù)克隆，并對(duì)其功能進(jìn)行了研究，發(fā)現(xiàn)來(lái)自不同植物的同源基因或相同植株的不同基因，其功能具有一定差異。如，煙草（）具有對(duì)SA 敏感的反式激活結(jié)構(gòu)域，在水稻（）中過(guò)表達(dá)基因會(huì)自發(fā)激活防御基因的表達(dá)，在水稻和草莓（）中防御反應(yīng)的激活可能會(huì)對(duì)植物生長(zhǎng)發(fā)育產(chǎn)生負(fù)面影響，影響作物生產(chǎn)相關(guān)的重要過(guò)程，但在擬南芥（）和許多其它植物中不存在這種現(xiàn)象。在大量植物中，過(guò)表達(dá)及其同源基因被證實(shí)能夠增強(qiáng)植物的抗病性，說(shuō)明了1在抗病中的保守作用?？梢?jiàn)，對(duì)未被報(bào)道的基因進(jìn)行研究，一方面可以探究基因應(yīng)對(duì)植株病害的功能，另一方面又可能挖掘基因的新功能。

核桃（）屬胡桃科（Juglandaceae）胡桃屬（）多年生落葉喬木，享有“世界四大堅(jiān)果”之一的美譽(yù)，也是我國(guó)主要經(jīng)濟(jì)樹(shù)種之一。隨著核桃栽培面積不斷擴(kuò)大，核桃產(chǎn)業(yè)鏈逐漸延長(zhǎng)，為區(qū)域經(jīng)濟(jì)發(fā)展做出了重要貢獻(xiàn)。但近年，核桃病蟲(chóng)害問(wèn)題日益突出，對(duì)核桃生長(zhǎng)及產(chǎn)量產(chǎn)生了嚴(yán)重的制約，尤其是核桃炭疽病、枝枯病和黑斑病對(duì)核桃植株的危害較大，嚴(yán)重時(shí)可能導(dǎo)致絕收，甚至樹(shù)木死亡。當(dāng)前，防治核桃此類(lèi)病害主要依靠噴灑化學(xué)殺菌劑?；瘜W(xué)農(nóng)藥存在殘留問(wèn)題，對(duì)食品安全、環(huán)境保護(hù)等都會(huì)造成不良影響，這急切需要掌握核桃的病害調(diào)控機(jī)制，進(jìn)而為尋求高效安全的病害防治措施提供依據(jù)。因此，本研究從核桃中鑒定獲得1條基因（命名），通過(guò)基本生物信息及病害響應(yīng)表達(dá)分析其生物學(xué)功能，為揭示核桃抗病響應(yīng)機(jī)制提供候選基因。

1 材料與方法

1.1 材料及處理

植物材料：2年生‘香玲’核桃嫁接苗。

病原：①膠孢炭疽菌（）（簡(jiǎn)寫(xiě)為CgTJ）：引起核桃炭疽病的病原菌。②矩圓黑盤(pán)孢菌（）（簡(jiǎn)寫(xiě)為MoZK）：引起核桃枝枯病的病原菌。③黃單胞菌（）（簡(jiǎn)寫(xiě)為XcHB）：引起核桃細(xì)菌性黑斑病的病原菌。

處理：分別制備濃度為1×10個(gè)·L的膠孢炭疽 菌 分 生 孢 子 菌 液，配 置100 μmol·L水 楊 酸（SA）溶液。對(duì)植株分別進(jìn)行SA、CgTJ、CgTJ+SA、MoZK、MoZK+SA、XcHB、XcHB+SA 噴灑處理。以噴水處理為對(duì)照。在處理后第0、1、6、9、12 d 等時(shí)間分別取葉和莖，用液氮速凍后保存于-80 ℃冰箱備用。每個(gè)處理包含6棵植株。

1.2 JrNPR1基因的克隆與分析

以“nonexpressor of pathogenesis related genes 1”為關(guān)鍵詞在‘香玲’核桃轉(zhuǎn)錄組中篩選基因，通過(guò)BLAST 分析對(duì)比后，選擇其中1 條基因（命名為）進(jìn)行后續(xù)分析。用ORF finder（http：∕∕www.ncbi.nlm.nih.gov∕gorf∕gorf.Html）確 定基因開(kāi)放讀碼框（ORF），通過(guò)ExPASy（http：∕∕web.expasy.org∕protparam∕）分 析基因的序列特征。采用CD-Search（https：∕∕www.ncbi.nlm.nih.gov∕Structure∕cdd∕wrpsb.cgi）對(duì)的保守結(jié)構(gòu)域進(jìn)行分析。利用BLASTP（http：∕∕blast.ncbi.nlm.nih.gov∕Blast）對(duì)同源性蛋白序列進(jìn)行搜索。利用MEGA7 構(gòu)建系統(tǒng)進(jìn)化樹(shù)。使用NEW PLACE（https：∕∕www.dna.affrc.go.jp∕PLACE∕）分析啟動(dòng)子順式作用元件。

1.3 JrNPR1基因的表達(dá)分析

樣品總RNA 采用十六烷基三甲基溴化銨（CTAB）方法提取。RNA 反轉(zhuǎn)錄cDNA 使用Prime-ScriptRT Reagent Kit（CWBIO，康為世紀(jì)，中國(guó)）試劑盒進(jìn)行。cDNA 稀釋10 倍用作實(shí)時(shí)熒光定量PCR（qRT-PCR）模板。qRT-PCR 反應(yīng)體系參照SYBR Green Real time PCR Master mix（CWBIO）制備。定量反應(yīng)儀器為Step OneReal-Time PCR System（Applied Biosystems 生產(chǎn)）。反應(yīng)程序?yàn)椋?4 ℃預(yù)變性30 s；94 ℃變性12 s，60 ℃退火45 s，72 ℃延伸45 s，45個(gè)循環(huán)；81 ℃讀板1 s，重復(fù)3次?；蚨恳餅椋?′-CCGATGTATGTTGGCTGT-3′∕5′-GGTCTACATGCATCATGG-3′。 內(nèi)參基因?yàn)?8S rRNA（HE574850）基因，引物為：5′-GGTCAATCTTCTCGTTCCCTT-3′∕5′-TCGCATTTCGCTACGTTCTT-3′。相對(duì)表達(dá)結(jié)果用2法進(jìn)行分析。數(shù)據(jù)使用SPSS 軟件包（SPSS，Chicago，Illinois，USA）處理。樣品變異性用標(biāo)準(zhǔn)偏差表示。不同時(shí)間點(diǎn)與0 d 之間的表達(dá)差異用檢驗(yàn)分析（＜0.05）。

2 結(jié)果與分析

2.1 JrNPR1全長(zhǎng)cDNA序列基本信息

通過(guò)查找核桃轉(zhuǎn)錄組數(shù)據(jù)獲得1 條基因，經(jīng)分析發(fā)現(xiàn)該基因ORF 長(zhǎng)1 782 bp（GeneBank登錄號(hào)：MW202261），擬推導(dǎo)的蛋白含有593 個(gè)氨基酸（見(jiàn)圖1），理論等電點(diǎn)為6.40。保守結(jié)構(gòu)域分析發(fā)現(xiàn)含有NPR1-like-C 結(jié)構(gòu)域（見(jiàn)圖1），多序列比對(duì)發(fā)現(xiàn)該蛋白與其它物種的NPR1 蛋白相似性很高（見(jiàn)圖2），表明屬于NPR1-like基因家族成員。構(gòu)建進(jìn)化樹(shù)分析發(fā)現(xiàn)JrNPR1 蛋白與楊梅（）、歐洲栓皮櫟（）以及加州白櫟（）等的進(jìn)化關(guān)系較近（見(jiàn)圖3）。

圖1 JrNPR1蛋白的NPR結(jié)構(gòu)域分布Fig.1 The NPR domain distribution of JrNPR1 protein

圖2 JrNPR1蛋白與其同源蛋白的NPR結(jié)構(gòu)域比對(duì)Jr.核桃；Mr.楊梅；Rc.蓖麻；Pd.美洲黑楊；Ss.密花豆Fig.2 Multiple sequence alignment of NPR domains from JrNPR1 and its homologous proteinsJr.Juglans regia；Mr.Morella rubra；Rc.Ricinus communis；Pd.Populus deltoids；Ss.Spatholobus suberectus

圖3 JrNPR1蛋白的進(jìn)化分析Jr.核桃；Mr.楊梅；Rc.蓖麻；Pd.美洲黑楊；Ss.密花豆；Ql.加州白櫟；Qs.歐洲栓皮櫟；Cc.木豆；Ap.非洲相思子；Pe.胡楊；Pa.銀白楊；Pt.毛果楊Fig.3 Phylogenetic tree analysis of JrNPR1 proteinJr.Juglans regia；Mr.Morella rubra；Rc.Ricinus communis；Pd.Populus deltoids；Ss.Spatholobus suberectus；Ql.Quercus lobata；Qs.Quercus suber；Cc.Cajanus cajan；Ap.Abrus precatorius；Pe.Populus euphratica；Pa.Populus alba；Pt.Populus trichocarpa

2.2 JrNPR1啟動(dòng)子元件特征

經(jīng)BLASTX 比對(duì)發(fā)現(xiàn)與GeneBank 中的XP_018838913.1 序列一致，以此鑒定其上游1 233 bp 的DNA 序列設(shè)計(jì)引物，以‘香玲’核桃DNA 為模板擴(kuò)增獲得該基因啟動(dòng)子序列進(jìn)行生物信息學(xué)分析，NEW PLACE 預(yù)測(cè)顯示，啟動(dòng)子中含有多種與抗病性相關(guān)順式作用元件，如WRKY 轉(zhuǎn)錄因子識(shí)別的WRKY71OS 元件，RAV1轉(zhuǎn)錄因子識(shí)別的RAV1AAT元件，組織特異表達(dá)的CAATBOX1 元件（見(jiàn)表1）。推測(cè)其具有調(diào)控參與抗病的能力。

表1 JrNPR1啟動(dòng)子正向包含的主要順式作用元件Table 1 The cis-acting elements contained in JrNPR1 promoter

2.3 JrNPR1基因響應(yīng)不同病害的表達(dá)

對(duì)試材分別進(jìn)行SA、CgTJ、CgTJ+SA、MoZK、MoZK+SA、XcHB、XcHB+SA 等處理，分析處理后第0、1、6、9、12 d 時(shí)葉和莖中基因的表達(dá)（見(jiàn)圖4～5）。結(jié)果發(fā)現(xiàn)，在葉中炭疽病脅迫下，脅迫1 d 時(shí)的表達(dá)達(dá)到最大值，為對(duì)照的9.96 倍，之后的表達(dá)逐漸降低，與SA 處理下的變化趨勢(shì)一致。在進(jìn)行CgTJ+SA處理后，基因隨脅迫后時(shí)間延長(zhǎng)的變化趨勢(shì)與單獨(dú)進(jìn)行SA和CgTJ處理的變化趨勢(shì)一致，但其表達(dá)水平高于單一處理下的表達(dá)，在0～12 d，分別為單一SA和CgTJ 處 理 下 的1.15～3.78、1.31～7.66 倍。在MoZK、MoZK+SA 脅迫下，基因的表達(dá)變化與SA、CgTJ、CgTJ+SA 相似，在1 d 的表達(dá)值最大，分別為對(duì)照的7.00、12.41 倍。XcHB、XcHB+SA 處理下的表達(dá)最大值出現(xiàn)在6 d，為對(duì)照的7.91、17.86 倍。在SA 存在情況下，基因響應(yīng)XcHB 表 達(dá) 是 單 獨(dú)XcHB 處 理 的2.21～178.89 倍（見(jiàn)圖4）。

圖4 不同處理下JrNPR1基因在葉中的表達(dá)*表示相同處理下不同時(shí)間點(diǎn)與0 d之間的差異顯著（P＜0.05），下同F(xiàn)ig.4 The expression of JrNPR1 gene in leaf under different treatments*Indicates that the difference between different time points and 0 d under the same treatment is significant（P＜0.05），the same as below

在莖中，基因在SA、CgTJ、CgTJ+SA、XcHB、XcHB+SA 處理下隨時(shí)間的變化趨勢(shì)相似，都是在6 d 達(dá)到最大值，分別為對(duì)照的12.06、13.62、20.19、8.10、17.06 倍。在MoZK、MoZK+SA處理下的最大表達(dá)值在9 d，分別為對(duì)照的9.40、12.41 倍。SA 存在情況下，基因響應(yīng)CgTJ、MoZK、XcHB 的轉(zhuǎn)錄水平均高于對(duì)應(yīng)的單獨(dú)處理下的表達(dá)，即，CgTJ+SA、MoZK+SA、XcHB+SA處理下的表達(dá)分別為CgTJ、MoZK、XcHB 處理下轉(zhuǎn)錄水平的1.67～12.70、1.32～2.07、1.53～3.59倍（見(jiàn)圖5）。

圖5 不同處理下JrNPR1基因在莖中的表達(dá)Fig.5 The expression of JrNPR1 gene in stem under different treatments

3 討論

核桃作為我國(guó)廣大山區(qū)百姓脫貧致富、實(shí)現(xiàn)國(guó)家鄉(xiāng)村振興的重要木本油料樹(shù)種，種植面積廣泛，形成了可明顯改善種植戶(hù)經(jīng)濟(jì)收入、推動(dòng)區(qū)域經(jīng)濟(jì)發(fā)展的重要產(chǎn)業(yè)鏈。但核桃產(chǎn)業(yè)的穩(wěn)產(chǎn)、增收必須要有科學(xué)的綜合管理，包括樹(shù)體、營(yíng)養(yǎng)、病蟲(chóng)危害防控等各個(gè)方面。特別是近年時(shí)常爆發(fā)的病蟲(chóng)害，值得關(guān)注。而這就要求必須掌握核桃主要病蟲(chóng)害發(fā)生規(guī)律及防御機(jī)制，選育抗病優(yōu)良品種。NPR1 是植物病害響應(yīng)中的重要因子，隨著在較多草本植物中的功能解析被報(bào)道，基因得到了較多的關(guān)注，認(rèn)為其表達(dá)可以增強(qiáng)植株的抗病能力，是植物抗病育種的重要候選基因。但在核桃中，基因尚未被鑒定，這對(duì)揭示核桃抗病機(jī)制不利。因此，本研究中我們克隆獲得一條核桃基因（命名為），通過(guò)保守結(jié)構(gòu)域、多序列比對(duì)、同源進(jìn)化等分析發(fā)現(xiàn)與其它物種的NPR 蛋白具有較高的相似性和較近的親緣關(guān)系（見(jiàn)圖2～3），在此基礎(chǔ)上，探究基因響應(yīng)核桃主要病害的潛力。

由于啟動(dòng)子包含較多的可起轉(zhuǎn)錄調(diào)節(jié)并可預(yù)測(cè)基因潛在功能的順式作用元件，首先，本研究鑒定獲得了基因上游1 233 bp啟動(dòng)子進(jìn)行順式作用元件分析。結(jié)果發(fā)現(xiàn)，該啟動(dòng)子包含WRKY、ABRE、Dof 等轉(zhuǎn)錄因子識(shí)別的順式作用元件，且有組織特異性表達(dá)、病害、干旱等逆境響應(yīng)相關(guān)原件（見(jiàn)表1）。由于其他較多啟動(dòng)子含有相似元件，均作為預(yù)測(cè)目標(biāo)基因潛在功能的重要依據(jù)，如，海蓬子（）基因啟動(dòng)子含有MYB、WRKY、MYC等轉(zhuǎn)錄因子識(shí)別的元件，且涉及病害、激素響應(yīng)，后被證實(shí)與基因抵抗?jié)B透脅迫相關(guān)。核桃基因的啟動(dòng)子包含Dof、MYB、MYC、WRKY 等順式作用元件和激素響應(yīng)相關(guān)的ARFAT 等元件，且受NaCl、PEG及ABA 誘導(dǎo)，認(rèn)為是重要的逆境調(diào)控候選基因。同時(shí)，考慮到基因是病害響應(yīng)重要基因，可推測(cè)基因很可能參與核桃抗病響應(yīng)過(guò)程。

由于大多數(shù)基因的表達(dá)水平與植物的抗性水平密切相關(guān)。如，擬南芥基因過(guò)表達(dá)對(duì)霜霉病菌（）和丁香假單胞菌（）的抗性增強(qiáng)；將基因轉(zhuǎn)入番茄（），不僅能對(duì)番茄花葉病毒產(chǎn)生抗性，而且對(duì)多種真菌性和細(xì)菌性病害表現(xiàn)出廣譜抗性；在水稻中組織特異性表達(dá)基因提高了對(duì)稻紋枯病的抗性。過(guò)表達(dá)蘋(píng)果（）基因可增強(qiáng)蘋(píng)果對(duì)火傷病和真菌病害的抗性，并激活一系列下游抗病相關(guān)基因的表達(dá)。在葡萄（）中超表達(dá)基因能夠增強(qiáng)對(duì)白粉病的抗性。因此，為了進(jìn)一步明確的病害響應(yīng)功能潛力，對(duì)核桃進(jìn)行了主要病害—炭疽病、枝枯病、細(xì)菌性黑斑病—脅迫處理，分析在接種后不同時(shí)間葉和莖中的表達(dá)。結(jié)果顯示，基因可被膠孢炭疽菌、矩圓黑盤(pán)孢菌、黃單胞菌誘導(dǎo)表達(dá)，且體現(xiàn)了一定的病害響應(yīng)特異性和組織差異性（見(jiàn)圖4～5）。綜上可得知，可作為抵抗核桃炭疽病、枝枯病、細(xì)菌性黑斑病的重要候選基因。

SA 在植物生長(zhǎng)發(fā)育等過(guò)程中具有重要的作用。內(nèi)源SA是誘導(dǎo)植物產(chǎn)生系統(tǒng)獲得抗性（SAR）的信號(hào)物質(zhì)，可調(diào)節(jié)植物體內(nèi)多種生理代謝。外源SA 可誘導(dǎo)植株產(chǎn)生生物與非生物脅迫耐受性，其中病害抗性更為突出。如，外源噴施SA 和吸附SA 插簽兩種給藥方式均能誘導(dǎo)煙草產(chǎn)生對(duì)煙草番茄斑萎病、黑脛病、普通花葉病、野火病和赤星病的抗性。在采收前對(duì)葡萄果穗噴施SA，可顯著降低果實(shí)生長(zhǎng)期間潛伏真菌率，提高紅地球葡萄對(duì)潛伏侵染病害的抗性，抑制采后灰霉病效果明 顯。因 此，對(duì) 核 桃 進(jìn) 行 了SA、CgTJ+SA、MoZK+SA、XcHB+SA 處理，結(jié)果發(fā)現(xiàn)，基因的表達(dá)可被SA 顯著誘導(dǎo)，且基因在CgTJ+SA、MoZK+SA、XcHB+SA 處理下的表達(dá)水平要明顯高于單一的CgTJ、MoZK、XcHB 處理下的表達(dá)（見(jiàn)圖4～5）。這與其他基因響應(yīng)激素和逆境表達(dá)推測(cè)的功能潛力一致。如，核桃MYB 轉(zhuǎn)錄因子在干旱、SA 等處理下的表達(dá)認(rèn)為可作為核桃逆境響應(yīng)機(jī)制研究及抗逆育種的優(yōu)良候選基因。番茄、、基因在SA 存在情況下，NaCl誘導(dǎo)的葉中表達(dá)被顯著增強(qiáng)?？梢?jiàn)，SA是介導(dǎo)基因響應(yīng)核桃炭疽病、枝枯病、細(xì)菌性黑斑病的重要因子，基因是核桃抗病的重要候選基因。在后續(xù)研究中，將通過(guò)在木本和草本植株中過(guò)量表達(dá)來(lái)全面解析的抗病功能。