轉(zhuǎn)錄組的黃山欒樹MYB基因家族鑒定及表達1)

呂運舟 董筱昀 梁珍海 孫海楠 黃利斌 蔣澤平

(江蘇省林業(yè)科學研究院，南京，211153)

MYB轉(zhuǎn)錄因子家族廣泛參與植物次生代謝、細胞分化、抗逆性等生理生化過程，其成員都具有高度保守的DNA結(jié)合結(jié)構(gòu)域，由52個氨基酸構(gòu)成的重復基序組成[1-3]。根據(jù)重復基序數(shù)目(R)可以分為1R-MYB、R2R3-MYB、3R-MYB、4R-MYB4種類型[3]。目前，MYB轉(zhuǎn)錄因子已在擬南芥(Arabidopsisthaliana)等模式物種中得到分離和功能鑒定[2-3]，在毛果楊(Populustrichocarpa)[4]、杉木(Cunninghamialanceolata)[5]、楠木(Phoebezhennan)[6]、桑樹(Morusnotabilis)[7]等喬木樹種基因組或轉(zhuǎn)錄組水平得到鑒定。

植物葉片顏色是由葉片中色素的組分、含量及分布決定[8-9]，包括葉綠素、類黃酮類和類胡蘿卜素類色素。類黃酮類色素和類胡蘿卜素的種類和含量決定了植物界的繽紛色彩，許多研究結(jié)果已經(jīng)揭示了它們的生物合成過程[10]。MYB轉(zhuǎn)錄因子對黃酮類化合物和類胡蘿卜素的生物合成起到重要的調(diào)控作用，如擬南芥MYB基因家族的MYB11[11]、MYB12[12-13]、MYB112[14]通過調(diào)控CHS、CHI、F3H和FLS1基因的表達實現(xiàn)對花青素等類黃酮類化合物生物合成的調(diào)控作用；彩艷龍頭花(Mimuluslewisii)R2R3-MYB類基因RCP1可以激活所有類蘿卜素合成基因而激活花瓣中類胡蘿卜素的生物合成[15]。此外，Amponmah et al.[16]發(fā)現(xiàn)獼猴桃(Actinidiadeliciosa)MYB7基因在轉(zhuǎn)錄水平正向調(diào)控類胡蘿卜素和葉綠素的生物合成。葉色變異植株是彩葉新品種和植物色素合成功能基因組學研究的重要材料[9，17]。‘金焰彩欒’(Koelreuteriabipinnata‘Jinyan’)是黃山欒樹(Koelreuteriabipinnatavar.integrifoliola)的天然黃化突變體，枝干呈橙黃色，春季葉片橙紅色，秋季全株葉片金黃，具有極高的觀賞價值[18-19]。本研究以‘金焰彩欒’及其原種黃山欒樹轉(zhuǎn)錄組數(shù)據(jù)入手，結(jié)合葉片發(fā)育不同時期高通量測序數(shù)據(jù)，挖掘鑒定黃山欒樹MYB轉(zhuǎn)錄因子家族基因，并對其預測蛋白特征及表達特征進行比較分析，旨在鑒別與類黃酮類或類胡蘿卜素調(diào)控相關(guān)MYB轉(zhuǎn)錄因子，探討‘金焰彩欒’葉片呈色的分子機制。

1 材料與方法

1.1 供試材料

供試黃山欒樹及其黃葉變異株‘金焰彩欒’種植于江蘇省林業(yè)科學研究院試驗苗圃內(nèi)。

1.2 RNA提取、cDNA文庫建立及測序

分別從固定樣株中采集春季新葉、夏季成熟葉及秋季變色期葉片，分別為T1(‘金焰彩欒’秋季葉片)、T2(實生對照秋季葉片)、T3(‘金焰彩欒’春季葉片)、T4(實生對照春季葉片)、T5(‘金焰彩欒’夏季葉片)及T6(實生對照夏季葉片)，液氮速凍；總RNA按照Trizol試劑盒(Thermo Fisher,USA)提取，測序文庫構(gòu)建及測序由北京百邁客生物科技有限公司完成，使用Illumina HiSeq 2500高通量測序技術(shù)進行測序[17]。

1.3 MYB基因鑒定

利用pfam數(shù)據(jù)庫PFAM(http://pfam.xfam.org/)的MYB基因家族結(jié)構(gòu)域hmmer隱馬爾科夫模型文件(PF00249)及hmmsearch軟件對轉(zhuǎn)錄組數(shù)據(jù)中所有的蛋白序列進行搜索，篩選期望值<0.001為候選基因家族基因；手動確認候選基因中的結(jié)構(gòu)域序列,利用關(guān)鍵“MYB”進行檢索，結(jié)合NCBI CDD(https://www.ncbi.nlm.nih.gov/Structure/bwrpsb/bwrpsb.cgi)工具，刪除不含MYB結(jié)構(gòu)域或重復序列。

1.4 蛋白序列及系統(tǒng)進化分析

從擬南芥數(shù)據(jù)庫TAIR(https://www.arabidopsis.org/)下載模式植物擬南芥的MYB基因蛋白序列。提交MYB蛋白序列至ExPASy(http://web.expasy.org/protparam/)，ProtParam在線分析蛋白序列性質(zhì)。利用多序列比對軟件muscle 3.8.31對97個黃山欒樹和130個擬南芥MYB蛋白進行多序列比對分析，然后用軟件MEGA7.0構(gòu)建鄰接法(neighbor joining，NJ)系統(tǒng)進化樹，使用p-distance距離模型構(gòu)建，選擇成對刪除選項，bootstrap值設置為1 000。

1.5 表達特征分析及驗證

利用97個MYB候選基因在黃山欒樹轉(zhuǎn)錄數(shù)據(jù)，包括‘金焰彩欒’和普通黃山欒樹的3個時期(春季新葉、夏季成熟葉和秋季衰老葉)的基因表達數(shù)據(jù)，并對表達量數(shù)據(jù)FPKM[20]進行歸一化轉(zhuǎn)換，然后繪制熱圖。

此外，選取10個與類黃酮類及苯丙氨酸類化合物代謝調(diào)控相關(guān)的MYB候選基因，選用看家基因KbGAPDH為內(nèi)參[21]。定量PCR引物利用Premier 5在非保守域設計，引物信息參見表1；RNB100 RNA提取試劑盒(SIGMA)提取以上6個組織總RNA，逆轉(zhuǎn)錄產(chǎn)物cDNA稀釋10倍，使用SYBR Premix Ex TaqTM(TAKARA)定量試劑盒，在CFX-96熒光定量PCR儀(Bio-Rad)上進行定量PCR。反應體系為20 μL，其中SYBR mix 10 μL，上下游引物(10 μmol·L-1)各0.8 μL，cDNA模板1.0 μL，并加6.4 μL ddH2O，每個樣本3次重復。PCR擴增反應程序:95 ℃預變性30 s，95 ℃變性5 s，60 ℃退火20 s，共40個循環(huán)。

表1 用于定量分析的引物序列

2 結(jié)果與分析

2.1 測序結(jié)果

經(jīng)過測序質(zhì)量控制，共得到31.44 Gb有效數(shù)據(jù)，各樣品Q30堿基百分比均不小于90.74%。經(jīng)Trinity軟件組裝后，共得到175 881條轉(zhuǎn)錄本和94 125條單基因簇，轉(zhuǎn)錄本與單基因簇的N50分別為2 791和1 340，組裝完整性較高。生物學分析比對發(fā)現(xiàn)，共有47 341條單基因簇在與NR、Swiss-Prot、GO、COG、KOG、KEGG數(shù)據(jù)庫得到注釋，這些序列信息為后續(xù)研究提供數(shù)據(jù)基礎。

2.2 黃山欒樹MYB基因家族的鑒定及編碼蛋白

基于黃山欒樹及其突變體葉片轉(zhuǎn)錄組數(shù)據(jù)的本地Blast比對檢索，在NR、NT、Swiss-Prot等數(shù)據(jù)庫進行注釋，再經(jīng)NCBI CDD工具檢驗，去除過短和重復冗余的序列后，共獲得97個黃山欒樹MYB轉(zhuǎn)錄因子，分別命名為KbMYB1～KbMYB97。從表2可以看出，包括48個1R-MYB類轉(zhuǎn)錄因子，47個R2R3-MYB類轉(zhuǎn)錄因子和2個3R-MYB類轉(zhuǎn)錄因子得到鑒定，預測編碼蛋白長度介于52(KbMYB89)～1780KbMYB25)個氨基酸，理論分子質(zhì)量大小介于52.0～193.35 ku，理論等電點范圍在4.51(KbMYB16)～11.87(KbMYB51)。

2.3 黃山欒樹MYB轉(zhuǎn)錄因子進化分析

利用MEGA7.0構(gòu)建的黃山欒樹與擬南芥MYB轉(zhuǎn)錄因子家族系統(tǒng)進化樹，拓撲結(jié)構(gòu)如圖1所示，可以分為23亞組，亞組G23僅包含22個黃山欒樹轉(zhuǎn)錄因子，亞組G3、G6和G21僅包含擬南芥MYB轉(zhuǎn)錄因子。19個亞組由黃山欒樹和擬南芥MYB轉(zhuǎn)錄因子共同組成，其中亞家族G2、G9、G14中只有1個黃山欒樹MYB基因家族成員。此外，在同一亞組中擬南芥基因和黃山欒樹基因之間彼此相互分開，說明兩者遺傳關(guān)系較遠。已有研究表明親緣關(guān)系較近的MYB轉(zhuǎn)錄因子具有相似生物學功能[2-3]，基于擬南芥中的研究，G1、G2、G4、G7、G8、G9和G13分別與次生生長、生長發(fā)育及脅迫應答等調(diào)控有關(guān)。

表2 黃山欒樹MYB基因家族編碼蛋白信息

圖1 黃山欒樹和擬南芥MYB類蛋白家族的進化分析

2.4 MYB轉(zhuǎn)錄因子在‘金焰彩欒’葉片呈色過程中的表達模式

利用黃山欒樹及其黃葉突變體‘金焰彩欒’春、夏、秋3個時期的葉片轉(zhuǎn)錄組數(shù)據(jù)庫，提取97個MYB基因?qū)腞PKM值(代表每百萬片段中來自于某基因每千堿基長度的片段數(shù))，并經(jīng)歸一化處理制作Heml熱圖(如圖2所示)。結(jié)果表明，MYB基因在‘金焰彩欒’及普通黃山欒樹不同時期葉片中表達差異明顯。根據(jù)表達模式特征可以將黃山欒樹MYB基因家族成員分為4類，I～III類在‘金焰彩欒’和普通欒樹中表達模式基本一致，其中I組為春季葉片中表達豐度較少，在夏季上調(diào)表達，秋季又下調(diào)表達；II組為春季和夏季葉片中表達量較少，在秋季上調(diào)表達；III組為春季葉片中表達豐度較高，在夏季和秋季葉片中表達下調(diào)。而IV組中的MYB基因在‘金焰彩欒’的3個時期葉片中表達量均較高，而在對照中呈現(xiàn)春季葉片表達豐度較高，其它兩個時期較低，這些MYB基因可能參與調(diào)控突變體葉片呈色過程。

2.5 類黃酮類等次盛化合物合成調(diào)控候選MYB基因表達

根據(jù)MYB轉(zhuǎn)錄因子表達特征及聚類分析結(jié)果，挑選與擬南芥類黃酮類、苯丙烷類等次生代謝合成調(diào)控相關(guān)MYB轉(zhuǎn)錄因子蛋白同源性較近的10個黃山欒樹MYB候選基因。實時定量分析結(jié)果表明(表3)，春季葉片相對高表達的MYB轉(zhuǎn)錄因子包括KbMYB40，KbMYB65，KbMYB80和KbMYB91，其中KbMYB80在‘金焰彩欒’表達量相對較高，而KbMYB40相反。候選基因KbMYB6和KbMYB8，預測與花青素等類黃酮類化合物合成正向調(diào)控相關(guān)，表達模式一致，在普通欒樹夏季成熟葉中表達較高，秋季表達下降。花青素合成負向調(diào)控候選基因KbMYB60在供試樣本中表達量較為一致，呈高水平狀態(tài)。而苯丙烷類合成調(diào)控候選基因KbMYB88在‘金焰彩欒’和普通欒樹的春、秋時期葉片中表達基本一致，在夏季葉片中表達量增加，且在普通欒樹中表達量最高。

T1、T3、T5與T2、T4、T6分別為‘金焰彩欒’與對照的秋季、春季及夏季葉片。

3 結(jié)論與討論

黃山欒樹是我國優(yōu)良的鄉(xiāng)土樹種，具有多種經(jīng)濟價值和園林觀賞價值，是園林綠化、庭院配植的理想樹種[18]。MYB轉(zhuǎn)錄因子在類黃酮類、類蘿卜素及苯丙烷類化合物代謝調(diào)控關(guān)系密切，在植物葉片、花瓣及果實呈色調(diào)控中發(fā)揮至關(guān)重要的作用[22-26]。通過對黃山欒樹葉片轉(zhuǎn)錄組數(shù)據(jù)挖掘，共鑒定出97個MYB基因，預測編碼蛋白長度介于52～1 780個氨基酸。通過與130個擬南芥MYB(AtMYB)蛋白聚類分析，預測到多個MYB基因可能參與類黃酮類、類胡蘿卜素類和苯丙烷類化合物代謝通路功能基因表達調(diào)控。Dubos et al.[3]根據(jù)結(jié)構(gòu)及功能將擬南芥MYB基因家族分為27個亞家族，其中AtMYB11、AtMYB12與AtMYB112[11-14]；AtMYB75[27]、AtMYB90[3]與AtMYB123[28]分別參與黃酮和花青素生物合成過程。共有13個黃山欒樹MYB基因與它們歸為4個亞組(G1、G2與G4)，涉及KbMYB5、KbMYB9、KbMYB80、KbMYB91，KbMYB6、KbMYB8、KbMYB44、KbMYB79、KbMYB93與KbMYB74、KbMYB40、KbMYB65，它們可能參與黃酮和花青素的生物合成調(diào)控。AtMYB58、AtMYB63[29]，AtMYB46[30]和AtMYB72[30]參與苯丙烷類化合物代謝而影響植物次生細胞壁合成，間接影響黃酮類代謝通路。由此推測，KbMYB16和KbMYB88參與黃山欒樹葉片細胞次生壁合成過程。

表3 候選MYB基因?qū)崟r熒光定量PCR分析

‘金焰彩欒’是黃山欒樹實生變異的黃化突變體，觀賞價值極高[18]。研究發(fā)現(xiàn)，金焰彩欒’整個生長期中的葉綠素a、葉綠素b、總?cè)~綠素含量和類胡蘿卜素、花色素苷含量均顯著低于原種，其光合能力減弱，生長較原種受到部分影響[19]。此外，研究發(fā)現(xiàn)脫落醛氧化酶基因的下調(diào)表達導致了類胡蘿卜素含量降低，葉綠素a酸酯氧化酶基因的下調(diào)表達與葉綠體合成受阻有關(guān)[17]，而在功能酶轉(zhuǎn)錄水平調(diào)控上也存在差異。表達模式分析發(fā)現(xiàn)，有25個MYB基因在春季葉片中表達較高，有21個在‘金焰彩欒’春、夏、秋葉片中表達量均高于普通欒樹，這些基因可能參與了‘金焰彩欒’葉片呈色調(diào)控。過表達AtMYB60可以抑制植物花青素生物合成[31-32]，KbMYB60與其進化關(guān)系最近，且在‘金焰彩欒’各時期葉片中表達均較高，預測其參與了花青素生物合成的負調(diào)控。進一步驗證此類MYB基因功能，有助于揭示‘金焰彩欒’葉片變異的調(diào)控機理，為后續(xù)彩葉植物新品種開發(fā)提供分子基礎。