代码库8-获得hg38启动子并添加gene symbol

library(TxDb.Hsapiens.UCSC.hg38.knownGene)
library(org.Hs.eg.db)

txdb <- TxDb.Hsapiens.UCSC.hg38.knownGene
gene_gr <- genes(txdb)
promoter_gr <- promoters(gene_gr)
promoter_gr$gene_symbol <- mapIds(x = org.Hs.eg.db,
                                  keys = promoter_gr$gene_id, 
                                  keytype = "ENTREZID", #需要转换的类型
                                  column = "SYMBOL") #需要转换为的类型
head(promoter_gr)
©著作权归作者所有,转载或内容合作请联系作者
【社区内容提示】社区部分内容疑似由AI辅助生成，浏览时请结合常识与多方信息审慎甄别。
平台声明：文章内容（如有图片或视频亦包括在内）由作者上传并发布，文章内容仅代表作者本人观点，简书系信息发布平台，仅提供信息存储服务。

相关阅读更多精彩内容

友情链接更多精彩内容