Question: how to obtain the locations of a exon and its corresponding peptide from a specific transcript in bioMart
0
gravatar for jjaeyun
4.0 years ago by
jjaeyun0
United States
jjaeyun0 wrote:

I tried to connect a transcript (ENSTxxx)'s exons to their corresponding peptides (residues) using bioMart in R.

library("biomaRt")
ensembl<-  useMart("ensembl")
ensembl<- useDataset( "ensembl", dataset="hsapiens_gene_ensembl")
getBM(c('peptide_location','exon_chrom_start','exon_chrom_end'), filters = c('ensembl_transcript_id','pdb'), values=list('ENST00000606293','2ZA5'), mart=ensembl)
Error in getBM(c("peptide_location", "exon_chrom_start", "exon_chrom_end"),  : 
  Query ERROR: caught BioMart::Exception::Usage: Attributes from multiple attribute pages are not allowed

I wish to get some help for this problem.
peptide exon biomart • 1.9k views
ADD COMMENTlink modified 4.0 years ago by informatics bot570 • written 4.0 years ago by jjaeyun0
1
gravatar for informatics bot
4.0 years ago by
United States
informatics bot570 wrote:

This has to do with the way the Biomart webserver/database is set up. Essentially, you can't do gene-based queries and transcript-based queries at one time. But you can do them sequentially and merge at the end.

transcript_data<-getBM(c('ensembl_transcript_id','exon_chrom_start','exon_chrom_end'), filters = c('ensembl_transcript_id','pdb'), values=list('ENST00000606293','2ZA5'), mart=ensembl)

pept_data<-getBM(c('ensembl_transcript_id','peptide_location'), filters = c('ensembl_transcript_id','pdb'), values=list('ENST00000606293','2ZA5'), mart=ensembl)

final_data<-merge(transcript_data, pept_data, by='ensembl_transcript_id')
ADD COMMENTlink modified 4.0 years ago • written 4.0 years ago by informatics bot570

Thank you so much! It worked. And, it answers for my rest of work. Thanks again!

BTW, this was my first use of Biostart. This is pretty nice

ADD REPLYlink written 4.0 years ago by jjaeyun0

You're very welcome :D

feel free to up-vote my "answer" and mark this question as solved.

ADD REPLYlink written 4.0 years ago by informatics bot570

Hi Lando, 

I thought your code was what I wanted, but actually, it is not.

In your merge function at the end, it was supposed to be exon range or exon id, other than transcript id. Could you correct this?

Thanks,

ADD REPLYlink written 4.0 years ago by jjaeyun0
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1935 users visited in the last hour