Ganoderma tsugae 发育阶段多糖积累及相关代谢重编程的多组学分析

《Frontiers in Microbiology》：Multi-omics analysis of polysaccharide accumulation and associated metabolic reprogramming across developmental stages of Ganoderma tsugae

　　

引言：多糖是Ganoderma tsugae子实体中重要的生物活性成分，然而其积累的发育动态及相关分子过程仍知之甚少。本研究旨在阐明在橡木仿生栽培条件下，子实体发育过程中与多糖相关的代谢重塑。方法：研究人员检测了五个发育阶段，包括原基（T1）、未成熟子实体（T2）、担孢子释放前（T3）、担孢子释放期（T4）和担孢子释放末期（T5）。测定了粗多糖含量，并整合了转录组学和非靶向代谢组学分析，以识别与多糖积累相关的关键通路、基因和代谢物。进一步通过qRT-PCR（定量逆转录聚合酶链式反应）验证了选定基因。结果：粗多糖含量从T1到T3逐渐增加，在T3达到最大值（0.90 g/100 g），随后下降。共鉴定出5,676个差异表达基因和911个差异积累代谢物，在转录组和代谢组水平均呈现明显的阶段依赖性分离。整合的KEGG（京都基因与基因组百科全书）分析表明，淀粉与蔗糖代谢是与多糖积累相关的核心通路，而半乳糖代谢作为补充模块，参与前体供应和结构多样化。FKS1（β-1,3-葡聚糖合酶相关基因）和多个糖苷水解酶相关基因的动态表达，以及糖磷酸酯和双糖代谢物的变化，提示底物降解、前体供应、聚糖延伸和碳再分配存在协调调控。qRT-PCR结果与RNA-seq数据一致。讨论：这些发现表明，G. tsugae子实体中多糖积累是一个发育协调的过程，而非孤立的单个生物合成步骤的结果。T3阶段似乎是多糖积累和碳代谢重编程的关键转折点。本研究为理解G. tsugae发育过程中多糖积累提供了分子和代谢框架，并为未来功能验证提供了潜在的候选基因和通路。

**论文解读文章**

**研究背景与问题**
Ganoderma tsugae Murrill 是一种珍稀药用真菌，主要分布于中国长白山地区，富含萜类、多糖、甾醇、腺苷等活性代谢物，具有抗氧化、抗肿瘤、增强记忆、止咳平喘等药理活性。其中，灵芝多糖被认为是该属最重要的功能成分之一，其生物活性与结构特征（如单糖组成、分子量、糖苷键类型、分支度及链构象）密切相关。然而，目前对G. tsugae多糖积累的发育动态及分子调控机制认识不足。尽管已有研究关注灵芝多糖的提取、结构解析及生物合成，但缺乏整合多组学方法对子实体发育全过程中碳源转化、前体供应、聚糖装配、细胞壁重塑及胞外输出等协同过程的系统解析。为此，研究人员以橡木仿生栽培方式，选取G. tsugae子实体发育的五个关键阶段（原基T1、未成熟子实体T2、担孢子释放前T3、担孢子释放期T4、担孢子释放末期T5），通过定量测定粗多糖含量，并整合转录组测序（RNA-seq）与非靶向液相色谱-串联质谱（LC-MS/MS）代谢组学分析，建立代谢组-转录组联合框架，揭示多糖生物合成及代谢重塑的分子基础。该论文发表在《Frontiers in Microbiology》。

**主要技术方法**
研究人员采用橡木仿生栽培体系，样本来源于中国吉林省长白县长白山志远生态农业有限公司栽培基地。核心方法包括：（1）苯酚-硫酸法测定粗多糖含量；（2）RNA提取、文库构建及Illumina测序获得转录组数据，用DESeq2鉴定差异表达基因（DEGs，|log₂ FC| ≥ 1, FDR ≤ 0.05）；（3）基于超高效液相色谱-四极杆静电场轨道阱质谱（UHPLC-Q Exactive HF-X）的非靶向代谢组学分析，用XCMS进行峰提取、对齐及定量，差异积累代谢物（DAMs）筛选标准为FDR ≤ 0.05、|log₂ FC| ≥ 1 且 VIP ≥ 1；（4）O₂PLS模型整合转录组-代谢组数据，Spearman相关网络分析；（5）qRT-PCR验证10个候选基因（以18S为内参，2^?ΔΔCt法计算相对表达）。

**研究结果**

**3.1 粗多糖含量的动态变化**
粗多糖含量从T1到T3逐渐升高，T3达到最大值（0.90 g/100 g），显著高于其他阶段（P < 0.01）；T4次之（0.84 g/100 g），T1（0.54 g/100 g）和T5（0.58 g/100 g）最低且无显著差异。结果表明，多糖积累在发育中期更为活跃。

**3.2 转录组与代谢组数据质量评估**
主成分分析（PCA）显示，转录组中T1与T2相近，T3、T4、T5另一簇，且T3与T5在PC1上接近；代谢组PCA中PC1和PC2解释总变异的46%，QC样本紧密聚集，表明分析稳定性良好。正交偏最小二乘判别分析（OPLS-DA）及200次置换检验证实模型稳健无过拟合。

**3.3 差异表达基因（DEGs）与差异积累代谢物（DAMs）分析**
共鉴定5,676个DEGs。相邻阶段比较中，T2 vs T3 DEGs最多（3,015个），T3 vs T4（588个）和T4 vs T5（339个）较少；交叉比较T5 vs T1 DEGs最多（3,675个）。UpSet分析显示，DEGs既有阶段特异性，也有多阶段共享（最大共享集909个）。共筛选出911个DAMs，化学分类以脂类及类脂分子（27.0%）、有机酸及衍生物（25.2%）、有机杂环化合物（15.4%）为主。层次聚类表明T1与T2代谢相似，T3与T4为一簇，T5独立。火山图显示T2 vs T3中下调代谢物多于上调，而T3 vs T4和T4 vs T5中上调多于下调。

**3.4 多糖生物合成相关DEGs与DAMs分析**
对KEGG数据库中淀粉与蔗糖代谢（ko00500）和半乳糖代谢（ko00052）两个通路分析。在基因层面，FKS1（β-1,3-葡聚糖合酶）及多个外切葡聚糖酶基因（CBH1、BGL1A、exgD）在后期表达升高；海藻糖合成基因（tpsA、TPS2）在早期高表达。代谢物层面，ko00500中大多数糖及磷酸化中间体（如麦芽糖、蔗糖、葡萄糖-1-磷酸、葡萄糖-6-磷酸）在早期丰富，后期下降；而海藻糖和异麦芽糖在后期积累升高。ko00052中，N-乙酰-D-半乳糖胺、半乳糖醇、二羟基丙酮磷酸等早期高，后期下降；棉子糖家族寡糖在T2阶段特异性增加。联合通路映射显示，ko00500中多个糖苷水解酶相关基因及糖磷酸酯代谢物的协调变化，表明G. tsugae发育过程中增强了对结构多糖和储存多糖的分解能力。

**3.5 qRT-PCR验证**
从ko00500和ko00052中选取10个基因（FKS1、exgD、celD、BGL1A、LKA1、amyB、agl1、SPAC2E11.16c、galE、GAL7）进行qRT-PCR验证，结果与RNA-seq表达趋势一致，证实了转录组数据的可靠性。

**3.6 代谢组与转录组相关性分析**
O₂PLS分析显示，基因与代谢物在第一联合轴正向区域中，关键候选基因（如aglB、CBH1、exgD、ARB_05372）与蔗糖、麦芽糖、乳糖、葡萄糖-1-磷酸紧密聚类；负向区域中，aglI、BGLU12、celD、bglL与异麦芽糖共分布。Spearman相关网络呈现正负相关的基因-代谢物对。

**3.7 多糖生物合成功能模块的整合分析**
将两个通路整合为七功能区域：红区（β-葡聚糖合成主链：FKS1编码的β-1,3-葡聚糖合酶及exo-1,3-β-葡萄糖苷酶）、绿区（纤维素降解：CBH1早期高表达，egl2和celD后期高表达）、粉区（淀粉/糖原降解：四个DEGs后期表达升高）、蓝区（海藻糖代谢：合成基因早期高，降解产物后期积累）、紫区（蔗糖降解：三个DEGs中后期高表达）、橙区（半乳糖转化：galK、GAL7、galE早期高表达）、黄区（水苏糖降解：agl1编码α-半乳糖苷酶）。整合表明，β-葡聚糖生物合成涉及35个DEGs、8个DAMs和20个酶。

**讨论与结论**
讨论指出，G. tsugae粗多糖含量呈现“先升后降”的发育模式，T3为峰值。这一趋势与碳分配策略相关：T1阶段碳资源主要用于菌丝生长和底物分解；T2–T3阶段碳流向细胞壁和胞外多糖合成；T4–T5阶段多糖部分降解以支持孢子形成和组织维持。机械论上，淀粉与蔗糖代谢（ko00500）构成核心碳通量路径，FKS1与多个糖苷水解酶基因的协同表达及糖磷酸酯、双糖的变化表明多糖积累是前体供应、聚糖延伸和细胞壁重塑的耦合过程。半乳糖代谢（ko00052）作为补充模块，在早期提供UDP-糖前体，且与底物中棉子糖家族糖类相关。海藻糖的后期积累反映了能量稳态和可快速动员碳库的建立。结论部分翻译如下：
本研究表明，G. tsugae子实体中粗多糖积累具有显著的阶段特异性，在担孢子释放前达到最高水平，随后在担孢子释放和成熟后期略有下降。整合多组学分析进一步指出，这一模式与转录重编程和代谢重塑密切相关。淀粉与蔗糖代谢（ko00500）是多糖生物合成的核心碳通量通路，而半乳糖代谢（ko00052）作为碳供应和结构多样化的补充模块。FKS1与多个糖苷水解酶相关基因的协调变化，以及糖磷酸酯和双糖的变化，表明G. tsugae多糖积累是一个涉及底物利用、前体生成、聚糖延伸和后期重塑的耦合过程。从应用角度看，这些发现将T3确定为阶段定向采收、质量评估及后续构效关系引导利用的有利发育窗口。然而，由于本研究聚焦于粗多糖含量及多组学关联，T3来源多糖的具体结构和生物活性仍有待确定。未来需结合多糖纯化、单糖组成分析、分子量测定、糖苷键分析、傅里叶变换红外光谱（FT-IR）/核磁共振（NMR）表征及抗氧化或免疫调节活性测定，以验证T3积累多糖组分的功能意义。总体而言，这些发现为理解G. tsugae多糖形成的发育调控及未来精准栽培、底物优化和候选基因功能验证提供了重要基础。