发酵诱导的非洲本土块茎分子重塑：木薯与cocoyam的非靶向代谢组学研究

《Food Chemistry》：Fermentation-induced molecular remodeling in African indigenous tubers: cassava and cocoyam

　　

木薯和cocoyam是撒哈拉以南非洲的主要膳食主食，在食用前通常经过自然发酵处理。尽管发酵能够减少抗营养化合物并改善食品质量，但其分子效应仍未被充分表征。研究人员采用基于非靶向质谱（mass spectrometry, MS）的代谢组学方法，结合计算注释流程，比较了这两种块茎中发酵诱导的分子重塑，结果显示二者呈现截然不同的响应。在木薯中，773个显著变化的特征中718个（92.9%）含量下降，表明其以分解代谢过程为主。在cocoyam中，响应更为平衡，1013个显著特征中有385个（38.0%）富集，包括与蛋白水解加工一致的三肽和二肽。类别分析（class analysis）、分子网络（molecular networking）及基于参考微生物模型的通路富集分析揭示了块茎特异性特征：在木薯中，嘌呤代谢是最富集的通路，脂肪酸是相对于其整体代谢组下降幅度最大的化学类别；而cocoyam则表现出更强的氨基酸通路富集。与生氰苷（cyanogenic glycoside）相关的特征含量下降，与解毒作用一致。生物转化预测还提示了当前数据库中不存在的推定发酵产物，突出表明这些块茎的化学特性尚未得到充分表征。

发酵诱导的分子重塑是食品科学中一个日益受到关注的领域，但非洲本土块茎在这一过程中的分子层面变化长期缺乏系统表征。木薯（Manihot esculenta）和cocoyam（Colocasia esculenta）是撒哈拉以南非洲重要的高碳水化合物块茎作物，为当地数亿人口提供主要膳食能量。二者均含有需经加工去除的抗营养因子，其中木薯含毒性氰化物，cocoyam含草酸钙晶体，而自然发酵是世代沿用的传统去毒与保藏手段。已有生化研究证实发酵能够提高抗氧化活性、改善蛋白质溶解度并降低氰化物含量，但这些结论大多来自近似组成、pH、可滴定酸度和微生物丰度等宏观指标，发酵如何在分子层面重塑这两种块茎的代谢物组，此前几乎无人涉及。高分辨质谱与计算代谢组学的快速发展使对复杂食品基质中数千种分子特征的系统追踪成为可能，而这类技术尚未充分应用于非洲传统发酵块茎。研究人员据此假设，组成和发酵条件相近的木薯与cocoyam在发酵后应具有相似的分子组成，并通过比较非靶向代谢组学对这一假设进行了检验，结果发现二者的分子响应实际上截然不同。该研究发表在《Food Chemistry》。

研究人员采集了购自南非林波波省Thohoyandou市场的木薯与cocoyam样品，在University of Venda食品科学与技术系中试工厂加工为粉末后进行自然（自发）发酵，发酵条件为室温（21–26?°C）厌氧浸没发酵72小时，每个条件设置三个生物学重复，以未发酵样品为对照。主要技术方法包括：基于液相色谱-四极杆飞行时间质谱（liquid chromatography-quadrupole time-of-flight mass spectrometry, LC-qTOF-MS/MS）的非靶向代谢组学数据采集；使用MZmine进行特征检测、对齐与缺口填充；使用SIRIUS（含CANOPUS）结合BioTransformer扩展候选数据库进行结构注释和化学类别归属；在GNPS2平台进行基于特征的分子网络（feature-based molecular networking, FBMN）分析；使用Mummichog算法基于Acetobacter、Lacticaseibacillus rhamnosus、Lactobacillus plantarum和Saccharomyces cerevisiae四种参考微生物模型进行通路富集分析；采用Welch's t检验与Benjamini-Hochberg错误发现率（false discovery rate, FDR）校正进行差异丰度分析，并采用Mann-Whitney U检验评估化学类别的非比例变化。

3.1 木薯与cocoyam在发酵中呈现不同的代谢组响应。主成分分析（principal component analysis, PCA）显示两种块茎的发酵与对照样品均明显分离，但分离模式相反：木薯发酵样品聚集紧密而对照分散，cocoyam则对照聚集紧密而发酵样品分散。差异丰度分析表明，木薯发酵以分解代谢为主，773个显著变化特征中718个（92.9%）含量下降，仅55个（7.1%）富集；cocoyam的响应更为平衡，1013个显著特征中385个（38.0%）富集、628个（62.0%）下降，富集特征中包含与蛋白水解加工一致的三肽和二肽。

3.2 发酵影响木薯和cocoyam中不同的化学类别。在NPClassifier通路水平，木薯中脂肪酸类相对于整体代谢组背景显著下降（校正后p=0.000576，n=283），是唯一显著偏离背景的通路；N-酰基胺和中性鞘糖脂的富集提示检测到的部分信号可能来自微生物生物量而非块茎自身转化。cocoyam中氨基酸和肽类显著区别于背景代谢组响应（校正后p=0.000005，n=158），糖基二酰甘油在发酵后丰度增加，与革兰氏阳性菌膜脂重塑相关；小肽显著富集而脂肪酰胺下降。生氰苷类特征在两种块茎中均呈下降趋势，与发酵脱毒作用一致。

3.3 发酵块茎的分子多样性在现有数据库中代表性不足。SIRIUS与GNPS2谱库匹配共获得218个推定注释特征，对应179个独特化合物，其中51个仅见于木薯、91个仅见于cocoyam、37个为二者共有。数据库交集分析显示，食品专用数据库FooDB对218个特征仅匹配99个；BioTransformer预测产生的两个推定结构未被任何参考数据库收录，可能是尚未被目录化的转化产物。代表性差异特征包括：甲基硫代腺苷在发酵木薯中显著富集（log₂FC=9.36，p<0.01），与微生物多胺代谢相关；阿魏酰胍丁胺在发酵cocoyam中显著富集（log₂FC=3.21，p<0.01），属于植物防御化合物；姜糖脂B在发酵木薯中下降（log₂FC=?4.83，p<0.05）；黄嘌呤在发酵cocoyam中下降（log₂FC=?5.92，p<0.01）。

3.4 分子网络将结构相关的发酵产物与底物相关联。两种块茎的分子网络均以鞘脂和甘油酯簇为主。木薯网络中观察到由BioTransformer预测的单酰甘油转化链（1-单亚油酸甘油酯、单油酸甘油酯、单硬脂酸甘油酯）和鞘脂系列（中性鞘糖脂到三种神经酰胺），且这些特征在发酵后均下降或不变，网络中也存在完全由下降特征组成的簇，与分解代谢特征一致。cocoyam网络包含更广泛的化合物家族，富集节点倾向于在特定分子家族内聚集；其中9个BioTransformer注释特征在发酵后富集，7个为推定注释的三肽（如Phe-Leu、Phe-Phe、Leu-Ala-Phe），指向活跃的蛋白水解活性，与乳酸菌蛋白水解系统作用于植物蛋白一致。相比之下，木薯中37个BioTransformer注释特征在发酵后无一富集。

3.5 通路富集揭示块茎特异性的发酵代谢路线。在四种参考微生物模型下，木薯中最清晰的信号为嘌呤代谢，在全部参考模型中显著富集（p=0.0015至0.0078），14个显著命中化合物在发酵后全部下降，包括脱氧腺苷和脱氧鸟苷；脂质和非脂质化学亚类模型显示溶血磷脂酰肌醇、溶血磷脂酰胆碱和溶血磷脂酰乙醇胺几乎完全消耗（98–100%），同时伴随类二十烷酸和类固醇下降，以及生氰苷耗尽。cocoyam中甘氨酸、丝氨酸和苏氨酸代谢为主导富集通路（p=0.0011至0.0163，富集比6至9），色氨酸代谢、半胱氨酸和甲硫氨酸代谢也显著富集；非脂质模型捕获了氨基酸消耗，但总命中数少于木薯。两种块茎之间不存在共享的参考模型富集通路。

讨论指出，木薯发酵以现有代谢物的分解为主，少数富集化合物的信号更可能来自微生物生物量；cocoyam发酵则在消耗原有代谢物的同时生成包括三肽和二肽在内的新产物，这与cocoyam相较于木薯更高的天然蛋白质含量一致。该研究在分子层面首次描绘了这两种非洲本土块茎发酵过程中截然不同的化学重塑模式，同时复现了生氰苷消耗等已知发酵效应，验证了分析流程的可靠性。研究还发现大量检测到的化合物未被现有食品或生物数据库充分收录，提示这类传统发酵食品的化学空间可能远未被认识。结论部分强调，若两种在近似组成层面看似营养相似的块茎均能呈现出如此不同的分子发酵谱，那么其他非洲传统发酵主食的分子多样性很可能远大于目前文献所记录的范围；对这些分子转化的理解将有助于改善发酵块茎产品的营养质量和安全性，并促进对具有文化重要性但科学上尚未充分开发的食品中生物活性化合物的发掘。研究局限包括：绝大多数注释为推定水平（MSI Level 2–3），未使用标准品验证；氰化物未直接定量，生氰苷结果仅作为质谱证据；实验采用实验室条件下面粉的浸没发酵，可能与传统整块块茎的家庭发酵不同；通路富集基于参考微生物模型而非实际表征的微生物群落；此外，实验设计无法区分微生物活动与加工过程中释放的内源植物酶（糖苷酶、酯酶、蛋白酶）各自对代谢重塑的贡献。