基因结构图从本质上就可以看成方块，直线，箭头的组合
解析GFF文件
bioconductor上有一个GenomicFeatures包，里面有一个makeTxDbFromGFF()函数可以解析GFF文件并构建TxDb对象，如何操作该对象见【Bioconductor系列】如何用Bioconductor对基因组注释.
以拟南芥TAIR10上的GFF为例(文件可以从TAIR10上下载),
library(GenomicFeatures)
GFF_file <- "C:/Users/DELL/Desktop/TAIR10_GFF3_genes.gff"
txdb <- makeTxDbFromGFF(GFF_file)
genes_df <- as.data.frame(genes(txdb))
exons_df <- as.data.frame(exons(txdb))
genes()和exons()函数分别用于获取所有基因，所有外显子的GRanges对象，之后用as.data.frame转成ggplot2使用的数据框格式
提取目标区间的基因结构信息
ggplot2作图
第一步，调整主题。如下的代码就是将画布清空
theme_syntenty <- theme_classic() +
theme(axis.line = element_blank(),
axis.text = element_blank(),
axis.ticks = element_blank(),
axis.title = element_blank())
theme)set()
第二步，画外显子的结构
外显子就是一个一个方块，我们可以用有一定厚度的线段来表示，geom_segment()
p1 <- ggplot(exons_in_gene_df,
aes(x=start, xend=end,y=0.5,yend=0.5)) +
geom_segment(size=4)
p1
外显子1
也可以用geom_rect()画框框
p1 <- ggplot(exons_in_gene_df) +
geom_rect(aes(xmin=start, xmax=end,ymin=-0.1,ymax=0.1),
colour="black", fill="white") +
ylim(c(-1,1))
p1
外显子2
外显子之间一般都是有连线的，所以下一步就把这些线加上去
p1 + geom_segment(data=gene_df, aes(x=start,xend=end,y=0,yend=0))
外显子3
然而，这个连线居然是在外显子的上面，一点都不美观。一种解决方式，就是分别在外显子间加线段，另一种解决方法就是先画线，然后用长方形的图层覆盖线。
p1 <- ggplot(exons_in_gene_df) +
geom_segment(data=gene_df, aes(x=start,xend=end,y=0,yend=0)) +
geom_rect(aes(xmin=start, xmax=end,ymin=-0.1,ymax=0.1),
colour="black", fill="white") +
ylim(c(-1,1))
p1
更好看的外显子
最后我们在加上这个基因的名字，以及把填一个比较好看的颜色
p1 <- ggplot(exons_in_gene_df) +
geom_segment(data=gene_df, aes(x=start,xend=end,y=0,yend=0)) +
geom_rect(aes(xmin=start, xmax=end,ymin=-0.1,ymax=0.1),
colour="#282a73", fill="#282a73") +
ylim(c(-1,1)) +
geom_text(data=gene_df, aes(x=(start + end)/2, y = 0.2,label=gene_id))
最后结果
