陈琪玉,李宇飞,2,徐 航,2,刘贤青,2,罗 杰,2
(1. 海南大学 热带作物学院/海南省热带生物资源可持续利用重点实验室,海口 570228;2. 海南省三亚崖州湾种子实验室,三亚 572000)
葡萄(grape )是葡萄科葡萄属的藤本类植物,不仅具有悠久的种植历史,也是目前全球种植面积最大以及产量最高的经济水果之一[1]。葡萄属可分为真葡萄和麝香葡萄两个亚属共计70多个种,但具有经济价值的驯化栽培种仅有20多个种且均属于真葡萄属,主要包括4大种群:欧亚种、美洲种、东亚种以及三大种群杂交选育出来的杂交种(以下简称种间杂交种)[2]。目前中国主要种植的类型为欧亚种与杂交种,其中杂交种的抗病性与环境适应性强于欧亚品种[3]。葡萄作为人们喜爱的水果之一,除了果肉的鲜美以外,以葡萄为原料的副产品如葡萄酒、葡萄汁和葡萄干等也受人们的欢迎[4-5]。在当今社会,消费者对食物的追求除最基本的口感外,对食物中的营养成分及含量水平也提出了更多的要求,因此探究葡萄中所含的营养物质的分布与含量成为研究热点[6-7]。
氨基酸、维生素、脂质以及类黄酮等物质作为四类常见营养元素,已有研究人员对其在葡萄中的分布情况开展研究。徐雯等[8]对中国河北省内主要种植的6个葡萄品种进行氨基酸的检测,发现葡萄中主要的氨基酸为谷氨酸、精氨酸、脯氨酸和天冬氨酸等。此外,Tassoni A等[9]通过研究抗病型和感病型葡萄品种的氨基酸水平差异,发现在抗病性强的葡萄品种中,谷氨酰胺、精氨酸、脯氨酸和苏氨酸等氨基酸的含量均高于感病型的葡萄品种。维生素是生物体维持正常生命活动所必需的营养物质,在生物体的生长发育等代谢过程中发挥重要作用。金洪艳等[10]对葡萄中所含的维生素进行检测,发现在葡萄中的维生素主要包括维生素B族和C族。脂质主要在葡萄籽中积累,Pérez-Navarro J等[11]对葡萄中所含的脂质进行检测,发现在果肉中广泛存在棕榈酸、硬脂酸、亚油酸和甘油磷脂等脂质,还发现在葡萄籽中不饱和脂肪酸的比例相对较高。类黄酮是葡萄中种类最多且含量最高的营养代谢物,其在葡萄中的分布与含量受到研究人员的广泛关注,葡萄中主要的类黄酮化合物包括花青素、黄酮醇和黄烷醇等[12]。研究人员还发现类黄酮在葡萄品种间存在很大的积累差异[13],例如,在大多数红葡萄品种中,黄烷醇和花青素是类黄酮化合物的主要组成成分,而在白葡萄品种中,酚类物质的含量相对较低,其中黄酮醇和黄烷醇占主导地位[14]。
代谢组学(metabonomics)是通过组群指标分析,进行高通量检测和数据处理,研究生物体整体或组织细胞系统的动态代谢变化,特别是对内源代谢、遗传变异、环境变化乃至各种物质进入代谢系统的特征和影响的学科。代谢组与基因组、转录组和蛋白组一起构成了系统生物学,代谢物是基因转录与蛋白表达的终端产物,承载着生物体众多的生命活动并体现生物的生命状态[15]。目前代谢领域常用的研究分析平台分别是液相色谱质谱联用系统(LC-MS)、气相色谱质谱联用系统(GC-MS)以及核磁共振(NMR)等。因LC-MS具有高灵敏和高重现性的双重优势,LC-MS是目前代谢组学研究中应用最广泛的检测手段。随着科学技术的不断进步与发展,液相色谱体系也从高效液相色谱体系(HPLC)发展到分离比更高的超高效液相色谱体系(UPLC),体系的升级使得物质分离效果和检测效率得到了大幅度的提升[16]。代谢组学常用物质检测方法有非靶向代谢检测[17]、靶向代谢检测[18]以及广泛靶向代谢检测方法[19],非靶向代谢检测和靶向代谢检测的检测方法各有优劣,而广泛靶向检测方法综合两者优势在代谢领域得到了广泛应用。代谢组学在植物相关研究中至关重要[20],在植物的品种鉴定、中医药学以及园艺作物的品质研究等领域均有广泛应用[21-23]。
葡萄作为高经济价值的水果之一,培育味道鲜美和营养价值高且具有广泛环境适应能力的葡萄品种是研究人员不断努力的方向。本研究收集了起源于欧亚种群(V.vinifera)或种间杂交种群两大种群共20份葡萄栽培品种,利用广泛靶向代谢检测方法对葡萄群体进行高通量代谢检测分析,建立了涵盖多类代谢物质葡萄代谢数据库,比较了几类营养代谢物在葡萄欧亚种和种间杂交种之间的分布情况,意为葡萄高营养品种的培育提供新的思路以及理论基础。
1.1 材 料 本研究使用的20份葡萄种质资源起源于葡萄欧亚种群(V.vinifera)或杂交种群两大种群,主要从北京、山西和山东等地的种质资源中心收集。在挑选材料时根据品种来源、果皮颜色以及品种用途等方面进行考量,所有的材料都使用统一的管理方法,包括栽培、灌溉、施肥、修剪和疾病控制等方面。
1.2 代谢样品提取 先利用真空低温冻干机将葡萄鲜样中的水分去除,将冻干后的样品磨成粉末装在2 mL离心管中。将冻干粉样按每份50 mg进行分装,将分装后的样品置于冰上按照W冻干粉样∶V提取液=1 mg∶8 mL加入由甲醇、乙腈和超纯水按V甲醇∶V乙腈∶V超纯水=2∶2∶1混合配制而成的提取液,再进行10 s涡旋混匀,然后放入加了冰块的超声仪中超声10 min,涡旋超声步骤共重复3次。超声结束后将样品置于预冻好的4℃高速离心机中离心10 min,吸取上清液打入装有0.22 μm尼龙滤膜的针筒进行过滤,将滤液装入含有内插管的棕色代谢瓶中等待上机检测,整个过程在冰上操作。代谢样品提取全过程中,除磨样和样品分装2个步骤外,其余操作均需在黑暗条件下进行。
1.3 代谢样品检测
1.3.1 液相分离参数 分离使用的液相仪器分别为UltiMate 3 000 高效液相色谱仪(赛默飞世尔科技公司)和DGU-20A5R高效液相色谱仪(岛津公司)。检测用的色谱柱为Shim-pack GISS C18(1.9 μm, 2.1 × 100 mm)反向色谱柱(岛津公司),柱温设为40℃。流动相体系为A相(含0.04 %乙酸水溶液);B相(含0.04 %乙酸的甲醇溶液)。液相洗脱时间为20 min,梯度为:0~5 min,流动相B相的比例由5 %逐渐增高至95 %,在B相占比达95 %时保持5 min,15.1 min时迅速降至B相比例5 %并保持至洗脱时间结束。在样品正式检测前必须使仪器达到平衡状态,平衡时间大约为3 min。
1.3.2 高分辨质谱检测参数 使用Q Exactive Plus Orbitrap HRMS仪器进行高分辨数据采集,扫描模式为Full MS/ddMS2,母离子扫描范围100~1 000 Da,二级谱图采集范围50~1 000 Da, 正离子模式下喷雾电压3.5 kV,毛细管温度325℃,加热器温度350℃,鞘气30 arb,辅助气10 arb。
1.3.3 低分辨质谱检测参数 Q-Trap 6500+MS用于进行代谢物的定量检测,检测模式为sMRM。ESI源的参数设置:温度500℃,喷雾气(GS1)50,辅助加热气(GS2)60,气帘气(CUR)35 psi。正离子扫描模式:离子喷雾电压(IS)5 000 V,去簇电压(DP)50 V,碰撞电压(CE)30 V。对代谢信号进行扫描时,分成多个文件,每个文件内放置1 000个离子对进行扫描,扫描窗口设定为70,总的循环时间为0.8 s。数据采集完后由Analyst软件导出然后导入MultiaQuant 软件进行后续的定量分析。
1.4 数据统计与分析
1.4.1 代谢数据的定性和定量 非靶检测所得的原始数据均由Compound Discoverer 3.1软件进行了峰对齐、反褶积和空白扣除等工作。代谢物定量数据使用 MultiQuant 3.0.3 软件进行处理。
1.4.2 代谢数据的处理由EXCEL完成 代谢物的聚类热图由R语言的pheatmap R包完成;主成分分析(PCA)使用线上网站HIPLOT (https://hiplot.com.cn/)完成。
2.1 葡萄代谢物数据库的建立 利用高分辨仪器UPLC-Q Exactive Plus Orbitrap对葡萄群体的代表性样本进行非靶向代谢组检测,并基于以下两个原则对代谢信号数据进行初步筛选:(1)是否拥有二级碎片;(2)前体离子(Q1)与产物离子(Q3)之差是否小于14,以此将质量较差以及冗余的信号去除。最终,通过非靶检测共计得到了2 353个代谢信号。为进一步对筛选得到的信号进行定量检测以及二次筛选,利用Q-Trap 6500+ MS对葡萄群体的代谢信号进行靶向检测。利用MultiQuant软件对靶向检测所得数据进行处理,去除峰型较差以及信噪比在10以下的代谢信号,最终得到了1 328个代谢信号,完成葡萄代谢物数据库的初步建立。
葡萄代谢物数据库初步建立后,利用多种方法对葡萄代谢物数据库中未知代谢信号进行物质注释。运用高分辨仪器配套的CD(Compound Discover)软件对采集到的代谢物质进行一个初步匹配,但因软件的局限性,匹配结果还需结合碎片信息以及出峰时间进行人工核对,再通过由已知到未知的顺序,结合实验室已有代谢数据库,根据Q1、Q3、碎片信息和保留时间(样品和检测仪器的不同会影响RT值,但出峰时间的差异范围大约在0.5~1 min)等信息注释了373个代谢物。再利用标准品核验,将标准品采用样品同样的检测方法进行代谢检测,根据所得的RT、Q1和Q3的值与已有数据进行匹配,通过该方法注释了247个物质。
此外,还对部分未知代谢信号进行了人工注释,如代谢物GRP0421的保留时间是4.97 min(图1-A),质荷比(m/z)为627.156 2,根据其二级谱图(图1-B),发现其主要的碎片为303.049 6,除此之外还有1个小的碎片465.102 6。通过分析该物质的质谱信息,其碎裂模式与糖基化类黄酮类物质的裂解方式非常相似,结合常见黄酮类物质的出峰时间大约在4~8 min,因此笔者推测该物质为类黄酮的糖基化修饰物。对该物质的二级谱图进行分析,303.049 6为该类黄酮物质母核的特征质荷比,其特征碎片303与另一个主要碎片465以及465与母离子627的中性丢失特征值均为162,恰好符合六碳糖丢失1个水(MW-H2O,162)的质荷比,因此推测在该物质中存在2个六碳糖和1个羟基,再结合CD软件的匹配结果人为核对后,将该物质注释为Quercetin 3-glucosylgalactoside,结构式以及物质断裂模式如图1-C所示。
图1 代谢物的鉴定,代谢数据库物质分类
本研究通过一系列代谢物注释方法后,在去除冗余值后,建立了1个包含768个已知代谢物的葡萄代谢数据库,涵盖了类黄酮、脂质、氨基酸、维生素、有机酸、萜类、多酚及酚胺等九类代谢物(图1-D),其中类黄酮和脂质化合物的占比较高。
2.2 葡萄代谢谱分析 基于所得代谢物数据库对葡萄群体进行定量分析后,本研究展开了初步的代谢轮廓分析。根据所得768个已知代谢物对20份葡萄品种利用R语言进行主成分分析(PCA)和聚类热图分析(图2)。由图2-A可知代谢物可以将葡萄欧亚种和杂交种区分开。为了探究是哪些化合物在葡萄不同种群之间存在差异,对80份栽培葡萄品种进行聚类热图分析,由图2-B可知类黄酮化合物、脂质、萜类化合物以及多酚化合物是造成葡萄品种存在代谢差异的几类主要化合物,其中萜类化合物在葡萄杂交种中的积累相对高于葡萄欧亚种而类黄酮物质和脂质在葡萄不同种群中的积累存在较大差异。此外,本研究统计了葡萄代谢库中的768个代谢物的变异系数,发现有482个代谢物的变异系数超过了60%,这说明在本研究中测得的代谢物在葡萄欧亚种和杂交种之间存在较大的积累差异(图2-C)。
图2 葡萄群体代谢组分析
2.3 葡萄欧亚种和杂交种差异代谢物分析 为明确在葡萄欧亚种和杂交种之间的代谢物的差异,本研究结合变量权重值(VIP值>1)和差异倍数的计算(Fold Chang<0.8或Fold Chang>1.2)两种分析代谢物的差异的方法,选取两种方法的交集,最终共筛选到了185个在葡萄欧亚种和种间杂交种之间的差异代谢物,其中在葡萄种间杂交种中呈现上调的有90种代谢物,有95种代谢物是在葡萄杂交种中出现下调(图3-A)。
图3 葡萄不同种群间差异代谢物分析
对185个差异的代谢物的注释信息进行统计,共有172个代谢物属于常见代谢物注释的类别,包含了14种氨基酸及其衍生物、85种类黄酮化合物、40种脂质、4种有机酸、3种酚胺类化合物、3种植物激素、8种萜类化合物以及15种维生素(图3-B)。这说明了在葡萄不同群体之间代谢物的差异主要由类黄酮、脂质和氨基酸等化合物引起。
2.4 葡萄欧亚种和杂交种营养代谢物的差异利用前述在葡萄中测得到的代谢数据库,对氨基酸、维生素、脂质和多酚等4类常见营养物质在葡萄的欧亚种和杂交种之间的含量差异进行分析。如图4-A所示,与欧亚种相比杂交种的总氨基酸的相对含量更高,其中脯氨酸(Proline)和谷氨酰胺(Glutamine)的相对含量分别高于葡萄欧亚种1.5倍和2.6倍,而酪氨酸(Tyrosine)和L-异亮氨酸(L-Isoleucine)在葡萄欧亚种与杂交种的积累含量无显著差异(图4-B),说明与欧亚种相比杂交种中氨基酸类化合物的积累相对更多。
图4 葡萄欧亚种和杂交种间营养代谢物含量差异
本研究发现,在葡萄中维生素B族的代谢物积累存在差异,维生素B3物质在葡萄杂交种中的积累极显著地高于欧亚种,但维生素B6化合物则在欧亚种中的积累更高(图4-C)。本研究还发现,维生素C化合物L-抗坏血酸(L-Ascorbic acid)在属于欧亚种的品种中的相对含量普遍高于属于杂交种的葡萄品种(图4-D)。甘油磷酯类化合物在葡萄不同种群中的积累也存在一定的差异,葡萄杂交种中总脂肪酸的相对含量显著高于欧亚种(图4-E),但亚麻酸(Linolenic acid)的相对含量在葡萄欧亚种中略高于杂交种(图4-F)。
类黄酮是葡萄中非常重要的一类营养物质,本研究通过靶向代谢检测方法共检测得到225个多酚类化合物,主要包含黄酮、黄酮醇、黄烷酮、黄烷醇、异黄酮、查尔酮、花青素和芪类化合物等九大类型。对检测到的多酚类物质进行聚类分析,如图5所示多酚类化合物在葡萄欧亚种和种间杂交种中存在较大的积累差异,其中黄酮醇类化合物在葡萄欧亚种和种间杂交种的积累模式与总的多酚类化合物在葡萄不同种群中的积累模式相似(图5-A),说明黄烷醇是形成多酚化合物在葡萄中存在显著积累差异的主要代谢物。与黄酮醇类物质不同,黄烷醇类化合物在葡萄欧亚种中的总相对含量极显著高于种间杂交种(图5-B)。对4种检测到的花青素在葡萄中的含量进行分析,发现花青素类物质在种间杂交种的品种中相对含量普遍高于葡萄欧亚种的品种(图5C-F),说明葡萄不同种群之间的杂交可以增加葡萄花青素的积累。
图5 葡萄欧亚种和杂交种多酚类物质含量差异
本研究以葡萄果实群体为材料,基于代谢组学策略建立了1个包含多类营养物质的葡萄代谢物数据库。尽管该数据库包含768个注释代谢物,但依旧有许多能被检测到的代谢信号仍未被注释出来。目前领域内常用的物质注释方式主要包括利用标准品比对,以及运用数据库匹配软件例如CD、MSDIAL和Sirius等通过导入已测得的高分辨数据进行匹配处理[24-25],这些物质注释方式方法多少都存在假阳性率过高的局限性。人工解谱也是研究人员会运用到的物质注释方法之一,但这一方法对相关专业知识要求较高因此使用范围更加狭窄。本研究主要是运用液相质谱的靶向代谢检测方式对葡萄群体进行研究,以及主要运用标准品匹配和软件匹配对代谢物进行注释,所运用的检测方法手段还有很大的提升空间,想要构建包含更多注释代谢物的葡萄代谢物数据库,还需进一步开发利用更多的方式策略加以完善。
本研究发现,脯氨酸、谷氨酰胺和维生素B3等代谢物在葡萄种间杂交种的相对含量显著高于葡萄欧亚种,这些代谢物已被报道与植物抗性相关,而种间杂交种群的培育是研究人员为了提高葡萄的抗性以及对于不同环境的适应性而产生的,本研究结果从代谢角度的方面验证了这一说法。本研究还发现几种参与葡萄果皮颜色形成的花青素在葡萄种间杂交种的品种中的相对含量高于欧亚品种,已有研究[26]表明运用嫁接砧木的育种方式培育出来的杂交葡萄品种,其果皮中花青素化合物的积累增加,促进了果皮着色,而葡萄颜色是迎合消费者审美的一种非常重要的农艺性状,这一结果为今后培育色泽鲜艳的葡萄品种提供了理论基础。通过利用营养成分不同的葡萄种群能为今后提高葡萄产业的价值品质提供代谢方面的物质基础。
随着代谢组学的高度发展,以代谢组学为基础的多组学研究已成为作物育种研究领域常用的手段之一,在蓝莓、咖啡和香菜等经济植物的研究上已被应用[27-29]。本研究建立的葡萄代谢物数据库为后续基于代谢组结果进一步开展葡萄基于代谢组的多组学分析奠定了代谢研究基础,有助于有针对性地对葡萄代谢品质等性质进行遗传机制的解析,为今后葡萄品种选育提供新方向与代谢基础。