10x 空間轉(zhuǎn)錄組數(shù)據(jù)跑完 Space Ranger 之后都會生成單個樣本的 cloupe.cloupe 文件,這個文件可以用 Loupe Browser 來對單個樣本進(jìn)行可視化,Loupe Browser 操作簡單、結(jié)果直觀,非常適合不太會寫代碼的老師或同學(xué)來自己挖掘數(shù)據(jù)。但是,做過單細(xì)胞或空間轉(zhuǎn)錄組項目的同學(xué)應(yīng)該都知道,實(shí)際上做數(shù)據(jù)分析的時候是需要多個樣本合并起來分析的,如果再用單個樣本進(jìn)行可視化就不太合適了,所以我們需要把多個樣本的 cloupe.cloupe 整合成一個合并的 cloupe 文件,這樣就方便自己用 Loupe Browser 來有效的挖掘有價值的信息。
這里我們來介紹一下,怎么用 10x spaceranger aggr 來合并多個空間轉(zhuǎn)錄組測序的樣本數(shù)據(jù)。
步驟一
確保單個樣本的 spaceranger count 已經(jīng)跑完了,假如我們前面運(yùn)行 spaceranger count 時生成了三個文件:
1 $ spaceranger count --id=LV123 ……
2 wait for pipeline to finish ……
3 $ spaceranger count --id=LB456 ……
4 wait for pipeline to finish ……
5 $ spaceranger count --id=LP789 ……
6 wait for pipeline to finish ……
步驟二
準(zhǔn)備樣本信息的 csv 文件,內(nèi)容格式如下(注意逗號分隔)。
列信息說明:
library_id: 樣本 id;
molecule_h5: 運(yùn)行 spaceranger count 生成的 molecule_info.h5 文件路徑;
cloupe_file: 運(yùn)行 spacerangercount 生成的 cloupe.cloupe 文件路徑;
spatial_folder: 運(yùn)行 spaceranger count 生成的 spatial 文件夾路徑。
除了上面指定的 4 列之外,也可以加入其它分類信息的列,之后會在 Loupe Browser 可視化的時候當(dāng)成一個分組方式來展示。比如下面,增加一列樣本分組信息,不過注意這里因?yàn)閷挾鹊脑虬?spatial_folder 列省略掉了,實(shí)際上 spatial_folder 這一列還是有的。
另外需要注意,spaceranger aggr 是不進(jìn)行批次校準(zhǔn)的,所以不同的染色方法的組織是不能合并到一起的(例如免疫熒光和 H & E 染色的組織切片)。
步驟三
運(yùn)行 spaceranger aggr:
1 $ spaceranger aggr --id=AGG123 \
2 -- csv=AGG123_libraries.csv \
3 -- normalize=mapped
參數(shù)說明:
--id: 輸出文件的 id
--csv: 樣本信息 csv 文件
--normalize:depth 歸一化的方法,默認(rèn)是 mapped,也可以選 none。
結(jié)果輸出:
成功運(yùn)行后會到的下面的信息(包括主要的輸出文件)
1 2018-10-04 13:36:33[runtime] (run:local)
ID.AGG123.SPATIAL_RNA_AGGREGATOR_CS.SPATIAL_RNA_AGGREGATOR.SUMMARIZE_AGGREGATED_REPORTS.fork0.join
2 2018-10-04 13:36:36 [runtime] (join_complete)
ID.AGG123.SPATIAL_RNA_AGGREGATOR_CS.SPATIAL_RNA_AGGREGATOR.SUMMARIZE_AGGREGATED_REPORTS
3 2018-10-04 13:36:45 [runtime] VDR killed 210files, 29MB.
4
5 Outputs:
6 - Aggregation metrics summary HTML: /home/jdoe/runs/AGG123/outs/web_summary.html
7 - Aggregation metrics summary JSON: /home/jdoe/runs/AGG123/outs/summary.json
8 - Secondary analysis output CSV: /home/jdoe/runs/AGG123/outs/analysis
9 - Filtered feature-barcode matrices MEX: /home/jdoe/runs/AGG123/outs/filtered_feature_bc_matrix
10 - Filtered feature-barcode matrices HDF5: /home/jdoe/runs/AGG123/outs/filtered_feature_bc_matrix.h5
11 - Unfiltered feature-barcode matrices MEX: /home/jdoe/runs/AGG123/outs/raw_feature_bc_matrix
12 - Unfiltered feature-barcode matrices HDF5: /home/jdoe/runs/AGG123/outs/raw_feature_bc_matrix.h5
13 - Copy of the input aggregation CSV: /home/jdoe/runs/AGG123/outs/aggregation.csv
14 - Loupe Browser file: /home/jdoe/runs/AGG123/outs/cloupe.cloupe
15 - Aggregated tissue positions list: /home/jdoe/runs/AGG123/outs/aggr_tissue_positions_list.csv
16 - Spatial folder containing spatial images and scalefactors: /home/jdoe/runs/AGG123/outs/spatial
17
18 Pipestance completed successfully!
伯豪生物 空間轉(zhuǎn)錄組測序 服務(wù)優(yōu)勢
高質(zhì)量切片:切片、貼片經(jīng)驗(yàn)豐富,針對不同組織優(yōu)化了解決方案。
流程化分析:完善的分析流程,準(zhǔn)確快速解析空間轉(zhuǎn)錄組數(shù)據(jù)。
標(biāo)準(zhǔn)化內(nèi)控:豐富的實(shí)操經(jīng)驗(yàn)構(gòu)建了標(biāo)準(zhǔn)化的內(nèi)控體系。
專業(yè)化團(tuán)隊:資深的技術(shù)團(tuán)隊具有多年項目方案設(shè)計、實(shí)驗(yàn)操作、售后分析等經(jīng)驗(yàn)。
全流程服務(wù):提供組織冷凍包埋、貼片、切片、透化、建庫測序及數(shù)據(jù)分析的全套服務(wù)。
更多伯豪生物人工服務(wù):