轉(zhuǎn)錄組及代謝組聯(lián)合解析茄子果色上位遺傳效應(yīng)

孫保娟，汪瑞，孫光聞，王益奎，李濤，宮超，衡周，游倩，李植良

轉(zhuǎn)錄組及代謝組聯(lián)合解析茄子果色上位遺傳效應(yīng)

孫保娟1，汪瑞2，孫光聞2，王益奎3，李濤1，宮超1，衡周1，游倩1，李植良1

1廣東省農(nóng)業(yè)科學(xué)院蔬菜研究所/廣東省蔬菜新技術(shù)研究重點(diǎn)實(shí)驗(yàn)室，廣州 510640；2華南農(nóng)業(yè)大學(xué)園藝學(xué)院，廣州 510642；3廣西壯族自治區(qū)農(nóng)業(yè)科學(xué)院蔬菜研究所，南寧 530027

【目的】茄子果色與商品果外觀品質(zhì)和價(jià)值密切相關(guān)，花青素是決定茄子果色的重要天然色素之一。通過(guò)基因表達(dá)和代謝物差異比較，解析上位基因互作調(diào)控茄子果皮花青素合成的作用機(jī)制，為不同果色茄子選育提供理論基礎(chǔ)。【方法】以花青素合成結(jié)構(gòu)基因突變的白果色母本19141、花青素合成關(guān)鍵調(diào)控基因突變的白果色父本19147及其紫紅果色F1代E3316茄子為試驗(yàn)材料，對(duì)商品果期果皮進(jìn)行轉(zhuǎn)錄組測(cè)序和廣靶代謝組分析?！窘Y(jié)果】轉(zhuǎn)錄組測(cè)序分析表明，19141_vs_19147的差異表達(dá)基因（DEGs）最多，其次是E3316_vs_19141，兩個(gè)比較組DEGs均在類(lèi)黃酮途徑富集程度最高；E3316_vs_19147篩選到的DEGs最少，未在類(lèi)黃酮途徑富集。廣靶代謝組共檢測(cè)分析到218個(gè)代謝物。E3316_vs_19141共檢測(cè)到差異代謝物（DAMs）113個(gè)，E3316_vs_19147共檢測(cè)到差異代謝物98個(gè)。轉(zhuǎn)錄組和代謝組聯(lián)合分析發(fā)現(xiàn)花青素生物合成關(guān)鍵結(jié)構(gòu)基因和，關(guān)鍵調(diào)節(jié)基因、和，修飾基因、、和，還有轉(zhuǎn)運(yùn)基因（）相對(duì)表達(dá)量均為19141＞E3316＞19147；果皮呈色相關(guān)的花青素類(lèi)代謝物矢車(chē)菊素（CAS：528-58-5）、矢車(chē)菊素3,5-二葡萄糖苷（CAS：20905-74-2）含量為E3316＞19147＞19141。上位基因調(diào)控下，茄子果皮中同時(shí)表達(dá)與，但是表達(dá)水平遠(yuǎn)高于。E3316花青素含量高于親本，綠原酸含量低于其親本?！窘Y(jié)論】上位性基因互作控制果色的遺傳背景下，親本突變基因類(lèi)型決定了花青素代謝通路基因表達(dá)趨勢(shì)和果皮呈色抑制方式。兩個(gè)白果色親本雜交F1代果皮呈紫紅色是由于花青素生物合成途徑兩個(gè)互作突變基因位點(diǎn)功能互補(bǔ)；高水平表達(dá)是合成矢車(chē)菊型花青素的主要原因；果皮花青素和綠原酸生物合成間存在競(jìng)爭(zhēng)關(guān)系。

茄子；果色；上位基因；花青素；轉(zhuǎn)錄組；代謝組

0 引言

【研究意義】茄子以果實(shí)為產(chǎn)品，是WHO推薦的世界十大健康蔬菜之一。茄子屬于茄科茄屬，在茄科蔬菜生產(chǎn)中占有重要的地位。茄子果皮顏色（以下統(tǒng)稱(chēng)果色）豐富多樣，有紫黑、紫紅、綠、白、粉紅等顏色，不同區(qū)域?qū)x擇具有不同的消費(fèi)習(xí)慣。茄子果色是決定商品外觀品質(zhì)和價(jià)值的重要因素，因此，果色是茄子育種過(guò)程中的一個(gè)重要選育目標(biāo)。筆者在育種過(guò)程中發(fā)現(xiàn)茄子果色上位性遺傳效應(yīng)，即2個(gè)白果色茄子親本雜交，F(xiàn)1為紫紅果，F(xiàn)2紫紅果色和白果色單株分離比符合9﹕7。在此基礎(chǔ)上，從轉(zhuǎn)錄組和代謝組水平對(duì)2個(gè)親本及其F1進(jìn)行比較分析，將有助于從基因表達(dá)和代謝產(chǎn)物差異方面解析茄子果色的上位遺傳效應(yīng)。【前人研究進(jìn)展】100多年前，英國(guó)生物學(xué)家Bateson和Punnett用2個(gè)純白花的豌豆品種雜交，F(xiàn)1后代均為紫花，F(xiàn)2后代中紫花和白花的分離比約9﹕7；為了解釋這種不符合孟德?tīng)柗蛛x比的現(xiàn)象，1909年Bateson提出了上位性的概念，來(lái)定義一個(gè)基因?qū)α硪粋€(gè)基因的遮蓋作用；許多動(dòng)植物重要性狀受2對(duì)以上的上位基因控制[1]。茄子果色遺傳是受多基因控制的復(fù)雜性狀?；ㄇ嗨厥菦Q定茄子果色的主要色素，關(guān)于茄子果色與花青素關(guān)系的遺傳研究中，TIGCHELAAR等[2]發(fā)現(xiàn)茄子果皮內(nèi)花青素由3個(gè)顯性基因（、和）決定，其濃度調(diào)節(jié)則由少數(shù)基因控制。NUNOME等[3]和DOGANLAR等[4]以數(shù)量性狀分析法和現(xiàn)代分子標(biāo)記技術(shù)結(jié)合，成功定位了少數(shù)與果色相關(guān)性較高以及與花青素在果實(shí)內(nèi)合成與含量相關(guān)的QTL，進(jìn)一步證實(shí)了茄子果色遺傳的多基因效應(yīng)。龐文龍等[5]的研究結(jié)果表明，主基因+多基因效應(yīng)是決定茄子果色表型變異的主要因素，環(huán)境因素起較小作用。LIAO等[6]開(kāi)發(fā)了茄子果皮紫紅和紫黑果色SCAR標(biāo)記?；ㄇ嗨亟?jīng)由苯丙烷途徑（phenylpropanoid pathway）和類(lèi)黃酮生物合成途徑（flavonoids biosynthetic pathway）生成[7]。花青素合成的關(guān)鍵結(jié)構(gòu)基因包括苯丙氨酸解氨酶基因、4-香豆酸-CoA連接酶、查爾酮合成酶基因、查爾酮異構(gòu)酶基因、黃烷酮3-羥化酶基因、黃酮3′-羥化酶基因′、黃酮3′,5′-羥化酶基因′′、二氫黃酮醇4-還原酶基因和花青素合成酶基因等。F3′H和F3′5′H決定了花青素合成的分支和類(lèi)型。黃烷酮3-羥化酶（F3H）將柚皮素轉(zhuǎn)化為二氫莰非醇（dihydrokaempferol，DHK），之后若被黃酮3′-羥化酶（F3′H）催化則生成二氫櫟皮黃酮（dihydroquercetin，DHQ），進(jìn)一步被DFR、ANS和UFGT催化生成矢車(chē)菊類(lèi)花青素；若被黃酮3′,5′-羥化酶（F3′5′H）催化則生成二氫楊梅黃酮（dihydromyricetin，DHH），進(jìn)一步被DFR、ANS和UFGT催化生成飛燕草類(lèi)花青素。除了結(jié)構(gòu)基因，轉(zhuǎn)錄調(diào)控因子（transcription factor，TF）也參與調(diào)控花青素合成，其中MYB轉(zhuǎn)錄因子、bHLH轉(zhuǎn)錄因子和WD40轉(zhuǎn)錄因子是發(fā)揮核心作用的TF，它們以MBW（MYB-bHLH-WD repeat）復(fù)合體的形式調(diào)控結(jié)構(gòu)基因的轉(zhuǎn)錄[8-9]。李翔等[10]通過(guò)同源克隆得到茄子花青素5-糖苷轉(zhuǎn)移酶基因（），發(fā)現(xiàn)其在云南紫長(zhǎng)茄和云南圓白茄所有檢測(cè)組織中都有表達(dá)，但在云南紫長(zhǎng)茄的紫色花瓣和商品期紫色果皮中的表達(dá)量略大于云南圓白茄。2016年WANG等[11]綜合生物信息學(xué)和轉(zhuǎn)錄組數(shù)據(jù)分析發(fā)現(xiàn)（Sme2.5_ 05099.1_g00002.1）為參與茄子花青素生物合成的關(guān)鍵轉(zhuǎn)錄因子，其在紫茄的莖、葉、花和果皮中的表達(dá)量顯著高于其在白茄中的表達(dá)量。2015年STOMMEL等[12]在茄子上的研究表明花青素生物合成途徑關(guān)鍵基因表達(dá)水平與花青素含量呈正相關(guān)，紫果皮中表達(dá)量比白果皮中高?；ㄇ嗨乜刂魄炎庸史奂t、淺紫、紫紅、紫黑、甚至黑色，茄子中主要花青素類(lèi)型為飛燕草素[13]?！颈狙芯壳腥朦c(diǎn)】筆者在茄子育種實(shí)踐中發(fā)現(xiàn)2對(duì)上位性效應(yīng)基因控制的茄子果色遺傳現(xiàn)象。白果色父本19147[14]和白果色母本19141[15]雜交F1代E3316果色為紫紅色（圖1），F(xiàn)2代紫紅果色和白果色植株分離比為9﹕7，說(shuō)明其果皮著色受到具有上位性互作的2個(gè)基因位點(diǎn)（和）控制，2個(gè)白果色茄子是由于控制花青素生物合成過(guò)程的2個(gè)獨(dú)立而又互補(bǔ)的上位基因位點(diǎn)分別失活引起。通過(guò)F2群體單株SLAF-Seq構(gòu)建遺傳圖譜，定位了控制茄子果色的上位基因位點(diǎn)[14]和[15]，初步確定和分別為[16]和位點(diǎn)候選基因[17]。然而，兩個(gè)上位基因互作對(duì)花青素代謝途徑基因表達(dá)模式以及代謝產(chǎn)物的影響尚不清楚。【擬解決的關(guān)鍵問(wèn)題】借助轉(zhuǎn)錄組學(xué)和代謝組學(xué)對(duì)白果色父本19147、白果色母本19141以及它們雜交產(chǎn)生的紫紅色果色F1代E3316商品果果皮進(jìn)行差異基因和差異代謝物分析，進(jìn)一步揭示茄子果皮著色的上位性遺傳效應(yīng)。

1 材料與方法

1.1 試驗(yàn)材料

茄子材料為白果色父本19147、白果色母本19141及其雜交得到的紫紅果色F1E3316（圖1）。父本來(lái)源于1998年粵北山區(qū)進(jìn)行蔬菜種質(zhì)資源考察時(shí)收集到的紫花白茄資源YC-212，經(jīng)6代自交純化獲得穩(wěn)定白茄自交系；母本來(lái)源為2009年初在廣州市廣東省農(nóng)業(yè)技術(shù)推廣中心柯木塱示范基地的茄子品種展示田中，收集到的茄子雜交一代組合‘優(yōu)美長(zhǎng)茄’的種果。2009年秋季，種植該雜交品種的種果后代群體共150株，編號(hào)為2911，發(fā)現(xiàn)其后代分離出現(xiàn)紫花紫黑皮、紫花紫褐皮、紫花紫紅皮、白花綠皮、白花白皮等多種類(lèi)型。其中1株綠皮茄子果實(shí)表現(xiàn)較突出，2010年春季種植該綠茄后代，從中分離獲得1個(gè)白花白果色單株，自2010年秋季起，對(duì)該白花白果色株系進(jìn)行定向選育，于2011年秋季獲得白色果皮穩(wěn)定自交系“2911-1-2”，后于2019年改編號(hào)為“19141”。父、母本及其F1于2019年7月育苗，8月種植于廣州市白云區(qū)鐘落潭鎮(zhèn)廣東省農(nóng)業(yè)科學(xué)院基地內(nèi)，10月取商品果期果實(shí)，刮取茄子果皮（厚度約為0.2 mm），稱(chēng)取2 g，用錫紙包裹并用液氮速凍，-80℃保存，用于后續(xù)轉(zhuǎn)錄組測(cè)序分析及代謝組檢測(cè)分析。共3組樣品：19141（白果皮）、19147（白果皮）、E3316（紫紅果皮）。轉(zhuǎn)錄組測(cè)序設(shè)3個(gè)生物學(xué)重復(fù)，共9個(gè)樣品。RNA提取完成并檢測(cè)合格后，委托北京諾禾致源生物信息科技有限公司進(jìn)行轉(zhuǎn)錄組測(cè)序。廣靶代謝組分析設(shè)5個(gè)生物學(xué)重復(fù)，共15個(gè)樣品，由蘇州帕諾米克生物科技有限公司提供廣靶代謝組檢測(cè)分析。

圖1 研究所用茄子材料果實(shí)圖片

1.2 轉(zhuǎn)錄組測(cè)序分析

1.2.1 cDNA文庫(kù)構(gòu)建及轉(zhuǎn)錄組測(cè)序 采用NEB Illumina測(cè)序文庫(kù)制備試劑盒（NEBNext? Ultra? RNA Library Prep Kit for Illumina?）構(gòu)建cDNA文庫(kù)。首先通過(guò)Oligo（dT）磁珠富集mRNA，隨后用二價(jià)陽(yáng)離子在NEB Fragmentation Buffer中將mRNA隨機(jī)打成短片段，以此為模板用隨機(jī)引物在M-MuLV逆轉(zhuǎn)錄酶體系中合成cDNA第一條鏈，然后在DNA polymerase催化下PCR合成cDNA第二條鏈，純化PCR產(chǎn)物，進(jìn)行末端修復(fù)，加A尾和測(cè)序接頭，然后利用AMPure XP beads篩選250—300 bp的片段，再進(jìn)行PCR擴(kuò)增、純化，得到cDNA文庫(kù)，最后使用IlluminaHiSeqTM2500測(cè)序儀進(jìn)行測(cè)序。

1.2.2 基因功能注釋 原始序列過(guò)濾掉接頭序列、低質(zhì)量序列、過(guò)短多N序列等，得到干凈序列（clean reads）。使用HISAT2軟件將clean reads與茄子基因組中數(shù)據(jù)（https://db.cngb.org/search/project/CNP0000734/：國(guó)家基因庫(kù)生命大數(shù)據(jù)平臺(tái)，accession number CNP0000734）進(jìn)行比對(duì)，獲得其在參考基因組定位信息。應(yīng)用BLAST軟件將基因數(shù)據(jù)庫(kù)序列與蛋白質(zhì)數(shù)據(jù)庫(kù)KEGG和GO進(jìn)行序列比對(duì)，得到Unigene的注釋信息。

1.2.3 基因表達(dá)量計(jì)算及差異基因篩選 用FPKM值（每100萬(wàn)條序列中每個(gè)基因以1 000個(gè)堿基為單位，與參考基因組序列比對(duì)上的序列數(shù)量）衡量基因表達(dá)水平。FPKM值越大，基因表達(dá)水平越高。通過(guò)FDR（false discovery rate）設(shè)定閾值的方法來(lái)篩選差異基因。兩組對(duì)比中，差異表達(dá)基因（DEGs）篩選標(biāo)準(zhǔn)為|log2fold change|＞0，F(xiàn)DR≤0.05。

1.3 果皮代謝組的測(cè)定及分析

1.3.1 代謝物的提取 精確稱(chēng)量樣本50 mg（±1%）于2 mL EP管中，準(zhǔn)確加入0.6 mL 2-氯苯丙氨酸（4 mg·L-1，用-20℃甲醇配制），渦旋振蕩30 s；加入100 mg玻璃珠，放入TissueLysisⅡ組織研磨儀中，25 Hz研磨60 s；室溫超聲15 min；12 000 r/min 25℃離心10 min，取上清液300 μL過(guò)0.22 μm過(guò)膜濾，濾液加入到檢測(cè)瓶中，進(jìn)行LC-MS上機(jī)檢測(cè)；從每個(gè)待測(cè)樣本各取20 μL混合成QC樣本，用來(lái)校正混合樣品分析結(jié)果的偏差以及由于分析儀器自身原因所造成的失誤；用剩余待測(cè)樣本進(jìn)行LC-MS檢測(cè)。

1.3.2 LC-MS液相質(zhì)譜分析 本試驗(yàn)采用ACQUITY UPLC? HSS T3（2.1 mm×150 mm，1.8 μm）色譜柱進(jìn)行分析，自動(dòng)進(jìn)樣器溫度設(shè)為8℃，以0.25 mL?min-1的流速，柱溫40℃，進(jìn)樣2 μL進(jìn)行梯度洗脫，流動(dòng)相組成為A：0.1%甲酸水+0.1%甲酸乙腈，B：5 mmol?L-1甲酸銨水+乙腈。質(zhì)譜分析使用Thermo Q Exactive Focus質(zhì)譜儀，采用正離子和負(fù)離子檢測(cè)模式的ESI電離方式，質(zhì)量掃描范圍為81—1 000 m/z。

1.3.3 質(zhì)譜數(shù)據(jù)分析 質(zhì)譜掃描獲得的原始數(shù)據(jù)通過(guò)ProteoWizard軟件轉(zhuǎn)換成mzXML格式，然后利用R（v3.3.2）語(yǔ)言的XCMS程序包進(jìn)行峰識(shí)別、峰過(guò)濾、峰對(duì)齊，得到質(zhì)核比、保留時(shí)間和峰面積等信息組成的數(shù)據(jù)矩陣；將正離子模式獲得的29 461個(gè)前體分子和負(fù)離子模式獲得的19 340個(gè)前體分子數(shù)據(jù)，導(dǎo)出至Excel進(jìn)行后續(xù)分析。對(duì)數(shù)據(jù)進(jìn)行峰面積的批次歸一化，以便不同量級(jí)的初始數(shù)據(jù)比較。利用SIMCA-P軟件進(jìn)行PCA-DA（偏最小二乘法判別分析）主成分分析。

1.3.4 數(shù)據(jù)分析和圖片處理 采用一維方差分析，-value≤0.05且VIP≥1相結(jié)合的方法篩選差異代謝物，并計(jì)算出差異代謝物在兩組間的差異變化倍數(shù)FC（fold change），通過(guò)差異變化倍數(shù)判斷差異代謝物的上調(diào)或者下調(diào)。用R語(yǔ)言（v3.3.2）中pheatmap程序包對(duì)數(shù)據(jù)集進(jìn)行縮放，得到代謝物相對(duì)定量值層次聚類(lèi)圖。借助MetPA數(shù)據(jù)庫(kù)分析，識(shí)別差異代謝物相關(guān)的KEGG通路，根據(jù)代謝通路中鑒定代謝物的相對(duì)響應(yīng)值和降維算法，獲得代謝通路的相對(duì)響應(yīng)值。

2 結(jié)果

2.1 親本及其F1代果皮轉(zhuǎn)錄組分析

2.1.1 親本及其F1代果皮差異表達(dá)基因分析 轉(zhuǎn)錄測(cè)序結(jié)果顯示，在供試茄子父、母本及其F1代之間的3個(gè)比較組中，19141_vs_19147即母本和父本相比，篩選到的差異基因最多，共3 064個(gè)，其中上調(diào)基因1 436個(gè)，下調(diào)基因1 628個(gè)；其次是E3316_vs_19141，即F1與母本相比，篩選到差異基因共2 531個(gè)，其中上調(diào)基因1 340個(gè)，下調(diào)基因1 191個(gè)；E3316_vs_19147，即F1代與父本相比，篩選到的差異基因最少，共1 213個(gè)，其中上調(diào)651個(gè)，下調(diào)562個(gè)（圖2）。

2.1.2 親本及其F1代果皮差異表達(dá)基因KEGG富集情況 對(duì)3個(gè)比較組19141_vs_19147、E3316_vs_19141和E3316_vs_19147的DEGs進(jìn)行KEGG富集分析發(fā)現(xiàn)，19141_vs_19147的DEGs顯著富集于類(lèi)黃酮生物合成、芪類(lèi)二苯基庚酮生物合成代謝途徑；E3316_ vs_19141的DEGs顯著富集于類(lèi)黃酮生物合成、芪類(lèi)二苯基庚酮生物合成、RNA聚合酶、苯丙烷生物合成代謝途徑；E3316_vs_19147的DEGs沒(méi)有顯著性富集的代謝途徑（圖3）。

圖2 不同對(duì)比組差異表達(dá)基因統(tǒng)計(jì)

2.1.3 類(lèi)黃酮代謝途徑相關(guān)基因在茄子親本及其F1果皮表達(dá)情況 類(lèi)黃酮生物合成途徑的一個(gè)重要分支是合成花青素，與植物各個(gè)組織部位呈色密切相關(guān)。由于E3316_vs_19147的DEGs沒(méi)有在類(lèi)黃酮途徑發(fā)生顯著富集，因此僅對(duì)E3316_vs_19141和19141_vs_19147兩個(gè)比較組在類(lèi)黃酮生物途徑顯著富集基因進(jìn)行分析（表1）。E3316_vs_19141在類(lèi)黃酮生物途徑富集了15個(gè)差異基因，其中4個(gè)基因上調(diào)，為綠原酸支路的羥基肉桂酰轉(zhuǎn)移酶基因（EGP13381、EGP07013、EGP32043和EGP10080）；11個(gè)基因下調(diào)，包括花青素代謝通路關(guān)鍵基因（EGP22200）、（EGP24232）、（EGP30923）、（EGP32037）、、（EGP18904）和花青素O-甲基轉(zhuǎn)移酶基因（EGP16482），綠原酸支路的基因（EGP05102、EGP30738、EGP01163）和羥基肉桂酰輔酶A奎尼羥基肉桂轉(zhuǎn)移酶基因（EGP02613）。由于花青素代謝通路的在轉(zhuǎn)錄組數(shù)據(jù)分析結(jié)果中未發(fā)生富集，因此根據(jù)已克隆的序列與參考基因組序列進(jìn)行對(duì)比，發(fā)現(xiàn)參考基因組注釋的（EGP31017）存在嵌合基因，影響了基因定量和富集結(jié)果，因此表1列出的數(shù)據(jù)為校正后結(jié)果。白茄親本比較組19141_vs_19147差異基因在類(lèi)黃酮生物途徑發(fā)生了顯著性富集，其中11個(gè)基因上調(diào)，包括、、、、、、和（EGP16482），綠原酸支路的反式肉桂酸4-單加氧酶基因（EGP24021）、（EGP05102和EGP30738）和（EGP16482）；4個(gè)基因下調(diào)，包括花青素代謝通路（EGP06309）和綠原酸支路的（EGP13381）、（EGP07013和EGP32043）。

表1 不同比較組類(lèi)黃酮生物合成途徑差異基因富集情況

比較組差異表達(dá)基因篩選標(biāo)準(zhǔn)為|log2倍數(shù)變化|＞0，F(xiàn)DR≤0.05；“-”表示該比較組中該基因未發(fā)生顯著富集；E3316_vs_19141比較組沒(méi)有顯著富集的差異基因，所以未列出

The screening criteria of DEGs in the comparison group were |log2 fold change|＞0, FDR≤0.05. “-” indicate that the gene was not significantly enriched in the comparison group. E3316_vs_19141 had no significantly enriched DEGs, so it was not listed

2.2 代謝組分析

2.2.1 代謝組的多元統(tǒng)計(jì)分析 利用SIMCA-P軟件對(duì)測(cè)得的白果色茄子母本19141、白果色茄子父本19147及其紫紅色F1代E3316的果皮代謝組數(shù)據(jù)主成分分析（principal component analysis，PCA）結(jié)果表明，無(wú)論是正離子模式下還是負(fù)離子模式下，白果色茄子母本19141、白果色茄子父本19147及其紫紅色F1代E3316等3組樣本數(shù)據(jù)產(chǎn)生明顯分離，且各組生物學(xué)重復(fù)樣本緊密地聚集在一起（圖4）。正離子模式下，第一主成分（PC1）和第二主成分（PC2）分別占總變量的36.6%和21.4%；負(fù)離子模式下，第一主成分（PC1）和第二主成分（PC2）分別占總變量的39.4%和21.4%。PC1值越大，說(shuō)明不同品種間的遺傳變異程度越高。主成分分析結(jié)果顯示3個(gè)茄子材料間存在顯著差異。

2.2.2 樣本間代謝物差異及代謝通路分析 本試驗(yàn)利用液相色譜-質(zhì)譜聯(lián)用技術(shù)，共檢測(cè)到茄子果皮代謝物218個(gè)。依據(jù)OPLS-DA模型第一主成分VIP值＞1、-value＜0.05，篩選得到不同樣本組間（E3316_vs_ 19141和E3316_vs_19147）的DAMs。E3316_vs_19141比較組共檢測(cè)到113個(gè)DAMs，其中54個(gè)上調(diào)，59個(gè)下調(diào)，注釋到42個(gè)代謝通路上，影響因子排名前5位的通路分別為-苯丙酸代謝、異喹啉生物堿生物合成、酪氨酸代謝、-亞麻酸代謝、泛酸和CoA生物合成（圖5）；E3316_vs_19147比較組共檢測(cè)到98個(gè)DAMs，其中24個(gè)上調(diào)，74個(gè)下調(diào)，注釋到49個(gè)代謝通路上，影響因子排名前5位的通路分別為苯丙氨酸代謝，單帖生物合成，異喹啉生物堿生物合成，丙氨酸、天冬氨酸和谷氨酸代謝及半胱氨酸和蛋氨酸代謝（圖5）。

2.2.3 樣本間差異代謝物層次聚類(lèi)分析及類(lèi)黃酮生物合成途徑差異代謝物 根據(jù)差異代謝物的相對(duì)含量，對(duì)其進(jìn)行層次聚類(lèi)分析。從19141、19147和E3316樣本間84個(gè)DAMs聚類(lèi)分析結(jié)果可見(jiàn)（圖6），白果色父本19147和紫紅果色E3316絕大多數(shù)DAMs含量更相近，它們要么顯著高于白果色母本19141，要么顯著低于白果色母本19141?？傮w來(lái)看，3組比較樣本代謝物差異趨勢(shì)與基因表達(dá)差異趨勢(shì)相似，都表現(xiàn)為母本比父本較E3316的差異大。與兩個(gè)白果色親本相比，紫紅果色F1E3316有19種代謝物上調(diào)，包括與果皮呈紫紅色相關(guān)的矢車(chē)菊素（cyanidin，CAS：528-58-5/13306-05-3）、矢車(chē)菊素3,5-二葡萄糖苷（cyanin，CAS：20905-74-2/2611-67-8），苯丙烷類(lèi)代謝物香豆素（coumarin）；-丙氨酸代謝相關(guān)的代謝物-丙氨酸、精胺（spermine）；生物堿假木賊堿（Anabasine）等；顯著下調(diào)的類(lèi)黃酮代謝物包括對(duì)香豆?？鼘幩幔╬-Coumaroyl quinic acid）和綠原酸（Chlorogenic acid，CGA，CAS：327-97-9）。

2.3 轉(zhuǎn)錄組和代謝組聯(lián)合分析上位性基因控制的果皮花青素合成

白果色父母本與其紫紅色F1果皮呈色不同主要是由于花青素引起的，因此，對(duì)花青素生物合成途徑相關(guān)基因表達(dá)和代謝物進(jìn)行聯(lián)合分析（圖7）。根據(jù)已知花青素生物合成途徑相關(guān)結(jié)構(gòu)基因、調(diào)控基因、修飾和轉(zhuǎn)運(yùn)基因（圖7-A），基于19141、19147及E3316果皮轉(zhuǎn)錄組相應(yīng)基因FPKM值構(gòu)建表達(dá)熱圖（圖7-B），由圖7-B可見(jiàn)，、、、、′′、和等花青素合成關(guān)鍵結(jié)構(gòu)基因，花青素3-糖苷轉(zhuǎn)移酶（anthocyanin 3-glycoside transferase，），花青素5-糖苷轉(zhuǎn)移酶（anthocyanin 5-glycoside transferase，），花青素?；D(zhuǎn)移酶（anthocyanin acyltransferase，），花青素-甲基轉(zhuǎn)移酶（anthocyanin-methyl transferase，）等花色苷形成和修飾基因，以及花青素轉(zhuǎn)運(yùn)基因（）的表達(dá)都為19141高于E3316，E3316高于19147。花青素合成途徑關(guān)鍵調(diào)節(jié)基因、和的表達(dá)模式與結(jié)構(gòu)基因、修飾基因及轉(zhuǎn)運(yùn)基因相似，表現(xiàn)為19141高于E3316，E3316高于19147。在3個(gè)樣品中相對(duì)表達(dá)量均較高，且顯著高于，這可能是上位性基因調(diào)控下合成的花青素是矢車(chē)菊素類(lèi)型的主要原因。

從圖7-C看，廣靶代謝組共檢測(cè)到3種花青素合成相關(guān)的代謝產(chǎn)物，分別為柚皮素、矢車(chē)菊素和矢車(chē)菊糖苷；與茄子果皮著色密切相關(guān)的代謝物矢車(chē)菊素和矢車(chē)菊素3,5-二葡萄糖苷相對(duì)含量表現(xiàn)為E3316＞19147＞19141。

綜合轉(zhuǎn)錄組和代謝組數(shù)據(jù)來(lái)看，兩個(gè)白果色親本雜交后，E3316花青素代謝通路相關(guān)基因表達(dá)水平居中，能夠進(jìn)行正常的轉(zhuǎn)錄和翻譯，催化花青素合成及轉(zhuǎn)運(yùn)，因此，E3316果皮呈紫紅色。

圖4 19141、19147和E3316茄子果皮代謝物PCA得分圖

圖5 差異代謝物的KEGG通路富集氣泡圖

紅色框內(nèi)代謝物為E3116較親本上調(diào)，綠色框內(nèi)代謝物為E3316較親本下調(diào)

A：花青素生物合成途徑圖解；B：根據(jù)花青素生物合成途徑相關(guān)基因FPKM值構(gòu)建的熱圖；C：廣靶代謝組檢測(cè)到的花青素生物合成相關(guān)代謝物結(jié)構(gòu)及相對(duì)含量

3 討論

3.1 果色上位性基因調(diào)控下基因表達(dá)與代謝產(chǎn)物間的關(guān)系

本研究采用的白果色父本材料為突變體[16]，白果色母本材料19141為突變體[17]。主要作用是將無(wú)色的原花色素氧化產(chǎn)生有色的花青素，該產(chǎn)物是花青素合成途徑中的第一個(gè)有色產(chǎn)物，對(duì)植物顯色形成至關(guān)重要的作用[18]。茄子上的研究表明表達(dá)水平與花青素含量呈正相關(guān)，紫果皮比綠果皮[19]和白皮[20]的表達(dá)量高。本研究上位基因互作背景下，白果色母本19141果皮中表達(dá)量顯著高于紫紅果色E3316。是控制茄子花青素合成的重要轉(zhuǎn)錄因子之一，對(duì)大多數(shù)促進(jìn)花青素生物合成的轉(zhuǎn)錄因子而言，都會(huì)上調(diào)結(jié)構(gòu)基因表達(dá)，從而導(dǎo)致大量花青素積累[21-22]。研究表明除調(diào)控花青素生物合成，還調(diào)控綠原酸CGA生物合成，刪除端的9個(gè)三聯(lián)體密碼子，突變體花青素合成途徑被破壞，CGA還可以積累[23]。

在2對(duì)雙隱性上位基因互作的遺傳背景下，花青素合成途徑基因的相對(duì)表達(dá)量與代謝產(chǎn)物含量間的關(guān)系非孟德?tīng)栠z傳條件下的簡(jiǎn)單對(duì)應(yīng)關(guān)系。ZHANG等[20]研究發(fā)現(xiàn)、、、、、、等7個(gè)結(jié)構(gòu)基因和、、、等4個(gè)轉(zhuǎn)錄因子在紫茄皮的表達(dá)量都顯著高于白茄皮。本研究中白果色母本19141果皮中花青素生物合成途徑相關(guān)基因表達(dá)量顯著高于紫紅果色F1E3316。與茄子果皮著色密切相關(guān)的代謝物氯化矢車(chē)菊素和矢車(chē)菊素3,5-二葡萄糖苷相對(duì)含量表現(xiàn)為E3316＞19147＞19141。這可能是由于19141突變方式為花青素合成關(guān)鍵結(jié)構(gòu)基因發(fā)生無(wú)義突變，終止密碼子提前，雖然整個(gè)通路的基因表達(dá)未受到影響，但由于基因翻譯提前終止，無(wú)法行使功能，導(dǎo)致代謝通路中斷，因此矢車(chē)菊色素和矢車(chē)菊雙甙含量最低；19147突變方式為花青素合成關(guān)鍵調(diào)控基因，其下調(diào)表達(dá)使花青素生物合成途徑結(jié)構(gòu)基因整體下調(diào)，所以各基因相對(duì)表達(dá)量均最低，但還有一定量的花青素合成，可能由于表達(dá)量也很低而無(wú)法行使把花色苷轉(zhuǎn)運(yùn)到液泡呈色的功能。E3316中和兩個(gè)上位性基因位點(diǎn)呈雜合，基因功能互補(bǔ)，花青素合成通路關(guān)鍵結(jié)構(gòu)和調(diào)控基因的表達(dá)量居中，花青素合成恢復(fù)，果皮呈紫紅色。

3.2 上位性基因調(diào)控下果皮合成花青素類(lèi)型

茄子果色與茄子果皮中的色素種類(lèi)和含量相關(guān)。茄子果皮一般含花青素和葉綠素兩種色素，花青素控制茄子果色呈粉紅、淺紫、紫紅、紫黑，甚至黑色，葉綠素控制茄子果色呈綠色[24]。紫色茄子中主要花青素為飛燕草素3-p-香豆酰蕓香糖苷5-葡萄糖苷（CAS: 28463-30-1，nasunin，茄色甙）；飛燕草素3-蕓香糖苷（CAS：58585-26-0，D3R，山慈姑花苷），少數(shù)茄子和部分野生茄子的花青素中具有飛燕草素-3-葡萄糖苷（CAS：6906-38-3）、矮牽牛素3-（p-香豆酰鼠李糖苷）-5-葡萄糖苷、飛燕草素3-咖啡?；罄钐擒?5-葡萄糖苷等[25-26]。本研究中檢測(cè)到的花青素類(lèi)代謝物為矢車(chē)菊素3,5-二葡萄糖苷（CAS：20905-74-2/2611- 67-8，cyanin），其使E3316果皮表現(xiàn)為紫紅色，這與以往報(bào)道中茄子主要花青素為飛燕草素衍生物有所不同。在葡萄上的研究表明，果皮中同時(shí)表達(dá)與，它們的比例控制著葡萄皮中花青素的組成，表達(dá)水平高的品種，果皮中矢車(chē)菊類(lèi)花青素累積多，表達(dá)水平高的品種，飛燕草素類(lèi)花青素累積多[27]。本研究茄子果皮相對(duì)表達(dá)水平高于，這可能是檢測(cè)到的花青素種類(lèi)為矢車(chē)菊素的重要原因。

3.3 果色上位性基因調(diào)控下花青素和綠原酸生物合成間的關(guān)系

花青素和綠原酸生物合成都是經(jīng)由苯丙素代謝途徑合成的，其在類(lèi)黃酮合成途徑的p-Coumaroyl-CoA處發(fā)生分支，分別被CHS和HCT酶催化進(jìn)入花青素和綠原酸生物合成途徑。綠原酸CGA是茄子果實(shí)中主要的抗氧化酚類(lèi)代謝物[28-29]，其在茄子中含量分別比番茄和馬鈴薯高出10倍和100倍[30]。在馬鈴薯上的研究發(fā)現(xiàn)通過(guò)RNAi抑制綠原酸代謝途徑的HQT基因降低綠原酸含量，增加了花青素含量[31]。本研究中與綠原酸合成相關(guān)的代謝物對(duì)香豆酰奎寧酸和綠原酸CGA含量趨勢(shì)為白果色親本高于紫紅果色E3316，這可能是由于CGA合成與花青素合成存在底物競(jìng)爭(zhēng)關(guān)系，當(dāng)紫紅果色E3316花青素合成恢復(fù)后，CGA合成受到競(jìng)爭(zhēng)性調(diào)控而含量下調(diào)。

4 結(jié)論

2對(duì)隱性上位基因互作控制茄子果色的遺傳背景下，親本突變基因的類(lèi)型決定了整個(gè)花青素代謝通路基因表達(dá)的趨勢(shì)和花青素類(lèi)代謝物合成及呈色的抑制方式。親本突變的基因位點(diǎn)為結(jié)構(gòu)基因，整個(gè)花青素合成通路基因表達(dá)不會(huì)受到影響，但關(guān)鍵酶功能缺失造成代謝通路中斷阻礙花青素積累；親本突變的基因位點(diǎn)為MYB轉(zhuǎn)錄因子基因，其幾乎下調(diào)整個(gè)花青素合成通路基因，雖然可以合成一定量的花色苷，但可能由于轉(zhuǎn)運(yùn)受阻而不能在液泡呈色。兩個(gè)白果色親本雜交后代果皮呈紫紅色可能是花青素生物合成途徑兩個(gè)不同突變基因位點(diǎn)功能互補(bǔ)的結(jié)果；表達(dá)水平顯著高于，促進(jìn)矢車(chē)菊型花青素合成；果皮花青素和綠原酸生物合成間存在競(jìng)爭(zhēng)關(guān)系。

[1] MIKO B I, WRITE P D, RIGHT S, EDUCATION N. Epistasis: Gene interaction and phenotype effects. Nature Education, 2011, 6(2008): 1-6.

[2] TIGCHELAAR E C, JANICK J, ERICHSON H T. The genetics of anthocyanin coloration in eggplant (L.). Genetics, 1968, 60: 475-491. doi: 10.1093/genetics/60.3.475.

[3] NUNOME T, ISHIGURO K, YOSHIDA T, HIRAI M. Mapping of fruit shape and color development traits in eggplant (L.) based on RAPD and AFLP markers. Breeding Science, 2011, 51: 19-26. doi: 10.1270/jsbbs.51.19.

[4] DOGANLAR S, FRARY A, DAUNAY M C, LESTER R N, TANKSLEY S D. Conservation of gene function in the Solanaceae as revealed by comparative mapping of domestication traits in eggplant. Genetics, 2002, 161(4): 1713-1726. doi: 10.1093/genetics/161.4.1713.

[5] 龐文龍, 劉富中, 陳鈺輝, 連勇. 茄子果色性狀的遺傳研究. 園藝學(xué)報(bào), 2008, 35(7): 979-986.

PANG W L, LIU F Z, CHEN Y H, LIAN Y. Genetic study on fruit color traits of eggplant. Acta Horticulturae Sinica, 2008, 35(7): 979-986. (in Chinese)

[6] LIAO Y, SUN B J, SUN G W, LIU H C, LI Z L, LI Z X, WANG G P, CHEN R Y. AFLP and SCAR markers associated with peel color in eggplant (L.). Scientia Agricultura Sinica, 2009, 8(12): 1466-1474.

[7] HOLTON T A, CORNISH E C. Genetics and biochemistry of anthocyanin biosynthesis. The Plant Cell, 1995, 7(7): 1071-1083. doi: 10.1105/tpc. 7.7.1071.

[8] 宮硤, 薛靜, 張曉東. 植物花青素合成途徑中的調(diào)控基因研究進(jìn)展. 生物技術(shù)進(jìn)展, 2011, 1(6): 381-390.

GONG X, XUE J, ZHANG X D. Regulation genes in plant anthocyanin synthesis pathway. Current Biotechnology, 2011, 1(6): 381-390. (in Chinese)

[9] ALBERT N W, DAVIES K M, LEWIS D H, ZHANG H B, MONTEFIORI M, BRENDOLISE C, BOASE M R, NGO H, JAMESON P E, SCHWINN K E. A conserved network of transcriptional activators and repressors regulates anthocyanin pigmentation in eudicots. The Plant Cell, 2014, 26(3): 962-980. doi: 10.1105/tpc.113.122069.

[10] 李翔, 劉楊, 李懷志, 王中彥, 陳火英. 茄花青素5-O糖基轉(zhuǎn)移酶基因克隆與表達(dá)特征分析. 上海交通大學(xué)學(xué)報(bào)(農(nóng)業(yè)科學(xué)版), 2011, 29(6): 1-5, 23.

LI X, LIU Y, LI H Z, WANG Z Y, CHEN H Y. Cloning and expression characterization of the anthocyanin 5-O glucosyltransferase in eggplant. Journal of Shanghai Jiao Tong University (Agricultural Science), 2011, 29(6): 1-5, 23. (in Chinese)

[11] WANG S J, CHEN Z, JI T, DI Q H, LI L J, WANG X F, WEI M, SHI Q H, LI Y, GONG B, YANG F J. Genome-wide identification and characterization of the R2R3MYB transcription factor superfamily in eggplant (L.). Agri Gene, 2016, 2: 38-52.doi: 10.1016/j.aggene.2016.09.006

[12] STOMMEL J R, DUMM J M. Coordinated regulation of biosynthetic and regulatory genes coincides with anthocyanin accumulation in developing eggplant fruit. Journal of the American Society for Horticultural Science, 2015, 140(2): 129-135.

[13] SAKAMURA S, WATANABE S, OBATA Y. The structure of the major anthoeyanin in eggplant. Agricultural and Biological Chemistry, 1963, 27(9): 663-665. doi: 10.1080/00021369.1963.10858161.

[14] 孫保娟, 李植良, 黎振興, 羅少波, 徐曉美. 茄子果色上位性基因的定位及其InDel分子標(biāo)記開(kāi)發(fā)與應(yīng)用. 中國(guó), ZL2016 1 1128459.5. 2018-03-02.

SUN B J, LI Z L, LI Z X, LUO S B, XU X M. Mapping of eggplant fruit color epistatic gene, and InDel molecular marker development and application. China, ZL2016 1 1128459.5. 2018-03-02. (in Chinese)

[15] 孫保娟, 李植良, 黎振興, 羅少波, 吳俊. 與茄子果色上位性基因緊密連鎖的InDel分子標(biāo)記及應(yīng)用.中國(guó), ZL201810018063.8. 2018-10-02.

SUN B J, LI Z L, LI Z X, LUO S B, WU J. InDel molecular marker closely linked with epistatic geneof eggplant fruit color and application of InDel molecular marker. China, ZL201810018063.8. 2018-10-02. (in Chinese)

[16] 李慧敏. 茄子果色上位D基因[D]. 廣州: 華南農(nóng)業(yè)大學(xué), 2020.

Li H M. Cloning of epigenetic D genecontrolling peel color and its function in anthocyanin synthesis in eggplant [D]. Guangzhou: South China Agricultural University, 2020. (in Chinese)

[17] 吳俊. 茄子果色上位基因基因分子標(biāo)記開(kāi)發(fā)及候選基因分析[D]. 廣州: 華南農(nóng)業(yè)大學(xué), 2020.

Wu J. Molecular marker development and candidate gene analysis of epistasisgene controlling fruit color in eggplant (L.) [D]. Guangzhou: South China Agricultural University, 2020. (in Chinese)

[18] SPRINGOB K, NAKAJIMA J, YAMAZAKI M, SAITO K. Recent advances in the biosynthesis and accumulation of anthocyanins. Natural Product Reports, 2003, 20(3): 288-303. doi: 10.1039/ b109542k.

[19] GISBERT C, DUMM J M, PROHENS J, VILANOVA S，STOMMEL J R. A spontaneous eggplant (L.) color mutant conditions anthocyanin-free fruit pigmentation. Hortscience, 2016, 51(7): 793-798.doi: 10.21273/HORTSCI.51.7.793.

[20] ZHANG Y J, HU Z L, CHU G H, HUANG C, TIAN S B, ZHAO Z P, CHEN G P. Anthocyanin accumulation and molecular analysis of anthocyanin biosynthesis-associated genes in eggplant (L.). Journal of Agricultural and Food Chemistry, 2014, 62(13): 2906-2912.

[21] 邵文婷, 劉楊, 韓洪強(qiáng), 陳火英. 茄子花青素合成相關(guān)基因SmMYB的克隆與表達(dá)分析. 園藝學(xué)報(bào), 2013, 40(3): 467-478.

SHAO W T, LIU Y, HAN H Q, CHEN H Y. Cloning and expression analysis of an anthocyanin-related transcription factor gene SmMYB in eggplant. Acta Horticulturae Sinica, 2013, 40(3): 467-478. (in Chinese)

[22] 王海竹, 曲紅云, 周婷婷, 徐啟江. 茄萼花色苷合成相關(guān)基因DFR和MYB克隆及表達(dá)分析. 中國(guó)農(nóng)業(yè)科學(xué), 2017, 50(14): 2781-2792.

WANG H Z, QU H Y, ZHOU T T, XU Q J. Cloning and expression analysis of anthocyanin biosynthesis-associated DFR and MYB genes inof eggplant (L.). Scientia Agricultura Sinica, 2017, 50(14): 2781-2792. (in Chinese)

[23] DOCIMO T, FRANCESE G, RUGGIERO A, BATELLI G, DE PALMA M, BASSOLINO L, TOPPINO L, ROTINO G L, MENNELLA G, TUCCI M. Phenylpropanoids accumulation in eggplant fruit: characterization of biosynthetic genes and regulation by a MYB transcription factor. Frontiers in Plant Science, 2015, 6: 1233. doi: 10.3389/fpls.2015.01233.

[24] DAUNAY M C, AUBERT S, FRARY A, DOGANLAR S, LESTER R N, BARENDSE G, VAN DER WEERDEN G, HENNART J W, HAANSTRA J, DAUPHIN F, JULLIAN E. Eggplant () fruit color: pigments, measurements and genetics// Proceeding of the 12 Eucarpia Meeting on Genetics and Breeding of Capsicum and Eggplant, 2004: 108-116.

[25] AZUMA K, OHYAMA A, IPPOUSHI K, ICHIYANAGI T, TAKEUCHI A, SAITO T, FUKUOKA H. Structures and antioxidant activity of anthocyanins in many accessions of eggplant and its related species. Journal of Agricultural and Food Chemistry, 2008, 56(21): 10154-10159. doi: 10.1021/jf801322m.

[26] 張映, 趙悅琪, 陳鈺輝, 連勇, 劉富中. 茄子紫色形成的分子研究進(jìn)展. 園藝學(xué)報(bào), 2019, 46(9): 1779-1796.

ZHANG Y, ZHAO Y Q, CHEN Y H, LIAN Y, LIU F Z. Progress in molecular research on purple formation of eggplant. Acta Horticulturae Sinica, 2019, 46(9): 1779-1796. (in Chinese)

[27] JEONG S T, GOTO-YAMAMOTO N, HSSHIZUME K, ESAKA M. Expression of the flavonoid 3’-hydroxylase and flavonoid 3’5’ hydroxylase genes and flavonoid composition in grape (). Plant Science, 2006, 17(1): 61-69. doi: 10.1016/j.plantsci.2005.07. 025.

[28] GRAMAZIO P, PROHENS J, PLAZAS M, ANDúJAR I, HERRAIZ F J, CASTILLO E, KNAPP S, MEYER R S, VILANOVA S. Location of chlorogenic acid biosynthesis pathway and polyphenol oxidase genes in a new interspecific anchored linkage map of eggplant. BMC Plant Biology, 2014, 14: 350. doi: 10.1186/s12870-014-0350-z.

[29] MENNELLA G, LO SCALZO R, FIBIANI M, D'ALESSANDRO A, FRANCESE G, TOPPINO L, ACCIARRI N, DE ALMEIDA A E, ROTINO G L. Chemical and bioactive quality traits during fruit ripening in eggplant (L.) and allied species. Journal of Agricultural and Food Chemistry, 2012, 60(47): 11821- 11831. doi: 10.1021/jf3037424.

[30] DOCIMO T, FRANCESE G, RUGGIERO A, BATELLI G, DE PALMA M, BASSOLINO L, TOPPINO L, ROTINO G L, MENNELLA G, TUCCI M. Phenylpropanoids accumulation in eggplant fruit: characterization of biosynthetic genes and regulation by a MYB transcription factor. Frontiers in Plant Science, 2015, 6: 1233. doi: 10.3389/fpls.2015.01233.

[31] PAYYAVULA R S, SHAKYA R, SENGODA V G, MUNYANEZA J E, SWAMY P, NAVARRE D A. Synthesis and regulation of chlorogenic acid in potato: Rerouting phenylpropanoid flux in HQT- silenced lines. Plant Biotechnology Journal, 2015, 13(4): 551-564. doi: 10.1111/pbi.12280.

Transcriptome and Metabolome Integrated Analysis of Epistatic Genetics Effects on Eggplant Peel Color

1Vegetable Research Institute, Guangdong Academy of Agricultural Sciences/Guangdong Key Laboratory for New Technology Research of Vegetables, Guangzhou 510640;2College of Horticulture, South China Agricultural University, Guangzhou 510642;3Institute of Vegetable, Guangxi Academy of Agricultural Sciences, Nanning 530007

【Objective】Peel color is closely related to the appearance and value of eggplant. Anthocyanin is one of the natural pigments that determines an eggplant’s peel color. Through the comparison of gene expression and metabolites as well as the analysis on the mechanism of epistatic gene interaction and regulation of anthocyanin synthesis in eggplant peel, this paper provided a theoretical basis for eggplant breeding with different peel colors.【Method】The white-peel female parent 19141 with mutation at the structural geneof anthocyanin biosynthesis pathway, white-peel male parent 19147 with mutation at regulatory geneof anthocyanin biosynthesis pathway, and their F1hybridE3316 with reddish-purple-peel were used as test materials, while transcriptome sequencing and wide-targeted metabolome analysis were performed on peels of commercial eggplant.【Result】The transcriptome sequencing analysis showed that: 19141_vs_19147 had the most differentially expressed genes (DEGs), followed by E3316_vs_19141. The two comparison groups both had the significant enrichment of DEGs in the flavonoid pathway. The DEGs of E3316_vs_19147 were the least and were not enriched in the flavonoid pathway. A total of 218 metabolites were detected by wide-targeted metabolomeanalysis. A total of 113 differentialaccumulated metabolites (DAMs) were detected in E3316_vs_19141, and total 98 DAMs were detected in E3316_vs_19147. The combined analysis of transcriptome and metabolome found that the relative expression levels of key structural genes,,,and, key regulatory genes,and, modifier genes,,and, and transporter gene() involving anthocyanin biosynthesis pathway had the relation of 19141＞E3316＞19147. The contents of anthocyanin metabolite cyanidin (CAS: 528-58-5) and cyanin (CAS: 20905-74-2)related to peel coloration had the relation of E3316＞19147＞19141. The combined analysis of transcriptome and metabolome bothshowed that the difference of E3316_vs_19141 was bigger than E3316_vs_19147. Under the regulation of epistatic genes,andgenes were expressed simultaneously in eggplant peel, but the expression level ofwas much higher than that of. The anthocyanin content of E3316 was higher than that of its parents, while the chlorogenic acid content was lower than that of its parents. 【Conclusion】Under the interaction of epistatic genes with controlling the peel color, the locus (type) and mutation mode of mutant gene in parent determined the trend of gene expression in the whole anthocyanin metabolism pathwayand the inhibition mode of peel pigmentation. The reddish-purple-peel F1hybrid of two white-peel parents was the result of the function complementation of two mutated gene involved in anthocyanin biosynthesis pathway.The highly expressedwas the main cause for synthesis of cyanidin, and the relationship between anthocyanin and chlorogenic acid biosynthesis was competitive in peels.

eggplant; peel color; epistatic genes; anthocyanins; transcriptome; metabolome

10.3864/j.issn.0578-1752.2022.20.011

2022-01-20；

2022-06-06

國(guó)家自然科學(xué)基金（31501755）、廣東省現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系（2021KJ110，2021KJ106）、廣東省自然科學(xué)基金（2021A1515012490）、廣東省農(nóng)業(yè)科學(xué)院農(nóng)業(yè)優(yōu)勢(shì)產(chǎn)業(yè)學(xué)科團(tuán)隊(duì)建設(shè)項(xiàng)目（202114TD，R2020PY-JX003）

孫保娟，Tel：020-38469456；E-mail：sunbaojuan@hotmail.com。通信作者李植良，E-mail：vri_li@163.com

（責(zé)任編輯 趙伶俐）