通过比较基因组学与转录组学分析筛选与香菇（Lentinula edodes）菌盖缺失表型相关的候选基因

《Journal of Fungi》：Identification of Candidate Genes Associated with the Pileus-Deficient Phenotype in Lentinula edodes Through Comparative Genomic and Transcriptomic Analyses Bo-Min Seo, Che-Hwon Park, Sung-Chul Lee, Rae-Won Kang and Young-Jin Park

本研究旨在阐明香菇（Lentinula edodes）两个菌株（Le_L和Le_S）间表型差异的分子机制，这两个菌株表现出不同的子实体形态。虽然菌株间的表型变异已被广泛观察到，但转录调控和结构基因组变异对这些差异的相对贡献仍知之甚少。比较转录组分析鉴定出8541个差异表达基因（DEGs），揭示了两个菌株之间清晰的功能分化。Le_S中上调的基因主要富集于核糖体成分和翻译相关过程，表明其蛋白质合成活性增强。相反，Le_L中上调的基因与转运蛋白、转录因子和多种代谢途径相关，暗示了其具有更广泛的调控和生理功能。蛋白质-蛋白质相互作用（PPI）网络分析进一步突出了不同的调控架构：在Le_S中，核糖体蛋白形成了高度互作的枢纽基因模块，而Le_L的枢纽基因功能多样，并包含多个主要易化子超家族（MFS）成员。对33个L. edodes菌株的同源基因分析表明，大多数枢纽基因是保守的，表明它们作为核心遗传组分发挥作用。尽管存在广泛的基因组水平变异（包括7931个SNPs和1149个INDELs），但枢纽基因内的序列变异有限，并且等位基因特异性表达（ASE）分析未显示显著的等位基因失衡。相比之下，存在-缺失变异（PAV）分析发现了影响MFS转运蛋白基因的结构差异，这些基因在Le_S中缺失，但在Le_L中存在并上调。总之，这些发现表明，在驱动表型分化方面，结构基因组变异（特别是涉及转运蛋白基因的变异）可能比序列水平变异起着更重要的作用。本研究为香菇菌株特异性性状的遗传基础提供了新见解，并强调了整合多水平基因组分析的重要性。

香菇菌盖缺失表型的遗传与分子机制解析

研究背景、问题与目的

香菇（Lentinula edodes）是全球范围内广泛栽培的食用菌之一，其子实体的正常形成，特别是菌盖（Pileus）的发育，直接关系到孢子的产生和商业价值。然而，在栽培过程中常观察到子实体形态异常，例如菌盖发育不全或缺失，这反映了发育调控的潜在紊乱。尽管先前研究已利用转录组学等方法探究了香菇子实体发育的分子基础，识别了与细胞壁重塑、碳水化合物代谢、信号转导等相关的基因表达变化，但对于驱动菌株间表型差异的遗传因素，尤其是转录调控、序列变异与结构基因组变异各自的作用，仍缺乏系统性的整合理解。

为了解决上述问题，研究人员选取了两个具有鲜明表型对比的香菇菌株：CBMLE44（Le_L）可产生具有发育良好菌盖的正常子实体，而KMCC05708（Le_S）则表现出菌盖显著减少或缺失的异常表型。本研究旨在通过整合多组学分析，探究这两个菌株表型差异的分子机制。具体目标包括：（1）表征两菌株间的转录组与功能差异；（2）评估枢纽基因在多个菌株间的保守性与变异性；（3）确定序列水平变异和结构基因组变异对表型分化的相对贡献。该研究为理解香菇菌株特异性性状的遗传基础提供了新的见解，相关成果发表在《Journal of Fungi》期刊上。

关键技术方法概述

本研究采用了多种生物信息学与基因组学技术进行整合分析。核心方法包括：1. 对Le_L和Le_S菌株的子实体组织进行RNA测序（RNA-seq），并利用DESeq2进行差异表达基因（DEG）分析；2. 对差异表达基因进行基因本体（GO）和京都基因与基因组百科全书（KEGG）通路富集分析，以明确功能分化；3. 基于STRING数据库构建蛋白质-蛋白质相互作用（PPI）网络，并使用Cytoscape的CytoHubba和MCODE插件识别枢纽基因及模块；4. 使用OrthoFinder对包括Le_L和Le_S在内的33个香菇菌株进行同源基因分析，评估枢纽基因的保守性；5. 基于全基因组测序数据，进行单核苷酸多态性（SNP）和插入缺失（INDEL）变异检测，并结合RNA-seq数据对枢纽基因内的变异位点进行等位基因特异性表达（ASE）分析；6. 开展存在-缺失变异（PAV）分析，以鉴定菌株间的结构基因组差异。

研究结果分析

3.1. Le_L与Le_S菌株的形态差异

研究人员观察到Le_L和Le_S菌株在子实体形态上存在显著差异。具体而言，Le_L菌株可产生具有发育良好菌盖的正常子实体，而Le_S菌株则表现出菌盖显著减少或几乎缺失的异常表型。这种差异表明Le_S的菌盖形成过程存在发育程序的破坏。

3.2. 转录组分析与差异表达基因的功能富集

转录组比较分析共鉴定出8541个差异表达基因。功能富集分析显示，Le_S中上调的基因主要富集于“核糖体结构成分”、“翻译”等GO术语，以及核糖体相关的KEGG通路，表明其蛋白质合成活性增强。相反，Le_L中上调的基因则显著富集于“跨膜转运”、“跨膜转运蛋白活性”、“DNA结合转录因子活性（RNA聚合酶II特异性）”以及葡聚糖代谢等功能类别。这表明Le_L与更广泛的运输、代谢和转录调控功能相关，而Le_S的转录程序则更集中于核糖体相关过程。

3.3. 蛋白质相互作用网络分析与枢纽基因鉴定

蛋白质-蛋白质相互作用网络分析进一步揭示了两菌株调控架构的差异。Le_L上调基因的网络中，枢纽基因功能多样，包含多个主要易化子超家族（MFS）转运蛋白、糖转运蛋白和转录因子。而Le_S上调基因的网络则高度紧密，核心枢纽基因几乎全部由核糖体蛋白组成。这表明Le_L的发育程序围绕转运和调控功能组织，而Le_S则表现出对翻译机器的高度协调上调。

3.4. 枢纽基因的表达模式

对鉴定出的枢纽基因进行表达模式验证，热图、火山图和箱线图均证实了这些基因在两菌株中存在显著的差异表达。具体而言，Le_L相关的枢纽基因（如MFS转运蛋白基因）在Le_L中表达更高，而Le_S相关的枢纽基因（如核糖体蛋白基因）在Le_S中表达更高。这进一步支持了枢纽基因是与表型差异相关的生物学上有意义的候选因子。

3.5. 枢纽基因在香菇菌株间的保守性与表达模式

同源基因分析表明，大多数枢纽基因在分析的33个香菇菌株中广泛存在，属于该物种的核心基因库。尽管存在广泛的基因组变异，但枢纽基因的拷贝数在两菌株间没有显著差异。这表明表型分化不太可能由大多数枢纽基因的大规模基因存在-缺失或拷贝数变异来解释，而更可能源于基因调控和网络水平表达动态的差异。

3.6. 全基因组变异分析与枢纽基因内变异的鉴定

全基因组变异分析在Le_S中鉴定了7931个SNP_s和1149个INDEL_s。然而，只有极少数的枢纽基因内部含有变异，且这些变异多为内含子替换或同义外显子替换，不太可能直接破坏蛋白质功能。这表明序列水平的变异本身不太可能是造成观察到的转录差异的主要机制。

3.7. 枢纽基因内变异的等位基因特异性表达分析

对枢纽基因内SNP_s位点的等位基因特异性表达分析显示，尽管Le_S在一些位点呈杂合状态，但等位基因表达比例接近0.5，未观察到一致的等位基因表达偏倚。因此，观察到的转录差异不太可能主要由这些SNP_s位点的顺式作用等位基因特异性效应驱动。

3.8. 枢纽基因的存在-缺失变异分析

存在-缺失变异分析揭示了一个关键发现：部分与膜转运相关的基因，特别是MFS转运蛋白基因，在Le_S菌株中缺失，而在Le_L及其他分析菌株中存在。这与转录组和网络分析中MFS基因在Le_L中上调并作为枢纽基因的结果相一致。因此，转运蛋白基因的结构变异被认为是导致表型分化的一个关键候选因素。

讨论与结论总结

讨论部分综合了上述发现，认为两菌株的表型差异并非主要由大规模的基因含量或编码序列变异驱动，而是更可能源于基因调控和网络组织的协同变化。其中，转运蛋白基因（尤其是MFS家族）的结构变异是一个关键贡献因素。这些基因在营养摄取、代谢物分配和发育稳态中起着关键作用，其缺失可能破坏正常菌盖发育所需的生理条件。研究也指出，尽管基于生物信息学分析提出了强有力的候选基因，但仍需未来的功能实验（如qPCR或功能测定）来验证其在菌盖形态发生中的具体作用。

结论翻译：总之，本研究对香菇表型分化的转录调控、网络结构和基因组变异进行了整合分析。研究结果表明，调控和转运蛋白相关基因网络的差异，而非单纯的序列水平变异，在决定菌盖形态中起着核心作用。重要的是，本研究中鉴定出的转运蛋白相关枢纽基因和结构变异位点，是未来功能验证的有希望候选基因。从应用角度看，这些发现为改善香菇的菌株选育、育种策略和栽培稳定性提供了潜在的框架。整合转录组学、网络和结构变异分析的方法，也可能有助于识别其他蘑菇物种发育性状的关键决定因素。