我们分享过单细胞转录组下游的降维聚类分群的很多例子，比如：[人人都能学会的单细胞聚类分群注释](https://mp.weixin.qq.com/s/1O1zuwLyM6_W0hZm5I26UA) ，就是使用了SingleR对细胞亚群进行生物学命名！
最近在这个教程下面看到了一下留言，说下载SingleR 的数据库失败。我帮忙看了看， 注意到 Revised: June 14th, 2020的时候SingleR 1.4.1 ，里面使用了下面的代码获取参考数据库:
library(celldex)
hpca.se <- HumanPrimaryCellAtlasData()
hpca.se
也就是说，SingleR 这个包本身并不会自带数据库啦，而是专门的把数据库文件丢给了celldex包。如果你不熟悉bioconductor包，可以先搜索了解一下。发布在bioconductor的包主要是生物信息学相关，在官方可以看到其主要是分成3类：
- 软件方面的包（包括各种芯片数据处理，NGS数据处理，差异分析等等！）
- 注释方面的包(第二类是一系列的基因组注释包，主要是各种ID的转换，kegg或者GO这样的功能注释，还有其它基因信息注释，转录本，外显子起始终止等等)
- 实验数据的包(每一个实验数据包都是一篇优秀的生物信息学分析文章，分析方法，思路都是值得学习的！)
不过，bioconductor除了罗列这3种包，还给了一些其它资源，比如：
-
S4对象的讲解(这个是综合性质的讲解，因为bioconductor系列的包的基础就是一系列对象及函数，需要细致的讲解)
-
AnnotationWorkflow (3) BasicWorkflow (5) EpigeneticsWorkflow (4) GeneExpressionWorkflow (11) GenomicVariantsWorkflow (2) ImmunoOncologyWorkflow (14) ProteomicsWorkflow (2) ResourceQueryingWorkflow (2) SingleCellWorkflow (2) -
其它学习资源收集与翻译( bioconductor牵头的一些网络公开课或者研讨会)
目前这个celldex包并没有发布在bioconductor
安装方法是：
install.packages("remotes")
remotes::install_github("LTLA/celldex")
celldex包里面的更多数据库如下所示：
library(celldex)
ref <- HumanPrimaryCellAtlasData()
ref <- BlueprintEncodeData()
ref <- MouseRNAseqData()
ref <- ImmGenData()
ref <- DatabaseImmuneCellExpressionData()
ref <- NovershternHematopoieticData()
ref <- MonacoImmuneData()
如果你下载失败，就需要考虑切换网络，或者求助海外朋友帮你下载。他需要使用下面的代码：
library(celldex)
ref <- HumanPrimaryCellAtlasData()
save(ref,file = 'HumanPrimaryCellAtlasData.Rdata')
ref <- BlueprintEncodeData()
save(ref,file = 'BlueprintEncodeData.Rdata')
ref <- MouseRNAseqData()
save(ref,file = 'MouseRNAseqData.Rdata')
ref <- ImmGenData()
save(ref,file = 'ImmGenData.Rdata')
ref <- DatabaseImmuneCellExpressionData()
save(ref,file = 'DatabaseImmuneCellExpressionData.Rdata')
ref <- NovershternHematopoieticData()
save(ref,file = 'NovershternHematopoieticData.Rdata')
ref <- MonacoImmuneData()
save(ref,file = 'MonacoImmuneData.Rdata')
如果你还使用的是旧版R包，会出现如下所示的提示：
> library(SingleR)
Attaching package: ‘SingleR’
The following objects are masked from ‘package:celldex’:
BlueprintEncodeData, DatabaseImmuneCellExpressionData, HumanPrimaryCellAtlasData,
ImmGenData, MonacoImmuneData, MouseRNAseqData, NovershternHematopoieticData
得到的文件大小如下：
14M Apr 14 21:51 BlueprintEncodeData.Rdata
292M Apr 14 21:52 DatabaseImmuneCellExpressionData.Rdata
42M Apr 14 21:51 HumanPrimaryCellAtlasData.Rdata
56M Apr 14 21:51 ImmGenData.Rdata
24M Apr 14 21:52 MonacoImmuneData.Rdata
16M Apr 14 21:51 MouseRNAseqData.Rdata
8.6M Apr 14 21:52 NovershternHematopoieticData.Rdata
这些Rdata文件只需要通过网络传递给你，你load进入即可，无需再使用函数实时联网下载数据库啦！