整合scRNA-seq和snRNA-seq与空间转录组学以解开木质部之谜

Integrating scRNA-seq and snRNA-seq with spatial transcriptomics to unlock the xylem puzzle

作者信息Mingke Wei, Jo-Wei Allison Hsieh, Jr-Fong Dang, Botong Tong, Hui Li, Ronald Sederoff, Ying-Chung Jimmy Lin, Quanzi Li
PMID41724975
发布时间2026-02-23
DOI10.1186/s13059-026-04007-z

实验完整度

高

包含snRNA-seq与scRNA-seq整合分析、解剖学验证、LCM RNA-seq相关性分析、RNA原位杂交及机器学习分类,多层级验证机制与轨迹。

主要模型

Populus alba × P. glandulosa cv. "84K"茎发育木质部(SDX) 原生质体scRNA-seq 单核RNA测序(snRNA-seq) 激光捕获显微切割(LCM)RNA-seq

重点核对

杨树品种(P. alba × P. glandulosa cv. "84K")和茎段位置(15节以下) 细胞捕获方法(原生质体分离、刮取、冻切) 整合分析参数(Seurat锚点K=5) 差异表达阈值(log2FC≥1,p<0.05) 用于相关性分析的LCM RNA-seq数据来源

摘要

背景:木质部发育是一个动态、连续的过程,对木本植物的次生生长和地球上的生物量积累至关重要。虽然单细胞RNA测序(scRNA-seq)能够重建早期木质部分化轨迹,但其依赖原生质体分离,排除了具有加厚次生细胞壁的晚期细胞,使得次生细胞壁沉积和程序性细胞死亡等关键阶段难以表征。结果:我们对杨树的茎发育木质部进行了单核RNA测序(snRNA-seq),并整合了之前的scRNA-seq数据集,以重建木质部形成的全面发育景观。解剖学验证证实,scRNA-seq图谱主要代表早期茎发育木质部,而snRNA-seq富集了更深的、次生细胞壁沉积层。整合分析揭示了一个在空间和转录上定义的次生细胞壁区域,该区域得到木质素自发荧光及其与激光捕获显微切割来源的转录组相关性的支持。差异表达和基因本体论分析揭示了木质素生物合成和程序性细胞死亡相关基因的富集,表明次生细胞壁形成和程序性细胞死亡启动在转录上是协调的。无监督聚类和支持向量机分类的机器学习进一步揭示了早期木质部细胞比晚期细胞具有更大的转录组异质性。结论:我们的研究结果表明scRNA-seq和snRNA-seq平台具有高度兼容性和互补性。这种整合方法不仅克服了单一技术的局限性,而且为解析植物复杂发育轨迹提供了一个实用且可扩展的框架。它成功重建了木质部发育的完整连续体,从形成层祖细胞到次生细胞壁形成再到程序性细胞死亡。

实验结论

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

研究问题
如何利用snRNA-seq和scRNA-seq的互补性,完整解析木质部发育轨迹,并明确次生细胞壁沉积和程序性细胞死亡的转录特征?
核心机制
次生细胞壁沉积和程序性细胞死亡的启动在转录上协调进行,相关基因(如木质素生物合成和PCD相关基因)在相同转录区域内共表达。
主要证据
整合scRNA-seq和snRNA-seq的UMAP分析显示两个平台互补,snRNA-seq捕获晚期SCW沉积细胞;LCM RNA-seq相关性、解剖学验证和RNA原位杂交支持SCW和PCD区域的存在;机器学习分析显示早期细胞异质性更高。
研究意义
该研究为植物单细胞研究提供了一个实用且可扩展的框架,通过整合多平台数据克服技术局限,并为全面理解木质部发育提供了新见解。

研究路径

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

1

解剖学分析验证采样互补性

验证去皮、刮取和冻切等处理对木质部细胞层捕获的影响,确认scRNA-seq和snRNA-seq采样区域。

对完整茎、去皮茎、原生质体分离、刮取及冻切样本进行甲苯胺蓝染色和自发荧光观察,比较不同处理下保留的细胞层。

2

单核RNA测序与数据整合

获得茎发育木质部单核转录组,并与已有scRNA-seq数据整合,重建连续发育轨迹。

从两个生物学重复的杨树茎发育木质部制备并纯化细胞核,构建snRNA-seq文库,测序后与scRNA-seq数据整合分析。

3

鉴定次生细胞壁沉积区域

在snRNA-seq数据中识别进行次生细胞壁沉积的细胞群体,并验证其空间位置。

通过LCM RNA-seq与snRNA-seq的Pearson相关性分析,并根据冷冻切片中保留的细胞层和木质素自发荧光,鉴定SCW区域。

4

SCW区域转录组特征分析

探究SCW沉积区域的基因表达特征,揭示与木质素生物合成和PCD相关的转录程序。

对SCW区域进行差异表达分析,并利用GO富集分析,识别富集的生物过程。

5

基因表达模式验证

验证通过转录组分析预测的候选基因在木质部中的空间表达模式。

选择SCW区域上调基因(如ACL5、FLA12、XCP2等)和早期簇基因(如RIN4等)进行RNA原位杂交。

6

转录异质性与分化轨迹分析

比较早期和晚期木质部细胞的转录组异质性,并推断分化轨迹。

对SDX-snRNAseq数据进行K-means聚类(K=2-20),通过采样归一化比较早晚期细胞聚类数,并使用SVM分类器验证分离;使用Slingshot推断轨迹。

研究方法

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

产品清单

实验环节名称品牌货号
CytoFLEX SRTBeckman--
Chromium单细胞3'试剂盒(V3.1)10x Genomics--
Illumina NovaSeq 6000Illumina--
Cell Ranger10x Genomicsv7.2.0
振动切片机LeicaVT1200S
冷冻切片机LeicaCM1950
OCT包埋剂Sakura Finetek--
光学显微镜OlympusBX53
旋转切片机LeicaRM2255
转录试剂盒Roche11175025910
NBT/BCIPThermo--

关键环节

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

环节核对要点
植物材料与生长条件
杨树品种(P. alba × P. glandulosa cv. '84K')、年龄(5个月或3个月)、茎段节位(15节以下)、生长条件
阅读提示:Methods(SnRNA-seq library construction and sequencing, Histological analyses)
细胞核分离与纯化
细胞核分离方法(参考文献26)、细胞核纯化仪器(CytoFLEX SRT)
阅读提示:Methods(SnRNA-seq library construction and sequencing)
文库构建与测序
文库试剂盒版本(Chromium Single Cell 3' Reagent Kits V3.1)、测序平台(Illumina NovaSeq 6000)、测序深度
阅读提示:Methods(SnRNA-seq library construction and sequencing)
生物信息学分析
比对参考基因组(P. trichocarpa)、Cell Ranger版本(v7.2.0)、Seurat参数(锚点K=5)、差异表达阈值(log2FC≥1, p<0.05)、GO数据库(TAIR)
阅读提示:Methods(SnRNA-seq analysis)
组织学分析
振动切片厚度(50μm)、OCT包埋、冷冻切片温度(-20℃)、TBO浓度(0.05%)
阅读提示:Methods(Histological analyses)
RNA原位杂交
固定液(FAA)、切片厚度(10μm)、探针合成方法(SP6/T7启动子)、检测试剂(地高辛、anti-DIG-AP、NBT/BCIP)
阅读提示:Methods(RNA in situ hybridization)
数据整合与聚类
整合方法(Seurat锚点)、K-means聚类范围(K=2-20)、采样归一化方法、SVM分类器参数
阅读提示:Results(Late-stage xylem cells exhibit less heterogeneity),Methods(SnRNA-seq analysis)