Biology Stack Exchange is a question and answer site for biology researchers, academics, and students. Join them; it only takes a minute:

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

I am working on a project and I'm trying to find data on the position of genes and heterochromatin binding sites such as HP1 in Drosophila melanogaster. Is this information available for the DGRP lines? Would I be able to measure the distance between a gene and its nearest binding site from such data?

share|improve this question

Have a look at these ChIP-seq data for HP1 in Drosophila: 1, 2 and 3. From ChIP-seq data you can find the distance between the TFBS peaks and the TSS of the gene.

You can also look for nucleosome positioning and DNAse hypersensitvity regions; for the former, I am sure that data is available for Drosophila.

share|improve this answer

As the previous answer said, there are some public HP1 ChIP-seq data from D. melanogaster if not from the DGRP, from modENCODE and maybe others. In the case of modENCODE, they've published not only the reads, but also their peak calls (mapping with Eland + calling with MACS).

BEDTools ( )is a nice command line tool for manipulating intervals with genomic coordinates, and the bedtools closest command can calculate the distance to the closest feature. It would look something like this:

bedtools closest -d -a genomic_features.gff3 -b peak_calls.gff3

Where the "genomic_features" part would be any subset of genome annotation that you're interested in.

I've worked through a small example with modENCODE data, Flybase transcript annotation and BEDTools:



share|improve this answer

Your Answer


By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.