低溫離體保存黃獨微型塊莖轉錄組、蛋白質(zhì)組和代謝組的關聯(lián)分析

洪森榮，吳夏俊鵬，徐文慧，占學林，謝妮妮，蔣 妍，汪金華，凌 飛，吳麗霞，萬 琳

(上饒師范學院 生命科學學院，江西 上饒 334001)

低溫離體保存黃獨微型塊莖轉錄組、蛋白質(zhì)組和代謝組的關聯(lián)分析

洪森榮，吳夏俊鵬，徐文慧，占學林，謝妮妮，蔣 妍，汪金華，凌 飛，吳麗霞，萬 琳

(上饒師范學院 生命科學學院，江西 上饒 334001)

為探究黃獨微型塊莖低溫離體保存的內(nèi)在機理，對其轉錄組、蛋白質(zhì)組和代謝組進行關聯(lián)分析。結果表明：黃獨微型塊莖低溫和常溫離體保存的轉錄本比例和蛋白比例的log2值均呈現(xiàn)正態(tài)分布，十分類似，符合真實生物樣品的要求；差異蛋白和差異轉錄本韋恩圖結果表明，84個轉錄本在兩者中均有差異，表明大多數(shù)差異的轉錄本在蛋白組數(shù)據(jù)均有差異；差異蛋白和差異轉錄本的log2對數(shù)熱圖分析結果表明，相對于蛋白質(zhì)組來說，轉錄組差異程度更大，兩者的差異轉錄本多數(shù)一致，少數(shù)不同；差異蛋白和差異轉錄本log2散點關聯(lián)密度圖分析結果表明，兩者的比例在0附近居多，左下角和右上角(轉錄本和蛋白的上下調(diào)關系一致)相對左上角和右下角(轉錄和蛋白上下調(diào)關系不一致)更多，表明兩者差異表達趨勢一致性的程度更高。低溫離體保存黃獨微型塊莖轉錄組、蛋白質(zhì)組和代謝組的關聯(lián)分析結果表明，低溫離體保存的黃獨微型塊莖主要涉及碳代謝、氨基酸生物合成、糖酵解途徑、淀粉蔗糖代謝和丙酮酸代謝等途徑，這為黃獨微型塊莖的低溫離體保存和低溫破除休眠提供了理論依據(jù)。

黃獨；微型塊莖；低溫離體保存；轉錄組；蛋白質(zhì)組；代謝組

黃獨(DioscoreabulbiferaL.)為薯蕷科薯蕷屬多年生藤本植物[1-2]，其地下塊莖，俗名黃藥子，性平、味苦、有毒，有散結消癭、清熱解毒和涼血止血的功效，在臨床上黃藥子還具有抗病毒、抗腫瘤等活性，常用于治療甲狀腺疾病[3]。黃獨長期進行營養(yǎng)繁殖，造成病毒感染嚴重，產(chǎn)量逐年下降，品質(zhì)不斷退化[4]。因此，保存黃獨種質(zhì)資源具有一定的現(xiàn)實意義。黃獨的腋芽處常會長出地上變態(tài)塊莖，俗稱零余子、珠芽等，在組織培養(yǎng)中特稱它為微型塊莖[5]。微型塊莖由于體積小、攜帶方便、活力久、成活率高，常可作為薯蕷屬植物快速繁殖的材料[6]。研究表明，黃獨微型塊莖能夠在黑暗和低溫條件下長期保存而不喪失活力，還有助于打破休眠[7]，但微型塊莖低溫離體保存的內(nèi)在原因尚無報道。高通量組學分析工具的發(fā)展，引導系統(tǒng)生物學進入大數(shù)據(jù)時代，在單組學研究顯現(xiàn)瓶頸之際，可通過多組學整合關聯(lián)分析，深入闡明生命活動的本質(zhì)和規(guī)律[8]。本研究將對黃獨微型塊莖低溫離體保存的轉錄組、蛋白質(zhì)組和代謝組三大組學進行關聯(lián)分析，旨在為黃獨微型塊莖低溫離體保存提供理論依據(jù)。

1 材料與方法

1.1 材料

黃獨微型塊莖由上饒師范學院生命科學學院植物組織培養(yǎng)室提供。

1.2 方法

1.2.1 轉錄組、代謝組和蛋白質(zhì)組檢測

利用黃獨試管苗的帶芽莖段，接種到MS + KT 2 mg·L-1+ NAA 0.5 mg·L-1+蔗糖60 g·L-1的液體培養(yǎng)基中，使其誘導出微型塊莖，然后挑選出直徑約0.5 cm的微型塊莖轉移到100 mL無菌三角瓶(空瓶)中，4 ℃低溫離體保存。試驗材料分成4 ℃處理組(以Con4表示)和25 ℃對照組(以Con25表示)，保存60 d后進行轉錄組學、代謝組學和蛋白質(zhì)組學檢測。每個溫度3個重復。基本方法如下。

轉錄組學：提取2組試驗樣品的總RNA后，通過Bioanalyzer 2100(Agilent，Germany)完成質(zhì)量控制。然后進行mRNA捕獲和片段化處理，隨后用SuperScriptIII cDNA Synthesis Kit (Life technologies，USA)進行cDNA反轉錄。根據(jù)Illumina公司HiSeq SBS Kit and Cluster Kit v4(Illumina，Sandiego)說明制備Total RNA測序文庫。最后，按照Illumina公司Hiseq 2500的操作說明對cDNA文庫進行2×125 bp的高通量測序。

代謝組學：2組試驗樣品在液氮中磨碎，并轉移到10 mL離心管，經(jīng)過渦旋振蕩、超聲波清洗機處理、氮氣吹干等處理，測定其次生代謝物含量，GC/MS檢測設備為Agilent 7890A/5975C氣-質(zhì)聯(lián)用儀(安捷倫，美國)。利用GC-MS預處理軟件XCMS(www.bioconductor.org/)對Agilent 7890A/5975C氣質(zhì)聯(lián)用儀檢測獲得的原始文件進行數(shù)據(jù)預處理。最后對樣本進行信息比對分析，并對得到的數(shù)據(jù)進行生物信息學分析。

蛋白質(zhì)組學：提取2組試驗樣品的蛋白，加入trypsin進行酶解和TMT標記，標記后的樣本分段是用Agilent 300Extend C18 column(5 μm particles， 4.6 mm ID， 250 mm length)在高pH進行反向HPLC分離，檢測用Thermo ScientificTM Q ExactiveTM Plus液質(zhì)聯(lián)用儀配備Thermo ScientificTM EASY-nLC 1000TM納升級UPHLC分離系統(tǒng)進行分析。

1.2.2 轉錄組、蛋白質(zhì)組和代謝組的關聯(lián)分析

密度圖分析：密度分布圖用于考察轉錄組和蛋白組Con4/Con25差異倍性變化的分布，選用軟件為R語言(www.r-project.org)。

Venn圖分析：韋恩圖用于考察轉錄組和蛋白組差異蛋白之間的異同。在本研究中，由于蛋白分析所用的參考庫為轉錄組庫，我們以所檢測蛋白為基礎，分析差異轉錄本和蛋白之間的異同，所用軟件為venny2.0.2(http://bioinfogp.cnb.csic.es/tools/venny/)。

熱圖分析：對轉錄組和蛋白組Con4/Con25差異倍性變化進行直觀的熱圖分析，考察轉錄本和蛋白差異的變化。所用方法為層次聚類，聚類軟件為R語言(www.r-project.org)包Pheatmap。

散點密度圖分析：對轉錄組和蛋白組Con4/Con25差異倍性變化進行直觀的散點密度圖分析，考察轉錄本和蛋白差異的變化。選用軟件為R語言(www.r-project.org)。

代謝通路分析：對差異的蛋白和轉錄本及差異的代謝產(chǎn)物同時向KEGG通路進行映射，從代謝通路圖上了解整體變化。

2 結果與分析

2.1黃獨低溫保存微型塊莖轉錄組、蛋白組和代謝組的簡要分析

樣品間基因表達水平相關性是檢驗實驗可靠性和樣本選擇是否合理的重要指標，每個樣品有3個重復，相關系數(shù)為0.885，表明樣品之間表達模式的相似度較高。在黃獨微型塊莖低溫離體保存中，共獲得164 145個差異表達基因，其中63 305個基因表達上調(diào)，100 840個基因表達下調(diào)，部分結果見表1。差異表達基因富集度極顯著(P<0.01)的部分相關GO term有液泡繼承、單鏈斷裂修復、紡錘體伸長、麥芽糖分解代謝過程、甘露糖基轉移酶活性、甘露糖磷酸轉移酶活性、細胞壁甘露糖蛋白的生物合成過程、G1期細胞有絲分裂周期早期細胞芽和有絲分裂紡錘體定位的建立等。Con4、Con25兩組樣本的差異蛋白有106個，上調(diào)表達蛋白61個，下調(diào)表達蛋白45個，其中，伸長因子3， 6-磷酸葡萄糖酸脫氫酶、磷酸甘油酸激酶、蔗糖合成酶、分子伴侶DNAK和S-腺苷甲硫氨酸合成酶等為在4 ℃離體保存黃獨微型塊莖中上調(diào)表達的蛋白，而α-淀粉酶、淀粉磷酸化酶等為在4 ℃離體保存黃獨微型塊莖中下調(diào)表達的蛋白。與25 ℃離體保存的黃獨微型塊莖相比較，4 ℃離體保存的黃獨微型塊莖差異性代謝物有丙氨酸、兒茶素、N， N-雙(2-羥乙基)甲胺、水楊酸和山梨糖等。

2.2黃獨低溫保存微型塊莖密度圖比較

對Con4/Con25各差異轉錄本和差異蛋白的log2進行柱狀密度曲線分析，考察轉錄本比例和蛋白比例的分布。結果(圖1)表明，大多數(shù)轉錄本比例和蛋白比例的log2位于0附近，轉錄本比例和蛋白比例的log2均呈現(xiàn)正態(tài)分布，兩者比較類似，符合真實生物樣品的要求。

圖1 黃獨低溫保存微型塊莖差異轉錄本和蛋白比例密度曲線Fig.1 Differential transcript and protein density profile of Dioscorea bulbifera L. microtuber conserved at low temperature

2.3 Venn圖比較分析

以蛋白組獲得的蛋白為基礎，同時對Con4和Con25中的差異蛋白和差異轉錄本進行韋恩圖分析。由圖2可知，Con4和Con25中共有84個轉錄本在轉錄組和蛋白質(zhì)組中的表達均有差異，表明大多數(shù)差異轉錄本在蛋白質(zhì)組中均有差異。

表1黃獨低溫保存微型塊莖轉錄組、蛋白組部分差異表達結果

Table1Partial differential expression in transcriptome and proteome ofDioscoreabulbiferaL. microtuber conserved at low temperature

基因名稱Genename轉錄組Transcriptome上/下調(diào)Up/down表達量倍性變化的對數(shù)Logarithmofgeneploidyvaria-tion蛋白質(zhì)組Proteome蛋白倍性變化Variationofproteincontent蛋白量倍性變化的對數(shù)Logarithmofproteincon-tentploidyvariation上/下調(diào)Up/downHEXA_B下調(diào)Down-1.662970.718-0.47794下調(diào)DownRFWD2.COP1下調(diào)Down-3.55020.639-0.64611下調(diào)DownE1.11.1.7下調(diào)Down-1.440570.535-0.90239下調(diào)DownGLGP.PYG下調(diào)Down-7.051420.386-1.37333下調(diào)DownGST.GST下調(diào)Down-1.459430.738-0.43831下調(diào)DownE3.2.1.86B.BGLA下調(diào)Down未獲得Notapplicable0.644-0.63487下調(diào)DownAMYA.MALS下調(diào)Down-1.725830.349-1.5187下調(diào)DownPLD1_2下調(diào)Down-1.834780.745-0.42469下調(diào)DownGST.GST下調(diào)Down-4.19640.652-0.61706下調(diào)DownFDH下調(diào)Down-2.397250.591-0.75877下調(diào)DownRHM下調(diào)Down未獲得Notapplicable1.9950.996389上調(diào)UpMETE下調(diào)Down-1.975751.5260.609755上調(diào)UpMETK下調(diào)Down-1.142792.741.454176上調(diào)UpATPF1A.ATPA上調(diào)Up2.8073551.5570.638769上調(diào)UpMETK下調(diào)Down-1.556042.0131.009347上調(diào)UpGDHA上調(diào)Up1.1890342.3691.244278上調(diào)UpMDH2上調(diào)Up未獲得Notapplicable1.6190.695103上調(diào)UpMETE下調(diào)Down-2.121021.70.765535上調(diào)UpHPPA下調(diào)Down-2.860461.6770.745883上調(diào)UpPDXS.PDX1下調(diào)Down-1.591241.5580.639695上調(diào)UpISPH.LYTB下調(diào)Down-1.284531.5820.66175上調(diào)UpPGD.GND下調(diào)Down未獲得Notapplicable4.5312.179829上調(diào)UpPGK.PGK上調(diào)Up未獲得Notapplicable3.9971.998918上調(diào)UpGAPDH.GAPA上調(diào)Up未獲得Notapplicable1.030.042644無變化NochangeALDO下調(diào)Down-1.619520.846-0.24127無變化NochangeMETE下調(diào)Down-2.536051.3230.403813無變化NochangeENO.ENO上調(diào)Up3.946050.839-0.25326無變化NochangeE3.6.3.6上調(diào)Up3.8073551.2560.328836無變化NochangeATPEV1A.ATP6A上調(diào)Up未獲得Notapplicable1.0610.085425無變化NochangeLOX1_5下調(diào)Down-4.481130.904-0.14561無變化NochangeATPEV1A.ATP6A下調(diào)Down-1.510470.95-0.074無變化NochangeASNB.ASNS上調(diào)Up未獲得Notapplicable1.1640.219091無變化NochangeNOP1.FBL下調(diào)Down-5.321931.3770.461529無變化NochangeALDO下調(diào)Down-1.046540.885-0.17625無變化NochangeMAEB下調(diào)Down-1.259060.824-0.27928無變化NochangePK.PYK上調(diào)Up2.6130730.989-0.01596無變化NochangeMETE下調(diào)Down-2.562161.4590.54498無變化NochangeCHS上調(diào)Up3.4059921.0350.049631無變化NochangeE3.2.1.21下調(diào)Down-2.169930.782-0.35476無變化NochangeALDO下調(diào)Down-2.165340.909-0.13765無變化NochangeASD上調(diào)Up2.1203520.772-0.37333無變化Nochange

續(xù)表1

圖2 黃獨低溫保存微型塊莖差異轉錄本和蛋白維恩圖分析Fig.2 Venn diagram analysis of differential transcripts and proteins of Dioscorea bulbifera L. microtuber conserved at low temperature

2.4 熱圖比較分析

以蛋白組獲得的蛋白為基礎，對差異蛋白和差異轉錄本的Con4/Con25各蛋白比例的log2進行熱圖分析，比較兩者之間的異同。從圖3可知，相對于蛋白組來說，轉錄組差異程度更大，兩者差異蛋白多數(shù)一致，少數(shù)不同。兩者均表現(xiàn)差異的基因主要有氨基己糖苷酶、過氧化物酶、淀粉磷酸化酶、谷胱甘肽S-轉移酶、6-磷酸-β-葡糖苷酶、α-淀粉酶 、高半胱氨酸甲基轉移酶、S-腺苷甲硫氨酸合成酶、F型H +轉運ATP酶α亞基、谷氨酸脫氫酶(NADP+)、蘋果酸脫氫酶酶、6-磷酸葡萄糖酸脫氫酶、磷酸甘油酸激酶等。

聚類軟件為R語言(www.r-project.org)包Pheatmap。紅色是表達顯著上調(diào)的蛋白，綠色表示表達顯著下調(diào)的蛋白The clustering software was R language pack (www.r-project.org) Pheatmap. Red meant significantly up-regulated proteins， and green meant significantly down-regulated proteins圖3 黃獨低溫保存微型塊莖差異轉錄本和蛋白比例熱圖Fig.3 Thermography of differential transcripts and proteins ratio of Dioscorea bulbifera L. microtuber conserved at low temperature

2.5 散點關聯(lián)密度圖比較

以蛋白組獲得的蛋白為基礎，同時對差異蛋白和差異轉錄本中各蛋白Con4/Con25比例的log2值進行散點關聯(lián)密度圖分析，比較兩者的異同。從圖4可知，兩者的比例在0附近居多，左下角和右上角(轉錄本和蛋白的上下調(diào)關系一致)相對左上角和右下角(轉錄本和蛋白的上下調(diào)關系不一致)更多，表明兩者差異表達趨勢一致性的程度更高。

圖4 黃獨低溫保存微型塊莖差異轉錄本和蛋白比例散點密度圖分析Fig.4 Density map analysis of differential transcripts and proteins of Dioscorea bulbifera L. microtuber conserved at low temperature

2.6 代謝通路綜合分析

在生物體內(nèi)，不同蛋白相互協(xié)調(diào)行使其生物學功能，通過通路顯著性富集能確定差異表達蛋白參與的最主要生化代謝途徑和信號轉導途徑，KEGG(Kyoto Encyclopedia of Genes and Genomes)是有關通路的主要公共數(shù)據(jù)庫。與常溫離體保存相比，低溫離體保存的黃獨微型塊莖轉錄組分析表明，樣本差異基因富集的KEGG途徑有鈣信號通路、檸檬酸循環(huán)(TCA循環(huán))、氧化磷酸化、黃酮和黃酮醇的生物合成等(表2)。蛋白質(zhì)組分析表明，黃獨低溫保存微型塊莖部分差異蛋白富集的GO term有谷胱甘肽轉移酶活性、高爾基液泡運輸和披網(wǎng)格蛋白小泡膜等(表3)。代謝組分析表明，差異性代謝物有丙氨酸、兒茶素和水楊酸等(表4)。其中，丙氨酸參與氰基氨基酸代謝，兒茶素參與次生代謝產(chǎn)物生物合成、黃酮類化合物的生物合成和苯丙素的生物合成；水楊酸參與多環(huán)芳烴降解、微生物在不同環(huán)境中的代謝，植物激素信號轉導，次生代謝產(chǎn)物生物合成，二惡英降解，苯丙氨酸代謝，芳烴降解，植物激素生物合成，鐵載體組非核糖體肽合成和苯丙素的生物合成等。

表2黃獨低溫保存微型塊莖部分差異基因富集的KEGG途徑

Table2KEGG pathways of part differentially expressed genes ofDioscoreabulbiferaL. microtuber conserved at low temperature

途徑PathwayKOKEGG注釋的轉錄本數(shù)量NumbersoftranscriptswithKEGGannotationKEGG注釋的DELs數(shù)量NumbersofDETswithKEGGannotationP值P-value校正P值CorrectedP-value鈣信號通路Calciumsignalingpathwayko0402064580.0001420.000589檸檬酸循環(huán)Citratecycleko000206084720.0054630.019306氧化磷酸化Oxidativephosphorylationko001904523530.0074930.025607黃酮和黃酮醇的生物合成Flavoneandflavonolbiosynthesisko0094477680.0003710.001500

表3黃獨低溫保存微型塊莖部分差異蛋白富集的GO terms

Table3Enriched GO terms of part differentially proteins ofDioscoreabulbiferaL. microtuber conserved at low temperature

類別CategoryP值P-value差異蛋白數(shù)量Numberofdifferentialproteins蛋白數(shù)量Numberofproteins功能描述信息Functiondescriptioninformation功能信息FunctionalinformationGO:00043640.01304748434谷胱甘肽轉移酶活性Glutathionetransferaseactivity分子功能MolecularfunctionGO:00068960.02403472422高爾基液泡運輸Golgitovacuoletransport生物學功能BiologicalprocessGO:00306650.02403472422披網(wǎng)格蛋白小泡膜Clathrin-Coatedvesiclemembrane細胞組件Cellularcomponent

表4黃獨低溫保存微型塊莖部分差異代謝物

Table4Differential metabolites inDioscoreabulbiferaL. microtuber conserved at low temperature

ID名稱NamesP值PvalueLog(Con4/Con25)VIP值VIPvalue34Alanine0.038054.1157861.9936197Catechin0.003892-1.635583.44036171Salicylicacid0.001429-2.139291.49439

VIP值，對模型貢獻度大(超過1)的化合物的值，一般根據(jù)PLS-DA模型的VIP值大于1、P<0.05的原則獲得差異性物質(zhì)。

The VIP value refered to the value of the compounds which were more than 1 and contributed greatly to the model. In generally， differential substances were obtained according to the principle that the VIP value of the PLS-DA model was greater than 1 and thePvalue was less than 0.05.

3 討論

系統(tǒng)生物學的發(fā)展產(chǎn)生了海量的基因組、轉錄組、蛋白質(zhì)組和代謝組等組學數(shù)據(jù)。基于多平臺組學數(shù)據(jù)整合的代謝通路分析為全面深刻了解生物系統(tǒng)提供了前提條件[9]。近年來，系統(tǒng)生物學已逐漸成為后基因組學時代生命科學研究的新熱點。轉錄組學、蛋白質(zhì)組學、代謝組學技術是目前系統(tǒng)生物學研究的重要方法。通過轉錄組學、蛋白質(zhì)組學和代謝組學數(shù)據(jù)關聯(lián)分析，可以深入理解基因、蛋白質(zhì)和代謝物之間的內(nèi)在聯(lián)系[10]。張振亞等[11]對正常培養(yǎng)和鹽脅迫條件下的2個南方型紫花苜蓿的根系進行轉錄組和蛋白質(zhì)組關聯(lián)分析，發(fā)現(xiàn)了紫花苜蓿耐鹽相關候選基因。馬進等[8]通過蛋白質(zhì)組和轉錄組關聯(lián)分析，揭示了南方型紫花苜蓿適應鹽脅迫環(huán)境的分子機制。張黛靜等[12]通過對銅脅迫下小麥幼根轉錄組學及蛋白質(zhì)組學的研究，發(fā)現(xiàn)小麥在重金屬逆境中轉錄水平和翻譯水平具有一致性，差異均表現(xiàn)在代謝、物質(zhì)合成、糖酵解等生命活動方面。本研究通過黃獨微型塊莖低溫離體保存轉錄組、蛋白質(zhì)組和代謝組的關聯(lián)分析，發(fā)現(xiàn)轉錄本和蛋白的比例均呈現(xiàn)正態(tài)分布，兩者比較類似，符合真實生物樣品的要求；黃獨低溫離體保存微型塊莖差異基因富集的KEGG途徑主要有鈣信號通路、檸檬酸循環(huán)、氧化磷酸化、黃酮和黃酮醇的生物合成等，差異蛋白富集的GO term主要有谷胱甘肽轉移酶活性、高爾基液泡運輸和披網(wǎng)格蛋白小泡膜等；差異代謝物主要有丙氨酸、兒茶素和水楊酸等。黃獨微型塊莖低溫離體保存的主要差異基因有過氧化物酶、谷胱甘肽S-轉移酶、淀粉磷酸化酶、S-腺苷甲硫氨酸合成酶、谷氨酸脫氫酶(NADP+)、蘋果酸脫氫酶、6-磷酸葡萄糖酸脫氫酶、磷酸甘油酸激酶等，這些基因大多與檸檬酸循環(huán)、氧化磷酸化相關。究其原因，推測與黃獨微型塊莖低溫保存需要較多能量供應有關。可見，黃獨微型塊莖低溫離體保存的轉錄水平、翻譯水平以及代謝水平具有一致性。

基因表達的主要環(huán)節(jié)包括轉錄、蛋白質(zhì)合成以及代謝物的形成。因此，從生長條件和狀態(tài)相同的細胞、組織或器官獲得的轉錄組、蛋白質(zhì)組和代謝物數(shù)據(jù)之間應該具有較高的相關性[13]。Koh等[14]比較了干旱脅迫下油菜的20個基因在轉錄和蛋白質(zhì)水平的表達，發(fā)現(xiàn)盡管有些基因轉錄本和蛋白的表達模式不盡相同，但二者之間仍存在顯著正相關。Budak等[15]發(fā)現(xiàn)，硬粒小麥品種Kiziltan和野生二粒小麥系TR39477、TTD22(T)在有些脅迫條件下，轉錄本的變化和蛋白水平的變化一致，但也有部分轉錄本的變化和蛋白水平的變化趨勢相反。本研究結果與這2項研究結果相一致。本研究中，黃獨微型塊莖低溫離體保存的差異蛋白和差異轉錄本韋恩圖分析結果表明，84個轉錄本在兩者中均為差異，表明大多數(shù)差異的轉錄本在蛋白組數(shù)據(jù)均為差異。差異蛋白和差異轉錄本的log2對數(shù)熱圖分析結果也表明，相對于蛋白組來說，轉錄組差異的程度更大，蛋白組數(shù)據(jù)差異的程度較小。兩者差異蛋白多數(shù)一致，少數(shù)有不同。差異蛋白和差異轉錄本log2散點關聯(lián)密度圖分析結果表明，兩者的比例在0附近居多，左下角和右上角(轉錄和蛋白上下調(diào)關系一致)相對左上角和右下角(轉錄和蛋白上下調(diào)關系不一致)更多，表明兩者差異表達趨勢一致性程度更高。但Peng等[16]在模擬干旱脅迫下，對體細胞雜交小麥品種Shanrong 3及其親本面包小麥品種Jinan 177苗期葉片和根進行蛋白質(zhì)組學和轉錄組學分析，發(fā)現(xiàn)僅有20個(27.0%)差異表達蛋白在轉錄水平與蛋白水平有相關性，指出其蛋白質(zhì)組和轉錄組之間的相關性不強。黃獨微型塊莖低溫離體保存的轉錄組、蛋白質(zhì)組和代謝組的關聯(lián)分析，可以為黃獨微型塊莖低溫離體保存提供理論依據(jù)。

[1] PASSAM H C， WICKHAM L D， WILSON L A. Comparative observations on the polarity of sprouting of bulbils ofDioscoreabulbiferaL. andDioscoreaalataL.[J].AnnalsofBotany， 1982， 49(3): 359-366.

[2] TEPONNO R B， PONOU B K， TAPONDJOU L A， et al. Bafoudiosbulbin H， a new clerodane diterpene from the flowers ofDioscoreabulbiferaL. varsativa[J].PhytochemistryLetters， 2013， 6(2): 310-314.

[3] 劉勁松， 高衛(wèi)娜， 鄭娟， 等. 黃獨鮮塊根化學成分研究[J]. 中國中藥雜志， 2017， 42(3): 510-516.

LIU J S， GAO W N， ZHENG J， et al. Chemical constituents from fresh tubers ofDioscoreabulbifera[J].ChinaJournalofChineseMateriaMedica， 2017， 42(3): 510-516. (in Chinese with English abstract)

[4] 尹明華， 洪森榮. 黃獨的莖尖培養(yǎng)和病毒檢測[J]. 中草藥， 2009， 40(9): 1462-1466.

YIN M H， HONG S R. Shoot tip culture and virus detection ofDioscoreabulbifera[J].ChineseTraditionalandHerbalDrugs， 2009， 40(9): 1462-1466. (in Chinese with English abstract)

[5] 王運英， 張曉麗， 白英豪， 等. 山藥微型塊莖萌發(fā)影響因素研究[J]. 北方園藝， 2015 (17): 194-196.

WANG Y Y， ZHANG X L， BAI Y H， et al. Impact factors in germination of microtubers fromDioscoreaopposite[J].NorthernHorticulture， 2015 (17): 194-196. (in Chinese with English abstract)

[6] 李明軍， 鄧麗， 劉欣英， 等. 生長素和細胞分裂素對懷山藥微型塊莖誘導形成的影響[J]. 河南農(nóng)業(yè)科學， 2008， 37(11): 102-106.

LI M J， DENG L， LIU X Y， et al. Effects of auxin and cytokinin on the induction of microtuber ofDisocoreaoppositaThunb.[J].JournalofHenanAgriculturalSciences， 2008， 37(11): 102-106. (in Chinese with English abstract)

[7] 尹明華， 洪森榮， 夏瑾華， 等. 黃獨微型塊莖低溫保存及其萌發(fā)苗遺傳穩(wěn)定性研究[J]. 廣西植物， 2015 (5): 733-740.

YIN M H， HONG S R， XIA J H， et al. Low temperature conservation ofDioscoreabulbiferamicrotuber and genetic stability of its germination seedling[J].Guihaia， 2015 (5): 733-740. (in Chinese with English abstract)

[8] 馬進， 鄭鋼， 裴翠明， 等. 南方型紫花苜蓿葉片響應鹽脅迫蛋白質(zhì)組和轉錄組關聯(lián)分析[J]. 核農(nóng)學報， 2016， 30(9): 1706-1715.

MA J， ZHENG G， PEI C M， et al. Correlation analysis on proteomic and transcriptomic of salt-response in leaves of southern typeAlfalfa[J].JournalofNuclearAgriculturalSciences， 2016， 30(9): 1706-1715. (in Chinese with English abstract)

[9] 周滔， 李靜宜， 馬毅， 等. 基于組學數(shù)據(jù)庫整合工具的代謝通路分析應用[J]. 國際藥學研究雜志， 2015， 42(5): 587-592.

ZHOU T， LI J Y， MA Y， et al. Metabolic pathway analysis and application based on the database and integrated tool of omics[J].JournalofInternationalPharmaceuticalResearch， 2015， 42(5): 587-592. (in Chinese with English abstract)

[10] WU H X， JIA H M， MA X W， et al. Transcriptome and proteomic analysis of mango (MangiferaindicaLinn) fruits[J].JournalofProteomics， 2014， 105: 19-30.

[11] 張振亞， 裴翠明， 馬進. 基于轉錄組和蛋白質(zhì)組關聯(lián)研究技術篩選紫花苜蓿耐鹽相關候選基因[J]. 植物生理學報， 2016， 52(3): 317-324.

ZHANG Z Y， PEI C M， MA J. Screening of candidate salt tolerance-related genes inalfalfabased on transcriptome-proteome correlation research techniques[J].PlantPhysiologyJournal， 2016， 52(3): 317-324. (in Chinese with English abstract)

[12] 張黛靜， 王多多， 董文，等. 銅脅迫下小麥幼根轉錄組學及蛋白質(zhì)組學研究[J]. 河南農(nóng)業(yè)科學， 2015， 44(4): 31-35.

ZHANG D J， WANG D D， DONG W， et al. Transcriptomics and proteomics analysis in root of wheat under copper stress[J].JournalofHenanAgriculturalSciences， 2015， 44(4): 31-35. (in Chinese with English abstract)

[13] 張紅亮， 王道文， 張正斌. 利用轉錄組學和蛋白質(zhì)組學技術揭示小麥抗旱分子機制的研究進展[J]. 麥類作物學報， 2016， 36(7): 878-887.

ZHANG H L， WANG D W， ZHANG Z B. Progress and prospects in the research on wheat drought response and resistance mechanisms using transcriptomic and proteomic approaches[J].JournalofTriticeaeCrops， 2016， 36(7): 878-887. (in Chinese with English abstract)

[14] KOH J， CHEN G， YOO M J， et al. Comparative proteomic analysis ofBrassicanapusin response to drought stress[J].JournalofProteomeResearch， 2015， 14(8): 3068-3081.

[15] BUDAK H， AKPINAR B A， UNVER T， et al. Proteome changes in wild and modern wheat leaves upon drought stress by two-dimensional electrophoresis and nanoLC-ESI-MS/MS[J].PlantMolecularBiology， 2013， 83(1): 89-103.

[16] PENG Z， WANG M， LI F， et al. A proteomic study of the response to salinity and drought stress in an introgression strain of bread wheat[J].Molecular&CellularProteomics， 2009， 8(12): 2676-2686.

(責任編輯侯春曉)

Correlationanalysisoftranscriptome，proteomeandmetabolomeofDioscoreabulbiferaL.microtubersconservedinvitroatlowtemperature

HONG Senrong， WU-XIA Junpeng， XU Wenhui， ZHAN Xuelin， XIE Ni’ni， JIANG Yan， WANG Jinhua， LING Fei， WU Lixia， WAN Lin

(CollegeofLifeSciences，ShangraoNormalUniversity，Shangrao334001，China)

In order to explore the internal mechanism ofDioscoreabulbiferaL. microtubers conservedinvitroat low temperature， the correlation of transcriptome， proteome and metabolome was analyzed. The results showed that Log2of transcript ratio and protein ratio ofD.bulbiferaL. microtubers conservedinvitroat low temperature and room temperature all showed normal distribution， which were very similar and consistent with the requirements of real biological samples; Wayne plot analysis of differential proteins and differential transcripts showed that the 84 transcripts were different in both groups， indicating that most of the transcripts were differentially expressed in proteome data; Log2logarithm thermal analysis of differential proteins and differential transcripts showed that the difference in transcriptome was greater than that in proteome， most of their differential transcripts were consistent， few were different; Analysis of log2logarithmic scatter correlation density plot showed that ratio of differential proteins and transcripts was in the vicinity of 0， the lower left corner and upper right corner that indicating the up/down regulation relationship of transcripts and proteins was consistent were more compared with the upper left corner and lower right corner that indicating the up/down regulation relationship of transcripts and proteins was not consistent， which indicated that consistency degree of the trend of differential expression was higher. Correlation analysis results of transcriptome， proteome and metabolome showed that conservationinvitroofD.bulbiferaL. microtubers at low temperature mainly involved the pathway， such as carbon metabolism， amino acid biosynthesis， glycolysis， starch and sucrose metabolism， and pyruvate metabolism etc.， which could provide a theoretical basis for conservationinvitroat low temperature and low temperature breaking dormancy forD.bulbiferaL. microtubers.

DioscoreabulbiferaL.; microtuber; conservationinvitroat low-temperature; transcriptome; proteome; metabolome

洪森榮，吳夏俊鵬，徐文慧， 等. 低溫離體保存黃獨微型塊莖轉錄組、蛋白質(zhì)組和代謝組的關聯(lián)分析[J]. 浙江農(nóng)業(yè)學報， 2017， 29(11): 1827-1834.

10.3969/j.issn.1004-1524.2017.11.08

2017-04-05

國家自然科學基金項目(31360072)

洪森榮(1974—)，男，江西永新人，碩士，教授，研究方向為藥用植物生物技術。E-mail: hongsenrong@163. com

S632.1

A

1004-1524(2017)11-1827-08