🔬 condensed matter

Multifractal Scaling in Hi-C Maps

This paper demonstrates that the multifractal scaling observed in Hi-C maps is a direct analytical consequence of the power-law contact probability $P(s)$ , establishing a universal physical link between chromatin organization's geometric competition and polymer contact statistics.

Original authors: Seong-Gyu Yang, Lucas Hedström, Jan Smrek, Ludvig Lizana

Published 2026-07-01

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Seong-Gyu Yang, Lucas Hedström, Jan Smrek, Ludvig Lizana

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine your DNA not as a long, straight string of beads, but as a giant, tangled ball of yarn inside a tiny box (the cell nucleus). Scientists have developed a way to take a "snapshot" of how different parts of this yarn touch each other. This snapshot is called a Hi-C map. It looks like a giant grid where bright spots show where two parts of the DNA are hugging each other.

For a long time, scientists noticed something strange about these maps: they have a "fractal" quality. This means the pattern looks complex and self-similar, no matter how much you zoom in or out. But why they looked this way was a mystery. Some thought it was due to a complex, hierarchical folding process, like a Russian nesting doll.

This paper solves that mystery with a simple, elegant idea: The complex pattern isn't caused by a complicated folding rule; it's just a natural side effect of how likely two pieces of yarn are to touch based on how far apart they are.

Here is the breakdown using simple analogies:

1. The Two Types of "Hugs"

The authors realized that the Hi-C map is actually a competition between two types of interactions, which they call intra-segment and inter-segment contacts.

The "Diagonal" Hug (Intra-segment): Imagine looking at the main diagonal line of the map. These are touches between DNA pieces that are neighbors on the string. They are close to each other, so they hug often.
The "Off-Diagonal" Hug (Inter-segment): These are touches between pieces of yarn that are far apart on the string but happen to be close in the 3D ball. These are rarer.

2. The Magic of "Zooming In" (The Box Size)

To measure the "fractal" nature, scientists use a grid (like a checkerboard) and count how many hugs happen inside each square. They change the size of the squares (the "box size") to see how the counts change.

The paper shows that the pattern you see depends entirely on which type of hug dominates your count:

When you look at the "big picture" (Small numbers): The "Off-Diagonal" hugs (the rare, far-away touches) dominate the math. This creates one specific pattern.
When you look at the "details" (Large numbers): The "Diagonal" hugs (the frequent, neighborly touches) take over. This creates a different pattern.

3. The "Crossover"

The most important discovery is that the map doesn't just switch abruptly from one pattern to another. Because the DNA is a continuous, tangled ball, there is a smooth crossover in the middle.

Think of it like listening to a song that slowly transitions from a drum solo to a violin solo. You don't hear a sudden cut; you hear the drums fade out as the violins fade in. The paper proves that this smooth transition is what creates the "multifractal" signature. It's not two separate patterns; it's one continuous geometric competition.

4. The "Universal" Rule

The authors tested this on humans, mice, and even a wide variety of other species (seals, swans, dragons, sharks, aphids, and plants).

The Finding: No matter the species, the "hug probability" follows a simple rule: the further apart two pieces are on the string, the less likely they are to touch, following a specific mathematical curve (a power law).
The Result: Because this rule is universal, the resulting "fractal pattern" in the map is also universal. It's a fundamental feature of how DNA is organized, not a quirk of a specific animal.

5. Why This Matters (Without the Jargon)

Before this paper, if you wanted to understand the 3D shape of DNA from these maps, you had to build a complex computer model and guess the rules.

This paper says: "You don't need to guess the rules."
If you just measure how often DNA touches itself at different distances (the power law), you can mathematically predict the entire fractal pattern of the map. The complex "fractal" look is just the shadow cast by a simple rule of distance.

The "Noise" Test

The scientists also added "static" (noise) to the data, like static on a radio. They found that even with a lot of noise, the fractal pattern remained clear. This proves the pattern is a real, robust feature of the DNA's organization, not just an artifact of messy data.

Summary

In short, the paper explains that the complex, fractal-looking maps of our DNA are not the result of a mysterious, multi-layered folding machine. Instead, they are the natural, geometric result of a simple rule: DNA pieces that are close together on the string touch more often than those far apart. The "multifractal" behavior is just the mathematical signature of this simple distance rule playing out across a 2D map.

Technical Summary: Multifractal Scaling in Hi-C Maps

Problem Statement
The three-dimensional organization of the genome exhibits scale-dependent structures observable via chromosome contact maps (Hi-C) and chromatin density microscopy. While recent studies have identified multifractal scaling in these data, the physical origin of this behavior remains unclear. Existing explanations rely on postulated models (e.g., hierarchical domain models) that reproduce multifractal signatures but fail to explain why a polymer with power-law contact probabilities should exhibit such behavior. Furthermore, the connection between multifractal formalism and fundamental polymer physics scaling laws has remained ambiguous.

Methodology
The authors bridge polymer physics and multifractal formalism by starting from a single empirical observable: the contact probability $P(s)$ , which decays as a power law $P(s) \sim s^{-\gamma}$ across genomic distances.

Analytical Derivation: Instead of assuming specific polymer conformations, the authors derive the mass exponent $\tau(q)$ analytically from $P(s)$ . They model the Hi-C map as a two-dimensional geometric object and decompose the generalized partition function $Z(q, l)$ into diagonal (intra-segment) and off-diagonal (inter-segment) contributions.
Model Variations:
- Single-Exponent Model: Assumes a single power-law decay $P(s) \sim s^{-\gamma}$ .
- Double-Exponent Model: Extends the analysis to $P(s)$ with two distinct exponents ( $\gamma$ for short distances and $\eta$ for long distances) separated by a transition scale $s_T$ , reflecting empirical observations of TADs and loop extrusion.
Validation: The theoretical predictions are validated against:
- Synthetic contact maps generated from the defined power-law models.
- Empirical Hi-C data from human (GM12878) and mouse (CH12.LX) cells.
- Hi-C data from six diverse Australian species (mammals, birds, reptiles, fish, insects, and plants).
Robustness Testing: The analysis includes the addition of additive white noise to empirical maps to test the stability of the multifractal signal against experimental perturbations.

Key Results

Origin of Multifractality: The study demonstrates that multifractal behavior in Hi-C maps is a direct geometric consequence of the power-law contact probability $P(s)$ and the competition between intra-segment (diagonal) and inter-segment (off-diagonal) contacts. No hierarchical cascade model is required to generate this behavior.
Analytical Form of $\tau(q)$ :
- For a single exponent $\gamma$ , the mass exponent $\tau(q)$ exhibits two asymptotic linear regimes separated by a smooth crossover.
- Small $q$ (dominated by off-diagonal/inter-segment contacts): The slope is $\tau(q) \approx 2(q-1)$ .
- Large $q$ (dominated by diagonal/intra-segment contacts):
  - If $\gamma \ge 1$ , the slope is $\tau(q) \approx q-1$ .
  - If $\gamma < 1$ , the slope is $\tau(q) \approx q(2-\gamma) - 1$ .
- Crucially, the large- $q$ slope for $\gamma < 1$ maps directly to the bulk contact exponent $\beta_b$ (where $\beta_b = \gamma$ in the context of the contact probability, leading to a slope of $2-\gamma$ ).
Double-Exponent Behavior: When $P(s)$ $P (s)$ has two regimes ( $\gamma$ $γ$ and $\eta$ $η$ ), the observed slope of $\tau(q)$ $τ (q)$ at large $q$ $q$ depends on the box size $l$ $l$ relative to the transition scale $s_T$ $s_{T}$ .
- For small $l$ ( $l < s_T$ ), the slope reflects the short-range exponent $\gamma$ .
- For large $l$ ( $l > s_T$ ), the slope reflects the long-range exponent $\eta$ , as the coarse-graining integrates over the short-range regime.
Universality and Robustness:
- The multifractal scaling is robust to additive white noise, provided the noise does not completely obscure the power-law decay of $P(s)$ .
- The behavior is consistent across diverse taxonomic groups (human, mouse, seal, swan, dragon, shark, aphid, clover), suggesting it is a universal feature of chromatin organization.
- The nonlinearity of $\tau(q)$ confirms true multifractality rather than bifractality, with a smooth crossover between regimes.

Significance and Claims
The paper claims to provide a physical explanation for multifractal scaling in Hi-C maps, establishing a direct link between the multifractal mass exponent $\tau(q)$ and polymer contact statistics.

Direct Extraction of Exponents: The authors propose a new method to extract the polymer contact exponent $\gamma$ (and $\eta$ ) directly from the slope of $\tau(q)$ without requiring structural modeling or fitting to specific polymer conformations.
Reduction to Two-Point Statistics: Unlike multifractal behaviors in turbulence or financial markets which often require irreducible hierarchical cascade models, the multifractality in Hi-C maps is shown to be "reducible" to two-point contact statistics ( $P(s)$ ) and the two-dimensional geometry of the map.
Physical Interpretation: The results clarify that the observed multifractality arises from the geometric competition between intra- and inter-segment contacts, rather than from spatial heterogeneity in chromatin density alone. The large- $q$ slope specifically encodes the bulk contact exponent, offering a route to probe polymer physics parameters directly from contact maps.