Longitudinal multi-omics characterization of the malignant evolution in multirelapsing glioblastoma

기술 요약: 다회 재발성 교모세포종의 종양 진화에 대한 종단적 멀티오믹스 특성 규명

문제 제기
교모세포종(GBM)은 2005년 이후 표준 치료법이 변하지 않은 채 여전히 극복하기 어려운 과제로 남아 있다. 치료 실패를 이해하고 신약 개발을 진전시키는 데 있어 주요한 장벽은 종양의 치료 과정 전반에 걸친 종양 조직의 종단적 분석이 부족하다는 점이다. 질병의 진행에 대한 포괄적인 분자 프로파일링을 시도했던 이전의 노력들은 단 한 번의 재발만을 샘플링하거나, 유전적 진화, 세포 가소성, 그리고 임상적 재발 사이의 역동적인 상호작용을 포착하지 못했다는 한계가 있었다. 더욱이, 교모세포종 환자의 짧은 생존 기간은 특히 뇌외 전이 및 말기 진행과 관련하여 전체적인 종양 궤적을 지도화하는 데 필요한 연속적인 샘플 수집을 어렵게 만든다.

연구 방법론
본 연구는 리-프라우메니 증후군(Li-Fraumeni syndrome, 생식선 TP53 돌연변이)을 가진 성인 남성 환자 1명의 사례를 대상으로 한 다차원적 특성 규명을 제시한다. 이 환자는 초기 진단부터 말기 뇌외 전이(림프절)에 이르기까지 31개월 동안 11회의 신경외과적 절제술을 받았다.

연구에는 다음과 같은 종단적 멀티오믹스 접근 방식이 사용되었다:

• 유전체 분석 (Genomic Analysis): 11개의 모든 종양 샘플(T1–T11)에 대해 벌크 DNA 시퀀싱을 수행하여 체세포 돌연변이, 복제수 변이(CNVs), 그리고 이형접합성 상실(LOH)을 식별하였다. 고신뢰도 변이를 사용하여 계통 발생 트리(Phylogenetic trees)를 재구성함으로써 종양 진화와 클론 분포를 추론하였다.
• 전사체 분석 (Transcriptomic Analysis): 11개 샘ール 모두에 대해 벌크 RNA 시퀀싱을 수행하여 유전자 발현 역동성을 평가하였다. 7개의 서로 다른 해부학적 위치를 포함하는 샘플 T4에서 T10까지는 단일 세포 RNA 시퀀싱(scRNA-seq)을 수행하였으며, 41,000개 이상의 고품질 세포를 확보하였다.
• 컴퓨터 통합 분석 (Computational Integration):
  • 디콘볼루션 (Deconvolution): 벌크 RNA-seq 데이터를 GBMDeconvolutR로 분석하여 세포 비율(종양, 면역, 기질 세포)을 추론하였다.
  • 세포 상태 매핑 (Cell State Mapping): scRNA-seq 데이터를 통합하고 Seurat을 사용하여 클러스터링함으로써 표현형 세포 상태(Mesenchymal [MES], Astrocyte-like [AC], Oligodendrocyte Progenitor-like [OPC], Neural Progenitor-like [NPC])를 식별하였다.
  • 궤적 분석 (Trajectory Analysis): Monocle을 사용하여 세포를 줄기세포 유사 상태에서 분화 상태로 순서를 매기는 의사시간(pseudotime) 분석을 수행하였다. RNA 속도(scVelo)를 활용하여 미래의 세포 상태와 전이 방향을 예측하였다.
  • 치료 반응 모델링 (Therapy Response Modeling): "Quadriga plot"을 개발하여 4가지 주요 세포 상태의 풍부도와 정밀 의료 타겟(예: EGFR, JAK1/2, mTOR)의 유전자 시그니처 강도를 특정 치료에 따른 반응과 함께 시각화하였다.

주요 결과

1. 유전적 진화 및 클론 구조:

   • 11개의 모든 종양 병변은 NF1, RB1, GPS2, APOB의 초기 돌연변이와 생식선 TP53 돌연변이에 의해 정의되는 공통된 유전적 조상을 공유했다.
   • 본 연구는 13개의 뚜렷한 종양 클론을 식별하였다. 초기 돌연변이는 보편적이었으나, 후기 단계의 진행은 TSC2, PBRM1, CD22, FANCI의 돌연변이 획득에 의해 주도되었다.
   • 클론들은 평행 진화를 보였다. 예를 들어, RB1과 TSC2 돌연변이가 독립적인 클론에서 나타났는데, 이는 적응도(fitness)의 우위를 시사한다.
   • 클론 $\epsilon$은 초기에 출현하여 말기 뇌외 전이(비강 및 림프절)를 포함한 후기 샘플들을 지배하였으며, 이는 광범위한 확산이 가능한 지속적이고 오래 살아남는 하위 집단의 존재를 나타낸다.

2. 전사체 역동성 및 세포 상태 가소성:

   • 치료 유도 변화: 종양은 정밀 의료에 반응하여 표현형 세포 상태의 가역적인 변화를 보였다. 예를 들어, EGFR 억제제인 오시머티닙(Osimertinib, Tagrisso) 처리는 MES 유사 시그니처의 감소와 보상적인 AC 유사 세포의 농축을 초래했으며, 이 세포들은 이후 EGFR을 상향 조절하였다. 마찬가지로, JAK/STAT 억제제(Jakafi)와 mTOR 억제제(Everolimus)는 MES 유사 인구 집단을 일시적으로 감소시켰고, 이는 AC 유사 세포로 대체되었다.
   • 면역 맥락: 초기 MES 유사 시그니처는 높은 골수구 세포 침윤과 상관관 $\��{e}$였으며, 이는 벌크 데이터에서의 MES 시그니처가 종양 자체의 상태보다는 면역 반응을 반영하는 경우가 많음을 시사한다. 후기 단계에서는 면역 소진과 pro-neural 시그니처로의 전환이 관찰되었다.
   • 전이 궤적: 질병이 말기 뇌외 전이로 진행됨에 따라, pro-neural 및 neuron-like GBM 세포 시그니처의 강력하고 일관된 농축이 나타났다. RNA 속도 분석은 세포가 MES/AC 유사 상태에서 침습적 특징을 가진 pro-neural 상태로 향하는 궤적을 따르고 있음을 보여주었다.

3. 뇌외 전이:

   • 본 연구는 뇌외 전이(비강 및 목 림프절)가 일차 뇌종양과 동일한 근원적 유전적 시그니처를 유지하면서도, 구별되는 pro-neural 및 neuronal 전사체 특징의 상향 조절을 특징으로 한다는 것을 확인하였다. 이는 뉴런 모방(neuronal mimicry)이 뇌내 전이와 뇌외 확산을 모두 주도하는 핵심 프로그램임을 시사한다.

주요 기여

• 전례 없는 종단적 규모: 본 연구는 11회의 수술 샘플을 통해 31개월 동안의 GBM 진화에 대한 고해상도 뷰를 제공하며, 일차 종양에서 말기 뇌외 전이로의 전이를 포착하였다.
• 멀티오믹스 상관관계: 특정 유전적 사건(예: TSC2, FANCI 돌연변이)을 전사체 변화 및 임상 결과와 성공적으로 연결함으로써, 유전적 진화가 어떻게 표현형 가소성을 유도하는지 입증하였다.
• 치료 반응 매핑: 단일 세포 데이터를 임상 치료 이력과 통합함으로써, 특정 표적 치료(Nivolumab, Palbociclib, Osimertinib, Everolimus, Jakafi)가 어떻게 특정적이고 종종 가역적인 세포 적응(예: MES-to-AC 전이)을 유도하는지 규명하였다.
• 자원 생성: 저자들은 통합된 단일 세포 데이터를 탐색할 수 있는 오픈 액세스 Shiny 앱인 "GBMaps"를 개발하였으며, 유전적, 분자적, 세포적, 임상적 진행을 연결하는 포괄적인 데이터셋을 제공하였다.

의의
본 논문은 GBM의 치료 반응과 실패의 복잡성을 이해하기 위한 "로드맵"을 제공한다고 주장한다. 이 연구의 주요 의의는 GBM 세포가 치료적 압박에 반응하여 표현형 상태(MES, AC, pro-neural) 사이를 이동하며 놀라운 적응력을 보인다는 점을 입증한 데 있다. 결정적으로, 본 연구는 특정 클론을 일시적으로 억제하더라도 종양이 대안 경로 또는 상태 전이를 통해 적응한다는 것을 보여준다.

저자들은 pro-neural 및 neuronal 특징의 농축이 뇌내 전이와 뇌외 전이를 모두 이끄는 근본적인 동력이라고 결론짓는다. 이러한 발견은 치료 패러다임의 재평가를 촉구하며, 재발을 방지하기 위해서는 세포 상태의 가소성과 GBM 세포의 뉴런 모방(neuronal mimicry)을 표적으로 삼는 것이 필요할 수 있음을 시사한다. 본 연구는 저항 기전을 해독하기 위해 연속적인 조직 샘플링이 필수적임을 강조하며, 정적인 단일 시점 분석으로는 GBM 진화의 역동적인 본질을 포착하기에 불충분하다고 주장한다.