Question: Overlap chipseq regions with repetative regions
0
gravatar for chipseq33
3.7 years ago by
chipseq3310
chipseq3310 wrote:

Hi guys,

I have a list of ChIP regions and I would like to know

1) which of them overlap with RepeatMask table that I took through UCSC?

chr1 10000 10468 + Simple_repeat chr1 10468 11447 - Satellite chr1 11503 11675 - LINE chr1 11677 11780 - DNA chr1 15264 15355 - SINE chr1 16712 16749 + Simple_repeat

2) how much overlap do you consider to say that one given region is located in repetitive regions?

Thanks in advance chip

chipseq • 1.1k views
ADD COMMENTlink modified 3.7 years ago by jotan1.2k • written 3.7 years ago by chipseq3310
1
gravatar for James Ashmore
3.7 years ago by
James Ashmore2.7k
UK/Edinburgh/MRC Centre for Regenerative Medicine
James Ashmore2.7k wrote:

You can use the command-line tool bedtools for this analysis. Use the intersect command to look for overlaps between your ChIP regions and the regions in the RepeatMask table. You can also specify a minimum overlap required for the region to count as an overlap, in the example below I've set the minimum overlap to be 50% of a ChIP regions length:

bedtools intersect -a ChIP_regions.bed -b RepeatMask_table.bed -f 0.5
ADD COMMENTlink modified 3.7 years ago • written 3.7 years ago by James Ashmore2.7k
1
gravatar for jotan
3.7 years ago by
jotan1.2k
Australia
jotan1.2k wrote:

If you are interested in identifying ChIP-seq enrichment at repeats, you should take a look at this paper from the Kharchenko lab:

http://genomebiology.biomedcentral.com/articles/10.1186/gb-2010-11-6-r69

Source code is available here:

http://compbio.med.harvard.edu/repeats/

This will give a more accurate idea of repeat enrichment compared to mapping reads across the genome.

ADD COMMENTlink written 3.7 years ago by jotan1.2k
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: 764 users visited in the last hour