TCGA 数据库中的基因编号采用的Esembl 的编号,但是有些分析软件,需要输入的基因编号是 gene symbol ,这就需要将Esemble 的ID 转换成gene symbol 。
今天介绍采用clusterProfiler 进行转换:
# 加载相关软件包 > library(clusterProfiler) > library(org.Hs.eg.db) # org.Hs.eg.db 包提供的ID转换类型 > keytypes(org.Hs.eg.db) [1] "ACCNUM" "ALIAS" "ENSEMBL" "ENSEMBLPROT" "ENSEMBLTRANS" "ENTREZID" "ENZYME" "EVIDENCE" "EVIDENCEALL" [10] "GENENAME" "GO" "GOALL" "IPI" "MAP" "OMIM" "ONTOLOGY" "ONTOLOGYALL" "PATH" [19] "PFAM" "PMID" "PROSITE" "REFSEQ" "SYMBOL" "UCSCKG" "UNIGENE" "UNIPROT" # 测试ID > test_id <- c("ENSG00000000971", "ENSG00000001084", "ENSG00000001460", "ENSG00000001461", "ENSG00000001626", "ENSG00000001630") # 采用bitr 命令进行ID的转换 > gene_ids <- bitr(test_id, fromType="ENSEMBL", toType=c("SYMBOL", "GENENAME"), OrgDb="org.Hs.eg.db") 'select()' returned 1:1 mapping between keys and columns # 查看转换的结果 > gene_ids ENSEMBL SYMBOL GENENAME 1 ENSG00000000971 CFH complement factor H 2 ENSG00000001084 GCLC glutamate-cysteine ligase catalytic subunit 3 ENSG00000001460 STPG1 sperm tail PG-rich repeat containing 1 4 ENSG00000001461 NIPAL3 NIPA like domain containing 3 5 ENSG00000001626 CFTR cystic fibrosis transmembrane conductance regulator 6 ENSG00000001630 CYP51A1 cytochrome P450 family 51 subfamily A member 1
如果您对TCGA数据挖掘感兴趣,请学习我的TCGA系列课程:
如果觉得我的文章对您有用,请随意打赏。你的支持将鼓励我继续创作!