如何将GEO下载的.rds文件转换为标准Seurat对象？

问题分析 #

你导入的GSE138826的.rds文件并非完整Seurat对象，而是该对象的RNA和SCT两个Assay的集合（以列表形式存储），直接调用head()自然无法得到正常输出。而macOS上出现的内存耗尽错误，大多是R默认内存限制未适配单细胞数据集的大小导致，和设备硬件配置无直接关联。

解决方案 #

1. 手动重建Seurat对象 #

分步提取列表中的Assay数据，重新构建完整Seurat对象：

# 导入列表数据
ref_list <- readRDS("GSE138826_reference.rds")

# 构建RNA Assay
rna_assay <- CreateAssayObject(
  counts = ref_list$RNA@counts,
  data = ref_list$RNA@data
)

# 构建SCT Assay（包含标准化后的scale.data）
sct_assay <- CreateAssayObject(
  counts = ref_list$SCT@counts,
  data = ref_list$SCT@data,
  scale.data = ref_list$SCT@scale.data
)

# 创建基础Seurat对象（以RNA为默认Assay）
ref_seurat <- CreateSeuratObject(counts = rna_assay)

# 将SCT Assay添加到对象中
ref_seurat[["SCT"]] <- sct_assay

# 若原列表包含元数据，可同步添加
# ref_seurat@meta.data <- ref_list$RNA@meta.data

2. 调整macOS的R内存限制 #

macOS下R默认的内存阈值较低，针对大单细胞数据集需要手动扩容：

• 终端启动前配置（推荐）：
  打开终端，先运行以下命令再启动R（将64G替换为你设备可用内存的70%-80%）：

  export R_MAX_VSIZE=64G
  

• R内动态调整：
  针对Seurat依赖的future包内存限制，运行：

  options(future.globals.maxSize = 64 * 1024^3)
  

3. 验证转换结果 #

运行以下命令确认对象构建成功：

# 检查对象类型（应返回Seurat）
class(ref_seurat)
# 查看已加载的Assay（应包含RNA和SCT）
Assays(ref_seurat)
# 查看表达矩阵前几行
head(ref_seurat@assays$RNA@data)

内容的提问来源于stack exchange，提问作者CroixJeremy_2