使用"biomaRt"包進行基因ID轉換

1、安裝"biomaRt"包

if (!requireNamespace("BiocManager", quietly = TRUE))
    install.packages("BiocManager")
BiocManager::install("biomaRt")
library("biomaRt")

2、"ensembl_transcript_id_version" 轉 "refseq_mrna", "ensembl_gene_id", "hgnc_symbol"

輸入文件file.txt如下所示：

轉換命令如下：

values<- read.table("file.txt",header=F)
ensembl<-  useMart("ensembl", dataset="hsapiens_gene_ensembl")
new=getBM(attributes=c("ensembl_transcript_id_version","refseq_mrna", "ensembl_gene_id", "hgnc_symbol"), filters = "ensembl_transcript_id_version", values = values, mart= ensembl)

輸出文件new如下所示：