美文网首页单细胞测序
单细胞转录组之cellranger探索——1.安装

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

作者: Dengfeng_Gao | 来源:发表于2019-07-31 17:43 被阅读0次

    实验室新服务器需要安装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探索与使用

    相关文章

      网友评论

        本文标题:单细胞转录组之cellranger探索——1.安装

        本文链接:https://www.haomeiwen.com/subject/qtlzrctx.html