找回“丢失的空间坐标”:SDAS 助力突破单细胞分析局限,解锁时空组学新视野

2026.08.07 内容来源:华大时空

一、SDAS 是什么

SDAS(Spatial Data Analysis Software  空间转录组数据分析软件)是华大时空为 Stereo-seq 时空组学数据开发的一款 "开箱即用、解压即跑" 的Linux 命令行分析工具包。SDAS覆盖原始数据预处理至图表输出的全流程分析模块,研究者不必手动拼接各类脚本工具,即可完成标准化、高质量空间组学分析。配套还有算法测评专栏(Analysis Guide)、时空高级分析 Pipeline。

  • SDAS 工具包:14 个分析模块、31 个经实测筛选的算法。模块以命令行方式调用,关键参数透出,可单独使用也可组合。

  • 算法测评专栏(Analysis Guide):针对空间组学常见分析任务(细胞注释、空间共表达等),从算法原理、参数影响、计算资源、运行结果四个维度做基准评测,给出 Stereo-seq 数据上的实测推荐。

  • 时空高级分析 Pipeline:把 SDAS 工具包中的模块按典型研究思路串成可执行的 Python 流程,配套示例数据与 conf 配置文件,不依赖 Snakemake / Nextflow。

简言之,算法测评回答"为什么用这个算法",工具包回答"怎么用",Pipeline 回答"按什么流程用"。


二、SDAS 解决什么问题

在和课题组的交流中,我们观察到空间转录组分析常面临三类实际问题。

分析思路局限。 最常见的情况:拿到空间数据后,仍按单细胞的惯性推进——降维、聚类、差异基因、通路富集。这些步骤在 scRNA-seq 阶段已做过一遍,缺少空间坐标的分析,最终图表里找不到"细胞在组织中的位置",与单细胞结果无本质差异,实验溢价难以体现。

算法选择成本高。 仅"细胞类型注释"一项,社区就有 cell2location、RCTD、Tangram、 等多种主流方法,精度、运行时间、参考集依赖、显存占用差异显著。多数研究者选算法凭文献频次或组内经验,缺乏在 Stereo-seq 数据上的实测依据,往往一个项目周期内只能试 1-2 种方法。

分析环境配置耗时。 cell2location 的 PyTorch 依赖与 Scanpy 冲突、CellChat 的 R 包与 SeuratWrappers 版本不兼容、不同算法对 CUDA 版本要求不一致——一个完整的空间分析栈手动配置通常需要 1-2 周,期间多在排查环境问题而非做分析。


三、SDAS 工具包

3.1 特点概览

  • 全流程覆盖 : 一个工具整合空间转录组领域 31 款主流算法,封装成  14 大分析模块,覆盖原始数据预处理至图表输出

  • 安装简便:Linux 解压即用,集成并预编译所需依赖,无需单独配置 Python、R、CUDA 等复杂依赖环境

  • 生态兼容 :读取STOmics分析软件SAW直接输出的.h5ad/.h5mu数据, 或者按照使用手册的说明准备AnnData/Seurat 格式的数据

3.2 14 大分析模块介绍

普通单细胞转录组告诉你"有哪些细胞、表达了什么基因",但 丢掉了"细胞在哪里"这个关键维度。空间转录组在保留单细胞所有信息的基础上,额外提供"空间坐标",从而解锁了一系列单细胞做不了的分析。下表是 14 个模块的完整功能说明,以及分析模块所选择的算法是否支持单细胞转录组数据,还是仅支持空间转录组数据。

分析模块

功能/生物学意义

适用范围

数据预处理

读取 SAW 输出,h5ad文件进行多片合并、格式转换、信息打印及子集提取等操作

单细胞/空间共用

细胞类型注释(非聚类)

*集成空间组学专用方法

破解原位 Spot 的混合信号,精细还原每个位点由哪些细胞类型组成及其比例。

空间特有

空间结构域识别

结合基因表达与空间位置,自动划定具有特定功能或病理特征的组织微结构区域

空间特有

空间基因共表达

寻找在空间上呈现协同高表达的基因模块,揭示组织局域的功能联动。

空间特有

细胞邻域分析

洞察细胞周围的“微环境生态位(Niche)”,揭示病理局部微环境的细胞构成模式。

空间特有

细胞空间距离关系

量化不同细胞类型之间的真实物理距离,评估细胞间的空间共定位或物理互斥倾向。

空间特有

细胞通讯(加空间)

在真实空间距离约束下,精确预测空间邻近细胞间发生的配体-受体信号互作。

空间特有

CNV 分析

从基因表达推断染色体拷贝数变异,在原位切片上高精度辨识肿瘤细胞与正常细胞。

单细胞/空间共用

差异基因分析

挖掘不同组织区域、病理分级或细胞群之间特异表达的关键分子标记。

单细胞/空间共用

基因集富集分析

将差异分子定位到已知生物学通路,快速解释该区域活跃的细胞功能与代谢状态。

单细胞/空间共用

轨迹分析

推算细胞在组织物理空间中的表达演变趋势,模拟细胞分化或病变演进的动态过程。

单细胞/空间共用

转录因子分析

识别调控特定空间区域或细胞状态的“主控”转录因子及其下游基因网络。

单细胞/空间共用

蛋白互作网络

将关键基因映射到蛋白互作网络,锁定驱动病理表型的核心蛋白节点。

单细胞/空间共用

公共数据库验证

将空转发现的基因或细胞特征放入 TCGA 等临床大队列,验证其临床预后与诊断价值。

-

3.3 模块算法清单

微信图片_20260807162144_11_2


四、典型案例:CRC 三级淋巴结构(TLS)分析

以一份结直肠癌(CRC)样本的 Stereo-seq bin20 数据为例,展示 SDAS 如何在肿瘤组织中识别三级淋巴结构(TLS),并探究 TLS 与肿瘤预后的关联。

4.1 时空高级分析 Pipeline 分析流程

微信图片_20260807162033_8_2

4.2 关键步骤与结论

步骤

生物学意义

注释

建立细胞图谱:识别出 B 细胞、T 细胞等免疫亚群

共表达

识别验证TLS结构:找到 3 个 module 与 TLS 空间模式相关,18 个基因落在 TLS_50 标志基因中

邻域

肿瘤微环境:5 个 CN:Epi 富集、免疫富集、Fibro/smoothmuscle、Plasma、Endo

DEG+富集

基因通路:CN2(免疫区)上调补体通路、炎症反应;EMT、糖酵解下调

CCI

细胞相互作用:B-T-NK-ILC-Myeloid 在 CN2 通讯增强;CD86-CD28 在 DC→T 显著

微信图片_20260807162033_9_2

微信图片_20260807162033_10_2


五、算法测评专栏

为帮助研究者在 SDAS 之外也能独立做方法学选择,我们围绕空间组学主流算法开设了"理论 + 实测 + 推荐"的算法测评专栏。已上线的两期:

5.1 细胞注释算法测评

对比 cell2location、RCTD、SPOTlight、Tangram、SCimilarity 等方法在 不同组织类型、不同分辨率下Stereo-seq 数据上的准确性与资源消耗情况。

5.2 空间基因共表达算法测评

对比 Hotspot、hdWGCNA、NeST 等方法在不同组织类型、不同分辨率下Stereo-seq 数据上的准确性与资源消耗情况。


六、下载获取与使用

6.1 SDAS 安装包与 Pipeline

6.2 SDAS 用户手册