Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

TAD missing #916

Open
80QvQ08 opened this issue Dec 1, 2024 · 0 comments
Open

TAD missing #916

80QvQ08 opened this issue Dec 1, 2024 · 0 comments

Comments

@80QvQ08
Copy link

80QvQ08 commented Dec 1, 2024

Hello, hicexplorer is amazing! I successfully used it to to analyze A/B compartment and call TADs, it was really helpful, but i encountered a problem recently. For example, as the pictures show, TADs were not called on some area of seveal chromsomes in different species.
f6a484b367445b5837e7bc1eb50a057
e0724e678295e3a488e3dcc4ea1bdc0

I think it is not the problem of pygenometrack, because I checked the TAD_boundaries.bed and TAD_domains.bed, the relevant information is indeed missing for these areas, for example, the length of chromsome in first picture is 19032923bp, the domain information of this chromsome is below:
4c46494027a71c5a2c0007cbeade885

The version of hicexplorer i used is the latest one 3.7.6 installed by conda, and the version of python is 3.8.
The code I used and the feedback are listed below:

hicFindTADs -m out.50k.cool.corrected.h5 --outPrefix out.50k.cool.corrected.h5.TAD --thresholdComparisons 0.05 --delta 0.01 --numberOfProcessors 50 --correctForMultipleTesting fdr
hicPlotTADs --tracks tracks.ini --region Chr11:0-19032923 -o Mpha_TAD_50X_chr11_all.png 
INFO:pygenometracks.tracksClass:initialize 1. [x-axis]
INFO:pygenometracks.tracksClass:initialize 2. [hic matrix]
WARNING:hicmatrix.HiCMatrix:Bin size is not homogeneous.                                       Median 50000.000000

INFO:pygenometracks.tracksClass:initialize 3. [tads]
100%|█████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 17/17 [00:00<00:00, 13850.65it/s]
INFO:pygenometracks.tracksClass:time initializing track(s):
INFO:pygenometracks.tracksClass:2.6904544830322266
DEBUG:pygenometracks.tracksClass:Figure size in cm is 40 x 2.6697970769549277. Dpi is set to 72

INFO:pygenometracks.tracksClass:plotting 1. [x-axis]
INFO:pygenometracks.tracksClass:plotting 2. [hic matrix]
WARNING:pygenometracks.tracks.GenomeTrack:*Warning*
The region to plot extends beyond the chromosome size. Please check.
Chr11 size: 19032923. Region to plot 0-49610884

INFO:pygenometracks.tracks.GenomeTrack:setting min, max values for track 2. [hic matrix] to: 1.0, 50.0

INFO:pygenometracks.tracksClass:plotting 3. [tads]
DEBUG:pygenometracks.tracks.GenomeTrack:ylim (0.0, 900000.0)

Meantime, i also found that althought other chromsomes seemed normal, but both ends of these are clear, Can hicexplorer not identify TADs at both ends of the chromosome?
c995d579896cd14ec6501dbb9b07df1
b5ee0dfa89746203d1660eeca9fec9e

Thank you very much for your time and assistance. I look forward to your guidance on this matter.

best

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant