空间剪接程序分析

目标：在真实 Stereo-seq 时空数据上开发「空间剪接程序分析」pipeline（对齐 Sandberg 实验室 Splice-rs 的方法学框架：split-reads→junction→PSI→伪聚合统计）。
数据：STDS0000058（MOSTA 图谱）E16.5 小鼠全脑，65,303 细胞 × 23,368 基因，单细胞分割（cell_bin），34 种细胞类型注释 + 17 个脑区 + 空间坐标。
方法：① 剪接因子表达提取 → ② 神经元剪接程序（Rbfox/Nova/Celf4）vs 前体剪接程序（Ptbp/Mbnl）双程序定义 → ③ ψ 伪聚合（Splice-rs 式）→ ④ 2×K χ² 统计检验 → ⑤ 空间可视化。

① 17/17 神经剪接因子全部检出，细胞类型特异性符合已知规律：Rbfox1 神经元/非神经元 = 3.8×，Nova1 = 2.9×。
② 剪接程序 ψ 得分呈清晰空间梯度：脑室生发区（VZ）ψ≈0.14（前体程序主导）→ 中脑/后脑分化区 ψ≈0.55-0.63（神经元程序主导）。
③ 统计显著：ψ 高/低 × 17 脑区 χ² = 7951.8（p≈0，df=16）；17 个剪接因子区域分布全部 padj≈0。