劉含,郭銀萍,穆興燕,劉曉燕,2*,岑順友
1(貴陽(yáng)學(xué)院 食品與制藥工程學(xué)院,貴州 貴陽(yáng),550005)2(貴州省果品加工工程技術(shù)研究中心,貴州 貴陽(yáng),550005) 3(貴州宏財(cái)聚農(nóng)投資有限責(zé)任公司,貴州 六盤水,553535)
刺梨(RosaroxburghiiTratt)系薔薇科薔薇屬藥食同源植物[1],主要分布在我國(guó)西南地區(qū)[2]。因其獨(dú)特的風(fēng)味而備受關(guān)注,并且富含多糖、維生素、黃酮、多酚等具有藥理功效的成分[3],使得刺梨果實(shí)在醫(yī)藥領(lǐng)域擁有較高的保健和藥用價(jià)值[4]。在藥理方面的研究顯示,刺梨具有抗氧化、抗衰老、降血糖[5]、治療胃潰瘍、促進(jìn)消化、抗動(dòng)脈粥樣硬化等多種藥理作用[6]。
國(guó)內(nèi)對(duì)刺梨的研究已經(jīng)十分火熱,但也局限于傳統(tǒng)方法的研究,如營(yíng)養(yǎng)特性及產(chǎn)品開(kāi)發(fā),抗氧化、降血糖、降血脂等。針對(duì)刺梨活性成分的研究大多是從多酚、多糖等單一化合物的角度進(jìn)行,且刺梨在代謝組學(xué)方面的研究鮮有報(bào)道。同一物質(zhì)用不同溶劑提取,其含量、特性都可能發(fā)生改變,非靶向代謝組學(xué)能系統(tǒng)、全面地分析,獲取大量代謝物的數(shù)據(jù),從而解析提取物中小分子的差異。代謝物提取是代謝組學(xué)研究中最關(guān)鍵的一步,選擇合適的提取方法為后續(xù)利用代謝組學(xué)尋找刺梨具有功能活性的關(guān)鍵代謝產(chǎn)物提供了有效的方法,有利于對(duì)刺梨中各代謝物進(jìn)行深入研究。
代謝組學(xué)(metabolomics)是在特定條件下對(duì)生物體所含代謝物整體進(jìn)行定性定量,可揭示特定時(shí)間和特定條件下代謝物的種類和含量,覆蓋范圍廣[7],對(duì)于20多萬(wàn)種的植物代謝物[8],代謝組學(xué)的思路對(duì)其具有較大的意義。依據(jù)檢測(cè)方式與研究目的的不同,代謝組學(xué)被分為非靶向和靶向[9]。非靶向旨在檢測(cè)出樣品中所有的代謝物,能發(fā)現(xiàn)一些新的化合物。由于數(shù)據(jù)庫(kù)數(shù)據(jù)不全,定性時(shí)較為困難,也只能相對(duì)定量[10]。但它的優(yōu)勢(shì)在于整體性的檢測(cè)出樣品中的化合物,以便于后續(xù)鑒定分析。靶向代謝組學(xué)能較為精確的定量出目標(biāo)代謝物[11]。本文采用非靶向代謝組學(xué)對(duì)刺梨不同提取物進(jìn)行分析,通過(guò)整體檢測(cè)刺梨中的化合物,結(jié)合多元統(tǒng)計(jì)方法篩選出目標(biāo)代謝物,對(duì)代謝物進(jìn)行分析,查閱文獻(xiàn)發(fā)現(xiàn)其具有的功能活性,為進(jìn)一步研究提供數(shù)據(jù)支持。
刺梨,貴州省龍里縣;甲醇、乙腈(純度≥99.9%),美國(guó)Thermo公司;2-氯苯丙氨酸(純度98.5%),上海阿拉丁生化科技股份有限公司;甲酸(純度≥ 98%),東京化成工業(yè)株式會(huì)社;甲酸銨(純度≥ 99%),美國(guó)Sigma-Aldrich公司;H2O(純度100%),美國(guó)Millipore公司。
Vanquish液相色譜儀、Q-Exactive質(zhì)譜儀,美國(guó)Thermo公司;ACQUITY UPLC?HSS T3 1.8 μm (2.1 mm×150 mm)色譜柱,美國(guó)Waters公司。
1.2.1 提取物制備
刺梨全果提取物水-乙醇(占比60%)溶液的制備:稱取刺梨樣品按料液比1∶10(g∶mL)與混合,45 ℃水浴2 h,45 ℃、功率200 W超聲30 min,過(guò)濾,收集濾渣,重復(fù)提取2次合并濾液,蒸發(fā)濃縮至浸膏狀,分別用水和60%乙醇定容得到25 mL刺梨粗提物,置于-20 ℃冰箱保存?zhèn)溆?6個(gè)生物學(xué)重復(fù)樣本)。
刺梨果渣(由同批次刺梨果壓榨干燥得到)提取物的處理同1.2.1節(jié)。
1.2.2 樣本預(yù)處理
精確稱量樣本200 mg(±1%)于2 mL EP管中,準(zhǔn)確加入0.6 mL 2-氯苯丙氨酸(4 mg/L)甲醇(-20 ℃)配制,渦旋振蕩30 s加入100 mg玻璃珠,放入組織研磨器中,55 Hz研磨60 s室溫超聲15 min;12 000 r/min,4 ℃離心10 min,取上清液300 μL 過(guò)0.22 μm膜過(guò)濾,過(guò)濾液加入到檢測(cè)瓶中用剩余待測(cè)樣本進(jìn)行超高效液相色譜-高分辨質(zhì)譜技術(shù)(ultra high performance liquid chromatography high resolution mass spectrometry, UPLC-HRMS)檢測(cè)。
1.2.3 色譜-質(zhì)譜分析
色譜條件:采用ACQUITY UPLC?HSS T3 1.8 μm (2.1 mm×150 mm)色譜柱,自動(dòng)進(jìn)樣器溫度設(shè)為8 ℃,以0.25 mL/min的流速,40 ℃的柱溫,進(jìn)樣2 μL進(jìn)行梯度洗脫,流動(dòng)相為正離子0.1%(體積分?jǐn)?shù))甲酸水(C)-0.1%甲酸乙腈(D);負(fù)離子5 mmol/L甲酸銨水(A)-乙腈(B)。梯度洗脫程序?yàn)?~1 min,2%B/D;1~9 min,2%~50% B/D;9~12 min,50%~98% B/D;12~13.5 min,98% B/D;13.5~14 min,98%~2% B/D;14~20 min,2% D-正模式(14~17 min,2% B-負(fù)模式)。
質(zhì)譜條件:儀器使用電噴霧離子源(electron spray ionization,ESI),正負(fù)離子電離模式,正離子噴霧電壓為3.50 kV,負(fù)離子噴霧電壓為2.50 kV,鞘氣30 arb,輔助氣10 arb。毛細(xì)管溫度325 ℃,以分辨率70 000進(jìn)行全掃描,掃描范圍81~1 000,并采用HCD進(jìn)行二級(jí)裂解,碰撞電壓為30 eV,同時(shí)采用動(dòng)態(tài)排除去除無(wú)必要的MS/MS信息[12]。
1.2.4 數(shù)據(jù)預(yù)處理
通過(guò)Proteowizard軟件對(duì)原始數(shù)據(jù)進(jìn)行分析及處理,利用R語(yǔ)言的XCMS程序包進(jìn)行峰識(shí)別、峰過(guò)濾、峰對(duì)齊;得到包括質(zhì)荷比、保留時(shí)間及峰面積等信息的數(shù)據(jù)矩陣。
1.2.5 物質(zhì)鑒定
代謝物的鑒定首先根據(jù)精確分子質(zhì)量進(jìn)行確認(rèn),后續(xù)根據(jù)MS/MS碎片模式對(duì)HMDB、METLIN、Massbank、LipidMaps以及mzClound標(biāo)準(zhǔn)品數(shù)據(jù)庫(kù)確認(rèn)注釋獲得代謝物。
利用SIMCA 14.1軟件進(jìn)行主成分分析(principal component analysis, PCA),可以判斷不同樣本生物學(xué)重復(fù)之間代謝組的數(shù)據(jù)可靠性和穩(wěn)定性[13]。正交偏最小二乘法判別分析(orthogonal partial least squares discriminant analysis, OPLS-DA)模型表明各提取物組與其他3個(gè)提取物組之間的差異性[14]。對(duì)差異代謝物的篩選依據(jù)是SIMCA 14.1軟件導(dǎo)出的變量重要性投影(variable importance projection value,VIP)值、P值和差異倍數(shù)(fold change,FC)。VIP值>1,P值<0.05及FC≥2和FC≤0.5為一般的差異代謝物篩選標(biāo)準(zhǔn)。在此基礎(chǔ)上,本實(shí)驗(yàn)設(shè)定VIP值>1.3、P值<0.01、FC≥3和FC≤0.1篩選最終的顯著差異代謝物。通過(guò)MBRole 2.0代謝通路分析功能,根據(jù)篩選出來(lái)的差異代謝物KEGG ID獲得刺梨不同提取物代謝物的代謝途徑[15]。
UPLC-HRMS有較高的靈敏性和較廣的分析范圍,樣本中含量極低的物質(zhì)也能檢測(cè)出,因此非靶向代謝組學(xué)常用UPLC-HRMS檢測(cè)樣品。在刺梨的提取實(shí)驗(yàn)中,常用的方法是水提、醇提。刺梨工業(yè)中,刺梨果渣產(chǎn)量巨大且利用率不高,但根據(jù)已有研究發(fā)現(xiàn)果渣中的活性成分十分豐富,研究前景可觀。故設(shè)計(jì)以下4種提取方法:刺梨全果水提、刺梨全果醇提、刺梨果渣水提、刺梨果渣醇提,這4種不同提取物的正負(fù)離子模式下重疊的總離子流圖(total ion chromatogram,TIC)如圖1所示。
ESI+模式獲得22 447個(gè)前體分子,ESI-模式獲得11 835個(gè)前體分子。在二級(jí)質(zhì)譜信息中,4種提取方式共得到316種代謝物。主要包括脂類、氨基酸類、糖類、酚類、維生素等,其中脂類68種、氨基酸52種、糖類48種、酚類28種、維生素14種。
A-ESI+;B-ESI-圖1 刺梨不同提取物總離子流圖Fig.1 Total ion flow diagram of different extracts of Rosa roxburghii Tratt
PCA結(jié)果(圖2)顯示,刺梨不同提取物的生物學(xué)重復(fù)都匯集在一起,且全部樣本均位于95%置信區(qū)間內(nèi)[16],表明不同提取物數(shù)據(jù)的可靠性與穩(wěn)定性較好,在測(cè)定過(guò)程中未出現(xiàn)明顯偏差[17]。每個(gè)圖形表示一個(gè)樣本,不同的組根據(jù)圖形的顏色和形狀來(lái)區(qū)別。圖形之間越近,樣品中代謝物越相似;反之,樣本越遠(yuǎn),其整體代謝物差異越大。通過(guò)觀察全部樣本的PCA得分圖,可以體現(xiàn)樣本的總體分布趨勢(shì)。PCA顯示出4個(gè)主成分,其中主成分1(PC1)的貢獻(xiàn)率為45.8%,主成分2(PC2)的貢獻(xiàn)率為23.2%,不同樣品組在兩個(gè)坐標(biāo)軸上彼此分離,說(shuō)明刺梨不同提取物的代謝物有一定的差異。
圖2 刺梨不同提取物代謝組的PCA得分Fig.2 PCA scores of different extracts of Rosa roxburghii Tratt
為了提高分類效能,在PCA模型初步處理樣本后,采用能降低系統(tǒng)噪聲干擾的OPLS-DA模型分析數(shù)據(jù)[18]。95%置信區(qū)間是衡量樣本可取性的重要指標(biāo)[19],本次實(shí)驗(yàn)的全部樣本處于區(qū)間內(nèi),表明本次樣本具有可取性。不同提取物OPLS-DA分析結(jié)果表明,每個(gè)提取物與其他3個(gè)提取物在第一主成分軸上明顯分開(kāi),說(shuō)明其代謝物之間代謝成分差異明顯(圖3)。在圖3四幅OPLS-DA得分圖中樣本間橫向距離越遠(yuǎn)說(shuō)明組間差異越大,縱向距離越近說(shuō)明組內(nèi)重復(fù)性越好。圖中數(shù)據(jù)表明,2組樣本存在明顯差異性。
置換檢驗(yàn)?zāi)軌蚺袛郞PLS-DA模型是否過(guò)擬合,對(duì)于OPLS-DA模型導(dǎo)出的VIP值、P值和FC才更具說(shuō)服力[20]。圖中橫坐標(biāo)為保留原始數(shù)據(jù)的百分比,縱坐標(biāo)為R2Y或Q2值[21]。最右邊橫縱坐標(biāo)均為1的點(diǎn)為真實(shí)值,左邊所有的模擬預(yù)測(cè)的R2(>0.9)和Q2(>0.9)都要低于真實(shí)值[22],并且R2總是大于Q2,Q2的回歸線截距要小于0.05。對(duì)OPLS-DA模型進(jìn)行200次排列實(shí)驗(yàn),得到模型評(píng)估指標(biāo)見(jiàn)表1。本次實(shí)驗(yàn)數(shù)據(jù)均滿足以上條件,表明模型可靠沒(méi)有過(guò)擬合[23]。
表1 OPLS-DA模型的主要評(píng)估指標(biāo)Table 1 Main evaluation indexes of OPLS-DA model
對(duì)316個(gè)物質(zhì)進(jìn)行差異代謝物篩選,根據(jù)VIP值>1、P值<0.05、FC≥2和FC≤0.5共篩選出168個(gè)差異代謝物。各提取組相較于其他3組分別篩選到30、42、75、21個(gè)差異代謝物。在此基礎(chǔ)上,對(duì)不同樣本組之間具有差異的代謝物進(jìn)行更嚴(yán)格的篩選,最終得到49種顯著差異代謝物(VIP>1.3、P<0.01、FC≥3和FC≤0.1)(表2)。
如表2所示,各提取組相較于其他3組篩選出的49種差異代謝物為5種、17種、21種、6種。其中脂類11種、維生素2種、糖類10種、酚類10種、氨基酸6種、核苷酸4種、醇2種、萜類4種。刺梨中的許多代謝產(chǎn)物,如萜類、脂類、酚類、多糖等,對(duì)人體具有促進(jìn)健康的功能。
A-全果水提與其他組對(duì)比;B-全果水提置換檢驗(yàn);C-全果醇提與其他組對(duì)比;D-全果醇提置換檢驗(yàn);E-果渣水提與其他組對(duì)比; F-果渣水提置換檢驗(yàn);G-果渣醇提與其他組對(duì)比;H-果渣醇提置換檢驗(yàn)圖3 刺梨不同提取物代謝物的OPLS-DA得分和置換檢驗(yàn)Fig.3 OPLS-DA score and displacement test of metabolites of different extracts of Rosa roxburghii Tratt
通過(guò)火山圖可以直觀地觀察到各組分與其他組分的篩選出的顯著差異代謝物(圖4),圖中橫坐標(biāo)表示log2FC縱坐標(biāo)表示-log10P圓點(diǎn)為不顯著代謝物,菱形為log2FC>0的代謝物,六邊形為log2FC<0的代謝物,篩選出的49種差異代謝物用注釋編號(hào)標(biāo)出。其中全果干粉水提上調(diào)代謝物有2個(gè),下調(diào)代謝物3個(gè)。全果干粉醇提上調(diào)代謝物有7個(gè),下調(diào)代謝物有10個(gè)。果渣干粉水提上調(diào)代謝物有19個(gè),下調(diào)代謝物有2個(gè)。果渣干粉醇提上調(diào)代謝物有1個(gè),下調(diào)代謝物有5個(gè)??傮w而言,上調(diào)差異代謝物多于下調(diào)差異代謝物。在本次實(shí)驗(yàn)中,上調(diào)差異代謝物黃酮類為最多。有研究發(fā)現(xiàn)在一定的條件下脂肪酸會(huì)降解為部分黃酮類化合物,這是因?yàn)橹悤?huì)被相應(yīng)的酶分解成脂肪酸,并且因?yàn)橹舅岬牟环€(wěn)定性,可在特定條件下分解為酮類和酚類[24]。
如圖5所示,聚類熱圖可將代謝組學(xué)數(shù)據(jù)可視化[25]。每行代表一個(gè)代謝物,每列代表一個(gè)樣本。4個(gè)不同提取物的代謝物水平存在顯著差異,這些特征被投射到熱圖上,并用于樣本聚類。圖中顏色表示相對(duì)含量,紅色塊表示含量高表達(dá),白色塊表示含量低表達(dá)。差異代謝物在4種不同提取物中的相對(duì)含量明顯區(qū)分。全果干粉水提中的積雪草酸、兒茶素、精氨基琥珀酸、富馬酸、柚皮素7-O-β-D-葡萄糖苷含量相對(duì)較高。全果干粉醇提中的海藻糖、柚皮素查爾酮、抗壞血酸鹽、蜜二糖含量相對(duì)較高。果渣干粉水提中的L-蘋果酸、D-木糖、杜松酸、3-甲基硫代丙酸、(-)-表沒(méi)食子兒茶素含量相對(duì)較高。果渣干粉醇提中的羥基異亮氨酸、異鼠李糖、青蟹肌醇、13-L-氫過(guò)氧基亞油酸含量相對(duì)較高。
表2 刺梨提取物差異代謝物的鑒定信息Table 2 Identification information of differential metabolites in Rosa roxburghii Tratt extract
A-全果干粉水提;B-全果干粉醇提;C-果渣干粉水提;D-果渣干粉醇提圖4 刺梨不同提取物差異代謝物的火山圖Fig.4 Volcanic map of different metabolites from different extracts of Rosa roxburghii Tratt
A-全果干粉水提;B-全果干粉醇提;C-果渣干粉水提;D-果渣干粉醇提圖5 不同提取物的聚類熱圖Fig.5 Cluster heat map of different extracts
利用篩選出的差異代謝物KEGG ID進(jìn)行MBRole 2.0的通路分析,得到P<0.01的代謝途徑有11條(P值越小,則該代謝通路的差異性越顯著)。這11條代謝途徑分別是淀粉與蔗糖的代謝、磷酸轉(zhuǎn)移酶系統(tǒng)、半乳糖代謝、谷胱甘肽代謝、酪氨酸代謝、類黃酮生物合成、植物激素的生物合成、β-丙氨酸代謝、鳥(niǎo)氨酸、賴氨酸和煙酸生物堿的生物合成、苯丙烷類化合物的生物合成、ABC轉(zhuǎn)運(yùn)蛋白生物合成。結(jié)合P值與刺梨功能活性物質(zhì)特征,本次實(shí)驗(yàn)的差異物代謝通路為類黃酮生物合成途徑。如圖6所示,參與類黃酮生物合成代謝通路的有4個(gè)代謝物,分別是有柚皮素查爾酮、兒茶素、柚皮素7-O-β-D-葡萄糖苷、(-)-表沒(méi)食子兒茶素[(-)-epigallocatechin,EGC]。
由類黃酮生物合成代謝通路圖可知(圖7),對(duì)香豆酰輔酶a(C6-C3)與3個(gè)丙二酰輔酶a(C3)分子的縮合成柚皮素查爾酮與二苯基丙烷(C6-C3-C6)單元,通過(guò)共軛環(huán)閉合將其轉(zhuǎn)化為柚皮素與黃酮(2-苯基色素-4-酮)主鏈,進(jìn)一步的修飾產(chǎn)生多種結(jié)構(gòu)形式,包括查爾酮、黃烷酮、黃酮以及異黃酮[26]。張帥[27]的研究表明柚皮素查爾酮具有抗腫瘤作用。
A-柚皮素查爾酮;B-(+)-兒茶素;C-(-)-表沒(méi)食子兒茶素; D-柚皮素7-O-β-D-葡萄糖苷?qǐng)D6 類黃酮生物合成代謝通路的4個(gè)代謝物Fig.6 Four metabolites of flavonoid biosynthesis and metabolism pathway
(+)-兒茶素的合成路徑為亮花青素+NADPH+H+生成(+)-兒茶素+NADP++H2O。KIM等[28]證明了植物源性兒茶素具有優(yōu)異的抗感染、抗炎和抗氧化活性。兒茶素和富含兒茶素的植物材料通過(guò)體內(nèi)研究中的炎癥機(jī)制在短期和長(zhǎng)期有效地抑制炎癥應(yīng)激。因此,兒茶素本身或含有兒茶素的營(yíng)養(yǎng)制劑可作為強(qiáng)效抗感染劑或具有優(yōu)異生理活性的功能性食品材料。
(-)-表沒(méi)食子兒茶素的合成路徑為delphinidin+2 NADPH+H+生成(-)-表沒(méi)食子兒茶素+2 NADP+。KIM等[29]的研究表明,用EGC與沒(méi)食子酸處理脂肪細(xì)胞可抑制脂肪細(xì)胞分化,表沒(méi)食子兒茶素可顯著調(diào)節(jié)成熟脂肪細(xì)胞的代謝基因轉(zhuǎn)錄,EGC確實(shí)在調(diào)節(jié)脂肪細(xì)胞分化和代謝方面發(fā)揮了有益的作用。在CHEN等[30]阿爾茲海默癥的研究中,證明了多酚化合物EGC和表兒茶素沒(méi)食子酸酯可以通過(guò)與Cu2+和Zn2+的螯合作用,有效緩解A40抗體聚集,減少活性氧的產(chǎn)生,從而降低Cu2+和Zn2+-A40誘導(dǎo)的神經(jīng)元毒性。
圖7 類黃酮生物合成代謝通路圖(來(lái)源:KEGG PATHWAY Database)Fig.7 Pathway map of flavonoids biosynthesis and metabolism (source:KEGG PATHWAY Database)注:方框?yàn)槊傅膰?guó)際命名編號(hào);普通圓圈為化合物;紅圓圈為本次的4個(gè)代謝物;實(shí)箭頭為反應(yīng)方向;虛箭頭為與其他代謝途徑的關(guān)系。
柚皮素7-O-β-D-葡萄糖苷的合成路徑為黃烷酮7-O-β-D-葡萄糖基轉(zhuǎn)移酶與柚皮素生成UDP和柚皮素7-O-β-D-葡萄糖苷[31],國(guó)內(nèi)外對(duì)其研究較少,但它同屬類黃酮一類,推測(cè)其也具備以上功能活性。在后續(xù)的研究中,可加強(qiáng)對(duì)刺梨抗腫瘤、抗氧化、抗感染、減肥等方面的研究。
代謝組學(xué)是一門交叉學(xué)科,對(duì)于各學(xué)科的綜合運(yùn)用要求較高,能更全面,更精準(zhǔn)地反饋物質(zhì)的整體信息。全面、細(xì)致地了解代謝異常和差異是代謝組學(xué)的一個(gè)重要特性。本研究基于UPLC-HRMS非靶向代謝組學(xué)分析刺梨不同提取物,PCA結(jié)果顯示,不同提取物的生物學(xué)重復(fù)都匯聚在一起,說(shuō)明刺梨不同提取物數(shù)據(jù)的穩(wěn)定性較好,在測(cè)定過(guò)程中未出現(xiàn)明顯偏差。組內(nèi)樣本點(diǎn)都很好地聚集在一起,證明代謝物種類和含量都相近。OPLS-DA模型分析表明各組與其他提取物代謝物差異性較大,并根據(jù)VIP、P、FC值篩選出49種顯著差異代謝物。總體而言,上調(diào)差異代謝物多于下調(diào)差異代謝物,上調(diào)差異代謝物黃酮類最多。經(jīng)過(guò)MBRole 2.0通路分析,這49種差異代謝物參與的代謝通路有11條,其中類黃酮生物合成途徑最為顯著。參與類黃酮生物合成代謝通路的有4個(gè)代謝物,分別是有柚皮素查爾酮、兒茶素、柚皮素7-O-β-D-葡萄糖苷、(-)-表沒(méi)食子兒茶素。兒茶素、柚皮素7-O-β-D-葡萄糖苷含量在全果水提中相對(duì)較高,柚皮素查爾酮含量在全果醇提中相對(duì)較高,(-)-表沒(méi)食子兒茶素含量在果渣水提中相對(duì)較高。這4種代謝物在國(guó)內(nèi)外的研究中展現(xiàn)出優(yōu)異的功能活性特點(diǎn)。本文探究不同提取物的代謝物差異,為以后以代謝組學(xué)為基礎(chǔ)的刺梨研究提供理論參考。在后續(xù)的實(shí)驗(yàn)中,可加強(qiáng)對(duì)刺梨抗腫瘤、抗氧化、抗感染、減肥等方面的研究。