BSTMatrix v_2.2流程說明

1.軟件依賴

1）python： 版本3.8及以上，安裝cv2模塊（4.0以上版本），matplotlib，seaborn，pandas
2）Star：版本2.6.1d及以上版本
3）perl：threads及threads::shared模塊
4）Rscript：需要以下R包，seurat、dplyr、tibble、ggplot2、broom、purrr、cowplot、cluster、ggpubr、plotly、htmlwidgets、kableextra、htmltools、shiny、knitr、rmarkdown、optparse。
依賴的軟件需要使用export添加到環(huán)境變量中，以實(shí)現(xiàn)流程的調(diào)用。

2.輸入數(shù)據(jù)準(zhǔn)備

1）測(cè)序數(shù)據(jù)：雙端測(cè)序fastq數(shù)據(jù)。
2）參考基因組數(shù)據(jù)：基因組序列文件，gtf文件（第3列需要包含exon），gff文件（可選）（第3列需要包含gene、exon）。
3）features.tsv文件：可使用gtf文件生成，參考命令：perl ./tools/features_generate.pl -i xxx.gtf -o features.tsv
4）STAR基因組索引文件：可使用基因組序列文件和gtf文件生成，參考命令：STAR –runThreadN 8 –runMode genomeGenerate –genomeDir star/ –genomeFastaFiles genome.fa –sjdbGTFfile gene.gtf
5）熒光解碼文件及HE圖片文件。

3.配置文件編寫

配置文件：
## fq測(cè)序數(shù)據(jù)文件路徑，支持.gz格式
FQ1 /path/to/read_1.fq.gz
FQ2 /path/to/read_2.fq.gz
## Flu info file 熒光解碼文件路徑
FLU /path/to/flu_info.txt
## 參考基因組STAR建庫(kù)目錄及gff/gtf文件路徑
INDEX /path/to/STAR/index/dir/
GFF /path/to/ref/gene/gff3/file #（也可使用gtf文件）
## 參考基因組features.tsv文件路徑
FEATURE /path/to/features.tsv
## HE染色圖片路徑
HE /path/to/HE.tif
## 輸出目錄及輸出文件前綴
OUTDIR /path/to/result/dir/
PREFIX outfile-prefix
### 程序參數(shù)
## fastq2BcUmi
BCType V2 #barcode 版本類型(一般為V2版本)
BCThreads 8 #線程數(shù)
## Umi2Gene
Sjdboverhang 100 #STAR建庫(kù)時(shí)使用的-sjdboverhang參數(shù)值，默認(rèn)100
STARThreads 8 #STAR比對(duì)線程數(shù)

ENV python和Rscript的路徑，如不提供則使用系統(tǒng)環(huán)境中的版本(不提供請(qǐng)注釋掉以下參數(shù))
PYTHON /path/to/python/dir/
Rscript /path/to/Rscript/dir/

4.流程運(yùn)行

1）流程說明：
流程分為6個(gè)步驟，如下所示：
A）步驟1：運(yùn)行fastq2BcUmi，識(shí)別fastq數(shù)據(jù)中的barcode、umi。
B）步驟2：運(yùn)行LinkBcChip，識(shí)別熒光數(shù)據(jù)的barcode信息并對(duì)應(yīng)到芯片上位置。
C）步驟3：運(yùn)行Umi2Gene，將reads與參考基因組比對(duì)，得到每個(gè)UMI對(duì)應(yīng)的基因信息。
D）步驟4：運(yùn)行MatrixMake，獲得基因表達(dá)矩陣。
E）步驟5：運(yùn)行AllheStat，處理HE圖片。
F）步驟6：運(yùn)行cluster.R，進(jìn)行聚類分析。
G）步驟7：運(yùn)行WebReport，得到網(wǎng)頁(yè)版報(bào)告。
2）流程參數(shù)：
-c config.txt 數(shù)據(jù)配置文件
-s 步驟選擇，0為運(yùn)行1-7所有步驟，也可選擇個(gè)別步驟單獨(dú)運(yùn)行，多個(gè)步驟中間使用“，”分割。
3）參考命令：
./BSTMatrix -c config.txt -s 0
./BSTMatrix -c config.txt -s 1,2,3,4,5,6,7
./BSTMatrix -c config.txt -s 1,3

5.結(jié)果文件說明

1）目錄結(jié)構(gòu)及結(jié)果說明：
outdir/
├── 01.fastq2BcUmi 步驟1運(yùn)行結(jié)果目錄
│?? ├── prefix.bc_dist 不同barcode檢測(cè)統(tǒng)計(jì)
│?? ├── prefix.bc_stat 不同barcode檢測(cè)統(tǒng)計(jì)
│?? ├── prefix.bc_umi_read.tsv barcode類型、對(duì)應(yīng)的umi及reads數(shù)統(tǒng)計(jì)文件
│?? ├── prefix.bc_umi_read.tsv.id barcode類型、對(duì)應(yīng)的umi及reads id文件
│?? ├── prefix.filter 沒有完整識(shí)別出來barcode的reads信息
│?? ├── prefix.full_stat barcode類型對(duì)應(yīng)的reads數(shù)、umi數(shù)
│?? ├── prefix.id_map id編號(hào)對(duì)應(yīng)關(guān)系文件
│?? ├── prefix.qual.stat reads統(tǒng)計(jì)文件
│?? ├── prefix.select_id 完整識(shí)別出barcode和UMI的reads id文件
│?? ├── prefix.stat barcode檢測(cè)統(tǒng)計(jì)
│?? ├── prefix.umi reads對(duì)應(yīng)的barcode類型及umi
│?? └── prefix.umi_cor.info umi校正信息
├── 02.LinkBcChip 步驟2運(yùn)行結(jié)果目錄
│?? ├── prefix.barcode_pos.tsv barcode類型對(duì)應(yīng)的芯片位置文件
│?? ├── prefix.barcode.tsv 芯片對(duì)應(yīng)的barcode類型文件
│?? ├── prefix.used.barcode
│?? └── prefix.pos reads對(duì)應(yīng)的芯片位置及barcode類型文件
├── 03.Umi2Gene 步驟3運(yùn)行結(jié)果目錄
│?? ├── prefixAligned.sortedByCoord.out.bam STAR軟件比對(duì)結(jié)果
│?? ├── prefix.final.transcript
│?? ├── prefixLog.final.out STAR軟件比對(duì)結(jié)果信息文件
│?? ├── prefix.cut90.fq 剪切成90bp長(zhǎng)度的reads文件
│?? ├── prefixLog.out
│?? ├── prefixLog.progress.out
│?? ├── prefix.map2exon
│?? ├── prefix.map2gene reads比對(duì)到的基因信息
│?? ├── prefix.map2intron
│?? ├── prefix.map2transcript
│?? ├── prefix.mapadjust.exon
│?? ├── prefix.mapadjust.intron
│?? ├── prefix.mapadjust.transcript
│?? ├── prefixSJ.out.tab
│?? ├── prefix_STARgenome
│?? ├── prefix_STARpass1
│?? ├── prefix_STARtmp
│?? ├── prefix.stat
│?? ├── prefix.total.stat reads比對(duì)結(jié)果統(tǒng)計(jì)文件
│?? └── prefix.umi_gene.tsv barcode對(duì)應(yīng)的umi及基因文件
├── 04.MatrixMake 步驟4運(yùn)行結(jié)果目錄
│?? ├── prefix.matrix.tsv 基因表達(dá)矩陣文件
│?? ├── prefix.matrix.tsv.filt
│?? ├── prefix.select.bc_umi_read.tsv
│?? ├── prefix.select.umi_gene.tsv
│?? ├── prefix.select.umi_gene.tsv.filter
│?? └── prefix.sequencing_saturation.stat
├── 05.AllheStat 步驟5運(yùn)行結(jié)果目錄
│?? ├── allhe
│?? ├── all_level_stat.txt 不同水平的spots統(tǒng)計(jì)
│?? ├── BSTViewer_project BSTViewer軟件輸入數(shù)據(jù)
│?? ├── heAuto_level_matrix 不同分辨率組織內(nèi)spot的矩陣目錄
│?? ├── level_matrix 不同分辨率的矩陣目錄
│?? ├── stat.txt
│?? └── umi_plot umi count圖片統(tǒng)計(jì)圖目錄
├── 06.Cluster 步驟6 聚類結(jié)果目錄
│?? ├── L13 L13水平聚類結(jié)果目錄
│?? ├── L3 L3水平聚類結(jié)果目錄
│?? ├── L4 L4水平聚類結(jié)果目錄
│?? ├── L5 L5水平聚類結(jié)果目錄
│?? ├── L6 L6水平聚類結(jié)果目錄
│?? └── L7 L7水平聚類結(jié)果目錄
├── 07.WebReport 步驟7運(yùn)行結(jié)果目錄
│?? ├── prefix.filelist
│?? ├── report.html 單文件版網(wǎng)頁(yè)版報(bào)告文件
│?? ├── index.html 網(wǎng)頁(yè)版報(bào)告html文件
│?? └── src 網(wǎng)頁(yè)版報(bào)告src目錄
└── prefix 收集的基因表達(dá)矩陣等文件目錄
├── barcode_pos.tsv barcode類型對(duì)應(yīng)的芯片位置文件
├── barcode.tsv 芯片對(duì)應(yīng)的barcode類型文件
├── bc_umi_read.tsv.gz barcode類型、對(duì)應(yīng)的umi及reads數(shù)統(tǒng)計(jì)文件
├── umi_gene.tsv.gz barcode對(duì)應(yīng)的umi及基因文件
├── features.tsv features.tsv文件
└── matrix.tsv 基因表達(dá)矩陣文件