植物生长、抗虫、抗旱、开花结实几乎都由生长素、细胞分裂素、赤霉素、脱落酸、乙烯、油菜素内酯、独脚金内酯、茉莉酸、水杨酸这九类经典激素协同调控。但传统研究有两个老问题一是bulk转录组把根、叶、维管、表皮混在一起细胞类型特异响应被“平均掉”二是植物单细胞多依赖原生质体酶解本身会诱发胁迫转录极早期激素信号容易失真。此外多数研究一次只做1-3种激素、单组织、单时间点九种激素在同一体系下怎么“跨细胞对话”一直缺全景图。近期南方科技大学翟继先、郭红卫团队联合东北师范大学刘智剑、云南大学姜凯等团队在Nature Communications发表题为“Cell-type specific early perception of nine phytohormones revealed by single-nucleus transcriptomics in Arabidopsis”的研究论文。该研究利用单核RNA-seq构建了拟南芥幼苗对九种植物激素早期响应的细胞类型图谱揭示了生长激素与防御激素截然不同的时空响应模式以及JA、SA、ABA三者之间非对称的交叉调控网络并发现了一个保卫细胞特异的SA/MYB60模块该模块通过借用ABA信号组分来精细调控气孔的开闭平衡。研究结果1.利用单核RNA测序构建全面的转录组图谱揭示九种植物激素的早期感知机制为构建拟南芥植物激素早期响应的综合时空转录组图谱研究收集了经九种植物激素外源处理的7日龄野生型幼苗处理时长分别为0.5h和3h同时设置相应的模拟处理mock作为对照用于制备snRNA-seq文库。质控降维注释获得35细胞簇其中31簇对应根/地上部已知类型根表皮、皮层、内皮层、维管、根冠、分生地上表皮、叶肉、保卫细胞、茎尖等。图1拟南芥幼苗早期植物激素响应单细胞核RNA-seq图谱2.不同激素处理观察到截然不同的早期反应趋势比较每种激素处理与对照确定每个时间点由激素诱导的差异表达基因DEGs。生长相关激素AUX、CK、GA、SL、BR、ET在0.5h处理时间点诱导的DEGs数量多于3h呈现早期偏向型模式防御/胁迫相关激素JA、SA、ABA在3h处理时间点诱导的DEGs数量更多、范围更广呈现持续或后期增加型模式。GO富集分析表明BR处理导致广泛的激素生物合成基因下调符合已知的负反馈调节机制。JA、SA、ABA处理后其自身生物合成通路相关基因在几乎所有细胞类型中均显著上调0.5h和3h表明这些激素在早期持续激活自身的合成程序。qRT-PCR验证通过0.5–24h的连续时间点检测发现激素响应并非简单线性。例如AUX/GH3.3、GA/GA20ox等在3h后逐渐减弱根部尤为明显而CK/ARR5、SL/BRC1在12h才达到峰值。即使是JA、SA、ABA在不同组织中也表现出非均匀的时间模式如JA在根中持续上升在地上部3h达峰后下降。将单细胞数据聚合成伪批量样本后超过30%的细胞类型特异性DEGs被掩盖尤其是ET、BR、SL、GA在3h时丢失超过50%的DEGs。功能富集分析显示伪批量分析会遗漏叶肉细胞中与叶绿体和光合作用相关的生物学程序凸显单细胞层面解析激素响应的必要性。生长激素如生长素在0.5h时DEGs主要在茎来源的细胞簇表皮、叶肉中共享根簇共享较少。防御激素ABA、SA、JA在0.5h时共享程度较低但在3h时扩展至多个地上部和根部簇提示其响应更具全局性。图2九种植物激素诱导的转录响应所呈现的时空异质性与组织内相干性3.激素应答相关的转录组重叠程度具有细胞类型与时间依赖性确定单个激素的动态变化后研究探讨了不同激素处理之间的转录响应重叠在不同细胞类型和时间点上如何变化。检测了所有激素两两比较中激素诱导的差异表达基因并区分处理0.5h和3h的上调和下调基因集涵盖了图谱中识别出来的聚类。用S-MOD超几何评分量化任意两激素DEG同向/反向重叠35簇聚出若干overlap module•同向重叠在0.5h整体最强处理3h整体信号减弱ABA-SA、ABA-JA等防御激素共同抑制。• 反向重叠多跨时间例如0.5h ET与3h ABA/SA/JA常呈拮抗提示早期的ET爆发可能会抑制、延迟或重新定向后续的防御程序。图3激素差异表达基因重叠的时空异质性4.JA、SA及ABA相关信号通路之间呈现出强烈但不对称的交叉调控作用JA、SA、ABA三者的响应基因重叠度显著高于其他激素对。例如3h约50%的ABA响应基因同时也对SA响应。然而这种高重叠并不意味着对称激活——通路层面的共调控是不对称的。SA和ABA处理均提高了JA生物合成和分解相关基因的表达特征但二者之间仍存在一些差异ABA处理会上调SA合成相关基因但对SA降解/失活基因也有上调作用。SA处理则导致ABA降解/失活基因在3h大幅下调这可能是ABA与SA之间高DEG重叠的部分原因。利用单细胞分辨率研究发现激素通路的交叉调控存在时空异质性。例如JA处理后ABA降解/失活基因在维管系统中下调但在茎和根的其他细胞中上调而ABA合成基因则相反。这表明激素交叉调控的模式是细胞类型依赖的。与公开数据集对比后发现同源激素通路的直接激活在不同数据集中基本保守但交叉响应高度依赖于处理条件尤其是激素浓度。该研究使用的SA浓度500μM高于公共Okamoto数据集300μM这解释了为何研究中SA-JA的交叉响应更强。通过单激素、两两组合和三激素组合的bulk RNA-seq实验发现大量基因表现出非加性转录响应。两两互作的强度不对称ABA×SA和SA×JA影响的基因集大于ABA×JA。GO富集分析显示两两互作基因不仅富集对应两条激素的通路还富集到第三条激素的通路。三重互作ABA×SA×JA特有基因在0.5h以上调为主富集防御/胁迫过程到3h转为下调为主呈现明显的时间转换特征。尽管全组织水平上尤其是地上部分存在明显的共向差异表达基因重叠模式但在保卫细胞中无论处理时长如何这三种激素诱导的DEG重叠度始终是茎部所有细胞类型中最低的。图4通路层级分析揭示了JA-SA-ABA信号通路的不对称响应重叠及交叉调控作用5.单细胞分析精准定位SA-ABA保卫细胞相互作用中一个与MYB60相关的调控节点保卫细胞在JA/SA/ABA中重叠最低且激素响应主要局限于这三种激素而对其他激素仅检测到轻微的转录变化。这表明保卫细胞在JA-SA-ABA激素交互调控网络可能与其他茎部细胞类型中的调控网络存在显著差异。在所有茎部细胞类型中保卫细胞对ABA、SA、JA三种激素的共方向DEG重叠度最低。保卫细胞的激素响应几乎仅限于这三种激素对其他六种生长类激素仅有微小转录变化。这说明保卫细胞中JA-SA-ABA的交互调控网络与其他茎部细胞类型显著不同。通过整合TF-motif富集与共表达网络发现一个MYB60为中心的网络仅在SA处理下的保卫细胞中被检测到其他激素或其他细胞类型中未见。MYB60是已知的保卫细胞特异性转录因子促进气孔开放而SA通常被认为诱导气孔关闭这一矛盾指向一种延迟负反馈机制SA初期诱导气孔关闭后通过上调MYB60及其连接的基因如IQD2、WRKY70、MPK12等借用ABA信号通路组分实现气孔的及时重新开放防止过度关闭。这个网络分析将MYB60界定为一种选择性SA枢纽分子——它可利用ABA相关信号组分以细胞类型特异性方式精细调控气孔动态。野生型WT在SA处理后3h气孔孔径减小24h恢复后显著重新开放。myb60突变体在SA处理后气孔重新开放显著延迟直到48小时才恢复而DMSO对照正常恢复。对突变体进行snRNA-seq发现SA诱导的MYB60调节子中7个ABA相关气孔调控基因的诱导倍数在myb60保卫细胞中减弱证实MYB60参与了SA响应性气孔调控程序的激活。图5一个由SA触发且特异性作用于保卫细胞的转录因子模块揭示MYB60是连接SA与ABA信号通路的枢纽节点总 结总之该研究构建的约50万个拟南芥细胞核的单核图谱揭示了九种植物激素如何以细胞类型分辨率快速且差异化地重塑转录调控网络。该图谱捕捉到以下关键现象由ABA、JA和SA驱动的、覆盖整个组织的主导性重编程激素信号交互中从早期协同作用到后期拮抗作用的转变以及保卫细胞特异性SA/MYB60模块——该模块通过重新利用ABA相关组分来精细调控气孔行为。这些发现共同揭示了植物激素共享的调控支架与定制化信号回路。
