摘要：NAC轉錄因子是植物中特有的轉錄因子，其在N端含有一段高度保守的特異結構域，對植物生長發(fā)育、脅迫反應和次生代謝產(chǎn)物生物合成等過程具有重要功能。研究基于高州油茶第三代全長轉錄組數(shù)據(jù)信息，通過生物信息學的方法篩選獲得51個NAC轉錄因子，為深入研究基因家族特征，對其進行理化性質分析、多列比對分析、蛋白系統(tǒng)進化樹、保守基序和基因結構分析、亞細胞定位，以及二級蛋白域結構預測、蛋白磷酸化位點分析。結果表明：NAC蛋白存在氨基酸數(shù)量差異，且氨基酸序列間疏水性呈現(xiàn)為親水性蛋白，其蛋白結構中具有5個保守的亞結構域。NAC轉錄因子系統(tǒng)進化樹分析表明116個基因可劃分為10個亞族群，部分高州油茶NAC蛋白具有特異性聚為單一亞族群，其他與擬南芥和水稻融合為不同亞族群。同時NAC蛋白均具有NAM超家族結構。亞細胞定位表明51個高州油茶NAC蛋白均位于細胞核，且二級結構中隨機卷曲結構為主，同時含有α-螺旋和β-轉角結構等。蛋白磷酸化位點都具有絲氨酸（Ser）、蘇氨酸（Thr）和酪氨酸（Tyr）。

關鍵詞：高州油茶；NAC轉錄因子；生物信息學分析

中圖分類號：S565.4; Q943.2 文獻標識碼：A 文章編號：1006-060X（2024）07-0007-08

Identification and Bioinformatics Analysis of NAC Transcription Factor Family in Camellia drupifera Loureiro

SHEN Hong-jian，DENG Jin-qing，LIAO Bo-yong，ZHANG Bi-pei，LI Yong-quan，WANG Yi

（Zhongkai University of Agriculture and Engineering， Guangzhou 510220， PRC）

Abstract： NAC transcription factors are unique transcription factors in plants， which contain a highly conserved segment of specific structural domains at the N-terminus and have important functions in the processes of plant growth and development， stress response， and secondary metabolite biosynthesis. In this study， based on the information of the third-generation full-length transcriptome data of Camellia drupifera Loureiro （C. drupifera）， 51 NAC transcription factors were obtained by bioinformatics screening. In order to study the characteristics of these gene families， their physicochemical properties were analyzed， and multilinear comparison， protein phylogenetic trees， conserved motifs， gene structure， and subcellular localization were investigated. The structure of the secondary protein domains was predicted， and the protein phosphorylation sites were analyzed. The results show that NAC proteins differ in the number of amino acids， and the hydrophobicity between amino acid sequences is presented as hydrophilic proteins. Their protein structures have five conserved sub-structural domains. The phylogenetic tree analysis of the NAC transcription factors shows that 116 genes can be classified into 10 sub-populations. Part of the NAC proteins of C. drupifera is clustered as a single sub-population with specificity， and the other parts of the NAC proteins are fused with Arabidopsis thaliana and rice as a mixed sub-population. Meanwhile， all the NAC proteins have NAM superfamily structures. The subcellular localization shows that all 51 NAC proteins of C. drupifera are located in the nucleus， and the secondary structure is dominated by a randomly curled structure， with the structures of α-helix and β-rotor present. All the protein phosphorylation sites contain serine （Ser）， threonine （Thr）， and tyrosine （Tyr）. The gene family characteristics of NAC transcription factors from C. drupifera provide a reference for the subsequent research on their gene functions.

Key words：C. drupifera; NAC transcription factors; bioinformatics analysis

引用格式：沈紅建，鄧錦清，廖柏勇，等. 高州油茶NAC轉錄因子家族鑒定及生物信息學分析[J]. 湖南農(nóng)業(yè)科學，2024（7）：7-14.

DOI：10.16498/j.cnki.hnnykx.2024.007.002

收稿日期：2024-02-23

基金項目：廣東省林業(yè)科技創(chuàng)新項目（2023KJCX006）；廣州市科技計劃項目（202201011754）；廣東省重點領域研發(fā)計劃項目（2020B

020215003）

作者簡介：沈紅建（1998—），男，湖南株洲市人，碩士研究生，專業(yè)為農(nóng)藝與種業(yè)。

通信作者：廖柏勇

NAC是植物中重要的一類轉錄因子[1]，其廣泛分布于陸生植物當中。NAC轉錄因子N端具有高度保守的NAC結構域，C端則呈現(xiàn)出高度多樣性[2]。

在擬南芥（Arabidopsis thaliana）、水稻（Oryza sativa L.）、

石榴（Punica granatum）等物種均有NAC轉錄因

子[3-5]。許多研究表明，NAC轉錄因子在植物的生長發(fā)育和生殖發(fā)育以及生物和非生物脅迫響應中都起著重要作用[6-9]。例如，在番茄中鑒定出的SlNAC35基因，對于干旱脅迫、鹽脅迫、細菌病原體均可誘導促進SlNAC35的表達；通過對SlNAC35的過表達提高了煙草干旱和鹽脅迫的耐受度，也提高了對病原體的抵抗能力性[10-11]。在植物信號轉導過程中，蛋白磷酸化與蛋白的翻譯和修飾有著重要關系，在參與調(diào)節(jié)酶的活性和植物的信號傳遞中起著重要作用[12]。蛋白質磷酸化可以影響蛋白質之間的相互作用，可以改變蛋白的結構，從而影響他們與其他蛋白相互作用，磷酸化修飾作為一種動態(tài)調(diào)節(jié)機制，可以被蛋白磷酸酶催化反應去除磷酸基團，從而實現(xiàn)關于細胞周期信號傳遞調(diào)控[13]。

油茶（Camellia oleifera）屬于山茶科（Theaceae）山茶屬（CamelliaL），為常綠闊葉樹種，適宜種植在丘陵、山區(qū)地帶[14]。干旱脅迫影響油茶樹體生長，造

成枝葉枯萎，葉片和果實脫落[10]。長期干旱脅迫容易

造成樹體衰老加速直至植株死亡[10，15]。目前，高州油茶NAC轉錄因子的研究未見報道，筆者基于高州

油茶全長轉錄組數(shù)據(jù)，利用生物信息學方法對高州油茶NAC轉錄因子結構進行分析[2]，為提高高州油茶抗脅迫相關的NAC基因功能研究提供理論參考。

1 材料與方法

1.1 植物材料

試驗樣品采自廣東省韶關市曲江區(qū)國有小坑林場（113°35′E，24°15′N）。分別在“B255”的3株無性系樹上采集花藥、葉片、葉芽、枝條、果皮、種仁鮮樣放入無紡布袋，再放入液氮冷凍帶回實驗室-80℃保存，稱取同樣重量的不同組織樣品進行混樣后，立即加入液氮磨樣，開始總RNA的提取。并將提取的RNA，分成３管送至中科新生命生物科技有限公司，基于第三代測序技術Pacbio[16]，采用單分子實時測序技術（SMRT）測序篩選從而得到較高質量高州油茶全長轉錄組數(shù)據(jù)。

1.2 NAC基因家族成員的鑒定

基于高質量高州油茶的全長轉錄組數(shù)據(jù)結果，通過cd–hit–v4.6.7軟件對下機處理好的高質量一致性序列進行去冗余，使用perl程序提其取蛋白序列，在TF數(shù)據(jù)庫（plantTFdb）中進行hmmscan比對，進行NAC成員鑒定，篩選獲得相關的NAC蛋白序列文件，利用PlantTFDB v5.0默認參數(shù)進行篩選分析最終獲得含保守結構域的NAC蛋白序列[17-18]。

1.3 NAC家族理化性質及序列比對分析

通過Expasy工具進行在線分析獲得NAC蛋白序列中所有相關理化性質進行分析處理。通過DNAMAN軟件中Multiple aligment sequence選項進行多序列比較，默認參數(shù)設置在Highlight homology level≥50%之間。

1.4 NAC蛋白進化樹分析

使用PlantTFDB v5.0數(shù)據(jù)庫，獲取擬南芥中關于NAC家族基因蛋白序列信息，在Ensembl Plants數(shù)據(jù)庫中挑選水稻數(shù)據(jù)庫獲得NAC家族基因的蛋白序列信息[19]。使用Mega11軟件將其匯總整理好的油茶、擬南芥以及水稻轉錄因子NAC蛋白序列文件進行系統(tǒng)進化樹構建[20]，先進行多序列比對，選用neighbor–joining算法，Bootstrap運行參數(shù)選擇

1 000，其余參數(shù)選擇默認值[21]。然后使用FigTreev 1.4.4軟件對構建好的進化樹進行優(yōu)化。

1.5 NAC基因家族的系統(tǒng)發(fā)育關系、保守基

序和基因結構分析

使用TBtools（v2.003）軟件拼接與繪制高州油茶51個NAC轉錄因子家族的系統(tǒng)發(fā)育關系、保守基序和基因結構圖[22]。使用MEME軟件，對高州油茶中NAC蛋白序列的保守基序進行分析[23-25]。

1.6 NAC轉錄因子亞細胞定位及跨膜分析

使用BUSCA亞細胞定位軟件[26-27]進行預測分析。通過TMHMM在線分析軟件，完成高州油茶NAC家族的跨膜分析[28]。

1.7 NAC蛋白二級結構預測分析

使用NPS@SOPMA軟件完成高州油茶NAC的蛋白質的預測與二級結構和分析[29-30]。

1.8 NAC蛋白磷酸化位點的分析

使用NetPhos 3.1 Server在線工具對51個高州油茶NAC家族中的蛋白氨基酸序列中可能存在的磷酸化位點進行統(tǒng)計分析[31]。

2 結果與分析

2.1 高州油茶NAC轉錄因子家族成員的鑒定

通過生物信息學篩選共獲得NAC轉錄因子基因家族51個成員（表1）。在ExPASy網(wǎng)站中預測結果表明，高州油茶NAC蛋白的氨基酸序列在250～756 aa之間，平均分子量為56 604.14 Da，其中CdNAC26氨基酸含量最多，其分子量為85 490.56 Da，而CdNAC51氨基酸數(shù)量的最少，其分子量只有27 611.84 Da。等電點范圍在4.65～10.13之間，其中CdNAC33最大，CdNAC15最小。NAC脂肪系數(shù)范圍在55.25～98.74之間，其中CdNAC48最大，而CdNAC21最小。經(jīng)過蛋白總平均疏水指數(shù)結果分析，得知高州油茶所有的NAC基因家族成員的平均疏水指數(shù)均為負值，表明其NAC蛋白為親水蛋白。

2.2 油茶NAC蛋白序列比對分析

通過對51個NAC蛋白多序列比對發(fā)現(xiàn)，NAC蛋白具有保守的5個亞結構域A～E[1]（圖1），但是

CdNAC3、CdNAC9、CdNAC13、CdNAC16、CdNAC25、

CdNAC27、CdNAC29和CdNAC43均缺少亞結構域A。5個結構域中C和D結構域相對比較穩(wěn)定，A、B、E結構域變異更多。

2.3 高州油茶NAC基因家族的系統(tǒng)發(fā)育關系、保守基序和基因結構分析

系統(tǒng)發(fā)育關系表明，116個NAC蛋白被劃分為10個亞族，其中CdNAC48和CdNAC49所在的亞族相對于其他的Ⅷ亞族最少，而Ⅹ亞族蛋白序列最多。對51個高州油茶蛋白序列進行保守基序和結構域分析（圖2A）發(fā)現(xiàn)，motif3和motif2主要構成了NAC蛋白的保守結構域（圖2B）。其中Ⅹ亞族全部都不含有motif2，CdNAC蛋白序列中保守基序最少的是motif7；其中有6種類型的motif是大多數(shù)的CdNAC蛋白序列中都含有的?；蚪Y構分析發(fā)現(xiàn)，NAM超家族結構域均存在于高州油茶的51個NAC蛋白序列中（圖2C）。

2.4 高州油茶、擬南芥和水稻NAC蛋白系統(tǒng)進化樹關系

為進一步研究NAC蛋白家族系統(tǒng)進化關系，對高州油茶51個，擬南芥36個和水稻的29個NAC蛋白序列進行匯總整理比對，同時采用鄰接法（neighbor–joining，NJ）方法構建系統(tǒng)進化樹[19-20]

（圖3）?？蓪⒏咧萦筒琛M南芥和水稻NAC蛋白劃分為10個亞族群，其中第Ⅹ號基因家族亞族群中未包含油茶蛋白只有OsNAC16和OsNAC17。OsNAC12

與其他蛋白亞群比較遠單獨形成一個分支。在Ⅲ亞族群中可能存在較遠的親緣關系CdNAC44、CdNAC45、CdNAC46、CdNAC47中未包含擬南芥和水稻蛋白。其余分支中既包含油茶蛋白序列，也包含擬南芥和水稻的蛋白序列。

2.5 油茶NAC轉錄因子亞細胞定位與跨膜預測分析

亞細胞定位分析表明，高州油茶中有23個NAC

蛋白均進行膜外運輸、跨膜運輸和膜內(nèi)運輸（表2），其中CdNAC11、CdNAC12、CdNAC13、CdNAC14位于細胞膜和細胞核，CdNAC22、CdNAC26位于細胞質和細胞核中，CdNAC48位于內(nèi)質網(wǎng)和細胞核中，其余蛋白均位于細胞核。

2.6 高州油茶NAC蛋白二級結構預測

表3是對高州油茶的NAC的蛋白序列二級結構進行預測結果，其中α–螺旋平均比例為25.36%，延伸鏈平均比例為19.30%，β–折疊平均比例為6.15%，

無規(guī)則卷曲平均比例為49.41%。其中無規(guī)卷曲占比最大，比例范圍達到29.49%～61.71%之間；β–折疊占比最小，比例范圍比較小只有3.16%～10.51%。

2.7 油茶NAC蛋白磷酸化位點分析

通過對51個高州油茶的蛋白序列磷酸化位點預測結果可知（表4），所有蛋白序列都包含絲氨酸（Ser）、蘇氨酸（Thr）和酪氨酸（Tyr）3個磷酸化位點，其中數(shù)目最多的是Ser位點，數(shù)目最少的是Tyr，CdNAC蛋白序列中磷酸化位點最多的是CdNAC19蛋白含有88個，最少的是CdNAC49蛋白只有23個。

3 討論與結論

基于高州油茶的三代全長轉錄組信息數(shù)據(jù)，共鑒定出51個高州油茶NAC家族成員，從蛋白理化性質分析來看，高州油茶NAC蛋白在氨基酸數(shù)、分子量、等電點、脂肪系數(shù)、總平均親水性、數(shù)量上存在較大的差異。高州油茶NAC蛋白的總平均疏水指數(shù)均為負數(shù)說明其蛋白多為親水性蛋白，這可能與該蛋白參與植物發(fā)育有關、基因表達調(diào)控、調(diào)控植物酶活、植物的信號傳遞等功能方面密切相關[32]。51個高州油茶NAC蛋白均具有NAM超家族結構域可能與NAC轉錄因子的DNA結合能力和轉錄活性有關[33]。51個高州油茶的蛋白結構域與擬南芥和水稻同時構建系統(tǒng)進化樹被分為10個亞群，其中Ⅲ亞簇群的CdNAC蛋白序列與其他兩個物種存在較大差異并未產(chǎn)生混合聚類，為高州油茶特異的NAC蛋白，后續(xù)可深入分析其在油茶中的特異功能。在51個油茶NAC蛋白中，有43個包含5個典型的亞結構域A～E，其中C和D在保守基序的組成上具有較高的相似性。

51個NAC蛋白均位于細胞核，其中還有23個蛋白參與了膜外與膜內(nèi)的跨膜運輸。蛋白結構預測結果中NAC蛋白最大的二級結構元件均為無規(guī)則卷曲結構，最小的為β–轉角結構，這與其他物種的NAC蛋白預測相似[19，28]。NAC轉錄因子是參與多種生物過程的一個家族，表明其作為植物育種的非生物脅迫耐受性具有巨大潛力[34]。從蛋白磷酸化位點預測結果中發(fā)現(xiàn)，Ser位點是高州油茶中蛋白中最多的磷酸化位點，這可能表明Ser磷酸化可能影響CdNAC蛋白的功能調(diào)控，可能參與植物細胞信號傳導相關[35]

油茶NAC轉錄因子在植物生長過程中具有重要的調(diào)控作用，但油茶NAC轉錄因子家族相關的研究比較少，研究首次利用高州油茶的全長轉錄組數(shù)據(jù)，利用生物信息學方法對高州油茶NAC轉錄因子家族進行生物信息分析，為后續(xù)相關基因功能的挖掘研究提供了參考。

參考文獻：

[1] 楊晨旭，哈斯巴根，杜超. 植物NAC轉錄因子結構和功能研究進展[J/OL]. 分子植物育種，1-16 [2024-07-29]. http：//kns.cnki.net/kcms/detail/46.1068.S.20230403.1809.023.html.

[2] 黎幫勇，胡尚連，曹穎，等. 毛竹NAC轉錄因子家族生物信息學分析[J]. 基因組學與應用生物學，2015，34（8）：1769-1777.

[3] RIECHMANN J L，HEARD J，MARTIN G，et al. Arabidopsis transcription factors：genome-wide comparative analysis among eukaryotes [J]. Science，2000，290（5499）：2105-2110.

[4] FANG Y J，YOU J，XIE K B，et al. Systematic sequence analysis and identification of tissue-specific or stress-responsive genes of NAC transcription factor family in rice [J]. Molecular Genetics and Genomics，2008，280（6）：547-563.

[5] 李圣龍，王傳銘，李曉靜. 石榴NAC轉錄因子家族的生物信息學分析[J]. 分子植物育種，2021，19（1）：88-99.

[6] DELESSERT C，KAZAN K，WILSON I W，et al. The transcription factor ATAF2 represses the expression of pathogenesis-related genes in Arabidopsis [J]. The Plant Journal，2005，43（5）：745-757.

[7] GUO Y F，GAN S S. AtNAP，a NAC family transcription factor，has an important role in leaf senescence [J]. The Plant Journal，2006，46（4）：601-612.

[8] SANJARI S，SHIRZADIAN-KHORRAMABAD R，SHOBBAR Z S，et al. Systematic analysis of NAC transcription factors' gene family and identification of post-flowering drought stress responsive members in sorghum [J]. Plant Cell Reports，2019，38（3）：361-376.

[9] 馬雪祺，陰艷紅，馮婧嫻，等. 植物NAC轉錄因子研究進展[J]. 植物生理學報，2021，57（12）：2225-2234.

[10] 董斌. 油茶抗旱性綜合評價及其干旱脅迫轉錄組研究[D]. 廣州：華南農(nóng)業(yè)大學，2018.

[11] WANG G D，ZHANG S，MA X C，et al. A stress-associated NAC transcription factor （SlNAC35）from tomato plays a positive role in biotic and abiotic stresses [J]. Physiologia Plantarum，2016，158（1）：45-64.

[12] 彭國穎，萬瑋，黃超，等. 紫鴨跖草NAC轉錄因子家族生物信息學分析[J]. 分子植物育種，2024，22（6）：1901-1909 .

[13] 吳軍，察艷艷，李霞，等.大豆中4個NIN轉錄因子的生物信息學及表達分析[J].華北農(nóng)學報，2022，37（S1）：26-34.

[14] 黎運. 油茶對干旱脅迫的生理響應及轉錄組分析[D]. 長沙：中南林業(yè)科技大學，2021.

[15] 陳永忠，羅健，陳隆升，等. 湖南省油茶旱害及防治對策[J]. 經(jīng)濟林研究，2014，32（3）：22-29，53.

[16] RHOADS A，AU K F. PacBio sequencing and its applications [J]. Genomics，Proteomics amp; Bioinformatics，2015，13（5）：278-289.

[17] GUO A Y，CHEN X，GAO G，et al. PlantTFDB：a comprehensive

plant transcription factor database [J]. Nucleic Acids Research，2008，36（Database issue）：D966-D969.

[18] POTTER S C，LUCIANI A，EDDY S R，et al. HMMER web server：2018 update [J]. Nucleic Acids Res，2018，46（w1）：W200-W204.

[19] 黎幫勇，胡尚連，曹穎，等. 毛竹NAC轉錄因子家族生物信息學分析[J]. 基因組學與應用生物學，2015，34（8）：1769-1777.

[20] 陳云霞，丁雨龍. 毛竹中NAC轉錄因子的克隆與序列分析[J]. 安徽農(nóng)業(yè)大學學報，2013，40（3）：372-377.

[21] SAITOU N，NEI M. The neighbor-joining method：a new method for reconstructing phylogenetic trees [J]. Molecular Biology and Evolution，1987，4（4）：406-425.

[22] CHEN C J，CHEN H，ZHANG Y，et al. TBtools：an integrative toolkit developed for interactive analyses of big biological data [J]. Molecular Plant，2020，13（8）：1194-1202.

[23] 郭怡婷，趙文菊，任延靖，等. 菊芋NAC轉錄因子家族基因的鑒定及分析[J]. 生物技術通報，2023，39（6）：217-232.

[24] 郭績寧. 亞麻薺 NAC 轉錄因子家族的生物信息學分析[D]. 唐山：華北理工大學，2022.

[25] 廖桂花，段鈺，王叢叢，等. 赪桐bHLH轉錄因子家族鑒定及生物信息學分析[J/OL]. 分子植物育種，1-17 [2024-07-26]. http：//kns.cnki.net/kcms/detail/46.1068.S.20230529.1303.022.html.

[26] CHOU K C，SHEN H B. Cell-PLoc：a package of Web servers for predicting subcellular localization of proteins in various organisms [J].

Nature Protocols，2008，3（2）：153-162.

[27] CHOU K C，SHEN H B. A new method for predicting the subcellular localization of eukaryotic proteins with both single and multiple sites：Euk-mPLoc 2.0 [J]. PLoS One，2010，5（4）：e9931.

[28] 段鈺，廖桂花，王叢叢，等. 赪桐NAC轉錄因子家族鑒定及生物信息學分析[J/OL]. 分子植物育種，1-18 [2024-07-26]. http：//kns.cnki.net/kcms/detail/46.1068.S.20230529.1302.020.html.

[29] GEOURJON C，DELéAGE G. SOPMA：significant improvements in protein secondary structure prediction by consensus prediction from multiple alignments [J]. Computer Applications in the Biosciences，1995，11（6）：681-684.

[30] 孫欣，上官凌飛，房經(jīng)貴，等. 葡萄NAC轉錄因子家族生物信息學分析[J]. 基因組學與應用生物學，2011，30（2）：229-242.

[31] 方志紅，劉建寧，張燕，等. 白羊草NAC轉錄因子家族生物信息學分析[J]. 草地學報，2020，28（5）：1266-1274.

[32] 陳倩，游雙梅，邢樂華，等. 果樹NAC轉錄因子的研究進展[J]. 分子植物育種，2021，19（19）：6396-6405.

[33] XU W J，DUBOS C，LEPINIEC L. Transcriptional control of flavonoid biosynthesis by MYB–bHLH–WDR complexes [J]. Trends in Plant Science，2015，20（3）：176-185.

[34] SUN H，HU M L，LI J Y，et al. Comprehensive analysis of NAC transcription factors uncovers their roles during fiber development and stress response in cotton [J]. BMC Plant Biology，2018，18（1）：150.

[35] HAN F X，WANG P Y，CHEN X，et al. An ethylene-induced NAC transcription factor acts as a multiple abiotic stress responsor in conifer [J]. Horticulture Research，2023，10（8）：uhad130.

（責任編輯：肖彥資）