推荐公众号奥智生物，分享前沿组学技术、实用生信技能、数据挖掘思路
全球气候变暖加剧高温-强光复合胁迫，严重威胁果树产量与品质。近年来，单细胞 / 单细胞核转录组（sc/snRNA-seq）、空间转录组等前沿技术极大推动植物逆境研究。其中 snRNA-seq 无需制备原生质体，可直接利用田间液氮冷冻样本，规避酶解解离带来的人工转录扰动，尤其适配木质化程度高的多年生果树材料。
但目前果树果实基于真实田间逆境的 snRNA-seq 研究仍然稀缺，能够将单细胞组学挖掘候选基因并完成多体系功能验证、形成完整证据链的工作更为少见。本篇发表于 Horticulture Research 的研究，以葡萄浆果日灼为对象，构建「田间生理-单核转录组-时序调控网络-基因功能验证」完整研究范式，为园艺作物非生物胁迫单细胞研究提供高质量参考案例。
奥智技术团队为本研究提供了单核转录组测序与分析服务。
● 发表期刊：Horticulture Research（IF=9.5）
日灼是高温-强光复合胁迫引发的葡萄重要生理性病害，会造成果实褐变坏死，显著降低果实商品价值。现有葡萄日灼研究多采用 Bulk 转录组，掩盖组织内部细胞异质性；且多数试验依靠离体果实或人工装置模拟胁迫，缺乏大田自然胁迫下细胞层面的应答解析，不同细胞如何分工协同抵御日灼的分子机制尚不明确。
snRNA-seq 可直接解析田间冷冻胁迫样本，规避原生质体制备带来的人为胁迫干扰。本研究以此技术为手段，结合生理生化、单细胞组学、时序网络与基因得失功能验证，解析葡萄浆果对日灼的细胞特异性应答机制。
大田条件下，随日灼胁迫程度加重，浆果细胞膜损伤持续加剧，MDA（丙二醛）含量不断上升。抗氧化酶 SOD、POD、CAT 活性在轻度日灼 SB1 达到峰值，中重度日灼（SB2、SB3）阶段酶活下降，反映抗氧化系统发生耗竭。 非酶抗氧化物质呈现双相响应特征：总酚、类黄酮在 SB1 阶段下调，中重度胁迫下积累；GSH（谷胱甘肽）持续升高；叶绿素、类胡萝卜素等光合色素随胁迫发生降解。
值得关注：尚未肉眼可见褐变的 SB1H（日灼病灶周边光晕组织）已经出现显著生理扰动，属于日灼早期潜在受损过渡组织。基于梯度生理结果，筛选 NE、SB1H、SB1 三个关键状态用于后续 snRNA-seq 测序。
构建葡萄浆果单细胞核转录组图谱，注释得到三大功能细胞类群
研究对 23,286 个高质量细胞核完成无监督分群，获得 12 个亚群，归并注释为 3 大细胞类型：EC（外果皮细胞）、MC（中果皮细胞）、VC（维管细胞）。
本研究采用多层交叉验证完成细胞注释：结合跨物种同源标记基因、亚群表达相关性、浆果分组织 RT-qPCR 验证 marker 基因的组织特异性。功能特征上，外果皮细胞富集角质层合成通路；中果皮高表达光合相关基因；维管细胞富集信号转导与物质转运通路；受标记基因资源限制，本研究未进一步拆分韧皮部、木质部亚群，统一注释为维管细胞 VC。
石蜡切片形态证据表明：日灼后浆果角质层、维管束面积显著缩减；胁迫样本中维管细胞占比下降，证实维管细胞对日灼胁迫高度敏感。
图2 葡萄浆果 snRNA-seq 细胞分群、细胞占比、组织石蜡切片及差异基因统计。
三类细胞对日灼胁迫既存在全局协同应答，又具备显著细胞类型特异性转录重编程
差异基因（DEGs，筛选标准：\(|log _{2} FC| ≥0.36\)，\(p<0.05\)）统计结果显示，无论 SB1H vs NE，还是 SB1 vs NE，维管细胞的差异基因数量显著高于外果皮、中果皮，发生最剧烈转录重编程。
GO 富集分析表明：全部细胞共享高温、高光、ROS 应答通路，同时抑制核糖体生物合成；而细胞应答存在特异性：维管细胞富集光受体、信号转导、物质转运；中果皮侧重碳水代谢与光保护；外果皮特异激活脂质、芪类次生代谢与损伤防御通路。
激素通路展现细胞分工：SA（水杨酸）、JA（茉莉酸）、ETH（乙烯）合成基因富集于外果皮与维管细胞；激素信号转导基因主要集中在维管细胞；激素响应基因在三类细胞广泛表达，不同细胞分别承担激素合成、信号转导、胁迫响应功能。
图3 植物激素相关差异基因的 UMAP 分布与表达气泡图。
伪时间轨迹与 hdWGCNA 网络共同证实 “维管细胞→中果皮细胞→外果皮细胞” 程序化级联应答
利用 Monocle 2 开展伪时间轨迹分析，细胞沿胁迫强弱顺序有序排布：伪时间起点以维管细胞为主，负责胁迫信号感知、物质转运；中段富集中果皮细胞，执行光保护、代谢底物供给；末端以外果皮细胞为主，启动角质、芪类物质合成，实现物理-化学防御。
hdWGCNA（高维加权基因共表达网络分析）构建得到 6 个共表达模块。模块激活时序与伪时间结果高度吻合：vvi3（维管主导）早期激活、vvi2（中果皮主导）中期激活、vvi4（外果皮主导）晚期激活，从网络层面证明葡萄浆果存在有序的细胞级联胁迫应答程序。
图4 高维加权基因共表达网络 hdWGCNA 分析结果。
外果皮特异性 VvVINST1 多维度正向调控葡萄浆果抗日灼能力
snRNA-seq 与 RT-qPCR 证实，芪合酶基因 VvVINST1 特异性在外果皮表达；基因转录水平与其催化产物反式白皮杉醇（trans-piceid）的含量均随日灼胁迫加重而上调。
多套得失功能实验表明：42 ℃、70,000 lx 强光胁迫条件下，异源过表达 VvVINST1 的 Micro-Tom 番茄抗日灼能力提升；葡萄浆果瞬时过表达可减轻日灼损伤，VIGS 沉默后日灼褐变皱缩表型加剧，证明 VvVINST1 是抗日灼正向调控基因。
机制层面，VvVINST1 促进外果皮积累 trans-piceid；提升 SOD/POD/CAT 酶活，缓解氧化损伤；维持 JA-ABA-SA 激素稳态；上调苯丙烷-芪类通路大部分结构基因（VvMYB14 反馈下调），多维度协同提升浆果日灼抗性。
图5 VvVINST1 调控芪类代谢通路的代谢物与基因表达检测。
建立葡萄浆果日灼 “细胞类型特异性程序化响应” 综合模型
整合生理表型、单核转录组、基因功能验证多维度证据，本研究提出葡萄浆果日灼细胞类型特异性程序化响应模型。高温-强光胁迫下，维管细胞最先感知并转导胁迫信号；中果皮细胞提供光保护与代谢物质；外果皮启动角质、芪类次生代谢，构建物理 + 化学防御屏障；不同细胞时序协同，共同抵御日灼伤害。
该模型解析急性复合胁迫下浆果器官空间-时序协同防御逻辑，为果树抗逆分子育种提供理论框架和候选基因资源。
图6 葡萄浆果日灼的细胞类型特异性程序化响应工作模型。
该研究充分发挥 snRNA-seq（单细胞核 RNA 测序）的技术优势：可直接使用大田冷冻样本，规避原生质体酶解带来的人为转录干扰，为多年生木本园艺作物开展田间逆境的单细胞解析提供可借鉴的研究范式。
面向未来，果树逆境组学研究正在从单纯的单细胞转录组，走向 snRNA-seq + 空间转录组 + 多组学（ATAC、代谢组）联合解析的新阶段。snRNA-seq 可以解析每个细胞核的转录状态，空间转录组能够还原细胞在组织内部原位空间位置，二者结合，将进一步解析器官内部细胞间通讯、信号分子长距离传递，真正还原果树果实应对田间复合胁迫的原位分子图景；而多组学联合，能够打通 “染色质开放-基因转录-代谢物积累” 完整调控链路，极大加速抗逆基因挖掘与分子育种进程。
然而木本果树材料本身存在细胞壁厚、木质化程度高、田间胁迫样本珍贵、取材难度大等现实痛点，对样本前处理、细胞核提取、建库质控、生物信息深度分析均提出极高门槛。奥智团队在植物单细胞、单细胞核转录组领域具备行业绝对实力：针对果树、林木等难处理木本材料建立专属细胞核提取优化体系，能够兼容液氮冷冻田间逆境样本；同时配套完整的高级生物信息分析管线，已经协助众多课题组产出园艺、植物领域高水平成果。
![]()
*未经许可，不得以任何方式复制或抄袭本篇文章之部分或全部内容。版权所有，侵权必究。
#单细胞转录组#单细胞核转录组#葡萄日灼研究#文献解读
