Integrative transcriptomic profiling of the phosphatome identifies DUSP15 as a chromophobe renal cell carcinoma-selective epithelial biomarker associated with an immune-metabolic tumour state

기술 요약: 인산가수분해효소(Phosphatome)의 통합적 전사체 프로파일링을 통한 색소성 신세포암 특이적 상피 바이오마커 DUSP15의 식별

문제 제기
색소성 신세포암(chRCC; KICH)은 신장 악성 종양의 약 5%를 차지하며, 특히 진행성 또는 전이성 질환에서 강력한 분자적 바이오마커와 작용 가능한 치료 표적이 부족하다는 특징이 있다. 투명세포형(KIRC) 및 유두형(KIRP) 신세포암과 달리, chRCC의 분자적 결정 요인은 아직 명확히 정의되지 않았다. 또한, 벌크 전사체 프로파일은 질병 선택적 발현 패턴을 식별할 수는 있지만, 해당 신호가 종양 상피에서 기원하는지, 정상 신장 계통에서 기원하는지, 혹은 침윤된 면역 및 기질 세포에서 기원하는지를 구별하지 못한다. 이러한 한계는 chRCC 특이적인 생물학적 취약성을 정밀하게 규명하고 정밀 치료 전략을 개발하는 데 걸림돌이 된다.

연구 방법론
본 연구는 인산가수분해효소(phosphatome)를 체계적으로 조사하고 후보 바이오마커의 세포 기원을 규명하기 위해 다층적인 통합 접근법을 채택하였다:

1. 범암 인산가수분해효소 스크리닝: 265개의 단백질 코딩 인산가수분해효소 유전자 목록을 선정하여, 21개 암종에 대한 TCGA(The Cancer Genome Atlas) 전사체 데이터를 사용하여 분석하였다. 조직 특이성은 Tau 지수를 통해 정량화하였다.
2. 벌크 전사체 검증: DUSP15 발현을 KICH, KIRC, KIRP의 조직학적 아형 주석이 명시된 253개의 신장 종양 샘플을 포함하는 10개의 독립적인 GEO(Gene Expression Omnibus) 마이크로어레이 코호트에서 검증하였다.
3. 진단 및 통계 분석: 진단 정확도와 아형 식별 능력을 평가하기 위해 ROC 분석, 로지스틱 회귀, 주성분 분석(PCA)을 수행하였다. 차등 발현은 Wilcoxon rank-sum test와 DESeq2를 사용하여 평가하였다.
4. 전사체 전체 상관관계 및 경로 분석: TCGA-KICH 내 DUSP15와 모든 단백질 코딩 유전자 간의 Spearman 상관관계를 계산하였다. 연관된 전사 프로그램 정의를 위해 Hallmark 유전자 세트를 이용한 GSEA(Gene Set Enrichment Analysis)와 Gene-Concept Network 시각화를 수행하였다.
5. 면역 탈포획(Immune Deconvolution): xCell 알고리즘을 적용하여 64가지 면역 및 기질 세포 유형의 상대적 농축도를 추정하고, DUSP15-high 그룹과 DUSP15-low 그룹으로 종양을 층화하였다.
6. 단일 세포 RNA 시퀀싱(scRNA-seq) 해상도: 두 개의 독립적인 scRNA-seq 데이터셋(GSE152938 및 GSE159115)을 분석하여 세포 해상도 수준에서 DUSP15 발현 위치를 국소화하였다. 여기에는 품질 관리, 클러스터링, 그리고 매칭된 정상 신장, ccRCC, pRCC 구획과 DUSP15 풍부도를 비교하기 위한 pseudobulk 분석이 포함되었다. 면역-대사 상태를 평가하기 위해 공여자별(Donor-resolved) 경로 분석을 수행하였다.

주요 결과

• 선택적 바이오마커로서의 DUSP15 식별: 265개의 인산가수분해효소 중 DUSP15는 0.977의 Tau 점수와 함께 가장 높은 chRCC 선택성을 나타냈다. DUSP15는 다른 종양 유형에 비해 KICH에서 매우 강력하고 독점적인 발현을 보였다(TCGA-KICH 샘플의 98.4%에서 검출).
• 진단 성능: TCGA 코호트에서 DUSP15는 0.976의 AUC(Area Under the Curve)로 KICH를 KIRC 및 KIRP와 구별하였다. GEO 마이크로어레이 코호트에서의 독립적 검증에서도 0.837의 AUC를 기록하였다. 로지스틱 회귀 분석 결과, DUSP15 발현량이 1단위 증가할 때마다 KICH 진단의 오즈비(odds)가 4.83배 증가하는 것으로 나타났다.
• 전사적 경관: 높은 DUSP15 발현은 산화적 인산화, MYC 타겟, DNA 복구 경로의 협응된 상향 조절과 관련이 있었다. 반대로, 인터페론-γ 반응, TNF-α/NF-κB 신호 전달 및 염증 반응 경로와는 음의 상관관계를 보였다.
• 면역-대사 상태: 면역 탈포획 결과, DUSP15-high 종양은 조절 T세포(Treg) 및 수지상 세포 신호가 감소하고, 전반적인 ImmuneScore 및 MicroenvironmentScore가 낮은 특징을 가진, 면역 미세환경이 현저히 결핍된 상태임을 밝혔다.
• 세포 국소화: scRNA-seq 분석을 통해 DUSP15 신호가 주로 chRCC 종양 상피에서 기원함을 확인하였다. GSE152931에서 DUSP15는 chRCC 종양 상피 세포의 72.2%(144.9 CPM)에서 검출되었으며, 면역, 기질 및 정상 신장 상피 구획에서는 발현이 현저히 낮았다. 이 패턴은 독립적인 GSE159115 데이터셋에서도 재현되었다(종양 세포에서 83.7% 검출; 204.3 CPM).
• 상피 맥락: scRNA-seq 데이터의 공여자별 경로 분석 결과, chRCC 상피는 정상 신장 상피와 비교했을 때의 보편적인 억제가 아니라, ccRCC 상피와 비교했을 때 특이적으로 대사/증식 프로그램의 농축과 염증 신호의 감쇄를 보이는 것으로 확인되었다.

의의 및 주장
본 논문은 DUSP15를 색소성 신세포암에 대한 고도로 선택적인 종양-상피 바이오마커 후보로 식별하였다. 본 연구는 DUSP15의 높은 벌크 종양 신호가 기질이나 면역 침윤이 아닌, 광범위한 종양 상피 내 발현에 기인한다는 점을 입증하였다.

저자들은 다음과 같은 특징을 가진 "DUSP15 관련 면역-대사 상태"를 정의하였다:

1. 대사 활성화: 산화적 인산화, MYC 타겟 및 DNA 복구 프로그램의 농축.
2. 면역 감쇄: 특히 ccRCC와의 비교 맥락에서 나타나는 염증 및 면역 세포 신호의 감소.

본 연구는 DUSP15를 높은 진단 정확도로 KICH를 다른 신장 아형과 구별할 수 있는 강력한 분자적 분류 도구로 위치시킨다. 이러한 결과는 DUSP15, 종양 대사, 면역 조절 및 chRCC 상피 정체성 사이의 관계에 대한 추가적인 기전 연구를 위한 집중적인 근거를 제공한다. 저자들은 결과가 여러 코호트와 플랫폼에 걸쳐 통계적으로 견고하지만, 임상적 유용성과 정밀한 종양 미세환경의 구조를 완전히 확립하기 위해서는 전향적 임상 검증 및 조직학적/공간적 검증이 필요하다고 언급하였다.