Inheritance Entropy: A Model-Independent Method to Probe the Hereditary Structure of Cell Lineage Trees

이 논문은 관찰된 클론 계통 나무가 비유전적 모델에 비해 유의미하게 낮은 엔트로피 값을 보인다는 점을 근거로, 인간 골수 기질 세포 집락의 이질성이 세포 주기 탈출을 조절하는 유전되지 않는 요인들로부터 기인함을 입증하기 위해 '상속 엔트로피(Inheritance Entropy)'라고 불리는 모델 독립적인 방법을 소개한다.

핵심

가족 계보를 상상해 보세요. 하지만 사람 대신, 배양 접시 안에서 자라나는 줄기세포 군집이 주인공입니다. 세포가 분열할 때마다 두 개의 새로운 가지가 생겨납니다. 보통 이 가지들은 계속 자라나며 나무를 채워 나갑니다. 하지만 가끔, 어떤 세포는 분열을 멈추기로 결정합니다. 이 세포는 '비활성' 상태(잠들거나 은퇴하는 것과 같은 상태)에 들어갑니다. 이렇게 되면 그 가족 나무의 해당 가지 전체가 성장을 멈춥니다.

과학자들이 던진 핵심 질문은 이것이었습니다: 왜 어떤 가지는 성장을 멈추고, 다른 가지는 계속 나아가는가?

이것은 단순히 무작위적인 우연일까요? 아니면 멈추라는 명령을 전달하는 숨겨진 '가족의 비밀'이 부모에서 자식으로 물려 내려오는 것일까요?

문제: 기계 속의 "유령"

이 현상을 연구하는 까다로운 점은, 세포가 분열을 멈추면 그 세포는 미래의 나무 구조에서 사실상 사라져 버린다는 것입니다. 자식이 없는 세포의 후손을 연구할 수는 없습니다. 이는 마치 존재하지 않는 잎사귀들을 보고 특정 나무 가지가 왜 죽었는지 알아내려는 것과 같습니다.

보통 돌연변이가 발생하면 그것은 파동처럼 퍼져 나갑니다. 하지만 '정지'의 경우, 그 파동은 오히려 '구멍'이 됩니다. 과학자들은 만약 정지 결정이 순수하게 무작위적이라면, 나무 곳곳에 구멍이 빗방창에 맺힌 물방울처럼 흩어져 있을 것이라고 생각했습니다. 하지만 만약 이 결정이 유전되는 것(가족의 특질처럼 물려 내려오는 것)이라면, 구멍들은 특정 가지들에 뭉쳐서 나타날 것입니다.

해결책: "유전 엔트로피(Inheritance Entropy)" 측정하기

과학자들은 유전 엔트로피라는 새로운 수학적 도구를 발명했습니다.

가족 나무를 지저집한 방이라고 생각해 보세요.

• 높은 엔트로피 (무질서/무작위): 만약 '정지된' 세포들이 나무 전체에 무작위로 흩어져 있다면, 방은 혼란스럽고 예측 불가능해 보일 것입니다. 이는 정지가 그저 운이 나빴거나 무작위적인 노이즈였음을 시사합니다.
• 낮은 엔트로피 (질서/유전): 만약 '정지된' 세포들이 단 하나 혹은 두 개의 특정 가지에 모여 있다면, 방은 이상하리만큼 조직화되어 보일 것입니다. 이는 특정 가지에 성장을 멈추라고 명령하는 어떤 '규칙'이 전달되었음을 시사합니다.

과학자들은 이 "무질서도"(엔트로피)를 32개의 서로 다른 인간 줄기세포 군집에 대해 계산했습니다.

발견: 그것은 무작위가 아니다

그들은 실제 세포 군집을 '뒤섞은(scrambled)' 버전과 비교했습니다. 실제 가족 나무를 가져와서 '정지된' 세포들을 모두 잘라낸 뒤, 이를 다른 가지들에 무작위로 다시 붙였다고 상상해 보세요. 이렇게 하면 유전이 불가능한 '가짜' 나무가 만들어집니다.

결과: 실제 군집의 75%에서 "무질서도"가 뒤섞인 가짜 나무보다 유의미하게 낮게 나타났습니다.

• 번역하자면: '정지된' 세포들은 무작위로 흩어져 있지 않았습니다. 그들은 특정 가지에 뭉쳐 있었습니다. 이는 정지 결정이 유전된다는 것을 증명합니다. 이것은 단순한 무작위적 우연이 아니라, 부모 세포에서 자식 세포로 물려 내려오는 하나의 형질입니다.

"지연(Lag)" 효과: 지연된 반응

과학자들은 또한 이것이 언제 일어나는지에 대해서도 흥ável한 사실을 발견했습니다. 그들은 "돌연변이 지연"을 발견했습니다.

어떤 세포가 '정지 신호'(돌연변이)를 받았다고 가정해 봅시다. 세포는 즉시 멈추지 않습니다. 세포는 몇 세대 동안 계속 분열하며 그 신호를 자식들에게 전달합니다. 몇 세대가 지난 후에야 세포들은 실제로 브레이크를 밟고 멈춥니다.

• 발견 내용: '정지 신호'가 생성된 순간과 세포가 실제로 분열을 멈추는 순간 사이에는 평균적으로 2세대의 지연이 존재합니다.
• 중요한 이유: 이 지연은 매우 중요합니다. 만약 세포가 즉각적으로 멈췄다면, 과학자들은 신호의 근원을 추적할 수 없었을 것입니다. 지연이 있었기 때문에, 그들은 '비어 있는' 가지들을 살펴보고 신호가 처음 생성된 정확한 시점까지 경로를 거슬러 올라갈 수 있었습니다.

"가족의 비밀" (후성유전학)

논문은 이 유전이 **후성유전적(epigenetic)**일 가능성이 높다고 결 결론짓습니다.

• 비유: DNA를 컴퓨터의 하드웨어(변하지 않는 코드)라고 한다면, 후성유전학은 소프트웨어 설정(볼륨, 밝기, 또는 '방해 금지 모드'와 같은 것)입니다.
• 세포들은 DNA 코드를 바꾸는 것이 아닙니다. 대신, 그들은 언제 일을 멈출지를 알려주는 "설정"을 물려받는 것입니다. 이것이 왜 동일한 유형의 세포에서 시작했음에도 불구하고 두 세포 군집이 서로 매우 다르게 보일 수 있는지를 설명해 줍니다.

요약

1. 방법: 그들은 세포의 죽음이 무작위인지 혹은 조직화되어 있는지 확인하기 위해 "무질서도"(엔트로피)를 사용했습니다.
2. 발견: 대부분의 군집에서 죽음은 조직화되어 있었으며, 이는 '정지' 신호가 유전된다는 것을 증명합니다.
3. 타이밍: 신호가 보내진 것과 세포가 실제로 멈추는 것 사이에는 2세대의 지연이 있습니다.
4. 원인: 이것은 유전적 돌연변이(하드웨어 변화)보다는 후성유전적 요인(소프트웨어 설정)에 의한 것입니다.

이 논문은 이것이 질병을 즉각 치료하거나 이식을 개선할 것이라고 주장하는 것이 아니라, 줄기세포 군집이 성장을 결정하고 멈추는 패턴을 결정하는 강력하고 유전적인 구조를 가지고 있음을 증명하는 새로운, 모델 프리(model-free) 방식의 방법을 제공한다는 점을 명시하고 있습니다.

기술 요약

기술 요약: 세포 계보 트리에서의 상속 엔트로피 (Inheritance Entropy)

1. 문제 정의

세포 계보 트리 내에서 유전적 전달(hereditary transmission)의 영향을 정량화하는 것은 생물학의 근본적인 과제로 남아 있다. 단일 세포 프로파일링은 세포 상태의 특성을 규명할 수 있게 해주지만, 이를 고전적인 계보 추적과 결합하여 유전적 파급 효과를 추적하는 것은 여전히 어렵다. 특히 비활성(세림, 분화 또는 휴지기로 인해 세포 분열이 멈춘 상태)과 관련하여 구체적인 문제가 발생한다. 활발한 증식의 경우 유전적 돌연변이가 자손에게 표현형을 증폭시키지만, 비활성은 자기 제한적이다. 즉, 일단 세포가 비활성화되면 해당 분지는 종료되며, 이는 하류(downstream)의 정보를 소멸시킨다. 이로 인해 비활성의 출현이 무작위적인 확률적 사건에 의한 것인지, 아니면 세포 주기 이탈을 조절하는 유전적(후성유전적) 요인에 의한 것인지 판단하기 어렵다. 또한, 줄기세포 집단에서 집단 간 이질성이 치료 잠재력에 영향을 미치는 상황에서, 유전적 돌연변이와 후성유전적 상속을 구분하는 것은 매우 중요하다.

2. 방법론

저자들은 계보 트리의 위상(topology)에 적용된 **샤논 엔트로피(Shannon entropy)**를 기반으로 한 모델 독립적 방법을 소개한다. 이 방법은 특히 비활성 세포의 분포에 초점을 맞춘다.

A. 데이터 수집

본 연구는 32개의 인간 골수 기질 세포(BMSC) 클론 집단의 고해상도 단일 세포 계보 추적 데이터를 활용한다. 실험은 모든 복제 세포가 7세대($k_{max}=7$)에 도달할 때까지 위상차 시간 경과 현미경(phase-contrast time-lapse microscopy)을 통해 수행되었다. 세포는 다음과 같이 분류되었다:

• 활성(Active): 분열 관찰됨.
• 비활성(Inactive, G0/사멸): 84시간 동안 분열이 없거나 세포 사멸.
• 리프(Leaves): 탄생 시점만 추적된 $k=7$ 단계의 세포.

B. 비활성 불균형 (Inactivity Imbalance)

비활성의 위상적 영향을 정량화하기 위해, 저자들은 트리 내의 모든 분열(mitosis) $m$에 대해 비활성 불균형($I_m$)을 다음과 같이 정의한다:
$$I_m = |M_{m}^{left} - M_{m}^{right}|$$
여기서 $M$은 왼쪽과 오른쪽 분지에서 누락된 $k_{max}$ 리프의 수를 나타낸다. 특정 노드에서 큰 불균형이 나타난다는 것은 두 딸 분지 사이의 비활성 확률이 갈라졌음을 시사하며, 이는 상류(upstream)에서의 돌연변이 이벤트를 나타낼 가능성이 있다.

C. 상속 엔트로피 (Inheritance Entropy)

핵심 지표인 상속 엔트로피($S$)는 계보 내의 모든 분열에서 발생한 불균형의 정규화된 가중치($w_m$)로부터 계산된다:
$$w_m = \frac{I_m}{\sum_n I_n}, \quad S = -\sum_m w_m \log w_m$$

• 높은 엔트로피: 불균형이 트리 전체에 균질하게 흩어져 있으며, 이는 무작위적이고 비유전적인 변동을 시사한다.
• 낮은 엔트로피: 불균형이 몇몇 노드에 집중되어 있으며, 이는 특정 하위 분지의 비활성 확률을 변화시킨 구체적인 유전적 이벤트(돌연변이)를 나타낸다.

D. 널 앙상블(Null Ensembles) 및 검증

통계적 유의성을 결정하기 위해, 생물학적 트리의 엔트로피를 세 가지 비유전적 널 앙상블과 비교한다:

1. 갈톤-왓슨(Galton-Watson, GW): 모든 세포에 대해 일정한 비활성 확률($q$)을 가짐.
2. 비균질 갈톤-왓슨(Non-Homogeneous GW, NHGW): 세대별로 다른 비활성 확률($q_k$)을 가짐.
3. 무작위 재배치(Random Reshuffling, RR): 선호되는 방법이다. 이는 각 세대 내에서 세포의 위치를 뒤섞되, 세대당 총 비활성 세포 수는 유지하는 방식이다. 이는 계보의 구조적 제약은 유지하면서 모든 모-자(mother-daughter) 유전적 연결을 끊어버린다.

**P-값(P-value)**은 널 앙상블 트리 중 $S \leq S_{biological}$을 만족하는 엔트로피의 비율로 정의된다. 낮은 P-값은 상속의 강력한 증거를 나타낸다.

E. 모델 피팅 (SKH)

저자들은 데이터를 설명하기 위해 확률적 커크우드-홀리데이(Stochastic Kirkwood-Holliday, SKH) 모델을 제안한다. 이 모델에서:

• 야생형(Wild) 세포는 확률 $p$로 "돌연변이" 상태로 변한다.
• 돌연변이 세포(및 그 자손)는 비활성 상태가 될 확률 $q$를 가진다.
• 결정적으로, 모델은 돌연변이 지연(mutation lag, $l$), 즉 돌연변이 이벤트와 비활성 표현형이 나타나는 시점 사이의 세대 수를 포함한다.

3. 주요 기여

1. 상속 엔트로피: 절단된(비활성) 분지의 위상만을 사용하여 계보 트리 내의 유전적 구조를 탐지하는 새로운 모델 불가지론적(model-agnostic) 지표를 제시하였다.
2. 널 모델의 검증: GW 및 NHGW 모델은 비균질한 비유전적 트리와 유전적 트리를 구분하지 못해 위양성(false positive)을 생성하는 반면, 무작위 재배치(RR) 방법은 견고하고 민감하며 모델 독립적임을 입증하는 엄격한 벤치마킹 연구를 수행하였다.
3. 돌연변이 지연 정량화: 비활성 확률이 변화한 특정 분열 지점을 찾아내고, 돌연변이 지연(돌연변이와 비활성 표현형 발현 사이의 지연)을 측정할 수 있는 능력을 제공한다.
4. 확률적 커크우드-홀리데이 모델: BMSC 계보의 유전적 특징(지연 포함)을 성공적으로 포착하는 일반화된 확률 모델을 구축하였다.

4. 결과

• 유전적 구조: 32개의 BMSC 집단 실험 데이터에서, **21개 집단(75%)**이 RR 테스트를 통해 유의미한 유전적 특성($P < 0.05$)을 보였다. 나머지 집단은 비활성 세포 수가 너무 적어 통계적 유의성을 확보하지 못한 것이지, 반드시 유전적 특성이 없음을 의미하는 것은 아니다.
• 엔트로피 분석: 생물학적 트리는 무작위로 재배치된 대조군보다 일관되게 유의미하게 낮은 엔트로피를 보였으며, 이는 비활성 세포의 분포가 무작위가 아니라 유전적 패턴을 따른다는 것을 증사한다.
• 돌연변이 지연: 평균 돌연변이 지연은 약 2세대로 나타났다. 이는 비활성을 유도하는 후성유전적 변화가 실제 표현형이 나타나기 전 상류에서 발생함을 확인해 준다.
• 모델 적합도: 데이터에 피팅된 SKH 모델은 $p \approx 0.1$(돌연변이 확률) 및 $q \approx 0.46$(돌연변이 세포의 비활성 확률)의 매개변수를 산출했다. 이 매개변수로 생성된 합성 트리는 실험적 BMSC 계보의 위상 구조와 매우 유사했다.
• 공간적 vs 위상적: 비활성 세포의 클러스터링이 배양 접시 상의 물리적 근접성(공간적 상관관계)에 의한 인위적 결과가 아니라, 계보의 위상적 특징임을 확인하였다.

5. 의의 및 주장

본 논문은 관찰된 인간 BMSC 집단의 이질성이 순수하게 확률적이거나 유전적인 메커니즘이 아닌, 후성유전적 상속 요인에 의해 조절된다고 주장한다.

• 견고성: 엔트로피 기반 테스트는 계보 트리 내의 차등 성장의 유전적 기원을 평가하는 단순하고 견고하며 모델 불가지론적인 방법을 제공한다.
• 임상적 관련성: 본 연구 결과는 BMSC 이식 시 발생하는 집단 간 이질성(골 재생 치료의 주요 장애물)이 후성유전적 상속 기원을 가지고 있음을 시사한다. 이러한 유전적 메커니즘을 이해하고 제어함으로써 줄기세포 치료의 일관성을 개선할 수 있다.
• 방법론적 발전: "돌연변이 지연"을 설정함으로써, 이 방법은 표현형(비활성)이 스스로 소멸되는 상황에서도 표현형 변화의 기원을 역추적할 수 있게 한다.
• 향후 방향: 저자들은 현재 방법이 비활성에 집중하고 있지만, 계보 위상과 세포 주기 시간 사이의 상관관계는 세포 주기 지속 시간이 후성유전적 상속 메커니즘에 의해 조절될 수 있음을 시사하며, 이에 대한 추가 연구가 필요하다고 언급했다.

저자들은 SKH 모델이 단순화된 모델임을 인정하며, 더 복적인 모델이 지연 시간의 변동성을 더 잘 포착할 수 있다고 밝히며 겸허한 태도를 유지했다. 또한, 이 방법은 변화를 탐지하는 것이므로, 만약 계보 내에 비활성 표현의 변이가 없다면 유전적 메커니즘이 존재하더라도 이를 탐지하지 못할 수 있음을 강조했다.