如何將TCGA的ENSG編號(hào)轉(zhuǎn)換成gene symbol

這篇文章給大家分享的是有關(guān)如何將TCGA的ENSG編號(hào)轉(zhuǎn)換成gene symbol的內(nèi)容。小編覺得挺實(shí)用的，因此分享給大家做個(gè)參考，一起跟隨小編過來看看吧。

TCGA的ENSG編號(hào)轉(zhuǎn)換成gene symbol

TCGA 數(shù)據(jù)庫中的基因編號(hào)采用的Esembl 的編號(hào)，但是有些分析軟件，需要輸入的基因編號(hào)是 gene symbol ，這就需要將Esemble 的ID 轉(zhuǎn)換成gene symbol  。 
采用clusterProfiler 進(jìn)行轉(zhuǎn)換：

# 加載相關(guān)軟件包
> library(clusterProfiler)
> library(org.Hs.eg.db)

# org.Hs.eg.db 包提供的ID轉(zhuǎn)換類型
> 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 命令進(jìn)行ID的轉(zhuǎn)換
> 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

# 查看轉(zhuǎn)換的結(jié)果
> 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

感謝各位的閱讀！關(guān)于“如何將TCGA的ENSG編號(hào)轉(zhuǎn)換成gene symbol”這篇文章就分享到這里了，希望以上內(nèi)容可以對大家有一定的幫助，讓大家可以學(xué)到更多知識(shí)，如果覺得文章不錯(cuò)，可以把它分享出去讓更多的人看到吧！