Multifractal Scaling in Hi-C Maps

기술 요약: Hi-C 맵의 다중 프랙탈 스케일링 (Multifractal Scaling)

문제 정의
게놈의 3차원 구조는 염색체 접촉 지도(Hi-C)와 크로마틴 밀도 현미경을 통해 관찰 가능한 스케일 의존적 구조를 나타낸다. 최근 연구들이 이러한 데이터에서 다중 프랙탈 스케일링(multifractal scaling)을 확인했으나, 이러한 거동의 물리적 기원은 여전히 불분명하다. 기존의 설명들은 다중 프랙탈 신호를 재현하는 가설적 모델(예: 계층적 도메인 모델)에 의존하고 있지만, 왜 거듭제곱 법칙(power-law)을 따르는 접촉 확률을 가진 폴리머가 이러한 거동을 보이는지에 대해서는 설명하지 못한다. 또한, 다중 프랙탈 형식론(multifractal formalism)과 근본적인 폴리머 물리 스케일링 법칙 사이의 연결 고리는 그동안 모호한 상태로 남아 있었다.

방법론
저자들은 단일 경험적 관측량인 접촉 확률 $P(s)$로부터 시작하여 폴리머 물리학과 다중 프록탈 형식론을 연결한다. 여기서 $P(s)$는 게놈 거리 $s$에 따라 $P(s) \sim s^{-\gamma}$로 감소하는 거듭제곱 법칙을 따른다.

1. 분석적 유도: 특정 폴리머 컨포메이션을 가정하는 대신, 저자들은 $P(s)$로부터 질량 지수(mass exponent) $\tau(q)$를 분석적으로 유도한다. 이들은 Hi-C 맵을 2차원 기하학적 객체로 모델링하고, 일반화된 분할 함수 $Z(q, l)$을 대각 성분(intra-segment)과 비대각 성분(inter-segment) 기여분으로 분해한다.
2. 모델 변형:
   • 단일 지수 모델 (Single-Exponent Model): 단일 거듭제곱 법칙 붕괴 $P(s) \sim s^{-\gamma}$를 가정한다.
   • 이중 지수 모델 (Double-Exponent Model): $P(s)$가 두 개의 뚜렷한 지수($s_T$에 의해 구분되는 단거리 $\gamma$와 장거리 $\eta$)를 갖도록 분석을 확장하며, 이는 TAD 및 루프 압출(loop extrusion)의 경험적 관찰을 반영한다.
3. 검증: 이론적 예측은 다음을 통해 검증된다:
   • 정의된 거듭제곱 법칙 모델로부터 생성된 합성 접촉 지도.
   • 인간(GM12878) 및 마우스(CH12.LX) 세포의 경험적 Hi-C 데이터.
   • 6종의 다양한 호주 생물 종(포유류, 조류, 파충류, 어류, 곤충, 식물)의 Hi-C 데이터.
4. 강건성 테스트: 실험적 섭동에 대한 다중 프랙탈 신호의 안정성을 테스트하기 위해 경험적 맵에 가법적 백색 잡음(additive white noise)을 추가하는 분석이 포함된다.

주요 결과

• 다중 프랙탈의 기원: 본 연구는 Hi-C 맵의 다중 프랙탈 거동이 접촉 확률 $P(s)$의 직접적인 기하학적 결과이며, 세그먼트 내부(diagonal) 접촉과 세그먼트 간(off-diagonal) 접촉 사이의 경쟁에서 비롯됨을 입증한다. 이러한 거동을 생성하기 위해 계층적 카스케이드(hierarchical cascade) 모델은 필요하지 않다.
• $\tau(q)$의 분석적 형태:
  • 단일 지수 $\gamma$에 대해, $\tau(q)$는 두 개의 점근적 선형 영역을 가지며 그 사이는 매끄러운 교차(crossover)로 연결된다.
  • 작은 $q$ (비대각/세그먼트 간 접촉이 지배적): 기울기는 $\tau(q) \approx 2(q-1)$이다.
  • 큰 $q$ (대각/세그먼트 내 접촉이 지배적):
    • 만약 $\gamma \ge 1$이면, 기울기는 $\tau(q) \approx q-1$이다.
    • 만약 $\gamma < 1$이면, 기울기는 $\tau(q) \approx q(2-\gamma) - 1$이다.
  • 결정적으로, $\gamma < 1$일 때의 큰 $q$ 기울기는 벌크 접촉 지수 $\beta_b$ (접촉 확률의 맥락에서 $\beta_b = \gamma$인 경우)와 직접 매핑된다.
• 이중 지수 거동: $P(s)$가 두 가지 영역($\gamma$와 $\eta$)을 가질 때, 관찰되는 $\tau(q)$의 기울기는 전이 스케일 $s_T$에 대한 박스 크기 $l$의 상대적 크기에 따라 달라진다.
  • 작은 $l$ ($l < s_T$)의 경우, 기울기는 단거리 지수 $\gamma$를 반영한다.
  • 큰 $l$ ($l > s_T$)의 경우, 거친 격자화(coarse-graining)가 단거리 영역을 통합함에 따라 기울기는 장거리 지수 $\eta$를 반영한다.
• 보편성 및 강건성:
  • 다중 프랙탈 스케일링은 $P(s)$의 거듭제곱 법칙 붕괴를 완전히 가리지 않는 한, 가법적 백색 잡음에 대해 강건하다.
  • 이 거동은 다양한 분류군(인간, 마우스, 물개, 백조, 드래곤, 상어, 진딧물, 클로버)에 걸쳐 일관되게 나타나며, 이는 크로마틴 구조의 보편적 특징임을 시사한다.
  • $\tau(q)$의 비선형성은 단순한 바이프랙탈(bifractality)이 아닌 진정한 다중 프랙탈임을 확인시켜 주며, 두 영역 사이의 매끄러운 교차를 보여준다.

의의 및 주장
본 논문은 Hi-C 맵의 다중 프랙탈 스케일링에 대한 물리적 설명을 제공하며, 다중 프랙탈 질량 지수 $\tau(q)$와 폴리머 접촉 통계 사이의 직접적인 연결 고리를 확립한다고 주장한다.

• 지수의 직접 추출: 저자들은 구조적 모델링이나 특정 폴리머 컨포메이션에 대한 피팅 없이도 $\tau(q)$의 기울기로부터 폴리머 접촉 지수 $\gamma$ (및 $\eta$)를 직접 추출할 수 있는 새로운 방법을 제안한다.
• 이점 통계로의 환원: 난류나 금융 시장의 다중 프랙탈과 달리, Hi-C 맵의 다중 프랙탈은 2점 접촉 통계($P(s)$)와 2차원 맵의 기하학으로 "환원 가능한" 것으로 나타났다.
• 물리적 해석: 연구 결과는 관찰된 다중 프랙탈이 단순히 크로마틴 밀도의 공간적 불균질성 때문이 아니라, 세그먼트 내부 및 세그서 간 접촉의 기하학적 경쟁에서 발생함을 명확히 한다. 특히 큰 $q$에서의 기울기는 벌크 접촉 지수를 인코딩하며, 이는 접촉 맵으로부터 폴리머 물리 파라미터를 직접 탐사할 수 있는 경로를 제공한다.