Biostar Test Site

This is site is used for testing only. Visit: https://www.biostars.org to ask a question.

add gene names from RefSeq instead of ENSEMBL using biomaRT for sleuth when reading kallisto output
1
0
Entering edit mode
3.3 years ago
Abiterkuile ▴ 20

Hi there,

I have reads that have been aligned using Kallisto constructed with a Refseq reference transcriptome (ftp://ftp.ncbi.nlm.nih.gov/refseq/H_sapiens/annotation/GRCh38_latest/refseq_identifiers/GRCh38_latest_rna.fna.gz). In order to include gene names into transcript-level analysis in Sleuth (an R package), I need to add the gene names using a Reqseq human trancriptome however in the Sleuth walk-through there are only instructions on how to do this using ENSEMBL (for an example data set constructed with the ENSEMBL human transcriptome);

mart <- biomaRt::useMart(biomart = "ENSEMBL_MART_ENSEMBL",
  dataset = "hsapiens_gene_ensembl",
  host = 'ensembl.org')
t2g <- biomaRt::getBM(attributes = c("ensembl_transcript_id", "ensembl_gene_id",
    "external_gene_name"), mart = mart)
t2g <- dplyr::rename(t2g, target_id = ensembl_transcript_id,
  ens_gene = ensembl_gene_id, ext_gene = external_gene_name)

from: https://pachterlab.github.io/sleuth_walkthroughs/trapnell/analysis.html

Does anyone know how i can do this but with RefSeq?

Thanks

RNA-Seq alignment sleuth biomaRt RefSeq • 2.1k views
ADD COMMENT
1
Entering edit mode
3.3 years ago

RefSeq don't really have gene names but you can get the NM transcript IDs with:

t2g <- biomaRt::getBM(attributes = c("refseq_mrna", "external_gene_name"), mart = mart)
ADD COMMENT

Login before adding your answer.

Traffic: 314 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6