能发CNS的单细胞转录组三大高阶分析

任何技术在流行的时候，都是一大堆人开发一大堆工具，但是到最后，大浪淘沙，只有一两个优秀工具被大家广为传颂！比如450K等甲基化芯片数据处理，就只剩下了ChAMP等：

只剩下了ChAMP

单细胞数据处理也是如此，在单细胞领域近5年超一千个软件工具算法被开发 ，详见近期Genome Biology (2021) 杂志的一个文章：《Over 1000 tools reveal trends in the single- cell RNA-seq analysis landscape》系统性的梳理了2016-2021这5年间，位于https://www.scrna-tools.org/网页的超一千个软件工具算法。

虽然说单细胞工具一直在增加，最基础的往往是降维聚类分群，参考前面的例子： 人人都能学会的单细胞聚类分群注释 ，这个领域基本上是Seurat的天下。

但是如果要提高文章档次，就需要一些高级分析了，其实它的三大高阶分析一直没有变化。

• 没想到自己会放弃conda（docker镜像的pyscenic做单细胞转录因子分析）
• 使用pyscenic做转录因子分析
• 使用基于python的velocyto软件做RNA速率分析
• CellPhoneDB的单细胞通讯结果的可视化之气泡图
• CellPhoneDB的单细胞通讯结果的理解
• 直接为CellPhoneDB创建一个独立的conda环境
• 把Seurat对象里面表达量矩阵和细胞表型信息输出给CellPhoneDB做细胞通讯
• 各个单细胞亚群特异性的转录因子热图
• 简单直接的拟时序分析方法，R包SCORPIUS推荐
• 把基因表达量画在拟时序结果图上
• 拟时序分析就是差异分析的细节剖析

比如张泽民课题组的泛癌层面NK单细胞水平分群，于2023年9月发表在CELL的《A pan-cancer single-cell panorama of human natural killer cells》，也是如此的单细胞转录组三大高阶分析：

单细胞转录组三大高阶分析