如何使用cell ranger拆分10X單細胞轉錄組原始數(shù)據(jù)，相信很多沒有經驗的人對此束手無策，為此本文總結了問題出現(xiàn)的原因和解決方法，通過這篇文章希望你能解決這個問題。

創(chuàng)新互聯(lián)建站是一家專業(yè)從事成都做網站、成都網站建設的網絡公司。作為專業(yè)網絡公司,創(chuàng)新互聯(lián)建站依托的技術實力、以及多年的網站運營經驗,為您提供專業(yè)的成都網站建設、全網營銷推廣及網站設計開發(fā)服務！

cell ranger是10X genomics公司提供的，專門用于分析10X 單細胞轉錄組數(shù)據(jù)的pipeline, 包含了原始數(shù)據(jù)拆分，表達定量，聚類分析等多個功能，本文主要介紹如何使用該軟件來拆分原始數(shù)據(jù)。

直接從官網下載最新版的軟件即可，網址如下

https://support.10xgenomics.com/single-cell-gene-expression/software/downloads/latest

該軟件由多個子命令構成，通過mkfastq命令拆分數(shù)據(jù)，流程示意如下

有以下兩種使用方式

cellranger mkfastq \
--id test \
--run run_directory \
--csv simple.csv

cellranger mkfastq \
--id test \
--run  run_directory \
--samplesheet samplesheet.csv

id參數(shù)指定輸出目錄的名字，run參數(shù)指定下機的原始bcl文件所在的目錄，該命令其實是對illumina提供的拆分數(shù)據(jù)的bcl2fastq命令的一個封裝，需要樣本名稱，index等信息，支持兩種格式，一種就是illlumina常規(guī)的samplesheet.csv文件，格式如下

另外一種是10X  genomics定制的一種簡化版的csv格式，內容如下

Lane,Sample,Index
1,test_sample,SI-GA-A3

只有3列，第一列指定lane ID, 第二列指定樣本名稱，第三列指定index的名稱，10X  genomics的每個index代表4條具體的oligo序列，示意如下

在根據(jù)index確定樣本時，允許1到2個堿基的錯配。在實際拆分數(shù)據(jù)時，更加推薦使用三列的CSV文件，因為samplesheet文件中需要根據(jù)不同版本的試劑盒修改對應的Reads信息。

V2試劑盒產生的文庫結構如下所示

V3試劑盒產生的文庫結構如下所示

和V2相比，V3試劑盒中所用的UMI和PolyT的長度都發(fā)生了變化，從而導致測序得到的R1和R2端的序列長度也不一致，V2試劑盒的R1端長度為26bp, 包含16bp的barcode和10bp的UMI序列，V3試劑盒的R1端長度為28bp, 包含16bp的barcode和12bp的UMI序列；V2試劑盒的R2端為98bp, V3試劑盒的R2端為91bp。

如果使用samplesheet文件，需要調整[Reads]中的序列長度，而使用簡化版的csv文件，cell ranger可以識別所用試劑盒版本，然后自動化的調整reads長度。
拆分好之后的目錄結構如下所示

├── fastq_path
│   ├── H35KCBCXY
│   │   └── test_sample
│   │       ├── test_sample_S1_L001_I1_001.fastq.gz
│   │       ├── test_sample_S1_L001_R1_001.fastq.gz
│   │       └── test_sample_S1_L001_R2_001.fastq.gz

對于每個樣本，除了常見的R1和R2端序列，還多出來一個I1序列文件，該文件中保存的是index序列，示意如下

@D00547:905:H35KCBCXY:1:1101:19188:87078 1:N:0:AGATCGGG
AGATCGGG
+
.<<....<

后續(xù)的子命令也是通過這種特定的目錄結構來進行分析，如果你有從其他地方下載的原始數(shù)據(jù)，也可以整理成這種目錄結構，方便后續(xù)使用cell ranger進行分析。

看完上述內容，你們掌握如何使用cell ranger拆分10X單細胞轉錄組原始數(shù)據(jù)的方法了嗎？如果還想學到更多技能或想了解更多相關內容，歡迎關注創(chuàng)新互聯(lián)行業(yè)資訊頻道，感謝各位的閱讀！

當前名稱：如何使用cellranger拆分10X單細胞轉錄組原始數(shù)據(jù)
文章位置：http://chinadenli.net/article36/jgigpg.html

成都網站建設公司_創(chuàng)新互聯(lián)，為您提供商城網站、做網站、手機網站建設、靜態(tài)網站、標簽優(yōu)化、自適應網站

廣告

聲明：本網站發(fā)布的內容（圖片、視頻和文字）以用戶投稿、用戶轉載內容為主，如果涉及侵權請盡快告知，我們將會在第一時間刪除。文章觀點不代表本網站立場，如需處理請聯(lián)系客服。電話：028-86922220；郵箱：631063699@qq.com。內容未經允許不得轉載，或轉載時需注明來源： 創(chuàng)新互聯(lián)