滇水金鳳花斑形成相關基因IuMYB114 和IuMYB36 的克隆及表達分析

陶宇蝶孟丹晨黃美娟魏春梅趙秋燕周敏瞿素萍黃海泉

(1.西南林業(yè)大學園林園藝學院/國家林業(yè)和草原局西南風景園林工程技術研究中心/云南省功能性花卉資源及產(chǎn)業(yè)化技術工程研究中心/西南林業(yè)大學園林園藝花卉研發(fā)中心，云南 昆明 650224；2.云南省農(nóng)業(yè)科學院花卉研究所，云南 昆明 650205)

花色是被子植物花器官的重要表型性狀之一，在提高植物抗逆性、吸引傳粉者、維持繁衍進化基石等方面發(fā)揮著重要作用[1，2]。 花斑一般指在花瓣或花萼上位置、大小和形態(tài)基本固定的斑塊，是被子植物一個獨特的表型特征，是人工選擇創(chuàng)新、植物進化和自然選擇的綜合結果，能極大地豐富植物的觀賞價值，是現(xiàn)代花卉業(yè)的重要研究領域[3]。 因此，揭示花斑形成的分子機理對花色改良具有重要意義[4，5]。

植物花斑的形成是由于花色素基因在花瓣或其他成花部位上差異表達，且調控基因可以調節(jié)花色素合成催化酶基因在不同區(qū)域的表達水平，進而導致花斑的形成[6]。MYB基因作為植物中最龐大的一個基因家族，在花斑及色素形成、生長發(fā)育等過程中發(fā)揮著重要作用。 已在多種植物中發(fā)現(xiàn)調節(jié)花色素苷合成及與色斑形成相關的MYB 轉錄因子，且以R2R3 類型為主，也有單一R重復結構的MYB 轉錄因子[7]。 如AtMYB75、At-MYB90、AtMYB113參與花青素的生物合成[8]；在蝴蝶蘭中，PeMYB2調控唇瓣的全紅著色，PeMYBI1調控唇瓣的紅色斑點，PeMYBI2調控唇瓣紋理的形成[9]；R2R3－MYB 轉錄因子LhMYB6、Lh-MYB12和LhMYB12－Lat分別在百合的凸起式斑點[10]、飛濺式斑點[11]和非斑區(qū)表達[12]；cjMYB1基因可能在紅色山茶品種的花色苷合成途徑中起關鍵作用[13]；FtMYB23提高了苦蕎原花青素的合成與積累[14]；在月季紅花突變體中，RhMYB基因正調控花青素的合成[15]。 劉盛雨[16]研究發(fā)現(xiàn)，MYB114轉錄因子可能在血橙果實成熟過程中介導藍光對多甲氧基黃酮和花青素的合成調控；Yao 等[17]提出PyMYB114與已報道的Py-MYB10有著色增強的效應，并在草莓果實和煙草葉片中得到驗證；趙英[18]提出SmMYB36的同源基因首次出現(xiàn)在洛磯山耬斗菜中， 過表達SmMYB36的轉基因材料可以提高丹參酮類物質而抑制酚酸類物質。 但關于MYB36基因調控花斑形成的相關研究尚未見報道。

滇水金鳳(Impatiens uliginosa)為中國特有的鳳仙花科鳳仙花屬一年或多年生草本植物，花色艷麗，以紅色為主，花形獨特，旗瓣圓形，花期長，具有較高的園林應用價值[19]。 目前關于滇水金鳳的研究主要集中于種子萌發(fā)特性、花發(fā)育調控基因克隆、花色調控基因及金屬元素測定等方面[20－23]，尚未見關于其花斑的研究報道。 本研究從斑色系滇水金鳳花瓣中克隆得到IuMYB114和IuMYB36基因，并分析了其序列結構、系統(tǒng)發(fā)育和表達模式，為研究滇水金鳳花斑形成的分子調控機理以及鳳仙花花色改良等提供一定的基礎數(shù)據(jù)和理論依據(jù)。

1 材料與方法

1.1 試驗材料

供試斑色系滇水金鳳于2020 年10 月取自昆明市尋甸縣清水海，于滇水金鳳盛花期采集樣品(圖1)，迅速放于液氮中，并帶回實驗室保存于－80℃冰箱中備用。

圖1 斑色系滇水金鳳花器官

1.2 總RNA 的提取與IuMYB114 和IuMYB36 基因的克隆

以滇水金鳳整個花器官為試驗材料，采用RNA 提取試劑盒(OMEGA)提取總RNA；RNA 逆轉錄采用逆轉錄試劑盒(北京全式金生物技術股份有限公司)，－20℃保存?zhèn)溆谩?根據(jù)轉錄組數(shù)據(jù)設計擴增全長的特異性引物(表1)。 以RNA 逆轉錄合成的cDNA 第一鏈作為模板進行Iu-MYB114和IuMYB36基因全長cDNA 擴增。 回收產(chǎn)物并進行克隆，送上海生工生物工程有限公司測序。

表1 滇水金鳳MYB 基因cDNA 序列擴增引物

1.3 滇水金鳳花斑IuMYB114 和IuMYB36 序列分析

采用ExPasy－ProtParam 在線軟件進行基因編碼蛋白質的基本理化性質分析；運用TMHMM 2.0 對蛋白進行跨膜結構域的預測分析；利用NCBI 的CDD 工具預測目的基因結構域；運用SMART 對蛋白結構域進行分析，蛋白質的二級結構預測采用SOPMA 在線軟件進行，蛋白質的三級結構預測采用SWISS－MODE 在線軟件進行；運用NCBI 中的Protein BLAST 工具尋找同源序列，借助NCBI BLAST、DNAMAN 和MEGA 11 進行序列分析和系統(tǒng)進化樹構建。

1.4 IuMYB114 和IuMYB36 基因的表達模式

按照前述方法分別提取滇水金鳳旗瓣斑區(qū)和非斑區(qū)2 個部位(圖1)的RNA，反轉錄合成cDNA，設計qRT－PCR 引物(表2)，以滇水金鳳Actin1作為內參基因，進行qRT－PCR 試驗。 每個樣本設置3 個生物學重復。 將非斑區(qū)的表達水平設為對照，與斑區(qū)的進行比較。

表2 滇水金鳳MYB 基因qRT－PCR 擴增相關引物序列

2 結果與分析

2.1 滇水金鳳IuMYB114 和IuMYB36 基因的克隆

通過PCR 擴增，獲得IuMYB114與IuMYB36基因的cDNA 全長，分別為315 bp 和876 bp(圖2)，分別編碼104 個和151 個氨基酸。 對這兩個MYB基因的內含子和外顯子進行分析發(fā)現(xiàn)兩者均不具有內含子。

圖2 IuMYB114 和IuMYB36 基因的PCR 擴增

2.2 IuMYB114 和IuMYB36 基因的序列結構分析

IuMYB114、IuMYB36 蛋白均為親水性不穩(wěn)定蛋白，均不具有跨膜結構域，均屬于MYB 超家族蛋白成員(圖3、表3)。 預測目的基因結構域發(fā)現(xiàn)IuMYB114具有MYB 轉錄因子典型的SANT 結 構域，IuMYB36不具有保守結構域(圖4)。

圖3 IuMYB114 (上)和IuMYB36 (下)基因的保守結構域

表3 IuMYB114 和IuMYB36 基因編碼蛋白的理化性質

圖4 IuMYB114 和IuMYB36 結構域

2.3 IuMYB114 和IuMYB36 基因編碼蛋白的二、三級結構預測

IuMYB114 和IuMYB36 蛋白的二級結構均以無規(guī)則卷曲為主，說明其結構不穩(wěn)定。 蛋白質三級結構預測結果(圖5)顯示，IuMYB114 是以轉錄因子WER 為模型，屬于R2R3 型亞類成員，覆蓋率為53.19%，序列相似性為47%；IuMYB36 是以涂層蛋白為模型，覆蓋率為27.08%，序列相似性為33%。

圖5 IuMYB114 和IuMYB36 蛋白三級結構預測

2.4 IuMYB114 和IuMYB36 的氨基酸序列同源性分析及系統(tǒng)進化分析

分別將IuMYB114 和IuMYB36 的氨基酸序列與其它植物的氨基酸序列進行對比，結果(圖6)顯示，IuMYB114 與芝麻(XP＿011100166.1)、歐洲李(XP＿034208044.1)等物種的MYB 同源性較高，同源性均達64%左右；IuMYB36 與常綠越橘(KAH7854734.1)、黃芩(AKA59789.1)和藍果樹(KAA8518714.1)等的同源性較高，同源性均在50%左右。 通過構建系統(tǒng)進化樹，發(fā)現(xiàn)IuMYB114和IuMYB36 均分別與各自的同源序列聚在一起，各處在兩個分支，與IuMYB114 親緣關系最近的是旱金蓮，而與IuMYB36 親緣關系最近的是除蟲菊(圖7)。

圖6 IuMYB114(A)和IuMYB36(B)的同源氨基酸序列比對

圖7 基于IuMYB114 和IuMYB36 氨基酸序列構建的系統(tǒng)發(fā)育樹

2.5 IuMYB114 和IuMYB36 基因在滇水金鳳旗瓣中的表達模式分析

通過對兩個基因在滇水金鳳旗瓣不同部位的表達分析(圖8)發(fā)現(xiàn)，IuMYB114與IuMYB36在斑區(qū)和非斑區(qū)均有表達，但斑區(qū)表達量顯著高于非斑區(qū)，分別為非斑區(qū)的12.83 倍和9.88 倍，與轉錄組分析結果相符。 推測IuMYB114 和IuMYB36 可能在盛花期存在一定的正調控作用從而使花斑顯現(xiàn)；另外，兩個基因的表達模式和表達量相似，推測兩者在斑區(qū)具有一定的協(xié)同作用。

圖8 IuMYB114 和IuMYB36 基因在滇水金鳳旗瓣2 個部位的表達分析

3 討論與結論

MYB 轉錄因子集特異性和多樣性于一體，在花色素中發(fā)揮著重要作用，是花斑形成的關鍵基因。 最早從玉米中克隆出第一個調節(jié)花色素合成的單一型MYB基因——C1基因[24]，隨后從巨峰葡萄中克隆出VlMYBA基因，能誘導葡萄果實上紅紫色斑點的產(chǎn)生[25]。 近幾年，關于MYB 轉錄因子組織表達特異性導致花斑出現(xiàn)的機理逐漸成為研究熱點，但未見其在滇水金鳳花斑形成中作用的相關研究報道。

本研究成功從滇水金鳳花器官中克隆得到2個cDNA 全長分別為315 bp 和876 bp、分別編碼104 個和151 個氨基酸的MYB基因，命名為Iu-MYB114和IuMYB36。 經(jīng)生物信息學分析，Iu-MYB114和IuMYB36均沒有內含子，這與從蓖麻莖稈中克隆的RcMYB114[26]、從丹參中克隆的SmMYB36[27]分析結果一致，且二者均屬于親水不穩(wěn)定蛋白，二級結構均以不穩(wěn)定的無規(guī)則卷曲為主；IuMYB114 三級結構以轉錄因子WER 為模型，屬于R2R3 型亞類，且具有典型的MYB 轉錄因子SANT 結構域，這與在甘藍型油菜中克隆的RFBnMYB114的分析結果一致[28]；IuMYB36不具有保守結構域，與丹參SmMYB36基因[27]的分析結果一致，推測其屬于新的R2R3－MYB 轉錄因子亞組[29]。 同源性分析發(fā)現(xiàn)，滇水金鳳的Iu-MYB114 與芝麻(Sesamum indicum) 和 歐洲李(Prunus dulcis)等的MYB 同源性均達64%左右，IuMYB36 與常綠越橘(Vaccinium darrowii)、黃芩(Scutellaria playfairii)等的MYB 同源性均達50%左右；系統(tǒng)進化分析發(fā)現(xiàn)IuMYB114 和IuMYB36各處在兩個分支。

植物花斑的出現(xiàn)主要依賴花青素在時間和空間上的受限沉積，由于MYB基因是調節(jié)花青素合成和積累的關鍵因子，研究MYB基因對花斑形成的作用具有重要意義。 研究表明，在牡丹中R2R3－MYB 轉錄因子PsMYB12 協(xié)同WD40 和bHLH 轉錄因子，共同調控牡丹花瓣底部花斑的形成[30]；文心蘭中得到的OgMYB1基因通過控制OgCHI和OgDFR兩個結構基因的表達進而調控花斑的形成[31]；外源αMYB基因對中國水仙花色素苷代謝途徑產(chǎn)生了影響，且誘導了結構基因的表達，使其花瓣發(fā)生了顏色的變化[32]；PcMYB114是導致梨果皮紅色突變的關鍵基因[33]；AT-MYB114作為激活劑參與了擬南芥中原花青素、類黃酮或花青素的生物合成[34]；PyMYB114與PybHLH3、PyERF3之間相互作用形成調控復合體，促進梨果實花青素的積累[35]；MdMYB114是蘋果花青素生物合成和運輸?shù)姆e極調節(jié)因子[36]。結合以往研究，推測IuMYB114可能與滇水金鳳的花色苷生物合成有關。 另外，有研究表明，CsMYB36可以調節(jié)黃瓜的黃綠色果皮著色[37]；SlMYB36－1基因過表達后對番茄根系內皮層凱氏帶形成具有影響[38]；MYB36控制擬南芥?zhèn)雀吔纾瑥亩{控根從增殖到分化的轉變[39，40]。本研究發(fā)現(xiàn)IuMYB114和IuMYB36在滇水金鳳旗瓣斑區(qū)和非斑區(qū)均有表達，但在斑區(qū)的表達量顯著高于非斑區(qū)，這與MYB114L在‘鳳丹’牡丹花色變化過程中的表達分析結果[41]一致。 推測Iu-MYB114可能在滇水金鳳盛花期的斑區(qū)中起一定的正調控作用，而IuMYB36可能在斑區(qū)顏色形成中發(fā)揮著一種新功能，從而對花斑的形成具有一定調控作用；因二者表達模式和表達量相似，推測它們具有一定的協(xié)同作用，但具體功能尚有待進一步試驗驗證。

本研究結果可為后續(xù)進行IuMYB114和Iu-MYB36基因的功能研究、探究其對滇水金鳳花斑形成的調控機制以及鳳仙花花色改良等提供一定的理論依據(jù)。