Microenvironmental Perfusion Gradients Shape the CpG Methylome to Drive Glioblastoma Intratumoral Heterogeneity

기술 요약: 미세환경 관류 구배가 CpG 메틸로옴을 형성하여 교모세포종의 종양 내 이질성을 유도함

문제 제기
교모세로종(GBM)은 극심한 종양 내 이질성(ITH)을 특징으로 하며, 이는 치료 저항성과 불량한 환자 생존율의 주요 동인이다. 유전적 및 후성유전적 변이가 이러한 ITH의 기여 요인으로 알려져 있으나, 역동적인 종양 미세환경(TME)과 안정적인 후성유전적 상태를 연결하는 상위 계층적 드라이버는 아직 해결되지 않은 과제로 남아 있다. 특히, 차별적인 조직 관류(perfusion)—이로 인해 형성되는 뚜렷한 혈관 주변, 중간, 저산소 니치(niche)—와 그에 따른 결과로서의 게놈 전반의 CpG 메틸화 지형 사이의 관계는 포괄적으로 매핑되지 않았다. 기존 연구들은 DNA 메틸화에 미치는 TME의 영향을 강조해 왔으나, 관류 구배가 어떻게 니치 특이적 후성유전적 서명을 각인하고 발암성 전사를 유도하는지에 대한 공간적으로 분해된 이해는 여전히 미지의 영역이다.

방법론
본 연구는 관류 구배와 후성유전적 및 전사 상태를 상관 분석하기 위해 다중 오믹스 및 공간 분해 접근법을 채택하였다:

1. 관류 기반 형광 염료 세포 표지(PFDLC): 저산소/혈관 분포를 모사하기 위해, 저자들은 NOD-SCID 마우스의 U87MG-GFP 제노그래프트 모델을 사용하여 꼬리 정맥 주사를 통해 Hoechst 33342 염료를 투여하였다. 이 염료는 기능적 혈관과의 근접성에 따라 세포를 표지한다.
2. 세포 분리(Cell Sorting): 유세포 분석법을 사용하여 염료 흡수 강도에 따라 세 가지 뚜렷한 하위 집단으로 구분된, 기질이 제거된 GFP+ 암세포를 분리하였다:
   • HH (Hoechst-High): 혈관 주변 니치 (혈류가 원활함).
   • HM (Hoechst-Medium): 중간 니치.
   • HL (Hoechst-Low): 저산소 니치 (혈류가 원활하지 않음).
3. 다중 오믹스 프로파일링:
   • 메틸로옴(Methylome): 세 분획 간의 차별적 메틸화 CpG 부위를 식별하기 위해(Δβ 임계값으로 정의) 고밀도 Illumina 850k EPIC 마이크로어레이 분석을 수행하였다.
   • 전사체(Transcriptome): 동일한 분획에서 차별 발현 유전자(DEG)를 식별하기 위해 RNA 시퀀싱(Cell-Seq2 프로토콜)을 수행하였다.
   • 통합(Integration): 메틸로옴과 전사체 데이터를 통합하여 저메틸화가 상향 조절(HoU)과 연관되고, 과메틸화가 하향 조절(HeD)과 연관되는 유전자를 식별하였으며, 특히 프로모터 영역과 CpG 아일랜드에 집중하였다.
4. 검증 및 공간 매핑:
   • 생체 내/임상(In Vivo/Clinical): 제노그래프트 및 인간 GBM 환자 조직(저산소 마커인 GLUT1과 공동 염색)에 대해 면역조직화학(IHC)을 수행하여 단백질 발현을 검증하였다.
   • 단일 세포 RNA 시퀀싱(scRNA-seq): 세 개의 공공 GBM 데이터셋을 통합 분석하고, 세포 저산소 예측 프레임워크(CHPF) 및 k-평균 클러스터링을 사용하여 연속적인 저산소 구배에 따른 유전자 발현을 매핑하였다.
   • 공간 전사체학(Spatial Transcriptomics): 인간 IDH-와일드타입 GBM 조직에 대한 Visium 기반 공간 전사체학을 사용하여, 이중 저산소 및 혈관 근접도 점수를 바탕으로 특정 조직학적 구역(예: 팔리세이딩 괴사, 저산소 핵심부)에 대한 유전자 발현을 매핑하였다.
   • 임상 상관관계: TC-GA 및 Ivy GAP 데이터셋의 발현 데이터를 분석하여 GBM 아형, G-CIMP 상태 및 환자 생존율과의 연관성을 조사하였다.

주요 결과

• 관류 주도적 후성유전적 구역화: 게놈 전반의 분석 결과, 가장 광범위한 CpG 메틸화 변이는 미세환경의 양 극단(HH vs. HL)에서 발생함을 밝혔다. 두 극단 사이에는 98,355개의 차별적 메틸화 CpG가 존재하였으며, 비무작위적인 염색체 핫스팟(특히 염색체 1번 및 2번)이 관찰되었다.
• 기능적 타겟팅: 저산소(HL) 니치에서는 프로모터 영역과 CpG 아일랜드에서 과메틸화가 우세하게 나타났다(각각 60.8% 및 63.9%). 경로 분석 결과, 과메틸화된 유전자들은 TGF-β, Rap-1, PD-L1/PD-1 면역 체크포인트 경로에 풍부하였고, 저메틸화된 유전자들은 MAPK, Wnt, Hedgehog 경로에 풍부하였다.
• 핵심 4-유전자 서명 식별: 메틸로옴과 전사체 데이터를 통합하여 저산소 니치에서 저메틸화 및 상향 조절되는 핵심 유전자 세트인 ABCA1, GBE1, GPNMB, CXCL14를 식별하였다.
  • 이 유전자들은 Ivy GAP 데이터셋에서 "팔리세이딩 괴사(palisading necrosis)" 및 "주변 괴사 구역(perinecrotic zones)"에서 가장 높은 발현을 보였다.
  • 이 서명의 높은 발현은 공격적인 중간엽(mesenchymal) 아형, 비-G-CIMP 표현형, 그리고 TCGA 데이터에서의 유의미한 환자 생生存 감소와 상관관계가 있었다.
• 다양한 양식 간의 검증:
  • IHC: 마우스 제노그래프트와 인간 환자 조직 모두에서 저산소(GLUT1-high) 구역 내 ABCA1, GBE1, GPNMB의 단백질 수준이 상승했음을 확인하였다. CXCL14는 상대적으로 덜 일관된 공간적 제한성을 보였다.
  • scRNA-seq: 타겟 유전자들(특히 ABCA1, GBE1, GPNMB)은 증가하는 저산소 점수에 따라 단계적인 상향 조절을 보였으며, 가장 심각한 저산소 하위 클러스터(H2 및 H3)에서 정점에 도달하였다.
  • 공간 전사체학: ABCA1, GBE1, GPNMB가 저산소 핵심부 및 팔리세이딩 구역에 공간적으로 제한되어 있음을 확인한 반면, CXCL14는 더 이질적인 분포를 보였다.

의의 및 주장
본 논문은 GBM에서 가장 공격적이고 치료 저항적인 구획을 재현성 있게 표시하는 "관류-메틸로옴-전사체 축(perfusion-methylome-transcriptome axis)"을 정의한다고 주장한다. 주요 기여는 다음과 같다:

1. 기전적 통찰: 본 연구는 미세환경의 관류 구배가 유전적 드라이버와 독립적으로, GBM에서 안정적인 니치 특이적 후성유전적 재프로그래밍을 유도하는 일차적 동인임을 확립한다.
2. 후성유전적 조절: 프로모터 저메틸화가 과메틸화보다 저산소 니치에서 활성 전사를 주도하는 데 있어 지배적인 역할을 하며, 특히 온코제닉(oncogenic) 캐스케이드(Wnt, MAPK, Hedgehog)를 활성화한다는 것을 입증한다.
3. 바이오마커 발견: ABCA1, GBE1, GPNMB를 저산소 정보를 반영한 환자 층화의 후보 바이오마커로 식별하였다. 이 유전자들은 후성유전적으로 활성화된 가장 악성인 저산소 종양 하위 집단의 마커로 제시된다.
4. 치료적 함의: 저자들은 이러한 악성 표현형이 TME에 의해 유도된 후성유전적 가소성을 통해 유지되기 때문에, 이들이 후성유전적 방향의 정밀 치료 전략을 위한 잠재적 타겟이 될 수 있다고 제안한다. 본 연구는 관류 구배가 이질성을 유도하는 다른 고형암에서도 니치 특이적 다중 오믹스 연구를 위한 프레임워크를 제공한다.

저자들은 CXCL14의 경우 메틸화 이외의 경로가 조절에 관여할 수 있음을 언급하며 신중한 어조를 유지하면서도, ABCA1, GBE1, GPNMB의 견고하고 공간적으로 제한된 특성을 검증된 타겟으로서 강조하였다.