This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How To Map Ucsc Mrna Isoforms To Gene ?

I have isoform expression file of human cancer from TCGA. I want to map all this isoform to gene, in other word, I want to group the isoform which belong to one gene together. the ID of isoforms belong to UCSC. Does anybody knows how can I do this ? 

isoform_id    normalized_count
uc011lsn.1    0.0000
uc010unu.1    20.1848
uc010uoa.1    7.1561
uc002bgz.2    36.1698
uc002bic.2    0.0000
uc010zzl.1    188.5822
uc001jiu.2    1085.9445
uc010qhg.1    0.0000
uc011krn.1
ngs bioinformatician
• 13,061 views • link
updated by Chun-Jie Liu Chun-Jie Liu • written by jack • jack

Normally one could use biomart (or the biomaRt Bioconductor package). However, none of these IDs seem to be included in that :(

• 310 views • link
written by Devon Ryan Devon Ryan

In UCSC genome browser itself u can output the ucsc_ids as well as gene names . try to see the tablebrowser and the output options :)

• 303 views • link
written by ancient_learner ancient_learner

3 answers

Go to UCSC Table Browser

Select:

group: Genes and Gene Predictions

track: UCSC Genes

table: hgFixed.transMapGeneUcscGenes

identifiers: paste list: (paste the list of isoform ids)

Click "get output"

screenshot

You will have something like:

#id    cds    db    geneName
uc001jiu.2    123..752    hg19    TIMM23
uc002bgz.2        hg19    UBE2QP2
uc002bic.2        hg19    UBE2QP2
uc010qhg.1    123..641    hg19    TIM23
uc010unu.1        hg19    UBE2QP2
uc010uoa.1        hg19    UBE2QP2
uc010zzl.1    1..633    hg19    HMGB1L1
uc011krn.1        hg19    MOXD2
uc011lsn.1        hg19    LOC100130426

You may want to look at the UCSC Table Browser Help or go to USCS genome support forum.

• 489 views • link
written by Zhaorong Zhaorong

Use biomaRt

library(biomaRt)

GENES = useMart("ENSEMBL_MART_ENSEMBL", dataset = "hsapiens_gene_ensembl")   

getBM(attributes = c('ensembl_gene_id','hgnc_symbol','ensembl_transcript_id','refseq_mrna','ucsc','chromosome_name','transcript_start','transcript_end'), mart = GENES)

This will list all Ensembl gene/transcript id, HGNC symbol, Refseq mRNA and ucsc id mapping. You can add filter and values into getBM to get mapping limited to values length.

The UCSC ID used by FIREHOSE isoform expression data is old version. You can not use the UCSC ID as filter in getBM.

• 0 views • link
written by Chun-Jie Liu Chun-Jie Liu

Log in to answer this question.

More posts like this