普通煙草熱激轉(zhuǎn)錄因子家族的全基因組鑒定與分析

郭 存，王 奇，李曉旭，張震彪，文利超，鄧智超，初雨蒙，劉 濤，崔萌萌，郭永峰*

普通煙草熱激轉(zhuǎn)錄因子家族的全基因組鑒定與分析

郭 存1,2，王 奇1，李曉旭3，張震彪1,2，文利超1,2，鄧智超1,2，初雨蒙1，劉 濤1,2，崔萌萌1，郭永峰1*

（1.中國農(nóng)業(yè)科學(xué)院煙草研究所，煙草行業(yè)基因資源利用重點(diǎn)實(shí)驗(yàn)室，青島 266101；2.中國農(nóng)業(yè)科學(xué)院研究生院，北京 100081；3.湖南中煙工業(yè)有限責(zé)任公司技術(shù)中心，長沙 410007）

熱激轉(zhuǎn)錄因子（heat shock transcription factor, Hsf）家族成員廣泛參與植物生長發(fā)育和非生物脅迫的調(diào)控。本研究應(yīng)用生物信息學(xué)和比較基因組學(xué)方法從普通煙草中鑒定了49個(gè)Hsf家族成員，發(fā)現(xiàn)這些Hsf家族成員均具有較為保守的寡聚化結(jié)構(gòu)域。通過系統(tǒng)進(jìn)化分析可將煙草和擬南芥的Hsf家族成員劃分成3個(gè)分支與14個(gè)亞家族。共線性分析顯示煙草Hsf基因與雙子葉植物（擬南芥，番茄和葡萄）基因組成共線基因?qū)Φ臄?shù)量要大于與單子葉植物（水稻和玉米）基因組成共線基因?qū)Φ臄?shù)量。全基因組復(fù)制預(yù)測共有12個(gè)基因可組成7條復(fù)制基因?qū)Α?dòng)子分析發(fā)現(xiàn)基因啟動(dòng)子區(qū)域存在著多個(gè)與非生物脅迫相關(guān)的順式作用元件。表達(dá)模式分析證實(shí)家族成員在不同組織的表達(dá)具有差異，并且大多數(shù)家族成員可響應(yīng)熱和干旱脅迫。研究結(jié)果可為煙草Hsf轉(zhuǎn)錄因子家族成員的功能研究奠定生物信息學(xué)基礎(chǔ)。

普通煙草；熱激轉(zhuǎn)錄因子；共線性；全基因組復(fù)制；表達(dá)模式

高溫是影響植物生長的一種非生物脅迫，高溫環(huán)境下作物產(chǎn)量和品質(zhì)都會(huì)受到嚴(yán)重影響。熱激轉(zhuǎn)錄因子（Heat shock transcription factor, Hsf）是植物體內(nèi)一類能夠調(diào)節(jié)高溫反應(yīng)的蛋白[1-2]。有研究表明，在高溫環(huán)境下，Hsf家族成員通過與下游相關(guān)基因的熱激元件結(jié)合進(jìn)而激活其轉(zhuǎn)錄，從而提高植物對高溫的抗性[3-4]。Hsf家族成員含有較為保守的寡聚化結(jié)構(gòu)域，根據(jù)寡聚化結(jié)構(gòu)的特征可將Hsf家族劃分成3個(gè)亞支（A、B和C）[5-7]。

目前，熱激轉(zhuǎn)錄因子已在多種植物基因組中發(fā)現(xiàn)，在擬南芥[5]、水稻[8]、馬鈴薯[9]、番茄[10]、小麥[11]、玉米[12]和葡萄[13]中分別鑒定出了21、25、27、24、82、25和19個(gè)Hsf成員。研究發(fā)現(xiàn)，擬南芥大部分與高溫脅迫相關(guān)的基因均受到A分支Hsf家族成員的調(diào)控，其中包括、、和。在三重突變體中，高溫響應(yīng)基因的轉(zhuǎn)錄水平顯著降低[4,14]。此外，與野生型相比，四重突變體表現(xiàn)出生長遲緩且對溫度更敏感的表型。在番茄和擬南芥中分別過表達(dá)和基因可以增強(qiáng)植物耐旱性和耐熱性[15-16]。此外在煙草中異源過表達(dá)或基因都能夠提高煙草的耐熱性[17-18]。據(jù)報(bào)道，和除了在植物的耐熱性方面起作用外，還參與調(diào)控植物耐缺氧性和氧化應(yīng)激耐受性反應(yīng)[19-20]。在高鹽和干旱環(huán)境下，基因的表達(dá)量顯著上調(diào)，并且在ABA介導(dǎo)的鹽、干旱和高溫脅迫中起正向調(diào)節(jié)作用[21]。

煙草是一種重要的經(jīng)濟(jì)作物，高溫、干旱等非生物脅迫影響其正常生長發(fā)育，并誘導(dǎo)多種根莖類病害發(fā)生，進(jìn)而影響煙葉的產(chǎn)量和品質(zhì)。Hsf轉(zhuǎn)錄因子參與多種非生物脅迫的響應(yīng)，而煙草中Hsf家族成員信息還未見報(bào)道。本研究首次從煙草基因組序列中鑒定出49個(gè)Hsf轉(zhuǎn)錄因子家族成員，并對其進(jìn)行了系統(tǒng)進(jìn)化、保守基序、共線性、復(fù)制事件、順式作用元件和表達(dá)模式等分析，為進(jìn)一步解析基因在煙草非生物脅迫響應(yīng)過程中的功能奠定基礎(chǔ)。

1 材料與方法

1.1 材料與儀器試劑

本試驗(yàn)所用的材料為普通煙草品種K326，保存于煙草行業(yè)煙草基因資源利用重點(diǎn)實(shí)驗(yàn)室。K326種子經(jīng)過滅菌后均勻地撒在MS培養(yǎng)基中，在光照培養(yǎng)箱中培養(yǎng)8 周（光周期：16 h光照/8 h黑暗；25 ℃）。挑選長勢一致的煙苗，一部分轉(zhuǎn)移到37 ℃培養(yǎng)箱中進(jìn)行高溫處理，另一部分吸干根部水分后轉(zhuǎn)移到濾紙上做干旱處理。在處理0、3、6 h時(shí)對煙苗整株取樣，并于液氮中保存。所有樣品取3次生物學(xué)重復(fù)。

本研究用的RNA提取試劑盒購自北京康為世紀(jì)公司，反轉(zhuǎn)錄試劑盒與SYBR?Green Premix Pro Taq HS qPCR Kit熒光定量試劑盒購自湖南艾科瑞生物工程有限公司。熒光定量儀器Roche LightCycler 480購自美國羅氏公司。

1.2 方法和條件

1.2.1 NtHsf轉(zhuǎn)錄因子家族成員的鑒定和序列分析 從擬南芥基因組數(shù)據(jù)庫（https://www.arabidopsis.org/）中下載已報(bào)道的AtHsf家族成員的蛋白序列[22]，并利用MAFFT工具在默認(rèn)參數(shù)下進(jìn)行序列比對，利用HMMER工具基于比對結(jié)果建立HMM文件。普通煙草K326蛋白序列數(shù)據(jù)庫從茄科基因組數(shù)據(jù)庫（https://solgenomics.net/organism，Sol Genomics Network）中下載得到，并用建立的HMM文件在默認(rèn)參數(shù)下對數(shù)據(jù)庫進(jìn)行檢索，初步得到煙草NtHsf序列。Pfam數(shù)據(jù)庫（https://pfam.xfam.org/）被用來分析所鑒定的蛋白序列，并去除不含有寡聚化結(jié)構(gòu)域的序列。利用ProtParam工具（https://web.expasy.org/protparam/）對新鑒定的NtHsf序列進(jìn)行理化性質(zhì)分析。

1.2.2 NtHsf轉(zhuǎn)錄因子家族成員的系統(tǒng)進(jìn)化分析 利用MAFFT在默認(rèn)參數(shù)下對已報(bào)道的擬南芥AtHsf成員和新鑒定的煙草NtHsf成員進(jìn)行蛋白全長序列的比對[23]。基于比對結(jié)果，利用MGEA X進(jìn)行鄰接樹的構(gòu)建。構(gòu)建鄰接樹的參數(shù)設(shè)置如下：采用Poisson Model和Pairwise Deletion，Bootstrap檢驗(yàn)設(shè)為1000次[24]。

1.2.3的共線性及復(fù)制事件分析 按照先前的研究，分別下載普通煙草和其他5個(gè)物種（擬南芥、番茄、普通、水稻和玉米）的基因組注釋文件，利用McScanX在默認(rèn)參數(shù)下進(jìn)行共線性分析，使用TBtools將結(jié)果進(jìn)行可視化展示[25]。使用ColinearScan軟件對煙草NtHsf基因的復(fù)制事件進(jìn)行預(yù)測和可視化[26]。

1.2.4 NtHsf的寡聚化結(jié)構(gòu)域序列分析 手動(dòng)提取NtHsf的寡聚化結(jié)構(gòu)域的序列，并利用MAFFT在默認(rèn)參數(shù)對其進(jìn)行比對與可視化[23]。

1.2.5 啟動(dòng)子分析 選取基因上游2000 bp區(qū)域[27]，并提交到PlantCare網(wǎng)站（http://bioinformatics.psb.ugent.be/webtools/plantcare/html/）進(jìn)行基因啟動(dòng)子順式作用元件的鑒定[28]。

1.2.6 轉(zhuǎn)錄組數(shù)據(jù)分析及可視化 從NCBI的GEO數(shù)據(jù)庫中下載得到普通煙草K326的根、莖和莖尖組織的轉(zhuǎn)錄組數(shù)據(jù)（https://www.ncbi.nlm.nih.gov，登錄號(hào)為GSE95717）[29]。從下載的數(shù)據(jù)中提取煙草基因的表達(dá)情況，利用TBtools軟件中的HeatMap程序?qū)虻谋磉_(dá)量進(jìn)行標(biāo)準(zhǔn)化分析和可視化展示[25]。

1.2.7 煙草樣品RNA的提取和反轉(zhuǎn)錄 利用康為世紀(jì)Ultrapure RNA Kit進(jìn)行樣品總RNA的提取，cDNA第一條鏈的合成在說明書參照下進(jìn)行，并將反轉(zhuǎn)錄的cDNA保存在?20 ℃待用。

1.2.8 熒光定量PCR 熒光定量PCR在Roche LightCycler 480上進(jìn)行，試驗(yàn)體系設(shè)置如下：2× SYBR?Green Pro Taq HS Premix 10mL, cDNA 1mL, Primer-F 0.4mL，Primer-R 0.4mL，無菌水8.2mL。PCR程序設(shè)置如下：95 ℃ 30s；95 ℃ 5s，60℃ 30s，40個(gè)循環(huán)。煙草基因（GenBank No. L18908）為內(nèi)參基因[30]，每個(gè)樣品進(jìn)行3次技術(shù)重復(fù)。試驗(yàn)數(shù)據(jù)采用2-△△Ct法進(jìn)行分析[31]。本試驗(yàn)所用到的引物如表1所示。

2 結(jié) 果

2.1 煙草NtHsf家族成員的鑒定和理化分析

本研究在煙草基因組中鑒定了49個(gè)家族成員（表2），并按照基因在染色體上的物理位置將新鑒定的煙草基因進(jìn)行命名（～）。新鑒定的49個(gè)基因中有24個(gè)成員錨定在15條煙草染色體上。利用ProtParam工具對NtHsf成員進(jìn)行理化性質(zhì)分析發(fā)現(xiàn)，49個(gè)NtHsf成員在理化性質(zhì)上存在較大差異，例如其氨基酸數(shù)目在146～504之間，相對分子量大小在16?521.7～55?842.1 Da之間，等電點(diǎn)范圍為4.45～10.12。

2.2 煙草NtHsf家族成員的進(jìn)化分析

為了進(jìn)一步研究煙草中NtHsf轉(zhuǎn)錄因子成員的進(jìn)化歷程，將新鑒定的49個(gè)煙草NtHsf成員和已報(bào)道的21個(gè)擬南芥AtHsf成員進(jìn)行系統(tǒng)進(jìn)化分析。如圖1所示，這些Hsf家族成員可劃分為A、B和C三個(gè)分支，其中A和B分支又可劃分為A1-A9和B1-B3亞家族。值得注意的是，每個(gè)亞家族中均含有來自擬南芥和煙草的Hsf家族成員，說明擬南芥和煙草的物種分化時(shí)間要晚于Hsf轉(zhuǎn)錄因子的分化時(shí)間。并且在各個(gè)亞家族中，煙草NtHsf家族成員要顯著多于擬南芥Hsf家族成員，例如A1、A4、A6、B1、B2、B3等亞家族，暗示在這兩個(gè)物種分化之后，各個(gè)亞家族的NtHsf成員可能又經(jīng)歷了復(fù)制事件，導(dǎo)致了基因的擴(kuò)張。

表1 熒光定量PCR引物信息

表2 普通煙草NtHsf轉(zhuǎn)錄因子家族成員鑒定及理化性質(zhì)分析

圖1 煙草和擬南芥Hsf家族成員的系統(tǒng)進(jìn)化分析

2.3 煙草NtHsf家族成員的寡聚化結(jié)構(gòu)域分析

Hsf家族成員的寡聚化結(jié)構(gòu)域由螺旋卷曲結(jié)構(gòu)和α螺旋結(jié)構(gòu)組成，其包含了大量的疏水氨基酸殘基[1]。為了驗(yàn)證煙草NtHsf家族成員的保守結(jié)構(gòu)域，選取NtHsf家族成員的寡聚化結(jié)構(gòu)域進(jìn)行多序列比對并進(jìn)行可視化。如圖2所示，新鑒定的NtHsf家族成員均具有較為保守的寡聚化結(jié)構(gòu)域，由HR-A和HR-B組成，其中A類Hsf成員有30個(gè)，C類2個(gè)，B類17個(gè)。

2.4 煙草NtHsf基因的共線性分析

為進(jìn)一步了解煙草家族成員的系統(tǒng)進(jìn)化關(guān)系，分析了煙草與其他5個(gè)物種（擬南芥、番茄、葡萄、水稻和玉米）的共線性關(guān)系。如圖3所示，灰色線條代表所有具有共線關(guān)系的基因，基因則用紅色線條標(biāo)記。其中，煙草中共22個(gè)基因與番茄基因成共線性關(guān)系，其次是葡萄（16個(gè)基因）、擬南芥（9個(gè)基因）、水稻（4個(gè)基因）和玉米（3個(gè)基因）。此外煙草基因與番茄、葡萄、擬南芥、水稻和玉米組成共線基因?qū)Φ臄?shù)量分別為30、19、11、11和7。值得注意的是，煙草中共有2個(gè)基因（和）與這5個(gè)物種的基因組成共線對，表明這2個(gè)基因可能在這些物種分化之前就已經(jīng)存在。

圖2 煙草NtHsf家族成員的寡聚化結(jié)構(gòu)域序列比對圖

2.5 煙草NtHsf基因家族的全基因組復(fù)制事件

系統(tǒng)進(jìn)化分析暗示煙草基因家族可能發(fā)生了全基因組復(fù)制事件，隨后利用McScan X對煙草家族成員進(jìn)行復(fù)制分析。如圖4所示，灰色線條代表煙草基因組中所有的復(fù)制基因?qū)Γt色線條表示基因。結(jié)果顯示，共有12個(gè)基因組成7條復(fù)制基因?qū)Γ凳具@些基因可能由全基因組復(fù)制事件產(chǎn)生，也說明全基因組復(fù)制事件在煙草基因家族擴(kuò)張過程中有一定的貢獻(xiàn)。

2.6 煙草NtHsf家族成員的啟動(dòng)子分析

利用PlantCARE網(wǎng)站對的啟動(dòng)子元件進(jìn)行探究，并選取12個(gè)順式作用元件進(jìn)行分析。如圖5所示，大多數(shù)基因啟動(dòng)子中存在激素響應(yīng)元件，例如響應(yīng)脫落酸的元件ABRE、響應(yīng)乙烯的ERE元件、響應(yīng)水楊酸的元件TCA-element和響應(yīng)茉莉酸甲酯的CGTCA-motif元件。值得注意的是，響應(yīng)生長素的AuxRR-core元件僅在個(gè)別基因的啟動(dòng)子區(qū)域中檢測到，例如、、、、、和。此外，厭氧誘導(dǎo)元件(ARE)、損傷響應(yīng)元件(WUN-motif)、脅迫響應(yīng)元件(TC-rich repeats)、WRKY轉(zhuǎn)錄因子結(jié)合位點(diǎn)(W-box)、MYB轉(zhuǎn)錄因子結(jié)合位點(diǎn)(MBS)和低溫響應(yīng)元件(LTR)等6種脅迫響應(yīng)元件也被檢測到廣泛地存在啟動(dòng)子區(qū)域中，暗示這些基因可能參與到不同脅迫的響應(yīng)過程中。同時(shí)，在一些基因的啟動(dòng)子區(qū)域檢測出與分生組織發(fā)育相關(guān)的順式作用元件(CAT-box)。

圖3 煙草和其他5個(gè)物種Hsf基因的共線性分析

圖4 煙草NtHsf基因家族的全基因組復(fù)制分析

圖5 煙草NtHsf基因啟動(dòng)子區(qū)順式作用元件分析

2.7 煙草NtHsf家族成員在不同組織中的表達(dá)分析

利用已發(fā)表的普通煙草轉(zhuǎn)錄組數(shù)據(jù)庫（https://www.ncbi.nlm.nih.gov，登錄號(hào)為GSE95717）對基因的表達(dá)模式進(jìn)行分析。如圖6所示，家族成員在煙草根、莖和莖尖組織中的表達(dá)差異較大。多數(shù)基因在3個(gè)煙草組織中均檢測到表達(dá)，如、和等，然而一些基因并未在這3個(gè)組織中檢測到表達(dá)量，如、、和。進(jìn)一步研究發(fā)現(xiàn)個(gè)別基因的表達(dá)模式具有組織特異性，例如、和在煙草的根和莖組織中特異高表達(dá)，在根和莖中檢測出表達(dá)量，而在莖尖中卻不表達(dá)。

2.8 煙草NtHsf家族成員在熱和干旱脅迫下的表達(dá)模式

利用熒光定量PCR的方法對煙草NtHsf家族成員在熱和干旱脅迫下的表達(dá)模式進(jìn)行分析。如圖7所示，一些基因可以被高溫誘導(dǎo)表達(dá)，尤其是、、、和基因的表達(dá)量在高溫處理3 h后顯著增加(>10倍)，而這些基因的表達(dá)水平在處理6 h后均下調(diào)，說明這些基因在高溫脅迫下的表達(dá)具有時(shí)間特異性。其余基因在高溫脅迫下的表達(dá)量均呈現(xiàn)出不同程度的上調(diào)表達(dá)，例如、和等基因。

圖6 煙草NtHsf家族成員的組織表達(dá)模式分析

圖8所示的為基因在干旱處理下的表達(dá)模式，結(jié)果發(fā)現(xiàn)和在干旱處理3 h后的表達(dá)量達(dá)到最高峰，分別為對照的6.5倍和8.4倍。且這兩個(gè)基因在干旱處理下的表達(dá)模式具有時(shí)間特異性，其表達(dá)量在處理3 h后達(dá)到頂峰，然后在處理6 h后下降，、、和等基因的表達(dá)模式也呈現(xiàn)出類似的趨勢。

圖7 普通煙草NtHsf家族成員在高溫處理下的表達(dá)模式分析

圖8 普通煙草NtHsf家族成員在干旱處理下的表達(dá)模式分析

3 討 論

煙草作為我國重要的經(jīng)濟(jì)作物，時(shí)常會(huì)受到高溫等環(huán)境因素的脅迫，對煙葉的品質(zhì)和生產(chǎn)造成負(fù)面影響。熱激轉(zhuǎn)錄因子被報(bào)道與植物的耐熱性和抗旱性有緊密的聯(lián)系，本研究采用比較基因組學(xué)的方法對煙草基因組中的NtHsf熱激轉(zhuǎn)錄因子進(jìn)行了鑒定。并從系統(tǒng)進(jìn)化、共線性、復(fù)制事件、順式作用元件和表達(dá)模式等方面對新鑒定的NtHsf轉(zhuǎn)錄因子進(jìn)行分析，為之后解析煙草NtHsf家族成員的生物學(xué)功能和調(diào)控機(jī)制提供基礎(chǔ)。

本研究共鑒定得到49個(gè)煙草NtHsf轉(zhuǎn)錄因子，進(jìn)一步分析發(fā)現(xiàn)其均含有較為保守的寡聚化結(jié)構(gòu)域。系統(tǒng)進(jìn)化分析表明煙草和擬南芥的Hsf家族成員可劃分為A、B和C三個(gè)分支和14個(gè)亞家族。有意思的是，每個(gè)亞家族中的NtHsf成員均比AtHsf成員多，暗示煙草中的NtHsf成員可能進(jìn)行一次或多次復(fù)制事件。隨后的全基因組復(fù)制事件顯示共有12個(gè)基因經(jīng)歷了復(fù)制過程，暗示這12個(gè)基因可能具有類似功能，也從側(cè)面證明復(fù)制事件在煙草NtHsf成員的擴(kuò)張過程中起著重要作用。

系統(tǒng)進(jìn)化和共線性的綜合分析是探索同源基因保守功能的有效手段之一，本研究通過分析煙草與擬南芥Hsf成員之間的系統(tǒng)進(jìn)化和共線性關(guān)系，初步推測和探討NtHsf成員的功能。在亞家族A1中，和與聚在一起，并且分別與這兩個(gè)煙草Hsf成員組成共線基因?qū)Γf明和與基因直系同源。進(jìn)一步分析發(fā)現(xiàn)和組成一對復(fù)制基因?qū)Γ凳酒涔δ芸赡芟嗨?。此外，在番茄中過表達(dá)可增強(qiáng)番茄幼苗的耐熱性[15]，而且表達(dá)模式發(fā)現(xiàn)在高溫脅迫下，和基因表達(dá)上調(diào)，暗示這兩個(gè)煙草基因可能參與到煙草的高溫脅迫響應(yīng)中。值得注意的是，基因在干旱脅迫下的表達(dá)量也顯著上調(diào)，暗示該基因可能參與到煙草多脅迫響應(yīng)調(diào)控中。

此外，被報(bào)道通過調(diào)控?zé)峒し磻?yīng)相關(guān)基因的表達(dá)來提高擬南芥耐熱性[16]。在本研究中，與聚在同一亞家族中，并且其編碼基因預(yù)測具有共線性關(guān)系，說明其直系同源。有趣的是，在的啟動(dòng)子上發(fā)現(xiàn)了多個(gè)與脅迫相關(guān)的順式作用元件，例如ABRE、G-box、W-box和MYB，且其轉(zhuǎn)錄本在高溫脅迫下明顯增加，暗示可能參與到煙草的高溫脅迫過程。此外，也與聚在同一亞家族中，且其編碼基因來自同一復(fù)制事件。在高溫脅迫下的表達(dá)量也有所增加，暗示可能與在煙草響應(yīng)高溫脅迫方面具有相似的功能。

4 結(jié) 論

本研究利用生物信息學(xué)等方法在普通煙草基因組中鑒定到49個(gè)熱激轉(zhuǎn)錄因子基因，它們編碼的蛋白均含有較為保守的寡聚化結(jié)構(gòu)域。此外，煙草NtHsf成員和擬南芥AtHsf成員被劃分成3個(gè)亞支和14個(gè)亞家族。共線性分析表明，基因和基因形成11對同源基因?qū)?，暗示這些共線基因可能具有保守的功能。啟動(dòng)子分析發(fā)現(xiàn)大量脅迫響應(yīng)相關(guān)的順式作用元件存在NtHsf成員的啟動(dòng)子上，并且多數(shù)基因能夠響應(yīng)高溫和干旱脅迫，說明家族成員可能在煙草的高溫和干旱脅迫過程中具有重要的功能。

[1] SCHARF K D, BERBERICH T, EBERSBERGER I, et al. The plant heat stress transcription factor (Hsf) family: structure, function and evolution[J]. Biochimica et Biophysica Acta (BBA)-Gene Regulatory Mechanisms, 2012, 1819(2): 104-119.

[2] VON KOSKULL-D?RING P, SCHARF K D, NOVER L. The diversity of plant heat stress transcription factors[J]. Trends in Plant Science, 2007, 12(10): 452-457.

[3] IKEDA M, MITSUDA N, OHME-TAKAGI M.HsfB1 and HsfB2b act as repressors of the expression of heat-inducible Hsfs but positively regulate the acquired thermotolerance[J]. Plant Physiology, 2011, 157(3): 1243-1254.

[4] YOSHIDA T, OHAMA N, NAKAJIMA J, et al.HsfA1 transcription factors function as the main positive regulators in heat shock-responsive gene expression[J]. Molecular Genetics and Genomics, 2011, 286(5): 321-332.

[5] GUO M, LIU JH, MA X, et al. The plant heat stress transcription factors (HSFs): structure, regulation, and function in response to abiotic stresses[J]. Frontiers in Plant Science, 2016, 7: 114-121.

[6] KOTAK S, PORT M, GANGULI A, et al. Characterization of C-terminal domains ofheat stress transcription factors (Hsfs) and identification of a new signature combination of plant class A Hsfs with AHA and NES motifs essential for activator function and intracellular localization[J]. The Plant Journal, 2004, 39(1): 98-112.

[7] YE J, YANG X, HU G, et al. Genome-wide investigation of heat shock transcription factor family in wheat (L.) and possible roles in anther development[J]. International Journal of Molecular Sciences, 2020, 21(2): 608-619.

[8] WANG C, ZHANG Q, SHOU H X. Identification and expression analysis ofin rice[J]. Journal of Zhejiang University Science B, 2009, 10(4): 291-300.

[9] TANG R, ZHU W, SONG X, et al. Genome-wide identification and function analyses of heat shock transcription factors in potato[J]. Frontiers in Plant Science, 2016, 7: 490-508.

[10] YANG X, ZHU W, ZHANG H, et al. Heat shock factors in tomatoes: genome-wide identification, phylogenetic analysis and expression profiling under development and heat stress[J]. PeerJ, 2016, 4: e1961.

[11] DUAN S, LIU B, ZHANG Y, et al. Genome-wide identification and abiotic stress-responsive pattern of heat shock transcription factor family inL.[J]. BMC Genomics, 2019, 20(1): 1-20.

[12] LIN Y X, JIANG H Y, CHU Z X, et al. Genome-wide identification, classification and analysis of heat shock transcription factor family in maize[J]. BMC Genomics, 2011, 12(1): 1-14.

[13] GUOTIAN L, FENGMEI C, YI W, et al. Genome-wide identification and classification of HSF family in grape, and their transcriptional analysis under heat acclimation and heat stress[J]. Horticultural Plant Journal, 2018, 4(4): 133-143.

[14] LIU HC, LIAO HT, CHARNG YY. The role of class A1 heat shock factors (HSFA1s) in response to heat and other stresses in[J]. Plant, Cell & Environment, 2011, 34(5): 738-751.

[15] LI H Y, CHANG C S, LU L S, et al. Over-expression ofheat shock factor gene () enhances chilling tolerance in transgenic tomato[J]. Botanical Bulletin of Academia Sinica, 2003, 44: 129-140.

[16] HIGASHI Y, OHAMA N, ISHIKAWA T, et al. HsfA1d, a protein identified via FOX hunting usingcDNAs improves heat tolerance by regulating heat-stress-responsive gene expression[J]. Molecular Plant, 2013, 6(2): 411-422.

[17] ZHU X, WANG Y, LIU Y, et al. Overexpression oftranscription factor fromimproved heat tolerance of transgenic tobacco[J]. PLoS One, 2018; 13(11): e0207277.

[18] XU L, GAO J, GUO L, et al. Heat shock protein 70 (HmHsp70) fromconfers thermotolerance to tobacco[J]. AMB Express, 2020; 10(1): 1-10.

[19] LI C, CHEN Q, GAO X, et al.modulates expression of stress responsive genes and enhances tolerance to heat and oxidative stress in[J]. Science in China Series C: Life Sciences, 2005, 48(6): 540-550.

[20] SCHRAMM F, LARKINDALE J, KIEHLMANN E, et al. A cascade of transcription factor DREB2A and heat stress transcription factor HsfA3 regulates the heat stress response of[J]. The Plant Journal, 2008, 53(2): 264-274.

[21] HUANG Y C, NIU C Y, YANG C R, et al. The heat stress factor HSFA6b connects ABA signaling and ABA-mediated heat responses[J]. Plant Physiology, 2016, 172(2): 1182-1199.

[22] 李曉旭，郭存，蒲文宣，等. 普通煙草WOX轉(zhuǎn)錄因子家族的全基因組鑒定及分析[J]. 中國煙草學(xué)報(bào)，2021，27（1）：90-100.

LI X X, GUO C, PU W X, et al. Genome-wide identification and systemic analysis of WOX family genes in tobacco[J]. Acta Tabacaria Sinica, 2021, 27 (1): 90-100.

[23] KATOH K, STANDLEY D M. MAFFT multiple sequence alignment software version 7: improvements in performance and usability[J]. Molecular Biology and Evolution, 2013, 30(4): 772-780.

[24] KUMAR S, STECHER G, LI M, et al. MEGA X: molecular evolutionary genetics analysis across computing platforms[J]. Molecular Biology and Evolution, 2018, 35(6): 1547-1549.

[25] CHEN C, 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.

[26] 李志遠(yuǎn)，宋青松，劉江，等. 普通煙草脂氧合酶基因家族鑒定及表達(dá)模式分析[J]. 中國煙草科學(xué)，2021，42（5）：86-94.

LI Z Y, SONG Q S, LIU J, et al. Identification and expression analysis of the lipoxygenase gene family in tobacco (L.) [J]. Chinese Tobacco Science, 2021, 42(5): 86-94.

[27] GUO C, WANG Q, LI Z Y, et al. Bioinformatics and expression analysis ofgenes reveal their potential functions in flower abscission and stress response in tobacco (L.)[J]. Frontiers in Genetics, 2021, 12: 670794.

[28] LESCOT M, DéHAIS P, THIJS G, et al. PlantCARE, a database of plant-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences[J]. Nucleic Acids Research, 2002, 30(1): 325-327.

[29] EDWARDS K, FERNANDEZ-POZO N, DRAKE-STOWE K, et al. A reference genome forenables map-based cloning of homeologous loci implicated in nitrogen utilization efficiency[J]. BMC Genomics, 2017, 18(1): 1-14.

[30] LI Z, CHAO J, LI X, et al. Systematic analysis of the bZIP family in tobacco and functional characterization ofinvolvement in salt stress[J]. Agronomy, 2021, 11(1): 148-165.

[31] LIVAK K J, SCHMITTGEN T D. Analysis of relative gene expression data using real-time quantitative PCR and the 2?ΔΔCTmethod[J]. Methods, 2001, 25(4): 402-408.

Genome-wide Identification and Systemic Analysis of the Hsf Gene Family inL.

GUO Cun1,2, WANG Qi1, LI Xiaoxu3, ZHANG Zhenbiao1,2, WEN Lichao1,2, DENG Zhichao1,2, CHU Yumeng1, LIU Tao1,2, CUI Mengmeng1, GUO Yongfeng1*

(1. Key Laboratory of Tobacco Gene Resources, Tobacco Research Institute of CAAS, Qingdao 266101, China; 2. Graduate School of Chinese Academy of Agricultural Science, Beijing 100081, China; 3. Technology Center, China Tobacco Hunan Industrial Co., Ltd., Changsha 410007, China)

Heat shock transcription factor (Hsf) family members are widely involved in the regulation of plant development and abiotic stress responses. In the current study, 49 Hsf members were identified from the tobacco genome by bioinformatics and comparative genomics methods, all of which contained conserved oligomerization domains. The Hsf members from tobacco andwere categorized into 3 clades and 14 subfamilies according to the phylogenetic analysis. The number of collinear gene pairs betweenandgenes from dicots (, tomato and grape) was more than the collinear gene pairs with thegenes from monocots (rice and maize). Besides, a total of 12genes were predicted to have arisen from duplication events. Promoter analysis showed that there were a variety of-elements related to abiotic stress responses on thepromoters. The various expression patterns ofin different tissues were detected, and mostgenes were induced by heat and drought stress. The results presented by this study might provide insights for the further studies on NtHsf members.

tobacco; heat shock transcription factor; collinearity; duplication events; expression pattern

S572.01

A

1007-5119（2022）03-0047-10

10.13496/j.issn.1007-5119.2022.03.008

中國農(nóng)業(yè)科學(xué)院科技創(chuàng)新工程（ASTIP-TRIC02）；湖南中煙工業(yè)有限責(zé)任公司科技項(xiàng)目（KY2022YC0010）；中央級(jí)公益性科研院所基本科研業(yè)務(wù)費(fèi)專項(xiàng)（1610232022009）

郭 存（1994-），男，在讀博士研究生，主要從事煙草育種研究。E-mail：82101191076@caas.cn。

，E-mail：guoyongfeng@caas.cn

2021-12-07

2022-04-15