蕓薹屬植物MYBL2基因的克隆及其在A、B、C基因組中的PCR鑒別

鄒婷，劉麗莉，向建華，周定港，吳金鋒，李莓，李寶，張大為，嚴明理,

蕓薹屬植物MYBL2基因的克隆及其在A、B、C基因組中的PCR鑒別

鄒婷1，劉麗莉1，向建華1，周定港1，吳金鋒1，李莓2，李寶2，張大為1，嚴明理1,2

1湖南科技大學生命科學與健康學院/經濟作物遺傳改良與綜合利用湖南省重點實驗室，湖南湘潭 410201；2湖南省作物研究所，長沙 410125

【目的】MYBL2負調控擬南芥花青素和原花青素的生物合成。從蕓薹屬6個物種的不同葉色材料中克隆MYBL2基因，分析其序列和表達模式，探究其在蕓薹屬植物花青素生物合成途徑中的功能，為油菜的品質、抗逆性、觀賞性等性狀改良提供參考。【方法】以蕓薹屬6個物種19份供試材料的總DNA為模板，同源克隆MYBL2基因，并進行多序列比對和進化樹分析；對白菜、甘藍型油菜、芥菜型油菜以及埃塞俄比亞芥的紫葉材料進行遮光處理，結合轉錄組和qRT-PCR分析MYBL2基因表達水平；對甘藍、甘藍型油菜、芥菜型油菜以及埃塞俄比亞芥紫、綠葉材料進行qRT-PCR分析MYBL2基因表達水平；根據克隆和序列的核苷酸變異位點設計特異性引物，開發能夠區分MYBL2基因組來源的PCR標記。【結果】克隆獲得和各9個同源基因共56個拷貝。其中，為首次獲得，編碼區序列全長為867 bp，包含2個內含子，分別為168和102 bp，編碼198個氨基酸，分子量為22.69 kD，等電點（pi）為8.72。序列比對和進化分析表明，來源于B基因組。蕓薹屬6個物種和同源基因中，僅、和在不同葉色材料中存在序列差異。經遮光處理后，紫葉材料葉色變淺，在白菜紫寶5號中，和表達量分別為未遮光部分的0.7和0.4倍；在紫葉白花甘藍型油菜中，、、和表達量分別為未遮光部分的0.4、0.5、0.4和0.4倍；在紫葉芥中，、、和表達量分別為未遮光部分的0.4、0.3、0.4和0.2倍；在紫稈埃芥中，、和表達量分別為未遮光部分的0.3、0.4和0.5倍，而表達量為未遮光部分的2.4倍。對比蕓薹屬不同葉色材料MYBL2基因表達情況，結果表明，除了羽衣甘藍，在紫葉材料中MYBL2基因大部分同源基因的表達量均高于綠葉。在羽衣甘藍中，綠葉羽衣甘藍和的表達量分別為紫葉羽衣甘藍的2.5和3.5倍；在甘藍型油菜中，紫葉白花、和的表達量分別為綠葉白花的7.5、8.6和26.0倍，而綠葉白花的表達量為紫葉白花的13.0倍；在芥菜型油菜中，紫葉芥、、和的表達量分別為四川黃籽的8.3、11.8、23.2和14.6倍；在埃塞俄比亞芥中，紫稈埃芥、的表達量分別為W-BCDH76的7.1和27.6倍，而W-BCDH76的和的表達量則分別為紫稈埃芥的2.8和5.0倍。依據克隆的基因序列設計出5對引物，可以有效鑒別蕓薹屬植物MYBL2基因的A、B、C基因組來源。【結論】蕓薹屬紫葉材料MYBL2基因的表達與光照密切相關，而且參與花青素生物合成的調控機制與擬南芥MYBL2負調控花青素生物合成的機制有所不同。

蕓薹屬植物；MYBL2基因；同源克隆；基因表達；基因組PCR鑒別

0 引言

【研究意義】油菜不僅是重要的油料作物，還是重要的膳食蔬菜，在全世界范圍內都有廣泛的栽培和食用，其豐富的花色和葉色具有觀賞價值。花青素（anthocyanidins）和原花青素（proanthocyanidins，PAs）是植物體內廣泛存在的色素，屬于類黃酮合成途徑中的次生代謝物，是影響植物花色、葉色、莖稈顏色和籽粒顏色形成的重要因素[1]，可以幫助植物本身應對干旱、蟲害、凍害和吸引昆蟲授粉等，蔬用和油用則具有抗氧化、預防心血管和腫瘤疾病等功效[2-3]。植物花青素和原花青素的生物合成途徑非常復雜，涉及很多關鍵酶、轉運蛋白和轉錄因子的參與[4-5]，探究并挖掘花青素和原花青素合成的關鍵基因，對油菜的品質、抗逆性、觀賞性等性狀的改良具有重要的意義。【前人研究進展】MYBL2（myeloblastosis family transcription factor like-2）轉錄因子屬于R3-MYB蛋白家族，參與植物抗逆境脅迫和類黃酮合成等重要生物過程的調控[6-7]。擬南芥負調控花青素的生物合成，且可能參與了花青素合成的反饋抑制，AtMYBL2通過與TT8（bHLH類轉錄因子）蛋白相互作用形成轉錄抑制復合物，該復合物能夠結合（編碼二氫黃酮醇4-還原酶，Dihydorflavonol 4-reductase）啟動子，抑制和的轉錄，從而抑制花色苷的積累[8-11]。在敲除系中，花青素合成關鍵基因和調控基因的表達水平上調，花青素含量顯著上升[12]。Song等[13]發現的表達與花青素含量呈負相關，紫甘藍（）葉片的紫色是由于啟動子替換或基因序列缺失導致的表達缺失所致。然而，Mushtaq等[14]結合qRT-PCR和轉錄組分析了紫、綠葉白菜（），結果表明，在紫葉中的表達量顯著高于綠葉。Rameneni等[15]結合qRT-PCR和轉錄組分析發現在紅紫色白菜（reddish purple Chinese cabbage，RPCC）中的轉錄水平高于綠色白菜。Ye等[16]發現與甘藍型油菜黃色和白色花瓣相比，、和是杏色花瓣和粉色花瓣中唯一顯著上調的3個花青素調節基因，這些結果與黃色和白色花瓣中花青素含量較低相一致，與花青素積累呈顯著正相關。因此，MYBL2基因可能是蕓薹屬植物花青素生物合成途徑中重要的調控基因，但MYBL2基因是發揮正調控還是負調控功能需要進一步驗證。有報道表明基因組原位雜交能夠辨別蕓薹屬雜種后代的A、B、C基因組[17-18]，也有根據基因上的酶切位點來鑒定A、C組的報道[19]。基因組原位雜交操作比較復雜，有多態性且能被酶切的位點相對有限，雖然有利用PCR和FISH等方法檢測B組染色體上的特異序列的報道，但這樣的特異序列也非常有限[20]。Drenkard等[21]以擬南芥為研究材料開發了一種改進的等位基因特異性PCR程序，設計引物時引入錯配堿基形成錯配，用于檢測單核苷酸多態性（single nucleotide polymorphism，SNP）位點。YAN等[22-25]根據數據庫已有和克隆得到的蕓薹屬植物、和的DNA序列分別在A、B、C基因組中特有的核苷酸多態位點進行引物設計，開發出能夠簡單、快速地檢測基因的多態性和等位基因的基因組來源的PCR標記。【本研究切入點】蕓薹屬植物基因的克隆和功能研究主要以白菜、甘藍、甘藍型油菜為主，在黑芥、芥菜型油菜和埃塞俄比亞芥中研究較少，蕓薹屬6個物種花青素生物合成途徑中MYBL2基因的功能尚未全部明確，關于蕓薹屬6個物種MYBL2基因序列的變異和多態性的研究也尚未見報道。【擬解決的關鍵問題】本研究基于已公布的和的DNA序列設計引物，以蕓薹屬6個物種不同葉色的材料為對象，同源克隆MYBL2基因，通過比較和在不同葉色材料中的序列差異和基因表達情況，以探討其與葉色之間的關系。通過分析蕓薹屬6個物種不同葉色材料中和特異的SNP位點，開發出能夠識別來自A、B、C基因組的MYBL2基因的PCR標記，為進一步了解其功能及起源進化關系奠定基礎。

1 材料與方法

1.1 試驗材料

蕓薹屬6個物種19份材料，包括白菜（，2n=20，AA）早熟40、綠葉白菜、白菜全紫株、紫寶5號；黑芥（，2n=16，BB）Giebra（由華中農業大學李再云課題組惠贈）；甘藍（，2n=18，CC）綠葉羽衣甘藍、紫葉羽衣甘藍；芥菜型油菜（，2n=36，AABB）四川黃籽、李G19、紫葉芥、Z44紫葉黃籽、紫葉cd30；甘藍型油菜（，2n=38，AACC）湘油787、ZSL409、綠葉白花和紫葉白花；埃塞俄比亞芥（，2n=34，BBCC）W-BCDH76、黃埃和紫稈埃芥，用途詳見電子附表1。

1.2 DNA提取和cDNA的合成

采用天根生化科技（北京）有限公司的試劑盒提取葉片的總DNA和總RNA，運用北京擎科生物科技有限公司的反轉錄試劑盒合成cDNA，1.0%瓊脂糖凝膠檢測DNA和RNA的完整性，核酸測定儀檢測其濃度和純度，將DNA和cDNA稀釋至50 ng·μL-1。

1.3 蕓薹屬MYBL2基因的克隆和序列分析

以19份供試材料葉片的總DNA為模板進行基因編碼區序列克隆，參照NCBI（https://www.ncbi.nlm.nih. gov/）、BnTIR（http://yanglab.hzau.edu.cn/BnTIR）和BRAD（http://brassicadb.cn）等數據庫已有的和序列設計引物（電子附表2），由北京擎科生物科技有限公司合成。PCR反應體系為2×Rapid Taq master mix 10 μL、10 mmol?L-1正向引物1 μL、10 mmol?L-1反向引物1 μL、50 ng·μL-1DNA模板2 μL，無菌去離子水補足至20 μL。PCR擴增程序為95℃ 5 min；95℃ 30 s，55℃ 30 s，72℃ 1 min，35個循環；72℃ 10 min。用1.0%瓊脂糖凝膠電泳檢測PCR產物，回收后，與pMD19-T載體連接，并轉化大腸桿菌DH5α感受態細胞，利用Amp抗性和X-gal/IPTG進行藍白斑篩選。選取6個陽性克隆送生工生物工程（上海）股份有限公司測序，采用DNAMAN軟件分析蕓薹屬不同葉色材料測序所得序列，結合數據庫已有的序列進行生物信息學分析，結合GenBank數據庫中下載與蕓薹屬親緣關系較近物種的MYBL2蛋白質序列進行比對，用MEGA11.0軟件構建系統發育樹。

1.4 蕓薹屬MYBL2基因轉錄水平的分析

將紫寶5號、紫葉白花甘藍型油菜、紫葉芥和紫稈埃芥材料的紫葉進行遮光處理（3個生物學重復）2周，取葉片。將紫葉羽衣甘藍、綠葉羽衣甘藍、紫葉白花甘藍型油菜、綠葉白花甘藍型油菜、紫葉芥、四川黃籽、W-BCDH76和紫稈埃芥材料進行土培，6葉期時取新鮮葉片。對上述葉片提取總RNA，反轉錄成cDNA。并將紫葉遮光和未遮光材料送上海美吉生物醫藥科技有限公司進行轉錄組測序分析。根據測序所得蕓薹屬6個物種MYBL2基因序列差異設計引物（電子附表2）。利用實時熒光定量PCR（qRT-PCR）檢測蕓薹屬MYBL2基因轉錄水平，以為內參基因，qRT-PCR反應體系為熒光染料SYBR 10 μL、10 mmol?L-1正向引物1 μL、10 mmol?L-1反向引物1 μL、50 ng·μL-1cDNA模板1 μL，無菌去離子水補足至20 μL。qRT-PCR擴增程序為95℃ 5 min；95℃ 10 s，55℃ 30 s，35個循環；65—95℃熔解曲線，每0.5℃讀板1次檢測擴增產物的特異性。用2-ΔΔCT算法分析目標基因的相對表達量。

1.5 MYBL2-1和MYBL2-2等位基因的特異PCR檢測

根據蕓薹屬6個物種和多態位點和基因特性，通過多次設計引物和優化PCR條件，設計5對特異性引物，能夠區分來自A、B、C基因組的MYBL2基因（電子附表2）。PCR反應體系為2×Rapid Taq master mix 10 μL、10 mmol?L-1正向引物1 μL、10 mmol?L-1反向引物1 μL、50 ng·μL-1DNA模板2 μL，無菌去離子水補足至20 μL。PCR擴增程序為95℃ 5 min；95℃ 30 s，52℃ 30 s，72℃ 1 min，33個循環；72℃ 5 min。1.5%瓊脂糖凝膠電泳檢測。

2 結果

2.1 蕓薹屬植物MYBL2基因的克隆

以蕓薹屬6個物種19份供試材料葉片總DNA為模板進行編碼區序列的克隆，共克隆獲得和各9個同源基因，共計56段序列（電子附表1，電子附圖1）。其中，用引物和從埃塞俄比亞芥W-BCDH76材料中擴增獲得1條約1 000 bp大小的片段。測序結果表明，該序列為867 bp（電子附圖2），包含2個內含子，分別為168和102 bp，且在參考基因組上未比對到相應基因。因此，該片段系首次獲得，命名為。對序列進行比對，發現該序列與參考基因組中和序列相似性最高，分別為100%和98.39%；與BRAD數據庫埃塞俄比亞芥另一個序列相似性為88.24%。因此，所獲得的確定為B基因組來源的。

2.2 蕓薹屬植物MYBL2-1氨基酸序列的推導與分子進化分析

根據埃塞俄比亞芥W-BCDH76材料的cDNA序列，用MEGA11.0和DNAMAN8.0軟件推導，BcaMYBL2-1共編碼198個氨基酸（圖1），分子量為22.69 kD，等電點（pI）為8.72，屬于R3-MYB家族。BLASTp分析結果表明，BcaMYBL2-1蛋白序列與黑芥BniB03.MYBL2-1的相似性為100%；與芥菜型油菜BjuB03.MYBL2-1的相似性為95.96%。經多序列比對，BcaMYBL2-1與擬南芥、蘿卜以及蕓薹屬6個物種MYBL2序列具有較高的保守性（電子附圖3）。利用MEGA11.0構建進化樹，發現來自不同材料的MYBL2-1與MYBL2-2分別聚在不同進化枝，且來源相同基因組（A、B或C基因組）的蛋白聚在一起，表明與的分化在蕓薹屬基因組分化和異源四倍化之前（圖2）。BcaMYBL2-1與來自B基因組的BniB03.MYBL2-1和BjuB03.MYBL2-1的進化關系較近，進一步確定了其來源于B基因組。

BniB03 G：黑芥（Giebra）；Bca G：WBCDH76；Bca P：紫稈埃芥；BjuB03 SY：四川黃籽；BjuB03 PM：紫葉芥

Bra G：綠葉白菜；Bra P1：白菜全紫株；Bra P2：紫寶5號；Bni G：Giebra；Bol G：綠葉羽衣甘藍；Bol P：紫葉羽衣甘藍；SY：四川黃籽；PM：紫葉芥；Bna G：綠葉白花；Bna P：紫葉白花；Bca G1：W-BCDH76；Bca G2：黃埃；Bca P：紫稈埃芥；At：擬南芥；Rs：蘿卜

2.3 蕓薹屬植物不同葉色材料MYBL2基因的序列差異

通過對克隆MYBL2基因進行分析，結果表明，大部分克隆基因與參考基因組中序列完全一致（50/56），僅在部分材料中、和與參考基因組存在差異（表1）。比較蕓薹屬植物不同葉色材料MYBL2基因序列，與綠葉白菜相比，白菜全紫株外顯子區域存在一段長為61 bp的片段插入（表1、電子附圖4）。在芥菜型油菜中，不同葉色材料的外顯子區域均有該61 bp片段的插入，而甘藍型油菜中則均無該片段的插入（表1）。比較綠葉埃芥（W-BCDH76）和紫稈埃芥的編碼區序列，共發現9個核苷酸差異位點差異（電子附圖2），導致4處氨基酸變化（表1和圖1）；比較綠葉和紫葉羽衣甘藍的序列，共發現11個核苷酸位點變異，導致3個氨基酸變化（表1，電子附圖5）。但和在不同葉色的甘藍型油菜和埃塞俄比亞芥材料中的均無差異（電子附圖5）。此外，8個拷貝（黑芥暫無紫葉材料）在不同葉色材料中編碼區序列均無差異（電子附圖6）。

表1 蕓薹屬不同葉色材料MYBL2基因序列差異位點

61 bp片段為TTGTTTCAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAGGGAATCGATCCAACT

The 61 bp fragment is TTGTTTCAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAGGGAATCGATCCAACT

2.4 遮光處理下蕓薹屬MYBL2基因的表達

白菜、甘藍型油菜、芥菜型油菜和埃塞俄比亞芥的紫葉材料經遮光處理2周后進行轉錄組測序與qRT-PCR分析，結果表明，紫葉材料經過遮光處理后葉色變淺，MYBL2基因表達普遍下調（圖3）。在紫寶5號、紫葉白花甘藍型油菜和紫葉芥中，與未遮光部分相比，MYBL2基因在遮光后均下調表達。在白菜紫寶5號中，和表達量分別為未遮光部分的0.7和0.4倍（圖3-B）；在紫葉白花甘藍型油菜中，、、和表達量分別為未遮光部分的0.4、0.5、0.4和0.4倍（圖3-C）；在紫葉芥中，、、和表達量分別為未遮光部分的0.4、0.3、0.4和0.2倍（圖3-D）。在紫稈埃芥中，、和表達量分別為未遮光部分的0.3、0.4和0.5倍，而表達量為未遮光部分的2.4倍（圖3-E）。對紫葉材料遮光前后轉錄組數據進行分析，結果表明，花青素代謝途徑中基因在遮光后平均表達量相比未遮光部分均下調（圖4）。其中，、、等基因在蕓薹屬植物多個材料遮光后均下調表達，因此，推測上述基因與MYBL2基因一樣，與光照密切相關，且遮光后被抑制表達。

A：蕓薹屬紫葉材料遮光處理后的表型；B：紫寶5號；C：紫葉白花甘藍型油菜；D：紫葉芥；E：紫稈埃芥；sP：紫葉未遮光部分；sG：紫葉遮光部分

2.5 蕓薹屬植物不同葉色材料MYBL2基因的表達

比較甘藍、甘藍型油菜、芥菜型油菜和埃塞俄比亞芥紫綠葉材料中MYBL2基因的表達，結果表明，除了羽衣甘藍，紫葉材料中MYBL2基因大部分同源基因均比綠葉上調表達。在羽衣甘藍中，綠葉羽衣甘藍和表達量分別為紫葉羽衣甘藍的2.5和3.5倍（圖5-B）；在甘藍型油菜中，紫葉白花、和表達量分別為綠葉白花的7.5、8.6和26.0倍，而綠葉白花表達量為紫葉白花的13.0倍（圖5-C）；在芥菜型油菜中，紫葉芥、、和表達量分別為四川黃籽的8.3、11.8、23.2和14.6倍（圖5-D）；在埃塞俄比亞芥中，紫稈埃芥、表達量分別為W-BCDH76的7.1和27.6倍，而W-BCDH76的和表達量則分別為紫稈埃芥的2.8和5.0倍（圖5-E）。對蕓薹屬紫綠葉材料中花青素合成相關（甘藍16個，甘藍型油菜31個，芥菜型油菜31個，埃芥32個）基因表達進行分析（電子附圖7），結果表明，甘藍、甘藍型油菜、芥菜型油菜和埃塞俄比亞芥紫綠葉材料中花青素合成相關基因表達模式存在差異，推測上述材料紫葉形成機制存在差異。

2.6 MYBL2-1和MYBL2-2等位基因的特異PCR檢測

根據蕓薹屬6個物種和序列的SNP位點，設計5對引物以區分來自A、B、C基因組MYBL2基因（圖6）。在9份供試材料中，引物可從具有A基因組材料中擴增出預期的目標條帶；引物B1、B2和B均可從B基因組材料中產生特異的目標條帶；引物對僅從具有C基因組的材料中產生預期的目標條帶。通過上述引物，可以有效地鑒別蕓薹屬植物MYBL2的基因組來源。

3 討論

3.1 羽衣甘藍不同葉色材料MYBL2基因的表達分析

比較不同葉色蕓薹屬植物MYBL2基因，結果表明，僅、和編碼區序列在紫、綠葉材料中存在差異。Song等[13]對紫葉甘藍位于C06染色體上的的啟動子區序列和編碼區序列研究，結果表明，無論是啟動子的替換，還是的全缺失，都可能導致甘藍紫葉性狀的出現。LIU等[26]在研究同一株觀賞性羽衣甘藍葉綠素、類胡蘿卜素和花青素色素沉著時，發現MYBL2基因在粉紅色葉片、淺粉紅色葉片和綠色粉色雜色葉片中依次上調表達。HAN等[27]在對甘藍花青素合成相關基因進行全基因組表征，發現和在紫葉甘藍（材料：15Z-P）中的表達量顯著低于綠葉羽衣甘藍（材料：120QY-G）。本研究綠葉羽衣甘藍和紫葉羽衣甘藍的編碼區序列存在11個核苷酸位點變異，會導致3個氨基酸變化，且紫葉材料中和表達量均相比綠葉低，但本研究所用甘藍材料編碼區序列的差異是否會導致葉色變化還有待進一步研究。

3.2 白菜、甘藍型油菜、芥菜型油菜和埃塞俄比亞芥MYBL2基因的表達分析

Zhao等[28]發現在MYBL2基因上游通過抑制MYBL2基因正向調控蔗糖誘導的擬南芥花青素生物合成，MYBL2基因過表達顯著抑制過表達幼苗的花青素合成。Xie等[29]在研究茉莉酸和赤霉素對花青素合成的影響過程中發現MBW復合物活性被MYBL2和JAZ家族蛋白抑制，它們分別與bHLH和MYB/bHLH競爭性結合，DELLA蛋白直接隔離MYBL2和JAZ抑制子，導致bHLH/MYB亞基的釋放，隨后形成活性MBW復合物，激活花青素生物合成途徑，推測JAZ-DELLA-MYBL2模塊在非生物脅迫誘導的花青素合成中也發揮了重要作用。Nguyen等[30]研究發現在高強度光照情況下，HY5（bZIP家族）直接與啟動子結合抑制擬南芥MYBL2基因的表達增加花青素的積累。然而，在蕓薹屬植物研究中，IMTIAZ等[31]發現在釩脅迫下紫葉芥和綠葉芥花青素含量降低，花青素生物合成基因在紫葉芥中高表達，尤其是釩脅迫下的、和，釩濃度與花青素生物合成之間呈負相關，結果表明，、和的上調可以觸發釩脅迫下花色苷的積累。KHUSNUTDINOV等[32]發現在芥菜型油菜的3個紫色品種（Vitamin、Red velvet和Freckle）中的表達量均顯著高于綠色品種 。Ye等[16]發現與甘藍型油菜黃色和白色花瓣相比，、和是杏色花瓣和粉色花瓣中唯一顯著上調的3個花青素調節基因，這些結果與黃色和白色花瓣中花青素含量較低相一致，與花青素積累呈顯著正相關。本研究白菜、甘藍型油菜、芥菜型油菜和埃塞俄比亞芥的紫葉材料經遮光處理后葉色變淺，MYBL2基因表達量普遍下降，表明MYBL2基因的表達水平與光照密切相關且可能與花青素的積累呈正相關。對甘藍型油菜、芥菜型油菜和埃塞俄比亞芥不同葉色材料MYBL2基因的表達進行分析發現，除了、和在紫葉中的表達量均低于綠葉，其余和在紫葉中的表達量均高于綠葉。因此，在蕓薹屬植物花青素合成途徑中MYBL2轉錄因子究竟發揮正調控、負調控還是拮抗功能？蕓薹屬植物與擬南芥MYBL2是否發揮著不同功能？仍需要進一步驗證。

A：紫寶5號；B：紫葉白花甘藍型油菜；C：紫葉芥；D：紫稈埃芥；sP：紫葉未遮光部分；sG：紫葉遮光部分

A：蕓薹屬不同葉色材料的表型；B（P）：紫葉羽衣甘藍；B（G）：綠葉羽衣甘藍；C（P）：紫葉白花甘藍型油菜；C（G）：甘藍型油菜綠葉白花；D（P）：紫葉芥；D（G）：四川黃籽；E（P）：紫稈埃芥；E（G）：W-BCDH76；P：紫葉材料；G：綠葉材料

A: Phenotypes ofwith different leaf colors; B(P): Green leaf kale; B(G): Purple leaf kale; C(P):with purple leaves and white flowers; C(G):with green leaves and white flowers; D(P): Purple leaf mustard; D(G): Sichuan Yellow; E(P):with purple leaf; E(G): W-BCDH76; P: Materials with purple leaf; G: Materials with green leaf

圖5 蕓薹屬不同葉色材料MYBL2基因的表達

Fig. 5 The expression ofwith different leaf color

M：標準分子量；A：白菜早熟40；B：黑芥；C：羽衣甘藍；AB1：紫葉芥；AB2：李G19；AC1：湘油787；AC2：紫葉白花甘藍型油菜；BC1：黃埃；BC2：紫稈埃芥；N：空白對照；ACT：內參基因對照

3.3 蕓薹屬植物MYBL2基因A、B、C基因組的鑒別

蕓薹屬植物中，高花青素含量的材料多見在白菜、甘藍和芥菜型油菜，而大規模種植的甘藍型油菜其高花青素含量材料相對較少，且部分為遠緣雜交獲得[33-35]。本研究基于蕓薹屬6個物種和序列多態性位點設計的特異性引物，通過等位基因特異PCR可以準確區分來自A、B基因組的基因和來自B、C基因組的。在解析蕓薹屬MYBL2基因功能后，上述標記可為利用遠緣雜交在A、B、C基因組間轉移MYBL2基因，創制高花青素含量的甘藍型油菜提供參考。

4 結論

蕓薹屬紫葉材料經遮光處理后，MYBL2基因表達量普遍下降，葉色變淺，表明該基因的表達與光照密切相關。除了羽衣甘藍，發現在甘藍型油菜、芥菜型油菜和埃塞俄比亞芥的紫葉材料中MYBL2基因的大部分同源基因的表達水平均高于綠葉，表明蕓薹屬植物中MYBL2基因參與花青素生物合成的調控機制與擬南芥有所不同。

[1] YAN M L, Ding S P, Liu L L, YIN X M, SHU J B. Cloning and expression analysis of an anthocyanidin synthase gene homologue fromJournal of genetics, 2014, 93(2): 513-516.

[2] IWASHINA T. Contribution to flower colors of flavonoids including anthocyanins: a review. Natural product communications, 2015, 10(3): 529-544.

[3] BAI S L, TAO R Y, TANG Y X, YIN L, MA Y J, NI J B, YAN X H, YANG Q S, WU Z Y, ZENG Y L, TENG Y W. BBX16, a B-box protein, positively regulates light-induced anthocyanin accumulation by activating

[4] HE Q, WU J Q, XUE Y H, ZHAO W B, LI R, ZHANG L G. The novel gene, located on chromosome A07, with a short intron 1 controls the purple-head trait of Chinese cabbage (L.). Horticulture research, 2020, 7: 97.

[5] ZHANG D W, LIU L L, ZHOU D G, LIU X J, LIU Z S, YAN M L. Genome-wide identification and expression analysis of anthocyanin biosynthetic genes inJournal of Integrative Agriculture, 2020, 19(5): 1250-1260.

[6] OGATA K, KANEI-ISHII C, SASAKI M, HATANAKA H, NAGADOI A, ENARI M, NAKAMURA H, NISHIMURA Y, ISHII S, SARAI A. The cavity in the hydrophobic core of MYB DNA-binding domain is reserved for DNA recognition and trans-activation. Nature structural biology, 1996, 3(2): 178-187.

[7] YE H X, LI L, GUO H Q, YIN Y H. MYBL2 is a substrate of GSK3-like kinase BIN2 and acts as a corepressor of BES1 in brassinosteroid signaling pathway inProceedings of the National Academy of Sciences of the United States of America, 2012, 109(49): 20142-20147.

[8] MATSUI K, UMEMURA Y, OHME-TAKAGI M. AtMYBL2, a protein with a single MYB domain, acts as a negative regulator of anthocyanin biosynthesis in. The Plant Journal, 2008, 55(6): 954-967.

[9] ZHANG X, ZHANG K, WU J, GUO N, LIANG J, WANG X, CHENG F. QTL-Seq and sequence assembly rapidly mapped the genefor the purple trait inica. Scientific reports, 2020, 10: 2328.

[10] WANG Y L, Wang Y Q, Song Z Q, Zhang H Y. Repression of

by both microRNA858a and HY5 leads to the activation of anthocyanin biosynthetic pathway in. Molecular plant, 2016, 9(10): 1395-1405.

[11] XU W J, GRAIN D, BOBET S, LE GOURRIEREC J, THéVENIN J, KELEMEN Z, LEPINIEC L, DUBOS C. Complexity and robustness of the flavonoid transcriptional regulatory network revealed by

comprehensive analyses of MYB-bHLH-WDR complexes and their targets inseed. New phytologist, 2014, 202(1): 132-144.

[12] DUBOS C, LE GOURRIEREC J, BAUDRY A, HUEP G, LANET E, DEBEAUJON I, ROUTABOUL J M, ALBORESI A, WEISSHAAR B, LEPINIEC L. MYBL2 is a new regulator of flavonoid biosynthesis in. The Plant Journal, 2008, 55(6): 940-953.

[13] SONG H, YI H, LEE M, HAN C T, LEE J, KIM H, PARK J I, NOU I S, KIM S J, HUR Y. Purpleis due to the loss ofexpression. BMC plant biology, 2018, 18(1): 82.

[14] MUSHTAQ M A, PAN Q, CHEN D Z, ZHANG Q H, GE X H, LI Z Y. Comparative leaves transcriptome analysis emphasizing on accumulation of anthocyanins in: molecular regulation and potential interaction with photosynthesis. Frontiers in plant science, 2016, 7: 311.

[15] RAMENENI J J, CHOI S R, CHHAPEKAR S S, KIM M S, SINGH S, YI S Y, OH S H, KIM H, LEE C Y, OH M H, LEE J, KWON O H, PARK S U, KIM S J, LIM Y P. Red Chinese cabbage transcriptome analysis reveals structural genes and multiple transcription factors regulating reddish purple color. International journal of molecular sciences, 2020, 21(8): 2901.

[16] YE S H, HUA S J, MA T T, MA X W, CHEN Y P, WU L M, ZHAO L, YI B, MA C Z, TU J X, SHEN J X, FU T D, WEN J. Genetic and multi-omics analyses revealas the key gene conferring anthocyanin-based color inflowers. Journal of Experimental Botany, 2022, 73(19):6630-6645.

[17] 孔芳, 蔣金金, 吳磊, 王幼平. 利用原位雜交及CAPS標記分析蕓薹屬A、B和C基因組間的關系. 作物學報,2008, 34(7): 1188-1192.

KONG F, JIANG J J, WU L, WANG Y P. relationships among genome A, B and C revealed by FISH and CAPS. Acta Agronomica Sinica, 2008, 34(7): 1188-1192. (in Chinese)

[18] 李宗蕓, 伍曉明, 王秀琴, 宋運淳. 甘藍與蕓薹屬5個近緣物種的基因組原位雜交分析. 中國油料作物學報, 2003, 25(4): 16-19, F003.

LI Z Y, WU X M, WANG X Q, SONG Y C. Genomic in situ hybridization analysis ofand 5 relatedspeciesChinese Journal of Oil Crop Scieves,2003, 25(4): 16-19, F003. (in Chinese)

[19] WU Y H, XIAO L, WU G, LU C M. Cloning of fatty acid elongase1 gene and molecular identification of A and C genome inspecies. Science in China. Series C, Life sciences, 2007, 50(3): 343-349.

[20] SCHELFHOUT C J, SNOWDON R, COWLING W A, WROTH J M. A PCR based B-genome-specific marker inspecies. Theoretical and applied genetics, 2004, 109(5): 917-921.

[21] DRENKARD E, RICHTER B G, ROZEN S, STUTIUS L M, ANGELL N A, MINDRINOS M, CHO R J, OEFNER P J, DAVIS R W, AUSUBEL F M. A simple procedure for the analysis of single nucleotide polymorphisms facilitates map-based cloning inPlant physiology, 2000, 124(4): 1483-1492.

[22] YAN M L, LIU X J, GUAN C Y, Liu L L, Lu Y, Liu Z S . Cloning and SNP analysis ofgene incaActa Agronomica Sinica, 2010, 36(10): 1634-1641.

[23] YAN M L, LIU X J, GUAN C Y, CHEN X B, LIU Z S. Cloning and expression analysis of an anthocyanidin synthase gene homolog from.Molecular Breeding, 2011, 28(3): 313-322.

[24] 嚴明理, 劉麗莉, 向建華, 丁素萍, 舒佳賓, 孫嬋. 黑芥基因的克隆和在蕓薹屬B基因組的PCR鑒別. 中國油料作物學報, 2013, 35(5): 484-490.

YAN M L, LIU L L, XIANG J H, DING S P, SHU J B, SUN C. Cloning ofgene inand PCR identification froma B-genomes. Chinese Journal of Oil Crop Sciences,2013, 35(5): 484-490. (in Chinese)

[25] YAN M L, Liu X J, Guan C Y, Liu L L, Xiang J H, Lu Y, Liu Z S. Cloning ofgene and PCR identification of genomes A, B and C inspecies. Genetica, 2014, 142(2): 169-176.

[26] Liu Y, Feng X, Zhang Y T, Zhou F H, Zhu P F. Simultaneous changes in anthocyanin, chlorophyll, and carotenoid contents produce green variegation in pink-leaved ornamental kale. BMC genomics, 2021, 22(1): 455.

[27] HAN F Q, ZHANG X L, YANG L M, ZHUANG M, ZHANG Y Y, LIU Y M, LI Z S, WANG Y, FANG Z Y, JI J L, LV H H. Genome-wide characterization and analysis of the anthocyanin biosynthetic genes in. Planta, 2021, 254(5): 92.

[28] ZHAO D M, ZHENG Y X, YANG L J, YAO Z Y, CHENG J F, ZHANG F, JIANG H Y, LIU D. The transcription factor AtGLK1 acts upstream of MYBL2 to genetically regulate sucrose-induced anthocyanin biosynthesis inBMC plant biology, 2021, 21(1): 242.

[29] XIE Y, TAN H J, MA Z X, HUANG J R. DELLA proteins promote anthocyanin biosynthesis via sequestering MYBL2 and JAZ suppressorsof the MYB/bHLH/WD40 complex inMolecularplant, 2016, 9(5): 711-721.

[30] NGUYEN N H, JEONG C Y, KANG G H, YOO S D, HONG S W, LEE H. MYBD employed by HY5 increases anthocyanin accumulation via repression ofThe Plant Journal, 2015, 84(6): 1192-1205.

[31] Imtiaz M, Mushtaq M A, Nawaz M A, Ashraf M, Rizwan M S, Mehmood S, Aziz O, Rizwan M, Virk M S, Shakeel Q, Ijaz R, Androutsopoulos V P, Tsatsakis A M, Coleman M D. Physiological and anthocyanin biosynthesis genes response induced by vanadium stress in mustard genotypes with distinct photosynthetic activity. Environmental toxicology and pharmacology, 2018, 62: 20-29.

[32] Khusnutdinov E, Artyukhin A, Sharifyanova Y, Mikhaylova E V. A mutation in thegene is associated with purple pigmentation inInternational journal of molecular sciences, 2022, 23(19): 11865.

[33] HU D D, JING J J, SNOWDON R J, MASON A S, SHEN J X, MENG J L, ZOU J. Exploring the gene pool ofby genomics-based approaches. Plant Biotechnology Journal, 2021, 19(9): 1693-1712.

[34] FU W, CHEN D, PAN Q, LI F, ZHAO Z, GE X, LI Z. Production of red-flowered oilseed rape via the ectopic expression ofOvPAP2.Plant biotechnology journal, 2018, 16(2): 367-380.

[35] GOSWAMI G, NATH U K, PARK J I, HOSSAIN M R, BISWAS M K, KIM H T, KIM H R, NOU I S. Transcriptional regulation of anthocyanin biosynthesis in a high-anthocyanin resynthesizedcultivar. Journal of Biological Research (Thessalonike , Greece), 2018, 25: 19.

Cloning of MYBL2 gene fromand its PCR identification in genomes A, B and C

ZOU Ting1, LIU LiLi1, XIANG JianHua1, ZHOU DingGang1, WU JinFeng1, LI Mei2, LI Bao2, ZHANG DaWei1, YAN MingLi1,2

1College of Life Science and Health, Hunan University of Science and Technology/Hunan Provincial Key Laboratory of Genetic Improvement and Comprehensive Utilization of Cash Crops, Xiangtan 410201, Hunan;2Hunan Crop Research Institute, Changsha 410125

【Objective】In,MYBL2 negatively regulates the biosynthesis of anthocyanins and proanthocyanidinsThe MYBL2 genes from sixspecies with different leaf colors were cloned. By analyzing the sequence and expression pattern of MYBL2, the function of MYBL2 in the biosynthesis of anthocyanins inspecies was explored. 【Method】The sequences of MYBL2 from the sixspecies with different leaf colors were obtained using homology-based cloning method and multi sequence alignment and phylogenetic tree analysis were performed. The purple leaf materials of,,andwere treated with shading, and the expression level of MYBL2 gene was analyzed by transcriptome sequencing and qRT-PCR. The qRT-PCR in,,andwith purple and green leaves were also performed to evaluate the expression level of MYBL2. Based on the nucleotide variation sites of the clonedandsequences, PCR markers which could distinguish the genomic origin of alleles were developed. 【Result】A total of 56 copies of 9 homologs ofandwere cloned from 19 samples of six species of. Thegene was obtained for the first time. The total length ofsequence was 867 bp, including two introns of 168 bp and 102 bp respectively, encoding 198 amino acids, with a molecular weight of 22.69 kD and an isoelectric point (pI) of 8.72. Sequence alignment and evolutionary analysis showed thatwas derived from B genome. Among theandcopies of six species in, only,andexhibited sequence differences in different leaf color materials. After shading treatment, the leaf color of purple leaf material becomes lighter than that of the unshaded part. In Chinese cabbage Zibao 5, the expression ofandwith purple leaves and white flowers, the expression of,,and,,andwith purple leaf, the expression of,and, the results showed that the expression of most of the homologous genes of MYBL2 gene in purple leaf materials was higher than that in green leaf materials except kale. In kale, the expression ofand, the expression of,and,,andwere 8.3 times, 11.8 times, 23.2 times and 14.6 times of those in Sichuan yellow respectively. Inwith purple leaf, BcaMYBL2-1 and BcaB03.MYBL2-1 were 7.1 and 27.6 times as much as W-BCDH76, respectively. However,andgenes of W-BCDH76 were 2.8 and 5.0 times as much as those ofwith purple leaf, respectively. Five pairs of primers were obtained, which can effectively identify thefrom A, B and C genomes of. 【Conclusion】After shading treatment, MYBL2 gene expression was closely related to light. the regulation mechanism of MYBL2gene inplants involved in anthocyanin biosynthesis was different from that ofplants in which MYBL2 gene negatively regulated anthocyanin biosynthesis.

species; MYBL2gene; homologous cloning; gene expression; genomic PCR identification

10.3864/j.issn.0578-1752.2023.03.002

2022-08-31；

2022-10-24

國家自然科學基金（31971980，U19A2029）

鄒婷，Tel：15297827231；E-mail：3392648002@qq.com。通信作者張大為，E-mail：zhangdawei.hnust@foxmail.com。通信作者嚴明理，E-mail：ymljack@126.com

（責任編輯 李莉）