植物生长、抗虫、抗旱、开花结实，几乎都由生长素、细胞分裂素、赤霉素、脱落酸、乙烯、油菜素内酯、独脚金内酯、茉莉酸、水杨酸这九类经典激素协同调控。但传统研究有两个老问题：一是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相关组分来精细调控气孔行为。这些发现共同揭示了植物激素共享的调控支架与定制化信号回路。
