10x 空間轉(zhuǎn)錄組數(shù)據(jù)跑完 Space Ranger 之后都會(huì)生成單個(gè)樣本的 cloupe.cloupe 文件，這個(gè)文件可以用 Loupe Browser 來(lái)對(duì)單個(gè)樣本進(jìn)行可視化，Loupe Browser 操作簡(jiǎn)單、結(jié)果直觀，非常適合不太會(huì)寫代碼的老師或同學(xué)來(lái)自己挖掘數(shù)據(jù)。但是，做過(guò)單細(xì)胞或空間轉(zhuǎn)錄組項(xiàng)目的同學(xué)應(yīng)該都知道，實(shí)際上做數(shù)據(jù)分析的時(shí)候是需要多個(gè)樣本合并起來(lái)分析的，如果再用單個(gè)樣本進(jìn)行可視化就不太合適了，所以我們需要把多個(gè)樣本的 cloupe.cloupe 整合成一個(gè)合并的 cloupe 文件，這樣就方便自己用 Loupe Browser 來(lái)有效的挖掘有價(jià)值的信息。
這里我們來(lái)介紹一下，怎么用 10x spaceranger aggr 來(lái)合并多個(gè)空間轉(zhuǎn)錄組測(cè)序的樣本數(shù)據(jù)。
步驟一
確保單個(gè)樣本的 spaceranger count 已經(jīng)跑完了，假如我們前面運(yùn)行 spaceranger count 時(shí)生成了三個(gè)文件：
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)容格式如下（注意逗號(hào)分隔）。
列信息說(shuō)明：
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 列之外，也可以加入其它分類信息的列，之后會(huì)在 Loupe Browser 可視化的時(shí)候當(dāng)成一個(gè)分組方式來(lái)展示。比如下面，增加一列樣本分組信息，不過(guò)注意這里因?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ù)說(shuō)明：
--id： 輸出文件的 id
--csv： 樣本信息 csv 文件
--normalize：depth 歸一化的方法，默認(rèn)是 mapped，也可以選 none。
結(jié)果輸出：
成功運(yùn)行后會(huì)到的下面的信息（包括主要的輸出文件）
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)錄組測(cè)序 服務(wù)優(yōu)勢(shì)
高質(zhì)量切片：切片、貼片經(jīng)驗(yàn)豐富，針對(duì)不同組織優(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)隊(duì)：資深的技術(shù)團(tuán)隊(duì)具有多年項(xiàng)目方案設(shè)計(jì)、實(shí)驗(yàn)操作、售后分析等經(jīng)驗(yàn)。
全流程服務(wù)：提供組織冷凍包埋、貼片、切片、透化、建庫(kù)測(cè)序及數(shù)據(jù)分析的全套服務(wù)。
更多伯豪生物人工服務(wù)：
