如何在bsseq的BSmooth结果图下方添加基因结构注释信息？

核心思路 #

利用bsseq本身的annoTracks参数，配合基因组注释工具包导入的基因、启动子、外显子区域数据，实现自动注释添加，无需手动绘图。

具体步骤 #

1. 导入并整理基因注释数据
   使用GenomicFeatures包从GTF/GFF文件提取所需的基因结构信息：

   library(GenomicFeatures)
   # 从GTF文件生成TxDb注释对象（替换为你的注释文件路径）
   txdb <- makeTxDbFromGFF("genome_annotation.gtf", format = "gtf")
   # 提取基因、外显子区域，自定义启动子范围（上游2000bp+下游200bp）
   genes <- genes(txdb)
   exons <- exons(txdb)
   promoters <- promoters(txdb, upstream = 2000, downstream = 200)
   

2. 通过annoTracks参数添加注释到plotRegion
   将整理好的注释区域打包成列表，传入plotRegion的annoTracks参数，同时可以通过annoTracksArgs定义各注释轨道的样式：

   library(bsseq)
   # 构建注释轨道列表，命名对应注释类型
   anno_tracks <- list(
     启动子 = promoters,
     基因 = genes,
     外显子 = exons
   )
   # 定义各轨道的绘图样式（颜色、线宽、填充色等）
   anno_styles <- list(
     启动子 = list(col = "#e74c3c", lwd = 2),
     基因 = list(col = "#2980b9", fill = "#3498db", alpha = 0.3),
     外显子 = list(col = "#2ecc71", fill = "#2ecc71", alpha = 0.6)
   )
   # 调用plotRegion生成带注释的图（替换为你的bsseq对象和目标DMR区域）
   plotRegion(
     bs = your_bsseq_object,
     region = target_dmr_grange,
     annoTracks = anno_tracks,
     annoTracksArgs = anno_styles,
     main = "DMR区域甲基化水平及基因注释",
     plotType = "p"
   )
   

3. 进阶：用Gviz实现更精细的基因结构展示
   如果需要显示内含子-外显子的连接结构，可结合Gviz包，在plotRegion绘图后叠加基因轨道：

   library(Gviz)
   # 创建基因结构轨道（自动解析外显子和内含子）
   gene_track <- GeneRegionTrack(
     txdb,
     chromosome = seqnames(target_dmr_grange),
     from = start(target_dmr_grange),
     to = end(target_dmr_grange),
     name = "基因结构",
     fill = "#2980b9",
     transcriptAnnotation = "transcript"
   )
   # 先绘制甲基化水平图
   plotRegion(your_bsseq_object, region = target_dmr_grange, plotType = "l")
   # 叠加基因结构轨道
   plotTracks(gene_track, add = TRUE, from = start(target_dmr_grange), to = end(target_dmr_grange))
   

内容的提问来源于stack exchange，提问作者Charité Learner