单细胞与空间RNA测序揭示果实衰老的时空轨迹

Single-cell and spatial RNA sequencing reveal the spatiotemporal trajectories of fruit senescence.

作者信息Xin Li, Bairu Li, Shaobin Gu, Xinyue Pang, Patrick Mason, Jiangfeng Yuan, Jingyu Jia, Jiaju Sun, Chunyan Zhao, Robert Henry
PMID38600080
发布时间2024-04-10
DOI10.1038/s41467-024-47329-x

实验完整度

高

包含单细胞RNA测序、空间转录组、RNA-FISH、VIGS功能验证及代谢组等多层证据

主要模型

火龙果(Hylocereus undatus)成熟与衰老期果皮组织 10× Genomics单细胞转录组 10× Visium空间转录组 VIGS沉默HuCMB1的果实

重点核对

成熟样品CK与贮藏后衰老样品Post的分组及衰老判定标准(第9天超氧阴离子爆发、类黄酮第二峰、表型皱缩或病斑) 单细胞悬液制备的酶解体系(1%纤维素酶、1%半纤维素酶、0.3%果胶酶等)及细胞活率(CK 80%、Post 77%) 空间转录组组织透化时间3 min及测序深度约250M read-pairs/样品 VIGS沉默HuCMB1的片段长度178 bp、农杆菌GV3101悬浮液OD600=0.6、5个生物学重复 RT-qPCR以β-actin为内参、2−∆∆Ct法计算、3个生物学重复

摘要

果实衰老是一个复杂的生理过程,果皮内存在多种细胞类型,使得阐明它们在果实衰老中的各自作用极具挑战性。本研究构建了火龙果(Hylocereus undatus)果皮的单细胞表达图谱,揭示外果皮和中果皮细胞在果实衰老过程中发生最显著的变化。拟时序分析建立了衰老过程中的细胞分化和基因表达轨迹。早期氧化应激失衡之后,外果皮细胞的抗性被激活,随后衰老相关蛋白在衰老晚期的中果皮细胞中积累。研究揭示了早期响应因子HuCMB1在衰老调控网络中的核心作用。本研究为深入理解植物动态衰老过程提供了时空视角。

实验结论

提炼研究问题、关键发现与证据,快速把握文章的核心贡献。

研究问题
果实衰老过程中果皮不同细胞类型如何参与以及其时空动态轨迹如何。
核心机制
衰老进程呈现ROS-抗性-衰老的时序:中果皮早期ROS积累,随后外果皮抗性基因(苯丙烷/类黄酮途径)激活,晚期中果皮衰老相关蛋白如HuSAG12积累;早期响应因子HuCMB1在调控网络中居核心,可能通过抑制类黄酮合成加速衰老。
主要证据
单细胞与空间转录组整合鉴定5种果皮细胞类型;拟时序与RNA velocity揭示分化轨迹;RNA-FISH验证HuSAG12在中果皮定位;VIGS沉默HuCMB1导致果实加速衰老且类黄酮早期合成显著降低;代谢组支持外果皮富集类黄酮。
研究意义
提供果实衰老的时空视角,揭示早期响应基因可作为监测果实衰老的候选报告基因,为减少采后损失和延长贮藏期提供思路与技术。

研究路径

按研究推进顺序梳理实验设计、验证步骤与关键观察。

1

构建果皮单细胞表达图谱

解析成熟与衰老期果皮细胞类型及表达异质性。

取成熟(CK)与贮藏后衰老(Post)果皮,酶解制备原生质体,经10× Genomics Chromium平台建库测序,Cell Ranger处理,Seurat聚类注释。

2

空间转录组定位细胞类型

将单细胞聚类对应到果皮空间区域。

对CK果皮进行10× Visium空间转录组测序,Seurat SNN聚类,鉴定区域标记基因,并用SingleR、SciBet、RCTD、CARD四种算法整合单细胞与空间数据。

3

拟时序与RNA velocity分析衰老轨迹

推断外果皮与中果皮细胞在衰老中的分化时序与状态转换。

对EX和ME细胞进行Monocle拟时序分析,构建四类功能基因集(衰老52、抗性642、ROS 1913、苯丙烷258),基因集评分筛选9019个高相关细胞,再进行拟时序与RNA velocity分析。

4

ROS与类黄酮表型验证

验证拟时序推断的衰老时序。

检测贮藏过程中果皮ROS(超氧阴离子、H2O2)和类黄酮水平,使用荧光探针、组织化学染色和分光光度法。

5

HuCMB1功能验证

验证HuCMB1在果实衰老中的调控作用。

构建pTRV2-HuCMB1载体,农杆菌GV3101侵染果实进行VIGS沉默,观察表型、失重率、类黄酮含量及RT-qPCR检测基因表达。

研究方法

按研究目的归类文中使用的方法,便于定位所需技术。

产品清单

实验环节名称品牌货号
纤维素酶----
半纤维素酶----
果胶酶----
甘露醇----
牛血清白蛋白----
台盼蓝----
Chromium单细胞3ʹ v3.1试剂盒10× Genomics--
Bioanalyzer 2100生物分析仪Agilent--
Qubit高灵敏度DNA定量仪Thermo Fisher Scientific--
NovaSeq6000测序仪Illumina--
Visium空间基因表达载玻片10× Genomics--
Visium组织优化载玻片10× Genomics--
NovaSeq S4试剂盒(200循环)Illumina20027466
DAPI----
蛋白酶K----
三氯化铈Sigma--
Spurr包埋剂----
H2DCFDA荧光探针----
盐酸羟胺----
芦丁标准品--153-18-4
混合球磨仪MM 400Retsch--
Shimadzu Nexera超高效液相色谱系统Shimadzu--
Phenomenex Kinetix C18色谱柱Phenomenex--
AB Sciex X500B质谱仪Sciex--
RNAprep Pure微量RNA提取试剂盒TIANGENDP432
TransScript一步法gDNA去除与cDNA合成SuperMix试剂盒TransGenAH311-02
pTRV2载体----
根癌农杆菌GV3101菌株----

关键环节

汇总复现实验时建议重点确认的条件及原文阅读提示。

环节核对要点
样品准备与衰老判定
成熟样品为授粉后35天果实;贮藏条件25°C、85%相对湿度;衰老样品判定为第9天出现超氧阴离子爆发、类黄酮第二峰及表型皱缩或病斑
阅读提示:Methods中“Preparation of mature and senescent samples”小节
单细胞悬液制备
酶解液组成(1%纤维素酶、1%半纤维素酶、0.3%果胶酶、0.6 M甘露醇、100 mM MES、20 mM KCl、20 mM CaCl2、10% BSA);酶解时间0.5 h;过滤孔径100 μm和40 μm;细胞浓度400–1200 cells/μl
阅读提示:Methods中“Sample collection and preparation of single-cell suspensions”小节
空间转录组
组织透化时间3 min;测序深度约250M read-pairs/样品;NovaSeq S4试剂盒200循环;Read 1 28 cycles、Read 2 91 cycles
阅读提示:Methods中“Spatial transcriptomic (stRNA-seq) experiments”小节
拟时序与网络分析
基因集评分筛选高相关细胞的四分位数法;SCODE运行50次取平均;MCODE参数Node Score Cutoff 0.2、K-Core 2、Max Depth 100
阅读提示:Methods中“Pseudotime analysis”和“Visualization of pseudotime-related gene regulatory networks”小节
VIGS功能验证
HuCMB1片段178 bp;农杆菌GV3101悬浮液OD600=0.6;25°C培养;5个生物学重复;贮藏9天后检测
阅读提示:Methods中“Virus-induced gene silencing (VIGS) of HuCMB1”小节
RT-qPCR
内参基因β-actin;2−∆∆Ct法;3个生物学重复;引物序列见Supplementary Data 10
阅读提示:Methods中“Gene expression analysis by RT-qPCR”小节及Supplementary Data 10
类黄酮与ROS检测
H2DCFDA浓度10 mM、25°C避光0.5 h;超氧阴离子羟胺氧化法530 nm;类黄酮氯化铝比色法510 nm、芦丁标准品CAS 153-18-4
阅读提示:Methods中“Flavonoid and ROS assessment”小节