R中CreateSeuratObject问题：首列为基因名时显示0个特征跨X样本

核心逻辑 #

CreateSeuratObject的counts参数有严格的输入要求：

• 必须是全数值型的矩阵/数据框，所有值对应基因表达计数
• 行对应基因，行名为基因ID/名称
• 列对应细胞，列名为细胞barcode/编号
  不需要将基因名作为单独的一列放在计数矩阵中。

操作问题点 #

你处理后的X数据框，第一列是字符型的基因名，其余列是数值型的细胞计数，同时merge操作会重置数据框的行名，此时X的行名是默认的1、2、3...数值索引，不符合counts参数的输入要求：

1. 存在非数值的基因名列，Seurat解析counts时会自动过滤非数值内容，无法识别到有效的表达特征，因此返回0 features的结果
2. 没有将基因名赋值给行名，即使后续删除基因名列能成功创建对象，基因名也会是默认的数值索引，不便于后续分析

正确操作 #

处理到X首列为基因名时，先将基因名赋值为行名，再删除基因名列即可：

# 将首列基因名设置为行名
rownames(X) <- X[,1]
# 删除首列基因名，保留纯数值计数矩阵
X <- X[,-1]
# 再创建Seurat对象即可同时保留正确的基因名
X_seurat <- CreateSeuratObject(counts = X, project = "X_seurat", min.cells = 5)

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