《Applied Microbiology and Biotechnology》：Transcriptomic insights into fruiting body malformations in Lentinula edodes
编辑推荐：
为应对香菇（Lentinula edodes）商业化栽培中频繁出现的子实体畸形问题，研究者对边缘凹陷、中心突起、孪生及不规则等四种畸形表型与正常样本进行了比较转录组学分析。研究揭示了表型特异性的转录特征，鉴定出跨表型一致的差异表达基因（DEGs），并通过功能富集与WGCNA阐明了细胞结构、代谢通路（如KEGG）失衡及环境应激响应在畸形发生中的作用，为定向育种和栽培管理提供了关键的分子靶点。
香菇，作为一种全球广泛栽培的食药用真菌，因其独特的风味和丰富的营养价值而备受青睐。然而，在商业化栽培过程中，子实体经常出现各种形态畸形，如菇盖边缘凹陷、中心异常突起、孪生融合或整体不规则皱缩等。这些畸形不仅降低了商品外观价值，更直接导致产量损失，对产业经济效益构成持续挑战。尽管环境因素（如温湿度波动）常被归咎于此，但在相同栽培条件下，畸形子实体与正常子实体常常同时出现，暗示着内在的分子调控紊乱可能是关键驱动因素。此前的研究多集中于子实体的发育过程或采后保鲜，对不同类型畸形之间的分子差异及其共同基础的系统性解析尚属空白。为深入阐明香菇子实体畸形的分子机制，并为精准预防和品种改良提供科学依据，一项针对四种明确畸形表型的比较转录组学研究得以开展。
该研究旨在通过高通量测序技术，系统比较四种畸形香菇子实体（边缘凹陷/H、中心突起/P、孪生/T、不规则/I）与其各自对应的正常对照（N1~N4）在相同栽培条件下的全基因组转录水平差异。研究最终揭示，香菇子实体畸形主要源于内在的转录组失调，而非病原体侵染。不同畸形表型展现出特异的转录特征，并共享一些关键的分子变化，这些发现为未来通过分子标记辅助育种和优化栽培管理策略来减少畸形发生率提供了潜在靶点和理论指导。相关研究成果已发表于《Applied Microbiology and Biotechnology》期刊。
本研究主要运用了几项关键技术方法。首先，研究人员于2024年5月在中国河北省阜平县的蘑菇农场，采集了四种畸形表型（H, P, T, I）及其对应的正常对照（N1~N4）的香菇（品种T2）成熟子实体的菌盖组织样本，每个类型设置三个生物学重复。随后，利用Illumina HiSeq平台对提取的总RNA进行高通量测序，获得转录组数据。数据分析方面，使用HISAT2软件将高质量测序片段比对到香菇参考基因组，并利用DESeq2软件包进行差异表达基因（DEGs）的鉴定。此外，还应用了基因本体（GO）和京都基因与基因组百科全书（KEGG）数据库对DEGs进行功能富集分析，并采用加权基因共表达网络分析（WGCNA）来识别与表型相关的基因共表达模块。
转录组测序、组装与注释
研究人员首先观察了畸形与正常香菇菌盖的形态差异，并通过切片排除了病原体感染的可能。随后，对24个样本（4组比较，各3重复）进行了RNA测序，获得了高质量的转录组数据。平均每个样本产出约6.47 Gb的clean reads，Q30值高于95%，其中约86.98%的reads能够比对到香菇参考基因组。共检测到9580个已知基因和1217个新基因的表达。这一部分工作为后续的差异分析奠定了坚实的数据基础。
差异表达基因的鉴定
通过比较分析，研究发现了不同畸形表型与正常对照之间存在显著差异的基因表达谱。其中，不规则（I vs N4）和孪生（T vs N3）表型表现出最严重的转录失调，分别有288和275个差异表达基因（DEGs）；中心突起（P vs N2）有185个DEGs；而边缘凹陷（H vs N1）的变化最小，仅61个DEGs。值得注意的是，两个基因（LENED_000547和 LENED_004546）在所有四种畸形类型的比较中都被一致地鉴定为差异表达。序列分析表明，尽管这两个蛋白共享多个保守序列基序，但它们在长度、理化性质、结构域组成和亚细胞定位预测上存在显著差异，前者可能是一个核定位的调控蛋白，而后者则可能是一个膜锚定的结构蛋白，暗示了它们在畸形发生中可能扮演不同的功能角色。
差异表达基因的功能注释与富集分析
为了解DEGs的生物学功能，研究进行了GO和KEGG富集分析。结果显示，在所有畸形类型中，与结构分子活性、核苷酸结合和转移酶活性相关的基因普遍下调，表明细胞完整性和基础代谢可能受损。相反，与碳水化合物代谢、跨膜运输和水解酶活性相关的基因则普遍上调，提示养分分配和细胞壁/膜动力学发生了改变。KEGG通路分析进一步指出，果糖/甘露糖代谢、谷胱甘肽代谢和甘油脂质通路在所有畸形类型中均受到干扰，揭示了代谢重编程和氧化应激反应是畸形的共同分子特征。此外，不同表型也显示出独特的富集模式，例如H型与碳水化合物代谢相关，P型与脂肪酸代谢和膜功能相关，T型与聚糖分解和细胞壁组织相关，I型则与应激响应和泛素化过程相关。
识别与四种表型异常相关的基因模块
通过WGCNA，研究构建了基因共表达网络，并识别出14个不同的基因模块。分析发现，绿色模块与中心突起（P）表型呈强正相关，而红色模块与边缘凹陷（H）表型呈强正相关。在关键的绿色模块中，多个应激响应基因（如谷胱甘肽转移酶、热激蛋白、DDR48、醛脱氢酶、胚胎发育晚期丰富蛋白（LEA）和酪氨酸酶）显著上调。这进一步支持了环境胁迫（如极端温度、低湿度或基质毒性）可能加剧了香菇的转录反应，从而诱发或加剧了子实体畸形。
讨论与结论
讨论部分对研究结果进行了深入整合与阐释。研究指出，在相同栽培条件下畸形与正常子实体的共存，强有力地表明内在的转录组失调是畸形的主要驱动因素，这与其他担子菌中的发现一致。研究提出了一个预测模型，即异常菌盖的形成是保守的代谢与氧化应激失衡，与表型特异的碳水化合物、脂质、细胞壁代谢、钙信号传导和应激防御通路紊乱相互作用的结果。其中，细胞壁动态重塑、钙信号传导（涉及钠/钙交换蛋白LeNCX和IQ基序钙调蛋白结合蛋白）以及强烈的应激防御机制（如LEA蛋白、热激蛋白的诱导）的激活，被认为是导致形态发生程序异常的关键环节。两个在所有畸形中一致的DEGs（LENED_000547和 LENED_004546）因其独特的结构特征被视为未来功能验证的核心候选基因。
综上所述，本研究通过对四种香菇子实体畸形表型的系统性比较转录组分析，初步揭示了畸形发生的分子基础是内在的转录调控紊乱。研究不仅刻画了不同畸形表型的特异性转录签名，也发现了跨表型共享的分子变化，包括代谢通路失调、应激响应激活和细胞结构相关基因的表达改变。这些发现构建了一个分子框架，将异常形态与特定的基因表达模式和生物学通路联系起来。尽管当前结果主要为相关性发现，但它们为后续通过基因敲除、过表达等实验进行因果验证提供了明确的方向和靶点。最终，这些分子见解有望应用于指导香菇的精准育种（如选育抗畸形品种）和优化栽培管理策略（如控制诱发应激的环境因子），从而减少经济损失，提升产业效益。研究也明确指出，需要未来在更广泛的遗传背景和栽培条件下进行验证，并通过功能实验将转录组关联转化为确凿的机制认知。
