How to annotate positions using biomaRt R package?
1
0
Entering edit mode
6.4 years ago
Ati ▴ 40

I have some genome positions and I want to annotate these positions (find Ensembl gene ID, the features like exonic, intronic, ...) using biomaRt R package.

part of my data

chr start stop strand
chr10 100572320 100572373 - 
chr10 100572649 100572658 +
r bioinformatics biomart • 3.0k views
ADD COMMENT
0
Entering edit mode

What do these regions represent?  Doing what you suggest with biomaRt is possible, but it is going to be very tedious.  Perhaps there is another tool that is a better fit, but we need to know a little more detail before suggesting.

ADD REPLY
0
Entering edit mode

There Polyadenylation sites and I want to know if they are intronic (annotate for introns)

ADD REPLY
0
Entering edit mode

Does this have to be done using biomaRt as part of a pipeline? If not you could look into the ChIPseeker package for straight-forward annotation.

ADD REPLY
0
Entering edit mode

No, just I need annotation base on Ensembl

ADD REPLY
0
Entering edit mode

This package need specific input format that I don't have

ADD REPLY
1
Entering edit mode

The input is a bedfile that can be read in by using the readPeakFile function, followed by the annotatePeak function. That should give you what you're looking for.

The example you've shown is basically a bedfile, so I'm not sure what the problem is.

ADD REPLY
0
Entering edit mode

Thanks, I face to this error

Error in file.exists(peak) : invalid 'file' argument
ADD REPLY
1
Entering edit mode
6.4 years ago
Emily 23k

Something like:

gene_list <-getBM(
              attributes=c(
                     'ensembl_gene_id',
                     'ensembl_transcript_id'
                     ),
              filters='chromosomal_region',
              values="10:100572320:100572373:-1,10:100572649:100572658:1",
              mart = ensembl
              )
ADD COMMENT
0
Entering edit mode

Thanks but for find intronic or exonic which attributes should be used?

ADD REPLY
1
Entering edit mode

biomaRt can't work at that level of detail.

ADD REPLY

Login before adding your answer.

Traffic: 1779 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