孫 葉 劉 紅 李風童 劉春貴 周春華 包建忠 陳秀蘭 趙國琦,*

(1 揚州大學動物科學與技術學院，江蘇 揚州 225009；2 江蘇里下河地區(qū)農業(yè)科學研究所，江蘇 揚州 225007)

春蘭(Cymbidiumgoeringii)是蘭科(Orchidaceae)蘭屬(Cymbidium)的原生種之一，大花蕙蘭(C.hybridum)是蘭屬的雜交種。蘭屬植物約有48個原生種，目前已知用來作為大花蕙蘭雜交親本的原生種有近20種，其中也包括春蘭。春蘭是傳統(tǒng)國蘭的重要代表種之一，自然花期正值我國傳統(tǒng)節(jié)日春節(jié)前后，在年銷花市場有很大的開發(fā)潛力。受傳統(tǒng)蘭文化的影響，春蘭優(yōu)良品種以瓣型綠花為主，紅花品種稀缺。大花蕙蘭是當今全球暢銷的主要蘭花種類之一，作為我國年銷花市場高檔盆栽花卉的明星品種，年銷售量在200萬盆以上[1]。大花蕙蘭品種的花形、花色極為豐富，在花色素、花色基因及轉基因調控花色等研究方面也取得了很大的進展[2-3]，但大多數品種無香味、葉形散亂，在我國達不到春化溫度的南方地區(qū)不易種植開花[4]。目前利用春蘭優(yōu)良品種與大花蕙蘭遠緣雜交選育小花型、花色豐富、有香味、耐寒的蘭屬新品種是市場開發(fā)的熱點[5]，同時這些新品種也是研究花色形成和呈色機理的理想材料。

高通量轉錄組測序技術彌補了非模式生物缺乏參考基因組信息的不足，已在觀賞植物領域廣泛應用。近年來春蘭、蕙蘭、建蘭、墨蘭、春劍等國蘭利用轉錄組測序技術，在花發(fā)育、花色、葉色、葉形等相關功能基因的挖掘和基因組轉錄調控分析方面取得了很大的進展[6-12]。蘭科植物的花色主要由類黃酮、類胡蘿卜素和葉綠素的含量及其比例決定，花青苷是紅花顯色最重要的色素。植物中類黃酮、類胡蘿卜素和葉綠素合成代謝途徑中的關鍵結構基因已研究得比較清晰，影響結構基因表達調控的轉錄因子正成為研究熱點，其中MYB(V-myb avain myeloblastosis viral ohcogene)轉錄因子的研究最為廣泛[13-15]。目前已在擬南芥中發(fā)現了339個MYB轉錄因子，位于美國斯坦福大學的擬南芥信息資源網站(TheArabidopsisInformation Resource, TAIR)已公布了133個MYB成員。MYB轉錄因子中 R2R3-MYB 的數量最多，在擬南芥和水稻已分別鑒定出126和109個R2R3-MYB家族成員[16]。根據蛋白C端保守的氨基酸基序的分析，R2R3-MYB轉錄因子劃分為 28 個不同的亞家族[17]，其中S7、S6 和 S2亞家族分別對黃酮醇[18]、花青素及原花青素[19]的調控起關鍵作用。植物花青苷的合成代謝途徑中MBW(MYB-bHLH-WD40)復合體調控模式處于核心地位[20]，其中R2R3-MYB轉錄因子起關鍵性作用[21]，Hsu等[22]認為PeMYB2、PeMYB11和PeMYB12三個 R2R3-MYB轉錄因子以不同比例的組合表達導致蝴蝶蘭中大量復雜的花色素沉著模式。

蘭科植物調控花色的結構基因和調節(jié)基因的研究已取得了很大的進展[23]，但相關基因的克隆和功能驗證主要集中在洋蘭上，關于國蘭的研究進展緩慢，春蘭花色形成的分子機理更鮮見報道。前期研究[24]發(fā)現綠花春蘭品種間雜交，F1出現了紅花新種質；綠花春蘭品種與黃花大花蕙蘭品種雜交，F1同樣也出現了紅花新種質，且花器官不同發(fā)育階段顯色差異較大。本研究采用綠花春蘭品種與黃花大花蕙蘭品種的紅花雜交后代黃金梅2個不同發(fā)育期的花瓣進行轉錄組測序，解析測序數據，系統(tǒng)挖掘春蘭與大花蕙蘭種間雜交種花色相關的關鍵結構基因和轉錄因子，探究基因的表達調控及不同發(fā)育期花瓣色素含量的變化，旨在為花色基因的克隆及功能分析，花色調控及呈色機理的解析和新種質的定向選育提供重要的參考依據。

1 材料與方法

1.1 試驗材料及轉錄組測序

供試材料黃金梅是母本綠花春蘭大宋梅(圖1-A)與父本黃花大花蕙蘭黃金虎(圖1-B)的雜交種，保存于江蘇里下河地區(qū)農業(yè)科學研究所中國蘭種質資源圃。該種質花器官遺傳了母本的瓣型和父本的著花數特征，花色發(fā)生了變異，小花蕾期為淡黃綠色，隨著花蕾的發(fā)育逐漸轉色，綻蕾后花色漸變?yōu)樽霞t色，與父、母有明顯差異。2019年1月選取黃金梅同一株系不同發(fā)育期小花蕾樣品(圖1-C)、黃綠色小花蕾花瓣(主瓣和側瓣)(圖1-D)、紫紅色始花期小花瓣(主瓣和側瓣)(圖1-E)，用于轉錄組測序，剩余花瓣同時采樣用于花色素測定。取樣前用無水乙醇輕擦拭外瓣，分別設3個生物學重復，取樣后立即用液氮凍存，測序材料用干冰保存送至上海為青生物技術有限公司進行RNA提取，構建cDNA文庫，在Illumina HiseqTM 2000測序平臺進行高通量測序；得到原始測序序列(sequenced reads)，即原始數據(raw data)。通過過濾去除N(N表示無法確定堿基信息)的比例大于10%、低質量序列、帶接頭(adapter)的reads，最終獲得Clean Data后，利用Trinity軟件對Clean Data進行序列組裝，得到Unigene庫。

注：A:母本大宋梅；B:父本黃金虎；C:雜交種黃金梅；D:小花蕾期；E：始花期。Note: A:Female parent Dasongmei. B:Male parent Huangjinhu.C:Hybrid Huangjinmei. D:Little flower Bud stage. E:Opening stage.圖1 取樣材料及其親本Fig.1 Sampling materials and their parents

1.2 差異表達基因功能注釋、富集分析及花色素基因分析

使用BLAST軟件將測序Unigene與Pfam、TrEMBL、eggNOG、Swiss-Prot、GO、KEGG數據庫進行序列比對，獲得Unigene的注釋信息，對差異表達基因(differential express genes, DEGs)做GO功能富集分析、KEGG注釋及通路富集分析。檢索類黃酮、類胡蘿卜素及葉綠色合成途徑中關鍵結構基因及花青素合成相關轉錄因子信息，獲取相關基因的核苷酸長度及2樣本差異表達的FPKM(fragments per kilobase per million)值。利用BLAST軟件在NCBI數)據庫中篩選花青苷合成途徑中下游關鍵基因和轉錄因子的同源基因。

1.3 基因轉錄模式驗證

選擇類黃酮及類胡蘿卜素形成關鍵通路14個表達差異較大基因設計引物(表1)，以甘油醛-3-磷酸脫氫酶(glyceraldehyde-3-phosphate dehydrogenase, GAPDH)為內參，利用轉錄組測序時提取的小花蕾、始花期花瓣RNA反轉錄獲得cDNA作為熒光定量PCR(quantitative real-time PCR, qPCR) 反應的模板，在LightCycler96實時熒光定量PCR儀(瑞士羅氏公司)中進行，驗證基因的表達。根據SYBR?PremixExTaqⅡ (Tli RNaseH Plus)(日本TaKaRa公司)說明書配置反應體系，每個樣品設3次重復。反應體系：SYBRPremixExTaqⅡ (Tli RNaseH Plus) (2×) 10.0 μL、10 μmol·L-1正反向引物各0.8 μL、cDNA模板50～100 ng、ddH2O補足至20 μL。PCR擴增程序：95℃預變性30 s；95℃變性5 s，60℃退火，延伸20 s， 40個循環(huán)；熔解曲線分析程序：95℃持續(xù)1 s，65℃持續(xù)60 s，97℃持續(xù)1 s。用2-ΔΔCt方法進行相對定量分析。

表1 基因及內參引物信息Table 1 The primer information of gene and internal reference

1.4 總花青素苷、總黃酮、類胡蘿卜素及葉綠素含量測定

花青素及黃酮提取參考Yang等[25]的方法，總花青素苷含量(total anthocyanin content, TAC)測定采用可見分光光度法[26]，總黃酮含量(total flavonoids content, TFC)的測定采用AlCl3顯色法[27],重復3次。參照 Lichtenthaler 等[28]的方法進行葉綠素和類胡蘿卜素含量的測定，3次重復。采用SPSS 17.0 進行t檢驗計算不同發(fā)育期色素含量的差異顯著性。

2 結果與分析

2.1 轉錄組測序數據分析與組裝

提取試驗材料黃金梅小花蕾期和始花期花瓣的總RNA，結果顯示，A樣品RNA濃度為276 μg·μL-1，總量為27.6 μg；B樣品RNA濃度為291 μg·μL-1，總量為29.1 μg；OD260/280>2.0，OD260/230>1.8，RIN>9.0，28S:18S>1.6；凝膠電泳結果顯示5S、18S、28S條帶完整，樣品質量符合建庫標準要求。建庫后對2個樣品進行轉錄組測序，測序數據經FASTX過濾后共獲得高質量堿基序列140 453 340個，堿基數為39 908 521 302 bp， 2個樣品Mapped Reads在Clean Reads中所占比例分別為86.83%(A)、88.29%(B)；質量值Q20和Q30均大于90%，GC含量均大于47.5%，證明測序質量好。Trinity合并組裝后分析統(tǒng)計拼接長度數據(表2)，得到 188 837條Transcript和113 780條Unigene，其中200～300 bp的Unigene比例最高有48 564 條，比例最少的為1 500～2 000 bp的Unigene 4 041條，僅占3.55%；組裝得到Transcript與Unigene的平均長度為725.64 bp和586.53 bp，N50分別為1 229 bp 和882 bp，組裝完整性較高，可進行下一步分析。

表2 轉錄組長度分布統(tǒng)計表Table 2 Statistical table of the length distribution of the transcriptome

注：1：胞外區(qū)；2：細胞；3：細胞核；4：細胞膜；5：細胞連接；6：胞外基質；7：膜封閉腔；8：復雜大分子；9：細胞器；10：細胞外基質部分；11：細胞外區(qū)域部分；12：細胞器組分；13：病毒粒子部分；14：膜結構；15：細胞組分；16：蛋白結合轉錄因子活性；17：核酸結合轉錄因子活性；18：催化活性；19：受體活性；20：鳥嘌呤核苷酸交換因子活性；21：結構分子活性；22：轉運活性；23：結合；24：電子載體活性；25：抗氧化活性；26：金屬伴隨子活性；27：酶調節(jié)活性；28：營養(yǎng)貯存活性；29：分子轉導活性；30：繁殖；31：細胞殺傷；32：免疫系統(tǒng)過程；33：代謝過程；34：細胞過程；35：生殖過程；36：信號傳導；37：多細胞生物過程；38：發(fā)育過程；39：生長；40：運動；41：單一生物過程；42：生物調節(jié)；43：刺激反 應；44：定位；45：多細胞生物過程；46：生物調節(jié)；47：細胞組成機制或生物合成。Tips: 1: Extracellular region. 2: Cell. 3: Nucleoid. 4: Membrane. 5: Cell junction. 6: Extracellular matrix. 7: Membrane-enclosed lumen. 8: Macromolecular complex. 9: Organelle. 10: Extracellular matrix part. 11: Extracellular region part. 12: Organelle part. 13: Virion part. 14: Membrane part. 15: Cell part. 16: Protein binding transcription factor activity. 17: Nucleic acid binding transcription factor activity. 18: Catalytic activity. 19: Receptor activity. 20: Guanyl-nucleotide exchange factor activity. 21: Structural molecule activity. 22: Transporter activity. 23: Binding. 24: Electron carrier activity. 25: Antioxidant activity. 26: Metallochaperone activity. 27: Enzyme regulator activity. 28: Nutrient reservoir activity. 29: Molecular transducer activity. 30: Reproduction. 31: Cell killing. 32: Immune system process. 33: Metabolic process. 34: Cellular process. 35: Reproductive process. 36: Signaling. 37: Multicellular organismal process. 38: Developmental process. 39: Growth. 40: Locomotion. 41: Single-organism process. 42: Rhythmic process. 43: Response to stimulus. 44: Localization. 45: Multi-organism process. 46: Biological regulation. 47: Cellular component organization or biogenesis.圖2 差異表達基因GO功能富集Fig.2 GO function enrichment of differentially expressed genes

2.2 轉錄組的功能注釋及轉錄組差異表達分析

基于基因在小花蕾期和始花期2個樣品中表達量，對識別到的DEGs進行功能注釋，使用BLAST軟件將Unigene與Pfam，TrEMBL，eggNOG，Swiss-Prot，GO，KEGG數據庫進行序列比對，44 088個Unigene獲得注釋信息，27 801個Unigene被GO數據庫注釋，32 528 個Unigene被eggNOG數據庫注釋，32 781個Unigene被KEGG數據庫注釋。組裝得到的Unigene進行FPKM值表達量統(tǒng)計,2個樣品測得4 349 個DEGs，較小花蕾期1 855個基因表達上調，較始花期2 494個基因表達下調。

2.3 差異表達基因GO功能富集

使用Blast2GO軟件對檢測到的DEGs進行GO功能注釋和富集分析。由圖2可知，DEGs劃分為47個功能組，包括生物學過程、細胞組分、分子功能三大類。其中6 658個GO條目屬于生物學過程，DEGs主要分布于細胞轉化(1 416，9.35%)、代謝過程(1 372，9.06%)、單細胞生物工程(1 101，7.27%)、生物調節(jié)(663，4.38%)、刺激反應(639，4.22%)；4 989個GO條目屬于細胞組成，注釋為細胞組分(1 651，10.90%)、細胞器(989，6.53%)、膜組分(759，5.01%)等相關的DEGs富集較多；3 494個GO條目屬于分子功能，注釋為結合(1 453，9.60%)、催化活性(13 83，9.13%)等功能相關的DEGs較多。

2.4 差異表達基因KEGG通路富集分析

本試驗測序數據中1 219個DEGs被KEGG數據庫注釋，涉及166條代謝途徑。DEGs富集排名前五的KEGG通路為苯丙醇生物合成(phenylpropanoid biosynthesis)(ko00940，68)、淀粉與蔗糖的代謝(starch and sucrose metabolism)(ko00500，41)、植物激素信號轉導(plant hormone signal transduction)(ko04075，40)、氰氨基酸代謝(cyanoamino acid metabolism)(ko00460，28)、α-亞麻酸代謝(alpha-linolenic acid metabolis)(ko00592，25)；與花色形成緊密相關的有類黃酮生物合成途徑(flavonoid biosynthesis)(ko00941，10)、黃酮和黃酮醇生物合成途徑(flavone and flavonol biosynthesis)(ko00944，2)、類胡蘿卜生物合成途徑(carotenoid biosynthesis)(ko00906，6)卟啉和葉綠素代謝(porphyrin and chlorophyll metabolism)(ko00860，1)。

2.5 花色相關基因表達分析

根據差異基因KEGG代謝途徑及基因注釋，從類黃酮、類胡蘿卜素、葉綠素合成代謝途徑，挖掘出與花色相關的54個關鍵結構基因(表3)和21個轉錄調控基因(表4)。由FPKM值可見小花蕾期和始花期2個樣本基因表達的差異，始花期樣品中結構基因有24條單基因較小花蕾期上調，30條單基因下調；類黃酮雙加氧酶基因(naringenin 3-dioxygenase,F3H)、DFR、ANS、矢車菊素 3-O-蕓香糖苷-5-O-葡萄糖苷轉移酶基因(cyanidin 3-O-rutinoside 5-O-glucosyltransferase, 3,5GT)、PDS、ZDS、類胡蘿卜素裂解雙加氧酶基因(cis-β-carotene 9′, 10′-cleaving dioxygenase,CCD7)、葉綠素a合成酶基因(chlorophyll a synthase, CHLG)、葉綠素b還原酶(chlorophyll b reductase，NOL)、葉綠素酶基因(chlorophyllase，CLH)、脫鎂葉綠酸加氧酶基因(pheophorbide a oxygenase，PAO)、順式-環(huán)氧類胡蘿卜素雙加氧酶基因(9-cis-epoxycarotenoid dioxygenase,NCED)均為上調表達；PAL、查耳酮異構酶基因 (chalcone isomerase,CHI)、類黃酮3′-羥化酶基因(flavonoid 3′-monooxygenase,F3′H)、類黃酮-3',5-羥化酶基因(flavonoid 3′,5′-hydroxylase, F3′5′H)、UFGT、花青素糖基轉移酶基因(Anthocyanidin 3-O-glucosyltransferase, 3GT)、PSY、BCH、ZEP、胡蘿卜素順反異構酶基因(carotene cis-trans isomerase, CRTISO)、鎂螯合酶基因(magnesium chelatase subunit H,CHLH)、葉綠素a加氧酶基因(chlorophyllide a oxygenase,CAO)均為下調表達；C4H、香豆酸：輔酶A連接酶基因(4-coumarate—CoA ligase, 4CL)、CHS、FLS、花青素雙糖基轉移酶基因(anthocyanidin5,3-O-glucosyltransferase, 5,3GT)的表達量既有上調也有下調；類胡蘿卜素合成代謝途徑的關鍵酶基因除ZDS和CCD7表達量較高外，其他基因表達量均較?。蝗~綠素合成酶基因僅CHLG表達量上調，葉綠素降解酶基因均上調。小花蕾期CHS(Unigene ID:G31409_c2_g1)基因表達量最高，始花期F3H(Unigene ID:G34942_c5_g1)基因表達量最高。

與小花蕾期樣品相比，始花期樣品的調節(jié)基因有8條單基因上調，13條單基因下調；通過基因注釋，在樣本中篩選出玉米糊粉層有色基因(coloured aleurone,C1)、玉米種皮有色基因(pericarp color,P)、玉米轉錄因子(zeamays1，Zm1)、MYB4、MYB6、外種皮透明基因2(transparent testa 2，TT2)、bHLH13同源性較高的轉錄因子，Zm1轉錄因子的表達量較高；小花蕾期轉錄因子MYB306(Unigene ID:G19067_c0_g1)表達量最高，始花期轉錄因子MYB31(Unigene ID:G12396_c0_g1)表達量最高。

與小花蕾期樣品相比，始花期樣品中花青苷合成下游關鍵結構基因DFR、ANS、3,5GT及轉錄因子Zm1、Hv1、MYB305、MYB31表達量整體較高，且在始花期均為上調。在NCBI數據庫BLAST花青苷下游合成基因及部分表達量較高的轉錄因子同源性物種，由表5可見，與DFR、ANS、3GT基因同源性最高的物種均為大花蕙蘭，與5,3GT、3,5GT、Zm1同源性最高的為蝴蝶蘭或石斛蘭，與轉錄因子MYB31同源性最高的為墨蘭，與轉錄因子Hv1、MYB305、MYB306同源性最高的均為建蘭。

表3 花色關鍵結構基因差異表達Table 3 Differential expression of key structural genes related flower color

表3(續(xù))

表4 花色相關轉錄因子差異表達Table 4 Differential expression of transcription factors related flower color

表5 部分花青苷合成關關鍵基因的最大同源性物種比對Table 5 The maximum identity of partial key genes of the synthesis of anthocyanins among species

2.6 部分差異結構基因轉錄模式驗證

對測序材料轉錄組中類黃酮、類胡蘿卜素代謝途徑的部分結構基因進行qPCR驗證，結果顯示所選基因qPCR結果與轉錄組測序表達結果的變化趨勢一致，說明轉錄組測序結果的可信度較高；基因PAL(Unigene ID:G34234_c1_g1)、CHS(Unigene ID:G25338_c1_g1)、F3H(Unigene ID:G34942_c5_g1)、FLS(Unigene ID:G33942_c5_g2)、DFR(Unigene ID:G31003_c2_g3)在小花蕾期和始花期表達量均較高，其中PAL(Unigene ID:G34234_c1_g1)基因在小花蕾期表達量最高，DFR(Unigene ID:G31003_c2_g3)基因在始花期表達量最高；類胡蘿卜素合成途徑的關鍵結構基因表達量整體明顯低于類黃酮合成途徑的關鍵結構基因表達量(圖3)。

圖3 花瓣轉錄組基因相對表達量的qPCR驗證Fig.3 qPCR verification of gene relative expression in petal transcription group

2.7 花瓣色素的檢測

對試驗材料不同發(fā)育期花瓣的色素含量進行檢測，結果顯示(表6)小花蕾期花瓣花青苷未檢出，始花期總花青苷含量為0.247 9 mg·g-1FW；不同花發(fā)育期各花色素含量均有極顯著差異；總黃酮的含量均高于類胡蘿卜素、葉綠素及總花青苷含量；與小花蕾期花瓣比較，始花期花瓣中類胡蘿卜素、葉綠素及總黃酮含量均降低，僅總花青素苷含量升高。

3 討論

本研究對2個不同發(fā)育期的黃金梅花瓣進行轉錄組測序，2個樣品Mapped Reads在Clean Reads中所占比例分別為86.83%(A)、88.29%(B)；質量值Q20和Q30均大于90%，GC含量均大于47.5%，可見測序數據質量較高。測序共獲得113 780條Unigene，組裝的Transcript與Unigene平均長度為725.64 bp和586.53 bp，N50分別為1 229 bp和882 bp。李文建等[9]對同株建蘭黃綠色和紅色花瓣轉錄組測序獲得106 479條Unigenes，200～300 bp的Unigene比例最高；Transcript與Unigenes的平均長度分別為1 150.14 bp和655.25 bp，N50分別為2 127 bp和1 087 bp，以上測序數據與本研究有一定差異，但GO功能富集分析中各類差異基因占比與本研究結果有一定趨同性。

表6 不同花發(fā)育期花瓣色素含量Table 6 The content of flower pigment in different flower development stages /(mg·g-1 FW)

花青苷與原花色素、黃酮醇等類黃酮物質共享底物及上游代謝路徑，CHS、CHI、F3H、DFR、ANS、FLS、GT基因是類黃酮合成代謝途徑的關鍵結構基因。本試驗測序數據中除ANS基因，其他關鍵酶基因的表達量均較高，推測ANS基因表達量較低是花瓣中花青苷含量大幅低于總黃酮含量的主要原因之一；糖基轉移酶是生成花色苷的關鍵酶[29]，本試驗篩出多個花青素單糖基轉移酶(3GT)及雙糖基轉移酶(3,5GT)基因，特別是矢車菊素 3-O-蕓香糖苷-5-O-葡萄糖苷轉移酶基因(3,5GT)在始花期花瓣中表達量上調明顯。Huai等[30]研究也認為類黃酮糖基轉移酶在開花后期表達強度與花色素合成呈正相關，推測黃金梅始花期花瓣為紫紅色，3,5GT酶基因起到了關鍵作用。本試驗測序數據中PSY、PDS、ZDS、BCH、ZEP、CCD7、CRTISO、NCED類胡蘿卜素合成代謝途徑的關鍵結構基因，除ZDS和CCD7表達量較高外整體表達量較低，CCD7為胡蘿卜素裂解酶，其高表達能導致類胡蘿卜素含量降低，始花期樣品中類胡蘿卜素含量較花青苷含量低，推測是花瓣表觀呈紫紅色的原因。自然界中綠色花較少，但國蘭特別是春蘭、蕙蘭以綠花為主，試驗材料黃金梅母本春蘭也為綠花。而葉綠素是綠色花的基礎，本試驗測序數據中CHLH、CHLG、CAO基因為葉綠素a(chlorophyll a)和葉綠素b(chlorophyll b)的合成酶基因，NOL、CLH、PAO基因為葉綠素降解酶基因，從轉錄組表達量結果來看，葉綠素合成基因僅CHLG表達量上調，而葉綠素降解酶基因均上調，推測從小蕾期到始花期花瓣葉綠素合成減緩、代謝增強，因此始花期葉綠素含量明顯降低。

本研究通過基因注釋篩選出了與C1、P、Zm1、MYB4、MYB6、TT2、bHLH13等轉錄因子序列,已有研究表明這些轉錄因子與花青苷的合成緊密相關[31-38]，C1、P、Zm1是玉米中與花青苷合成相關的轉錄因子，本研究中Zm1轉錄因子表達量較高；此外，與花色素合成代謝有關的Hv1、MYB31、MYB305轉錄因子表達量較高，且在始花期上調，張春玲等[39]認為萬壽菊中Hv1與番茄紅素β-環(huán)化酶基因，即類胡蘿卜素的代謝機制與葉黃素含量的調控有關，Fornalé等[40]研究認為玉米中MYB31參與了類苯丙烷基因的調控；Liu等[41]和Wang等[29]認為MYB305調節(jié)類黃酮代謝基因的表達；本試驗測序數據中MYB306轉錄因子表達水平較高，且在始花期下調，有研究指出，MYB306 的表達參與調控花青苷、總黃酮、酚類物質合成等生理代謝[42-44]，相關轉錄因子的調控功能將在后續(xù)研究進行重點驗證。

本研究材料黃金梅的母本為綠花品種大宋梅，父本為黃花品種黃金虎，雜交后代黃金梅花瓣顯紫紅色，分析其花青苷合成的下游關鍵結構基因和部分轉錄因子，BLAST分析發(fā)現與DFR、ANS、3GT基因同源性最高的物種均為大花蕙蘭，推測基因與父本同源性高；與花青苷雙糖基轉移酶基因(5,3GT、3,5GT)同源性最高的為蝴蝶蘭或石斛蘭，與轉錄因子MYB31同源性最高的為墨蘭，Hv1、MYB305、MYB306同源性最高的均為建蘭，推測部分基因與母本同源性高；黃金梅顯色的原因是否為父母花青苷結構基因和調控基因的互補或重組，可進一步通過春蘭全基因組測序結果作驗證。

4 結論

本研究通對黃金梅不同發(fā)育期花瓣轉錄組測序及數據分析，獲得44 088個Unigene注釋信息；ANS和GT基因的表達在花青苷合成代謝中起到了關鍵調控作用；紅花的呈色主要由花青苷決定，與黃金梅DFR、ANS、3GT基因同源性最高的物種均為大花蕙蘭，與花青苷雙糖基轉移酶基因(5,3GT、3,5GT)同源性最高的為蝴蝶蘭或石斛蘭，與轉錄因子MYB31同源性最高的為墨蘭，與Hv1、MYB305、MYB306同源性最高的均為建蘭。本研究為下一步花色基因的克隆、功能分析及分子調控機理的解析提供了一定參考。