Question: Mapping CLIP-seq Regions to 5'UTR, CDS, or 3'UTR
0
gravatar for mb314
10 months ago by
mb3140
mb3140 wrote:

Dear Community,

I have a list of peaks (bed format) from my CLIP-seq project and I'm trying to map where each of these peaks are located with respect to an mRNA (5' or 3' UTR, or CDS). Is there a program I can use to map this all at once rather than manually?

Thank you!

3'utr 5'utr clip-seq • 339 views
ADD COMMENTlink modified 10 months ago • written 10 months ago by mb3140
0
gravatar for Pierre Lindenbaum
10 months ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum126k wrote:

bedtools intersect with a gff file.

ADD COMMENTlink written 10 months ago by Pierre Lindenbaum126k

Thanks! I did that, however I only got the ID, geneID, and gene_name (for example: ID=NM_002032;geneID=FTH1;gene_name=FTH1). In your experience have you gotten the location of the peak within the coding transcript (i.e. 5'UTR/3'UTR etc). It could be an issue with my gff file.

ADD REPLYlink written 10 months ago by mb3140
0
gravatar for Giovanni Quinones Valdez
10 months ago by
UCLA

Annovar is a software that typically annotates variants, but you can use it to annotate your CLIP-Seqs. It's very smart when assigning ambiguous genomic regions.

ADD COMMENTlink written 10 months ago by Giovanni Quinones Valdez0
0
gravatar for mb314
10 months ago by
mb3140
mb3140 wrote:

Following up in case someone has this question, I ended up using HOMER's annotatePeaks.pl. It was very easy to use with my BED6 files.

More information here: http://homer.ucsd.edu/homer/ngs/annotation.html

ADD COMMENTlink written 10 months ago by mb3140
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1832 users visited in the last hour