Spatial mitochondrial lineage tracing uncovers a premetastatic niche and microenvironment programmed fate switching in osteosarcoma

기술 요약: 공간적 미토콘드리아 계통 추적을 통한 골육종의 전이 전 니치(Premetastatic Niche) 및 미세환경 유도형 운명 전환 규명

문제 제기
골육종의 진행은 복잡한 시공간적 역동성을 특징으로 하지만, 종양 세포의 운명 결정을 지배하는 계통 관계는 여전히 제대로 이해되지 않고 있다. 단일 세포 RNA 시퀀싱(scRNA-seq)은 세포 내 이질성을 밝혀냈으나, 이는 세포를 고유한 조직 구조로부터 분리하여 세포 근린(cellular neighborhoods)을 이해하는 데 필수적인 공간 좌표를 소실시킨다. 반대로 공간 전사체학(ST)은 구조를 보존하지만 정적인 스냅샷만을 포착하므로, 인접한 세포들이 공통된 클론 기원을 공유하는지 아니면 수렴 진화를 통해 유사한 상태로 수렴했는지를 구별할 수 없다. 또한, 기존의 계통 추적 방법들은 대개 형질전환 조작(예: CRISPR 바코딩)에 의존하므로, 일차 인간 임상 검체에는 적용하기 어렵다. 따라서 내인성 추적자를 사용하여 공간, 시간, 계통 역사를 통합함으로써, 전이 촉진 상태가 일시적인 미세환경 가소성에서 비롯된 것인지 아니면 결정론적인 클론 잠금(clonal locking)에 의한 것인지를 결정할 수 있는 기술적 필요성이 존재한다.

방법론
저자들은 처리되지 않은 상태의 말기 비전이성 인간 골육종 검체에서 세 가지 분석 차원—공간, 시간, 계통—을 통합하였다:

1. 단일 세포 및 공간 전사체학: 11,203개의 세포에 대해 scRNA-seq를 수행하여 세포 집단을 정의하고, 악성 세포와 비악성 세포를 구분하기 위해 복제수 변이(CNV)를 추론하였다. 동일한 샘플에 대해 공간 전사체학(10x Visium)을 수행하여 유전자 발현을 조직 구조 상에 매핑하였다.
2. 의사시간 재구성(Pseudotemporal Reconstruction): 저자들은 베이지안 RNA 속도 모델인 cell2fate를 적용하여 전사 역동성과 유전자 모듈을 추론하고, 종양 상태와 기질 리모델링의 시간적 진행 과정을 재구성하였다.
3. 내인성 계통 추적: 유전적 조작 없이 클론의 조상을 해결하기 위해, 저자들은 자연적으로 발생하는 체세포 미토콘드리아 DNA(mtDNA) 변이를 활용하였다. 이들은 MASTER 프로토콜을 사용하여 Visium cDNA 라이브러리로부터 미토콘드리아 전사체를 농축하였다.
4. 공간 및 계통 분석:
   • 공간적 일관성: 미토콘드리아 변이 대립유전자 빈도(VAF)에 Moran's I 기반의 공간 자기상관 필터를 적용하여, 클론의 공간적 일관성을 확인하고 깊이(depth)에 의한 인위적 오류를 제거하였다.
   • 계통 재구성: Wright-Fisher 표류 모델과 은닉 마르코프 트리(WF-HMT)를 결합한 MitoDrift 프레임워크를 사용하여 계통수를 재구성하고 확률적 미토콘드리아 분리를 모델링하였다.
   • 운명 전환 모델링: 이산 상태 연속 시간 마르코프 모델(ARD)을 적합시켜 세포 상태(COL3A1⁺, ALPL⁺, THY1⁺ 등) 간의 방향성 전환율을 추정하였다.
   • 상호작용 분석: 세포 간 통신을 분석하기 위해 CellChat과 SpaCET를 사용하여 공간적 맥락 내에서의 신호 전달 허브를 식별하였다.

주요 결과

• 종양 이질성 및 CNV 프로파일링: scRNA-seq는 12개의 뚜렷한 집단을 식별하였으며, 종양 세포(66.5%)가 경관을 지배하였다. 악성 하위 집단에는 ALPL⁺(골모세포), THY1⁺(중배엽/줄기세포 유사), COL3A1⁺(전구세포/매트릭스 리모델링) 세포가 포함되었다. CNV 프로파일링은 이 하위 집단들의 악성 기원을 확인하였으며, 특히 THY1⁺ 세포는 염색체 16의 독특한 획득을 나타냈다.
• 단계적 운명 분기: 의사시간 재구성은 COL3A1⁺ 전구세포가 ALPL⁺ 골모세포 계통과 THY1⁺ 중배엽 계통이라는 두 개의 뚜렷한 계통으로 분기되는 단계적 폭포 현상을 보여주었다. 이 분기는 골수 세포 유도 리모델링 모듈의 활동 정점과 시간적으로 일치한다.
• 분화된 전사 프로그램:
  • ALPL⁺ 계통으로의 전환은 EMT 프로그램과 산화적 인산화로의 대사 전환에 의해 유도된다.
  • THY1⁺ 계통으로의 전환은 NF-κB 매개 염증 재프로그래밍, 항원 처리 저하, 그리고 ALPL⁺ 경로와는 구별되는 대사 재편성을 특징으로 한다.
• THY1_Endo 니치의 발견: 공간 전사체학은 THY1⁺ 종양 세포가 PLVAP⁺ 기능 장애 내피세포와 공존하는 특정 전이 촉진 니치를 식별하였다. 이 THY1_Endo 니치는 PTN(pleiotrophin) 및 NOTCH 신호 전달이 풍부한 공간적으로 제한된 신호 허브를 형성한다. 세포 간 상호작용 분석 결과, THY1⁺ 세포는 다른 종양 아형에서는 관찰되지 않는 내피세포 및 특정 대식세포 아형과의 강력한 신호 상호작용을 수행함이 확인되었다.
• 계통-상태 결합 및 방향성 운명 잠금:
  • 공간적 일관성: 미토콘드리아 변이는 유의미한 공간 자기상관(Z = 10.07, p ≈ 4×10⁻²⁴)을 보였으며, 이는 mtDNA VAF가 Visium 해상도에서 충분한 클론 정보를 유지함을 입증한다.
  • 비가역성: WF-HMT 모델을 이용한 계통 재구성은 전환율의 현저한 비대칭성을 드러냈다. COL3A1⁺에서 THY1⁺로의 전환율은 724.08인 반면, 역방향 전환율은 4.09에 불과했다(약 177배의 비대칭성). 반면, COL3A1⁺에서 ALPL⁺로의 전환율은 0.47로 낮았다.
  • 클론 독립성: 낮은 Blomberg's K 값(<0.1)은 THY1⁺ 및 ALPL⁺ 상태가 서로 다른 COL3A1⁺ 전구세포로부터 여러 차례 독립적으로 발생했음을 나타내지만, 일단 확립되면 THY1⁺ 운명은 본질적으로 비가역적이다.
  • 돌연변이 시그니처: 고순도 THY1⁺ 및 ALPL⁺ 클론에서 클레이드 특이적 미토콘드리아 돌연변이가 식별되었으며, 이는 표현형 잠금이 일어나기 전 돌연변이 프라이밍이 선행됨을 시사한다.

의의 및 주장
본 논문은 골육종 진행에 대한 "시공간-계통(time-space-lineage)" 전달 모델을 제안한다. 저자들은 전이 능력을 갖춘 THY1⁺ 중배엽 운명이 세포 자율적 프로그램에 의해서만 결정되는 것이 아니라, 기능 장애 내피세포와의 직접적인 공존을 통한 엄격한 **공간적 라이선싱(spatial licensing)**을 필요로 하며, 이를 통해 THY1_Endo 니치를 형성한다고 주장한다.

주요 주장:

1. 방법론적 검증: 본 연구는 체세포 미토콘드리아 변이가 형질전환 조작 없이도 일차 인간 종양에서 강력한 내인성 계통 추적자로서 클론의 역사와 운명 결정을 해결할 수 있음을 입증한다.
2. 기전적 통찰: 본 연구는 THY1⁺ 상태가 COL3A1⁺ 전구세포로부터의 방향성 전환에 의해 유도되는, 일시적인 가소적 상태가 아닌 비가역적이고 클론적으로 잠긴 종말 분화임을 확립한다.
3. 니치 의존성: THY1⁺ 세포의 전이 능력은 THY1_Endo 니치(PTN/NOTCH)가 제공하는 공간적으로 제한된 신호 환경에 부수적이다.
4. 임상적 함의: THY1_Endo 니치는 실질적인 공간적 타겟으로 제안된다. 저자들은 전이 폭포를 차단하기 위해 계통 결정(비가역적 운명)과 니치 지원(공간적 신호)을 모두 저해하는 전략을 목표로 해야 한다고 제안한다.

결론적으로, 본 연구는 RNA 속도, 공간 전사체학, 미토콘드리아 계통 추적을 통합함으로써 단일 세포 상태, 공간적 조직, 시간적 역동성, 그리고 클론 조상을 연결하는 다차원적 프레임워크를 제공하며, 이는 암의 진화를 이해하고 잠재적으로 차단하는 데 필수적이다.