(文中引用了一些其他作者的教程,如果是不允许转载的请私信告知我,我会立即删除,谢谢)
velocyto
用于单细胞数据的分析(RNA速率分析)的工具。除了常用的拟时序轨迹分析,RNA速率分析也是一种细胞命运分析的工具。原理不再赘述(主要是咱也不太懂),简言之,细胞中的RNA在转录后会存在剪切成熟的过程,根据对不同成熟度的RNA的分析,确定细胞在发育过程中的先后。
velocyto
有pyhton版本和R版本,velocyto.R
只能用于python版分析后产生的.loom
格式文件的下游处理,并且只能在Linux系统运行(在Windows系统安装了多次失败后的血泪经验)。
1. python版velocyto
的安装
velocyto
python版使用conda
和pip
还是比较顺利的。
1.1 作用
pyhton版velocyto
的作用是将.bam
转换为.loom
格式,.loom
可以用python继续下游分析,也可以使用R进行分析。
1.2 输入输出文件
velocyto
的输入文件包括3个。其中,2个gtf
,rmsk
重复序列的.gtf
文件可以在UCSC下载,教程可以参考这个文章中的“下载特定物种的特殊gtf文件”使用基于python的velocyto软件做RNA速率分析 - 云+社区 - 腾讯云 (tencent.com);
另一个gtf
文件可以在cellranger官网下载。
同时还需要cellranger软件下机的outs
文件夹,该文件夹中包括.bam
文件和其他一些中间文件。
1.3 安装
python版velocyto
的安装根据velocyto官网的教程进行安装即可,Installation Guide — velocyto 0.17.16 documentation。
1.3.1 创建并激活子环境
conda create -n velocyto # python版本应该大于3.6
conda activate velocyto
1.3.2 安装依赖的包
conda install numpy scipy cython numba matplotlib scikit-learn h5py click
1.3.3 pip安装pysam、velocyto
pip install pysam
pip install velocyto
安装成功后,输入velocyto --help
,输出如下,表明安装成功:
~$ velocyto --help
Usage: velocyto [OPTIONS] COMMAND [ARGS]...
Options:
--version Show the version and exit.
--help Show this message and exit.
Commands:
run Runs the velocity analysis outputting a loom file
run10x Runs the velocity analysis for a Chromium Sample
run-dropest Runs the velocity analysis on DropEst preprocessed data
run-smartseq2 Runs the velocity analysis on SmartSeq2 data (independent bam file per cell)
tools helper tools for velocyto
2. R版velocyto.R
的安装
参考的教程:hdf5安装_成功安装velocity.R!
2.1 conda创建一个独立的环境
强烈建议创建一个独立的环境,因为各种依赖包版本不同的话就无法进行下去,要防止各种包和库文件之间的相互影响。
conda create -n R4velocyto
conda activate R4velocyto
2.2 conda安装R及相关R包、库
我使用conda安装了4.0版本的R,其他版本或者最新版本没有尝试。
conda install -c conda-forge R=4.0 # 搜索python的channel进行添加
经验表明,使用conda安装包括devtools
包在内的包是个不错的选择。缺什么依赖包,就安装好了
conda install r-curl
conda install r-rversions
conda install r-credentials
conda install r-httr
conda install r-gert
conda install r-gh
conda install r-rcmdcheck
conda install r-usethis
conda install r-hdf5r
conda install r-devtools
安装库文件,参考了该博主的教程hdf5安装_成功安装velocity.R!
conda install -c conda-forge openmp
conda install -c anaconda boost
2.3 安装velocyto.R
安装依赖包并进行下载velocyto.R。
install.packages("BiocManager")
BiocManager::install('pcaMethods') # BiocManager安装依赖包pcaMethods
library(devtools) # 加载devtools包
install_github("velocyto-team/velocyto.R") # 安装velocyto.R
启动R
,加载包,显示下面的界面,便是安装成功了。
> library(velocyto.R)
Loading required package: Matrix