Kinetics-seq在单细胞水平同时测量RNA合成与降解速率
单细胞 RNA 测序(scRNA-seq)有一个它自己解决不了的问题:裂解时刻的 RNA 存量是合成速率和降解速率共同塑造的,这两者在丰度快照里无法拆开。一条 RNA 多,可能因为正在高速转录,也可能因为降解极慢。数字相同,背后的调控程序未必一样。
Kinetics-seq 把体内代谢标记(in vivo metabolic labeling)嵌进 scRNA-seq 流程来打开这个黑箱。带标记原子的核苷酸前体引入生物体后,新合成的 RNA 分子会掺入标记,存量 RNA 不会;裂解后同时测总 RNA 和新合成 RNA,就能反推每条转录本的合成速率(synthesis rate)和降解速率(degradation rate)。加上原有的丰度(abundance)和周转率(turnover),Kinetics-seq 在单细胞分辨率上可以同时给出四个参数,传统 scRNA-seq 只给其中一个。(PMID: 42771508)
批量水平的代谢标记测序早已存在,Kinetics-seq 的推进点是把它落到单细胞分辨率且在组织原位完成。不过摘要没有披露标记试剂的具体类型、时间窗口长度和标记效率,也没有说用的是哪种肿瘤和多大样本量。动力学参数估计的精度和这些细节高度相关,方法的适用边界目前难以判断。
肿瘤细胞群内部的动力学层
丰度分析能区分细胞类型,这是 scRNA-seq 的基本功。加入合成速率和降解速率之后,在同一个按丰度分好的细胞群内部,还能再分出一层异质性。
两个细胞的某条基因丰度相近,一个是高合成 + 快降解,另一个是低合成 + 慢降解——丰度数字一样,调控状态是两套。通过对合成速率和降解速率的联合建模(joint modeling),研究者在肿瘤细胞之间观察到显著的动力学多样性(pronounced kinetic diversity),不同基因也倾向于采用不同策略:有的靠提高合成维持高表达,有的靠减缓降解达到同样结果。(PMID: 42771508)
这部分结论目前局限在肿瘤组织里。体内代谢标记在动物模型里操作相对直接,但能否适用于人类临床样本,摘要没有交代。放到正常组织或非肿瘤病理状态下,动力学图谱是否呈现类似规律,也有待独立数据。
合成速率筛出的三条通路
用 RNA 合成速率而非丰度作为排序依据,筛出转录活跃的基因子集再做富集,结果在三条通路上找到了更清晰的信号:雌激素早期应答(estrogen response early)、PI3K-AKT-mTOR 信号通路,以及上皮-间质转化(epithelial-mesenchymal transition,EMT)。
这三条通路与肿瘤进展的时间关联,在动力学分析里比丰度分析更强。(PMID: 42771508)
不是说丰度分析找不到这些通路,而是动力学指标的变化和肿瘤进展节奏更对齐。如果合成速率在某个阶段先响应、丰度后变,那基于丰度设计的分析就会错过这个窗口。这对早期生物标志物识别有直接含义,但具体影响要靠独立临床队列来验证。
另外两篇:
血清素(serotonin)神经元发育方面(PMID: 42771501),这项小鼠研究用深度学习流水线构建了出生后全脑血清素轴突密度的四维图谱,发现两层并存的调控:脑区特异性的轴突生长与精修,以及出生早期的全脑同步化浪涌。单细胞转录组揭示血清素神经元亚型多样性在胚胎期已确立,亚型分类在小鼠与人之间高度保守;出生后第 4 天存在一个跨亚型协同的接线分子(wiring molecule)表达峰,与随后全脑支配浪涌时间上吻合。人类部分为转录组比较,无人体发育影像验证。
黑色素瘤免疫亚型与药物重定向方面(PMID: 42777131),用 432 例 TCGA 皮肤黑色素瘤 RNA-seq 数据做免疫冷热分型,耗竭型(D,n = 169)与免疫富集型(IE + IE/F,n = 184)对比得到 300 上调 / 300 下调差异特征,经 LINCS L1000 连接作图(connectivity mapping)筛选候选化合物。DNA 甲基转移酶(DNMT)抑制剂信号最一致,fdcyd、阿扎胞苷(azacitidine)和 RG-108 三个化合物在两折交叉验证里均保持较高排名(decitabine 不稳健),但类别富集经 FDR 校正后未达显著(q = 0.126)。IDO 抑制剂趋于负向(q = 0.095),MEK 抑制剂无明显模式(q = 0.698)。冷型特异上调模块的折间重叠仅 25.7%,远低于热型的 88.7%,说明耗竭表型的转录特征本身不够稳定。整体属于探索性信号,现有一项阿扎胞苷 + 卡铂联合抗 PD-L1 的 II 期试验数据与此收敛。
- System-level synchronization and subtype-specific modules govern brain-wide serotonin axon innervation.,Proceedings of the National Academy of Sciences of the United States of America(PMID: 42771501)
- Kinetics-seq enables comprehensive profiling of single-cell RNA kinetics in vivo to reveal dynamic tumor heterogeneity.,Proceedings of the National Academy of Sciences of the United States of America(PMID: 42771508)
- Transcriptomic immune subtyping and connectivity mapping nominate DNA methyltransferase inhibitors for immunedepleted melanoma.,Melanoma research(PMID: 42777131)