栽培方式与生长阶段对灵芝子实体生物活性成分谱及抗氧化能力的影响：正常与鹿角状子实体的比较分析

《Food Chemistry: X》：Effects of cultivation methods and growth stages on bioactive compound profiles and antioxidant capacity in Ganoderma lucidum: a comparative analysis of normal versus antler-shaped fruiting bodies

　　

：0.5% CO2胁迫诱导灵芝（Ganoderma lucidum）产生鹿角状形态并抑制子实体生物量。鹿角状灵芝（AGL）在第7天表现出更强的抗氧化能力，而正常栽培灵芝（CK）在成熟期（第40天）抗氧化能力更优。酚酸和类黄酮可能是两种表型抗氧化差异的关键代谢物。加权基因共表达网络分析（WGCNA）筛选出与抗氧化相关代谢物显著富集的棕色和绿色模块。网络药理学和分子对接预测了其多靶点协同抗氧化机制。

论文解读：

一、研究背景与问题

灵芝（Ganoderma lucidum）作为一种具有药用价值的大型真菌，富含多糖、三萜类等生物活性成分，具有抗氧化、免疫调节、抗肿瘤及保肝等药理作用，在现代功能食品工业中得到广泛应用。灵芝还含有酚类化合物和类黄酮等抗氧化成分，是开发功能性膳食补充剂和天然抗氧化食品的优质原料。然而，在人工大规模栽培过程中，较高的二氧化碳浓度会显著影响子实体的分化和发育，诱导灵芝产生鹿角状的形态异常（AGL）。与正常栽培的灵芝（CK）相比，这些由CO₂诱导的形态变异体在生长动力学和生物量积累模式上表现出显著差异。研究表明，高浓度CO₂会抑制细胞周期调控、DNA复制和蛋白酶体功能相关基因的表达，进而抑制细胞增殖并损害细胞壁生物合成，最终导致鹿角状形态的形成。同时，菌丝体从幼龄到成熟期的发育进程会动态重塑次生代谢途径，直接影响生物活性化合物的合成与积累，从而调节灵芝的药用和营养特性。

目前，不同灵芝表型在发育各阶段酚酸和类黄酮的积累动态仍缺乏系统表征，支撑阶段依赖性抗氧化活性差异的关键代谢物也尚未被系统鉴定。这一知识空白不仅制约了灵芝栽培实践的标准化和优化，也阻碍了其深加工产业的高值化发展。此外，生物体内活性氧（ROS）等自由基的过度积累会诱导氧化应激，促进炎症反应、肿瘤细胞增殖以及糖尿病、骨关节炎等慢性疾病的发生发展。与合成抗氧化剂相比，天然抗氧化剂具有更优的安全性和更低的副作用风险。酚酸和类黄酮作为灵芝中具有强抗氧化活性的天然次生代谢产物，通过直接清除自由基和抑制脂质过氧化发挥作用，同时还能调节microRNA（miRNA）的表达，间接调控细胞内抗氧化酶活性。

二、研究内容与意义

针对上述研究空白，研究人员选取两种灵芝表型，在三个关键发育阶段采集样本，测定了生长参数、体外抗氧化能力、总酚含量（TPC）、总黄酮含量（TFC）、总多糖含量（TPS）和可溶性糖含量（SSC）。通过整合广泛靶向代谢组学、加权基因共表达网络分析（WGCNA）、网络药理学和分子对接技术，系统阐明了培养条件和生长阶段对灵芝代谢谱的影响，预测了调控抗氧化活性的潜在酚酸和类黄酮单体，并从化合物-靶点相互作用、分子靶点参与和通路水平调节三个互补视角阐释了其协同抗氧化机制。该研究为灵芝栽培的定向优化和天然抗氧化食品配料的开发提供了理论基础。

三、主要技术方法

研究人员采用的主要技术方法包括：广泛靶向代谢组学分析（基于UPLC-MS/MS技术，由武汉迈维代谢生物科技有限公司完成）；加权基因共表达网络分析（WGCNA，用于构建代谢物共表达网络并分析模块与生理性状的关联）；网络药理学分析（整合SwissTargetPrediction靶点预测、GeneCards疾病靶点筛选、STRING蛋白互作网络构建及Metascape富集分析）；分子对接（采用AutoDock Vina进行配体-受体结合能计算）。样本来源为新疆石河子农垦科学院提供的灵芝菌株。

四、研究结果

4.1 表型与生理学分析

在第7天，两组均处于菌丝萌发阶段，形态无差异；第20天，CK形成肾形小子实体，AGL则出现分枝状鹿角结构；第40天，CK菌盖完全展开，AGL分枝伸长变细。鲜重（FW）在两组中均随时间递增，第20天和第40天CK显著高于AGL（第40天为AGL的2.31倍）。抗氧化能力测定显示，AGL在第7天和第20天的DPPH清除活性显著高于CK；ABTS测定中AGL在第7天和第20天活性更高，但第40天CK超过AGL约28%；FRAP结果呈现一致的时间模式。TPC和TFC的动态变化证实了这一趋势：第7天两者均以AGL为高，第40天CK显著高于AGL。TPS在各时间点均以CK为高，MDA含量在CK中持续偏高，表明AGL的氧化损伤相对较低。双因素方差分析显示，培养条件对FRAP和TPC的影响最强（η²p分别为0.930和0.924），表明抗氧化能力和酚类积累主要由CO₂浓度主导。

4.2 灵芝代谢物变化

代谢组学分析共鉴定出867种代谢物，分为12大类。主成分分析（PCA）显示PC1和PC2累计贡献率达72.63%，CK与AGL在代谢层面显著分离。正交偏最小二乘判别分析（OPLS-DA）模型验证R²Y值为1.000，Q²值均超过0.9，模型稳健可靠。差异代谢物（DAMs）筛选结果显示，三个发育阶段分别鉴定出213、275和340个差异代谢物，且各阶段下调代谢物数量均多于上调代谢物。

4.2.1 差异代谢物分析

酚酸分析显示，AGL_1中龙胆酸（Gentisic Acid）含量为CK_1的2.16倍，其对位二羟基取代结构赋予其接近没食子酸的强抗氧化活性。第40天，CK_3中富集了托品酸、3-(4-羟基苯基)丙酸、DL-3-苯乳酸和1-O-咖啡酰-β-D-木糖等高度生物活性的酚酸单体。黄酮分析显示，CK在所有阶段黄酮含量均高于AGL；AGL_1中wistin和槲皮素-3-O-(6"-O-乙酰基)葡萄糖苷的倍数变化分别达20.99和19.62，支撑了AGL早期的抗氧化潜力。

4.2.2 KEGG富集分析

第7天，CK中核苷酸代谢和类黄酮生物合成途径活性更高；第20天，AGL中苯丙烷代谢途径逐渐激活，驱动酚酸和黄酮糖苷的持续合成；第40天，长期高CO₂胁迫显著抑制AGL的初级代谢途径，而苯丙烷途径活性达到峰值，使类黄酮逐渐成为AGL的优势活性化合物。

4.3 生理与代谢物的WGCNA分析

WGCNA将代谢物分为九个模块。棕色和绿色模块与抗氧化相关参数（DPPH、ABTS、FRAP、TPC、TFC）呈正相关，其中棕色模块与TPC（r=0.79）、FRAP（r=0.82）和DPPH（r=0.66）显著正相关。棕色模块富集酚酸、类黄酮、氨基酸衍生物和不饱和脂肪酸，枢纽代谢物包括2,4,6-三羟基苯甲酸、龙胆酸和槲皮素糖苷。韦恩图分析鉴定出90个重叠代谢物，占棕色和绿色模块代谢物的48.65%，其中鉴定出4种重叠酚酸，表明酚酸是灵芝抗氧化能力差异的潜在驱动因素。

4.4 网络药理学结果与分析

酚酸网络药理学分析鉴定出56个共同靶点，核心靶点包括PTGS2、PPARG、EGFR、ESR1和HSP90AA1。排名前五的酚酸代谢物为3-(4-羟基苯基)丙酸、香草酸、DL-3-苯乳酸、托品酸和2,4,6-三羟基苯甲酸甲酯。GO富集分析显示靶点显著富集于炎症反应、氧化应激调控和凋亡等生物过程；KEGG富集分析显示靶点主要集中于NF-κB和PI3K-Akt通路。黄酮网络药理学分析鉴定出79个共同靶点，核心靶点包括TNF、PTGS2、PPARG、EGFR和BCL2。黄酮类化合物如二氢山柰酚、西里辛醇、没食子儿茶素3-O-没食子酸酯等表现出多靶点特征。KEGG富集分析显示靶点显著富集于PI3K-Akt、MAPK、HIF-1和p53通路。两类化合物共享PTGS2、PPARG和EGFR等核心靶点，构成协同抗氧化效应的分子基础。

4.5 分子对接结果与分析

分子对接结果显示，五种酚酸与靶蛋白的结合能均低于?5.5 kcal/mol，其中3-(4-羟基苯基)丙酸和DL-3-苯乳酸与HSP90AA1的结合能分别为?7.1和?7.6 kcal/mol。五种黄酮与靶蛋白的结合能均低于?6.8 kcal/mol，其中山柰酚-4'-O-葡萄糖苷和二氢山柰酚与PTGS2的结合能分别为?9.9和?10.0 kcal/mol，黄酮类化合物整体结合能力优于酚酸。但考虑到酚酸在灵芝中的总浓度更高且具有更优的细胞通透性和体内生物利用度，其对整体抗氧化活性的贡献更为显著。

五、结论

本研究通过整合培养条件与生长阶段，系统研究了常规栽培和0.5% CO₂胁迫对灵芝表型、生物量、酚酸和类黄酮积累及体外抗氧化活性的影响。结果表明，高浓度CO₂诱导鹿角状畸形子实体的形成并在早期增强抗氧化能力；然而，长期暴露会扰乱初级和次生代谢途径，抑制化合物合成和生物量积累。相比之下，常规栽培的灵芝表现出规则的形态发育和稳定的代谢稳态，能够在成熟期大量积累高生物活性的酚酸和类黄酮，从而获得更优的综合抗氧化品质。研究预测酚酸可能是灵芝抗氧化活性的主要贡献者，而黄酮发挥辅助协同作用。在实际生产中，应采用正向通风方式生产高抗氧化品质的灵芝，避免长期暴露于高浓度CO₂。本研究建立的综合分析系统可为其他食药用菌的相关研究提供参考。