新型紫色大白菜15NG28后代的鑒別及其轉錄組分析

張德雙，趙 泓，張鳳蘭，余陽俊，趙岫云，于拴倉，汪維紅，蘇同兵，盧桂香

(北京市農林科學院 蔬菜研究中心，北京 100097)

大白菜(Brassica.rapaL.ssp.pekinensis)起源于中國，種質資源豐富、生態類型多樣，是名副其實的大眾化蔬菜。隨著居民生活質量的提高和家庭消費模式、人口結構的改變，大白菜富含的蛋白質、糖、氨基酸、維生素C、膳食纖維、胡蘿卜素以及花青苷等營養品質日益受到國內外學者的重視。大白菜品種類型多樣，繼黃心、桔紅心大白菜后，近幾年育種家又成功選育出紫色葉球大白菜。紫色葉球大白菜色澤艷麗、風味獨特，富含花青苷等營養成分，逐漸被消費者接納和喜愛。植物體內的花青素為水溶性色素，在自然條件下不能穩定存在，常常結合一個或多個糖苷，形成花青苷。花青苷是黃酮類次生代謝產物，在植物細胞的內質網上合成，然后運輸到液泡并積累，可以引起植物的葉片、花薹、根莖和種莢等呈現紫色、紅色、藍色和粉色等顏色，如甘藍類的紫紅色甘藍、紫色羽衣甘藍、紫色苤藍和白菜類的紅菜薹、紫色小白菜和紫紅色蕪菁等。花青苷具有吸引昆蟲傳粉和抵御低溫、紫外線傷害以及防治植物病害等多種生理功能；也具有抗氧化衰老、抗癌、抗動脈硬化，減少夜盲癥和預防近視，防止神經紊亂以及增強人體免疫力等重要的醫療保健功能；花青素還是一種安全的天然食用色素。花青素是迄今為止所發現的最強效的自由基清除劑，抗自由基氧化能力是VC的20倍、VE的50倍，并且對光、熱和氧等穩定性好。因此，花青素在園藝、醫藥、化妝及食品等方面具有重要的應用價值和研究意義[1]。

目前，紫色大白菜來源主要有4種：紫紅色芥菜[2]、紫色小白菜[3-6]、紫菜薹[7]和紫紅色蕪菁[8]。在國內，已育成紫羅蘭、紫冠一號、紫貴妃等小白菜品種和京研紫快菜等大白菜品種，但是紫色性狀均來源于紫色小白菜，其紫色僅分布在葉片的正面，而在葉片的反面和葉柄處沒有分布，且葉球剖面顏色不均勻，表現為雜色[3-4]，因此，該紫色性狀存在一定的局限性，無法廣泛利用。當前，國內還沒有利用另外3種紫色性狀育成白菜品種的報道。近年來，在國際上韓國首次培育出紫色大白菜品種紫寶(http：//www.1688.com/chanpin/D7CFB0D7B2CB.html)，該品種為細胞質不育系配制的雜交F1，紫色均勻分布在葉片的正反面、葉脈和葉柄處，且葉球剖面紫色分布均勻，顏色美觀，但是由于細胞質雄性不育系具有自我保護的生物學功能，即細胞質母性遺傳的特點，不能產生花粉，無法自交留種，所以該紫色資源也無法應用。目前，紫寶由國內代理公司銷售，品種名稱也不同，且種子價格昂貴，每粒種子平均售價為0.8～1.0元。2014年Ahmed等[9]報道了全紫色大白菜材料Palgangbaechu和部分紫色材料8053。2016年，韓國又推出了大白菜紫色品種紫裔。當前，有關紫寶、紫裔、Palgangbaechu和8053等紫色性狀的來源仍不清楚，猜測有可能來源于紫色芥菜。

轉錄組是基因組遺傳信息與生物功能蛋白質組相連接的必然紐帶，也是研究基因功能及結構的基礎，已廣泛應用于基礎研究等多個領域。轉錄組測序的研究對象為特定細胞在某一功能狀態下所能轉錄出來的所有RNA的總和，目前主要是mRNA。通過新一代高通量測序技術，能夠全面快速地獲得紫色大白菜葉片在特定狀態下的幾乎全部轉錄本序列信息，特別是有關類黃酮、黃酮、黃酮醇和花青素合成途徑相關基因的轉錄本信息，有利于獲得差異表達的基因，闡明調控葉片花青苷合成相關基因的分子機制[10-11]。

2015年春季，筆者在1份繁育的小白菜雜交一代中發現了1株紫色材料，不僅紫色性狀突出，而且葉片正反面、葉脈和葉柄均表現為紫色，命名為15NG28。隨后，將15NG28與大白菜雜交、回交，繼續轉育該紫色性狀，已初步創制了一批大白菜新型紫色資源材料。本研究擬開展15NG28及其雜交F1、回交一代BC1F1和自交等后代的植物學性狀、基因組的倍性和染色體數、轉錄組學等研究，為拓寬新的紫色大白菜資源，最終選育出具有我國自主知識產權的深紫色白菜新品種并進一步獲得控制15NG28紫色葉片花青苷合成途徑候選基因奠定基礎。

1 材料和方法

1.1 試驗材料

2015年4月，選用2份大白菜(1份春白菜15綠1-1-26，1份秋白菜16-1081)，分別與15NG28雜交，配制3組F1：15NG28×16-1081、15NG28×15綠1-1-26和15綠1-1-26×15NG28，同時在蕾期對15NG28進行人工剝蕾、自交授粉，收獲全部種子。12月18日，在溫室中播種15NG28自交種子、3組F1材料。2016年1月13日開始，觀察全部材料葉片在幼苗期和成株期的紫色性狀，并拍照。4月，配制回交一代BC1F1：16M-170-22×15綠1-1-26，收獲全部種子。12月20日，播種2組BC1F1種子，2017年4月，觀察2組BC1F1葉片紫色性狀，并拍照，同時在蕾期對紫色BC1F1進行人工剝蕾、自交授粉，紫色BC1F1繼續與白菜進行回交轉育。試驗材料均來源于北京市農林科學院蔬菜研究中心大白菜課題。所用試驗材料的世代、名稱和特性見表1。

表1 材料的世代、名稱和特性Tab.1 Generation，name and characteristics of materials

1.2 觀察葉片中花青素的分布情況

采用徒手切片的方法，將2把刀片并排重疊，快速切取紫色和綠色葉片，放于清水中，并用毛筆撈取，置于載玻片上，在顯微鏡下觀察葉片中花青素的分布情況。

1.3 觀察基因組倍性

參考饒琳莉[12]的方法并進行優化，觀察樣品的基因組倍性，具體操作為：選取樣品新鮮嫩葉100 mg，加入1 mL預冷的緩沖液Galbraith(45 mmol/L MgCl2，30 mmol/L檸檬酸納，20 mmol/L MOPS，0.1% Triton X，pH=7)于培養皿中，用鋒利的刀片快速切碎嫩葉，再用50 μm孔徑的尼龍膜過濾，收集濾液，1 000 r/min離心5 min，倒掉上清液，加入200 μL PI染液，振蕩混勻，避光靜置20 min后，用流式細胞儀進行檢測分析并拍照。試驗在北京市農林科學院蔬菜研究中心生物技術研究室完成。

1.4 染色體數觀察

在晴天上午(9：00-11：00)，選取植株盛花期花序上端的幼小花蕾，用卡諾氏液(V(無水酒精)∶V(冰乙酸)=3∶1)固定24 h，再換至70%乙醇中，置于4 ℃冰箱保存備用。以改良品紅溶液染色，常規壓片，在顯微鏡下進行觀察，統計染色體數目并拍照。試驗完成地點同1.3。

1.5 轉錄組RNA-Seq分析

選取F1紫色植株16M-170-22、BC1F1紫色植株17M-245-21和17M-247-22、F1綠色植株16M-170-23的葉片，3次生物學重復，送交杭州壹基因公司進行轉錄組分析，并獲得測序結果報告。

獲得差異表達基因的方法：以紫色單株BCZ為參考，對綠色單株BCL與紫色單株BCZ基因表達量FPKM(Fragments Per Kilobase of transcript per Million fragments mapped)進行比較，即BCL-VS-BCZ，再依據錯誤發現率FDR(False Discovery Rate)篩除假陽性，最終獲得BCZ、BCL 2個樣品之間的差異表達基因。

2 結果與分析

2.1 15NG28植物學性狀

觀察15NG28葉片顏色性狀，結果表明，在幼苗期，15NG28的子葉和真葉正反面，葉柄、葉脈均為紫色；在成株期，15NG28葉片為長圓形，葉緣淺鋸齒，葉柄窄，葉柄正面為紫色、背面為綠色，葉柄沒有葉翼(圖1)。由于在繁育的普通小白菜雜交種中發現15NG28材料，所以可以肯定15NG28是以小白菜為母本、父本未知的F1雜交種，因此，15NG28表現為中間類型，偏向小白菜，但不表現白菜的抱球、小白菜的束腰等性狀(圖1)，也沒有典型剖面特征；在花期，15NG28花瓣顏色為黃色，花蕾小，花瓣、花藥和柱頭等花器官也小，花粉偏少，但腋芽萌生能力特別旺盛，花枝細長且軟。15NG28自交結實結果表明：15NG28自交結實率極低。

A.全株；B.葉柄。A.Whole plant；B.Leaf blade.

2.2 15NG28與白菜雜交的F1性狀

以15NG28為母本，分別與大白菜15綠1-1-26和16-1081雜交，均可以獲得少部分F1種子，而以15NG28為父本與15綠1-1-26雜交，后代結實率很高。比較15NG28及其3組F1植株葉片性狀，結果表明：以15NG28為母本的2組雜交F1(15NG28×16-1081、15NG28×15綠1-1-26)，其葉片顏色出現了分離，表現為紫色和綠色2類(表2)。其中紫色F1的子葉、真葉正反面均為紫色，特別是葉脈和葉柄也為紫色，且顏色深(圖2-4)。另外，調查了雜交F1(15綠1-1-26×15NG28)的性狀，結果表明：9株F1植株全部表現為綠色，而且植物學性狀相似，因此，這9株F1應該為自交種。同時，調查了紫色F1植株16M-170-21現蕾期、抽薹期和開花期的性狀，結果表明：隨著植株逐漸經歷現蕾期、抽薹期和開花期的過程，葉片紫色逐漸褪色、變淺和變綠，在花期，只有葉脈部分呈現出零星紫色，因此，紫色F1植株葉片紫色由深到淺的顯色順序為：苗期>現蕾期>抽薹期>開花期(圖2，5)。可見，利用常規雜交的方法，可以獲得少部分15NG28與大白菜的雜交F1種子，且F1后代分離出部分紫色植株。

表2 15NG28及其雜交F1的性狀Tab.2 Characteristics of 15NG28 and it′s F1 hybrids

A.16M-170-21；B.16M-170-22；C.16M-171-24；D.16M-171-26。

圖3 16M-170-22苗期真葉反面性狀Fig.3 Abaxial leaf characteristics of 16M-170-22 at seedling stage

A.16M-171-24；B.16M-170-21；C.16M-170-22；D.16M-171-26。

圖5 16M-170-21在現蕾期、抽薹期和開花期的性狀Fig.5 Characteristics of 16M-170-21 at budding，bolting and flowering stage

2.3 紫色F1的自交、回交一代BC1F1植物學性狀

2017年4月，重點調查了F116M-170-21的自交后代以及回交一代BC1F1(16M-170-22×15綠1-1-26)的性狀。結果表明：16M-170-21自交后代表現出多種類型，其中8株的主要性狀為：芥菜型，葉片顏色深紫色，圓形葉，葉緣鋸齒，深淺不等，葉柄窄、細長，無葉翼；另外7株的主要性狀為：白菜型，葉片顏色深紫色，圓形葉，葉柄長，有葉翼(圖6)。31株BC1F1性狀全部表現為白菜類型，其中14株葉片為紫色，17株葉片為綠色。10株BC1F1后代性狀見圖7。2017年5月8日，調查BC1F1紫色植株自交結實情況，結果表明：14株紫色植株結實性很好。

圖6 16M-170-21自交后代性狀Fig.6 Characteristics of self-inbred plants of 16M-170-21

2.4 葉片中花青素的分布情況

在顯微鏡下觀察葉片花青素的分布情況，結果表明：F1紫色大白菜16M-170-22的花青素主要分布在葉片的上表皮和下表皮，下表皮少于上表皮，而在綠色大白菜16-1081的上、下表皮中均沒有觀察到花青素(圖8)。

A.紫色株；B.綠色株。A.Purple plants；B.Green plants.

A.16-1081；B.16M-170-22。

2.5 植株的倍性分析

2.5.1 紫色F1的倍性 經流式細胞儀鑒定，二倍體白菜2n熒光強度高峰位于48(圖9-A)，異源四倍體芥菜4n熒光強度高峰位于96(圖9-B)。在F1紫色植株中，16M-170-21熒光強度高峰位于105，倍性為4n(圖10-A)，16M-170-22熒光強度高峰位于80(圖10-B)，推測為3n，即2n+n，16M-171-24熒光強度高峰位于150，倍性為6n(圖10-C)，16M-171-26熒光強度高峰位于70(圖10-D)，推測為3n，即2n+n。

2.5.2 綠色F1的倍性 與紫色F1植株一樣，綠色F1植株也存在復雜的倍性變化。圖11顯示為四倍體和混倍體。

2.6 染色體數分析

經觀察發現，在F1紫色大白菜中，16M-170-21染色體數為36條(圖12)，基因型可能為芥菜型AABB，即IIA+IIB；16M-170-22染色體數為28條(圖13)，基因型可能為AAB，即IIA+IB；16M-171-24染色體數為24條(圖14)，基因型可能為ABB，即IIB+IA；16M-171-26染色體為28條(圖15)，基因型可能為AAB，即IIA+IB。

A.白菜二倍體；B.芥菜異源四倍體。A.B.rapa diploid；B.B.juncea heterotetraploid.

A.16M-170-21；B.16M-170-22；C.16M-171-24；D.16M-171-26。

A.四倍體F1植株；B.混倍體(二倍與四倍F1植株)。A.4n；B.Mix-ploid (2n and 4n).

圖12 16M-170-21染色體(4n=36，100×)Fig.12 Chromosome of 16M-170-21(4n=36，100×)

圖13 16M-170-22染色體(3n=28，100×)Fig.13 Chromosome of 16M-170-22(3n=28，100×)

圖14 16M-171-24染色體(3n=24，100×)Fig.14 Chromosome of 16M-171-24(3n=24，100×)

圖15 16M-171-26染色體(3n=28，100×)Fig.15 Chromosome of 16M-171-26(3n=28，100×)

2.7 轉錄組測序分析

F1紫色植株16M-170-22、BC1F1紫色植株17M-245-21和17M-247-22、F1綠色植株16M-170-23的植物學性狀特征見圖16。利用Illumina HiSeqTM 4000，對紫色單株BCZ(16M-170-22、17M-245-21、17M-247-22)和綠色單株BCL(16M-170-23)進行初步的轉錄組測序分析。將原始序列數據(Raw data)進行過濾處理，去除接頭污染Reads、低質量Reads和含N比例大于5%的Reads，最后得到Clean data。

使用比對軟件Tophat(http：//ccb.jhu.edu/software/tophat/index.shtml)，將每個樣品的Clean data與白菜參考基因組(http：//brassicadb.org/brad/)進行比對。由表3可見，紫色單株BCZ的定位總個數百分比、雙末端定位個數百分比、左端定位個數百分比、右端定位個數百分比均低于綠色單株BCL，體現了紫色單株BCZ基因組的復雜性，進一步證明15NG28的紫色性狀可能并非來源于大白菜的事實。

以紫色單株BCZ為參考，經進一步比較紫色單株16M-170-22、17M-245-21、17M-247-22和綠色單株16M-170-23差異表達基因的表達量FPKM，并篩除假陽性，結果表明：明顯上調的基因共有140個，明顯下調的基因共有16個。在差異表達基因分析結果中，選取參與黃酮類合成代謝通路的相關基因，如在紫色單株明顯上調表達的15個差異表達基因，進行聚類分析(圖17)。

A.16M-170-22；B.16M-170-23；C.17M-245-21；D.17M-247-22。

樣品Samples讀碼總數Totalreads定位總個數百分比/%Percentoftotalmappedreads雙末端定位個數百分比/%Percentofmappedpairreads左端定位個數(百分比)Percentofmappedleftreads右端定位個數(百分比)PercentofmappedrightreadsBCL4787061484.975.920554669(85.9%)20066163(83.8%)BCZ4659573873.564.417337028(74.4%)16929714(72.7%)

每行代表一個基因，不同顏色代表基因的表達變化倍數或表達量。對于差異表達基因聚類：淺顏色表示表達量上調，深顏色表示表達量下調。對于樣品聚類：黑色到白色表示表達量從小到大，顏色越淺表示差異表達基因的表達量越高。

Each row presented one gene，and different colors meant expression amount or expression level.For differentially expressed genes in each row，light color meant up-regulation and dark color meant down-regulation.For clustering chart in each column，colors from dark to light meant low to high gene expression level.

圖1715個差異表達基因聚類模式圖
Fig.17Clusteringchartoffifteendifferentiallyexpressedgenes

進一步對紫色單株和綠色單株差異表達基因的GO功能注釋進行分類作圖，結果表明：3種GO功能注釋分別描述為基因的分子功能(Molecular function)、細胞成分(Cellular component)和參與的生物過程(Biological process)(圖18)。圖18解釋了在差異基因背景和全部基因背景下，GO各二級功能的18個生物過程、10個細胞成分和7個分子功能差異基因富集情況，體現了具有明顯比例差異的二級功能差異基因富集與全部基因富集的不同趨勢。

同時，對紫色單株BCZ中參與花青素合成相關代謝通路及其基因進行分析，發現參與類黃酮、黃酮和黃酮醇、花青素生物合成等3個代謝通路相關的差異表達基因共12個：7個類黃酮生物合成途徑基因生物，4個黃酮和黃酮醇生物合成途徑基因，1個花青素合成途徑基因(表4)。另外，12個差異表達基因在紫色材料BCZ中差異表達量最高的2個基因均為結構基因：Bra013652(AT4G22880)為無色花青素雙加氧酶(LDOX)，Bra027457(AT5G42800)為二氫黃酮醇還原酶(DFR)。12個差異表達基因在白菜和擬南芥中比較和功能見表5，將是今后的重點研究內容。

3 討論與結論

有關大白菜紫色資源的創制以及與花青苷合成相關的基因表達、克隆等內容是目前國內外研究的熱點，如紫色大白菜[13-19]、紫色小白菜[20-22]、紫色菜薹[22]等。

目前，國內普遍利用紫色小白菜與大白菜雜交，獲得葉片紫色的大白菜材料，育成紫色白菜品種，如京研紫快菜等，但是該紫色性狀僅在葉片正面表現出深紫色，而子葉、真葉和心葉的背面均為綠色，尤其是葉球剖面呈現出紫色、白色、黃色等多種顏色，即紫色表現不均勻，因此，該材料在大白菜育種中存在一定的局限性，該紫色品種的推廣也受到一定限制[3-4]。西北農林科技大學選育出高花青素含量大白菜純系11S91，該紫色性狀來源于紫色菜薹，但其外葉、葉柄和葉脈均為綠色，只有葉球剖面的心葉部分表現為紫色，花青素主要分布在葉片表皮下2～3層細胞[15，19]。紫色性狀表現優異的大白菜資源主要有2種：紫紅色大白菜新種質B90335-5[2，10-11]和紫色大白菜品種，如紫寶、紫裔。研究表明，紫紅色大白菜B90335-5紫紅色性狀不符合孟德爾遺傳規律[10]。當前國內還沒有育成類似于紫寶的紫色大白菜品種，因此，大白菜紫色資源的原創工作非常迫切和重要。喬海云等[23]獲得了大白菜與紫色甘藍的種間雜種。15NG28與大白菜雜交F1、回交BC1F1的紫色植株性狀不僅葉片正反面、葉脈全部為紫色，而且葉柄正反面也為紫色，與紫紅色芥菜、紫寶的紫色性狀較一致，與來源于紫色小白菜、紫色菜薹的紫色性狀明顯不同。初步認為，新型紫色材料15NG28的紫色性狀可能來源于紫色芥菜。同時，僅用常規方法，15NG28可以成功與白菜進行雜交、回交，轉育紫色性狀，且后代較容易獲得種子。雖然雜交F1、BC1F1后代的倍性、染色體行為存在復雜性，這也說明芥菜的紫色性狀整合到白菜基因組是一個復雜的、漫長的生物學過程，但是可喜的是，本研究僅在31株BC1F1(16M-170-22×15綠1-1-26)后代中就發現了14個紫色植株，而且全部植株均為白菜，結實性好，植物學性狀、染色體數目等都具有白菜基因組AA的特征特性。下一步，筆者將利用2株BC1F1(17M-245-21、17M-247-22)分別與綠色大白菜骨干親本系雜交，繼續轉育紫色性狀，創制大白菜骨干親本系的紫色資源，為選育出紫色大白菜新品種奠定基礎。因此，15NG28、雜交F1，尤其是BC1F1后代中的紫色植株具有潛在的應用價值。

橫坐標為GO 各分類內容，縱坐標左側為基因數目所占百分比，右側為基因數目。矩形直方圖代表上調基因，橢圓圖代表下調基因。

通路Pathway注釋通路的DEG基因DEGswithannotationpathway(107)注釋通路的全部注釋基因Allgeneswithannotationpathway(22358)P值PvalueQ值Qvalue通路路徑PathwayID差異表達的基因Differentiallyexpressedgenes類黃酮生物合成Flavonoidbiosynthesis7(6.54%)285(1.27%)0.00044723050.013864146ko00941Bra027457、Bra030550、Bra013652、Bra036208、Bra019350、Bra007142、Bra009312黃酮和黃酮醇生物合成Flavoneandflavonolbiosynthesis4(3.74%)127(0.57%)0.0032644060.050598293ko00944Bra034677、Bra009312、Bra008316、Bra035271花青素生物合成Anthocyaninbiosynthesis1(0.93%)13(0.06%)0.060475140.288419898ko00942Bra038445

表5 白菜差異表達基因在擬南芥中名稱和功能注釋Tab.5 Names and annotations of differentially expressed genes in B.rapa and A.thaliana

蕓薹屬禹氏三角(U′s triangle)包括3個二倍體基本種：白菜AA(B.rapa，2n = 20)、甘藍CC(B.oleracea，2n = 18)、黑芥BB(B.nigra，2n = 16)和3個四倍體復合種：甘藍型油菜AACC(B.napus，2n=38)、芥菜型油菜AABB(B.juncea，2n=36)和埃塞俄比亞芥BBCC(B.carinata，2n=34)。測序結果表明，白菜基因組大小為485 Mb，黑芥基因組大小為591 Mb，甘藍基因組大小為630 Mb，而芥菜基因組大小為922 Mb。本研究中，以觀察15NG28、雜交F1、回交BC1F1后代的倍性和染色體數目為輔助手段，重點選擇具有更多AA基因組特征的后代進行人工自交、回交授粉，以減少工作量，進而快速獲得具有大白菜性狀的紫色中間材料，如17M-245-21、17M-247-22等。在試驗中，流式細胞儀峰值的位置一方面受植物染色體組成決定，另一方面還有2.5%～5.0%的系統誤差。通過流式細胞儀，很難確定某個材料的基因組組成。下一步，將利用黑芥的特異標記等，輔助篩選具有芥菜紫色性狀的大白菜中間材料。同時，我們還發現，16M-171-24染色體數為24條，可能為三倍體，但流式細胞儀檢測的結果為六倍體，二者結果不同，可能原因是16M-171-24同株上同時出現了三倍體和六倍體的枝條，而筆者選取了三倍體花蕾觀察染色體數目。在白菜和芥藍等蕓薹屬小孢子培養過程中，經常可以發現在同一DH植株上生長不同倍性枝條的情況。同樣，即使經過多次人工自交和與白菜回交，仍沒能獲得16M-171-24自交種子以及回交后代種子，說明16M-171-24具有復雜的染色體和基因組。

目前，有關大白菜花青苷結構基因和轉錄因子等有關研究已見報道。Xie等[11]對紫紅色和綠色大白菜進行轉錄組學研究，共獲得可定位在白菜參考基因組上的差異表達基因(DEGs)930個，其中389個基因上調表達，541個基因下調表達。通過進一步重新組裝未定位的reads，又發現了2 031個Unigenes。作者認為芥菜的R2R3-MYB轉錄因子c3563g1i2被導入到大白菜基因組中，可能參與大白菜花青素合成。本研究對F1的1個綠色植株16M-170-23、1個紫色植株16M-170-22和BC1F1的2個紫色植株17M-245-21、17M-247-22分別進行轉錄組測序，結果表明：在紫色材料中差異表達量最高的2個基因均為結構基因：Bra013652(LDOX)，Bra027457(DFR)，該結論與Xie等[11]報道控制紫紅色大白菜的R2R3-MYB轉錄因子c3563g1i2完全不同，但與Xie在比對花青素合成途徑的轉錄子后發現花青素合成途徑中Bra013652和Bra027457基因的上調表達量最高及這2個基因的Log2FC值10.05和10.07也最高的結論相吻合。由此可見，15NG28及其雜交、回交后代可能為新的大白菜紫色資源材料。

2017年秋季，將進一步調查14株BC1F1紫色白菜的自交后代和回交二代BC2F1的植物學性狀和葉片顏色分離情況，重點研究17M-245-21和17M-247-22的自交后代和BC2F1后代的性狀，以期獲得一批紫色大白菜中間材料，再通過人工加代等多代自交和游離小孢子技術獲得紫色大白菜的純合系和DH系，最后配制紫色與綠色的F1，構建葉片顏色分離群體，如F2和DH群體，進一步開展紫色性狀的遺傳規律、染色體定位等相關研究。同時，將利用實時定量熒光PCR(qRT-PCR)技術，對12個差異表達基因進行分析，進一步研究控制15NG28紫色性狀花青素合成途徑的相關基因。

[1] 張德雙，張鳳蘭，于拴倉，等.與紫色大白菜花青苷合成相關的基因[J].分子植物育種，2015，13(11)：2619-2625.

[2] 孫日飛，張淑江，章時蕃，等.紫紅色大白菜種質的創新研究[J].園藝學報，2006，33(5)：1032-1032.

[3] 張德雙，張鳳蘭，余陽俊，等.紫色大白菜育種材料的創造[J].長江蔬菜，2007(11)：52-53.

[4] 張德雙，張鳳蘭，余陽俊，等.紫色大白菜育種思路初探[J].長江蔬菜，2008(22)：14-17.

[5] 劉 瑾，汪維紅，張德雙，等.控制白菜葉片紫色的pur基因初步定位[J].華北農學報，2013，28(1)：49-53.

[6] Wang W H，Zhang D S，Yu S C，et al.Mapping theBrPurgene for purple leaf color on linkage group A03 ofBrassicarapa[J].Euphytica，2014，199(3)：293-302.

[7] 張明科，張魯剛，鞏振輝，等.白菜紫色性狀RAPD連鎖標記的篩選與染色體定位研究[J].西北植物學報，2008，28(5)：901-906.

[8] Hayashi K，Matsumoto S，Tsukazaki H，et al.Mapping of a locus regulating anthocyanin pigmentation inBrassicarapa[J].Breeding Science，2010，60(1)：76-80.

[9] Ahmed N U，Park J I，Jung H J，et al.Anthocyanin biosynthesis for cold and freezing stress tolerance and desirable color inBrassicarapa[J].Functional & Integrative Genomics，2015，15(4)：383-394.

[10] 張淑江．大白菜紫紅色基因定位和候選基因分析[D].北京：中國農業科學院研究生院，2014.

[11] Xie L L，Li F，Zhang S F，et al.Mining for candidate genes in an introgression line by using RNA sequencing：the anthocyanin over accumulation phenotype inBrassica[J].Frontiers in Plant Science，2016，7：1245.

[12] 饒琳莉.榨菜與紫甘藍種間異源六倍體的植物與分子特性初探[D].杭州：浙江大學，2015.

[13] Burdzinski C，Wendell D.Mapping the anthocyaninless (anl) locus in rapid-cyclingBrassicarapa(RBr)to linkage group R9[J].BMC Genetics，2007，8：64.

[14] Zhang J，Lu Y，Yuan Y，et al.Map-based cloning and characterization of a gene controlling hairiness and seed coat color traits inBrassicarapa[J].Plant Molecular Biology，2009，69(5)：553-563.

[15] 段巖嬌，張魯剛，何 瓊，等.紫心大白菜花青素積累特性及相關基因表達分析[J].園藝學報，2012，39(11)：2159-2167.

[16] Guo N，Cheng F，Wu J，et al.Anthocyanin biosynthetic genes inBrassicarapa[J].BMC Genomics，2014，15(1)：426.

[17] Kim C K，Kim J A，Kikuchi S，et al.Computational identification of Chinese cabbage anthocyanin-specific genes[J].BioChip Journal，2011，5(2)：184-192.

[18] Wang Z，Tang J，Hu R，et al.Genome-wide analysis of the R2R3-MYB transcription factor genes in Chinese cabbage (Brassicarapassp.pekinensis) reveals their stress and hormone responsive patterns[J].BMC Genomics，2015，16：17.

[19] He Q，Zhang Z F，Zhang L G.Anthocyanin accumulation，antioxidant ability and stability，and a transcriptional analysis of anthocyanin biosynthesis in purple heading Chinese cabbage (BrassicarapaL.ssp.pekinensis)[J].Journal of Agricultural and Food Chemistry，2016，64(1)：132-145.

[20] 張 彬.蕓薹屬植物花青素生物合成代謝途徑調控機制的研究[D].重慶：重慶大學，2011.

[21] Zhang Y J，Chen G P，Dong T T，et al.Anthocyanin accumulation and transcriptional regulation of anthocyanin biosynthesis in purple Bok-Choy (Brassicarapavar.chinensis)[J].Journal of Agricultural and Food Chemistry，2014，62(51)：12366-12376.

[22] 郭 寧.白菜花青苷和黃酮醇苷自然變異及遺傳機制研究[D].北京：中國農業科學院，2014.

[23] 喬海云，李 菲，張淑江，等.大白菜與紫甘藍種間雜種的獲得與鑒定[J].植物科學學報，2012，30(4)：407-414.