Hi!
I have a list of peaks and I would like to annotate them to the nearest genes using ensemble build say 69 or 70.
The desired output is to have the nearest gene and also to know if the peak falls in exon, intron, intergenic or promoter region.
How can I do it?
Kindly help
Thank you
4 answers
Different tools, check out Galaxy Web Interface, If you wanna do offline, try CEAS.
What I do is pull the tables from UCSC(see posts like A: How to download all the introns from UCSC, Format RefSeq output obtained from UCSC Table Browser?) and then use intersectBed from Bedtools. check the usage here A: how to determine overlaps from coordinates. There are lot of posts discussing about this, please skim through archive quickly.
•
written
by
Sukhi Singh
If you are R/Bioconductor inclined, the locateVariants function in the VariantAnnotation package does something similar to what you are trying to do.
I suggest using, ChIPpeakAnno a package in Bioconductor, R. Here is a link to the paper which talks about the package. It gives you a visual analysis of all the regions near to the peak as well.
Log in to answer this question.
-
I am getting more annotation hits from HOMER than I have actual peaks. Is that normal?
written by Ronin •
-
How to annotate probe to nearest gene promoter?
written by BioNoobie •
-
How to filter out the promoter-TSS peaks from the `.narrowPeak` file?
written by Dan •
-
How to map query coordinates on to genome gff annotation files ?
written by isha.lily20 •
-
annotationPeaks.pl algorithm for defining promoter-Tss
written by yahelsal •
-
Defining exons, UTRs, introns and intergenic regions using "overlap" in Ensemble REST API
written by Maria •
-
Tool: Tool to annotate ChIP-Seq peaks (hg19 or mm10) and find neighboring peaks from multiple exper…
written by abhishek.mitra •
-
Generate Random Genomic Regions
written by ChIP
-
Annotate Chip-Peaks With Gene Symbol/Name
written by GeneInfo •
-
Annotatepeaks Function From Homer
written by Dataminer
