小麥NPR1基因家族鑒定及其表達分析

李雅倩，王林，宋婧含，牛紅莉，朱光，高春保，方正武，劉易科

摘要：【目的】對小麥NPR1基因家族成員進行鑒定及表達分析，為探究該家族基因的作用機制及小麥遺傳改良提供理論參考。【方法】以擬南芥NPR1家族蛋白序列為參考序列，從小麥基因組中鑒定出小麥NPR1基因家族成員，利用生物信息學軟件對其序列特征進行分析，并分別利用RNA-Seq原始數據和實時熒光定量PCR（qRT-PCR）分析小麥NPR1家族基因在不同組織及不同脅迫下的表達水平。通過STRING在線網站構建TaNPR1s蛋白的互作網絡。【結果】共鑒定獲得20個小麥NPR1基因家族成員，其編碼蛋白的不穩定指數均大于40，為不穩定蛋白;平均總親水性值（GRAVY）均為負值（除TaNPR5-D為正值外），為親水蛋白;主要分布于細胞核內，在葉綠體、線粒體、內質網和細胞質等部位也有分布;二級結構均由α-螺旋、延伸鏈、β-轉角和無規卷曲組成，以α-螺旋和無規則卷曲為主;對應的三級結構模型有10種。20個TaNPR1s蛋白均可與轉錄因子HBP-1b及未知蛋白A、B、C、D發生相互作用，這5種蛋白均含有bZIP結構域（含TGACG基序）和種子休眠特異基因結構域（DOG1）。TaNPR1s基因在不同組織中的表達模式不同，可分為在多個組織中表達、在特定的組織或發育階段表達和在不同組織發育階段均低表達或不表達，共三大類。隨機挑選的8個TaNPR1s基因中，有3個基因在禾谷鐮刀菌脅迫下表達量降低，但在白粉病菌脅迫下表達量升高;有2個基因在這兩種菌脅迫下表達量均升高，有2個基因在兩種菌脅迫下表達量均下降。TaNPR1s基因對6種非生物脅迫處理均有響應，但表達模式存在差異。【結論】小麥NPR1基因家族成員在不同組織生長發育過程和生物和非生物脅迫響應中發揮重要調控作用，且基因的可變剪接體也表現出不同組織表達特性，豐富了NPR1s蛋白功能。TaNPR1s蛋白可能通過與bZIP和DOG1結構域結合發揮其生物學功能。

關鍵詞： 小麥;NPR1基因家族;鑒定;生物信息學分析;表達分析

中圖分類號： S512.103.53? ? ? ? ? ? ? ? ? ? ? ? ? 文獻標志碼： A 文章編號：2095-1191（2021）09-2339-11

Genome-wide identification and expression analysis of

wheat NPR1 gene family

LI Ya-qian1，2， WANG Lin1，2， SONG Jing-han1，2， NIU Hong-li1， ZHU Guang2，

GAO Chun-bao1，2， FANG Zheng-wu1*， LIU Yi-ke2*

（1College of Agronomy， Yangtze University/Hubei Collaborative Innovation Center for Grain Industry， Jingzhou， Hubei? 434025， China; 2Food Crops Institute，Hubei Academy of Agricultural Sciences/Wheat Disease Biology

Research Station on Central China，Ministry of Agriculture and Rural Affairs/Hubei Engineering and

Technology Research Center of Wheat， Wuhan? 430064， China）

Abstract：【Objective】The identification and expression analysis of NPR1 gene family members in wheat provided theoretical reference for exploring the mechanism of this family gene and wheat genetic improvement. 【Method】Wheat NPR1 gene family members from wheat genome were identified? based on NPR1 family protein sequence from arabidopsis. The sequence characteristics were analyzed by bioinformatics software， and the expression levels of wheat NPR1 family genes under different tissues and different stresses were analyzed by RNA-Seq data and real-time fluorescence quantitative PCR（qRT-PCR）. The online website STRING was used to predict and analyze the interaction relationship between related proteins of TaNPR1 family members. 【Result】A total of 20 members of wheat NPR1 gene family were identified， and the TaNPR1s encoded by them were unstable proteins with instability index greater than 40. The average total hydrophilic va-lue （GRAVY） was negative （except TaNPR5-D was positive）， which was hydrophilic protein. It was mainly distributed in the nucleus， chloroplast， mitochondria， endoplasmic reticulum and cytoplasm. The secondary structure was composed of α-helix， extended strand， β-turn and random coil， and α-helix and random coil were the main components. There were 10 kinds of tertiary structure models. All 20 TaNPR1s proteins could interact with the transcription factor HBP-1b and unknown proteins A， B， C and D， the five proteins all contained bZIP domain（containing TGACG motif） and seed dormancy specific gene domain（DOG1）. The expression pattern of TaNPR1s gene in different tissues was different， which could be divided into three categories：expressed in multiple tissues， expressed in specific tissues or developmental stages， and all low expression or no expression in different tissue development stages. Among the 8 randomly selected TaNPR1s genes， the expression levels of 3 genes decreased under Fusarium graminearum stress， but increased under powdery mildew stress. The expression levels of two genes increased under both strains， and the expression levels of two genes decreased under both strains. TaNPR1s gene responded to all six abiotic stress treatments， but their expression patterns were different. 【Conclusion】Members of wheat NPR1 gene family play an important role in regulating the growth and development of different tissues and the response to biological and abiotic stresses， and the variable spliceosomes of genes also show different tissue expression characteristics， enriching the function of NPR1s protein. TaNPR1s protein may play its bio-logical function by binding to bZIP and DOG1 domains.

Key words： wheat; NPR1 gene family; identification; bioinformatics analysis; expression analysis

Foundation item： Hubei Provincial Special Project of Central Government Guides Local Science and Technology Development（2020ZYYD011）;National Wheat Industry Technology System Construction Project（CARS-03）

0 引言

【研究意義】小麥是世界廣泛種植的主要農作物之一，也是人類的主要口糧來源，其生長發育過程中易遭到各種病原微生物的侵擾，從而影響其產量和品質（Wang et al.，2018;陳文燁等，2020）。植物在生物體進化過程中為了抵抗各種生物和非生物脅迫形成了各種免疫防御機制（Jones and Dangl，2006;陳文燁等，2020），其中，植物的系統性抗病機理存在兩種發生類型：一種是致病菌所誘導的系統獲得性抗性（Systemic acquired resistance，SAR）;另一種為非致病性微生物所誘導的系統誘導性抗性（Induced systemic resistance，ISR）（Kogel and Langen，2005;鄭世仲等，2020）。非表達病程相關蛋白（Non-expresser of pathogenesis related genes 1，NPR1）是植物防御信號中的一個重要調節因子，廣泛參與植物體內的免疫防御反應（Mhamdi，2019），其在SAR反應過程中是水楊酸（Salicylic acid，SA）信號通路的重要組成部分，也是植物病程相關蛋白（Pathogenesis-rela-ted proteins，PRs）基因表達和SAR的激活子（Cao et al.，1997;Zhou et al.，2000;Backer et al.，2019）;其在茉莉酸（Jasmonic acid，JA）介導的ISR反應過程中是JA信號途徑的負調控因子（Pieterse et al.，1998;Spoel et al.，2003）。此外，NPR1是植物響應病原菌侵染的核心基因，參與調控植物體內SA與JA的平衡（韓永光等，2018）。在植物對冷環境適應過程中NPR1基因參與調節冷誘導基因表達（Olate et al.，2018）。因此，對小麥基因組中的NPR1基因家族成員進行鑒定及分析，對探究NPR1基因在小麥抗逆過程中的作用機制及培育小麥抗性品種具有重要意義。【前人研究進展】NPR1基因首次從擬南芥中被克隆，該基因啟動子區域的W-box元件和WRKY轉錄因子結合進而調控NPR1基因的轉錄，其編碼蛋白的C端含有錨蛋白重復序列（Ank）和NPR1-like C-terminal結構域，N端有一個BTB/POZ（Broad complex Tramtrack and Bric-a-brac/Pox virus and zinc finger）結構域（Cao et al.，1994;Cao et al.，1997;Aravind and Koonin，1999）。目前擬南芥基因組測序已完成，共鑒定出6個AtNPR1s基因（Initiative，2000），其中AtNPR1和AtNPR2是SA受體，在植物免疫中起轉錄共激活因子的作用（Castelló et al.，2018）;AtNPR3和AtNPR4也是SA受體，在植物防御中起轉錄輔抑制因子的作用（Zhang et al.，2006）;AtNPR5和AtNPR6基因又被稱為AtBOP1和AtBOP2基因，因缺乏NPR1-like C-terminal結構域，在系統發育進化樹中分支差異較大，二者主要在擬南芥生長發育過程和形態建成中發揮重要作用（Hepworth et al.，2005;Mikael et al.，2005）。目前大量研究證實，NPR1基因對植物的抗病性中發揮重要作用，如在蘋果中過量表達MpNPR1基因能加強植株對火疫病、黑星病菌及膠銹菌的抗病性（Malnoy et al.，2007）;將擬南芥AtNPR1基因轉入草莓中可增強草莓對炭疽病、白粉病及葉斑病等的抗性（Silva et al.，2015）;水稻中過表達OsNPR1基因能增強植株對葉枯病菌的抗性（Yuan et al.，2007）。抑制煙草NPR1基因表達可導致植株喪失煙草花葉病毒（TMV）抗性（Liu et al.，2002）。保守的條銹菌蛋白可與小麥NPR1蛋白發生相互作用，并降低病原體對致病相關基因（Pathogenesis-related genes，PR）的誘導（Wang et al.，2016）。NPR1基因在受到病原物侵染的條件下，才能激活下游防御基因表達，在未被誘導物或病原體激活時，PR基因的表達并未顯著提高（蔡韡韡等，2019;Tripathi et al.，2019）。【本研究切入點】目前，NPR1基因家族已在番木瓜（Peraza-Echeverria et al.，2012）、蘋果（焦鵬，2016）、香蕉（任陪娣等，2019）、小麥（Liu et al.，2019） 和油菜（Wang et al.，2020）等多種植物中鑒定分析。但針對小麥NPR1基因在生物和非生物脅迫下的表達模式分析較少且缺乏驗證，同時對該家族基因的蛋白特性還缺乏系統研究。【擬解決的關鍵問題】以擬南芥NPR1家族蛋白序列為參考序列，從小麥基因組中鑒定出小麥NPR1基因家族成員，利用生物信息學軟件對其序列特征進行分析，并分別利用RNA-Seq原始數據和實時熒光定量PCR（qRT-PCR）分析小麥NPR1家族基因在不同組織及不同脅迫下的表達水平;通過STRING在線網站構建TaNPR1s蛋白的互作網絡，為探究該家族基因的作用機制及小麥遺傳改良提供理論參考。

1 材料與方法

1. 1 試驗材料

供試的小麥品種為鄂麥170，由湖北省農業科學院糧食作物研究所提供。RNA提取試劑TRIzol購自賽默飛世爾科技（中國）有限公司。反轉錄試劑盒（Prime ScriptTM RT reagent Kit with gDNA Eraser）和實時熒光定量PCR試劑盒（SYBR Premix Ex Taq）購自TaKaRa公司。主要設備儀器：人工氣候箱（PQX智能型，寧波萊福科技有限公司）、凝膠成像系統（Bio-Rad，美國）和CFX Connect Real-Time System（Bio-Rad，美國）。

1. 2 樣品處理及采集

參考Jiang等（2019）的方法對小麥種子消毒處理后進行無菌培養，生長出來的小麥幼苗在人工氣候箱中培養至兩葉一心時期。對小麥幼苗分別進行生物脅迫[禾谷鐮刀菌（PH-1）、白粉病菌（Bgt strain E09）]和非生物脅迫[熱（42 ℃）、冷（4 ℃）、脫落酸（ABA，100 μmol/L）、水楊酸（SA，100 μmol/L）、茉莉酸甲酯（MeJA，100 μmol/L）和氯化鈉（NaCl，150 μmol/L）]處理24 h，對照組（CK）不作脅迫處理。收集幼苗根和葉部位組織，每3株混作一個樣，液氮速凍處理后，-80 ℃保存備用，用于分析不同逆境脅迫下TaNPR1基因的表達特性。

1. 3 小麥NPR1基因家族成員鑒定

以擬南芥AtNPR1（AT1G64280.1）、AtNPR2（AT4G26120.1）、AtNPR3（AT5G45110.1和AT5G451 10.2）、AtNPR4（AT4G19660.1和AT4G19660.2）、AtNPR5（AT2G41370.1）和AtNPR6（AT3G57130.1）的蛋白序列為查詢序列，在小麥基因組數據庫中進行BLASTp比對（E值<1×10-10）（Altschul et al.，1997）。將檢索到的序列提交到Pfam（http：//pfam.xfam.org/）（Finn et al.，2006）和InterProScan（http：//www.ebi.ac.uk/interpro/）（Mulder and Apweiler，2007）網站，篩選出同時含有BTB/POZ和Ankyrin結構域的蛋白，即為小麥NPR1基因家族成員。

1. 4 生物信息學分析及互作預測

小麥NPR1蛋白的理化性質通過ExPASy Ser-ver10（https：//prosite.espasyp.org/）（Wilkins et al.，1999）進行預測。利用WoLF PSORT：Protein Subcellular Localization Prediction（https：//www.genscript.com/wolf-psort.html）（Horton et al.，2007）在線預測亞細胞定位。利用SOPMA（https：//npsa-prabi.ibcp.fr/cgi-bin/npsa_automat.pl？page=npsa_sopma.html）預測蛋白二級結構。使用SWISS-MODEL（https：//swissmodel.expasy.org/interactive#alignment）（Schwede et al.，2003）對蛋白的三級結構進行預測。通過功能蛋白關聯網絡在線網站STRING（https：//stringdb.org）（Szklarczyk et al.，2018）預測小麥NPR1基因家族成員相關蛋白的互作關系，可信度分數為0.70。

1. 5 轉錄組數據分析

從NCBI的Short Read Archive（SRA）數據庫中下載關于不同組織生長發育相關的RNA-Seq原始數據（SBA編號：PRJEB25639），以及不同逆境脅迫處理的RNA-Seq原始數據（SBA編號：條銹菌PRJEB12497、禾谷鐮刀菌PRJEB12358、干旱PRJNA306536和高溫─干旱共處理PRJNA257938），并通過HISAT2將其映射到小麥基因組中。轉錄組數據采用TPM值（Transcripts Per Million）來表示基因的表達量。TPM值低于1的看作基因不表達或表達量很低。利用R語言“pheatmap”包，基于TPM值繪制基因表達譜，以Log2（TPM+1）值繪制NPR1基因差異表達熱圖。

1. 6 實時熒光定量PCR（qRT-PCR）表達分析

利用反轉錄試劑盒將RNA反轉錄合成cDNA。采用SYBR Premix Ex Taq進行qRT-PCR檢測。反應體系（20.0 μL）：cDNA 2.0 μL，正、反向引物各0.8 μL，50×ROX Reference Dye II 0.4 μL，TB Green Premix Ex Taq II 10.0 μL，滅菌水補足至20.0 μL。擴增程序：95 ℃ 3 min;95 ℃ 10 s，60 ℃ 30 s（同時收集信號值），共進行40個循環，溶解曲線分析為65 ℃。β-Actin作為相對定量的內參基因。利用Primer Premier 5.0設計內參基因和小麥NPR1基因家族成員的引物序列，如表1所示。采用2-ΔΔCt方法進行數據計算，每個樣本技術重復3次，基因的相對表達量用平均值±標準差表示。對照的相對表達量為1。

2 結果與分析

2. 1 小麥NPR1基因家族成員鑒定及序列特征分析結果

以8個擬南芥NPR1家族蛋白序列為參考序列，在小麥基因組數據庫進行BLASTp比對分析，并結合Pfam和SMART網站去除不具有BTB/POZ和Ank保守結構域的序列，共鑒定獲得20個小麥NPR1基因家族成員（其中包含3個可變剪接體：TaNPR3.1-A、TaNPR3.1-B、TaNPR4.1-B），并根據其與擬南芥NPR1基因家族的進化關系進行分組及命名（表2）。20個小麥NPR1家族蛋白（TaNPR1s）的不穩定指數均大于40，為不穩定蛋白;平均總親水性值（GRAVY）均為負值（除TaNPR5-D為正值外），為親水蛋白;主要分布于細胞核內，在葉綠體、線粒體、內質網和細胞質等部位也有分布。

2. 2 TaNPR1s蛋白結構預測結果

由表3可知，TaNPR1s蛋白二級結構由α-螺旋、延伸鏈、β-轉角和無規卷曲4種元件組成，以α-螺旋和無規則卷曲為主，其中TaNPR1-A、TaNPR1-B、TaNPR1-D、TaNPR2-A、TaNPR2-D、TaNPR3-A、TaNPR3-B、TaNPR3-D、TaNPR4-A、TaNPR4-B和TaNPR4-D蛋白以α-螺旋比例最大，TaNPR5-A、TaNPR5-B、TaNPR5-D、TaNPR6-A、TaNPR6-B和TaNPR6-D蛋白中α-螺旋和無規卷曲所占比例較相近。此外，TaNPR5-A、TaNPR5-B、TaNPR5-D、TaNPR6-A、TaNPR6-B和TaNPR6-D蛋白結構不含NPR1-like C-terminal結構域，與該家族其他蛋白差異較大。

利用SWISS-MODEL構建小麥NPR1蛋白的三級結構模型，共有10種，如表3所示。TaNPR1-A、TaNPR1-B、TaNPR1-D、TaNPR2-A和TaNPR2-D蛋白的二級結構相似，但三級結構差異較大，除TaNPR1-B和TaNPR1-D對應的結構模型相同外，其余蛋白對應的結構模型均不相同;TaNPR3-A、TaNPR3.1-A、TaNPR4-A、TaNPR4.1-B和TaNPR4-D蛋白的二級結構相似，對應結構模型均為5y4d.1。5y4d.1三維結構在I和II組均有分布，其原因可能是與I和II組的遺傳關系較近（Liu et al.，2019）。第III組內除TaNPR6-A和TaNPR6-B蛋白三級結構差異較大外，其余蛋白三級結構相似（圖1）。TaNPR3-A和TaNPR3-B蛋白分別與其對應的可變剪接體TaNPR3.1-A和TaNPR3.1-B三級結構類似，而TaNPR4-B蛋白與其對應的可變剪接體TaNPR4.1-B三級結構差異相對較大。

2. 3 TaNPR1s蛋白互作網絡構建

為了更好地理解基因功能，使用STRING在線數據庫構建TaNPR1s蛋白的互作網絡（圖2）。所有TaNPR1s蛋白均可與轉錄因子HBP-1b及未知蛋白A、B、C、D發生相互作用。通過Pfam檢測發現，這5種蛋白均含有bZIP結構域（含TGACG基序）和種子休眠特異基因結構域（DOG1）。TaNPR1可能通過與這兩個結構域發生相互作用從而發揮其功能。

2. 4 TaNPR1s基因在不同組織中的表達分析結果

從NCBI下載與中國春小麥籽粒、穗、葉、莖和根5個組織在生長發育過程中的RNA-Seq原始數據（SBA編號：PRJEB25639），對TaNPR1s基因在小麥生長發育過程中表達水平進行分析。在20個TaNPR1s基因中，大部分基因在小麥生長發育過程中均有表達，TPM值為1.00～31.32（圖3）。按照在不同組織的表達情況，TaNPR1s基因可分為三大類：第一類在多個組織中表達，包括TaNPR3.1-A、TaNPR3-B、TaNPR3-D、TaNPR4-A、TaNPR4-B、TaNPR4-D和TaNPR4.1-B共7個基因;第二類僅在特定的組織或發育階段表達，包括TaNPR1-A、TaNPR1-B、TaNPR1-D、TaNPR2-A、TaNPR2-D、TaNPR5-A、TaNPR5-B、TaNPR5-D、TaNPR6-A、TaNPR6-B和TaNPR6-D共11個基因;第三類在不同組織發育階段均低表達或不表達，包括TaNPR3-A和TaNPR3.1-B。以第二類中的TaNPR5-B基因為例，利用BAR在線軟件繪制其基因表達的電子熒光象形圖（圖4），該圖直觀顯示了TaNPR5-B基因在小麥早、中、晚3個生長發育階段在不同組織的表達情況。由上述結果推測小麥NPR1基因家族成員在小麥不同組織生長發育過程中發揮重要作用。

2. 5 TaNPR1s基因在不同脅迫條件下的表達分析結果

從NCBI下載與生物脅迫（SBA編號：條銹菌PRJEB12497;禾谷鐮刀菌PRJEB12358）和非生物脅迫（SBA編號：干旱PRJNA306536;高溫─干旱共處理PRJNA257938）相關的RNA-Seq原始數據，對TaNPR1s基因在不同逆境脅迫下的表達模式進行分析，結果圖5所示。TaNPR1s基因在不同逆境脅迫條件下表達模式不同。TaNPR1-A、TaNPR1-B、TaNPR1-D和TaNPR4.1-B基因在接種條銹菌后1 d時表達量上升，而后表達量下降;TaNPR3.1-A和TaNPR3-D基因在接種后3 d時表達量下降，而在5 d時表達量上升。禾谷鐮刀菌侵染后，TaNPR1-A、TaNPR1-B、TaNPR1-D、TaNPR5-A、TaNPR5-B和TaNPR6-B基因的表達量總體上呈上升趨勢，TaNPR4-B、TaNPR4.1-B和TaNPR5-D基因的表達量總體呈下降趨勢。在干旱脅迫處理條件下，TaNPR1-A、TaNPR1-B、TaNPR1-D、TaNPR4.1-B和TaNPR4-D基因表達量下降。此外，TaNPR1-A、TaNPR1-B、TaNPR1-D、TaNPR4.1-B和TaNPR4-D基因在高溫或高溫—干旱共處理條件下表達量變化不顯著，而TaNPR1-A、TaNPR1-B和TaNPR1-D基因在干旱處理1 h時表達量上升，隨后下降，TaNPR4.1-B、TaNPR4-A和TaNPR4-D基因在干旱處理6 h時表達量下降。TaNPR3-B基因在干旱處理條件下表達量下降，在高溫或高溫—干旱共處理條件下表達量上升。TaNPR4-B干旱處理條件下表達量下降，在高溫或高溫—干旱共處理條件下表達量先上升后下降。綜上所述，TaNPR1s基因參與逆境環境的響應及信號轉導，且不同NPR1基因在不同脅迫下基因的表達量變化不同。

隨機挑選8個TaNPR1s基因，利用qRT-PCR檢測其在生物脅迫（禾谷鐮刀菌和白粉病菌）下的表達情況，結果如圖6-A和圖6-B所示。與對照相比，TaNPR2-D、TaNPR3.1-A和TaNPR3-D基因在禾谷鐮刀菌脅迫下的情況下相對表達量降低，而在白粉病菌脅迫下相對表達量升高;TaNPR3-B和TaNPR3.1-B基因在兩種菌脅迫下相對表達量均升高，TaNPR2-A和TaNPR4-A基因在兩種菌的脅迫下相對表達量均下降。

利用qRT-PCR檢測上述隨機挑選的8個TaNPR1s基因在非生物脅迫（冷、熱、NaCl、SA、MeJA和ABA）下的表達情況，結果如圖6-C～圖6-H所示。TaNPR1s基因對6種非生物脅迫處理均有響應，但表達模式存在明顯差異。與對照相比，在冷脅迫下，除TaNPR3-B和TaNPR3.1-A基因的相對表達量上升外，其余基因相對表達量均下降;在熱和SA脅迫下，所有基因的相對表達量均下降;在NaCl脅迫下，TaNPR2-A、TaNPR2-D和TaNPR3.1-B基因的相對表達量上升，其余基因相對表達量下降;在MeJA脅迫下，除TaNPR2-A基因外，其他基因的相對表達量下降;在ABA脅迫下，除TaNPR3-A基因外，其余基因的相對表達量下降。

3 討論

NPR1是植物SAR及PRs基因表達的激活子，也參與SA與JA/乙烯信號途徑的相互作用，是調節植物整體抗病性的重要作用因子，參與植物多種抗性代謝通路，在植物生長發育和抗病性研究中具有重要作用（Cao et al.，1997），因而植物NPR1基因的研究受到國內外學者的普遍關注。本研究在Liu等（2019）對于NPR1基因系統分類分析的基礎上，對小麥NPR1家族蛋白的理化性質、結構預測、互作網絡及在不同組織和不同逆境脅迫下的表達情況進行系統分析，初步揭示了小麥NPR1基因家族成員的結構和功能，為該基因家族在小麥遺傳改良中的應用提供了理論參考。前人研究發現，NPR1蛋白的BTB結構域與TGA2的抑制區結合會導致TGA2抑制PRs基因表達的功能缺失，從而提升植物抗病性，可見，NPR1典型特征結構域對其功能發揮極為重要（Rochon et al.，2006;Boyle et al.，2009;李微巍，2014）;NPR1蛋白的NPR1-like C-terminal結構域含有核定位信號，對定位于細胞核內的NPR1參與轉錄共激活發揮重要作用（Kinkema et al.，2000）。本研究通過構建蛋白互作網絡發現，所有TaNPR1s蛋白均與轉錄因子HBP-1b及未知蛋白A、B、C、D發生相互作用，經Pfam檢測發現這5種蛋白均含有bZIP結構域和DOG1結構域，這兩類結構域均屬于TGA轉錄因子。在細胞核內，NPR1蛋白能與TGA亞家族成員（含有TGACG基序，屬于bZIP轉錄因子家族）結合，這些轉錄因子參與SAR反應，在SA的誘導下能激活下游防衛基因PR1的表達（Zhang et al.，1999;Fan and Dong，2002）。而DOG1結構域可以通過結合ABA信號傳遞的負調控因子，增加對ABA的敏感性從而引起種子休眠，在種子成熟和發育過程中發揮作用（宋松泉等，2020）。說明TaNPR1蛋白可能通過與bZIP或DOG1結構域結合，進而在生物和非生物脅迫過程中發揮功能。

本研究發現，部分TaNPR1s基因（TaNPR3.1-A、TaNPR3-B、TaNPR3-D、TaNPR4.1-B和TaNPR4-D）在不同組織中均有較高的相對表達量，即使在不同脅迫下，這些基因仍保持較高的相對表達量，表明這些基因在小麥不同生長階段發揮重要的調控作用，可能直接或間接參與小麥的生長發育過程和形態建成。此外，本研究還發現，盡管TaNPR1s基因在不同生物脅迫和非生物脅迫下的相對表達量不同，但大多數TaNPR1s基因的相對表達量明顯上調或下調，表明這些基因參與逆境環境的響應及信號轉導（Ma and Bohnert，2007），且TaNPR3.1-A、TaNPR3-B和TaNPR4.1-B基因在不同發育時期均具有較高的相對表達量，但其對應的可變剪接體TaNPR3-A、TaNPR3.1-B和TaNPR4-B基因的相對表達量較低。可變剪切體可使同一個基因翻譯成不同的蛋白，進而表現出不同的功能，可變剪接是基因在轉錄過程中普遍存在的現象，是導致蛋白質功能多樣性的重要原因之一（Wang et al.，2015）。小麥的3個NPR1基因與其對應的可變剪切體在逆境脅迫下表達情況不同，說明可變剪切增加了小麥NPR1家族的轉錄本，豐富了其蛋白的遺傳多樣性。NPR1介導的逆境脅迫下植物信號通路非常復雜，今后應進一步比較TaNPR1s基因可變剪切體之間的差異，并探究其具體作用機制，對挖掘NPR1基因在響應逆境脅迫中的關鍵位點具有重大意義。

4 結論

小麥NPR1基因家族成員在不同組織生長發育過程和生物和非生物脅迫響應中發揮重要調控作用，且基因的可變剪接體也表現出不同組織表達特性，豐富了NPR1s蛋白功能。TaNPR1s蛋白可能通過與bZIP和DOG1結構域結合發揮其生物學功能。

參考文獻：

蔡韡韡，曾志芳，魏春，楊華麗，佘文琴，陳桂信. 2019. ?NPR1同源基因全長cDNA的分離與表達分析[J]. 園藝學報，46（3）：567-576. [Cai W W，Zeng Z F，Wei C，Yang H L，She W Q，Chen G X. 2019. Isolation and expression analysis of full-length cDNA of NPR1 homologous genes in Nai（Prunus salicina var. cordata）[J]. Acta Horticulturae Sinica，46（3）：567-576.] doi：10.16420/j.issn.0513-353x.2019-0030.

陳文燁，楊帆，劉永偉，董福雙，趙和，柴建芳，呂孟雨，周碩. 2020. 小麥TaNAC-B072基因的克隆和表達分析[J]. 江西農業學報，32（8）：1-7. [Chen W Y，Yang F，Liu Y W，Dong F S，Zhao H，Chai J F，Lü M Y，Zhou S. 2020 Cloning and expression analysis of TaNAC-B072 in wheat[J]. Acta Agriculturae Jiangxi，32（8）：1-7.] doi：10.19386/ j.cnki.jxnyxb.2020.08.01.

韓永光，馬利剛，趙樂，馮衛生，鄭曉珂. 2018. 植物抗性基因NPR1研究進展[J]. 安徽農業科學，46（26）：18-20. [Han Y G，Ma L G，Zhao L，Feng W S，Zheng X K. 2018. Research progress on resistance gene NPR1 in plants[J]. Journal of Anhui Agricultural Sciences，46（26）：18-20.] doi：10.13989/j.cnki.0517-6611.2018.26.006.

焦鵬. 2016. 蘋果NPR1基因家族的全基因組鑒定與表達分析[D]. 楊凌：西北農林科技大學. [Jiao P. 2016. Genome-wide identification and expression analysis of NPR1 gene family in apple[D]. Yangling：Northwest A & F University.]

李微巍. 2014. 楊樹NPR-like基因家族抗病相關基因的鑒定及功能研究[D]. 合肥：安徽農業大學. [Li W W. 2014. Identification of disease-resistance cancer-critical genes and functional research of NPR-Like family in poplar[D]. Hefei：Anhui Agricultural University.]

任陪娣，謝碧玉，陳翠，賈彩紅，王靜毅，張建斌. 2019. 香蕉 NPR1基因家族的鑒定及在枯萎病菌脅迫下表達分析[J]. 植物遺傳資源學報，20（6）：1621-1629. [Ren P D，Xie B Y，Chen C，Jia C H，Wang J Y，Zhang J B. 2019. Genome-wide survey of the NPR1 gene family in banana（Musa spp.） and expression pattern analysis upon infection with Foc TR4[J]. Journal of Plant Genetic Resources，20（6）：1621-1629.] doi：10.13430/j.cnki.jpgr.20190416 002.

宋松泉，劉軍，徐恒恒，劉旭，黃薈. 2020. 脫落酸代謝與信號傳遞及其調控種子休眠與萌發的分子機制[J]. 中國農業科學，53（5）：857-873. [Song S Q，Liu J，Xu H H，Liu X，Huang H. 2020. ABA metabolism and signaling and their molecular mechanism regulating seed dormancy and germination[J]. Scientia Agricultura Sinica，53（5）：857-873.] doi：10.3864/j.issn.0578-1752.2020.05.001.

鄭世仲，王培育，張春柳，江勝滔，江金蘭，顏沛沛，葉煒，賴鐘雄. 2020. 文心蘭NPR1基因的克隆與誘導抗性過程中的表達分析[J]. 福建農業學報，35（2）：140-149. [Zheng S Z，Wang P Y，Zhang C L，Jiang S T，Jiang J L，Yan P P，Ye W，Lai Z X. 2020. NPR1 Cloning and expression du-ring development of induced systemic resistance in Oncidium hybridum[J]. Fujian Journal of Agricultural Scien-ces，35（2）：140-149.] doi：10.19303/j.issn.1008-0384.2020. 02.004.

Altschul S F，Madden T L，Sch?ffer A A，Zhang J，Zhang Z，Miller W，Lipman D J. 1997. Gapped BLAST and PSI-BLAST：A new generation of protein database search programs[J]. Nucleic Acids Research，25（17）：3389-3402. doi：10.1093/nar/25.17.3389.

Aravind L，Koonin E V. 1999. Fold prediction and evolutiona-ry analysis of the POZ domain：Structural and evolutio-nary relationship with the potassium channel tetramerization domain[J]. Journal of Molecular Biology，285（4）：1353-1361. doi：10.1006/jmbi.1998.2394.

Backer R，Naidoo S，Berg N. 2019. The NONEXPRESSOR OF PATHOGENESIS-RELATED GENES 1（NPR1） and related family：mechanistic insights in plant disease resistance[J]. Frontiers in Plant Science，10：102. doi：10.3389/fpls. 2019.00102.

Boyle P，Su E L，Rochon A，Shearer H L，Murmu J，Chu J Y，Fobert P R，Després C. 2009. The BTB/POZ domain of the Arabidopsis disease resistance protein NPR1 interacts with the repression domain of TGA2 to negate its function[J]. The Plant Cell，21（11）：3700-3713. doi：10.1105/tpc.109.069971.

Cao H，Bowling S A，Gordon A S，Dong X. 1994. Characteri-zation of an Arabidopsis mutant that is nonresponsive to inducers of systemic acquired resistance[J]. The Plant Cell，6（11）：1583. doi：10.1105/tpc.6.11.1583.

Cao H，Glazebrook J，Clarke J D，Volko S，Dong X. 1997. The Arabidopsis NPR1 gene that controls systemic acquired resistance encodes a novel protein containing ankyrin repeats[J]. The Cell，88（1）：57-63. doi：10.1016/S0092-8674（00）81858-9.

Castelló M J，Medina-puche L，Lamilla J，Tornero P. 2018. NPR1 paralogs of Arabidopsis and their role in salicylic acid perception[J]. PLoS One，13（12）：e0209835. doi：10. 1371/journal.pone.0209835.

Fan W H，Dong X N. 2002. In vivo interaction between NPR1 and transcription factor TGA2 leads to salicylic acid-mediated gene activation in Arabidopsis[J]. The Plant Cell，14（6）：1377-1389. doi：10.1105/tpc.001628.

Finn R D，Mistry J，Schuster-b?ckler B，Griffiths-jones S，Hollich V，Lassmann T，Moxon S，Marshall M，Khanna A，Durbin R，Eddy S R，Sonnhammer E L L，Bateman A. 2006. Pfam：Clans，web tools and services[J]. Nucleic Acids Research，34（S1）：247-251. doi：10.1093/nar/gkj149.

Hepworth S R，Zhang Y，Mckim S，Li X，Haughn G W. 2005. BLADE-ON-PETIOLE-dependent signaling controls leaf and floral patterning in Arabidopsis[J]. The Plant Cell，17（5）：1434-1448. doi：10.1105/tpc.104.030536.

Horton P，Park K J，Obayashi T，Fujita N，Harada H，Adams-collier C J，Nakai K. 2007. WoLF PSORT：Protein loca-lization predictor[J]. Nucleic Acids Research，35（S2）：585-587. doi：10.1093/nar/gkm259.

Initiative T A G. 2000. Analysis of the genome sequence of the fowering plant Arabidopsis thaliana[J]. Nature，408（6814）：796-815. doi：10.1038/35048692.

Jiang W Q，Yang L，He Y Q，Zhang H T，Li W，Chen H G，Ma D F，Yin J L. 2019. Genome-wide identification and transcriptional expression analysis of superoxide dismutase（SOD） family in wheat[J]. PeerJ，7（2）：e8062. doi：10.7287/peerj.preprints.27402.

Jones J D G，Dangl J L. 2006. The plant immune system[J]. Nature，444（7117）：323-329. doi：10.1038/nature05286.

Kinkema M，Fan W，Dong X. 2000. Nuclear localization of NPR1 is required for activation of PR gene expression[J]. The Plant Cell，12（12）：2339-2350. doi：10.1105/tpc. 12.12.2339.

Kogel K，Langen G. 2005. Induced disease resistance and gene expression in cereals[J]. Cellular Microbiology，7（11）：1555-1564. doi：10.1111/j.1462-5822.2005.00592.x.

Liu X，Liu Z G，Niu X H，Xu Q，Yang L. 2019. Genome-wide identification and analysis of the NPR1-Like gene family in bread wheat and its relatives[J]. International Journal of Molecular Sciences，20（23）：5974. doi：10.3390/ ijms20235974.

Liu Y，Schiff M，Marathe R，Dinesh-kumar S P. 2002. Tobacco Rar1，EDS1 and NPR1/NIM1 like genes are required for N-mediated resistance to tobacco mosaic virus[J]. The Plant Journal，30（4）：415-429. doi：10.1046/j.1365-313x. 2002.01297.x.

Ma S，Bohnert H J. 2007. Integration of Arabidopsis thaliana stress-related transcript profiles，promoter structures，and cell-specific expression[J]. Genome Biology，8（4）：R49. doi：10.1186/gb-2007-8-4-r49.

Malnoy M，Jin Q，Borejsza-wysocka E E，He S Y，Aldwinckle H S. 2007. Overexpression of the apple MpNPR1 gene confers increased disease resistance in Malus x domestica[J]. Molecular Plant-Microbe Interactions，20（12）：1568-1580. doi：10.1094/MPMI-20-12-1568.

Mhamdi A. 2019. NPR1 has everything under control[J]. Plant Physiology，181（1）：6-7. doi：10.1104/pp.19.00890.

Mikael N，Mattias H，Ove N. 2005. The BLADE ON PETIOLE genes act redundantly to control the growth and development of lateral organs[J]. Development，132（9）：2203-2213. doi：10.1242/dev.01815.

Mulder N，Apweiler R. 2007. InterPro and InterProScan：Tools for protein sequence classification and comparison[J]. Methods in Molecular Biology，396：59-70. doi：10. 1007/978-1-59745-515-2_5.

Olate E，Jiménez-gómez J M，Holuigue L，Salinas J. 2018. NPR1 mediates a novel regulatory pathway in cold acclimation by interacting with HSFA1 factors[J]. Nature Plants，4（10）：811-23. doi：10.1038/s41477-018-0254-2.

Peraza-echeverria S，Santamaría J M，Fuentes G，Menéndez-cerón M，Vallejo-reyna M，Herrera-valencia V A. 2012. The NPR1 family of transcription cofactors in papaya：Insights into its structure，phylogeny and expression[J]. Genes & Genomics，34（4）：379-390. doi：10.1007/s13258- 011-0218-7.

Pieterse C M J，Van Wees S C M，Van Pelt J A，Knoester M，Laan R，Gerrits H，Weisbeek P J，Van Loon L C. 1998. A novel signaling pathway controlling induced systemic resistance in Arabidopsis[J]. The Plant Cell，10（9）：1571-1580. doi：10.1105/tpc.10.9.1571.

Rochon A，Boyle P，Wignes T，Fobert P R，Despres C. 2006. The coactivator function of Arabidopsis NPR1 requires the core of its BTB/POZ domain and the oxidation of C-terminal cysteines[J]. The Plant Cell，18（12）：3670-3685. doi：10.1105/tpc.106.046953.

Schwede T，Kopp J，Guex N，Peitsch M C. 2003. SWISS-MODEL：An automated protein homology-modeling ser-ver[J]. Nucleic Acids Research，31（13）：3381-3385. doi：10.1093/nar/gkg520.

Silva K J P，Brunings A，Peres N A，Mou Z，Folta K M. 2015. The Arabidopsis NPR1 gene confers broad-spectrum disease resistance in strawberry[J]. Transgenic Research，24（4）：693-704. doi：10.1007/s11248-015-9869-5.

Spoel S H，Koornneef A，Claessens S M C，Korzelius J P，Van Pelt J A，Mueller M J，Buchala A J，Métraux J，Brown R，Kazan K，Van Loon L C，Dong X，Pieterse C. 2003. NPR1 modulates cross-talk between salicylate- and jasmonate-dependent defense pathways through a novel function in the cytosol[J]. The Plant Cell，15（3）：760-770. doi：10.1105/tpc.009159.

Szklarczyk D，Gable A L，Lyon D，Junge A，Wyder S，Huerta-cepas J，Simonovic M，Doncheva N T，Morris J H，Bork P，Jensen L J，Mering C. 2018. STRING v11：Protein-protein association networks with increased coverage，supporting functional discovery in genome-wide experimental datasets[J]. Nucleic Acids Research，47（D1）：D607-D613. doi：10.1093/nar/gky1131.

Tripathi D，Raikhy G，Kumar D. 2019. Chemical elicitors of systemic acquired resistance—Salicylic acid and its functional analogs[J]. Current Plant Biology，17：48-59. doi：10.1016/j.cpb.2019.03.002.

Wang X D，Bi W S，Gao J，Yu X M，Wang H Y，Liu D Q. 2018. Systemic acquired resistance，NPR1，and pathoge-nesis-related genes in wheat and barley[J]. Journal of Integrative Agriculture，17（11）：60345-60347. doi：10.1016/S2095-3119（17）61852-5.

Wang X D，Yang B J，Li K，Kang Z S，Cantu D，Dubcovsky J. 2016. A conserved Puccinia striiformis protein intera-cts with wheat NPR1 and reduces induction of pathogenesis-related genes in response to pathogens[J]. Molecular Plant-Microbe Interactions，29（12）：977-989. doi：10.1094/MPMI-10-16-0207-R.

Wang Y，Liu J，Huang B O，Xu Y M，Li J，Huang L F，Lin J，Zhang J，Min Q H，Yang W M，Wang X Z. 2015. Mechanism of alternative splicing and its regulation[J]. Biome-dical Reports，3（2）：152-158. doi：10.3892/br.2014.407.

Wang Z，Ma L Y，Li X，Zhao F Y，Sarwar R，Cao J，Li Y L，Ding L N，Zhu K M，Yang Y H，Tan X L. 2020. Genome?wide identification of the NPR1 like gene family in Brassica napus and functional characterization of BnaNPR1 in resistance to Sclerotinia sclerotiorum[J]. Plant Cell Reports，39：709-722. doi：10.1007/s00299-020-02525-z.

Wilkins M R，Gasteiger E，Bairoch A，Sanchez J C，Williams K L，Appel R D. 1999. Protein identification and analysis tools in the ExPASy Server[J]. Methods in Molecular Biology，112（112）：531-552. doi：10.1385/1-59259-584-7：531.

Yuan Y X，Zhong S H，Li Q，Zhu Z R，Lou Y G，Wang L Y，Wang J J，Wang M Y，Li Q L，Yang D L，He Z H. 2007. Functional analysis of rice NPR1-like genes reveals that OsNPR1/NH1 is the rice orthologue conferring disease resistance with enhanced herbivore susceptibility[J]. Plant Biotechnology Journal，5（2）：313-324. doi：10.1111/j.1467- 7652.2007.00243.x.

Zhang Y L，Cheng Y T，Qu N，Zhao Q G，Bi D L，Li X. 2006. Negative regulation of defense responses in Arabidopsis by two NPR1 paralogs[J]. The Plant Journal，48（5）：647-656. doi：10.1111/j.1365-313X.2006.02903.x.

Zhang Y L，Fan W H，Kinkema M，Li X，Dong X N. 1999. Interaction of NPR1 with basic leucine zipper protein transcription factors that bind sequences required for salicylic acid induction of the PR-1 gene[J]. Proceedings of the National Academy of Sciences，96（11）：6523-6528. doi：10.1073/pnas.96.11.6523.

Zhou J，Trifa Y，Silva H，Pontier D，Lam E，Jyoti S，Klessig D F，Zhou J M，Shah J. 2000. NPR1 differentially interacts with members of the TGA/OBF family of transcription factors that bind an element of the PR-1 gene required for induction by salicylic acid[J]. Molecular Plant-Microbe Interactions，13（2）：191-202. doi：10.1016/S1369-5266（00）80027-8.

（責任編輯 陳 燕）