番茄組蛋白乙酰轉(zhuǎn)移酶(HAT)的全基因組鑒定與分析

史建磊 熊自立 李濤 張海利 宰文珊

摘要：組蛋白乙?；揎椩诨虮磉_(dá)調(diào)控方面扮演著重要角色。為深入挖掘番茄組蛋白乙酰轉(zhuǎn)移酶（HAT）基因，本研究運(yùn)用生物信息學(xué)方法，共鑒定出26個(gè)番茄HAT。聚類分析發(fā)現(xiàn)，這些HAT可以分為7組，分別為HAG、HAG1、HAG2、MCC1、HAM、HAF和HAC。其中，GNAT家族（GCN5、HAT1、MCC1） 18個(gè)、p300/CBP家族4個(gè)、MYST和TAFII250家族各1個(gè)。6對同源基因Ka/Ks值均小于1，說明其進(jìn)化中經(jīng)歷了純化選擇。不同家族HAT理化性質(zhì)存在差異，但均為親水蛋白質(zhì)，主要定位于細(xì)胞核，蛋白質(zhì)二級結(jié)構(gòu)主要包括α螺旋和無規(guī)則卷曲，具有典型的保守基序和特征結(jié)構(gòu)域。染色體定位發(fā)現(xiàn)，番茄HAT基因不均勻地散布于12條染色體上，且呈兩端分布，其中1號和9號染色體最多，部分基因經(jīng)串聯(lián)重復(fù)形成基因簇。基因表達(dá)分析發(fā)現(xiàn)，GNAT家族番茄特有的HAT和HAG2組基因在病菌、鹽和熱脅迫下具有較高的表達(dá)水平，一些基因具有病菌誘導(dǎo)特異性和耐鹽（熱）與鹽（熱）敏材料上表達(dá)差異性?？梢?，番茄HAT在結(jié)構(gòu)、性質(zhì)和功能上具有多樣性，其在全基因組水平上的鑒定與分析為相關(guān)基因克隆利用提供了依據(jù)。

關(guān)鍵詞：番茄;組蛋白乙酰轉(zhuǎn)移酶;生物信息學(xué)

中圖分類號：Q786文獻(xiàn)標(biāo)識碼：A文章編號：1000-4440（2020）03-0666-09

Genome-wide identification and analysis of histone acetyltransferase （HAT） in tomato

SHI Jian-lei1，XIONG Zi-li1，LI Tao2，ZHANG Hai-li1，ZAI Wen-shan1

（1.Wenzhou Vocational College of Science and Technology， Wenzhou 325006， China;2.Vegetable Research Institute， Guangdong Academy of Agricultural Sciences， Guangzhou 510640， China）

Abstract： Histone acetylation plays an important role in the regulation of gene expression. In order to further explore histone acetyltransferase （HAT） gene in tomato， 26 HAT genes were identified by using bioinformatics method in this study. The results of cluster analysis showed that these HAT could be divided into seven groups， namely HAG， HAG1， HAG2， MCC1， HAM， HAF and HAC. Among them， 18 HAT genes were GNAT family members （GCN5， HAT1， MCC1）， four HAT genes were p300/CBP family members， and one was MYST or TAFII250 family member. The Ka/Ks values of six pairs of paralogous genes were all less than one， indicating that purifying selection had occurred in evolution. The physicochemical properties of HAT in different families were different，but all of them were hydrophilic proteins， mainly located in the nucleus. Furthermore， the secondary structure of proteins mainly included alpha helix and random coil， and had typical conserved motifs and domains. Chromosome localization results showed that tomato HAT gene were unevenly distributed on 12 chromosomes and at both ends. In addition， the number of chromosome 1 and 9 was the most， and some genes formed gene clusters by tandem duplication. Results of gene expression analysis indicated that the HAT genes， which were unique in tomato， and HAG2 gene， had higher expression levels under pathogen， salt and heat stress. Also， some genes had pathogen-inducible specificity and different expression levels between resistant and sensitive materials. It can be seen that HAT gene in tomato? has diversity in structure， character and function， and the genome-wide identification and analysis results providea basis for cloning and utilization of related genes.

Key words：tomato;histone acetyltransferase （HAT）;bioinformatics

植物在生長發(fā)育過程中面臨著各種生物與非生物脅迫，同時(shí)，在長期歷史進(jìn)化中植物本身也形成了多種調(diào)控與防御機(jī)制，來保障個(gè)體正常生長發(fā)育與群體繁衍。組蛋白修飾在基因表達(dá)調(diào)控方面扮演著重要角色，通過組蛋白乙酰轉(zhuǎn)移酶（HAT）促進(jìn)基因表達(dá)和組蛋白去乙?；福℉DAC）抑制基因表達(dá)。組蛋白乙?；且环N重要的翻譯后修飾，通常與基因的轉(zhuǎn)錄激活有關(guān)，根據(jù)HAT在細(xì)胞內(nèi)分布及誘導(dǎo)乙?；蟮男?yīng)分為A型（核內(nèi)）和B型（胞質(zhì)）兩類[1]。目前已知的A型HAT種類較多，如GNAT （Gcn5相關(guān)N-末端乙酰轉(zhuǎn)移酶）、MYST （MOZ，Ybf2/Sas3，Sas2和Tip60）、p300/CBP （cAMP響應(yīng)原件結(jié)合蛋白）、TAFII250 （TATA結(jié)合蛋白相關(guān)因子）等[2-4]。

前人已對擬南芥（Arabidopsis thaliana） [5]、水稻（Oryza sativa） [6]和葡萄（Vitis vinifera） [7]組蛋白乙?；蚣易暹M(jìn)行了生物信息和基因功能研究。組蛋白乙?；腿ヒ阴；揎棾齾⑴c植物生長發(fā)育外，也參與多種生物與非生物脅迫響應(yīng)[8-10]。番茄（Solanum lycopersicum）是世界上重要的園藝和經(jīng)濟(jì)作物，也是重要的科研模式作物，盡管Cigliano等[11] 2013年在番茄中發(fā)現(xiàn)了32個(gè)HAT，但與HDAC相比，相關(guān)研究仍非常有限。同時(shí)，番茄全基因組測序的完成[12]，為在全基因組水平上鑒定番茄HAT提供了條件。鑒于HAT基因家族的重要生理功能，本研究利用生物信息學(xué)方法對番茄HAT進(jìn)行鑒定與分析，旨在為其進(jìn)一步功能研究和分子育種應(yīng)用提供基礎(chǔ)信息。

1材料與方法

1.1番茄組蛋白乙酰轉(zhuǎn)移酶（HAT）全基因組鑒定

番茄全基因組序列及注釋信息（ITAG3.2 release）下載自茄科基因組網(wǎng)站SGN （https：//solgenomics.net/），并利用BioEdit 7.2.5構(gòu)建本地?cái)?shù)據(jù)庫。采用2種方法鑒定番茄組蛋白乙酰轉(zhuǎn)移酶（HAT）：從Pfam 32.0 （http：//pfam.xfam.org/）數(shù)據(jù)庫中下載HAT隱馬爾科夫模型（Hidden Markov Model，HMM）特征序列，結(jié)合擬南芥HAT，作為問詢序列Blastp （E-value<1）構(gòu)建本地?cái)?shù)據(jù)庫;同時(shí)以“histone acetyltransferase”為關(guān)鍵詞在SGN中進(jìn)行檢索，去除冗余基因，即獲得候選基因。然后通過Pfam和SMART （http：//smart.embl-heidelberg.de/smart/batch.pl）鑒定HAT特征結(jié)構(gòu)域，將不含目的結(jié)構(gòu)域的候選基因排除。

1.2番茄組蛋白乙酰轉(zhuǎn)移酶（HAT）系統(tǒng)樹的構(gòu)建

提取番茄HAT氨基酸序列，使用ClustalX 1.83進(jìn)行序列聯(lián)配，然后利用MEGA 7.0的鄰接法（Neighbor-Joining，NJ） （參數(shù)設(shè)置：Bootstrap method 1 000，p-distance and Pairwise deletion）進(jìn)行系統(tǒng)發(fā)育樹構(gòu)建，去除Bootstrap支持率低于50%的節(jié)點(diǎn)，顯示各分支。

1.3番茄組蛋白乙酰轉(zhuǎn)移酶（HAT）理化參數(shù)分析

利用在線工具ExPASy的ProtParam程序（https：//web.expasy.org/protparam/）進(jìn)行蛋白質(zhì)理化參數(shù)預(yù)測，SOPMA （https：//npsa-prabi.ibcp.fr/cgi-bin/npsa_automat.pl？page=npsa_sopma.html）預(yù)測蛋白質(zhì)二級結(jié)構(gòu)，利用SWISS-MODEL （https：//swissmodel.expasy.org/）進(jìn)行蛋白質(zhì)三維結(jié)構(gòu)同源建模，利用WoLF PSORT （https：//wolfpsort.hgc.jp/）進(jìn)行蛋白質(zhì)亞細(xì)胞定位預(yù)測。

1.4番茄組蛋白乙酰轉(zhuǎn)移酶（HAT）結(jié)構(gòu)特征分析

利用在線工具GSDS 2.0 （http：//gsds.cbi.pku.edu.cn/）繪制基因的外顯子-內(nèi)含子結(jié)構(gòu)模式圖。利用數(shù)據(jù)庫Pfam 32.0和在線工具M(jìn)EME 5.1.0 （http：//meme-suite.org/tools/meme）分別分析蛋白質(zhì)保守結(jié)構(gòu)域和基序（Motifs），并使用TBtools軟件[13]對其進(jìn)行可視化。

1.5番茄組蛋白乙酰轉(zhuǎn)移酶（HAT）的染色體定位

結(jié)合SGN中基因位置信息（Tomato SL3.0），利用MapDraw 2.1對番茄HAT進(jìn)行染色體定位，繪制其染色體分布圖。同時(shí)進(jìn)行基因簇和串聯(lián)重復(fù)分析，其中劃分基因簇的原則為：2個(gè)相鄰目的基因間距小于200 kb且其他基因不得多于8個(gè);判斷串聯(lián)重復(fù)的原則為：相鄰目的基因間距小于100 kb且序列相似度大于70% [14]。用KaKs_Calculator 2.0軟件[15]計(jì)算基因間堿基替換率（Ka/Ks）。

1.6番茄組蛋白乙酰轉(zhuǎn)移酶（HAT）基因表達(dá)分析

基于已發(fā)表的Affymetrix番茄基因組和TOM2芯片數(shù)據(jù)，下載8張芯片表達(dá)結(jié)果，分別為E023 （真菌激發(fā)子）、E024 （番茄細(xì)菌性潰瘍病菌）、E031 （灰霉菌）、E048 （馬鈴薯紡錘塊莖類病毒）、E046 （番茄斑萎病毒）、E047 （叢枝菌根真菌球霉菌）、E051 （鹽脅迫）和E053 （干旱和熱脅迫），利用MeV 4.9.0繪制基因表達(dá)譜。

2結(jié)果與分析

2.1番茄組蛋白乙酰轉(zhuǎn)移酶（HAT）鑒定及系統(tǒng)聚類

通過Blastp和關(guān)鍵詞檢索，共鑒定出26個(gè)番茄組蛋白乙酰轉(zhuǎn)移酶（HAT）。經(jīng)多重序列比對后利用鄰接法構(gòu)建系統(tǒng)發(fā)生樹（圖1）。38個(gè)番茄和擬南芥HAT可聚為8組（Ⅰ～Ⅷ）。根據(jù)擬南芥HAT分類標(biāo)準(zhǔn)，第Ⅰ組8個(gè)HAT為番茄所特有;第Ⅱ組1個(gè)番茄HAT，且與擬南芥AtHAG2同源;第Ⅳ組1個(gè)番茄HAT，且與AtHAG1同源;第Ⅵ組4個(gè)番茄HAT，且與擬南芥AtMCC1同源，另有4個(gè)AtMCC1同源HAT散落于樹中，這8個(gè)HAT屬于GNAT家族;第Ⅲ組1個(gè)番茄HAT，且與AtHAM同源，屬于MYST家族;第Ⅴ組1個(gè)番茄HAT，且與AtHAF同源，屬于TAFII250家族;第Ⅷ組4個(gè)HAT，且與AtHAC同源，屬于p300/CBP家族。整體來看，GNAT家族的18個(gè)HAT大部分（77.78%）聚在第Ⅰ、Ⅱ、Ⅳ和Ⅵ組，但仍有少數(shù)HAT散落于其他組。同時(shí)，番茄中未發(fā)現(xiàn)與擬南芥AtHAG3同源的基因（Ⅶ）。

結(jié)合序列相似度比較，發(fā)現(xiàn)6對（13個(gè)）同源基因（Solyc01g008120.3.1和Solyc04g008610.3.1，Solyc02g089790.3.1和Solyc05g006180.2.1，Solyc06g065570.3.1和Solyc08g067440.2.1，Solyc09g082240.3.1和Solyc09g082250.2.1，Solyc09g072580.3.1和Solyc09g090120.3.1，Solyc01g097990.3.1、Solyc01g098010.3.1和Solyc01g098020.2.1），且對應(yīng)基因間Ka/Ks值均小于1，說明其在自然進(jìn)化過程中經(jīng)歷了純化選擇。

2.2番茄組蛋白乙酰轉(zhuǎn)移酶（HAT）的理化特征分析

從表1可見，26個(gè)番茄HAT氨基酸數(shù)為80～1 856，相對分子質(zhì)量為9 200～210 630 ，其中TAFII250和p300/CBP家族成員所含氨基酸較多、相對分子質(zhì)量較大;理論等電點(diǎn)為4.62～9.58，平均7.40，整體接近中性;TAFII250和p300/CBP為不穩(wěn)定蛋白質(zhì)，而多數(shù)GNAT和MYST家族成員為穩(wěn)定蛋白質(zhì);所有鑒定成員均為親水性蛋白質(zhì)，主要定位于細(xì)胞核，其次定位于細(xì)胞質(zhì)和葉綠體;蛋白質(zhì)二級結(jié)構(gòu)主要包括α螺旋和無規(guī)則卷曲。結(jié)合系統(tǒng)聚類結(jié)果（圖1），選取有代表性的8個(gè)基因（G1 Solyc01g098020.2.1， G2 Solyc08g067440.2.1， G3 Solyc11g013520.2.1， G4 Solyc10g045390.2.1， G5 Solyc02g068580.1.1， G6 Solyc07g006820.3.1， G7 Solyc09g072580.3.1和G8 Solyc05g006180.2.1）編碼的蛋白質(zhì)進(jìn)行蛋白質(zhì)三維建模（圖2），發(fā)現(xiàn)GNAT （G1、G2、G4、G5）和MYST （G3）家族除了α螺旋，還存在β折疊，但TAFII250 （G6）和p300/CBP （G8）家族無β折疊。同時(shí)，發(fā)現(xiàn)部分GNAT和MYST家族成員具有活性位點(diǎn)。

基因結(jié)構(gòu)分析發(fā)現(xiàn)，GNAT家族基因整體較小、外顯子數(shù)較少，Solyc10g045390.2.1基因較大但氨基酸數(shù)較少，推測主要由非編碼區(qū)組成;而MYST、TAFII250和p300/CBP家族基因較大、外顯子數(shù)較多。GNAT家族基因主要是1/2位內(nèi)含子，MYST是0/1位內(nèi)含子，TAFII250和p300/CBP是0/1/2位內(nèi)含子，同時(shí)，部分基因的部分外顯子為非對稱外顯子（圖3）?；蚍治霭l(fā)現(xiàn)，20個(gè)蛋白質(zhì)保守基序大小為21～50，不同家族HAT所含基序種類不同，如GNAT組1主要含Motif 1、Motif 7、Motif 8和Motif 17，GNAT組2主要含Motif 14、Motif 18和Motif 19;MYST和TAFII250家族與GNAT組2相似;p300/CBP家族主要含Motif2～Motif6、Motif 9～Motif 13、Motif 15～Motif 16和Motif 20，同一家族成員存在基序丟失情況。另外，不同家族HAT具有典型的特征結(jié)構(gòu)域，如GNAT家族的Acetyltransf_1和Hat1_N，MYST家族的zf-MYST和MOZ_SAS，p300/CBP家族的HAT_KAT11和zf-TAZ。

G1～G8為8個(gè)代表性HAT。

2.3番茄組蛋白乙酰轉(zhuǎn)移酶基因（HAT）染色體定位

26個(gè)番茄組蛋白乙酰轉(zhuǎn)移酶基因（HAT）不均勻地散布于番茄12條染色體上，且存在于兩端，其中9號染色體分布有7個(gè)HAT，1號染色體為5個(gè)HAT，其他染色體為1～3個(gè)HAT （圖4）。其中，9號染色體和1號染色體分別有4個(gè)和3個(gè)基因形成單基因簇。聚類分析結(jié)合序列相似度比較獲得的8對同源基因中，3對位于不同染色體上，可能是經(jīng)片段重復(fù)而來;4對分別位于同一染色體基因簇內(nèi)，為串聯(lián)重復(fù)而來。

實(shí)線方框表示基因簇，虛線方框表示串聯(lián)重復(fù)，虛線連接表示片段重復(fù)。

2.4番茄組蛋白乙酰轉(zhuǎn)移酶基因（HAT）表達(dá)譜分析

番茄生物和非生物脅迫（鹽、干旱和熱）芯片數(shù)據(jù)表達(dá)譜見圖5 （部分基因未篩選到相應(yīng)探針）。由圖5A和5B可知，GNAT家族番茄特有HAT （Ⅰ）和AtHAG2同源基因Solyc02g068580.1.1（Ⅱ）在所有病菌侵染下，都具有較高的表達(dá)水平，部分AtMCC1同源基因在灰霉菌和番茄斑萎病毒（TSWV）侵染下具有較高的表達(dá)水平，p300/CBP家族番茄HAT （Ⅷ）在TSWV、叢枝菌根真菌球霉菌（MYC）和馬鈴薯紡錘塊莖類病毒（PSTVd）侵染下具有較高的表達(dá)水平，MYST家族的Solyc11g013520.2.1 （Ⅲ）在灰霉菌和PSTVd侵染下具有較高的表達(dá)水平，TAFII250家族的Solyc07g006820.3.1 （Ⅴ）在TSWV和MYC侵染下具有較高的表達(dá)水平，其他AtMCC1同源基因和AtHAG1同源基因Solyc10g045390.2.1 （Ⅳ）表達(dá)水平較低。

由圖5C可知，GNAT家族番茄特有HAT （Ⅰ）在鹽和熱脅迫下均有較高的表達(dá)水平，但耐鹽（熱）和鹽（熱）敏材料區(qū)分不明顯;p300/CBP家族的Solyc01g008120.3.1、Solyc02g089790.3.1 （Ⅷ）和AtMCC1同源基因Solyc05g025890.2.1在干旱脅迫下具有較高的表達(dá)水平，同樣與對照材料區(qū)分不明顯，但后兩者能夠區(qū)分耐熱和熱敏材料;AtMCC1同源基因Solyc03g097140.3.1 （Ⅵ）在熱脅迫下具有較高的表達(dá)水平，與對照材料區(qū)分不明顯，但可以區(qū)分耐鹽/干旱和不耐鹽/干旱材料;除AtMCC1同源基因Solyc06g065570.3.1、Solyc08g067440.2.1 （Ⅵ，鹽脅迫下表達(dá)水平較高）和Solyc05g005900.3 （表達(dá)水平較低）外，其他基因在至少一種非生物脅迫下與對照材料間表達(dá)差異明顯。

A和B：番茄HAT在6種生物脅迫下的基因表達(dá);C：番茄HAT在3種非生物脅迫下的基因表達(dá)。a1：真菌激發(fā)子誘導(dǎo)4 h;a2：真菌激發(fā)子誘導(dǎo)8 h;a3：番茄細(xì)菌性潰瘍病菌接種4 d;a4：番茄細(xì)菌性潰瘍病菌接種8 d;a5：接種灰霉菌的成熟綠果;a6：接種灰霉菌的成熟紅果;a7：接種馬鈴薯紡錘塊莖類病毒的栽培種Rutgers;a8：接種馬鈴薯紡錘塊莖類病毒的栽培種Moneymaker。b1：番茄斑萎病毒侵染51 d的根;b2：番茄斑萎病毒侵染51 d的葉;b3：叢枝菌根真菌球霉菌處理51 d的根;b4：叢枝菌根真菌球霉菌處理51 d的葉。c1：鹽脅迫下的栽培種Moneymaker;c2：鹽脅迫下的野生醋栗番茄PI365967;c3干旱脅迫下的敏感材料;c4：干旱脅迫下的耐旱材料;c5：熱脅迫下的敏感材料;c6：熱脅迫的耐熱材料。

3討論

植物表觀遺傳學(xué)是當(dāng)今研究熱點(diǎn)之一，而組蛋白修飾調(diào)控是表觀遺傳學(xué)的重要研究領(lǐng)域。組蛋白乙?；鳛榻M蛋白修飾的一種，由組蛋白乙酰化和去乙?；餐瓿蒣16]。在真核生物中，組蛋白乙酰化是染色質(zhì)中的可逆生物學(xué)過程，是促進(jìn)染色質(zhì)松弛和基因轉(zhuǎn)錄調(diào)控的主要因素之一[17]。

3.1番茄組蛋白乙酰轉(zhuǎn)移酶（HAT）鑒定與聚類

本研究結(jié)合Blastp和關(guān)鍵詞檢索，共鑒定出26個(gè)番茄組蛋白乙酰轉(zhuǎn)移酶（HAT）。通過與擬南芥HAT共同構(gòu)建系統(tǒng)發(fā)生樹，可將番茄HAT分為7組，第Ⅰ組為番茄所特有，第Ⅱ組與擬南芥AtHAG2同源，第Ⅳ組與AtHAG1同源，第Ⅵ組與AtMCC1同源，這4組屬于GNAT家族;第Ⅲ組與AtHAM同源，屬于MYST家族;第Ⅴ組與AtHAF同源，屬于TAFII250家族;第Ⅷ組與AtHAC同源，屬于p300/CBP家族。GNAT家族擴(kuò)張程度較大，散布于多個(gè)組中，但番茄中未發(fā)現(xiàn)AtHAG3同源基因。所建發(fā)生樹包含基因染色體來源不同的混合進(jìn)化枝，反映出其在染色體上的擴(kuò)張情況，可能是產(chǎn)生具有新功能基因的源泉[18]。同一進(jìn)化枝包含不同物種基因，說明其可能按照物種特異性的方式進(jìn)行了擴(kuò)張。A型HAT僅存在于核內(nèi)，并與基因轉(zhuǎn)錄密切相關(guān);B型HAT通常只能使新合成的H3和H4組蛋白N端尾部特定賴氨酸殘基乙酰化，并影響其在核小體組裝中的定位[1]。本研究發(fā)現(xiàn)多數(shù)番茄HAT為A型，但也有B型存在，如第Ⅱ組的AtHAG2同源基因Solyc02g068580.1.1。

3.2番茄組蛋白乙酰轉(zhuǎn)移酶（HAT）的理化特征

26個(gè)番茄HAT中，8個(gè)GNAT家族HAT為番茄所特有，其他3個(gè)家族HAT均比擬南芥要少[5]，也比Cigliano等[11]發(fā)現(xiàn)的番茄HAT要少，可能與物種本身、所用數(shù)據(jù)庫和鑒定標(biāo)準(zhǔn)有關(guān)，今后需進(jìn)一步完善數(shù)據(jù)庫和標(biāo)準(zhǔn)，使鑒定種類和數(shù)量更科學(xué)。這26個(gè)HAT基因大小及其編碼的蛋白質(zhì)穩(wěn)定性等存在明顯差異，且具有家族特異性，但均為親水蛋白質(zhì)，整體接近中性，主要定位于細(xì)胞核，包括α螺旋和無規(guī)則卷曲，GNAT/MYST家族成員還存在β折疊?；蚪Y(jié)構(gòu)分析發(fā)現(xiàn)不同家族基因內(nèi)含子相位不同，可能發(fā)生了內(nèi)含子單獨(dú)獲取或丟失[19]，且部分基因的部分外顯子為非對稱外顯子，容易造成閱讀框的推移，可能是造成基因多樣性的原因之一。不同家族蛋白質(zhì)具有典型的保守基序和特征結(jié)構(gòu)域，是彼此區(qū)分的重要依據(jù)。

3.3番茄組蛋白乙酰轉(zhuǎn)移酶基因（HAT）染色體定位

26個(gè)番茄組蛋白乙酰轉(zhuǎn)移酶基因（HAT）不均勻地散布于番茄12條染色體上，且主要（46.15%）位于1號和9號染色體上，呈兩端分布，這種基因在進(jìn)化過程中可能更易變異?；虼睾痛?lián)重復(fù)分析發(fā)現(xiàn)，1號和9號染色體上各有一個(gè)基因簇，簇內(nèi)分別有3個(gè)和2個(gè)基因形成串聯(lián)重復(fù)，說明串聯(lián)重復(fù)是基因擴(kuò)張和形成基因簇的重要機(jī)制[20]，在其他家族基因中同樣如此[19，21]?；蛲葱苑治霭l(fā)現(xiàn)，有3對基因來自不同染色體，可能經(jīng)片段重復(fù)而來，其中1對基因在染色體上的排列順序相反，暗示這些區(qū)段可能發(fā)生了倒位。同時(shí)，同源基因間Ka/Ks均小于1，說明其在自然進(jìn)化過程中經(jīng)歷了純化選擇。

3.4番茄組蛋白乙酰轉(zhuǎn)移酶基因（HAT）表達(dá)

通過番茄生物和非生物脅迫（鹽、干旱和熱）芯片表達(dá)數(shù)據(jù)分析，發(fā)現(xiàn)GNAT家族番茄特有HAT和AtHAG2同源基因Solyc02g068580.1.1在所有病菌侵染下，都具有較高的表達(dá)水平，部分AtMCC1同源基因和AtHAG1同源基因Solyc10g045390.2.1表達(dá)水平較低，其他家族基因具有病菌誘導(dǎo)特異性。GNAT家族番茄特有HAT在鹽和熱脅迫下均有較高的表達(dá)水平，耐鹽（熱）和鹽（熱）敏材料區(qū)分不明顯，但仍有超過50%的基因在至少一種非生物脅迫下與對照材料間存在差異表達(dá)。同樣，GNAT家族基因、尤其是AtMCC1同源基因具有表達(dá)多樣性，當(dāng)然與其數(shù)量較多有關(guān)，也從側(cè)面反映出其家族特性?；蚋弑磉_(dá)雖能響應(yīng)不同生物與非生物脅迫，但同時(shí)增加了植物負(fù)擔(dān)，植物通過假基因化、丟失和相關(guān)表達(dá)調(diào)控機(jī)制進(jìn)行自我保護(hù)，維持生命活動物質(zhì)和能量代謝的平衡。

4結(jié)論

本研究運(yùn)用生物信息學(xué)方法，共鑒定出26個(gè)番茄組蛋白乙酰轉(zhuǎn)移酶（HAT），根據(jù)序列特征可以分為7組，分別為HAG、HAG1、HAG2、MCC1、HAM、HAF和HAC。其中，GNAT家族（GCN5、HAT1、MCC1） 18個(gè)，p300/CBP家族4個(gè)，MYST和TAFII250家族各1個(gè)。同時(shí)，發(fā)現(xiàn)8對同源基因且經(jīng)歷了純化選擇。這些HAT均為親水蛋白質(zhì)，整體接近中性，主要定位于細(xì)胞核，包括α螺旋和無規(guī)則卷曲。不同家族酶具有典型的保守基序和特征結(jié)構(gòu)域。染色體定位發(fā)現(xiàn)，番茄HAT不均勻地散布于12條染色體上，且呈兩端分布，部分基因經(jīng)串聯(lián)重復(fù)形成基因簇。GNAT家族番茄特有HAT和HAG2組基因在病菌、鹽和熱脅迫下具有較高的表達(dá)水平，一些基因具有病菌誘導(dǎo)特異性和耐鹽（熱）、鹽（熱）敏材料表達(dá)差異性。

參考文獻(xiàn)：

[1]MIZZEN C A，ALLIS C D. Linking histone acetylation to transcriptional regulation [J]. Cell Metabolism Life Science，1998，54： 6-20.

[2]STERNER D E，BERGER S L. Acetylation of histones and transcription-related factors [J]. Microbiology and Molecular Biology Reviews，2000，64（2）： 435-459.

[3]MIZZEN C A，YANG X J，KOKUBO T，et al. The TAFII250 subunit of TFIID has histone acetyltransferase activity [J]. Cell，1996，87： 1261-1270.

[4]BANNISTER A J，KOUZARIDES T. The CBP co-activator is a histone acetyltransferase [J]. Nature，1996，384： 641-643.

[5]PANDEY R，MULLER A，NAPOLI C A，et al. Analysis of histone acetyltransferase and histone deacetylase families of Arabidopsis thaliana suggests functional diversification of chromatin modification among multicellular eukaryotes [J]. Nucleic Acids Research，2002，30（23）： 5036-5055.

[6]LIU X，LUO M，ZHANG W，et al. Histone acetyltransferases in rice （Oryza sativa L.）： phylogenetic analysis，subcellular localization and expression [J]. BMC Plant Biology，2012，12： 145.

[7]AQUEA F，TIMMERMANN T，ARCE-JOHNSON P. Analysis of histone acetyltransferase and deacetylase families of Vitis vinifera [J]. Plant Physiology and Biochemistry，2010，48（2/3）： 194-199.

[8]KORNET N，SCHERES B. Members of the GCN5 histone acetyltransferase complex regulate PLETHORA-mediated root stem cell niche maintenance and transit amplifying cell proliferation in Arabidopsis [J]. The Plant Cell，2009，21（4）：1070-1079.

[9]LATRASSE D，BENHAMED M，HENRY Y，et al. The MYST histone acetyltransferases are essential for gametophyte development in Arabidopsis [J]. BMC Plant Biology，2008，8： 121.

[10]HAN S K，SONG J D，NOH Y S，et al. Role of plant CBP/p300-like genes in the regulation of flowering time [J]. The Plant Journal，2007，49（1）：103-114.

[11]CIGLIANO R A，SANSEVERINO W，CREMONA G，et al. Genome-wide analysis of histone modifiers in tomato： gaining an insight into their developmental roles [J]. BMC Genomics，2013，14： 57.

[12]SATO S，TABATA S，HIRAKAWA H，et al. The tomato genome sequence provides insights into fleshy fruit evolution [J]. Nature，2012，485（7400）： 635-641.

[13]CHEN C J，XIA R，CHEN H，et al. TBtools，a Toolkit for Biologists integrating various HTS-data handling tools with a user-friendly interface [J]. BioRxiv，2018，3. doi： https：//doi.org/10.1101/289660.

[14]HUANG S，GAO Y，LIU J，et al. Genome-wide analysis of WRKY transcription factors in Solanum lycopersicum [J]. Molecular Genetics and Genomics，2012，287（6）： 495-513.

[15]ZHANG Z，LI J，ZHAO X Q，et al. KaKs_calculator： calculating Ka and Ks through model selection and model averaging [J]. Genomics Proteomics Bioinformatics，2006，4： 259-263.

[16]李濤，蘇慧慧，李植良，等.番茄HDACs家族基因在脅迫條件下的表達(dá)分析[J].熱帶作物學(xué)報(bào)，2015，36（11）：1994-2001.

[17]WATERBORG J H. Plant histone acetylation： in the beginning [J]. Biochimica et Biophysica Acta，2011，1809： 353-359.

[18]劉云飛，萬紅建，韋艷萍，等.番茄NBS-LRR抗病基因家族全基因組分析[J].核農(nóng)學(xué)報(bào)，2014，28（5）：790-799.

[19]劉潮，韓利紅，宋培兵，等.辣椒類甜蛋白基因家族鑒定及表達(dá)分析[J].江蘇農(nóng)業(yè)學(xué)報(bào)，2018，34（1）：122-129.

[20]YANG Z F，ZHOU Y，WANG X F，et al. Genome-wide comparative phylogenetic and molecular evolutionary analysis of tubby-like protein family in Arabidopsis，rice，and poplar [J]. Genomics，2008，92（4）： 246-253.

[21]向貴生，王開錦，晏慧君，等.薔薇科植物MLO蛋白家族的生物信息學(xué)分析[J].基因組學(xué)與應(yīng)用生物學(xué)，2018，37（5）：2043-2059.

（責(zé)任編輯：陳海霞）