空间质量指标为什么也要画回组织坐标
质量指标不是脱离位置的一列数字。把UMI、检测基因、线粒体比例、阈值命中与配准结果画回组织,才能看见边缘、批次、图像和真实组织梯度。
空间转录组分析常先查看每个点位的UMI数、检测基因数和线粒体比例,再用统一阈值标记低质量点。若这些数字只留在直方图和表格中,分析人员看不到被标记点是否刚好集中在真实组织区域、切片边缘、拼接缝或配准错误附近。
空间数据多了一层不能丢掉的证据:位置。质量指标既可能反映捕获、测序、分割和成像问题,也可能随细胞密度、组织结构与转录活动自然变化。删除前把指标画回组织坐标,是区分技术异常和生物梯度的必要步骤。
质量检查从图像与坐标开始
表达矩阵带有坐标,不代表坐标一定正确。10x Visium成像指南要求组织和定位框清晰可见、对焦正确、曝光适当,拼接不能扭曲定位框;定位框随后用于Space Ranger对齐图像。
欠曝会隐藏组织细节,过曝可能让信号饱和;低分辨率、失焦、白平衡或阴影修正问题都会改变形态。拼接缝如果拉伸定位框,点位与组织的对应可能整体偏移。此时后续任何“空间异常”都可能建立在错误坐标上。
所以第一张质量图不是UMI直方图,而是原始组织图、定位框与点位覆盖。确认图像文件对应正确玻片和捕获区,方向一致,组织边界与点位对齐。发现明显配准问题时,应先修正图像和坐标流程,不能靠表达阈值删除错位区域。
文件身份要先于图形解释
10x成像指南建议文件名同时包含玻片序号和捕获区标识。这个细节防止图像交换:表达文件与另一张切片图像即使都能打开,叠加后也会产生看似复杂的空间不一致。
分析清单应记录玻片、捕获区、图像、FASTQ、Space Ranger运行和输出目录的稳定关系,并保存哈希。图像重新导出或手动配准时生成新版本,不覆盖原图。
若多个荧光通道分别导出,还要确认位深、尺寸、对齐与格式一致。合并RGB图便于查看,却会失去独立调整每个通道的能力;原始单通道资料应按实验设计保留。
全局分布只能发现候选异常
OSTA质量控制章节列出常见点位指标:文库大小、检测基因数、线粒体转录比例和每点细胞数。低文库或低检测基因可能与捕获效率、细胞损伤或反应效率有关;高线粒体比例可能提示细胞损伤;异常细胞数可能涉及分割。
这些解释都应写成“可能”。同一个低文库值可以来自技术损失,也可以来自本来转录较少或细胞稀疏的区域。直方图能显示尾部和离群值,却没有组织语境。
第一轮可用宽松规则标记候选,不立即删除。保存每项指标原值、变换方式、阈值、规则版本和命中原因。多个指标共同异常与单一指标轻微越界,应保持不同标签。
组织结构会混杂质量指标
OSTA以人脑灰质和白质为例:神经元胞体较多的灰质可以检测到更多基因和转录本,白质则可能有较低文库和更高线粒体比例。如果不了解组织结构,统一阈值可能把白质误判为低质量。
这个例子不能外推为所有组织中的固定方向,却揭示一般机制:细胞组成、密度和形态区域会共同改变QC指标。空间点不是从同一分布随机抽出的独立样本,它们邻近并共享组织背景。
因此,阈值选择前要查看已知形态或注释区域。若被标记点沿完整解剖结构连续分布,而不是孤立散点,生物混杂的可能性上升。此时更适合保留并在下游模型中处理,或使用区域内比较,而不是全局删除。
把每项指标画回组织
至少绘制文库大小、检测基因数、线粒体比例、组织内外标记和细胞数的空间图。连续色标显示梯度,另用二值图显示阈值命中。两种图不能互相替代:色标保留程度,二值图显示实际删除边界。
颜色尺度必须固定并记录。每个样本自动拉伸到自身最小与最大值,会让跨切片图看起来相似,却隐藏绝对差异。极端值可以为显示而截断,但原值和截断分位必须保存。
将阈值命中与组织图并排,而不是用彩色点完全遮住形态。检查异常是否集中在组织边缘、裂隙、折叠、脱落、染色不均、拼接缝、定位框附近或某一批次区域。
Space Ranger汇总值也要回到位置
10x的Space Ranger Web Summary技术说明把“组织下点位读取比例”与图像方向、组织识别和环境RNA等问题联系起来。单一百分比告诉整体偏离,却不能说明问题位于哪一处。
同一文档提供按UMI着色的组织点位图,并建议观察UMI与组织形态是否一致。组织细胞密集或高表达区域出现较高UMI可能合理;若信号跨出组织、集中在不相符区域或呈现异常条带,则需要回查试剂、组织、图像和配准。
中位UMI、每点检测基因和测序饱和度也受样本类型、RNA质量、深度和实验目标影响。把厂商示例中的期望值当成跨组织硬阈值,会忽略这些条件。
聚类与形态不一致是调查入口
Space Ranger文档指出,UMI或聚类空间分布与组织形态不一致,可能与试剂不均、组织脱落、图像对齐或图像互换有关。这里的关键词是“可能”:形态不一致提供调查方向,不是自动失败判定。
先检查文件身份与配准,再看异常是否沿试剂覆盖、切片损伤或拼接区域出现。之后检查表达标志与已知组织结构。若形态与表达确实反映不同层面,生物差异也可能合理。
聚类本身来自表达相似性和参数设置,不能反过来作为图像绝对真相。多个参数下都稳定且与独立形态证据一致,才提高解释可信度。
全局阈值为什么会误删整片区域
OSTA示例中,文库阈值设得较高后,被标记点沿已知皮层层次集中。若只看总体直方图,阈值似乎合理;画回坐标后,才看见规则正在选择一个生物区域。
这类误删会改变下游组成。某一层或细胞群被系统性移除,差异表达、空间域和邻域关系都会随之变化。删除点数量不多,也可能因为位置集中而产生很大偏差。
所以每次改阈值,都要报告命中总数和空间分布。只说“删除3%低质量点”不够,还要说明这3%是随机散布、集中边缘,还是覆盖特定组织结构。
局部邻域提供另一种受控比较
全局阈值把每个点与整个样本比较,隐含各处共享分布的假设。OSTA还展示局部离群方法:把每个点的文库、检测基因和线粒体比例与空间邻居比较,从而寻找相对周围异常的点。
如果一个区域整体较低,但内部相互一致,局部方法可能保留它;若单个点明显偏离邻居,则更容易被标记。这适合组织结构强、全局差异大的场景。
局部方法也不是答案。邻居数量、六角或方格点阵、组织边界和分布变换都会影响结果。跨越不同组织区域的邻域可能再次混入生物差异。应把全局与局部结果并排,检查两者为何不同。
切片边缘既可能有技术问题,也可能有真实结构
边缘容易发生组织损伤、脱落、低细胞密度或配准不稳定,因此常出现低计数。边缘同时也可能是有意义的解剖界面。不能建立“距离边缘小于某值就删除”的通用规则。
先看原始图像是否有折叠、破裂、缺失或焦点问题,再看多个QC指标是否共同异常。若只有文库较低而形态完整、邻近点连续且标志表达合理,更适合保留或单独分析。
对于成像型平台,还要检查视野边界和分割。OSTA成像QC章节说明视野差异、边界细胞和分割可能影响计数与形态指标。不同平台产生异常的机制不同,不能把Visium点位规则直接套到单细胞分割数据。
阈值应与样本和问题一起记录
固定阈值便于重复,但如果来自另一组织或平台,就可能错误。数据驱动阈值也会随样本组成改变。无论采用哪种方式,都要记录依据、变换、指标方向和预期影响。
若多张切片构成实验,先逐样本查看,再决定是否共享规则。统一规则便于比较,却可能让质量较弱样本过度删除;每片单独调参又可能引入人为差异。可以先设共同原则,再为特殊样本记录偏离理由。
阈值不是最终生物结论。被标记点保留在对象中,用状态字段区分,而不是立刻从所有文件永久删除。这样可以重跑并检查结果对规则是否敏感。
删除前做过滤敏感性分析
至少保留三个版本:未过滤、仅去除强技术失败、按完整规则过滤。比较每个版本的点位数、组织覆盖、主要空间域、标志基因和关键结论。若结论只在最严格阈值下出现,应把依赖关系明确报告。
检查被删除点在组织图中的覆盖面积,而不只看数量。少量连续点可能切断空间邻接,影响边界或轨迹。下游方法使用邻接图时,更要确认过滤没有制造人工孔洞。
所有版本使用相同坐标基线和图像,过滤规则通过状态字段重建。不要为每版重新导出后丢失点位对应,避免后续无法判断差异来自数据还是坐标变化。
先做坐标与矩阵的一致性验收
在计算阈值前,检查表达矩阵中的点位标识能否一一映射到坐标表,坐标数量、组织内标记和图像尺寸是否符合本次运行记录。重复标识、缺失坐标或未知捕获区都应先解决,不能在绘图函数中静默丢弃。
随机抽取若干定位框附近、组织中心和边缘点,确认图上的位置与表中坐标方向一致。水平翻转、上下颠倒或缩放比例错误有时仍会产生一张“像组织”的图,但与真实位置不对应。抽样核对能在复杂分析前发现这种系统性错误。
所有修正都生成新坐标版本并保存变换参数。手动拖动图像后只保留最终截图,会让分析无法重现,也无法判断某个异常图案是在配准前还是配准后出现。
为每一种异常图案指定调查顺序
孤立低文库点先比较邻域和原图局部,再查对应孔洞、气泡或组织损伤;连续边缘带先检查组织厚度、焦点和组织内识别;跨越整片的条带或拼接缝则优先回查扫描、曝光、试剂分布和图像拼接。
若高线粒体比例与低文库、低检测基因同时出现,技术损伤解释通常比单一指标更强,但仍要观察其是否覆盖完整组织区域。若只有一个指标异常,应先确认计算分母、基因注释和平台输出定义。
UMI或聚类跨出组织边界时,先查组织掩膜、图像方向和环境RNA,而不是立刻给外部点赋予新生物含义。相反,形态完整区域出现稳定表达结构时,也不应仅因与预期图谱不同而删除。
批次比较要保留各自的空间背景
多张切片的总体分布不同,可能来自测序深度、组织组成、处理时间或真实条件。先在每片内部画图与标记,再建立跨片汇总。把所有点合并后计算一个阈值,容易让高深度切片成为标准,系统性移除另一片的正常区域。
跨片比较可报告每片中位数、分位数、组织覆盖和异常空间形态,同时保留统一的计算定义。若必须使用共同门槛,说明门槛如何从实验设计或验证样本得到,并展示每片命中区域。
技术重复应检查异常图案是否复现。相同解剖区域在重复中出现类似梯度,支持生物或稳定流程解释;仅一片沿拼接缝或试剂方向出现,更值得调查技术过程。复现仍不是因果证明,但能缩小排查范围。
过滤决策需要机器可读的审计表
每个点保留原始指标、全局规则命中、局部规则命中、图像观察、最终状态和理由代码。理由使用有限词表,例如图像错位、组织外、强多指标异常、待复核或保留生物区域,同时允许附加人工说明。
审计表与图像、参数和软件版本一起进入结果包。下游分析只按最终状态选择点位,不直接重算另一套隐含阈值。这样同一过滤版本能在差异表达、聚类和空间邻接中一致复用。
人工复核者若改变状态,要记录前后值和证据。质量控制不是无法追踪的鼠标选择,而是一组可以重放、比较和撤销的决策。
用负控与实验记录增强解释
空白区域、无组织点或实验设计中的负控可以帮助观察背景,但其定义依平台而异。背景升高若与组织外区域、特定试剂方向或某批次同时出现,比单点高值更能支持流程调查。
实验记录要与空间图对时:染色、成像、透化、文库制备和测序中出现的延迟、设备告警或样本损伤,都可成为异常图案的候选解释。记录用于提出和排除假设,不可为了符合日志而忽略不一致数据。
没有负控或完整日志时,应降低结论强度,保留问题点并报告不确定性。删掉无法解释的数据会让结果看起来更整洁,却不会补回缺失证据。
让图表能够复核,而不只是好看
每张空间质量图都标注样本、指标、单位、变换、颜色范围、阈值和点位状态。图例不使用“好”“坏”代替数值,避免观看者把连续差异误解成已经验证的类别。用于显示的截断与用于判定的原始值必须分开。
导出静态图时同时保存生成代码和数据摘要。报告中的图若经过裁切、翻转或缩放,应注明与分析坐标之间的变换;否则下一位复核者可能在另一方向寻找同一异常。
关键结论至少配一张未过滤图、一张标记图和一张过滤后图,并在相同视野与颜色范围下比较。这样可以直接看见规则移除了什么空间结构,而不是只在表格中看到点位数量减少。
把不确定点留给后续证据
并非所有候选异常都要在第一次分析中二选一。证据冲突的点可以标为待复核,在主要分析中做包含与排除两套结果。若核心结论在两套结果中一致,过滤争议的影响较小;若结论改变,就应把不确定性写入报告。
后续可利用重复切片、独立染色、病理复核或平台建议补证。新增证据只更新状态与理由,不改变原始指标和坐标。保留这条审计链,才能知道结论变化来自新证据,而不是无记录地调高或调低阈值。
一套删除前审查顺序
先冻结图像、坐标、表达矩阵和文件关系。检查玻片与捕获区身份、方向、焦点、曝光、拼接和定位框。然后计算每点QC指标,查看全局分布,但只做候选标记。
接着把原值、连续色标和阈值命中分别画回组织,与形态、边缘、批次和已知区域比较。对明显空间混杂的指标,尝试区域内或局部邻域方法,并记录邻域定义。
最后检查多指标一致性与实验记录,决定保留、标记、重算配准或删除。运行过滤前后敏感性分析,保存所有状态、阈值和软件版本。删除不应是清理表格的第一步,而是证据整合后的最后决定。
空间点文库大小、检测基因数和线粒体比例既可能反映技术质量,也会随真实组织区域变化。10x要求把UMI与聚类叠加到组织图像,并检查定位框、焦点、拼接和配准。OSTA则展示了全局阈值误删和局部邻域比较的边界。
空间图案不能单独证明技术失败或生物机制。把指标画回坐标的价值,是让错误的统一阈值无处隐藏,也让真正局部异常获得可调查的位置。质量控制保留空间证据,后续生物解释才不会从一张被静默裁剪的组织图开始。
本文核对资料
- OSTA作者组/Bioconductor生态:《Orchestrating Spatially Resolved Transcriptomics Analysis: Quality control》,当前页面交叉核对
- 10x Genomics:《Interpreting Space Ranger Web Summary Files for Visium Spatial Gene Expression for FFPE Assay》,2022-02
- 10x Genomics:《Visium Spatial Gene Expression Imaging Guidelines》,2023-09-26
资料来源
- OSTA作者组/Bioconductor生态:《Orchestrating Spatially Resolved Transcriptomics Analysis: Quality control》,发布或更新于 2026-08-12
- 10x Genomics:《Interpreting Space Ranger Web Summary Files for Visium Spatial Gene Expression for FFPE Assay》,发布或更新于 2022-02-01
- 10x Genomics:《Visium Spatial Gene Expression Imaging Guidelines》,发布或更新于 2023-09-26