Biostars beta testing.
Question: How to annotate positions using biomaRt R package?
0
Entering edit mode

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 +
ADD COMMENTlink 4.0 years ago abc • 30 • updated 17 months ago RamRS 21k
Entering edit mode
0

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 REPLYlink 4.0 years ago
Sean Davis
25k
Entering edit mode
0

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

ADD REPLYlink 4.0 years ago
abc
• 30
Entering edit mode
0

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 REPLYlink 4.0 years ago
cbio
• 410
• updated 17 months ago
RamRS
21k
Entering edit mode
0

No, just I need annotation base on Ensembl

ADD REPLYlink 4.0 years ago
abc
• 30
Entering edit mode
0

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

ADD REPLYlink 4.0 years ago
abc
• 30
Entering edit mode
1

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 REPLYlink 4.0 years ago
cbio
• 410
• updated 17 months ago
RamRS
21k
Entering edit mode
0

Thanks, I face to this error

Error in file.exists(peak) : invalid 'file' argument
ADD REPLYlink 4.0 years ago
abc
• 30
• updated 17 months ago
RamRS
21k
1
Entering edit mode

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 COMMENTlink 4.0 years ago Emily_Ensembl 18k • updated 17 months ago RamRS 21k
Entering edit mode
0

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

ADD REPLYlink 4.0 years ago
abc
• 30
Entering edit mode
1

biomaRt can't work at that level of detail.

ADD REPLYlink 4.0 years ago
Emily_Ensembl
18k

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.0