在微生物的研究中,常需要展示不同物种在样品中的分布情况,常规的展示图表有如堆叠图、聚类热图等。在大量高分文章中,也常看到使用circos来展示关系型数据,同样地,我们也可以使用circos图对物种与样本关系可视化。
( Nature ,2011)
绘制Circos图的方法有繁有简。如果有一定的生信基础,掌握perl语言或者R语言,可以安装circos软件绘制;当然,也可以使用OmicShare Tools的circos物种关系图进行绘制,只需上传OTU丰度表即可绘制,大大节省了精力和时间。
工具链接:
https://www.omicshare.com/tools/Home/Soft/circos
数据准备
工具可以选用物种分类水平作图,因此上传的丰度表格必须带有物种注释。
OTU丰度表
如果样本中有分组,还可以选择上传分组文件,左列样本名称,右列为分组名。
分组文件
参数选择
Circos物种关系图可以选择门(phylum)/纲(class)/目(order)/科(family)/属(genus)/种(species),六个分类水平进行绘图。
最多只能选择分类水平下丰度前10的物种进行展示,超过10个,则其他归为Others;不足10个物种,则自动使用全部物种作图,也可以自定义展示的物种数(不超过10)。
Circos图结果
Circos图中左侧为物种,右侧则为样本,左侧为物种,各用不同颜色表示。连线则为样本与物种的关系,线段越粗,表明该物种在该样本中越丰富。
如下图中,深紫色代表的Pontibacter,其跟TC-3的连线最粗,CK-1连线最细,则Pontibacter这一物种在TC-3中是更丰富的。
而同一个样本TC-3中,浅蓝色占比最大,则说明TC-3这个样本中,Sphingomonas最丰富。
除了物种与样本关系,Circos图还适用于转录组、微生物、代谢物等数据分布展示,如分组+基因/功能+表达量;Go term+上下调+基因数量;分组+富集的代谢通路+通路的基因数量等。
绘制也可以使用OmicShare Tools的另一个工具:动态Circos图,欢迎大家探索绘制~
参考文献
[1]Berger MF, Lawrence MS, Demichelis F ,et al.The genomic complexity of primary human prostate cancer .Nature 470:214-220(2011).
*未经许可,不得以任何方式复制或抄袭本篇文章之部分或全部内容。版权所有,侵权必究。
基迪奥生物|专业定制测序服务
联系方式:020-39341079;service@genedenovo.com返回搜狐,查看更多
责任编辑: