velocyto-scRNAseq速率的计算
velocyto-scRNAseq速率的计算
介绍：https://cloud.tencent.com/developer/article/1607921
安装：注意：velocyto 需要 python 版本>=3.6.0
git clone https://github.com/velocyto-team/velocyto.py.git
cd velocyto.py
pip install -e . # note the trailing dot其他安装方式https://velocyto.org/velocyto.py/install/index.html#install
使用方法以10X为例
1.Running the CLI
Usage: velocyto run10x [OPTIONS] SAMPLEFOLDER GTFFILE
Runs the velocity analysis for a Chromium 10X Sample
10XSAMPLEFOLDER specifies the cellranger sample folder
GTFFILE genome annotation file
Options:
-s, --metadatatable FILE Table containing metadata of the various samples (csv fortmated rows are samples and cols are entries)
-m, --mask FILE .gtf file containing intervals to mask
-l, --logic TEXT The logic to use for the filtering (default: Default)
-M, --multimap Consider not unique mappings (not reccomended)
-@, --samtools-threads INTEGER The number of threads to use to sort the bam by cellID file using samtools
--samtools-memory INTEGER The number of MB used for every thread by samtools to sort the bam file
-t, --dtype TEXT The dtype of the loom file layers - if more than 6000 molecules/reads per gene per cell are expected set uint32 to avoid truncation (default run_10x: uint16)
-d, --dump TEXT For debugging purposes only: it will dump a molecular mapping report to hdf5. --dump N, saves a cell every N cells. If p is prepended a more complete (but huge) pickle report is printed (default: 0)
-v, --verbose Set the vebosity level: -v (only warinings) -vv (warinings and info) -vvv (warinings, info and debug)
--help Show this message and exit.需要准备的文件：
1.下载genome annotation file
Download a genome annotation (.gtf file) for example from GENCODE or Ensembl. If you use the cellranger pipeline, you should download the gtf that comes prepackaged with it here.
2.下载 expressed repeats annotation(可选）
You might want to mask expressed repetitive elements, since those count could constitute a confounding factor in the downstream analysis. To do so you would need to download an appropriate expressed repeat annotation (for example from UCSC genome browser and make sure to select GTF as output format).
我的数据是大鼠的，所以选择rn6
image.png
所以输入关键三个要素：
1.路径：outs的上级目录 (e.g.包含 outs, outs/analys and outs/filtered_gene_bc_matrices的总目录).
2.genes.gtf is the genome annotation file provided with the cellranger pipeline.
3.repeat_msk.gtf is the repeat masker file described in the Preparation section above.
例子：
velocyto run10x -m yourpath/repeat_msk.gtf yourpath_to_outs/NRVC yourpath/cellranger_rn6/genes/genes.gtf如果出现samtools的报错：可以参考以下解决方案
https://github.com/velocyto-team/velocyto.py/issues/192#issuecomment-545349435
2.Estimating RNA velocity
This guide covers the analysis and assumes that you have produced a .loom file using the velocyto CLI (follow the guide above).
关于velocyto run10x输入的loom文件的介绍
A valid .loom file is simply an HDF5 file that contains specific groups representing the main matrix as well as row and column attributes. Because of this, .loom files can be created and read by any language that supports HDF5.
.loom files can be easily handled using the loompy package.
Get started with the analysis
At this point you are ready to start analyzing your .loom file. To get started read our analysis tutorial and have a look at the notebooks examples.
library(devtools)
install_github("velocyto-team/velocyto.R")
library(velocyto.R)未完待续
欢迎关注~
参考：
https://velocyto.org/velocyto.py/tutorial/cli.html#running-velocyto
https://doi.org/10.1038/s41586-018-0414-6
http://pklab.med.harvard.edu/velocyto/notebooks/R/SCG71.nb.html