单细胞转录组之cellranger探索——1.安装

实验室新服务器需要安装cellranger,发现最新版(3.1.0)与旧版(1.2.0)有点差别,故整理流程如下:

1. 下载安装

1.1下载地址

[https://support.10xgenomics.com/single-cell-gene-expression/software/downloads/latest](https://support.10xgenomics.com/single-cell-gene-expression/software/downloads/latest)
wget -O cellranger-3.1.0.tar.gz "http://cf.10xgenomics.com/releases/cell-exp/cellranger-3.1.0.tar.gz?Expires=1564516562&Policy=eyJTdGF0ZW1lbnQiOlt7IlJlc291cmNlIjoiaHR0cDovL2NmLjEweGdlbm9taWNzLmNvbS9yZWxlYXNlcy9jZWxsLWV4cC9jZWxscmFuZ2VyLTMuMS4wLnRhci5neiIsIkNvbmRpdGlvbiI6eyJEYXRlTGVzc1RoYW4iOnsiQVdTOkVwb2NoVGltZSI6MTU2NDUxNjU2Mn19fV19&Signature=LjEM6kHZLAa-5QpBoR8adaynXbxxMczkx9hYBprJqYBTeM1duSTqNQ-kTe-qvNWY18C-mrjX0cql9jETFVUHIgQzh3UILSdzScAH82dk7iCaDPcsjpwkajQpR9gQf5Es2laj5yvRDw0vgd1OsifTn5XwDxc557Bba58mTxz5TChZsh1Apf~2sgQkDqh4toBRJXXfn7EzX~VbxI1x3ZOnyAqS0D1h54T5K6mL1ez3LPIRfk7WWifNs0dT4B6qcBHZ~vTjwo6kU2A2azM0EzUdnJvutdO1BoIaq4ymEjdvug4r8AEhAjnM7knc9CeGlTp9Pd7Fl5rneafCPw5sYfqfSw__&Key-Pair-Id=APKAI7S6A5RYOXBWRPDA"

1.2 安装指南

[https://support.10xgenomics.com/single-cell-gene-expression/software/pipelines/latest/installation](https://support.10xgenomics.com/single-cell-gene-expression/software/pipelines/latest/installation)

1.2.1安装cellranger

cd biosoft/ #(文件目录)
tar -xzvf cellranger-3.1.0.tar.gz

1.2.2 添加环境变量便于调用

#临时变量
export PATH=/biosoft/cellranger-3.1.0:$PATH
#当前用户永久变量
vim ~/.bashrc
export PATH=“~/software/biosoft/cellranger-3.1.0:$PATH”

1.2.3基本命令

    cellranger mkfastq #拆分

    cellranger count #定量
    cellranger aggr #标准化,见仁见智
    cellranger reanalyze #分析
    cellranger mat2csv

    cellranger mkgtf #gtf索引
    cellranger mkref # 索引

    cellranger vdj

    cellranger mkvdjref

    cellranger testrun #测试
    cellranger upload #问题反馈
    cellranger sitecheck #系统检测

1.2.4 references

cellranger 提供了集成human+mouse的参考索引。
下载安装:

下载地址 wget http://cf.10xgenomics.com/supp/cell-exp/refdata-cellranger-GRCh38-and-mm10-3.1.0.tar.gz
tar -xzvf refdata-cellranger-GRCh38-3.0.0.tar.gz

自建索引:#针对3.1.0,其他物种可以参考这种方式

参考:[https://support.10xgenomics.com/single-cell-gene-expression/software/release-notes/build#grch38mm10_3.1.0](https://support.10xgenomics.com/single-cell-gene-expression/software/release-notes/build#grch38mm10_3.1.0)
#以人为例
1.获取fasta参考基因组
wget ftp://ftp.ensembl.org/pub/release-93/fasta/homo_sapiens/dna/Homo_sapiens.GRCh38.dna.primary_assembly.fa.gz
gunzip Homo_sapiens.GRCh38.dna.primary_assembly.fa.gz
2.获取gtf注释文件
wget ftp://ftp.ensembl.org/pub/release-93/gtf/homo_sapiens/Homo_sapiens.GRCh38.97.chr.gtf.gz
gunzip Homo_sapiens.GRCh38.97.chr.gtf.gz
3.构建cellranger调用gtf
cellranger mkgtf Homo_sapiens.GRCh38.97.chr.gtf Homo_sapiens.GRCh38.97.chr.filtered.gtf \
                 --attribute=gene_biotype:protein_coding \
                 --attribute=gene_biotype:lincRNA \
                 --attribute=gene_biotype:antisense \
                 --attribute=gene_biotype:IG_LV_gene \
                 --attribute=gene_biotype:IG_V_gene \
                 --attribute=gene_biotype:IG_V_pseudogene \
                 --attribute=gene_biotype:IG_D_gene \
                 --attribute=gene_biotype:IG_J_gene \
                 --attribute=gene_biotype:IG_J_pseudogene \
                 --attribute=gene_biotype:IG_C_gene \
                 --attribute=gene_biotype:IG_C_pseudogene \
                 --attribute=gene_biotype:TR_V_gene \
                 --attribute=gene_biotype:TR_V_pseudogene \
                 --attribute=gene_biotype:TR_D_gene \
                 --attribute=gene_biotype:TR_J_gene \
                 --attribute=gene_biotype:TR_J_pseudogene \
                 --attribute=gene_biotype:TR_C_gene
根据服务器配置时间长短略有不同,我这里大概一分钟。
4.构建参考索引
cellranger mkref --genome=GRCh38 \
                 --fasta=Homo_sapiens.GRCh38.dna.primary_assembly.fa \
                 --genes=Homo_sapiens.GRCh38.97.chr.filtered.gtf \
                 --ref-version=3.1.0
这里cellranger用STAR构建ref,时长约半个小时。
也可以利用 cellranger mkref 将两个物种ref合并。3.10版本建议直接采用合并方式,以前均有单独ref方式。
具体区别还得仔细看文档。
此外,旧服务器1.2版本还有ERCC reference,3.1中已被取消。10x不建议加ERCC,Smart-seq2可以加。

1.2.5 测试

cellranger建议运行以下命令来看系统能否支持运算:
cellranger sitecheck > sitecheck.txt
cellranger upload your@email.edu sitecheck.txt
生成一个sitecheck.txt文件,记录一些系统配置的东西,如果配置不满足,官网会发信给你。

1.2.6预运行

cellranger提供预运行检查流程能否跑的通。
cellranger testrun --id=tiny
会在当前目录下生成tidy文件,cellranger3.1下生成cellranger-tiny-fastq与cellranger-tiny-ref文件。
tiny下有tiny.mri.tgz,可以作为问题诊断指南。
如果预运行失败:
cellranger upload your@email.edu tiny/tiny.mri.tgz
发邮件给官网即可。

明天更新四个pipeline使用方法。接下来计划为:
1.软件探索
2.PBMC数据集使用
3.seurat在10x上的应用
4.long ranger探索与使用

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 212,686评论 6 492
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 90,668评论 3 385
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 158,160评论 0 348
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 56,736评论 1 284
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 65,847评论 6 386
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 50,043评论 1 291
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 39,129评论 3 410
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 37,872评论 0 268
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 44,318评论 1 303
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 36,645评论 2 327
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 38,777评论 1 341
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 34,470评论 4 333
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 40,126评论 3 317
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 30,861评论 0 21
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,095评论 1 267
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 46,589评论 2 362
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 43,687评论 2 351

推荐阅读更多精彩内容