摘要

从单细胞微生物到复杂的多细胞生物，所有生物在生命过程中都会不可避免地经历胁迫。从胁迫中恢复的能力是决定生物整体抗逆性的基本特征，但目前对胁迫恢复的研究仍不够深入。为探究植物如何从干旱中恢复，研究者对中度干旱后复水15 min开始的过程进行了精细时间序列的RNA测序分析。结果发现，干旱恢复是一个快速过程，涉及数千个恢复特异性基因的激活。为捕捉拟南芥（Arabidopsis thaliana，A. thaliana）叶片不同细胞类型中的这些快速恢复响应，研究团队在干旱恢复起始阶段进行了单细胞核转录组分析，鉴定出一种在不同细胞类型中独立形成的细胞类型特异性转录状态。为进一步验证干旱恢复过程中观察到的细胞类型特异性转录变化，研究者采用基于多重抗错荧光原位杂交（MERFISH）的空间转录组技术，明确了拟南芥叶片组织中恢复诱导型基因表达程序的解剖学定位。此外，研究还发现，恢复过程会自主诱导免疫系统激活，且这种激活能在体内增强拟南芥、野生番茄（Solanum pennellii）和栽培番茄（Solanum lycopersicum cv. M82）的病原菌抗性。由于复水会促进微生物增殖，进而增加感染风险，因此干旱恢复诱导的免疫激活可能对植物在自然环境中的存活至关重要。这些发现表明，干旱恢复与预防性防御响应同时发生，揭示了在植物不同细胞类型中调控胁迫恢复的复杂机制。

主要结果

图1. 干早恢复的精细时间分辨率RNA-seq揭示恢复特异性基因

首先，研究者使用精细时间分辨率RNA-seq技术，以蛭石培养的拟南芥莲座叶为研究对象，在中度干旱（土壤相对含水量30%）及复水后6 h内的7个时间点采集样本，并以同期正常浇水植株作为对照以排除昼夜转录变化干扰；通过差异基因表达分析，将基因分为干旱响应基因和恢复特异性基因两类，结合K-means聚类分析基因表达模式。结果显示，共鉴定出1248个干旱响应基因，其中53%在恢复阶段仍有差异表达，且97%在复水6 h后恢复正常表达水平；同时鉴定出超3000个恢复特异性基因，这些基因在干旱期无差异表达但在复水后不同时间点表现出差异，且呈现早期恢复诱导、晚期恢复诱导、复水下调三种表达模式，证明干旱恢复不仅是干旱诱导变化的逆转，还涉及大量恢复特异性基因的激活。

图2. 复水15 min后叶片相对含水量和气孔导度部分恢复

针对拟南芥叶片水分生理指标，分别使用叶片相对含水量（RWC）测定法和叶片孔隙计测定气孔导度（gsw）技术，对比正常浇水（WW）、中度干旱（30%的土壤相对含水量，D）、复水15 min（R15）、复水1 h（R60）四种处理下的叶片状态。结果显示，中度干旱时叶片RWC显著降低，复水15 min后无明显变化，复水1 h后则显著回升；而气孔导度在复水15 min后就已显著向正常水平恢复，复水1 h后完全恢复，且气孔导度的恢复速度独立于叶片RWC的缓慢变化，表明在中度干旱胁迫下，叶片水分含量和气孔导度均会下降，且气孔导度对复水的响应更迅速，1 h内即可完全恢复。

图3. 拟南芥叶片干旱复水后的细胞类型特异性转录响应

结合单细胞核RNA测序（snRNA-seq） 与MERFISH空间转录组技术，设置正常浇水、中度干旱、复水15 min、正常浇水+额外灌溉15 min（W15，排除根系机械刺激干扰）四种处理，对拟南芥叶片细胞核进行测序分析并验证细胞类型注释的空间一致性。结果显示，snRNA-seq数据经质控和无监督聚类得到27个代表不同细胞身份的簇，通过已知细胞类型标记基因注释及与已发表拟南芥叶片单细胞数据集投影比对，明确了22个簇的细胞类型（包括表皮细胞、叶肉细胞、维管细胞等），剩余5个未注释簇经GO富集分析显示独特功能特征（如免疫活性、代谢应激等）；MERFISH进一步验证了簇标记基因的解剖学定位与细胞类型注释一致，且snRNA-seq检测到的差异表达基因数量是RNA-seq的10倍，同时发现复水快速诱导多种细胞类型中乙烯响应因子（ERF）、WRKY等转录因子的表达，部分转录因子在分裂细胞中优先上调。

图4. 独特的叶肉细胞亚群在干旱恢复期间表现出特定基因表达程序

随后，基于snRNA-seq数据的亚聚类分析、基因共表达网络分析（WGCNA）及MERFISH结果，对拟南芥叶片16种主要细胞类型中的细胞进行亚聚类，重点分析叶肉细胞亚群（Mes 2）在不同水分处理组（W、D、W15、R15）中的基因表达特征。结果显示，在表皮细胞、叶肉细胞、排水器细胞、筛管分子细胞中均鉴定出高度富集（细胞比例超过50%）于干旱恢复植株的亚群，其中叶肉细胞亚群 Mes 2的亚簇4在复水组（R15）中高度富集，代表进入恢复特异性细胞状态（RcS）；对这些恢复富集亚群进行基因共表达网络分析，鉴定出富集基因模块及核心枢纽基因，其中212个枢纽基因中有50个在至少两个模块中共享，且70%的免疫活性细胞簇标记基因与这些共享枢纽基因重叠；MERFISH空间转录组进一步验证了叶肉细胞中SRC2、AT4G29780、SZF1等枢纽基因在复水早期的特异性富集，明确了RcS的组织定位和时间动态。

图5. 拟南芥从中度干旱恢复后增强了对病原体的抗性

将RNA-seq数据与生物诱导子转录响应数据关联分析，通过无菌低水分琼脂平板病原菌接种实验、土壤栽培植株病原菌喷雾接种实验，探究拟南芥干旱恢复后的免疫功能。结果显示，复水15 min后上调的恢复特异性基因中，50% 与已知生物诱导子（如Flg22、EF-Tu）触发的模式触发免疫（PTI）核心响应基因重叠，且免疫相关基因表达在复水早期达到峰值后随时间逐渐下降；无菌平板实验中，中度干旱（R-MD）和重度干旱（R-SD）恢复的植株接种丁香假单胞菌番茄致病变种（Pst DC3000）后，细菌载量显著低于正常浇水对照，且中度干旱恢复植株的抗病效果更优；土壤栽培实验中，复水90 min的植株接种Pst DC3000后，24 h后的细菌载量显著低于对照，而4 h初始接种时细菌载量无差异，证明干旱恢复诱导的免疫（DRII）能抑制病原菌增殖，且该响应不依赖于ABA信号。

图6. 干旱恢复诱导的免疫增强了栽培番茄S. lycopersicum cv. M82对穿孔黄单胞菌和Pst DC3000 的抗性

最后，为了探究DRII是否也能增加栽培番茄的病原菌抗性，研究团队采用病原菌注射器浸润接种法、病害严重度图像分析及细菌载量计数法，以Solanum lycopersicum cv. M82为研究对象，接种Pst DC3000和番茄斑点病病原菌穿孔黄单胞菌（Xanthomonas perforans 97-2），对比正常浇水（W）和中度干旱恢复（R）植株的抗病表现。结果显示，复水90 min后接种病原菌的番茄叶片，5 天后的病斑面积百分比显著低于正常浇水对照；初始接种4 h后各组叶片细菌载量无差异，但5天后恢复组叶片的细菌载量显著低于对照；两种病原菌处理下结果一致，证明干旱恢复诱导的免疫（DRII）在栽培番茄中依然保留，且对不同病原菌均有抗性增强效果。