PGK1 expression differs across breast cancer molecular subtypes and is associated with overall survival in the TCGA BRCA cohort

기술 요약: 유방암 아형에서의 PGK1 발현 및 생존율

문제 제로 (Problem Statement)
Phosphoglycerate kinase 1 (PGK1)은 암 대사와 유방암 진행에 관여하는 당분해 효소이다. 이전 연구들은 PGK1 발현 증가와 불량한 임상 결과 사이의 연관성을 시사해 왔으나, 이러한 결과들은 일관되지 않았다. 저자들은 이전 보고들이 유방암의 고유한 분자적 아형(예: PAM50 분류)과 질병 단계를 고려하지 못함으로써 발생한 혼란 때문일 수 있다고 가정한다. 또한, 이진화된(dichotomized) 발현 데이터 사용, 불완전한 임상적 조정, 그리고 이질적인 데이터 파이프라인과 같은 이전 연구들의 방법론적 한계는 PGK1-생존율 연관성에 대한 불확실성을 야기했다.

방법론 (Methodology)
본 연구는 cBioPortal DataHub에서 확보한 공식 TCGA-BRCA PanCancer Atlas 2018 데이터셋을 사용하여 회고적 단일 코호트 분석을 수행하였다. 분석은 RNA-seq V2 RSEM 발현 데이터를 가진 1,082개의 원발성 종양 샘플에 집중하였다.

• 데이터 전처리: 발현 값은 log2 변환되었다. 임상 데이터에는 연령, AJCC 병리적 단계, 그리고 전체 생존(OS) 상태가 포함되었다. PAM50 아형은 Luminal A, Luminal B, HER2-enriched, Basal-like, Normal-like 범주로 매핑되었다.
• 생존 분석: 915명의 환자(118명 사망)로 구성된 완전 사례(complete-case) 코호트가 생존 모델링에 활용되었다. PGK1 발현은 표준화된 연속 변수로 취급되었다. 세 가지 Cox 비례 위험 모델이 적합되었다:
  1. 단변량 (PGK1만 포함).
  2. 연령 조정 모델.
  3. 주요 모델: AJCC 단계 및 PAM50 아형에 따라 기저 위험(baseline hazard)을 층화하여 연령을 조정한 모델. 이러한 층화는 단계 및 아형의 범주형 모델에서 관찰된 비비례 위험(non-proportional hazards) 문제를 해결하기 위해 선택되었다.
• 공통 발현 분석: 큐레이션된 119개 유전자 패널을 사용하여 당분해, HIF-1 조절, 상피-간엽 이행(EMT), 세포 주기에 대한 복합 시그니처를 계산하였다. PGK1과 이러한 시그니처 간의 피어슨(Pearson) 및 편 상관관계(partial correlation, 아형 조정 후)가 계산되었다.
• 통계적 엄밀성: 분석은 Python을 사용하여 수행되었다. 본 연구는 REMARK 가이드라인을 준수하였으며, 양측 검정을 사용하였고, 완전한 재현성 스크립트와 체크섬을 제공하였다.

주요 결과 (Key Results)

• 아형별 발현: PGK1 발현은 PAM50 아형에 따라 유의미하게 달랐다 ($P = 4.35 \times 10^{-57}$). 발현은 HER2-enriched 종양(평균 13.76)에서 가장 높았으며, Basal-like (13.41), Luminal B (13.16), Normal-like (12.82), Luminal A (12.76) 순이었다.
• 생존 연관성: 높은 PGK1 발현은 유의미하게 짧은 전체 생존율과 연관되었다.
  • 단변량 모델에서, 1-SD 증가당 위험비(hazard ratio, HR)는 1.55 ($P = 7.10 \times 10^{-7}$)였다.
  • 주요 층화 모델(연령, 단계 및 아형을 조정한 모델)에서도 연관성은 통계적으로 유의미하게 유지되었으며, HR은 1.41 (95% CI, 1.13–1.76; $P = 0.002$)이었다. 이는 단계 및 분자적 아형과 독립적으로, PGK1 발현의 1 표준 편차당 위험이 41% 증가함을 나타낸다.
• 공통 발현 패턴: PGK1은 당분해 복합 점수와 가장 강한 양의 상관관계를 보였다 ($r = 0.675$; 아형 조정 후 편 상관관계 $r = 0.575$). 또한 세포 주기 및 HIF-1 조절 점수와 양의 상관관계를 보였다. 반대로, PGK-1은 루미날 호르몬 수용체 마커(ESR1, PGR) 및 BCL2와 음의 상관관계를 보였다.

의의 및 주장 (Significance and Claims)
저자들은 PGK1 발현이 유방암 전반에 걸쳐 균일하지 않으며, 특히 공격적인 HER2-enriched 및 Basal-like 종양에서 가장 높은 수준을 보이는 등 분자적 아형에 따라 실질적으로 달라진다고 결론지었다. 결정적으로, 본 연구는 연령, AJCC 단계 및 PAM50 아형에 대한 엄격한 조정을 거친 후에도 높은 PGK1 발현과 불량한 전체 생존율 사이의 연관성이 지속됨을 입증하였다.

본 논문은 이러한 발견을 **탐색적(exploratory)**인 것으로 규정한다. 저자들은 본 연구가 PGK1을 위험 층화의 잠재적 바이오마커로 뒷받침하고 있으나, 인과적 조절, 임상적 유용성 또는 치료적 관련성을 확립하는 것은 아니라고 명시적으로 밝히고 있다. 저자들은 이러한 단일 코호트 연구 결과가 임상적 결론을 도출하기 전에 독립적인 임상 검증과 기능적 실험 검증을 필요로 한다는 점을 강조한다. 본 연구의 주요 기여는 공식 소스 파일을 사용하고, 연속적 발현 모델링을 적용하며, 비비례 위험을 명시적으로 처리함으로써 방법론적 엄밀성을 확보하여 유방암의 이질성 맥락 내에서 PGK1의 독립적인 예후 가치를 명확히 한 데 있다.