Single-cell and multi-omic analyses reveal a mito_lysis-like malignant program and a TNFRSF12A-MAP3K14-VCAM1 inflammatory remodeling axis in clear cell renal cell carcinoma

기술 요약: 단일 세포 및 멀티오믹스 분석을 통한 ccRCC에서의 mito_lysis 유사 악성 프로그램과 TNFRSF12A-MAP3K14-VCAM1 염증 재구축 축의 규명

문제 제기
투명 세포 신세포암(ccRCC)은 상당한 악성 이질성과 고도로 재구성된 종양 미세환경(TME)을 특징으로 한다. 면역 관문 억제제와 항혈관 생성 치료가 일부 환자의 예후를 개선했으나, 치료 저항성은 여전히 흔하게 나타난다. 기존의 벌크 수준 분석은 특정 악성 세포 상태와 종양 내재적 프로그램과 TME를 연결하는 복잡한 회로를 포착하지 못하는 경우가 많다. 구체적으로, 악성 상피 세포 내 미토콘드리아 스트레스, 산화 불균형, 그리고 염증 재구축 사이의 상호작용은 아직 명확히 정의되지 않았다. 또한, 골수구 풍부 니치(niche)와 스트레스 관련 악성 프로그램 사이의 리간드-수용체 신호, 그리고 이러한 신호를 지속적인 염증 및 부착 출력으로 변환하는 하위 수신 측 경로에 대해서도 완전히 이해되지 않은 상태이다.

방법론
본 연구는 단일 세포 RNA 시퀀싱(scRNA-seq), 벌크 전사체, 단백질체, 그리고 공간 전사체학을 결합한 통합 멀티오믹스 프레임워크를 활용하였다:

1. 데이터 통합 및 악성 세포 정제: 세 개의 공공 ccRCC scRNA-seq 데이터셋(GSE152938, GSE159115, GSE171306)을 Seurat과 Harmony를 사용하여 배치 교정(batch correction) 후 통합하였다. 정상 네프론 유사 세포로부터 악성 상피 세포를 구분하기 위해, 저자들은 inferCNV를 적용하여 복제수 변이(CNV)를 추론하고, 상피 세포를 이배체(diploid), 모호한(ambiguous), 고신뢰도 이배성(aneuploid) 악성 그룹으로 분류하였다.
2. 프로그램 특성 규명: 악성 세포를 대상으로 미토콘드리아 스트레스, 산화 불균형 및 염증 활성화와 관련된 유전자를 기반으로 합성된 "mito_lysis-like" 점수를 계산하였다. 고점수 세포를 특징짓기 위해 차별 발현 및 경로 농축 분석(GO, KEGG)을 수행하였다.
3. 리간드-수용체 추론: NicheNet을 사용하여 고(high) mito_lysis-like 프로그램을 조절하는 골수구 유래 송신 세포로부터의 상위 리간드를 추론하였다.
4. 가상 섭동(Virtual Perturbation): scTenifoldKnk를 활용하여 주요 수신 측 노드(IL1A, TNFRSF12A, MAP3K14, VCAM1)의 가상 녹아웃(knockout)을 수행함으로써, 하위 염증 및 부착 출력에 미치는 기능적 영향을 평가하였다.
5. 멀티오믹스 검증:
   • TCGA-KIRC: 벌크 전사체 데이터를 사용하여 단일 세포 유래 축 점수와 TNF/NF-κB 활성, 골수구 침윤(xCell 이용) 및 환자 생존율 간의 상관관계를 분석하였다.
   • CPTAC: 단백질 수준에서 VCAM1과 염증/골수구 마커 간의 상관관계를 검증하기 위해 단백질체 데이터를 분석하였다.
   • 공간 전사체학: 공공 공간 데이터(GSE210041)를 분석하여 골수구 니치의 공간적 분포와 수신/출력 염증-부착 프로그램을 매핑하였다.

주요 기여 및 결과

• 연속적인 Mito_lysis-like 프로그램의 식별: 본 연구는 ccRCC의 악성 세포에서 나타나는 연속적이고 등급이 매겨진 전사 상태인 "mito_lysis-like"를 정의한다. 이는 불연속적인 서브클론이 아니라 산화 스트레스 적응, TNF/NF-κB 신호 전달, IL-17 신호 전달, 단백질 처리 스트레스(ER 스트레스) 및 부착 관련 염증 출력이 풍부한 스펙트럼이다. 고점수 세포는 높은 수준의 스트레스 반응 유전자(예: STAT1, HMOX1, DDIT3)와 염증 매개체를 발현한다.
• 골수구 유래 염증 입력: NicheNet 분석을 통해 골수 유래 TNFSF12(TWEAK)와 IL1A를 최상위 후보 리간드로 식별하였다. TNFSF12는 골수 세포에 대한 가장 강력한 특이성을 보였으며, 그 수용체인 TNFRSF12A(Fn14)는 고(high) 점수 악성 세포에서 높게 발현되었다.
• TNFRSF12A-MAP3K14-VCAM1 축: 가상 섭동 및 네트워크 분석을 통해 수신 측 염증 재구축 축을 확립하였다:
  • TNFRSF12A: 수용체 진입점으로 작용하며, 이의 섭동은 백혈구 부착 및 활성화 유전자에 유의미한 영향을 미쳤다.
  • MAP3K14: 비정형(non-canonical) NF-κB 신호 전달 경로로 수렴하는 신호 증폭기 역할을 한다.
  • VCAM1: 원위부(distal) 염증-부착 출력 역할을 한다.
  • 이 축은 mito_lysis-like 프로그램을 ICAM1, PTGS2, NFKB2, CEBPB와 같은 하위 출력으로 연결한다.
• 임상 및 공간적 상관관계:
  • TCGA-KIRC에서 해당 축 점수는 mito_lysis-like 점수, TNF/NF-κB 활성 및 골수구 침윤과 강하게 상관되었다. TNFRSF12A와 VCAM1의 높은 발현은 낮은 전체 생존율을 예측했다.
  • CPTAC 단백질체 데이터는 VCAM1 단백질 양이 염증(STAT1, ICAM1) 및 골수구(TYROBP, AIF1) 마커와 상관관계가 있음을 확인하였다.
  • 공간 전사체학 분석 결과, 수신/출력 염증-부착 프로그램은 균일하게 분포하기보다 골수구가 풍부한 조직 니치에 농축되어 있음이 밝혀졌다.

의의
본 논문은 ccRCC에서 스트레스와 염증의 결합을 이해하기 위한 멀티오믹스 프레임워크를 제공한다고 주장한다. "mito_lysis-like" 상태를 정의함으로써, 본 연구는 ccRCC에서의 미토콘드리아 스트레스가 단순한 대사적 특징이 아니라 염증 적응과 밀접하게 연관되어 있음을 시사한다. TNFRSF12A-MAP3K14-VCAM1 축의 식별은 골수 유래 신호(특히 TNFSF12)가 어떻게 악성 상피 세포의 염증 재구축과 부착 출력을 유도하는지에 대한 구체적인 기전적 가설을 제시한다. 이 축은 불량한 임상 결과의 잠재적 바이오마커이자 종양-골수구 상호작용을 이해하기 위한 표적으로 제시된다.

한계 및 주장의 겸손함
저자들은 자신들의 발견이 일차적으로 연관성(associative)에 기반하고 있음을 명시적으로 밝히고 있다. 저자들은 다음과 같이 언급하였다:

• "mito_lysis-측" 상태는 전사체 특징으로부터 추론된 것이며, 특정 세포 사멸 과정에 대한 직접적인 증거를 구성하지 않는다.
• NicheNet 및 scTenifoldKnk를 통해 생성된 기전적 가설은 발현 기반 네트워크에 근거하며, 리간드-수용체 신호 전달이나 유전자 기능에 대한 실험적 증거가 부족하다.
• 본 연구는 소급적 공공 데이터셋에 의존하므로, 골수 세포, 악성 상피 세포, 그리고 혈관 인터페이스 사이의 인과관계를 확립할 수 없다.
• 공간 분석은 스팟(spot) 수준의 해상도로 제한되었으므로, 단일 유전자 공간 추론보다는 모듈 단위의 해석이 더 신뢰할 만하다.
  저자들은 이 축의 생물학적 및 중개적 관련성을 명확히 하기 위해 향-후 전향적 임상 검증과 실험적 섭동(in vitro/in vivo)이 필요하다고 결론지었다.