写在前面
最近在作图,一直迷恋上TBtools了,真是感觉上手后,都不想自己写代码来分析或画图了。TBtools确实牛X,再一次发自内心的!!!!昨晚,使用TBtools来做共线分析后做Circos图。关于使用TBtools做Circos图,CJ也做了很多的教程,以及也有很的同学或公众号推了相关的教程。具体可以看,定制Circos图:使用TBtools,从数据准备到可视化,视频手把手 | 用户教你Advanced Circos 画最好看的象形图~,还有很多。这个教程也只是自己粗糙的记录自己操作过程(是真的很粗......)。
--Du
![](https://img.haomeiwen.com/i18865546/b48d0004045bdbb7.jpg)
所需文件
-
基因组fa文件
-
基因组注释文件(gtf or gff)
一、获得染色体长度信息
![](https://img.haomeiwen.com/i18865546/877ada702af2a37a.png)
获得结果
![](https://img.haomeiwen.com/i18865546/6486da39aa570b7f.png)
二、获得基因位置信息文件
![](https://img.haomeiwen.com/i18865546/47334867455d38de.png)
拖拽gff
文件到TBtools软件,点击Start
。
![](https://img.haomeiwen.com/i18865546/f0780015f63a9cc9.png)
删除不必要的信息
![](https://img.haomeiwen.com/i18865546/b65a70887d53dedb.png)
三、提取目标基因的信息
![](https://img.haomeiwen.com/i18865546/adb2dbf1c8f2d71e.png)
![](https://img.haomeiwen.com/i18865546/f505f225706b5412.png)
拖拽文件文件时,
Selected Coluumn
选择我们要match的列![](https://img.haomeiwen.com/i18865546/7f52123b65caedcc.png)
四、基因组内的同源复制信息制作
使用MCScanx,分别输入物种的fa和gtf信息。
![](https://img.haomeiwen.com/i18865546/98fcae58d182bfee.png)
结束以后,可以获得geneLinkedRegion.tab.xls
的文件,
![](https://img.haomeiwen.com/i18865546/ae02478a0eea1eaa.png)
五、选择目的基因关联信息
使用Table Row Extract or Filter
进行match。
![](https://img.haomeiwen.com/i18865546/9c6828b12952f29b.png)
下面的操作为重点,不然匹配出来的文件是空。
![](https://img.haomeiwen.com/i18865546/f48990ded550dad4.png)
这样就可以获得次作物中,基因组复制的信息位点。
将获得的信息复制到geneLinkedRegion.tab.xls
的文件中,并修改颜色。
![](https://img.haomeiwen.com/i18865546/e7eef6c59dc727fc.png)
六、出图
将以上的信息文件输入到Advanced Circos
中即可出图
![](https://img.haomeiwen.com/i18865546/37081ef4c3f61588.png)
七、调整
出图后,我们自己需要调整很多的参数,如颜色,字体大小,距离......等等,具体可以根据自己的需求来调整。
![](https://img.haomeiwen.com/i18865546/b3f5a2ee76638a4b.png)
在这里可以调整图形整体的参数,如chr、距离等。
![](https://img.haomeiwen.com/i18865546/4e38846b4353b3c0.png)
比如,调整
Feature Label
可以调整显示基因的间距。![](https://img.haomeiwen.com/i18865546/304425920111325a.png)
调整染色体的颜色
更多的参数调整,我们可以更具那些优秀的教程一步一步的操作。以及,需要我们自己结合自己实际进行操作。
![](https://img.haomeiwen.com/i18865546/cbe6fd55bd9e3b69.png)
对于第一次绘制类似图形的童鞋,我还是建议找一张自己看着好看的图形作为参考模仿着绘制。自己作图时一直这样做的,模仿着并优化着,你可以将多张图形整合到自己图形中。
ENDING !!!
往期文章:
1. 最全WGCNA教程(替换数据即可出全部结果与图形)
2. 精美图形绘制教程
3. 转录组分析教程
- 1.课程介绍
- 2.第一章 Linux基础
- 3.第一章 生信软件安装
- 4.第二章 转录组数据的下载
- 5.第三章 参考基因组序列和注释文件的下载
- 6.转录组上游分析教程 | 第四章 数据过滤
- 7.第四章 Hisat2进行数据比对
- 8.第四章 Bowtie2进行数据比对
- 9.第四章 BWA进行数据比对
- 10.第四章 Tophat2比对
- 11.第五章 无参考基因组的转录组分析
- 12.第六章 转录本定量分析
小杜的生信筆記,主要发表或收录生物信息学的教程,以及基于R的分析和可视化(包括数据分析,图形绘制等);分享感兴趣的文献和学习资料!!
网友评论