比较代谢组学揭示不同生长期筇竹笋的代谢物差异.pdf

成分分析 食品科学 2023,Vol.44,No.14 237比较代谢组学揭示不同生长期筇竹笋的 代谢物差异王卫华1,2，杨晓琴2，桑正林1，杨顺强1，赵 平2，刘 云2，朱国磊2，叶夏英1,*（1.昭通学院农学与生命科学学院，滇东北高原特色农业研究中心，云南 昭通 657000；2.西南林业大学 西南山地森林资源保育与利用教育部重点实验室，云南 昆明 650224）摘 要：为探究筇竹笋在生长过程中代谢物的动态变化和阶段特异性，采用基于超高效液相色谱-串联质谱的广泛靶向代谢组学，对2 个不同生长期出土高度分别为10cm和25cm的筇竹鲜笋样本进行分析。结果表明，2 个生长期的筇竹笋样本共鉴定出13 类875 个代谢物，其中糖及醇类、黄酮类等成分的相对含量在生长期为25cm样本中明显上升。多元统计分析结果显示，2 个不同生长期的筇竹笋代谢物存在明显阶段特异性，生长期10cm与25cm的筇竹笋共存在283 个显著差异代谢物，其中102 个代谢物被注释到82 条代谢通路中。通过综合分析得出，20 种竹笋主要营养成分在2 个不同生长期的笋中无显著差异，而与口感风味、活性功能有关的代谢物在生长期为25cm的样本中相对含量更高，且25cm筇竹笋样本中苦味因子相对含量更低。因此与生长期为10cm的筇竹笋相比，生长期为25cm的筇竹笋综合经济价值更高，可能是更加适合的采收规格。本研究结果可为筇竹资源的综合开发和可持续利用提供理论参考。关键词：筇竹笋；不同生长期；比较代谢组学；差异代谢物；超高效液相色谱-质谱联用Comparative Metabolomics Revealed Metabolite Differences in Bamboo Shoots(Chimonobambusa tumidissinoda Hsueh&T.P.Yi ex Ohrnberger)at Different Growth StagesWANG Weihua1,2,YANG Xiaoqin2,SANG Zhenglin1,YANG Shunqiang1,ZHAO Ping2,LIU Yun2,ZHU Guolei2,YE Xiaying1,*(1.Research Center for Plateau Characteristic Agriculture in Northeast Yunnan,Department of Agronomy and Life Science,Zhaotong University,Zhaotong 657000,China;2.Key Laboratory for Forest Resources Conservation and Utilization in the Southwest Mountains of China,Ministry of Education,Southwest Forestry University,Kunming 650224,China)Abstract:In order to explore the dynamic changes and stage specificity of metabolites in bamboo shoots during its growth,widely targeted metabolomics based on ultra-high performance liquid chromatography tandem mass spectrometry(UPLC-MS/MS)was used to analyze fresh bamboo shoots harvested at heights of 10 and 25 cm.The results showed that a total of 875 metabolites belonging to 13 categories were identified from bamboo shoots,including sugars,alcohols and flavonoids,whose relative contents significantly increased in 25 cm-high bamboo shoots.Moreover,the results of multivariate statistical analyses showed that the metabolites in bamboo shoots had an obvious stage specificity.Totally 283 significantly differential metabolites were identified between the two growth stages,of which 102 metabolites were annotated to 82 metabolic pathways.The contents of 20 major nutrients did not significantly differ between the two growth stages;however,the contents of metabolites related to taste,flavor and bioactive function were higher,while the content of bitter compounds was lower in the 25 cm-high bamboo shoots.Therefore,compared with the 10 cm-high bamboo shoots,the 25 cm-high bamboo shoots have a higher comprehensive economic value and are more suitable for harvesting.The 收稿日期：2022-09-26基金项目：国家自然科学基金地区科学基金项目（32260307）；云南省应用基础研究计划青年项目（2019FD059）；昭通市凤凰计划项目第一作者简介：王卫华（1989）（ORCID:0000-0002-6503-2164），男，主管药师，硕士，研究方向为天然产物化学。E-mail:*通信作者简介：叶夏英（1990）（ORCID:0000-0003-2066-945X），女，副教授，博士，研究方向为竹类系统生物学。E-mail: 238 2023,Vol.44,No.14 食品科学 成分分析results of this study can provide a theoretical reference for the comprehensive development and sustainable utilization of Chimonobambusa tumidissinoda.Keywords:bamboo shoots;different growth stages;comparative metabolomics;differential metabolites;ultra-high performance liquid chromatography-mass spectrometry DOI:10.7506/spkx1002-6630-20220926-278中图分类号：S644.2 文献标志码：A文章编号：1002-6630（2023）14-0237-08引文格式：王卫华,杨晓琴,桑正林,等.比较代谢组学揭示不同生长期筇竹笋的代谢物差异J.食品科学,2023,44(14):237-244.DOI:10.7506/spkx1002-6630-20220926-278.http:/WANG Weihua,YANG Xiaoqin,SANG Zhenglin,et al.Comparative metabolomics revealed metabolite differences in bamboo shoots(Chimonobambusa tumidissinoda Hsueh&T.P.Yi ex Ohrnberger)at different growth stagesJ.Food Science,2023,44(14):237-244.(in Chinese with English abstract)DOI:10.7506/spkx1002-6630-20220926-278.http:/筇竹（Chimonobambusa tumidissinodaHsueh&T.P.YiexOhrnberger）为地下茎复轴混生中小型竹类，隶属禾本科（Poaceae），竹亚科（Bambusoideae）、方竹属（ChimonobambusaMakino）1，是收录于中国珍稀濒危保护植物目录的保护竹种，也是我国特有种质资源2。筇竹自然分布于我国西南金沙江下游地区，是云南昭通重要的生物资源，昭通天然筇竹林面积占全国天然筇竹林面积的70%以上3。筇竹笋具有肉厚、质脆、味美之特性，是深受消费者青睐的食品。果蔬的口感风味与其代谢物的种类、含量及这些代谢物在发育成熟过程中的动态变化等密切相关，例如番茄的口感风味与番茄中的有机酸、可溶性糖类等代谢物有关，其中的33 种代谢物与消费者的喜好存在明显关联4-5。WangRuochen等6对生长成熟过程中猕猴桃的风味物质变化进行了研究，结果显示，猕猴桃风味相关代谢物在其成熟过程中具有明显的阶段特异性。GaoQuan等7研究发现了新鲜竹笋中的苦涩味主要来源于竹笋的代谢产物L-苯丙氨酸。植物次生代谢产物是植物体内合成的一类与生长发育不相关的小分子物质，它们通常在植物适应环境和应对胁迫的过程中起到关键作用8，植物来源的食品也具有种类丰富的次生代谢产物（如酚酸类、黄酮类、生物碱类等），而这些代谢物与食品发挥抗氧化、抗炎、和改善代谢综合征的生物活性与功能有关9-12。目前，对筇竹笋营养成分的研究结果表明，筇竹鲜笋富含蛋白质、氨基酸、可溶性糖和粗脂肪等营养物质，必需氨基酸比例在34%以上13。杨奕等14通过分析测定3 个不同出土高度筇竹笋主要营养成分，认为出土高度为1020cm的筇竹笋营养品质和经济价值最高。然而现有的研究仅局限于笋中的营养物质，对筇竹笋在生长过程中相关代谢产物的动态变化和阶段特异性，以及与筇竹笋口感风味相关的其他物质和次生代谢产物等的研究尚不够深入。广泛靶向代谢组学是一种基于高通量检测和多元数据处理，能够全面、快速对样品中小分子代谢物进行定性、定量分析，从而系统分析样品代谢物差异水平和种类变化的新兴技术15-18。因此，本研究利用基于超高效液相色谱-串联质谱（ultra-highperformanceliquidchromatography-tandemmassspectrometry，UPLC-MS/MS）的广泛靶向代谢组学方法对2 个不同生长期（出土高度为10cm和25cm）筇竹鲜笋进行分析，进一步探究筇竹笋代谢物以及其在笋生长过程中的动态变化和阶段特异性，为筇竹资源的综合开发和可持续利用提供参考。1 材料与方法1.1 材料与试剂筇竹笋样品于2021年5月在云南昭通市大关县三江口采集，经昭通学院农学与生命科学学院生物科学教研室鉴定确认为禾本科植物筇竹（C.tumidissinodaHsueh&T.P.YiexOhrnberger）的嫩苗。分别量取出土高度为10cm和25cm的无虫害竹笋，脱壳后在顶部、中部和基部各取5 g，混合后放置于冻存管中用液氮速冻后放置于80 冰箱待测。水、甲醇、乙腈、甲酸（均为色谱纯）美国Merck公司；标准品二甲基亚砜（色谱纯）云南西力生物技术股份有限公司。1.2 仪器与设备NexeraX2UPLC仪 岛津（中国）有限公司；4500QTRAP串联质谱仪 美国AppliedBiosystems公司；SB-C18色谱柱（2.1mm100mm，1.8m）美国安捷伦公司；Scientz-100F普通型冷冻干燥机 宁波新芝生物科技股份有限公司；MM400冷冻混合研磨仪 德国莱驰公司；PTX-FA210S电子天平 华志（福建）电子科技有限公司；XH-D型旋涡振荡器 美国赛默飞世尔科技公司；JIDI-16D台式高速离心机 广州吉迪仪器有限 公司。成分分析 食品科学 2023,Vol.44,No.14 2391.3 方法1.3.1 样品提取样品放置于冻干机中真空冷冻干燥后利用研磨仪研磨（30 Hz，1.5min）至粉末状，称取100mg的粉末，溶解于1.2mL70%甲醇提取液中，后置于旋涡振荡器，每30min涡旋一次，每次持续30s，共涡旋6 次，样本置于4 冰箱过夜，离心（12 000 r/min，10min）后，吸取上清，用微孔滤膜（0.22m）过滤样品，并保存于进样瓶中，2 个不同生长期的筇竹笋样本，分别重复取样3 次，用于UPLC-MS/MS分析。质控样本由不同组的筇竹笋样本提取物等量混合制备而成，与分析样本采用相同处理和检测方法，重复测定3 次。仪器分析过程中，每10 个检测分析样本中插入一个质控样本，以监测分析过程的重复性。1.3.2 UPLC-MS/MS采集条件UPLC条件：柱温40，流速0.35mL/min，进样量4L；流动相：A相为超纯水（加入0.1%甲酸），B相为乙腈（加入0.1%甲酸）；洗脱梯度：0.00min，95%A，5%B；0.009.00min，95%5%A，5%95%B，保持1min；10.0011.10min，5%95%A，95%5%B，并以5%B平衡至14min。MS条件：电喷雾离子源；离子源温度550；离子喷雾电压5 500 V（正离子模式）/4 500 V（负离子模式）。三重四极杆扫描使用多反应监测模式，根据每个时期内洗脱的代谢物，在每个时期监测一组特定的多反应监测模式离子对。1.4 数据处理利用软件Analyst1.6.3处理质谱数据，用MultiaQuant软件，基于自建数据库（metwaredatabase，MWDB），对样本的代谢物进行质谱定性定量分析，对所有代谢物质谱峰进行积分，并对同一物质质谱峰进行积分校正19。多元统计包括无监督模式识别的主成分分析（principalcomponentanalysis，PCA）和有监督模式识别的正交偏最小二乘判别分析（orthogonalpartialleastsquares-discriminantanalysis，OPLS-DA）。OPLS-DA在原始数据进行log2转换后，再进行中心化处理，然后利用R软件进行分析。差异代谢物筛选采取将OPLS-DA模型的变量重要性投影（variableimportanceinprojection，VIP）值、差异倍数（fold change，FC）值相结合的方法。显著差异代谢物基于2个筛选标准：1）FC2和FC0.5的代谢物；2）VIP1的代谢物。2 结果与分析2.1 总代谢物分析经过代谢组分析，2 个生长时期的筇竹笋样本中共检测鉴定出875 种代谢物（图1），按代谢物类型将代谢物分为13 类（图2）。利用面积归一法获得各类物质的相对含量（表1），25cm筇竹笋较10cm样本相比，糖及醇类、黄酮类、生物碱类相对含量上浮明显，脂质、氨基酸及衍生物占比有所下降。0.010.841.0923456789 10 11 12 13 140.41.2A1.00.80.60.2?107?/min0.014.1323456789 10 11 12 13 140.41.2B1.00.80.60.2?107?/minA.正离子模式；B.负离子模式。图 1 筇竹笋样本代谢物多反应监测多峰图Fig.1 MRM chromatograms of bamboo shoot samples 160140120100806040200?1561471019794776756231917417?图 2 筇竹笋样本中检测到的代谢物分类Fig.2 Classification of metabolites detected in bamboo shoot samples根据代谢组分析结果发现，在875 种代谢物中，筇竹笋样品中的脂质成分种类最多，相对丰度最高，鉴定出脂质156 种，包括74 种游离脂肪酸，18 种甘油酯及94 种其他脂类等。其中亚油酸和-亚麻酸在875 个鉴定代谢物的相对含量分别达3.94%和2.90%。亚油酸和-亚麻酸是人体必需的不饱和脂肪酸，具有降脂和保护心脑血管等作用20-22。本研究表明筇竹笋中不饱和脂肪酸占脂肪酸的主导地位，这与其他竹笋研究结论23基本一致。另有研究显示，竹笋膳食纤维能够降低人体低240 2023,Vol.44,No.14 食品科学 成分分析血清总胆固醇、低密度脂蛋白胆固醇和动脉粥样硬化 指数24-25。根据本研究结果推测，筇竹笋中丰富的不饱和脂肪酸等脂质成分可能是筇竹笋具备改善代谢综合征（血脂异常）作用的重要物质基础。表 1 2 个不同生长期筇竹笋样本各类代谢物相对含量Table 1 Relative contents of metabolites in bamboo shoots at two different growth stages物质类别相对含量/%10cm样本25cm样本脂质38.1532.68酚酸类5.986.40黄酮类1.843.18氨基酸及其衍生物10.638.10生物碱13.1916.80有机酸9.608.98糖和醇类4.7910.13核苷酸及其衍生物12.4410.62萜类1.140.65维生素类1.211.27木脂素和香豆素0.150.16醌类0.080.07其他类0.790.95从筇竹笋中鉴定出氨基酸及其衍生物共97 种，其中包含蛋白氨基酸17 种，必需氨基酸8 种。利用测出的17 种蛋白氨基酸的相对含量与17 种氨基酸相对含量总和进行比较，求出各氨基酸的相对含量（表2）。从表2可以看出，筇竹笋样本中L-谷氨酸的相对含量最高，达23.59%；必需氨基酸在17 种氨基酸中的相对含量达34.13%。根据游离氨基酸的呈味特征26，从表2还可以分析得出，筇竹笋中鲜味氨基酸如L-谷氨酸、L-天冬氨酸、L-赖氨酸，甜味氨基酸如L-脯氨酸相对含量较高，而苦味氨基酸如竹笋中最主要的苦味因子L-苯丙氨酸相对含量较低。这可能是筇竹笋具有鲜甜口感，受消费者喜好的关键原因和重要物质基础。表 2 筇竹笋样本中氨基酸的相对含量Table 2 Relative contents of amino acids in bamboo shoots序号氨基酸相对含量/%相对含量在本类别中排名1L-谷氨酸23.591/972L-天冬酰胺11.784/973L-脯氨酸9.855/974L-缬氨酸8.106/975L-谷氨酰胺7.268/976L-赖氨酸6.1610/977L-天冬氨酸5.5711/978L-亮氨酸4.9912/979L-异亮氨酸4.9113/9710L-甲硫氨酸4.4915/9711L-精氨酸2.7019/9712L-酪氨酸2.4320/9713L-苏氨酸2.2821/9714L-色氨酸2.1523/9715L-丝氨酸1.6927/9716L-苯丙氨酸1.0634/9717L-组氨酸0.9936/9718总必需氨基酸34.13从筇竹笋中还鉴定出有机酸77 种（柠檬酸、L-苹果酸、茉莉酸等）、糖及醇类67 种（3-脱氢-L-苏糖酸、D-蔗糖、D-果糖等）、维生素19 种（烟酰胺、D-泛酸、抗环血酸等）、核苷酸及其衍生物56 种（腺嘌呤、鸟苷等）、其他类17 种（4-胍基丁醛等）。据此可以看出筇竹笋中维生素种类丰富，营养价值高，其中有机酸和可溶性糖种类丰富、相对含量较高，可能也是筇竹笋风味口感形成的重要物质基础。此外，从筇竹笋中鉴定出比较丰富的次生代谢产物，包括黄酮类101 种（4-羟基-5,7-二甲氧基黄酮、芹菜素-8-C-(2”-葡萄糖基)阿拉伯糖苷、槲皮素-3-O-芸香糖苷等）、酚酸类147 种（3-O-阿魏酰奎宁酸、邻苯二甲酸二(2-乙基己基)酯、4-羟基苯甲酸等）、生物碱94 种（2-羟基-7-甲氧基-1,4-苯并噁嗪-3(2H)-酮、五羟色胺、哌啶等）、萜类23 种（贝壳杉烯酸等）、木脂素和香豆素17 种（东莨菪内酯-7-O-葡萄糖醛酸苷等）、醌类4 种（大黄酚-9-蒽酮等）。由于植物次生代谢产物通常具备一定生物活性，因此这些物质成分可能是筇竹笋发挥活性功能，适应环境的重要物质基础。2.2PCA结果为了解各分组样本之间的主要代谢物和组内样本之间的差异，运用PCA对2 个不同生长期的筇竹笋样本和质控样本（mix）进行分析。由图3可知，PC1贡献率为44.40%，PC2贡献率为40.07%，前5 个PC的累计可解释变异为93.86%，所有样品坐标点均分布于95%置信区域，说明该模型可分析不同生长期筇竹笋代谢物的差异情况。质控样本和各组内样本基本重合，表明检测的稳定性和可靠性良好。生长期10cm的筇竹笋样本与生长期25cm的筇竹笋样本组间差异大，分离趋势明显，说明2 个不同生长时期的样本代谢物存在明显差异。?20?40?200204020400PC2?40.07%?PC1?44.40%?25 cm C25 cm A25 cm Bmix 01mix 03mix 0210 cm A10 cm C10 cm B图 3 2 个不同生长期筇竹笋样本与质控样本质谱数据的PCA得分图Fig.3 PCA score plot of MS data of bamboo shoots at two different growth stages and those of quality control samples2.3 OPLS-DA结果为进一步分析2 个不同生长期筇竹笋的代谢物差异，筛选2 个组别的差异代谢物，对样本进行OPLS-DA。从图4A可知，生长期为10cm的筇竹笋与25cm的筇竹笋在代谢物上存在显著差异，这一结果与PCA的成分分析 食品科学 2023,Vol.44,No.14 241结果相似。为验证OPLS-DA的准确性，本模型对数据进行了200 次随机排列组合，模型的R2X=0.801、R2Y=1、Q2=0.992，表明所构建的OPLS-DA模型稳定可靠，具备良好预测性，适合进一步筛选2 个不同生长期筇竹笋的差异代谢物。?20?30?20?100201030010?10PC2?7.64%?PC1?72.4%?10 cm B10 cm C10 cm A25 cm A25 cm B25 cm C图 4 2 个不同生长期筇竹笋样本OPLS-DA得分图Fig.4 OPLS-DA plot of bamboo shoots at two different growth stages2.4 2 个不同生长期筇竹笋差异代谢物的分析基于OPLS-DA模型的VIP、单变量分析的FC和相应的筛选标准，对分组样本的差异代谢物进行筛选。如 图5A所示，生长期为10cm的筇竹笋与生长期为25cm筇竹笋共存在283 个显著差异代谢物，2 个不同生长期的筇竹笋显著差异代谢物个数占所有检测出代谢物总数的32.34%，表明2 个不同生长期筇竹笋的代谢物差异较为显著。其中生长期25cm的样本相较于10cm样本上调代谢物154 个（包括酚酸类44 个、黄酮类32 个、糖及醇类19 个等），下调代谢物129 个（包括生物碱类22 个、氨基酸及衍生物15 个、脂质类15 个等）。将衡量代谢成分显著差异的定量信息FC通过数学变换为log2FC，按log2FC值筛选出排名前20的显著差异代谢物（图5B）。从图5B可知，在25cm的筇竹笋样本中，柠檬酸、没食子酸、绿原酸等有机酸，酚酸类和多种黄酮类成分出现显著上调，而N-阿魏酸五羟色胺、吲哚等成分出现显著下调。为直观反映出筇竹笋中主要营养成分和风味物质的阶段特异性，将2.1节分析得到的20 种主要营养成分的具体差异情况列出（表3），结合2 个生长期主要营养成分的VIP、FC值和显著差异代谢物的筛选标准，从表3可知，亚油酸、-亚麻酸及多个氨基酸在筇竹笋2 个不同生长期的差异并不显著，而在25cm的筇竹笋样本中，苦味因子L-苯丙氨酸比10cm样本的相对含量略低。0.0?16?8?4?2?1 0 1 241680.90.61.2A0.3VIPlog2FC?129?576?154B2020202020log2FC?-7-O-?-4-O-?(?)?-7-O-?2,3-?-3-?2-?-1-?-8-O-?3-?-1-?-4-O-?4-?-6-C-(2”-?)?(3-O-?)?-2”-O-?N-?C?-N-?1,3-O-?1-O-?-3-O-?-4-O-?2-O-?6-O-?-7-O-?-4-O-?-3-O-?N-?16:0(2n?)3,5,6,7,8,3,4-?13-?-6Z,9Z-11E-?9,10-?-12,13-?19.2416.1815.5114.9314.4114.0313.8213.6012.8512.4812.4112.1511.6411.1710.7710.7410.6710.6410.3410.22?3.48?3.66?4.18?4.19?4.76?5.49?8.42?8.48?10.25?10.63?11.20?11.68?11.97?12.23?12.69?13.03?13.11?13.13?14.84?15.28?图 5 2 个不同生长期筇竹笋样本差异代谢物火山图（A）和排名前20 差异代谢物FC柱状图（B）Fig.5 Volcano map(A)and histogram(B)of top 20 most significantly differential metabolites from bamboo shoots at two different growth stages表 3 2 个不同生长期筇竹笋样本中20 种主要营养成分VIP值和FC值Table 3 VIP and FC of 20 major nutrient components in bamboo shoots at two different growth stages代谢物VIPFC差异性亚油酸1.020.91不显著-亚麻酸1.081.09不显著-亚麻酸0.811.06不显著L-谷氨酸1.170.73不显著L-天冬酰胺1.150.60不显著L-脯氨酸1.170.62不显著L-缬氨酸0.601.05不显著L-谷氨酰胺1.170.57不显著L-赖氨酸1.170.59不显著L-天冬氨酸1.170.40显著下调L-亮氨酸1.170.74不显著L-异亮氨酸1.170.75不显著L-甲硫氨酸1.170.52不显著L-精氨酸1.170.44显著下调L-酪氨酸1.070.82不显著L-苏氨酸1.160.53不显著L-色氨酸1.170.21显著下调L-丝氨酸0.031.00不显著L-苯丙氨酸1.170.66不显著L-组氨酸1.170.43不显著注：单变量分析的FC值为生长期为25cm样本中代谢物相对含量与10cm样本代谢物相对含量的比值。242 2023,Vol.44,No.14 食品科学 成分分析2.5 2 个不同生长期筇竹笋差异代谢通路的分析为进一步探索2 个不同生长期筇竹笋中代谢通路的差异，利用京都基因和基因组百科全书（KyotoEncyclopediaofGenesandGenomes，KEGG）对统计筛选出的显著差异代谢物进行功能注释。283 个显著差异代谢物中有102 个代谢物被注释到82 条代谢通路中，其中被注释的差异代谢物数目排前3的KEGG通路类型分别为代谢途径、次生代谢物的生物合成和ABC转运蛋白，分别占总数的75.49%、43.14%和14.71%。如图6所示，生长期为25cm的筇竹笋样本与10cm的样本相比，黄酮与黄酮醇的生物合成、半乳糖代谢、-亚麻酸代谢等通路富集明显且呈现较明显的上调趋势。植物激素信号转导、苯并噁嗪类生物合成等通路呈现出明显上调或下调趋势，但富集程度不高。与竹笋木质化进程关系密切的苯丙素类的生物合成通路富集程度和上调趋势尚不明显。0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0?/?A?P?0.3020400.20.11.00.50.00.51.0B/P1.000.750.500.250.0010203040图 6 2 个不同生长期筇竹笋样本差异代谢物KEGG富集图（A）及 差异丰度得分图（B）Fig.6 KEGG enrichment map(A)and differential abundance score map(B)of differential metabolites in bamboo shoots at two different growth stages2.6 2 个不同生长期筇竹笋代谢物阶段特异性综合分析已有研究结果表明，竹笋中的黄酮类成分具有抗氧化、降血糖等多种活性功能27-28。综合2 个不同生长期筇竹笋的差异代谢物和代谢通路结果发现，在生长期为25cm的筇竹笋样中，包括槲皮素3-O-新橘皮糖苷、槲皮素3-O-芸香糖苷（芦丁）、槲皮素7-O-芸香糖苷、荭草素-7-O-葡萄糖苷、槲皮素-4-O-葡萄糖苷（绣线菊苷）等27 个黄酮类代谢物的相对含量较生长期为10cm的筇竹笋样出现显著增加。图7显示了25cm筇竹笋样中鉴定出的两个相对含量最高的黄酮类代谢物：4-羟基-5,7-二甲基黄酮和槲皮素-3-O-芸香糖苷（芦丁），在总代谢物相对含量分别为0.95%和0.43%。4-羟基-5,7-二甲基黄酮是典型的二氢黄酮，具有抗晚期糖基化终末产物（advancedglycationendproducts，AGEs）的活性功能，AGEs则与多种退行性疾病有关29，而芦丁是黄酮醇二糖苷，具有抗氧化和抗糖尿病的活性功能 30-31。据此可以分析得出，筇竹笋中含有能够抗氧化和改善代谢综合征（糖尿病）的黄酮类成分，且具有活性功能的黄酮类成分在25cm的筇竹笋中较10cm的含量显著增加。茉莉酸（图7）是由-亚麻酸合成的脂类派生植物激素，是植物多个生长发育过程的关键调节信号 分子32，同时可驱动多个次生代谢产物（包括黄酮类、花青素等）在植物中的积累33。在麦草中的一项研究表明，外源茉莉酸还能诱导麦草中的酚类物质含量增多，从而增加麦草抗氧化和抗炎的活性功能34。在本研究中，2 个不同生长期的筇竹笋样本都检测到了茉莉酸这一代谢物，与10cm样本相比，生长期为25cm的筇竹笋样本中茉莉酸相对含量出现显著上调（VIP=1.17、FC=2.35）。同时，代谢途径结果（图6）显示-亚麻酸代谢和植物激素信号转导通路出现了富集与上调。据此综合分析推测，筇竹笋生长过程中，茉莉酸由-亚麻酸经生物合成途径快速富集产生并在筇竹笋生长过程中起到关键作用，可能是黄酮类代谢物富集上调的关键调控因子。能否通过外源性的茉莉酸诱导筇竹笋活性功能成分的增加从而增强筇竹笋活性功能有待进一步的研究探索。苯并噁嗪类是存在于禾本科植物中的一类次生代谢产物，也是一种重要的化感物质，具有抗虫、抗病的防御功能35。在本次筇竹笋代谢组学研究中，发现包括2-羟基-7-甲氧基-1,4-苯并噁嗪-3(2H)-酮、2-羟基-1,4-苯并噁嗪-3(2H)-酮葡萄糖苷等多个苯并噁嗪类生物碱（图7），代谢组学结果还显示，生长期25cm的筇竹笋样本相比10cm样本，苯并噁嗪类生物合成通路出现下调（图6），苯并噁嗪生物碱相对含量下降，呈现出阶段特异性，这可能与筇竹笋出笋早期需要更多化感物质对抗生物胁迫有关。筇竹笋是滇东北高原特色农产品，也是产区农民重要的经济收入来源。然而，目前在采笋季节，筇竹笋成分分析 食品科学 2023,Vol.44,No.14 243采收的规格参差不齐，部分地区采收者甚至会“见笋就打”，这对筇竹资源的可持续开发利用十分不利。本研究结果显示，相对于生长期为10cm的筇竹笋，生长期为25cm筇竹笋中的20 种主要营养成分变化并不显著，而与口感风味密切相关的有机酸、糖与醇类、鲜甜味氨基酸等代谢物相对含量上升，苦味因子L-苯丙氨酸相对含量下降，与抗氧化和改善代谢综合征等活性功能密切相关的黄酮类代谢物显著上升。因此与生长期为10cm的筇竹笋相比，生长期为25cm的筇竹笋具有更高的综合经济价值，可能是更加适合的采收规格。NHO2-?-7-?-1,4-?-3(2H)-?OOHONHOOO2-?-1,4-?-3(2H)-?OHOHOHOHHOOOOHHOGlcOHOHRha?-3-O-?OOOOHO4-?-5,7-?HOOOHH?图 7 筇竹笋中部分代谢物的结构式Fig.7 Structures of some metabolites in bamboo shoots3 结 论运用基于UPLC-MS/MS的广泛靶向代谢组学对昭通产筇竹笋进行分析。从2 个不同生长期（出土高度分别为10cm和25cm）筇竹笋样本中共鉴定出13 类875 个代谢物，表明筇竹笋是营养成分丰富的食品原料。多元统计分析显示，2 个不同生长期筇竹笋代谢物存在明显的阶段特异性，生长期为10cm的筇竹笋与生长期为25cm的筇竹笋共存在283 个显著差异代谢物，其中102 个代谢物被注释到了82 条代谢通路中。通过综合分析得出，20 种主要营养成分在生长期为10cm和25cm的筇竹笋样本中差异并不显著，而与口感风味、活性功能有关的代谢物在生长期为25cm的样本中相对含量更高，且25cm筇竹笋样本中苦味因子相对含量更低。因此与生长期为10cm的筇竹笋相比，生长期为25cm的筇竹笋综合经济价值更高，可能是更加适合的采收规格。本研究结果可为筇竹资源的综合开发和可持续利用提供重要的理论参考。参考文献：1 朱石麟,马乃训,傅懋毅.中国竹类植物图志M.北京:中国林业出版社,1994.2 傅立国,金鉴明.中国植物红皮书:稀有濒危植物:第1册M.北京:科学出版社,1992:318-319.3 董文渊,邱月群,王逸之,等.天然筇竹居群形态遗传多样性J.东北林业大学学报,2016,44(5):101-103.DOI:10.3969/j.issn.1000-5382.2016.05.023.4ZHANGB,TIEMAND,JIAOC,etal.Chilling-inducedtomatoflavorlossisassociatedwithalteredvolatilesynthesisandtransientchangesinDNAmethylationJ.ProceedingsoftheNationalAcademyofSciences,2016,113(44):12580-12585.DOI:10.1073/pnas.161391011.5TIEMAND,ZHUGT,RESENDEM,etal.AchemicalgeneticroadmaptoimprovedtomatoflavorJ.Science,2017,35:391-394.DOI:10.1126/science.aal1556.6WANGRC,SHUP,ZHANGC,etal.Integrativeanalysesofmetabolomeandgenome-widetranscriptomerevealtheregulatorynetworkgoverningflavorformationinkiwifruit(Actinidia chinensis)J.NewPhytologist,2022,233(1):373-389.DOI:10.1111/nph.17618.7GAOQ,JIANGH,TANGF,etal.Evaluationofthebittercomponentsofbambooshootsusingametabolomicsapproach