Multifractal Scaling of Hi-C Maps: From Contact Probability to the Mass Exponent ✛
Context
The spatial organization of chromatin is one of the key unsolved problems in molecular biology. Hi-C experiments Hi-C, akin to a kind of interaction spectroscopy, measure contact frequencies between loci. The observed power-law decay P(s)∼$s^{−γ}$ with γ≈1 indicates a 'crumpled globule' conformation. Simultaneously, multifractal scaling was discovered in the same maps, but its origin remained a mystery. Our study fills this gap, showing that multifractality naturally emerges from the geometry of the contact map and the form of P(s), without invoking cascade processes.
Methods
The authors used numerical simulations of synthetic Hi-C maps based on a standard model with one or two power-law exponents, and analyzed experimental data from human, mouse, and six Australian species. Multifractal analysis followed standard procedure: the map was covered with squares of size l, the mass $μ_l$ in each square was computed, and the generalized partition sum Z(q,l)=∑$μ_l^q$. The exponent τ(q) was determined from the scaling Z(q,l)∼$l^{τ(q)}$ using a finite-difference derivative. The method also yields generalized dimensions $D_q$, related to Rényi entropy. Analytical contributions from diagonal (intra-segment) and off-diagonal (inter-segment) contacts were derived, leading to a piecewise linear τ(q) and linking its slopes to γ. Robustness was tested by adding white noise.
Results
Experimental τ(q) for human chromosome 1 (γ≈1.09) showed two linear regimes: a slope of 2 for small q and 1 for large q, exactly matching the theoretical prediction τ(q)=q−1 for γ≥1. For a smaller window l=200 kb, where effective γ≈0.7, the slope was ∼1.3 (=2−γ), confirming the formula τ(q)=q(2−γ)−1. Analysis of maps from mouse and six other species, including the New Zealand fur seal, black swan, and subterranean clover, showed the same pattern, underscoring the universality of multifractal scaling. Noise did not disrupt the behavior as long as the power-law P(s) held.
Implications
For the first time, a direct link is established between multifractal analysis of Hi-C maps and fundamental contact exponents of polymer physics. The slope of τ(q) at large q directly yields the bulk contact exponent $β_b$, previously extracted only from fits to P(s). This provides a new, noise-robust method for estimating polymer exponents that does not require structural models. Moreover, it is shown that multifractality here is reducible to two-point statistics, unlike in turbulent systems, and is not a signature of hidden cascades.
Future development
The proposed approach can be extended to chromatin dynamics, for example, during the cell cycle, and to joint analysis with spatial microscopy data. An intriguing direction remains extracting the surface contact exponent $β_s$ from τ(q), currently obscured by the map geometry.
Impact
These results impact chromatin biophysics, offering a new tool for interpreting Hi-C data. The method can be adapted for analyzing other network or matrix data with power-law correlations, including social and cosmological networks.
Next steps
The analysis should be applied to higher-resolution Hi-C maps and under various biological conditions, such as knockout of genes involved in loop extrusion.
Key open problems
This study connects the problem of extracting three-dimensional chromatin organization from two-point contacts with universal scaling laws. It also raises the question of the exact relationship between γ and the spatial fractal dimension d_f, especially in the regime γ<1, where d_f>3.
The knowledge machine's view
This section is addressed to the author of the paper and their colleagues. Neighbouring works were found by meaning-based search across our archive, not by word overlap. We do not review: we show the field around the work.
What the work does
The paper derives a connection between the multifractal scaling of Hi-C maps and the power-law contact probability P(s)~s^{-γ}. The authors show that multifractality arises from the geometry of the contact map and the form of P(s), not from cascade processes. Tests on synthetic and experimental data (human, mouse, six Australian species) confirm universality: the slopes τ(q) are determined by the exponent γ.
Where its strength is
A direct analytical link is established between multifractal analysis and fundamental polymer exponents, providing a new, noise-robust method for estimating the bulk contact exponent without structural models.
Where it could go next
- It might be interesting to compare the reduction of multifractality to two-point statistics with the analysis of chaos in disordered systems: the work on quantum unpredictability in crystals examines how disorder generates complex behavior, and perhaps similar ideas could help understand under what violations of the power law P(s) multifractality ceases to be trivial.based on: Chaos in a Crystal: How Disorder Awakens Quantum Unpredictability arXiv:2607.03294
- A further development could be applying the proposed method to time series or spectra where multifractality is also observed: the work on the atomic nucleus as a choir of voices discusses multifractal analysis of energy spectra, and transferring your approach, which links multifractality to simple statistical laws, could provide a new interpretation.based on: The Atomic Nucleus: Not a Soloist, but a Choir of Voices arXiv:2605.01308
- Since your work emphasizes the distinction from turbulent cascades, it might be worth considering a connection with gravitational wave signals: the analysis of the four layers of black holes uses hierarchical structures, and comparing multifractal spectra could reveal common features in the organization of complex systems.based on: Four Layers of Black Holes: Gravitational Waves Reveal Merger History arXiv:2607.22011v1
Worth trying
- Test whether the predicted piecewise linearity of τ(q) holds when realistic noise, such as Poisson noise, is added instead of just white noise.
- Apply the method to Hi-C maps with known structural anomalies (e.g., in cancer translocations) to assess the sensitivity of the bulk contact exponent.
- Compare the obtained exponents with results from other methods, such as loop or TAD analysis, for cross-validation.
Why this may matter
Directly extracting polymer exponents from multifractal analysis could simplify and unify Hi-C data processing, and open the way to studying the universality of chromatin folding across species.
Nearby works in our archive
- How magnetic particles get entangled on schedulesimilarity 0.57
- Twisted Solar Loop Reveals the Secret of Magnetic Explosionssimilarity 0.57
- Quantum Map Detectives Find Hidden Connectionssimilarity 0.56
- The Atomic Nucleus: Not a Soloist, but a Choir of Voicessimilarity 0.56
- Weak Molecular Bonds Step Out of the Shadowssimilarity 0.56
- Four Layers of Black Holes: Gravitational Waves Reveal Merger Historysimilarity 0.55
- Chaos in a Crystal: How Disorder Awakens Quantum Unpredictabilitysimilarity 0.55
The archive is sparse around this work: the closest review is only 0.57, against a usual 0.64. It looks like few have drilled nearby.
Original abstract
as written by the authors arXiv:2606.31997
The three-dimensional organization of the genome exhibits rich, scale-dependent structure, as revealed by both chromosome contact maps (e.g., Hi-C maps) and chromatin density measured by microscopy. Recent studies have reported multifractal scaling in these data. Yet, the origin of this scaling behavior remains unclear: existing efforts describe it through postulated models. Here, we show that the multifractal structure of Hi-C maps is a direct consequence of the power-law contact probability $$P(s)$$, which is itself an empirical observable measured from Hi-C maps. Starting from $$P(s)$$ with a single exponent $$\gamma$$, we analytically derive the mass exponent $$\tau(q)$$, which characterizes how the $$q$$-th moment of contact density scales with box size $$l$$ used to coarse-grain the genomic coordinate. This multifractal behavior reflects the geometric competition between intra- and inter-segment contacts. We find that the slope of $$\tau(q)$$ at large $$q$$ is given by $$2 -\gamma$$ when $$\gamma <1$$, and by $$1$$ when $$\gamma \geq 1$$. We further show that this behavior is robust to noise and consistent across diverse organisms, indicating that it is a universal feature of chromatin organization. We extend our analysis into double-exponent $$P(s)$$, and show the $$l$$ dependence in multifractal behavior. Taken together, these results provide a physical explanation for multifractal scaling and establish a direct link between the multifractality in Hi-C maps and polymer contact statistics, with the large-$$q$$ slope of $$\tau(q)$$ mapping onto a known polymer contact exponent.
The multifractal formalism, born in turbulence theory, unexpectedly became a key to understanding DNA folding. In the genome, however, multifractality turned out to be simpler: it is entirely determined by two-point correlations, without a cascade hierarchy.
Key numbers
- γ exponent human Chr1: 1.09
- γ exponent mouse Chr1: 1.1
- τ(q) slope at small q: 2
- Hi-C map resolution: 25 kb
- number of species studied: 8