一张热图同时压缩了多种条件
Hi-C矩阵把两个基因组区间的接触计数放进二维坐标。对角线附近通常更亮,是因为线性距离较近的片段更容易发生接触;远离对角线的信号更稀疏。颜色映射还会受到测序深度、分箱大小和显示阈值影响,因此两张颜色相似的图不一定来自可比较的数据。
比较之前应先确认参考基因组、染色体、分辨率、有效读段数量和标准化方法。若一张图按40 kb分箱,另一张按1 Mb分箱,视觉纹理的差异可能主要来自尺度,而不是细胞结构。
原始计数不能直接回答结构问题
接触频率随着基因组距离增加而下降,限制性酶切位点、GC含量、可比对性和片段长度也会制造系统偏差。原始矩阵适合检查数据规模和异常条带,却不适合直接比较不同区域的结构强弱。
O/E矩阵把观察值与距离背景下的期望值比较,矩阵平衡则尝试让行列贡献更一致。两种方法解决的问题不同,结果也不能混用。页面或论文只写“normalized”而不说明方法时,结论应保持谨慎。
颜色需要回到数值和区间
可靠阅读会同时保存颜色范围、缺失值规则和矩阵分辨率。为了突出局部结构而截断极端值是合理的可视化选择,但必须说明,否则读者会误以为不同样本使用相同尺度。
HK-BEUP传输矩阵时应把绘图参数、处理版本和原始文件摘要一起交付。接收端重新绘图后若出现差异,团队才能判断问题来自文件、标准化还是颜色设置。
一个更可靠的对照方式
若要比较两种条件,可以先在相同参考版本与分辨率下重算矩阵,再分别查看有效配对数量、距离衰减曲线和矩阵稀疏度。随后固定颜色范围绘图,并在代表区间之外抽查多个染色体。只有局部结构在重复样本中稳定出现,才值得进入后续机制解释。
侧重展示的局部热图还应配套全局质量摘要。否则一个视觉鲜明的区域可能掩盖整份矩阵覆盖不足、重复率过高或染色体命名错误。
什么情况下不宜继续比较
如果原始读取、参考版本或标准化记录缺失,重新绘图只能提高视觉一致性,不能补回方法依据。此时更适合把结果标记为探索性资料,并优先寻找原始数据或独立实验。
矩阵阅读的目标不是找到最亮的方块,而是判断信号在什么条件下产生、能否重复,以及换一种处理方法后是否仍然存在。