香菇菌丝体生理成熟与ROS、氨基酸及碳通量互作网络研究

为解决香菇菌丝体在“转色期”生理成熟机制不清，影响子实体优质高产的问题，研究人员通过多组学技术，解析了菌丝成熟过程的动态分子事件。研究揭示了活性氧(ROS)信号、氨基酸稳态和碳通量三者间的时序协同变化，并构建了一个整合工作模型，为理解从菌丝生长到子实体发育的转变提供了可验证的框架，对优化香菇栽培实践具有重要意义。

香菇，作为一种营养丰富、风味独特的食用菌，在全球范围内被广泛栽培和消费。然而，在香菇的工厂化栽培中，有一个关键而神秘的阶段——“转色期”（又称菌丝体生理成熟期）。菌丝体在完全“吃透”培养基料后，并不会立即长出子实体，而是需要一段时间的“酝酿”。在此期间，菌丝表面会形成一层褐色菌膜。只有当菌丝达到充分成熟的生理状态，才能在适宜的环境条件下，高效地扭结分化，形成产量高、品质优的子实体。反之，如果菌丝未成熟或过度成熟，都会导致减产、畸形菇甚至出菇失败。因此，深入理解菌丝体在转色期内部发生的生理和分子变化，对于精准判断成熟度、实现香菇的稳产高产至关重要。尽管科学家们推测此阶段可能涉及次生代谢物合成、黑色素积累、活性氧(ROS)调控以及营养能量稳态等一系列复杂过程，但究竟是哪些关键通路在起主导作用？它们之间又如何相互作用？这些问题一直是香菇生物学和栽培技术中的“黑箱”。为了解开这些谜团，宋晓霞、陈雨佳等研究人员在《Scientia Horticulturae》上发表了一项研究，为我们揭示了香菇菌丝体生理成熟背后一个潜在的分子互作网络。

为了系统解析菌丝体成熟的分子机制，研究团队采用了多时间点动态分析结合多组学整合的研究策略。他们从同一批次接种的KS11香菇菌棒中，选取了八个不同培养时间点（D35至D105）的样品进行研究。研究首先通过观察菌棒的转色特征和测定出菇率、子实体产量与形态等农艺性状，确定了菌丝体的最佳成熟时间点。在此基础上，研究人员对来自最佳成熟时间点前后的菌丝样品，同时进行了转录组学（分析基因表达）和代谢组学（分析小分子代谢物）的高通量测序与分析。通过生物信息学方法，他们筛选了不同时间点间的差异表达基因(DEGs)和差异积累代谢物(DAMs)，并进行了功能富集分析，以揭示在成熟过程中被激活或抑制的生物学通路。此外，研究还结合两种方法（MⅠ：基因与代谢物共表达分析；MⅡ：差异基因与差异代谢物共同富集的通路分析）对转录组和代谢组数据进行了整合分析，以挖掘更深层次的调控网络。同时，为了验证组学分析的发现，研究检测了与ROS代谢和纤维素降解相关的关键生化指标，如过氧化氢(H₂O₂)含量以及内切-1,4-β-葡聚糖酶(EC 3.2.1.4)和纤维二糖水解酶(EC 3.2.1.91)的活性。所有分析的样本均来自上海成英农业发展有限公司提供的KS11菌棒。

3.1. 八个处理阶段木屑菌棒的动态转色和出菇特性

研究发现，菌棒表面在光照培养下从D35开始形成褐色菌膜，到D65时完全覆盖。出菇率从D45到D65持续上升，但未达到100%。从D75开始，所有菌棒均能出菇，但D75菌棒出菇不均匀，且单棒平均出菇数显著低于D85、D95和D105。在D85至D105中，D95在大多数性状上表现最佳，包括单棒产量、单菇重量、菌盖直径、菌柄直径和长度等。因此，D95被确定为KS11菌丝体完全成熟的最佳时间点。

3.2. 八个处理阶段菌丝体的动态转录组特征

转录组分析揭示了在七个相邻时间点比较对中，差异表达基因显著富集于三大类功能：氧化还原酶活性、蛋白质代谢和碳水化合物代谢。与氧化还原酶活性相关的术语在不同比较对中各不相同，涉及锰过氧化物酶(MnP)、细胞色素P450 (P450)、单加氧酶、谷氧还蛋白(Grx)、硫氧还蛋白(Trx)和过氧化还原酶(Prdx)等多种基因的动态变化。蛋白质代谢相关术语在多个比较对中富集，涉及丝氨酸蛋白酶抑制剂、26S蛋白酶体（包括核心颗粒和调控颗粒相关基因）以及多种氨基酸（如缬氨酸、亮氨酸、异亮氨酸、色氨酸等）的代谢通路。碳水化合物代谢相关术语同样在多个比较对中富集，主要集中在纤维素分解过程（涉及碳水化合物结合模块CBM家族1和糖苷水解酶GH家族基因），以及果胶降解、淀粉和蔗糖代谢、糖酵解/糖异生、丙酮酸代谢等基础代谢途径。

3.3. 八个处理阶段菌丝体的动态代谢组特征

代谢组分析在所有七个比较对中鉴定出12个共同富集的KEGG通路术语，主要属于氨基酸代谢。其中鉴定出三个共同的差异代谢物（MP46、MN189和MN276），它们的代谢途径集中于核心碳代谢通路（糖酵解和三羧酸循环）。此外，不同比较对还独特富集了与其他氨基酸、碳水化合物代谢、辅因子和维生素以及能量代谢相关的通路术语。

3.4. 转录组和代谢组的整合分析

通过两种整合分析方法（MⅠ和MⅡ）发现，富集的KEGG通路主要涉及氨基酸、其他氨基酸、碳水化合物以及辅因子和维生素的代谢。其中，四个通路被两种方法共同鉴定出：D55 vs. D65中的β-丙氨酸代谢；D65 vs. D75和D85 vs. D95中的丙酮酸代谢；D85 vs. D95中的柠檬酸循环(TCA循环)；D95 vs. D105中的丙氨酸、天冬氨酸和谷氨酸代谢。研究重点构建了两个在多个比较对中富集的通路（map00250 丙氨酸、天冬氨酸和谷氨酸代谢；map00620 丙酮酸代谢）的网络图，发现其中一些基因（如LENED_003951天冬酰胺合成酶、LENED_005767谷氨酸脱羧酶等）的表达水平在所有检测的时间点比较对中都发生了动态变化。

3.5. 八个处理阶段生化指标检测的动态特征

H₂O₂含量、内切-1,4-β-葡聚糖酶活性和纤维二糖水解酶活性在八个处理阶段均呈现动态变化。H₂O₂含量在D55和D95达到峰值，而两种纤维素酶活性均在D95达到最高值。

研究的讨论部分构建了一个关于香菇菌丝体生理成熟的潜在整合工作模型。该模型的核心是一个涉及活性氧(ROS)代谢、氨基酸稳态和碳通量的相互作用网络。ROS动态可能充当信号枢纽：在光照诱导下，菌丝体经历活跃的ROS代谢。细胞外H₂O₂可能由GMC氧化还原酶家族产生，并被锰过氧化物酶(MnP)等用于木质素降解，从而暴露出更多植物多糖供菌丝利用。细胞内ROS则可能在线粒体和内质网中生成，并受到谷胱甘肽、Grx、Trx和Prdx等抗氧化系统的清除。ROS水平的波动（如H₂O₂在D55和D95的峰值）可能作为关键信号。氨基酸代谢充当信号-代谢界面：未被清除的ROS可能氧化蛋白质，这些受损蛋白质通过泛素-蛋白酶体系统（26S蛋白酶体）被降解回收氨基酸。代谢组数据显示，在D85左右，多种氨基酸含量处于低点，提示菌丝可能处于一种氨基酸匮乏状态，而氮饥饿已被证实是子实体诱导的重要因素之一。特别是丙氨酸、天冬氨酸和谷氨酸的代谢(map00250)可能扮演关键角色，它们连接着氮资源调度与核心碳代谢。碳通量为发育提供能量和储备：菌丝体通过分泌纤维素酶等（如CBM1、GH家族基因）降解培养基中的纤维素，获取碳源。碳通量一方面通过糖酵解、丙酮酸代谢和TCA循环产生能量(ATP)和碳骨架；另一方面，碳可能被转化为主要的可溶性储存和运输形式——甘露醇。代谢组数据显示甘露醇含量在D35较高，随后下降，并在D95前又回升至峰值，暗示其作为碳储备和潜在抗氧化剂的作用。

综上所述，这项研究首次通过系统的多组学分析，揭示了香菇菌丝体在转色期生理成熟过程中，ROS信号、氨基酸循环和碳通量动员三者之间存在的时序协同与潜在互作关系。研究所提出的整合模型为理解菌丝体从营养生长向生殖生长转变的复杂分子机制提供了一个可验证的框架。其重要意义在于，不仅深化了我们对香菇生物学基础的认识，还为栽培实践提供了潜在的优化靶点。例如，通过调节氧气供应来干预ROS信号，或通过优化培养基配方来调控碳/氮通量，都可能成为优化菌丝成熟时机、提升整体生产效能的实用杠杆。该研究标志着我们在揭开香菇“转色”奥秘的道路上迈出了关键一步。