最近有粉丝提到是否可以把多个样本混杂到一起建立一个10X单细胞转录组库进行测序后数据分析,的确是有这样的例子,比如我前些天在Twitter看到的发表在Nephrology Dialysis Transplantation, 的文章, https://doi.org/10.1093/ndt/gfz227 题目是:A single-cell map for the transcriptomic signatures of peripheral blood mononuclear cells in end-stage renal disease 就是这样。
就是10个ESRD病人,10个正常人志愿者的血液,提取PBMC进行10X仪器的单细胞转录组数据而已。值得注意的是,其中研究团队就做了两个10X哦,也就是说他们把 10个ESRD病人混合成为一个样品,10个正常人志愿者也混合成为一个样品。
下面的流程图写的很清楚具体细胞数量,平均检测到的基因数量是1000,每个样品检测到1万个左右的细胞数量,都符合10X仪器的技术水平。
而且,可以看到,分析其实仍然是我们一直讲解的R包及基础流程,分别是: scater,monocle,Seurat,scran,M3Drop 需要熟练掌握它们的对象,:一些单细胞转录组R包的对象 流程也大同小异:
其中一个步骤是判断重要的基因,我分享过:比较5种scRNA鉴定HVGs方法 里面就提到了 。
我昨天就在生信技能树分享了这个观点;表达矩阵数据的免疫环境组成分析可以是标配 反正对bulk的转录组表达矩阵,是可以通过算法,比如CIBERSORT来推断免疫组成的。而本研究也就是对细胞分群后进行注释,然后看看不同亚群的免疫细胞的比例变化。
又或者,既然是纯粹的看免疫细胞的各种比例,其实跑一个流式细胞仪看看就好了,单细胞转录组有点高射炮打蚊子的感觉。不过,反正就2个10X样品,成本也很低,无所谓啦!
因为本研究设计非常简单,就一个变量,就是10 healthy volunteers and 10 patients with ESRD ,所以主要的分析都是围绕这个变量来进行。既然已经定位到了健康人和疾病组的差异细胞群,就可以继续把该群细胞细分,然后具体看某个亚群的表达差异。
其实,也许,早点发出去是最好的策略吧。总比很多文章最后都无法发出去的要好!