Integrated Analysis of Bulk and Single Cell Transcriptomic Sequencing Reveals Immune Patterns of T Cell Regulation and Prognostic Biomarkers in Glioblastoma

기술 요약: 교모세포종의 벌크 및 단일 세포 전사체 시퀀싱 통합 분석

문제 제기
교모세포종(GBM)은 성인 중추 신경계에서 발생하는 가장 공격적인 원발성 악성 종양으로, 신규 진단 환자의 5년 생존율이 10% 미만인 불량한 예후를 특징으로 합니다. 수술적 절제, 방사선 치료, 테모졸로마이드 화학요법을 포함한 표준 중재에도 불구하고, 면역 억제성 종양 미세환경(TME)은 면역 관문 억제제(ICI) 및 CAR-T 세포와 같은 면역 요법의 효능을 심각하게 제한합니다. 교모세포종 미세환경 내 T 세포의 이질성, 다른 세포군과의 상호작용 패턴, 그리고 기능 부전을 유도하는 구체적인 분자 기전을 이해하는 데에는 결정적인 공백이 존재합니다. 따라서 강력한 예후 바이오마커를 식별하고 개인화된 임상 중재를 가이드할 수 있는 모델을 개발하는 것이 시급합니다.

연구 방법론
본 연구는 공공 저장소(TCGA, GEO)의 벌크 RNA 시퀀싱(bulk RNA-seq) 및 단일 세포 RNA 시퀀싱(scRNA-seq) 데이터를 통합하는 멀티 오믹스 접근 방식을 채택하였습니다. 분석 워크플로우는 다음과 같이 진행되었습니다:

1. WGCNA를 통한 핵심 유전자 식별: 벌크 전사체 데이터셋(GSE4290 및 GSE50161)에 대해 가중 유전자 공발현 네트워크 분석(WGCNA)을 수행하여 공발현 네트워크를 구축했습니다. GBM 상태와 높은 상관관계를 보이는 모듈을 식별하였으며, 이를 차등 발현 분석과 교차 검증하여 47개의 GBM 관련 핵심 유전자 세트를 도출했습니다.
2. 단일 세포 해상도 및 주석 달기: 4개의 정상 조직과 4개의 GBM 뇌 조직으로 구성된 scRNA-seq 데이터(GSE162631)를 Seurat 패키지를 사용하여 처리했습니다. 엄격한 품질 관리와 차원 축소(PCA, UMAP, t-SNE)를 거친 후, 세포를 클러스터링하고 내피세포, 섬유아세포, 대식세포, 악성 세포, 미세아교세포, T 세포의 6가지 주요 하위 유형으로 주석을 달았습니다.
3. T 세포 특이적 분석: 47개의 핵심 유전자를 UCell 알고리즘을 사용하여 단일 세포 데이터에 매핑하여 농축도를 결정했습니다. T 세포가 이 핵심 유전자들의 주요 운반체임을 확인했습니다. 이후 분석에는 다음이 포함되었습니다:
   • 차등 발현: 100개의 T 세포 특이적 차등 발현 유전자(DEG) 식별.
   • 궤적 및 통신: T 세포 상태 전이를 매핑하기 위한 의사시간(Pseudotime) 궤적 분석(Monocle 2) 및 T 세포와 기타 기질 구성 요소 간의 리간드-수용체 상호작용을 정량화하기 위한 CellChat 분석.
   • 기능적 농축: 생물학적 기능을 규명하기 위한 GO 및 KEGG 경로 분석.
4. 하위 유형 분류 및 예후 모델링:
   • 컨센서스 클러스터링: T 세포 DEG 발현을 기준으로 TCGA-GBM 샘플을 두 가지 분자적 하위 유형(C1 및 C2)으로 계층화했습니다.
   • 면역 및 약물 프로파일링: ESTIMATE, 다양한 면역 침윤 알고리즘(CIBERSORT, XCELL 등) 및 약물 민감도 예측(oncoPredict)을 사용하여 하위 유형을 평가했습니다.
   • 리스크 모델 구축: T 세포 DEG를 사용하여 5개 유전자로 구성된 예후 시그니처를 도출하기 위해 LASSO-Cox 회귀 모델을 구축했습니다. 모델의 성능은 Kaplan-Meier 생존 분석, ROC 곡선 및 Human Protein Atlas(HPA) 데이터베이스를 이용한 독립적인 단백질 수준 검증을 통해 확인되었습니다.

주요 기여 및 결과

• T 세포 중심 핵심 유전자 식별: 본 연구는 GBM과 관련된 47개 유전자의 핵심 세트를 성공적으로 식별했습니다. 단일 세포 분석 결과, T 세포가 모든 세포 성분 중 이 유전자 세트에 대해 가장 높은 농축 점수를 나타내어, GBM 면역 환경에서 T 세포의 중심적인 역할을 강조했습니다.
• T 세포 고갈의 기전적 통찰: 의사시간 분석은 GBM 내 T 세포가 세포 독성 감소와 면역글로불린 유전자(예: IGHA2, IGHG2)의 비정상적 상향 조절을 특징으로 하는 말기 고갈 상태로 진화함을 보여주었습니다. 세포 간 통신 분석 결과, 섬유아세포가 T 세포로 유입되는 신호의 주요 원천이며, COL 패밀리(콜라겐)가 T 세포의 CD44 및 인테그린 수용체를 조절하는 강력한 조절자로 작용하여 기질 매개 면역 억제 기전을 시사함을 밝혔습니다.
• 분자적 하위 유형(C1 vs. C2): 본 연구는 T 세포 유전자 발현을 기반으로 두 가지 뚜렷한 GBM 하위 유형을 정의했습니다:
  • C1 하위 유형: 면역 관문 유전자의 높은 발현과 PI3K 억제제(예: Alpelisib, Dactolisib) 및 알킬화제(Carmustine)에 대한 민감성을 특징으로 합니다.
  • C2 하위 유형: 더 높은 면역 침윤(대식세포, Treg)과 높은 면역 점수를 보였으며, 후성유전학적 조절제(HDAC 억제제인 Vorinostat 등) 및 CDK4/6 억제제(Ribociclib)에 대한 더 높은 민감도를 나타냈습니다.
• 예후 5-유전자 시그니처: GNG11, CD5, COL1A1, CAMK4, BCAS4의 5개 핵심 유전자를 사용하여 견고한 예후 모델을 구축했습니다.
  • 위험 인자: GNG11과 CD5의 높은 발현은 불량한 예후와 상관관계가 있었습니다.
  • 보호 인자: COL1A1, CAMK4, BCAS4의 높은 발현은 더 나은 생존 결과와 상관관계가 있었습니다.
  • 검증: 모델은 AUC 값 0.59(1년), 0.64(3년), 0.69(5년)로 예측 정확도를 입증했습니다. HPA 데이터베이스의 면역조직화학 분석은 5개 유전자의 단백질 발현 수준이 정상 뇌 조직에 비해 GBM 조직에서 유의하게 상향 조절되었음을 확인하여, 전사체 결과를 단백질 수준에서 검증했습니다.

의의 및 주장
저자들은 본 연구가 교모세포종의 T 세포 매개 면역 미세환경의 이질성을 체계적으로 규명했다고 주장합니다. 벌크 및 단일 세포 데이터를 통합함으로써, 본 연구는 종양 기질(특히 섬유아세포 및 ECM)과 T 세포 고갈 사이의 복잡한 상호작용을 이해하기 위한 이론적 토대를 제공합니다. 5-유전자 예후 모델의 개발 및 검증은 임상적 예후 평가를 위한 실질적인 도구를 제공합니다. 또한, 하위 유형별 약물 민감도 식별은 교모세포종 환자를 위한 개인화된 임상 중재, 진단, 예후 평가 및 표적 면역 요법 전략을 개선할 수 있는 경로를 제시합니다. 본 연구는 회고적 공공 데이터에 의존한다는 한계와 향후 인비보(in vivo)/인비트로(in vitro) 검증의 필요성을 인정하면서도, 본 연구의 결과가 GBM 바이오마커 및 치료 표적 탐색을 위한 기초적인 프레임워크를 제공한다고 주장합니다.