大头菜又名芥菜，是一种十字花科芸苔属草本植物，常用来制成酱腌菜，因其脆嫩爽口而深受消费者喜爱。大头菜含有丰富的蛋白质、糖类、维生素以及膳食纤维。通常，大头菜经过3～4个月短期发酵即可成熟，作为佐餐菜食用。

为了研究不同超长周期发酵大头菜代谢物差异和主要差异代谢物及其关联代谢通路，成都大学的张钰麟、赵志平*和四川省食品检验研究院的徐跃成*等人从自贡市某大头菜加工企业采集了工业化生产的3种不同发酵年限的大头菜，利用LC-MS/MS技术分析其差异代谢物，以期为定向调控大头菜特征代谢物的生物合成提供一定的理论依据。

1、大头菜代谢物的PCA

如图1所示。0 a组、5 a组与10 a组均分布在不同的区域内，说明3种不同发酵周期大头菜代谢物之间存在显著差异，因此基于PCA得分图确定3个处理组：0 a与5 a为A1，5 a与10 a为A2，0 a与10 a为A3。 和 Q 2 是评判PCA模型的主要参数，其中 、 Q 2 ＝0.64表明所建立的PCA模型具有较好的可解释度。

2 大头菜代谢物的OPLS-DA

对3个处理组，采用OPLS-DA进一步分析均得到了较好的分离，如图2A 1～C 1所示。每个处理组中两两组别均显著分离，说明发酵周期的不同导致大头菜代谢物种类及含量产生显著差异，这可能与长期发酵大头菜中耐盐微生物活力一直处于较高水平有关，如乳酸菌可产生多种氨基酸，通过厌氧发酵产生丙酸、乳酸；酵母菌可产生有机醇，基于利用还原糖能力强的特点，使其更能适应发酵过程的高盐、高酸环境，加之代谢物之间的相互作用使得组别间差异较明显。OPLS-DA模型参数为 均大于0.5，且 ＝0.988说明该模型可以解释和预测两种大头菜代谢物之间的差异。

为了进一步验证模型是否出现过拟合现象，采用200次循环迭代置换检验，结果如图2A 2～C 2所示。 Q 2回归线与 Y轴交点都在负半轴，表明OPLS-DA模型稳定可靠，不存在过拟合现象。

3、大头菜差异代谢物的筛选与鉴定

本研究以OPLS-DA模型变量投影重要性（VIP）＞1.3且P＜0.05为筛选标准。如图3A所示，3个处理组差异代谢物主要以有机酸及其衍生物、氨基酸及其衍生物和核酸及其衍生物为主，发酵时间的延长使得微生物通过酶促反应利用大头菜中底物能力增强。研究显示，随着发酵的进行，盐厌氧菌属、盐单胞菌属和色盐杆菌属逐渐成为相对丰度较高的优势菌群；在高盐的环境下，优势耐盐菌群充分利用能源物质，使得大分子营养物质被降解成小分子物质，如氨基酸、脂肪酸等，这体现在A3较A2及A1组差异代谢物种类及个数均更加丰富。利用Venn图对各处理组的差异代谢物分析表明，A1、A2、A3共鉴定到53、113种和174种差异代谢物，如图3B所示。取3个处理组共有的19种差异代谢物定义为主要差异代谢物，包括3种氨基酸及其衍生物（L-苯丙氨酸、N-乙酰亮氨酸、4-乙酰氨基丁酸），5种有机酸（6-羟基烟酸、3-(2-羟基苯基)丙酸、2-氨基苯磺酸、L-3-苯乳酸、烟酸），3种核酸及其衍生物（N-D-核糖嘌呤、鸟嘌呤、7-甲基腺嘌呤），2种脂肪酸（辛酸、癸二酸），2种酮糖（β-D-果糖-6-磷酸、D-核酮糖），1种黄酮类物质（根皮苷），1种酚类物质（苯酚），1种酚酸（水杨酸），1种生物碱（β-咔啉）。

氨基酸及其衍生物在植物中发挥各种功能，除了参与合成蛋白质外，还与能量代谢、植物的胁迫反应等有关。发酵周期的不同使得微生物代谢功能基因表达程度不同，尤其是氨基酸的相关代谢活动。氨基酸的种类直接决定发酵蔬菜的滋味，Xiao Yangsheng等通过对四川泡菜、江西盐菜和东北酸菜研究发现，发酵蔬菜都富含谷氨酸、苯丙氨酸、丙氨酸、天冬氨酸等氨基酸，均在本研究中被检出。徐亚洲等研究证实，泡菜发酵过程中氨基酸含量变化趋势并非一致，这在本研究中体现在天冬氨酸为上调的差异代谢物，而丝氨酸与丙氨酸为下调的差异代谢物。多种氨基酸共同作用有助于大头菜滋味的形成，赋予大头菜鲜辣等味觉特征。

本研究鉴定的大部分有机酸含量随着发酵周期延长而增加。研究显示，发酵蔬菜中有机酸主要来源于蔬菜原料及发酵过程中微生物的代谢作用。随着发酵周期的增加，乳酸及乙酸等有机酸由于糖代谢作用逐渐积累增加 。 次生代谢物以产量较大的初生代谢物为前体，通过糖代谢、莽草酸延伸、氨基酸延伸、乙酸延伸途径所合成的结构复杂的化合物。次生代谢物种类很多，包括生物碱、色素、抗生素等。由于色素含量的不同，0、5 a和10 a大头菜的色泽表现出显著差异，如图4所示。酚酸及其衍生物广泛分布在植物中 ，发酵过程中酚酸的积累有助于大头菜口感、酸度、风味的形成 ，如在芥菜和客家擦菜中可检测出40余种酚酸 ，与本研究检出的大多数酚酸变化趋势一致，但主要差异代谢物水杨酸与传统发酵的芥菜变化趋势不一致。传统芥菜水杨酸含量随发酵周期延长含量逐渐下降，而本研究水杨酸含量呈上升趋势。与水杨酸代谢有关的酶为苯丙氨酸解氨酶，其通过介导肉桂酸途径从而合成水杨酸。邹同雷等 发现龙须菜中苯丙氨酸解氨酶在高盐浓度下活性增加，而5 a和10 a大头菜的食盐含量高于0 a，因此5 a和10 a大头菜水杨酸含量显著高于0 a。 本研究鉴定到的主要差异生物碱为 β -咔啉，其含量呈先上升后下降的趋势， β -咔啉生物碱在新药开发上具有广阔的前景，基于其在长期发酵大头菜中的含量变化规律，后续可针对5 a大头菜进行相关功能性代谢物的生物合成调控。

4、大头菜主要差异代谢物热图分析

如图5所示，主要差异代谢物在各个组别中含量分布具有一定的规律特征。0 a中3种下调的主要差异代谢物与其他两组差别明显，其中以3-(2-羟基苯基)丙酸和 β - D -果糖-6-磷酸变化最为显著。5 a组中含量较高的有 β -咔啉、癸二酸、辛酸、苯酚，它们多为脂肪酸、生物碱和酚类物质，可能对大头菜发酵前期褐色的形成起重要作用。10 a中含量较高的有根皮苷、 L -苯丙氨酸、2-氨基苯磺酸、 L -3-苯乳酸、烟酸、6-羟基烟酸、 N - D -核糖嘌呤。

5、主要差异代谢物通路分析

对主要差异代谢物所富集的代谢通路进行分析，将主要差异代谢物映射到KEGG数据库中，主要富集到苯丙氨酸代谢（Impact＝0.47）、磷酸戊糖途径（Impact＝0.13）、果糖和甘露糖代谢（Impact＝0.21）3 条代谢通路上，如图6所示。芳香族氨基酸除了作为蛋白质的组成成分外，还是许多次生代谢物生物合成的前体 ，通过莽草酸途径在质体内部合成。由于缺乏完整的糖酵解过程，使其只能通过磷酸烯醇式丙酮酸特异性的磷酸转运体从胞浆导入。 L -苯丙氨酸是许多酚类化合物的前体物质，包括木质素、缩合单宁以及苯丙烷类物质。在一些植物中，苯丙氨酸生物合成的主要途径是芳香酸途径 。 β - D -果糖-6-磷酸作为连接磷酸戊糖途径与果糖和甘露糖代谢的主要代谢产物，其含量随着发酵周期的延长而下降，可能是由于6-磷酸葡萄糖异构酶活性减弱，使来源于糖酵解过程的 α - D -葡萄糖-6-磷酸向其转化能力减弱；作为呼吸作用的底物，通过三羧酸循环、磷酸戊糖途径及果糖和甘露糖代谢途径转化为三磷酸腺苷以维持正常的生命活动 。

结 论

基于LC-MS/MS探究了工业化生产的发酵0、5 a和10 a的传统大头菜代谢物的差异，差异代谢物主要以有机酸、氨基酸和核苷酸衍生物为主。基于Venn图得到3个处理组共有的19种主要差异代谢物，将主要差异代谢物通过KEGG数据库检索出3 条最显著的代谢通路，为基于代谢通路定向调控大头菜代谢物合成提供一定的理论基础。

大头菜在发酵过程中发生复杂的生物化学反应，代谢物种类繁多。深入解析长期发酵大头菜代谢物组成及变化规律是大头菜后续产品开发的前提，也是未来的主要研究方向之一。此外，微生物是影响发酵大头菜代谢物组成最为关键的因素之一，利用宏基因组学、宏转录组学、蛋白质组学及代谢组学等多组学技术深度揭示工业化生产长期发酵大头菜的微生物多样性及群落结构，明确功能微生物及活性基因，阐明微生物调控长期发酵大头菜代谢物合成的机理是未来的研究方向。同时，基于代谢通路利用定向调控技术加速特征代谢产物生物合成进程也是工业化生产亟需解决的技术问题。

本文《长期发酵大头菜差异代谢物及其关联代谢通路分析》来源于《食品科学》2022年43卷22期192-198页，作者：张钰麟，陈泓帆，赵志平*，王 卫，聂 鑫，程 杰，张 崟，刘达玉，徐跃成*，罗淮良。DOI:10.7506/spkx1002-6630-20211220-236。点击下方 阅读原文即可查看文章相关信息。

图片来源于文章原文及摄图网。