這里是佳奧搪桂!多一個(gè)工具多一種方法透敌,除了R語(yǔ)言,下游分析的方法還有很多很多踢械。
這里記錄一些我使用TBtools的心得酗电。
1 、安裝
https://github.com/CJ-Chen/TBtools
2内列、手冊(cè)
https://www.yuque.com/cjchen/hirv8i
3撵术、熱圖
按住Alt水平或豎直方向移動(dòng)
4、Fasta Tools——Fasta Stats .fq信息統(tǒng)計(jì)
5话瞧、Fasta Tools——Fasta Extract (Recommended) 根據(jù)ID提取序列
每次使用需要點(diǎn)擊Initialize初始化
6嫩与、Fasta Tools——Sequence Manipulate (Rev&Comp) 對(duì)序列進(jìn)行操作
顯示互補(bǔ)、反向交排、只看ID等
7划滋、Fasta Tools——ID Simplify 簡(jiǎn)化基因ID
8、ORF Prediction——Batch Translate CDS to Protein 根據(jù)CDS預(yù)測(cè)蛋白序列
Input/Output
9埃篓、ORF Prediction——Complete/Partial ORF Predict (Single Mode) 翻譯序列古毛,單一序列
10、GFF3/GTF Manipulate——Gtf/Gff3 Sequence Extract 提取序列
Feature Tag:選取特征都许,CDS,exon嫂冻,protein等
Up Stream Bases:選取上游2000bp序列提取
Retain Only UpStream or DownStream Bases:只保留上下游序列
11胶征、BLAST GUI Wrapper——Blast Zone 本地比對(duì)
11.1 NCBI Genome下載參考基因組
比如擬南芥的:
https://www.ncbi.nlm.nih.gov/genome/4
Arabidopsis thaliana (thale cress) Reference genome: Arabidopsis thaliana (assembly TAIR10.1)
[下載這里]Download sequences in FASTA format for genome, transcript, protein
Download genome annotation in GFF, GenBank or tabular format BLAST against Arabidopsis thaliana genome, transcript, protein
下載基因組和蛋白組文件
左下config設(shè)置路徑
+添加文件
-刪除文件
11.2 本地比對(duì)
選擇基因.fa序列
Paste Input/File Input
Output需要給文件命名:D:/TBtools/Result/1
選擇左側(cè)數(shù)據(jù)庫(kù),Start
可視化 Visulize——Alignment Shower
根據(jù)ID找序列:右鍵數(shù)據(jù)庫(kù) Fetch Sequences
一鍵建樹 Build a Phylo Tree
12桨仿、BLAST GUI Wrapper——Several Sequences to a Big File [Commonly Used] 序列比對(duì)參考基因組
Outfmt:Pairwise或默認(rèn)值BlastXML
Visualize:可視化結(jié)果
13睛低、BLAST XML Visualize——Alignment View 比對(duì)XML格式文件
14、BLAST XML Visualize——Pileup Graph 檢測(cè)高通量數(shù)據(jù)覆蓋率
15、Gene Ontology——GO Enrichment 富集分析
需要:
go-basic.obo文件
注釋Background文件 一列基因ID钱雷,一列GO號(hào)
DEG差異基因列表
得到excel輸出文件
16骂铁、Common Tools——Enrichment Bar Plot GO/KEGG結(jié)果可視化
需要:
GO.Enrichment.final.xls
其他軟件得到的結(jié)果:Other
17、HeatMap llustrator——HeatMap 熱圖繪制
參數(shù)很多罩抗,多嘗試
還可添加行拉庵、列的group名
疊加熱圖:右鍵Panel Editor——Copy Panel——Paste Panel
18、HeatMap llustrator——Super HeatMap Browser——eFP Browser 熱圖疊加在解剖圖上
19套蒂、Venn and Upset Plot——Venn (Up to Six Sets) 韋恩圖
輸入.ids文件
Region Overlap Mode:查看區(qū)間之間交集
20钞支、Venn and Upset Plot——UpSet Plot (Up to Any Sets)
21、Venn and Upset Plot——UpSet Plot (Batch) 直接拖入文件做交集
22操刀、Advanced Circos
需要:
染色體長(zhǎng)度文件:Chr1 23000 可用Fasta state提取染色體長(zhǎng)度
Discrete Color Scheme Generator:輸入染色體長(zhǎng)度文件烁挟,制作染色體著色文件
LinkRegion文件:共線性分析 Chr1 59560 63890 Chr1 26217 27221 240,240,240,1.0CRLF
GenePos文件:Chr1 ARF3 59560 63890 240,240,240CRLF
GeneDensity文件:Chr1 0 10000 26 繪制成基因密度熱圖/線條/柱形圖
QTLs文件:定位
HighLight文件:標(biāo)注,調(diào)整透明度
調(diào)整LinkPos