ARTICLE DETAIL

资讯详情

深耕网站建设与运营推广的一线实战洞察。

Pangolin-NPU 输出解读:12 通道剪接位点评分背后的生物学含义

Pangolin-NPU 输出解读:12 通道剪接位点评分背后的生物学含义 Pangolin-NPU 输出解读12 通道剪接位点评分背后的生物学含义【免费下载链接】pangolin-npu项目地址: https://ai.gitcode.com/atlasleong/pangolin-npu拿到 Pangolin-NPU 的推理结果第一眼看到LOGITS_SHAPE(1, 64, 12)很多人会疑惑这个 12 通道剪接位点评分到底代表什么Pangolin-NPU 是运行在昇腾AscendNPU 上的 RNA 剪接预测模型核心是膨胀残差 1D-CNN 集成输入一段 pre-mRNA 序列A/C/G/U为每一个核苷酸位置输出 12 个数值。读懂这 12 个通道就等于读懂了模型在哪些位置、哪些组织里、以多大强度参与剪接这件生物学大事。本文用真实输出带你逐通道拆解。为什么是 12 个通道从4 个组织 × 3 类评分说起剪接splicing是 pre-mRNA 变成成熟 mRNA 的关键步骤而剪接强度具有明显的组织特异性——同一个位点在心脏里可能是主要剪接位在肝脏里可能几乎不被使用。Pangolin 正是为此而生它同时预测 heart、liver、brain、testis 四个组织的剪接行为见 config.json 中的tissue_names。每个组织输出3 个通道四个组织合计正好 4 × 3 12 个通道通道索引通道名含义取值范围0 / 3 / 6 / 9{tissue}_no_splice该位点在该组织中不参与剪接0~1softmax1 / 4 / 7 / 10{tissue}_splice_site该位点是剪接位点供体/受体的评分0~1softmax2 / 5 / 8 / 11{tissue}_usage该剪接位点的实际使用率0~1sigmoid通道的命名与顺序由output_channels属性明确定义见 modeling_pangolin.py。三类通道的生物学含义评分、分类与使用率这 12 个通道并不是各自独立的随机数它们背后对应着三种不同的生物学问题splice_site与no_splice是一对 softmax 通道两者之和恒为 1。数值接近 1 的splice_site表示模型认为该位置大概率是真实的剪接供体/受体位点。这部分逻辑来自 PangolinPredictionHead 中的 2 通道 softmax 头。usage是 sigmoid 通道即使某个位置是剪接位点也不代表它每次都被使用。usage衡量的是被转录机器实际选中的比例越接近 1 说明该位点越是这条组织里主力的剪接位点。组织间对比的价值同一个位置把 heart、liver、brain、testis 的usage放在一起就能看出组织特异性剪接的完整图景——这正是 Pangolin 相比传统 SpliceAI 最大的进步。从真实输出看 argmax 类别为什么全是 0、6、9在项目交付的确定性测试中输入是 64 个核苷酸的ACGU重复序列输出ARGMAX_CLASS_IDS[6, 0, 6, 6, ...]。如果你把 argmax 当作模型认为最可能的类别会发现结果里几乎只有 0、6、9三种类别 0 heart_no_splice心脏、非剪接类别 6 brain_no_splice大脑、非剪接类别 9 testis_no_splice睾丸、非剪接这其实是一个非常合理的生物学信号ACGUACGU...是一段人工合成的重复序列不含任何真实剪接信号因此四个组织里不参与剪接的通道全部占优。换句话说即使序列是噪声模型也没有胡乱给高分——它知道这里没有剪接位点。输出的经典用法算变异效应 Δ 评分理解了 12 通道的语义Pangolin-NPU 最常见的用途就顺理成章了比较参考序列与变异序列的输出差Δ 评分。输入参考序列得到 12 通道剪接位点评分把某个位置的碱基替换成变异版本如 A→G再跑一次逐通道相减Δ 0表示该变异增强了剪接位点或使用率Δ 0表示削弱。由于 12 个通道覆盖了四个组织的完整剪接行为你一次就能回答这个突变会不会改变某组织中的剪接改变多大这正是临床变异解读中判断剪接致病变异的标准流程。在昇腾 NPU 上复现这 12 个通道Pangolin-NPU 的交付入口是 inference.py它直接在逻辑设备npu:0上完成前向禁止 CPU 回退。核心就三步分词 →PangolinModel前向 → 从probabilities取(1, L, 12)张量见 inference.py。有一个细节值得注意CANN 默认会把 fp32 卷积放到 HF32 降精度通路在近并列位置上可能让 argmax 翻转例如位置 12 的 brain 与 heart 评分只差 4.5e-5。交付代码在首个 NPU 计算前执行torch_npu._C._npu_setOption({ALLOW_CONV_HF32: disable})见 inference.py把误差从 9.3e-5 压到 1.8e-7保证 64 个位置的类别 100% 与 CPU 基线一致。对剪接位点评分这种差之毫厘、谬以千里的任务精度修复不是可选项。常见疑问速答 12 个通道的分数能直接相加吗不能。splice_site/no_splice是成对 softmax组内和为 1usage是独立 sigmoid三者的数值语义不同。ARGMAX_CLASS_IDS只有 12 种取值吗是的它就是把 12 个通道里分数最高的那个通道索引取出来范围恒为 0~11。序列两端信息不足怎么办模型原生支持 10000 核苷酸的上下文转录本两端会以 N 补足不会破坏单点预测。总结Pangolin-NPU 的 12 通道剪接位点评分并不神秘4 个组织 ×剪接位点评分 非剪接评分 使用率。看懂通道顺序就能直接解读 argmax 类别理解 softmax/sigmoid 的差别就能正确比较组织间的剪接强度掌握 Δ 评分法就能把输出用于真实的变异致病性分析。下次再看到(1, L, 12)你就能一眼读出背后的生物学含义了。【免费下载链接】pangolin-npu项目地址: https://ai.gitcode.com/atlasleong/pangolin-npu创作声明:本文部分内容由AI辅助生成(AIGC),仅供参考
返回列表