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
ADD COMMENT
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)

ADD COMMENT
0
Entering edit mode

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

Yogi

ADD REPLY
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

ADD COMMENT
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.

ADD REPLY

Login before adding your answer.

Traffic: 1563 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6