This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Best way to get gene IDs for Salmon transcript output

I generated TPM counts from fastq data using salmon. This leaves me with the NM_transcript IDs. I would like to generate the gene symbols from the transcript IDs . Biomart does not recognize transcripts, NCBI Datasets produces an error when I run the entire transcriptome. I have been exploring tximport and tximeta, however, I have run into numerous issues particularly with tximeta not detecting my ref file. Any advice would be greatly appreciated.

Update: I have txiimport and tximeta now running, however they create S4 objects and I am unsure how to make these readable

rna-seq r
• 2,793 views • link
written by dk0319 • dk0319

2 answers

tximeta was able to compile all my quant.sf files and summarize to gene level

• 0 views • link
written by dk0319 • dk0319

I would like to generate the gene symbols from the transcript IDs

If you need just the gene symbols, and not the sequence, you can parse the gene2refseq.gz file. Depending on the age of your input set of accessions, you may not find information for all of them. That's because gene2refseq.gz file is regularly updated and any NM_ accessions that are no longer latest will be absent. You can download the gene2refseq.gz file from here: https://ftp.ncbi.nlm.nih.gov/gene/DATA/

As a first pass, you can extract information for as many NM_ accessions as you can from this file and then use Entrez Direct or NCBI Datasets to get the information for the remaining ones.

• 0 views • link
written by vkkodali_ncbi vkkodali_ncbi

Log in to answer this question.

More posts like this