circos个人使用教程_circos怎么加上边框

karyotype

作图的必须文件之一，这个文件可以看作是一个坐标文件，圈图一般利用基因组为坐标进行作图，自然也可以用其他数据进行作图。
假如是利用conda安装的circos,那么你的karyotype需要位于我们的perl所在目录，比如利用which perl来确定perl所在目录，比如
(perl) [fanxuezhe@192 karyotype]$ which perl
~/anaconda3/envs/perl/bin/perl
那么就在~/anaconda3/envs/perl/bin中找到data目录，该目录就是将来karyotype文件存放的目录，到该目录下的karyotype文件夹下可以看到已经预定义好的karyotype文件 
(perl) [fanxuezhe@192 karyotype]$ pwd
/home/fanxuezhe/anaconda3/envs/perl/data/karyotype
(perl) [fanxuezhe@192 karyotype]$ ls
assembly                          karyotype.chimp.pt4.txt              karyotype.human.hg16.txt  karyotype.mouse.mm9.txt    karyotype.sorghum.txt
chromosome.band.hg19.txt          karyotype.chimp.txt                  karyotype.human.hg17.txt  karyotype.mouse.txt        karyotype.yeast.txt
chromosome.band.hg38.txt          karyotype.drosophila.dm6.hires.txt   karyotype.human.hg18.txt  karyotype.oryzasativa.txt  karyotype.zeamays.txt
dm6.hires.txt                     karyotype.drosophila.dm6.lowres.txt  karyotype.human.hg19.txt  karyotype.rat.rn4.txt      parse.karyotype
dm6.lowres.txt                    karyotype.drosophila.hires.dm3.txt   karyotype.human.hg38.txt  karyotype.rat.txt          README
karyotype.arabidopsis.tair10.txt  karyotype.drosophila.lowres.dm3.txt  karyotype.human.txt       karyotype.rm.3.txt
karyotype.arabidopsis.txt         karyotype.drosophila.txt             karyotype.mouse.mm10.txt  karyotype.rm.txt
格式类似于 
(perl) [fanxuezhe@192 karyotype]$ vim karyotype.human.txt
chr - hs1 1 0 249250621 chr1
chr - hs2 2 0 243199373 chr2
chr - hs3 3 0 198022430 chr3
chr - hs4 4 0 191154276 chr4
chr - hs5 5 0 180915260 chr5
chr - hs6 6 0 171115067 chr6
chr - hs7 7 0 159138663 chr7
chr - hs8 8 0 146364022 chr8
chr - hs9 9 0 141213431 chr9
chr - hs10 10 0 135534747 chr10
chr - hs11 11 0 135006516 chr11
chr - hs12 12 0 133851895 chr12
chr - hs13 13 0 115169878 chr13
chr - hs14 14 0 107349540 chr14
chr - hs15 15 0 102531392 chr15
chr - hs16 16 0 90354753 chr16
chr - hs17 17 0 81195210 chr17
chr - hs18 18 0 78077248 chr18
chr - hs19 19 0 59128983 chr19
chr - hs20 20 0 63025520 chr20
chr - hs21 21 0 48129895 chr21
chr - hs22 22 0 51304566 chr22
chr - hsX x 0 155270560 chrx
chr - hsY y 0 59373566 chry
band hs1 p36.33 p36.33 0 2300000 gneg
band hs1 p36.32 p36.32 2300000 5400000 gpos25
band hs1 p36.31 p36.31 5400000 7200000 gneg
band hs1 p36.23 p36.23 7200000 9200000 gpos25
band hs1 p36.22 p36.22 9200000 12700000 gneg
band hs1 p36.21 p36.21 12700000 16200000 gpos50
band hs1 p36.13 p36.13 16200000 20400000 gneg
band hs1 p36.12 p36.12 20400000 23900000 gpos25
band hs1 p36.11 p36.11 23900000 28000000 gneg
band hs1 p35.3 p35.3 28000000 30200000 gpos25
band hs1 p35.2 p35.2 30200000 32400000 gneg
band hs1 p35.1 p35.1 32400000 34600000 gpos25
band hs1 p34.3 p34.3 34600000 40100000 gneg
band hs1 p34.2 p34.2 40100000 44100000 gpos25
数据内容可以粗略分为两类，第一类为chr为起始的数据，这些数据表示染色体的数据。第二类为band表示染色体上的特定区域。对两类数据进行解释（文件是以空格为分隔符进行分隔的）：
 chr - hs1 1 0 249250621 chr1
 chr表示染色体
 hs1表示物种缩写加上染色体的编号（尽量用物种缩写，有多个物种便于区分，即使有一个用缩写也合理）
 1表示染色体的编号
 0表示染色体的起始位置
 249250621表示染色体的结束位置
 chr1表示染色体的染色(chr1对应的颜色在etc/colors.ucsc.conf中有定义，格式为rgb的染色)
 例如： 
chr1  = 153,102,0
chr2  = 102,102,0
chr3  = 153,153,30
chr4  = 204,0,0
chr5  = 255,0,0
chr6  = 255,0,204
chr7  = 255,204,204
chr8  = 255,153,0
chr9  = 255,204,0
chr10 = 255,255,0
chr11 = 204,255,0
chr12 = 0,255,0
chr13 = 53,128,0
chr14 = 0,0,204
chr15 = 102,153,255
chr16 = 153,204,255
chr17 = 0,255,255
chr18 = 204,255,255
chr19 = 153,0,204
chr20 = 204,51,255
chr21 = 204,153,255
chr22 = 102,102,102
chr23 = 153,153,153
chrX  = 153,153,153
chr24 = 204,204,204
chrY  = 204,204,204
chrM  = 204,204,153
chr0  = 204,204,153
chrUn = 121,204,61
chrNA = 255,255,255
我们也可以照葫芦画瓢定义自己的颜色 
对band数据的解释
 band hs1 p36.33 p36.33 0 2300000 gneg
 band表示当前的数据是一个band，也就是染色体中的一部分
 hs1表示当前band位于chr中定义的hs1染色体上
 p36.33表示当前的染色体位置，将来可以选择是否显示在圈图上
 p36.33表示当前的染色体位置，将来可以选择是否显示在圈图上
 0表示当前band的起始位置
 23000000表示当前band的结束位置
 gneg表示颜色，同样在color.ucsc.conf中定义
 注意
 1.band之间不能发生重叠
 2.band加和必须能够覆盖完整的染色体
 比如
 chr为0-10
 band为0-5、5-10即可 
conf文件
 
conf文件是用来对图形的各个元素以及形式进行控制的文件，conf文件可以实现对圈图更细粒度的控制。示例 
# 1.1 MINIMUM CIRCOS CONFIGURATION 
# This is a 'hello world' Circos tutorial. 
# Only required configuration elements are included.
# Subsequent tutorials in this section build on this example to
# generate a representative image with common elements found in Circos
# figures in the literature.
# Chromosome name, size and color definition
karyotype = data/karyotype/karyotype.human.txt
# The <ideogram> block defines the position, size, labels and other
# properties of the segments on which data are drawn. These segments
# are usually chromosomes, but can be any integer axis.
<ideogram>
<spacing>
# Spacing between ideograms. Suffix "r" denotes a relative value. It
# is relative to circle circumference (e.g. space is 0.5% of
# circumference).
default = 0.005r
# You can increase the spacing between specific ideograms.
#<pairwise hsY;hs1>
#spacing = 20r
#</pairwise>
</spacing>
# Ideogram position, thickness and fill. 
# Radial position within the image of the ideograms. This value is
# usually relative ("r" suffix).
radius           = 0.90r
# Thickness of ideograms, which can be absolute (e.g. pixels, "p"
# suffix) or relative ("r" suffix). When relative, it is a fraction of
# image radius.
thickness        = 20p
# Ideograms can be drawn as filled, outlined, or both. When filled,
# the color will be taken from the last field in the karyotype file,
# or set by chromosomes_colors. Color names are discussed in
# http://www.circos.ca/documentation/tutorials/configuration/configuration_files
# When stroke_thickness=0p or if the parameter is missing, the ideogram is
# has no outline and the value of stroke_color is not used.
fill             = yes
stroke_color     = dgrey
stroke_thickness = 2p
</ideogram>
################################################################
# The remaining content is standard and required. It is imported from
# default files in the Circos distribution.
# These should be present in every Circos configuration file and
# overridden as required. To see the content of these files, 
# look in etc/ in the Circos distribution.
# It's best to include these files using relative paths. This way, the
# files if not found under your current directory will be drawn from
# the Circos distribution. 
# As always, centralize all your inputs as much as possible.
<image>
# Included from Circos distribution.
<<include etc/image.conf>>
</image>
# RGB/HSV/LCH color definitions, color lists, location of fonts, fill
# patterns. Included from Circos distribution.
# In older versions of Circos, colors, fonts and patterns were
# included individually. Now, this is done from a central file. Make
# sure that you're not importing these values twice.
# *** DO NOT DO THIS ***
# <colors>
# <<include etc/colors.conf>>
# <colors>
# **********************
<<include etc/colors_fonts_patterns.conf>>
# Debugging, I/O an dother system parameters
# Included from Circos distribution.
<<include etc/housekeeping.conf>>
解释：
 1.利用karyotype=data/karyotype/*****.txt的方式来指定，data目录就是我们在karyotype中提到的data目录
 2.<ideogram></ideogram>是一个xml格式的数据，利用这个区块来划分出图形的设定参数，比如这里设定了<spacing></spacing>也就是染色体之间的距离，0.005r表示利用周长的0.005倍作为间距。radius表示圈型的内径，thickness表示圈的厚度，厚度可以用p作为单位，也可用r半径作为相对单位。fill表示是否有填充颜色，stroke_color表示边框的颜色，stroke_thickness表示边框的厚度。
 3.<image></image>表示图形的设置，比如图形的输出名称，图形的格式等等，这里引出一个重要的引用方法–include，这相当于将另一个文件中的信息放在此处。，另外几个include表达的也是这种含义。 
一开始跟着tutorial作图，能够很快的看到效果，提升成就感，并且在tutorial中的conf文件基础之上去修改能够简化很多设置。 
获取tutorial压缩包
 
#获取最新的tutorial文件
wget  http://circos.ca/distribution/circos-tutorials-current.tgz
tar -zxvf circos-tutorials-current.tgz
#转移工作目录到tutorial之下
cd circos-tutorials-0.67/tutorials/
#查看所有文件
(perl) [fanxuezhe@192 tutorials]$ ls
0  1  10  2  3  4  5  6  7  8  9  circos.png  connectogram.tgz  makeimages
#这里的0，1，2，3....与官网的教程相对应
制作自己的第一个图
 
#转移工作目录到tutorial之下的1/1中
cd circos-tutorials-0.67/tutorials/1/1
circos -conf circos.conf
在执行命令的文件夹下会看到这样的图形，这个图形是人类的23个染色体的按顺序排列的形状，至此已经完成了第一张圈图
 
 更加详细的排布后续在其他博客中继续更新
				1 Circos简介
Circos是一个非常酷炫的基因组数据可视化软件，能够将基因组数据映射到环形的基因组坐标上，用相互嵌套的环道来展示基因组数据，还可以通过连接线来呈现基因组区块之间的关系。Circos的开发者是一位加拿大的生物信息学家Martin Krzywinski，值得一提的是，他还是一位专业的摄影师。于是，科学和艺术的交融，催生了Circos这款强大的可视化软件，果真不懂艺术的码农不是好的生物学家。
自从Circos诞生以后，圈图就成了生物信息可视化的明星图形，不时见诸顶尖杂志的文章。尽管有意见
官方安装指导资料 
寂寞先生___Circos教程系列
由于官方文档基于 Ubuntu和windows给出安装指导，而我平时使用的系统CentOS6.5，所以不得不进行如下操作：
perl -v
### This is perl, v5.10.1 (*) 
 which perl
### /usr/bin/perl
sudo yum install gd gd-de...
				一．安装circos：1.从官网下载circos-0.68.tar（此时文件在Downloads中）2.在Downloads中建立文件夹，名为circos，将circos-0.68.tar移至circos文件夹中并按以下命令解压缩（首先要切换目录到Downloads/circos）：$ tar xvf circos-0.68.tar
3.生成链接ln -s circos-0.68 current
				练习一下大名鼎鼎的circos，画图风头正劲，各种CNNS引用。话说是子江兄引进国内，但是从安装到调试，从语法到逻辑，真心说比较复杂，不过看到结果的话心理面就好受点。。。
这个博主不知道是算原创还是算翻译，不过由于大部分都是自己在做测试，顺便就算翻译吧。
karyotype = data/karyotype/karyotype.human.txt #核型文件，包括染色体具体信息和颜色配置，