An Integrative Transcriptomic Analysis of a Stemness-Associated Gene Set on Ovarian Cancer Recurrence

기술 요약: 난소암 재발에 대한 줄기세포 관련 유전자 세트의 통합 전사체 분석

문제 제기
난소암(OC)은 진행된 단계에서 5년 생존율이 30% 미만인 매우 치명적인 부인과 암으로 남아 있다. 백금 기반 화학요법에 대한 초기 반응은 양호함에도 불구하고, 환자의 70~80%가 질병 재발을 경험한다. 이러한 재발은 빈번하게 화학요법에 내성을 가진 암 줄기세포(CSC)의 잔존에 기인한다. 특정 전사 인자(예: NANOG, SOX2) 및 대사 조절 인자(예: ALDH 패밀리)가 줄기세포성(stemness)을 유지하는 것으로 알려져 있으나, 일차 질환에서 재발성 질환으로의 전이를 구체적으로 특징짓는 일관된 유전자 발현 패턴을 식별할 필요가 있다. 본 연구는 정의된 6개의 줄기세 세포 관련 유전자 패널이 일차 난소암(POC) 대비 재발성 난소암(ROC)에서 재현 가능한 발현 변화를 보이는지 결정하여, 종양 재발의 기저에 깔린 전사적 및 대사적 재프로그래밍을 규명하는 것을 목표로 한다.

연구 방법론
본 연구는 대규모 공공 전사체 데이터와 실험적 검증을 결합한 통합적 접근 방식을 채택하였다:

1. 데이터 수집 및 코호트:

   • TCGA 코호트: 421개의 일차 난소암(POC) 샘플과 62개의 재발성 난소암(ROC) 샘플로부터 얻은 RNA-seq 데이터(FPKM 정규화).
   • GEO 참조 데이터: 고등급 장액성 난소암(HGSOC) 일차 종양에서 유래하여 실험적으로 검증된 난소암 줄기세포(OCSC)의 데이터셋(GSM4684836).
   • 임상 검증 코호트: 이스탄불 대학교 의과대학에서 수집된 8명의 POC 환자 및 8명의 ROC 환자의 신선 조직 샘플(각각 1차 화학요법 전 및 문서화된 재발 후 수집).

2. 대상 유전자 패널:
   분석은 줄기세포성, 자기 복제 및 화학요법 내성과 관련이 있는 것으로 알려진 6개의 유전자, 즉 ALDH1A1, ALDH1A2, ALDH1A3, NANOG, PROM1 (CD133), 및 SOX2에 초점을 맞추었다.

3. 컴퓨터 분석:

   • 전처리: FPKM 값을 필터링하고(샘플의 80% 이상에서 FPKK < 10인 유전자 제거) $\log_2$ 변환 $[\log_2(\text{FPKM}+1)]$을 수행하였다.
   • 차등 발현 분석: POC와 ROC를 비교하기 위해 허위 발견율(FDR) 보정이 적용된 모수적 통계 테스트를 사용하였다. 유의성은 $|\log_2 \text{fold change}| > 1$ 및 조정된 $p < 0.05$로 정의되었다.
   • 상관관계 분석: POC와 ROC 전사체가 OCSC 참조 시그니처와 얼마나 유사한지 평가하기 위해 피어슨 상관계수를 계산하였다.
   • 생존 분석: TCGA 환자를 대상으로 합성 줄기세포 점수(6개 유전자의 평균 z-score)를 산출하였다. 전체 생존율(OS)을 평가하기 위해 Kaplan-Meier 곡선과 Cox 비례 위험 회귀 모델을 사용하였다.

4. 실험적 검증:

   • 임상 샘플로부터 총 RNA를 분리하고 이를 cDNA로 전환하였다.
   • GAPDH를 참조 유전자로 사용하여 qRT-PCR을 통해 유전자 발현을 정량화하였다.
   • 상대적 발현량은 $\Delta\Delta\text{Cq}$ 방법을 사용하여 계산하였으며, qRT-PCR과 RNA-seq의 fold-change 간의 일치도를 Spearman 상관계수로 평가하였다.

주요 결과

• TCGA 데이터에서의 차등 발현:

  • ROC에서 상향 조절됨: SOX2와 PROM1은 일차 종양에 비해 재발성 종양에서 유의미한 상향 조절을 보였다 (조정된 $p < 0.01$).
  • ROC에서 하향 조절됨: ALDH1A1 발현은 POC에 비해 ROC에서 유의하게 감소하였다 ($p < 0.01$).
  • 유의한 변화 없음: ALDH1A2, ALDH1A3, 및 NANOG는 TCGA 코호트에서 통계적으로 유의미한 차등 발현을 보이지 않았다.

• 줄기세포 시그니처와의 상관관계:

  • 일차 종양은 OCSC 시그니처와 음의 상관관계($r = -0.59$)를 보였으며, 이는 줄기세포와 구별되는 발현 프로파일을 나타낸다.
  • 재발성 종양은 거의 중립적인 상관관계($r = -0.04$)를 보였으며, 이는 재발 과정에서 OCSC와 유사한 전사 상태로 수렴함을 시사한다.

• 생존 분석:

  • 높은 합성 줄기세포 점수를 가진 환자들은 더 나쁜 전체 생존율을 보이는 경향이 있었으나, 이 차이는 통계적으로 유의하지 않았다 ($\text{log-rank } p = 0.11$).

• qRT-PCR 검증:

  • 일치하는 결과: SOX2와 PROM1은 재발 임상 코호트에서도 유의하게 상향 조절되어 ($p < 0.01$ 및 $p < 0.05$), TCGA의 결과를 뒷받침하였다.
  • 상충하는 결과: ALDH1A1은 재발 qRT-PCR 코호트에서 통계적으로 유의한 증가를 보여 ($p < 0.05$), TCGA 데이터셋에서 관찰된 하향 조절과 대조를 이루었다.
  • 기타 유전자: ALDH1A3는 qRT-PCR 코호트에서 유의하게 상승했으나 ($p < 0.05$) TCGA에서는 그렇지 않았다. ALDH1A2는 유의하지 않은 증가를 보였고, NANOG는 두 데이터셋 모두에서 중간 정도의 비유의적인 증가를 보였다.
  • 일치도: RNA-seq와 qRT-PCR 플랫폼의 fold-change 값 사이에서 약하거나 중간 정도의 양의 상관관계 ($\text{Spearman } r = 0.42$)가 관찰되었다.

의의 및 주장
저자들은 난소암 재발이 특정 하위 형태의 다능성 조절 인자, 특히 SOX2–PROM1 축의 재활성화에 중심을 둔 특정한 전사적 재프로그래밍에 의해 유도된다고 주장한다.

• 강건한 바이오마커: 두 대규모 RNA-seq 및 소규모 qRT-PCR 데이터셋 모두에서 일관되게 상향 조절된 SOX2와 PROM1은 재발 관련 줄기세포 활성화의 신뢰할 수 있는 바이오마커로서의 위치를 확립한다.
• 대사 조절 인자의 이질성: 본 연구는 ALDH 이소형들이 이질적인 조절 양상(예: 서로 다른 코호트에서 반대되는 경향을 보이는 ALDH1A1)을 보임을 강조하며, 이는 이들의 발현이 균일한 재발 마커라기보다 화학요법 노출, 종양 아형 및 미세환경 요인에 의해 크게 영향을 받음을 시사한다.
• 기전적 통찰: (POC의) 음의 상관관계에서 (ROC의) 중립적 상관관계로의 변화는 재발이 별개의 계통에서 발생하는 것이 아니라, 종양 세포의 일부가 점진적으로 줄기세포와 유사한 전사적 특성을 다시 획득함으로써 발생함을 시사한다.
• 치료적 함의: 저자들은 SOX2–PROM1 전사 줄기세포 시그니처를 표적으로 삼는 것이 내성 세포 집단을 제거하고 장기적인 결과를 개선하는 효과적인 전략이 될 수 있다고 제안하며, 현재의 결과가 예비적 단계이며 추가적인 기능적 검증이 필요함을 인정한다.

본 논문은 TCGA의 ROC 코호트 샘플 크기가 작고 검증된 임상 샘플 수가 제한적이라는 점을 언급하며, 이것이 생존 분석의 일반화 및 미세한 전사적 차이 탐지에 영향을 미칠 수 있음을 명시하며 신중한 어조를 유지한다.