This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to distinguish protein isoforms using biomaRt?

I want to retrieve the canonical sequence for some proteins from ensembl, based upon their UniProt IDs.

library(biomaRt)

ensemblMart <- useMart(
  biomart = "ensembl",
  dataset = "hsapiens_gene_ensembl"
)

(result <- getBM(
  attributes = c("uniprot_swissprot", "peptide"),
  filters    = "uniprot_swissprot",
  values     = "P31749",
  mart       = ensemblMart,
  bmHeader   = TRUE
))

This particular protein has two isoforms, both of which are returned. I only care about the canonical form, so how can I distinguish the isotopes?

I suspect that I need to return another of the attributes from the mart

listAttributes(ensemblMart)

but I don't know which one.

biomart r
• 4,156 views • link
updated by Ram Ram • written by richierocks • richierocks

2 answers

There is no canonical defined in BioMart.

• 0 views • link
written by Emily Emily

Thanks. Just to be clear: are there any other ways of distinguishing isoforms?

• 0 views • link
updated by Ram Ram • written by richierocks • richierocks

If you are trying to retrieve protein sequences using UniProt IDs, why not just fetch them directly from UniProt using something like urllib2 on Python (I'm sure there's a way to do this in R as well)? Or even just wget on linux, i.e., wget -nv http://www.uniprot.org/uniprot/O53166.fasta. UniProt will return the canonical sequence.

• 0 views • link
updated by Ram Ram • written by shelly.deforte shelly.deforte

Log in to answer this question.

More posts like this