Scanpy项目中的Visium HD空间转录组数据支持解析

背景介绍

Scanpy作为单细胞转录组分析的重要工具，近期在社区讨论中关于其对10X Genomics Visium HD空间转录组数据的支持问题引起了广泛关注。随着10X Genomics Space Ranger软件升级至v3.0版本，其输出文件格式发生了变化，导致Scanpy原有的读取功能不再适用。

技术解决方案

针对这一问题，scverse生态系统中的spatialdata-io项目已经实现了对Visium HD数据的读取支持。该解决方案通过以下方式工作：

1. 数据读取：spatialdata-io提供了专门的visium_hd读取器，能够解析Space Ranger v3.0输出的原生格式文件（包含多种文件扩展名，其中部分为Zarr格式）

2. 数据结构转换：读取后的数据会被转换为SpatialData对象，用户可以直接从中获取AnnData对象进行后续分析

3. 与Scanpy集成：获得的AnnData对象可以无缝接入Scanpy的分析流程，进行预处理、质量控制和聚类分析等标准操作

使用建议

对于需要使用Visium HD数据的用户，推荐采用以下工作流程：

1. 通过spatialdata-io的visium_hd函数读取原始数据
2. 将返回的SpatialData对象转换为AnnData格式
3. 使用Scanpy进行后续的转录组分析
4. 如需保存处理后的数据，可将其转换为SpatialData的Zarr格式进行存储

注意事项

1. 该解决方案目前不支持直接从Zarr格式读取Visium HD数据，而是处理Space Ranger的原生输出格式

2. 在某些情况下可能需要明确指定dataset_id参数，特别是当自动推断失败时

3. 对于只想进行基本转录组分析而不需要复杂空间信息的用户，可以直接使用提取出的AnnData对象

未来展望

随着空间转录组技术的快速发展，Scanpy团队正在考虑将原有的空间数据IO功能迁移到专门的spatialdata项目中，以提供更专业、更灵活的空间数据分析支持。这种模块化的设计将使Scanpy能够更专注于核心的单细胞分析功能，同时通过与其他专业工具的集成来扩展其应用范围。

对于需要同时分析基因表达和空间信息的用户，建议关注scverse生态系统中spatialdata和spatialdata-io项目的持续发展，这些工具将提供更全面的空间数据分析解决方案。