單細胞轉(zhuǎn)錄組之cellranger探索——1.安裝

實驗室新服務(wù)器需要安裝cellranger酣衷,發(fā)現(xiàn)最新版(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 添加環(huán)境變量便于調(diào)用

#臨時變量
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 #系統(tǒng)檢測

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.構(gòu)建cellranger調(diào)用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
根據(jù)服務(wù)器配置時間長短略有不同狂打,我這里大概一分鐘。
4.構(gòu)建參考索引
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構(gòu)建ref混弥,時長約半個小時趴乡。
也可以利用 cellranger mkref 將兩個物種ref合并。3.10版本建議直接采用合并方式蝗拿,以前均有單獨ref方式晾捏。
具體區(qū)別還得仔細看文檔。
此外哀托,舊服務(wù)器1.2版本還有ERCC reference惦辛,3.1中已被取消。10x不建議加ERCC仓手,Smart-seq2可以加胖齐。

1.2.5 測試

cellranger建議運行以下命令來看系統(tǒng)能否支持運算:
cellranger sitecheck > sitecheck.txt
cellranger upload your@email.edu sitecheck.txt
生成一個sitecheck.txt文件,記錄一些系統(tǒng)配置的東西嗽冒,如果配置不滿足呀伙,官網(wǎng)會發(fā)信給你。

1.2.6預(yù)運行

cellranger提供預(yù)運行檢查流程能否跑的通添坊。
cellranger testrun --id=tiny
會在當前目錄下生成tidy文件剿另,cellranger3.1下生成cellranger-tiny-fastq與cellranger-tiny-ref文件。
tiny下有tiny.mri.tgz,可以作為問題診斷指南雨女。
如果預(yù)運行失敗:
cellranger upload your@email.edu tiny/tiny.mri.tgz
發(fā)郵件給官網(wǎng)即可谚攒。

明天更新四個pipeline使用方法。接下來計劃為:
1.軟件探索
2.PBMC數(shù)據(jù)集使用
3.seurat在10x上的應(yīng)用
4.long ranger探索與使用

?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末氛堕,一起剝皮案震驚了整個濱河市馏臭,隨后出現(xiàn)的幾起案子,更是在濱河造成了極大的恐慌讼稚,老刑警劉巖位喂,帶你破解...
    沈念sama閱讀 206,013評論 6 481
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件,死亡現(xiàn)場離奇詭異乱灵,居然都是意外死亡塑崖,警方通過查閱死者的電腦和手機,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 88,205評論 2 382
  • 文/潘曉璐 我一進店門痛倚,熙熙樓的掌柜王于貴愁眉苦臉地迎上來规婆,“玉大人,你說我怎么就攤上這事蝉稳∈阊粒” “怎么了?”我有些...
    開封第一講書人閱讀 152,370評論 0 342
  • 文/不壞的土叔 我叫張陵耘戚,是天一觀的道長嗡髓。 經(jīng)常有香客問我,道長收津,這世上最難降的妖魔是什么饿这? 我笑而不...
    開封第一講書人閱讀 55,168評論 1 278
  • 正文 為了忘掉前任,我火速辦了婚禮撞秋,結(jié)果婚禮上长捧,老公的妹妹穿的比我還像新娘。我一直安慰自己吻贿,他們只是感情好串结,可當我...
    茶點故事閱讀 64,153評論 5 371
  • 文/花漫 我一把揭開白布。 她就那樣靜靜地躺著舅列,像睡著了一般肌割。 火紅的嫁衣襯著肌膚如雪。 梳的紋絲不亂的頭發(fā)上帐要,一...
    開封第一講書人閱讀 48,954評論 1 283
  • 那天把敞,我揣著相機與錄音,去河邊找鬼宠叼。 笑死先巴,一個胖子當著我的面吹牛其爵,可吹牛的內(nèi)容都是我干的冒冬。 我是一名探鬼主播伸蚯,決...
    沈念sama閱讀 38,271評論 3 399
  • 文/蒼蘭香墨 我猛地睜開眼,長吁一口氣:“原來是場噩夢啊……” “哼简烤!你這毒婦竟也來了剂邮?” 一聲冷哼從身側(cè)響起,我...
    開封第一講書人閱讀 36,916評論 0 259
  • 序言:老撾萬榮一對情侶失蹤横侦,失蹤者是張志新(化名)和其女友劉穎挥萌,沒想到半個月后,有當?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體枉侧,經(jīng)...
    沈念sama閱讀 43,382評論 1 300
  • 正文 獨居荒郊野嶺守林人離奇死亡引瀑,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點故事閱讀 35,877評論 2 323
  • 正文 我和宋清朗相戀三年,在試婚紗的時候發(fā)現(xiàn)自己被綠了榨馁。 大學(xué)時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片憨栽。...
    茶點故事閱讀 37,989評論 1 333
  • 序言:一個原本活蹦亂跳的男人離奇死亡,死狀恐怖翼虫,靈堂內(nèi)的尸體忽然破棺而出屑柔,到底是詐尸還是另有隱情,我是刑警寧澤珍剑,帶...
    沈念sama閱讀 33,624評論 4 322
  • 正文 年R本政府宣布掸宛,位于F島的核電站,受9級特大地震影響招拙,放射性物質(zhì)發(fā)生泄漏唧瘾。R本人自食惡果不足惜,卻給世界環(huán)境...
    茶點故事閱讀 39,209評論 3 307
  • 文/蒙蒙 一别凤、第九天 我趴在偏房一處隱蔽的房頂上張望劈愚。 院中可真熱鬧,春花似錦闻妓、人聲如沸菌羽。這莊子的主人今日做“春日...
    開封第一講書人閱讀 30,199評論 0 19
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽注祖。三九已至,卻和暖如春均唉,著一層夾襖步出監(jiān)牢的瞬間是晨,已是汗流浹背。 一陣腳步聲響...
    開封第一講書人閱讀 31,418評論 1 260
  • 我被黑心中介騙來泰國打工舔箭, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留罩缴,地道東北人蚊逢。 一個月前我還...
    沈念sama閱讀 45,401評論 2 352
  • 正文 我出身青樓,卻偏偏與公主長得像箫章,于是被迫代替她去往敵國和親烙荷。 傳聞我的和親對象是個殘疾皇子,可洞房花燭夜當晚...
    茶點故事閱讀 42,700評論 2 345

推薦閱讀更多精彩內(nèi)容