《Frontiers in Cell Death》：Cross-stress transcriptomic integration in Marchantia polymorpha reveals stress-specific expression of regulated cell-death-associated genes
编辑推荐：
调控性细胞死亡（regulated cell death, RCD）参与植物发育、胁迫适应和免疫过程，但RCD相关基因在不同胁迫下是否表现出保守的转录响应仍不清楚。研究人员对公开可用的地钱（Marchantia polymorpha）RNA-seq数据集进行了
调控性细胞死亡（regulated cell death, RCD）参与植物发育、胁迫适应和免疫过程，但RCD相关基因在不同胁迫下是否表现出保守的转录响应仍不清楚。研究人员对公开可用的地钱（Marchantia polymorpha）RNA-seq数据集进行了比较性再分析，以探究RCD相关基因的重复性及胁迫特异性表达。研究人员统一处理了来自三项独立研究的36个RNA-seq样本，得到9组胁迫与对照的比较：8组为非生物胁迫（其中7组来自同一项研究，另有一项独立的高光实验）以及1组病原菌侵染。针对9组胁迫数据集中每一组与其各自对照相比，均鉴定了差异表达基因（differentially expressed genes, DEGs），随后检查了不同胁迫与对照比较组合之间的共享DEGs。研究人员使用此前在地钱中鉴定的14个RCD相关基因的精选列表，检验转录响应的保守性。此外，在数据集批次校正后，利用加权基因共表达网络分析（weighted gene co-expression network analysis, WGCNA）对整合表达数据进行检测。在所有9组比较共有的34个DEGs中，仅12个基因在每次比较中均以相同方向变化。这34个基因包含两个转录因子：一个AP2/ERF家族因子（Mp5g22160）和一个bHLH因子（Mp2g00890），它们在侵染和Schr?der高光处理中上调，但在Tan等人（Tan et al.）的所有7组比较中下调。在14个精选RCD相关基因中，10个在至少一组比较中差异表达，其中MpZOU2是唯一在多于一种胁迫中受调控的基因（在冷、黑暗和侵染中上调，在Tan高光中下调）。WGCNA鉴定了64个共表达模块，包括一个与侵染相关的模块，该模块富集蛋白质泛素化、蛋白酶体介导的降解、内质网相关降解、自噬和膜运输；以及一个与缺氮和渗透胁迫相关的模块，该模块富集氧化应激反应和细胞壁组织。综上，这些发现表明地钱在多种胁迫下并未表现出统一的转录RCD程序。相反，RCD相关基因似乎在更广泛的、依赖于特定胁迫背景的胁迫响应网络中被选择性招募。MpZOU2的重复调控使其成为功能研究的候选基因，同时强调仅凭转录水平变化并不能确证RCD的执行。
研究背景与问题：调控性细胞死亡（regulated cell death, RCD）是有机体发育、组织重塑、胁迫适应和免疫的重要组成部分，由遗传程序精确调控以清除多余或受损细胞。植物RCD在孢子体和配子体阶段均发挥关键作用，涉及大孢子发生、管胞分化、胚胎发生、叶片形态建成，以及盐碱、干旱、紫外线、重金属和极端温度等逆境下的受损细胞清除；在免疫中，过敏性反应（hypersensitive response, HR）相关的局部细胞死亡可限制病原增殖。然而，植物RCD的分子调控机制仍 largely 未知，现有知识主要来自被子植物（拟南芥、烟草、百日草等）。地钱（Marchantia polymorpha）作为苔藓植物，具有相对简单的调控架构和已充分研究的胁迫响应，为探究RCD相关基因在进化中的保守性与分化提供了合适体系。近期研究已在地钱中鉴定出遗传调控的发育性细胞死亡，如MpZOU1/MpICE1调控模块促进根椐细胞（pegged rhizoids）的终末分化与程序性死亡，MpMLO1依赖的钙信号调控精子释放过程中的程序性死亡；MpZOU2被鉴定为MpZOU1的古老旁系同源基因，但其细胞死亡功能未知。目前尚不清楚地钱中多种胁迫是否会激活共同的RCD转录响应，还是以胁迫特异性组合招募RCD相关基因。已有地钱转录组研究多聚焦于单一胁迫或特定生物学过程，未直接回答RCD相关转录响应是否在独立实验中重复出现，也未应用共表达网络分析整合非生物与病原数据。为此，研究人员对公共RNA-seq数据进行了跨胁迫整合再分析。
研究人员开展的研究：研究人员筛选公共数据库中的地钱RNA-seq数据集，依据野生型材料接受胁迫处理且存在匹配对照的标准，选取3项已发表研究共36个样本，构成9组胁迫与对照比较：8组非生物胁迫（来自Tan等人（Tan et al.）研究的冷、黑暗、热、高光、缺氮、渗透和盐胁迫，以及来自Schr?der等人（Schr?der et al.）研究的独立高光实验）和1组病原菌侵染（来自Grenz等人（Grenz et al.）研究的丁香假单胞菌（Pseudomonas syringae）侵染）。研究人员对原始测序数据统一进行质量控制和定量分析，使用Salmon进行转录本定量，以地钱v6.1注释汇总至基因水平；对每组比较独立使用DESeq2鉴定差异表达基因（differentially expressed genes, DEGs），阈值设为padj < 0.05且|log2倍数变化| ≥ 1。随后检查9组比较的共享DEGs，并评估调控方向一致性。研究人员整合全部样本，使用ComBat-seq进行批次校正，经方差稳定变换后开展主成分分析（PCA）和加权基因共表达网络分析（WGCNA），鉴定与特定胁迫相关的共表达模块。同时，研究人员从文献和MarpolBase注释中整理14个地钱RCD相关基因的精选列表（涵盖硫堇、NADPH氧化酶、metacaspase、自噬相关蛋白、ZOU/ICE转录因子等），分析其在9组比较中的差异表达模式。最后对显著胁迫相关模块进行基因本体（Gene Ontology, GO）生物学过程富集分析。
研究结果：(1) 整合转录组分析捕获多种胁迫响应——9组比较的DEGs数量从渗透胁迫的1252个到侵染的4789个不等。两独立高光实验共享1211个DEGs，表明高光转录响应具有高度可重复性；盐与渗透胁迫共享741个DEGs，低于高光间的共享数。(2) 9组比较共有的DEGs仅34个，其中仅12个在所有比较中方向一致（3个上调，包括线粒体导入基因Mp7g17400；9个下调，包括脂氧合酶Mp2g16000）；其余22个包括两个转录因子（AP2/ERF家族因子Mp5g22160和bHLH因子Mp2g00890）表现出方向相反的变化。8组非生物胁迫共有的DEGs为51个，少于9组共有的结果说明生物与非生物胁迫部分重叠但各有特异程序。(3) 14个RCD相关基因中，10个在至少一组比较中差异表达；MpZOU2是唯一在多于一种胁迫中受调控的基因，在冷、黑暗和Tan高光及侵染中均显著变化（冷、黑暗、侵染中上调，Tan高光中下调）；其余9个仅在单一比较中显著，其中侵染比较中就有5个基因上调（MpRBOHA、MpATG7、MpMCA-I、MpATG1等）。Schr?der高光比较中没有RCD基因达到显著阈值。14个基因中无一出现在9组共有DEGs中。(4) 批次校正有效减少研究来源相关变异，PCA显示校正后不同研究的对照样本彼此重叠，两独立高光样本聚在一起，盐和渗透胁迫聚类相近，而黑暗条件仍为明显离群。(5) WGCNA鉴定出64个共表达模块；其中turquoise模块（1931个基因）与侵染显著相关，富集蛋白质泛素化、蛋白酶体介导降解、内质网相关降解（ER-associated degradation, ERAD）、自噬和膜运输；green模块与缺氮和渗透胁迫显著相关，富集氧化应激反应、过氧化氢分解代谢和细胞壁组织。(6) GO富集分析显示turquoise模块涉及蛋白稳定性调控和细胞内囊泡运输等21个显著生物学过程，green模块则富集氧化应激和细胞壁重塑相关过程。
讨论与结论总结：讨论部分指出，地钱中RCD相关转录响应主要表现为胁迫特异性，而非统一程序。9组比较仅共享34个DEGs，其中方向一致者仅12个，支持植物胁迫响应由部分重叠但相对独立的网络组成这一观点。MpZOU2在四种独立条件下重复调控，使其成为功能研究的有力候选，但转录变化本身不能证明RCD执行。侵染比较中5个RCD基因上调与原始研究中观察到的组织坏死和毒素作用一致，但不能将病原诱导的坏死直接等同于典型RCD。Schr?der高光比较中RCD基因无显著变化，而Tan高光中MpZOU2显著调控，说明名义相同的胁迫类别不一定产生相同细胞结局。WGCNA鉴定的侵染相关模块揭示病原胁迫伴随蛋白质质量控制、自噬和膜运输的协同激活；缺氮/渗透相关模块则关联氧化还原稳态和细胞壁重塑。研究结论总结：地钱对环境及病原胁迫的响应依赖于高度背景特异的转录网络；少量RCD相关基因出现重复性响应，而大多数变化为胁迫特异。侵染引发最强的RCD基因协同诱导，并伴随富含蛋白质降解、ERAD、自噬和膜运输的大型共表达模块；高光、渗透和缺氮虽引起广泛转录变化，但RCD基因响应有限或缺失。这些发现强化了胁迫驯化与RCD执行之间的区别，提示保守的细胞死亡相关组分可能是选择性招募而非通过通用途径激活。未来需对MpZOU2等重复候选进行功能研究，并结合细胞水平的RCD执行测量，以确定这些转录特征是否代表古老植物细胞死亡网络的保守组分，还是仅与RCD部分交叉的胁迫响应基因。
