单细胞数据“质检员”指南：拿到表达矩阵后，你的第一件事应该是检查这些

单细胞数据质检实战指南：从表达矩阵到可靠分析的五大检查点

当你第一次拿到单细胞RNA测序的表达矩阵时，那种兴奋感可能让你想立刻开始聚类分析和可视化。但作为一名严谨的研究者，按下暂停键进行系统质检(QC)才是明智之举。我曾见过太多案例因为跳过QC步骤而导致后续分析出现偏差——某个异常细胞群后来被证明是死细胞聚集，或者"差异表达基因"其实是线粒体污染的信号。这份指南将带你像专业质检员一样，用Seurat工具系统检查表达矩阵中的关键指标。

1. 理解表达矩阵的基本结构

在开始任何分析前，我们需要明确表达矩阵的构成。典型的单细胞表达矩阵是一个基因×细胞的二维表格，其中的数值代表每个基因在每个细胞中的UMI计数。这种稀疏矩阵通常有以下特征：

• 行(rows)：代表基因(通常用官方基因符号如TP53)
• 列(columns)：代表细胞(通常用细胞条形码如AAACCTGAGCTAACGT-1)
• 值(values)：UMI计数(通常为整数，0表示未检测到)

# 查看表达矩阵基本信息的R代码示例 dim(expression_matrix) # 显示基因和细胞数量 head(rownames(expression_matrix)) # 显示前几个基因名 head(colnames(expression_matrix)) # 显示前几个细胞条形码

常见问题警示：

• 如果矩阵维度显示只有几十个基因或细胞，可能数据加载出错
• 如果基因名显示为ENSG...编号，需要考虑是否转换为更易读的基因符号
• 如果数值包含小数，可能不是原始计数矩阵而是经过标准化处理

专业提示：始终保留原始计数矩阵的备份，所有转换和标准化都应在新对象中进行

2. 细胞层面的基础QC指标

2.1 检测每个细胞的基因数量

每个细胞检测到的基因数(nFeature_RNA)是评估细胞质量的首要指标。健康的细胞通常检测到数千个基因，而低质量细胞或空微滴(empty droplets)则基因数明显偏少。

典型阈值参考：

# 计算并可视化每个细胞的基因数量 library(Seurat) seurat_obj <- CreateSeuratObject(counts = expression_matrix) hist(seurat_obj$nFeature_RNA, breaks = 50, main = "Distribution of Genes per Cell", xlab = "Number of Genes")

2.2 检测每个细胞的UMI总数

每个细胞的总UMI计数(nCount_RNA)反映测序深度。异常高或低的UMI数可能分别表示双细胞(doublets)或低质量细胞。

异常模式诊断：

• UMI数极低：可能来自死细胞或空微滴
• UMI数异常高：可能是多个细胞被捕获在一个微滴中(doublets)
• 双峰分布：可能混合了不同细胞类型或质量群体

2.3 线粒体基因比例

线粒体基因占比(percent.mt)是评估细胞完整性的关键指标。高比例通常表明细胞膜受损，RNA从线粒体泄漏。

# 计算线粒体基因比例 seurat_obj[["percent.mt"]] <- PercentageFeatureSet( seurat_obj, pattern = "^MT-") # 人类用MT-，小鼠用mt- # 可视化QC指标关系 FeatureScatter(seurat_obj, feature1 = "nCount_RNA", feature2 = "percent.mt")

线粒体基因阈值建议：

• 一般样本：<10-20%
• 敏感样本(如神经元)：可能需要放宽至25%

3. 基因层面的QC检查

3.1 检测基因在细胞中的表达频率

有些基因在极少数细胞中表达，可能是测序噪音或低质量转录本。我们可以计算每个基因在多少细胞中表达(>0 UMI)。

# 计算基因表达频率 gene_detection_rate <- rowSums(expression_matrix > 0) / ncol(expression_matrix) # 可视化 hist(gene_detection_rate, breaks = 50, main = "Gene Detection Rate Across Cells", xlab = "Fraction of Cells Expressing the Gene")

值得关注的基因类型：

• 在<5个细胞中表达的基因：考虑过滤
• 在>95%细胞中表达的基因：可能是管家基因

3.2 识别可能的技术干扰基因

某些高表达基因可能来自实验技术而非生物学信号：

• 核糖体蛋白基因：RPL/RPS开头的基因
• 热休克蛋白基因：HSP家族
• 线粒体基因：MT-开头的基因

# 检查核糖体蛋白基因表达比例 seurat_obj[["percent.rb"]] <- PercentageFeatureSet( seurat_obj, pattern = "^RP[SL]") # 匹配RPL或RPS开头的基因

4. 样本层面的QC与批次检查

4.1 检查样本间细胞数量平衡

如果实验包含多个样本，需要确认各样本贡献的细胞数量是否均衡。

# 假设metadata中有sample_id列 table(seurat_obj@meta.data$sample_id) # 可视化 library(ggplot2) ggplot(seurat_obj@meta.data, aes(x = sample_id)) + geom_bar() + labs(title = "Cell Counts per Sample")

异常情况处理：

• 某个样本细胞数极少：考虑技术问题或重新处理
• 样本间细胞数差异极大：可能需要下采样平衡

4.2 检测批次效应

即使使用相同protocol，不同批次的数据也可能存在系统性差异。早期发现批次效应有助于后续校正。

# 使用PCA初步检查批次效应 seurat_obj <- NormalizeData(seurat_obj) seurat_obj <- FindVariableFeatures(seurat_obj) seurat_obj <- ScaleData(seurat_obj) seurat_obj <- RunPCA(seurat_obj) DimPlot(seurat_obj, reduction = "pca", group.by = "batch_id") # 假设metadata中有batch_id列

5. 综合QC与数据过滤策略

5.1 建立多指标过滤标准

结合前述指标，制定适合你数据集的过滤标准。例如：

# 创建过滤逻辑 keep_cells <- seurat_obj$nFeature_RNA > 500 & seurat_obj$nFeature_RNA < 6000 & seurat_obj$percent.mt < 15 # 应用过滤 seurat_obj_filtered <- subset(seurat_obj, subset = keep_cells)

5.2 过滤前后的比较

记录过滤前后的数据变化至关重要：

5.3 处理特殊情况的技巧

• 冷冻样本：线粒体基因比例可能更高，需调整阈值
• 小型细胞(如血小板)：基因检测数可能较低
• 高转录活性细胞(如浆细胞)：UMI数可能异常高

# 针对特殊细胞类型的自适应过滤 if(sample_type == "frozen"){ mt_threshold <- 25 } else { mt_threshold <- 15 }

完成这些QC步骤后，你的表达矩阵已经准备好进行后续的标准化和分析了。记住，好的QC不是寻找"完美"的数据，而是理解数据的局限并在分析中适当考虑这些因素。