Can I use two restriction enzymes to analyze HiC data using HiCExplorer
2
0
Entering edit mode
2.5 years ago

Hi My HiC library was prepared using two restriction enzymes. The HiCExplorere tool findRestSites seems to take only one enzyme site as input. Is there any way I can use two restriction enzymes and analyze my data using HiCExplorer.

Thank you!

HiC HiCExplorer Arima HiC Multiple enzyme HiC • 2.0k views
1
Entering edit mode
2.5 years ago
harish ▴ 350

I've had this same issue a while back. I ended up using genome_disgest.py from HiC-Pro to do that. (https://github.com/nservant/HiC-Pro)

0
Entering edit mode

Thank you Harish! I was able to generate a .bed file using digest_genome.py.

Yogi

0
Entering edit mode
2.5 years ago

I was able to generate the list of restriction fragments (bed file) after genome digestion using two restriction enzymes used in Arima HiC kit. I used HiC-Pro digest_genome.py to do that (https://github.com/nservant/HiC-Pro). Arima kit uses two restriction enzymes ^GATC and G^ANTC. Since digest_genome.py does not allow for "N" in the cut sites. I used following command to get around this problem:

digest_genome.py -r ^GATC G^AATC G^ATTC G^AGTC G^ACTC -o your_output.bed your_genome.fa


Yogi Bordiya

1
Entering edit mode

Ahh, the latest update does support the ^N sites, I found about it on github.

However you will definitely need to process your files in a compatible format if you are using Juicer etc.