CIRCOS增加热图、点图、线图和区块属性

评论6,635

CIRCOS图在有了染色体信息界定绘图区域后，就可以向里面添加离散数据如标记特定的区域或连续数据如展示修饰的丰度等。

经过前面部分对CIRCOS基本安装，最简单绘图和调整各部分属性的熟悉后，这些基本配置都可以写入单独的文件，供多次使用，就像下面主配置文件中的include所示。

而我们每次绘图主要修改的部分就在主配置文件里面，更换下对应圈的文件名、半径和绘制属性即可。

main circos configure 

(test.circos.conf)

  1. <<include etc/colors_fonts_patterns.conf>>
  2.  
  3. # # 必须的部分，控制染色体信息显示
  4. # # 之前两篇教程为了方便理解ideogram和ticks都是直接写在主配置文件的
  5. # # 但通常实际使用时会拆分成不同文件，方便管理，所以这里采用了引用方式
  6. # # 具体解释都见前面的文章
  7. <<include ideogram.conf>>
  8. <<include ticks.conf>>
  9.  
  10. # 定义输出
  11. <image>
  12. <<include etc/image.conf>>
  13. file*=test.circos.png
  14. dir*=./
  15. </image>
  16.  
  17. # karyotype定义染色体的名字、ID、起始位置信息，是绘制图的根本
  18. # 具体解释看前面2篇文章
  19. karyotype = test.chromsomes.circos_input.txt
  20.  
  21. # `chromosome_units`定义染色体一个单位的大小，缩写为`u`。若`chromosome_units=1000000`, 则`10u=10000000`。
  22. # 后面会用到这个单位，尤其是ticks中
  23. chromosomes_units =1000000
  24. chromosomes_display_default = yes
  25.  
  26. # 大标签highlights，复数定义里面有多个highlight
  27. <highlights>
  28.  
  29. # 高亮的区域放在 <highlight>标签中，配置简单，给定文件名和内外半径就好
  30. # 文件格式见下面解释
  31. <highlight>
  32. file=HL1.bed.circos_input.txt
  33. r0=0.914285714286r
  34. r1=0.991428571429r
  35. </highlight>
  36.  
  37. # 高亮的区域放在 <highlight>标签中，配置简单，给定文件名和内外半径就好
  38. # 文件格式见下面解释
  39. <highlight>
  40. file=HL2.bed.circos_input.txt
  41. r0=0.828571428571r
  42. r1=0.905714285714r
  43. </highlight>
  44.  
  45. # 高亮的区域放在 <highlight>标签中，配置简单，给定文件名和内外半径就好
  46. # 文件格式见下面解释
  47. <highlight>
  48. file=HL3.bed.circos_input.txt
  49. r0=0.742857142857r
  50. r1=0.82r
  51. </highlight>
  52.  
  53. </highlights>
  54.  
  55. # 定义绘制的线图、点图、热图、直方图等，复数plots
  56. <plots>
  57.  
  58. # 全局属性定义
  59. color = spectral-7-div-rev
  60. stroke_thickness = 1
  61. stroke_color = black
  62. scale_log_base = 0.5
  63.  
  64. # 其中一个子图，给定文件名，内外半径，和绘制类型
  65. # 文件格式见下面解释
  66. <plot>
  67. file=Heatmap1.bed.circos_input.txt
  68. type=line
  69. r0=0.657142857143r
  70. r1=0.734285714286r
  71. </plot>
  72.  
  73. # 其中一个子图，给定文件名，内外半径，和绘制类型
  74. # 文件格式见下面解释
  75. <plot>
  76. file=Heatmap2.bed.circos_input.txt
  77. type=scatter
  78. r0=0.571428571429r
  79. r1=0.648571428571r
  80. </plot>
  81.  
  82. # 其中一个子图，给定文件名，内外半径，和绘制类型
  83. # 文件格式见下面解释
  84. <plot>
  85. file=Heatmap3.bed.circos_input.txt
  86. type=histogram
  87. r0=0.485714285714r
  88. r1=0.562857142857r
  89. </plot>
  90.  
  91. # 其中一个子图，给定文件名，内外半径，和绘制类型
  92. # 文件格式见下面解释
  93. <plot>
  94. file=Heatmap4.bed.circos_input.txt
  95. type=heatmap
  96. r0=0.4r
  97. r1=0.477142857143r
  98. </plot>
  99.  
  100. </plots>
  101.  
  102. <<include etc/housekeeping.conf>>
  103. max_points_per_track*  =  2500000
  104. data_out_of_range* = trim

下面看下文件格式的要求

test.chromsomes.circos_input.txt

  1. ##chr - ID LABEL START END COLOR
  2. chr - chr1 hs1 0 248956422 chr1
  3. chr - chr2 hs2 0 242193529 chr2
  4. chr - chr3 hs3 0 198295559 chr3
  5. chr - chr4 hs4 0 190214555 chr4
  6. chr - chr5 hs5 0 181538259 chr5

高亮文件的格式一致，选其中一个做例子 HL1.bed.circos_input.txt

  1. #ParentID    START    END    ATTRIBUTE
  2. #ParentID对应于karyotype文件的ID（第3列)
  3. #ATTRIBUTE列颜色的获取看前面两篇的介绍以获知有多少可用颜色
  4. chr1 0 6422000 fill_color=vlyellow
  5. chr1 12447822 24895644 fill_color=vlpred
  6. chr1 136926042 149373864 fill_color=vlppurple
  7. chr1 236508618 248956422 fill_color=vlporange
  8. chr2 0 2421000 fill_color=vlpred
  9. chr2 36329031 48438708 fill_color=vlpgreen
  10. chr2 193754832 205864509 fill_color=vlpred
  11. chr2 133206447 145316124 fill_color=vlyellow
  12. chr3 59488668 69403446 fill_color=vlpred
  13. chr3 128892114 138806892 fill_color=vlpblue
  14. chr3 1800000 198295559 fill_color=vlpred
  15. chr4 0 2145550 fill_color=vlpgreen
  16. chr4 95107280 104618008 fill_color=vlpblue
  17. chr4 133150192 142660920 fill_color=vlppurple
  18. chr4 190000000 190214555 fill_color=vlpgreen
  19. chr5 0 1538259 fill_color=vlpblue
  20. chr5 27230739 36307652 fill_color=vvlpgreen
  21. chr5 72615304 81692217 fill_color=vvlyellow
  22. chr5 127076782 136153695 fill_color=vvlporange
  23. chr5 15382590 15392590 fill_color=vlpblue

数值文件格式也一致，只是给定不同的type绘制不同属性的图，如线图、热图等。

Heatmap1.bed.circos_input.txt

  1. #ParentID    START    END    Value
  2. #ParentID对应于karyotype文件的ID (第3列)
  3. #Value列为我们想展示的值
  4. chr1 0 12447822 1
  5. chr1 12447822 24895644 2
  6. chr1 211612974 224060796 18
  7. chr1 224060796 236508618 19
  8. chr1 236508618 248956422 20
  9. chr2 0 12109677 1
  10. chr2 12109677 24219354 2
  11. chr2 24219354 36329031 3
  12. chr2 36329031 48438708 4
  13. chr2 48438708 60548385 5
  14. chr2 205864509 217974186 18
  15. chr2 217974186 230083863 19
  16. chr2 230083863 242193529 20
  17. chr3 0 9914778 1
  18. chr3 9914778 19829556 2
  19. chr3 19829556 29744334 3
  20. chr3 29744334 39659112 4
  21. chr3 39659112 49573890 5
  22. chr4 0 9510728 1
  23. chr4 9510728 19021456 2
  24. chr4 19021456 28532184 3
  25. chr4 28532184 38042912 4
  26. chr4 133150192 142660920 15
  27. chr4 171193104 180703832 19
  28. chr4 180703832 190214555 20
  29. chr5 0 9076913 1
  30. chr5 9076913 18153826 2
  31. chr5 18153826 27230739 3
  32. chr5 154307521 163384434 18
  33. chr5 163384434 172461347 19
  34. chr5 172461347 181538259 20

所有文件都准备好之后，运行circos -conf test.circos.conf就可以获得test.circos.svg和test.circos.png两幅图了。

CIRCOS增加热图、点图、线图和区块属性

如果有属性不理解或还需修改属性，请参考前两篇及前两篇提供的属性修改文档方便查找。

发表评论

匿名网友

拖动滑块以完成验证
加载失败