Seurat空间转录组分析中Overlay功能在Harmony整合后的异常处理

背景介绍

在单细胞和空间转录组数据分析中，Seurat是一个广泛使用的R语言工具包。其中Overlay功能常用于对空间数据进行区域裁剪，这在分析特定解剖区域时非常有用。然而，当数据集经过不同整合方法处理后，这一功能的稳定性表现存在差异。

问题现象

用户在使用Seurat V5进行空间转录组数据分析时，发现Overlay功能在原始数据上工作正常，但在经过Harmony整合后出现异常。具体表现为：

1. 原始数据上成功执行：

mouse_905_orig[["AC_left"]] <- Overlay(mouse_905_orig_orig[["x95"]], lAC_to_crop, invert = FALSE)

2. Harmony整合后失败：

mouse_905_harmony[["AC_left"]] <- Overlay(mouse_905_harmony[["x95"]], lAC_to_crop, invert = FALSE)

错误信息提示点数组必须包含0或多个元素。

问题排查与解决方案

1. 不同整合方法的比较

测试发现，当使用RPCA方法进行整合时，Overlay功能仍能正常工作：

mouse_905_integrated_rpca[["AC_left"]] <- Overlay(mouse_905_integrated_rpca[["x95"]], lAC_to_crop, invert = FALSE)

这表明问题可能与Harmony整合过程的特定实现有关。

2. 关键发现

通过调整Harmony整合的参数配置，发现问题可以得到解决：

• 问题配置：使用assay = "SCT"
• 解决方案：改用normalization.method = "SCT"

这一调整使得Overlay功能在Harmony整合后的数据上恢复正常。

技术原理分析

Harmony整合与数据结构的改变

Harmony整合过程可能会改变原始数据的空间坐标信息结构，特别是：

1. 坐标精度变化：整合算法可能对坐标值进行归一化或缩放
2. 数据结构重组：整合后的对象可能采用不同的内部存储格式
3. 元数据处理：空间多边形信息在整合过程中可能未被正确保留

SCT标准化方法的差异

assay = "SCT"和normalization.method = "SCT"虽然都使用SCTransform方法，但实现方式不同：

1. assay参数：直接指定使用现有的SCT assay
2. normalization.method参数：在整合过程中重新计算SCT标准化

后者更完整地保持了数据的空间属性信息。

最佳实践建议

1. 整合方法选择：

   • 优先使用normalization.method参数而非直接指定assay
   • 对于空间数据，推荐在整合前测试关键功能

2. 工作流程优化：

   # 推荐方式
   integrated_obj <- IntegrateData(
     anchorset = anchors,
     normalization.method = "SCT",
     verbose = FALSE
   )
   
   # 替代方案：在整合前完成空间操作
   orig_obj[["region"]] <- Overlay(orig_obj[["fov"]], polygon)
   integrated_obj[["region"]] <- orig_obj[["region"]]
   

3. 错误处理：

   • 遇到类似错误时可检查空间坐标的完整性
   • 验证多边形对象与空间数据的坐标系是否匹配

总结

在Seurat空间转录组分析中，整合方法的选择和参数配置会显著影响后续空间操作的功能。特别是使用Harmony整合时，采用normalization.method = "SCT"而非直接指定assay能够更好地保持数据的空间属性，确保Overlay等空间操作正常执行。这一发现为处理类似问题提供了明确的解决方案，也提醒我们在数据分析流程中需要注意方法间的兼容性。