《World Journal of Microbiology and Biotechnology》：Multi-omics analysis reveals Protein Kinase A-associated regulatory remodeling during adaptation of Trichoderma reesei to lignocellulosic substrate
编辑推荐：
丝状真菌里氏木霉（*Trichoderma reesei*）是全纤维素降解酶的主要工业来源，其对复杂碳源的响应受营养感应机制调控，包括环磷酸腺苷（cAMP）–蛋白激酶A（PKA）信号通路。本研究整合转录组学、定量蛋白质组学和磷酸化蛋白质组学，分析了在葡萄糖或甘
丝状真菌里氏木霉（*Trichoderma reesei*）是全纤维素降解酶的主要工业来源，其对复杂碳源的响应受营养感应机制调控，包括环磷酸腺苷（cAMP）–蛋白激酶A（PKA）信号通路。本研究整合转录组学、定量蛋白质组学和磷酸化蛋白质组学，分析了在葡萄糖或甘蔗渣条件下培养的亲本菌株QM9414和Δpkac1菌株中PKAc1相关响应。*pkac1*缺失与生长相关表型改变及部分生物质解聚酶胞外活性降低相关。多组学分析揭示了条件依赖性变化，涉及碳水化合物活性酶（CAZymes）基因和蛋白质亚群、营养转运蛋白、应激相关蛋白及调控因子。磷酸化蛋白质组学鉴定了与*pkac1*缺失相关的磷酸化状态变化，包括富含PKA共有基序位点的磷酸化水平降低。采用计算机模拟肽段对接对候选PKAc1关联底物进行优先级排序以供后续验证，其中源自Sec7的肽段相较对照肽段表现出更优的对接行为。综上，这些数据支持一个工作模型，即PKAc1在里氏木霉适应甘蔗渣过程中参与调控重塑及磷酸化状态重塑，并影响所选CAZyme相关输出的幅度和/或时序。
**多组学解析PKAc1在里氏木霉木质纤维素适应中的调控重塑**
**1 研究背景与科学问题**
木质纤维素生物燃料是化石燃料的可持续替代方案，但植物生物质的顽固性使高效糖化仍面临技术挑战。第二代（纤维素）乙醇生产需要强大的酶系统将纤维素和半纤维素分解为可发酵糖。丝状真菌里氏木霉（*Trichoderma reesei*）因其数十年来经菌株改良获得的卓越分泌能力，是这些纤维素酶和半纤维素酶的主要工业来源。然而，里氏木霉木质纤维素酶的天然生产受到营养可利用性、分解代谢物阻遏机制及光照等环境信号的严格调控，尤其Cre1介导的阻遏作用甚至可调节对槐糖等强诱导物的响应。在简单碳源（如葡萄糖）存在时，高效酶表达被抑制，仅当真菌检测到纤维素等复杂多糖时才被诱导，引发与葡萄糖分解代谢不同的代谢重构。该调控由复杂转录因子网络（如主调控因子Xyr1、阻遏因子Cre1）以及感应碳源的信号通路共同介导，其中cAMP–PKA通路对分泌至关重要，MAPK通路在代谢和纤维素酶调控中发挥双重作用，液泡钙转运系统协调生长信号与产物合成。
PKA是真菌营养响应信号的主要调控因子，位于葡萄糖感应系统下游，根据营养状态协调生长、发育和代谢。在模型酵母和丝状真菌中，PKA活性影响从孢子形成、应激响应到水解酶表达等多种过程。此前在粗糙脉孢菌和烟曲霉中的研究表明PKA突变体表现出生长减缓及发育异常。在里氏木霉中，间接证据提示PKA参与纤维素酶调控和次级代谢；cAMP升高可激活钙信号并增强纤维素酶基因表达；磷酸化蛋白质组分析提示PKA可能翻译后调控分泌酶；敲除PKA催化亚基编码基因*pkac1*（QM6a背景）导致纤维素酶产生受损和生长缺陷。尽管PKA/PKAc1先前已被认为参与里氏木霉的生长、纤维素酶产生、光依赖性调控和次级代谢，但其在转录组、蛋白质组和磷酸化蛋白质组层面的效应尚未在整合实验框架中系统检验。
**2 研究方法概述**
研究人员选用工业相关的木质纤维素底物甘蔗渣（由纤维素、半纤维素、木质素及次要组分构成），以葡萄糖作为碳阻遏对照条件。以里氏木霉QM9414（ATCC26921）为亲本菌株，使用此前报道的*pkac1*缺失突变体Δpkac1，通过诊断PCR和RT-qPCR确认缺失状态。在含2%葡萄糖或1%甘蔗渣的液体Mandels培养基中培养24小时。主要技术方法包括：①表型与酶活分析（菌落形态、径向生长、生物量干重，以及纤维素内切酶活性（CMCase）、滤纸酶活（FPase）、β-葡萄糖苷酶、β-木糖苷酶、纤维二糖水解酶和木聚糖酶活性测定）；②转录组测序（RNA-seq，Illumina平台150 bp双端测序，HISAT2比对至里氏木霉参考基因组v2.0，DESeq2鉴定差异表达基因）；③定量蛋白质组学（LC-MS/MS，Orbitrap Exploris 240，MaxQuant标记-free定量，Perseus统计检验）；④磷酸化蛋白质组学（TiO?亲和富集磷酸肽，MaxQuant磷酸化位点定位，基于强度的标记-free定量）；⑤生物信息学分析（GO和KEGG富集、PKA共有基序扫描、HPEPDOCK肽段-蛋白质对接）。所有组学测量均使用每种菌株和碳源条件下三个生物学重复，在标准化时间点（24小时）采集样本。蛋白质组学数据已提交至ProteomeXchange（PXD072673），RNA-seq数据已提交至NCBI SRA（PRJNA1398409）。
**3 研究结果**
**3.1 Δpkac1在甘蔗渣上表现出生长表型改变和胞外酶活性降低**
PKA催化亚基缺失与里氏木霉生长、菌落形态和色素沉着改变相关。Δpkac1菌株径向生长略减，形成更小更致密的菌落，并丧失QM9414特有的绿色色素。在木糖、纤维二糖、乳糖和甘油为唯一碳源时Δpkac1径向生长显著降低，但葡萄糖和CMC条件下无统计差异。在甘蔗渣条件下，Δpkac1的纤维二糖水解酶和FPase胞外活性按培养物体积计算分别显著降低约61%（p=0.009）和36%（p=0.0003），木聚糖酶活性也降低；而CMCase、β-葡萄糖苷酶和β-木糖苷酶无显著差异。由于24小时干生物量在菌株间相当（QM9414: 0.051±0.004 g/L；Δpkac1: 0.050±0.004 g/L；p=0.67），这些降低也对应于单位生物量活性的降低。但甘蔗渣上的酶活性仍高于葡萄糖对照（<0.1 U/g）。
**3.2 转录组分析揭示Δpkac1中CAZyme相关表达模式改变**
主成分分析显示生物重复聚类良好，碳源（PC1，76%方差）和基因型（PC2，12%方差）是转录变异的主要驱动因素。差异表达分析表明*pkac1*缺失在两种条件下均引发广泛转录变化。在甘蔗渣条件下，许多正常诱导的CAZyme基因在Δpkac1中表达水平降低，包括木聚糖酶基因*xyn1*（log?FC=-1.28）和转录因子*ace3*（log?FC=-1.22）。主要纤维素酶基因*cel7a/cbh1*和*cel6a/cbh2*在Δpkac1中转录水平较低但未达统计显著，提示诱导减弱而非丧失。关键调控因子*xyr1*有数值下降但未达显著，而阻遏因子*cre1*转录略升，细胞糊精转运蛋白*stp1*显著下调（log?FC=-1.11）。Δpkac1中上调最高的基因富集于应激响应因子和替代营养转运蛋白，而非纤维素分解系统成分。这些结果支持PKAc1贡献于部分半纤维素分解响应的转录幅度。
**3.3 蛋白质组学揭示部分CAZyme蛋白丰度降低和应激相关通路变化**
在甘蔗渣条件下，Δpkac1中若干CAZyme相关蛋白丰度改变，但Cel7A接近中性；菌株内的甘蔗渣对葡萄糖比较显示两菌株均强烈诱导Cel7A和其他CAZyme蛋白，说明Δpkac1保留了木质纤维素响应性CAZyme程序，仅特定蛋白丰度模式改变。GO富集分析显示，Δpkac1在甘蔗渣条件下高丰度蛋白主要映射到三羧酸（TCA）循环、糖酵解和氨基酸代谢等核心初生及能量代谢通路，以及线粒体呼吸体和质子转运ATP合酶复合物等细胞组分；而亲本菌株高丰度蛋白富集于翻译起始、孔复合物和内质网动力学。葡萄糖条件下，Δpkac1富集到D-苏式-醛糖1-脱氢酶和CoA转移酶活性，且两种碳源条件下均有过氧化物还原酶活性等类别过度表达；而核糖体结构成分在葡萄糖条件下显著耗竭，提示突变体翻译准备受损。
**3.4 磷酸化蛋白质组学揭示PKA相关信号网络**
在葡萄糖条件下，Δpkac1中检出的磷酸肽数量（702）少于亲本（1211），且低磷酸化位点富集PKA共有基序（R-R/K-x-S/T）。GO富集显示差异磷酸化蛋白涉及染色质组织和细胞运输。在甘蔗渣条件下，多数位点仍显示低磷酸化，但部分位点显示高磷酸化；GO富集涉及信号转导、转录调控和运输。具体候选调控蛋白中，转录因子VIB1（一个锌簇蛋白，参与半纤维素酶基因诱导）和CLP1（一个假定的PHD指蛋白，与Xyr1协同激活纤维素酶基因）在Δpkac1中磷酸化显著降低；糖转运蛋白STP1存在差异磷酸化位点；高尔基蛋白Sec7（ARF GTP酶鸟嘌呤核苷酸交换因子）在Δpkac1中低磷酸化，与突变体胞外酶积累降低一致。
**3.5 对接分析优先排序候选PKAc1底物**
对含PKA共有基序且在Δpkac1中低磷酸化的9残基肽段进行HPEPDOCK对接，以PKAr1衍生肽为对照。排名靠前的候选肽段来自Tr76037、Tr120819、Tr120411和Sec7/Tr43877。其中Tr76037和Tr120819的ΔScore最接近PKAr1对照，Tr120411（序列HRRKSSTHH）是甘蔗渣条件下排名最高的候选。这些结果提供了候选PKAc1关联底物的优先验证列表，但无法确立直接激酶-底物关系。
**3.6 多组学整合揭示Δpkac1中转录-蛋白不一致性**
跨组学比较显示，转录水平变化与蛋白质或磷酸化位点丰度变化之间的决定系数极低（R2≤0.04），回归斜率平缓，大量特征点位于1:1对角线下方，表明许多基因的蛋白或磷酸化变化比转录变化更负或更正。在葡萄糖条件下，70个靶标在mRNA和总蛋白水平协调上调（Up/Up象限），显著多于甘蔗渣条件下的42个。而转录-磷酸化位点关联在甘蔗渣条件下有30个同步Up/Up靶标，葡萄糖条件下仅有5个Down/Down小簇。该结果支持PKAc1效应涉及转录、转录后、翻译后和分泌相关过程的多层调控。
**3.7 讨论与结论**
讨论指出，*pkac1*缺失并未完全丧失木质纤维素分解诱导，而是削弱和重新分配了特定CAZyme相关输出。残余响应并非由旁系同源基因*pkac2*的代偿性上调驱动（其表达未改变）。PKA可能作用于Xyr1上游或平行通路，通过影响GPCR（如PTH11家族）感应葡萄糖耗竭和替代碳源信号。磷酸化蛋白质组数据提示PKAc1可能通过磷酸化VIB1和CLP1调节转录因子活性，并通过Sec7相关高尔基运输影响分泌能力。突变体中特定Rab GTP酶和SNARE蛋白丰度增加可能反映囊泡运输的代偿变化。与先前在QM6a背景中的研究相比，本研究中较弱生长表型和保留的诱导能力可能源于遗传背景、培养基、培养方式等差异。结论部分指出：PKAc1有助于里氏木霉对木质纤维素底物的调控适应；在转录组、蛋白质组和磷酸化蛋白质组层面，*pkac1*缺失与甘蔗渣上生长相关表型改变、所选木质纤维素酶胞外活性降低以及富集PKA共有基序的广泛磷酸化状态变化相关；这些多组学模式支持一个工作模型，即PKAc1相关磷酸化可能有助于碳源感应、调控输出和分泌相关过程；数据集支持多种CAZyme相关响应诱导减弱，但生长阶段效应、转录调控、分泌相关过程和直接激酶-底物关系的相对贡献需要通过时间进程实验和定向验证来阐明。
