This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Biomart Refseq Data Problem

Hello,

I want to extract all the mRNA sequence that have RefSeq ID with the coding sequence, and start and end positions of 3 and 5 prime UTRs. I have tried;

EM = useMart("ensembl", dataset = "hsapiensgeneensembl")

attr = c("ensembltranscriptid", "cdna", "cdnacodingstart", "cdnacodingend", "5utrstart", "5utrend", "3utrstart", "3utrend")

Refseq = getBM(attributes = attr, filters = "withoxrefseq_mrna", values = TRUE, mart = EM, uniqueRows = TRUE)

and keep getting

Error in scan(file, what, nmax, sep, dec, quote, skip, nlines, na.strings, : line x did not have y elements (depending on variables its different lines etc)

I have tried getting ride of different variables but even just the ID, cDNA sequence, coding start and end it is still the same. What is the problem?

Secondly I am developing an R package that needs a lot of background biological data. Is it a wise move to incorporate biomaRt section in the code to get all the necessary data so the user does not have to provide it by himself? I am not sure how future proof it will be with potential problems even a minor change might cause.

Thanks for your help.

biomart
• 3,341 views • link
updated by Biojl Biojl • written by Pppp • Pppp

2 answers

The problem is that you are trying to force sequence data ("cdna") and other types of data into a data frame - which is not going to work.

So just omit "cdna" from attr and it will work fine (I just tried it).

As for the R package, you should read up on how other people do it. For example, many of the Bioconductor packages depend on one another but you don't necessarily "incorporate" code from one in another; you just ensure that it loads when required. The developers section at the Bioconductor website should have more details of best practices.

• 2 views • link
written by Neilfws Neilfws

That solves the problems, thanks very much.

• 0 views • link
written by Pppp • Pppp

And just to clarify, if you do want to fetch sequences, use the getSequence() method - it's described in the biomaRt PDF. Then you can merge the results of getSequence() and getBM() if required.

• 0 views • link
written by Neilfws Neilfws

Regarding the second part of your question: Take into account that you can only access through an API-database to the last version of ENSEMBL. Previous versions are only accessible through the website. This may be a problem, because you normally don't finish a project within a single release or have previous work with another release version.

• 1 views • link
written by Biojl Biojl

Log in to answer this question.

More posts like this