馬尾松NAC轉錄因子基因PmNAC8的克隆及表達分析

鎬青青 姚圣 劉佳禾 陳佩珍 張夢洋 季孔庶

（南京林業大學 南方現代林業協同創新中心 林木遺傳與生物技術教育部重點實驗室，南京 210037）

植物經常遭受自然界中有害的環境脅迫（如干旱、高溫、冷害、凍害以及高鹽等非生物脅迫），對植物的生長和發育產生不利影響，對農林業生產和生態環境造成嚴重的損失［1-2］。為了應對不利的環境條件，植物通過調節自身基因的表達來抵御脅迫，從而獲得脅迫耐受性［3-4］。AP2/EREBP、WRKY、bZIP、MYB和NAC等類型的轉錄因子（transcription factors，TFs）已被證明在植物非生物脅迫反應中起著關鍵作用［5］。NAC轉錄因子家族是植物特有的一類轉錄因子，它的命名取自最早發現的矮牽牛（Petunia hybrida Vilm.）NAM 基因、擬南芥（Arabidopsis thaliana）ATAF1基因和CUC2基因的首字母［6］。

在擬南芥中首次報道了NAC轉錄因子參與干旱反應的調控，據報道，過表達ANAC019、ANAC055和ANAC072可以調控擬南芥干旱應答基因的表達，提高轉基因植物的耐旱性［7］。水稻 OsNAP［8］、小麥 TaNAC69［9］和 玉 米 ZmSNAC1［10］等 NAC 轉 錄因子均參與干旱反應的調控。此外，一些干旱相關的NAC基因也參與植物激素介導的信號通路，如脫落酸、茉莉酸、水楊酸和乙烯的信號通路［11］。Mahmood等［12-13］研究表明ANAC032促進植株葉片衰老的同時，在非生物脅迫（如蔗糖、鹽和氧化）下，又能抑制花青素的合成；而擬南芥ANC046既能促進葉片衰老也能促進花青素的合成［14］。在馬鈴薯（Solanum tuberosum）中，StNAC2對鎘脅迫有一定的應答和調節作用［15］。在水稻（Oryza sativa L.）中，過表達OsNAC5可以增大水稻根系的直徑，提高對干旱的耐受性，從而增加水稻產量［16］。在豆科植物中間錦雞兒（Caragana intermedia）中，過表達CiNAC1可促進擬南芥葉片衰老，并且在乙烯處理后與葉綠素降解和衰老相關的基因表達量明顯提高［17］。在山葡萄（Vitis amurensis）中，VaNAC26 可以通過上調茉莉酸合成基因和茉莉酸信號通路中的基因表達，提高擬南芥體內的茉莉酸含量來增強對干旱的耐受力［18］。在桃樹（Amygdalus persica）中，PpNAC72在桃樹果實中高表達，推測該基因可能參與果實成熟［19］。

馬尾松（Pinus massoniana）是松科（Pinaceae）松屬常綠喬木，具有速生、豐產、適應性強、經濟價值高等優良特性，是我國南方最主要的優質針葉用材樹種之一［20］。由于我國南方紅壤中有效磷含量低，頻繁發生干旱、低溫等非生物脅迫，使得馬尾松造林存活率和生產力降低，造成了嚴重的經濟損失。NAC作為植物特有的一類轉錄因子，在植物生物和非生物脅迫反應的調控中起著關鍵作用。而且，現階段對NAC基因的研究主要集中在模式植物，如擬南芥、水稻和大豆等，對木本植物的研究較為薄弱，尤其是馬尾松等針葉樹種。因此，對馬尾松NAC轉錄因子進行功能鑒定是必要的。

本研究運用RT-PCR的方法成功克隆了馬尾松PmNAC8，對其進行生物信息學分析、亞細胞定位，組織特異性以及在干旱、損傷等非生物脅迫和外源激素方面的響應表達模式分析，以期為馬尾松抗逆分子育種提供候選基因資源和理論依據。

1 材料與方法

1.1 材料

馬尾松植物材料取自南京林業大學樹木園中15年生實生苗和實驗室盆栽的30 d幼苗；本氏煙草在光照培養箱中培養，24℃長日照（18 h光照/6 h黑暗），28 d后取長勢較好且一致的植株作為試驗材料。

多糖多酚植物總RNA提取試劑盒（離心柱型）購自天根生化科技（北京）有限公司；大腸桿菌菌株Trans1-T1感受態、克隆載體pEASY-Blunt Simple購自北京全式金生物技術有限公司；PrimeScriptTMRT Master Mix（Perfect Real Time）、SMARTer RACE5′/3′Kit和 3′-Full RACE Core Set with PrimeScrip RTase Kit購 自 TaKaRa公 司；DNA marker、Green Taq Mix、Phanta Max Super-Fridelity DNA Polymerase均購自南京諾唯贊生物科技有限公司；本研究所有引物均由Primer 5軟件設計，并由擎科生物科技有限公司合成。

1.2 方法

1.2.1 RNA的提取和cDNA合成 參照TIANGEN多糖多酚植物總RNA提取試劑盒說明書，提取30 d盆栽馬尾松幼苗全株RNA，所得產物經1.2%凝膠電泳檢測，保證RNA完整性，同時，使用超微量紫外分光光度計（NANODROP 2000，ThermoFisher公司）對RNA濃度進行測量，若OD260/280、OD260/230比值為1.8-2.1，說明RNA質量較好。選取2 μg濃度較高且純度較好的RNA，利用全式金反轉錄試劑盒反轉錄合成cDNA第一條鏈，作為基因克隆的模板。根 據 SMARTer?RACE 5′Kit 和 3′-Full Core Set with PrimeScriptTM RTase試劑盒說明書，合成5′和3′末端模板。

1.2.2 PmNAC8的克隆 從NCBI（national center for biotechnology information）搜索不同植物的NAC序列，與馬尾松當年生嫩枝轉錄組（PRJNA655997）數據［21］相比較，篩選出目的基因，利用Primer Premier 5.0軟件設計中間片段引物PmNAC8-Mid（表1）。以馬尾松幼苗的cDNA為模板進行PCR反應，總反應體系為 cDNA（50 ng/μL）2 μL、PmNAC8-Mid F 1 μL、PmNAC8-Mid R 1 μL、5×GXL Buffer（Mg2+Plus）10 μL、DNTP Mix-ture 5 μL，ddH2O 補至 50 μL ；反應程序為 98℃ 3 min ；98℃ 30 s，58℃ 30 s，72℃2 min，35個循環；72℃ 7 min，4℃保存。

以測序獲得的PmNAC8中間片段的序列為參照，利用Primer 5.0軟件分別設計3′RACE特異性引物 ：3′RACE-PmNAC8-F、3′RACE-PmNAC8-R 以 及5′RACE 特異 性引物 5′RACE-PmNAC8-F、5′RACEPmNAC8-R（表1），分別以 3′-Full Core Set with Prime-ScriptTMRTase 試劑盒和 SMARTer?RACE5′ Kit試劑盒反轉錄得到的cDNA為模板，分別與接頭引物 3′RACE Outer、3′RACE Inner 和 5′RACE Outer、5′RACE Inner進行巢式PCR，各反應均遵循試劑盒說明。將反應產物經過1.2%瓊脂糖凝膠電泳驗證、回收，與pEASY-Blunt載體連接，轉化到大腸桿菌感受態細胞Trans1-T1，挑取陽性單克隆送上海杰李生物公司進行測序。用Seqman軟件拼接序列獲得基因全長。

為驗證拼接結果的準確性，在序列兩端設計引物，分別為：PmNAC8-F和PmNAC8-R（表1），以馬尾松cDNA為模板進行PCR反應，產物經連接轉化后，送至上海杰李公司進行測序。

1.2.3 PmNAC8的生物信息學分析 利用ORF Finder（https://www.ncbi.nlm.nih.gov/orffinder/） 在線預測開放閱讀框，確定基因編碼序列；利用ExPASy（https://web.expasy.org/protparam/） 在 線 軟件分析氨基酸序列的基本理化性質；利用ExPASy TMpred（https://embnet.vital-it.ch/software/TMPRED_form.html） 和 MHMM 2.0（http://www.cbs.dtu.dk/services/TMHMM/）在線分析蛋白質的跨膜區和擴膜方向；利用SOPMA（https://npsa-prabi.ibcp.fr/cgibin/npsa_automat.pl?page=npsa_sopma.html）在線軟件分析預測蛋白質二級結構；利用SWISS-MODEL（https://swissmodel.expasy.org/interactive）在線預測蛋白質三級結構；利用SignalP-4.0 Server（http://www.cbs.dtu.dk/services/SignalP-4.0/）在線軟件預測蛋白的信號肽；利用Motif Scan（https://myhits.isb-sib.ch/cgibin/motif_scan）在線軟件分析蛋白的結構域。

將PmNAC8編碼氨基酸序列在NCBI數據庫中進行blastp比對，下載blast score較高的氨基酸序列，利用 MEGA X［22］，Evolview（https://www.evolgenius.info/ evolview/） 和 iTOL5.0（https://itol.embl.de/） 在線軟件構建系統進化樹。

1.2.4 PmNAC8同義密碼子偏好性分析 密碼子是生物體內遺傳信息從DNA到蛋白質準確傳遞的載體［23-24］。不同生物對密碼子的選擇具有不同的偏好性，當宿主不具備或僅含有稀少的外源基因密碼子時，外源基因的表達可能降低或終止。因此，對密碼子使用偏好性進行分析，有助于選擇合適的外源表達宿主，保證外源基因順利表達。利用CodonW軟 件（http://codonw.sourceforge.net/）計 算 PmNAC8密碼子的使用特性參數，包括A3s、C3s、U3s、GC含量、有效密碼子數（effective number of codons，ENc）、密碼子中第3位堿基的GC含量（GC3s）、相對同義密碼子使用度（relative synonymous codon usage，RSCU）和密碼子適應指數（codon adaptation index，CAI）等。運用EMBOSS中的CUSP在線程序計算密碼子的使用頻率。

擬南芥、煙草（Nicotiana tabacum）、歐洲山楊（Populus tremula）、大腸桿菌（Escherichia coli）和釀酒酵母（Saccharomyces cerevisiae）的基因組密碼子使用頻率數據來自密碼子使用數據庫（http://www.kazusa.or.jp/codon）。

1.2.5 馬尾松PmNAC8亞細胞定位載體的構建及農桿菌轉化 使用Xba I和BamH I酶切瞬時表達載體PBI121-GEP，同時利用高保真酶擴增含有酶切位點的目的基因，上下游引物分別為：PmNAC8-GFP-F和PmNAC8-GFP-R（表1），將PCR擴增產物與雙酶切產物回收純化后重組并轉化大腸桿菌，選取陽性菌液送至杰李公司測序。將構建成功的融合表達載體質粒（35S∷PmNAC8-GFP）轉化農桿菌GV3101，以只含有GFP標簽的質粒農桿菌菌株為對照。將P19菌液與目標菌液等比例混合后按照李尊強等［25］方法將農桿菌重懸液注射于生長狀態較好的煙草葉片中，在人工氣候培養箱中避光培養48 h后，利用激光共聚焦顯微鏡LSM710（Zeiss，Jena，Germany）觀察GFP融合蛋白的位置。

1.2.6 qRT-PCR分析 為了研究馬尾松PmNAC8的表達水平，根據天根生化科技有限公司的多糖多酚植物總RNA提取試劑盒說明書，分別提取了15年生馬尾松不同組織及8種處理下30 d幼苗的RNA。利用超微量分光光度計測定總RNA濃度與OD260/280、OD260/230的值并進行凝膠電泳檢測其彌散程度（方法同1.2.1），將檢測結果符合要求的RNA作為模板進行反轉錄。8種處理包括非生物脅迫和激素處理，分別是機械損傷、15% PEG6000、10 mmol/L MeJA、50 μmol/L ETH、10 mmol/L H2O2、2 mmol/L GA、1 mmol/L SA和400 μmol/L ABA，其中機械損傷的處理方法是通過切斷松針的上半部分、15% PEG6000處理方法是將植物浸泡在溶液中模擬干旱，其他脅迫和處理均采用噴淋植物表面的方式進行［26］。每個處理選擇3株長勢均勻的幼苗作為3個生物重復，每個植株收集3個樣品作為3個技術重復，每隔0、3、6、12和24 h取幼苗針葉作為樣品，經液氮速凍后，-80℃保存備用，其中，在0 h收集的樣品未經任何處理用作對照組。利用Primer 5.0軟件設計實時定量PCR特異性引物Q-PmNAC8-F和Q-PmNAC8-R，以Zhu等［27］篩選出的馬尾松較穩定的Tubulin alpha（GenBank登錄號：KM496535.1）作為內參基因，并設計內參引物（表1），進行定量PCR反應。將得到的cDNA模板稀釋20倍，備用。反應體系和反應程序參考 Hieff UNICON?qPCR SYBR Green Master Mix試劑盒說明書。反應結束后，將數據導出，參照2-ΔΔCT法計算該基因的相對表達量，使用Microsoft Office Excel 2019和GraphPad Prism 8.0等軟件進行數據分析及基因表達圖的繪制。

1.2.7 啟動子克隆和序列分析 根據在馬尾松近緣物種火炬松基因組數據庫檢索到的PmNAC8啟動子序列，設計2條特異性引物GSP1和GSP2（表1），利用天根植物基因組提取試劑盒提取15年生馬尾松基因組DNA，按照TaKaRa生物公司染色體步移技術試劑盒進行PCR擴增，將PCR產物經電泳檢測，回收清晰條帶，連接克隆載體pEASY-Blunt，轉化大腸桿菌，挑選正確的單克隆擴大培養后送至上海杰李生物公司測序。根據測序序列設計特異性引物PmNAC8-Pro（表1），以馬尾松基因組DNA為模板克隆啟動子全長。

表1 本研究中使用的引物序列Table 1 Sequences of primers used in this study

利 用 PlantCARE（http://bioinformatics.psb.ugent.be/webtools/plantcare/html/）在線軟件預測PmNAC8的啟動子內部順式作用元件。

2 結果

2.1 PmNAC8的克隆與序列分析

根據馬尾松轉錄組數據庫中找到的目的基因的EST序列［21］，設計特異性引物后，依次進行中間片段擴增、3′RACE擴增、5′RACE擴增，之后利用SeqMan軟件拼接獲得該基因的全長cDNA序列1 726 bp，為驗證拼接結果的正確性，對其全長及ORF區域進行擴增（圖1）。測序結果顯示拼接結果正確，該基因的ORF為1 209 bp，3′端序列148 bp，5′端序列369 bp，預測其可編碼402個氨基酸，且具有NAC蛋白特有的NAM結構域，因此將其命名為PmNAC8，基因登錄號為MZ291447。

圖1 馬尾松PmNAC8的全長及ORF的擴增Fig.1 Full-length and ORF amplification of PmNAC8 gene in P.massoniana

2.2 PmNAC8的生物信息學分析

2.2.1 蛋白理化性質及跨膜區預測分析 將PmNAC8序列通過NCBI在線軟件進行ORF預測發現，其含有開放閱讀框1 209 bp，編碼402個氨基酸，以ATG為起始密碼子，TGA為終止密碼子。通過Motif Scan在線分析PmNAC8蛋白的結構域，顯示在其58-191位氨基酸序列是NAC蛋白特有的NAM結構域，其E-value值為2.80E-26，說明其屬于NAC轉錄因子（圖2-A）。

利用ExPASy ProtParam在線軟件預測蛋白質的理化性質，結果表明，PmNAC8蛋白由碳、氫、氧、氮、硫5種元素組成，總分子式為C1941H3011N559O639S20，分子量為45 042.82 Da，由402個氨基酸殘基組成，理論等電點為5.11，其中帶正電荷（Arg+Lys）的氨基酸殘基數為43，帶負電荷（Asp+Glu）的氨基酸殘基數為65。利用ProtScale在線軟件進行蛋白疏水性分析，其中橫坐標為蛋白質氨基酸殘基，縱坐標為該位點氨基酸殘基疏水性值，正值表示疏水性氨基酸，負值表示親水性氨基酸。結果顯示，分值為負的氨基酸多于分值為正的氨基酸，即親水大于疏水，因此，推測PmNAC8蛋白為親水性蛋白（圖2-B）。通過在線軟件SignalP預測PmNAC8蛋白信號肽，結果表明，該蛋白不存在信號肽（圖2-C）。

圖2 PmNAC8蛋白的生物信息學分析Fig.2 Bioinformatics analysis of PmNAC8 protein

通過ExPASy TMpred和TMHMM在線軟件預測分析蛋白跨膜區、擴膜方向。結果（圖3-A、B）表明，PmNAC8蛋白為非跨膜蛋白。

圖3 PmNAC8蛋白跨膜區的預測分析Fig.3 Transmembrane region prediction of PmNAC8 protein

2.2.2 蛋白二級結構、三級結構預測 利用SOPMA在線預測分析蛋白的二級結構。結果（圖4-A）表明，PmNAC8蛋白的二級結構由無規則卷曲（Cc）、α-螺旋（Hh）、延伸鏈（Ee）和β-折疊構成。其中，無規則卷曲有251個氨基酸，占序列的62.44%；α-螺旋有70個氨基酸，占序列的17.41%；延伸鏈有61個氨基酸，占序列的15.17%；β-折疊有20個氨基酸，占序列的4.98%。通過SWISS-MODEL網站使用同源建模法在線預測蛋白質的三級結構（圖4-B）。

圖4 PmNAC8蛋白二級（A）、三級（B）結構預測Fig.4 Secondary（A）and tertiary（B）structure predictions of PmNAC8 protein

2.2.3 系統進化樹分析 利用進化樹設計軟件MEGA5.0采用鄰接法（neighbor-joining，NJ），設置參數自展值（Bootstrap值）為1 000次進行進化分析，結果（圖5）表明，馬尾松的NAC蛋白與火炬松、北美云杉的遺傳距離較近，說明它們在進化的過程中親緣關系較為密切，而與其他雙子葉植物的遺傳距離較遠，說明與雙子葉植物的NAC蛋白在進化過程中親緣關系較遠。

圖5 PmNAC8蛋白系統發育進化樹Fig.5 Phylogenetic analysis of PmNAC8 protein

2.3 PmNAC8密碼子偏好性分析

2.3.1 GC、GC3s、ENc和CAI分析 利用CondonW程序分析PmNAC8密碼子的CAI值、ENc值，同時運用EMBOSS中的CUSP在線程序計算總GC含量以及GC3s有效密碼子數［28］。結果（表2）可知，PmNAC8基因的ENc值為50.19，ENc值偏大，說明PmNAC8在編碼氨基酸時密碼子使用偏好性較弱，表達水平相對較低。CAI值為0.187，進一步表明PmNAC8對密碼子的選擇偏好性較弱。GC含量和GC3s值分別為31.10%和42.50%，說明密碼子偏好以A/T結尾。

表2 PmNAC8密碼子選擇偏好性參數Table 2 Codon selection preference parameters of PmNAC8 gene

2.3.2 相對同義密碼子使用度 在馬尾松PmNAC8的密碼子中，有29個密碼子的RSCU值大于1（表3），為馬尾松PmNAC8偏好密碼子，并且其中有23個密碼子是以A/T結尾；其中偏好性較強的有UCA、GCA和UGA（RSCU＞2.0）；然而編碼甲硫氨酸（Met）的密碼子ATG、編碼色氨酸（Trp）的密碼子TGG、編碼亮氨酸（Leu）的密碼子CTC以及編碼異亮氨酸（Ile）的ATA，其RSCU值均等于1，表明馬尾松PmNAC8基因對Met、Trp、Leu和Ile的使用沒有偏好性。大多數G或C堿基結尾的密碼子的RSCU值和Fraction值均較低，表明這些密碼子在該基因中的使用頻率較低。

表3 PmNAC8同義密碼子相對使用度Table 3 RSCU of PmNAC8 gene

2.4 PmNAC8外源宿主的選擇

將PmNAC8密碼子使用頻率分別與擬南芥、煙草、釀酒酵母、大腸桿菌和歐洲山楊基因組密碼子的使用頻率進行分析比較（表4），將比值大于或等于2.0和小于或等于0.5的密碼子作為篩選標準統計差異個數。結果表明，PmNAC8與煙草、擬南芥和歐洲山楊3種模式植物之間的絕大部分密碼子使用偏好性差異較小，其中，與煙草、擬南芥、歐洲山楊基因組密碼子的使用頻率差異較大的分別有12、14和13個，說明在PmNAC8遺傳轉化試驗中，煙草更適合作為異源表達的受體。與3種植物相比，PmNAC8與大腸桿菌和釀酒酵母2種微生物的基因組密碼子偏好性差異較大，尤其與大腸桿菌密碼子偏好性差異很大，而與釀酒酵母密碼子偏好性差異小于大腸桿菌，可能由于釀酒酵母為真核生物，而大腸桿菌為原核生物，可見，相較于大腸桿菌，酵母更適合PmNAC8微生物內異源表達。

表4 PmNAC8與部分模式生物基因組密碼子使用偏好性比較Table 4 Comparison of codon usage preference between PmNAC8 and other model organisms

續表Continued

2.5 PmNAC8亞細胞定位試驗結果

構建PBI121-PmNAC8-GFP基因的融合表達載體，以空載體PBI121-GFP為對照，利用煙草瞬時表達系統，通過農桿菌將載體轉入煙草葉片中，在人工氣候培養箱中避光培養48 h后，利用激光共聚焦顯微鏡LSM710觀察GFP融合蛋白的位置。由圖6可知，PmNAC8蛋白質主要定位在細胞核中，而對照組細胞核、細胞膜和細胞質中均能觀察到綠色熒光信號。

圖6 PmNAC8亞細胞定位結果Fig.6 Subcellular localization results of PmNAC8

2.6 PmNAC8的組織特異性表達分析

利用實時熒光定量PCR技術分析了PmNAC8在15年生馬尾松不同組織中的表達情況。結果（圖7）表明，PmNAC8在花、幼葉、老葉、幼莖、老莖、根、木質部和韌皮部中均有表達，但表達水平不同，在根中的表達水平顯著高于在其他組織，是花的21.4倍，是幼葉和老葉的9.0倍和7.6倍。因此，PmNAC8主要在馬尾松的根中表達。

圖7 馬尾松PmNAC8在不同組織中的表達水平Fig.7 PmNAC8 expression level in the different tissues of P.massoniana

2.7 PmNAC8在非生物脅迫下的表達分析

以 10 mmol/L MeJA、50 μmol/L ETH、10 mmol/L H2O2、2 mmol/L GA、1 mmol/L SA、15% PEG6000、機械損傷和400 μmol/L ABA處理后的30 d馬尾松盆栽幼苗cDNA為模板，檢測PmNAC8在脅迫處理下表達情況（圖8）。在10 mmol/L MeJA和15%PEG6000處理下PmNAC8表達量均受到顯著抑制，在3 h時表達量約為對照組的1/2和1/4，脅迫至6h后表達量幾乎完全被抑制。在50 μmol/L ETH、10 mmol/L H2O2、2 mmol/L GA、1 mmol/L SA、損傷和400 μmol/L ABA處理下，PmNAC8表達量總體呈現先升高后降低的趨勢。在ETH、H2O2、GA、SA處理下，PmNAC8均在3 h時表達量最大，與對照組差異顯著，分別是對照組的1.13、1.86、6.89和1.15倍，其中在2 mmol/L GA處理下PmNAC8的表達量顯著高于對照組和其他脅迫組。在損傷和ABA處理下，PmNAC8的表達量分別在6和12 h時最高，與對照組差異顯著，分別是對照組的1.62和1.30倍。

圖8 馬尾松PmNAC8在不同脅迫下的表達水平Fig.8 PmNAC8 expression of P.massonniana under different stress

2.8 PmNAC8啟動子順式作用元件分析

為明確PmNAC8啟動子是否具有活性，克隆了PmNAC8上游1 492 bp啟動子序列（圖9），命名為PmNAC8-Pro。PlantCARE分析顯示，該啟動子除含有TATA-box（21個）、CAAT-box（34個）等核心啟動子元件和光響應元件（6個）以外，還包含與非生物脅迫相關的元件，如參與干旱誘導的MYB結合位點MBS（1個）、抗氧化反應元件ARE（2個），以及植物激素相關響應元件，如赤霉素響應元件（1個）（表5）。基于上述結果，推測馬尾松PmNAC8啟動子可能參與調控非生物脅迫以及響應植物激素脅迫的功能。

表5 PmNAC8啟動子中的順式作用元件Table 5 Cis-acting elements in the PmNAC8 promoter

圖9 PmNAC8啟動子擴增電泳Fig.9 Amplification electrophoresis of PmNAC8 promoter

3 討論

NAC是植物特有的一類轉錄因子，是植物長發育及適應環境脅迫的重要轉錄調節因子。許多研究表明，NAC轉錄因子家族在植物中具有不可替代的生物學功能，如葉片衰老［29］、果實成熟［30］、胚珠發育［31］、逆境脅迫應答［32］和激素信號轉導［33］等。目前，對NAC轉錄因子的報道主要集中在模式植物中，如擬南芥、水稻、玉米等，有關馬尾松NAC家族基因功能的研究和報道較少。本研究從馬尾松中克隆了PmNAC8，將其編碼的氨基酸序列與其他物種NAC進行比對，結果表明，PmNAC8蛋白具有NAC轉錄因子家族典型的NAM結構域，N端具有150 bp組成的高度保守序列，C端具有高度特異性。另外，蛋白質的二級結構是其整個空間結構的基礎，在蛋白質折疊過程的早期形成，為蛋白質三級結構的形成奠定基礎，對亞細胞定位有一定作用［34］，預測PmNAC8蛋白的二級結構發現主要以無規則卷曲結構為主，這與大多數NAC蛋白二級結構相似。系統進化樹結果表明其編碼的蛋白與火炬松、北美云杉高度同源。對PmNAC8蛋白結構分析發現，其為親水性蛋白，且不存在信號肽，與大多數植物NAC蛋白相似［35］。亞細胞定位結果表明，PmNAC8定位于細胞核，屬于典型的核蛋白，表明PmNAC8可能在細胞核中發揮重要功能。

探究PmNAC8基因密碼子偏好性，可以揭示其基因的密碼子組成特性，為深入研究PmNAC8的功能提供依據。本研究利用生物信息學的方法對PmNAC8的密碼子偏好性進行分析，結果表明，該基因有29個密碼子的RSCU值大于1，其中23個密碼子是以A/T結尾，且該基因GC含量和GC3s值分別為31.10%和42.50%，這一結果均符合馬尾松基因總體上偏好第3位堿基為A/T密碼子的特征［36］。煙草、擬南芥、歐洲山楊等模式植物常用于基因功能的探究，大腸桿菌常作為基因的原核表達系統受體，而酵母菌常被用于真核系統受體，由于馬尾松轉基因平臺尚未構建成功，因此，需要借助相應遺傳轉化體系進行異源表達。在異源表達過程中，為實現外源基因的成功表達并提高其表達量，應盡量選擇密碼子使用偏好性差異較小的作為受體，本研究通過比較PmNAC8與擬南芥、煙草、歐洲山楊、釀酒酵母、大腸桿菌的基因組密碼子使用頻率的差異，發現酵母真核表達系統更適合作為PmNAC8基因的表達系統，而煙草相較于擬南芥和歐洲山楊更適合作為PmNAC8的異源表達受體。

本研究對馬尾松PmNAC8進行qRT-PCR分析發現，該基因在植物不同組織中均有表達，但在根中的表達水平明顯高于其他組織，其次是老莖和幼莖。同時，PmNAC8受損傷及植物激素GA、ABA的誘導表達，推測該基因可能參與植物非生物脅迫和激素的應答調控。啟動子影響外源基因的表達水平，是基因工程表達載體的一個重要元件，對基因的調控具有重要意義［37］。如水稻Os08PTS啟動子序列能夠調控Os08PTS基因在水稻莖和種子胚中表達，驅動下游基因轉錄［38］；谷子膜結合NAC家族基因SiNAC啟動子參與生長素、甲基茉莉酸和氧化脅迫應答［39］等。本研究以馬尾松為材料經染色體步移法克隆了馬尾松PmNAC8啟動子，長度為1 492 bp。序列分析發現PmNAC8-Pro啟動子序列中包括多種光響應元件、生物脅迫響應元件及植物激素相關響應元件，如赤霉素（GA）響應元件，且序列分析結果與qRT-PCR結果一致，推測PmNAC8可能受到非生物脅迫和激素等的誘導響應。

4 結論

從馬尾松中克隆出一個新的NAC基因PmNAC8，全長1 726 bp，編碼402個氨基酸，其編碼蛋白質分子量為45 042.82 Da，等電點為5.11，屬于親水性蛋白，且不存在信號肽，有明顯的NAM結構域，屬于典型的NAC轉錄因子家族成員；密碼子使用偏性較弱，偏好使用A/T結尾的密碼子；酵母真核表達系統較大腸桿菌原核表達系統更適合PmNAC8異源表達，模式植物煙草較擬南芥和歐洲山楊更適合做PmNAC8的遺傳轉化受體。該蛋白定位在細胞核內，與火炬松、北美云杉同源性較高；PmNAC8在植物不同組織中均有表達，可能參與植物非生物脅迫和激素的應答調控。