Single-section spatial hypoxia-cytotoxic associations do not consistently reproduce across breast cancer patients

기술 요약: 단일 섹션의 공간적 저산소증-세포독성 연관성은 유방암 환자들 사이에서 일관되게 재현되지 않음

문제 제격(Problem Statement)
공간 전사체학은 온전한 조직 내의 국소적인 종양-면역 관계를 연구할 수 있게 해준다. 그러나 중요한 통계적 과제가 존재한다: 단일 조직 섹션 내의 수천 개 측정 지점(spot)들은 단일 생물학적 표본으로부터 나온 공간적으로 의존적인 관측치이며, 독립적인 생물학적 반복구가 아니다. 이 지점들을 환자 수준의 주장을 뒷받침하는 독립적인 증거로 취급하는 것은 "의사 반복(pseudoreplication)"에 해당한다. 결과적으로, 하나의 섹션 내에서 지점 수준(spot-level)에서 관찰된 강력하고 통계적으로 유의미한 연관성이 서로 다른 환자들 사이에서 일관된 생물학적 효과를 반영하지 못할 수 있다. 본 연구는 이러한 단일 섹션 내의 공간적 연관성과 환자 수준의 재현성 사이의 간극을 유방암 연구를 통해 다룬다.

방법론(Methodology)
저자들은 공개된 유방암 데이터셋을 사용하여 2단계 이차 분석 설계를 채택하였다:

1. 발견 단계 (Discovery Phase):

   • 데이터셋: 단일 10x Genomics Visium 유방암 섹션 (3,798개 지점).
   • 유전자 프로그램: 6개의 사전 지정된 전사 프로그램이 계산되었다: 종양 상피, 저산소증(6개 유전자), 세포독성 림프구(8개 유전자), 골수구, 섬유아세포/ECM, 그리고 EMT-조절 프로그램.
   • 분석: 저자들은 공간적 이웃을 정의하기 위해 대칭화된 6-최근접 이웃 그래프를 구축하였다. 저산소증 관련 전사가 공간적 자기상관(Moran's I)을 평가하고, 저산소증 점수와 이웃 지점들의 평균 세포독성 점수 사이의 연관성을 Spearman 상관계수와 순위-이분 상관 효과 크기(high-hypoxia 지점과 low-hypoxia 지점 비교)를 사용하여 평가하였다.

2. 외부 검증 단계 (External Validation Phase):

   • 데이터셋: 8명의 서로 다른 환자로부터 얻은 36개의 섹션(총 13,619개 지점)으로 구성된 독립적인 HER2 양성 유방암 코호트 (Andersson et al.).
   • 프로토콜: 각 섹션은 독립적으로 처리되었다 (라이브러리 크기 정규화, 로그 변환, 유전자별 표준화). 저산소증과 이웃 세포독성 사이의 방향성 연관성이 각 섹션별로 계산되었다.
   • 집계: 섹션 수준의 상관관계는 환자 수준으로 중앙값(median)을 사용하여 집계되었다.
     통계적 검정: 일방향 exact Wilcoxon 부호 순위 검정을 사용하여 환자 수준의 중앙값 상관관계가 0보다 유의하게 작은지 결정하였다.
   • 민감도 분석: 저자들은 이웃 크기 변화(4, 8, 10개 이웃), effector-only 세포독성 시그니처 적용, 그리고 컴팩트한 저산소증 시그니처를 200개 유전자로 구성된 MSigDB Hallmark 저산소증 세트로 대체함으로써 견고성을 테스트하였다.

주요 결과 (Key Results)

• 발견 섹션 (Discovery Section): 명확한 역상관 관계가 관찰되었다. 저산소증 관련 전사는 이웃한 세포독성 유전자 활성과 음의 상관관관계(Spearman $\rho = -0.202$)를 보였다. 고저산소증 지점은 저산소증 지점보다 유의하게 낮은 이웃 세포독성 점수를 나타냈다 (rank-biserial effect = $-0.286$).
• 환자 수준 검증 (Patient-Level Validation): 역상관 관계가 외부 코호트에서 일관되게 재현되지 않았다.
  • 36개 섹션 중 19개, 8명 환자 중 5명만이 음의 연관성을 보였다.
  • 환자 수준의 중앙값 상관관계는 $-0.043$이었으며, 이는 0과 통계적으로 구별되지 않았다 (one-sided exact Wilcoxon $P = 0.473$).
  • 환자 간 이질성이 상당했다. 예를 들어, 환자 B는 균일하게 양의 섹션 수준 연관성을 보인 반면, 환자 A는 균일하게 음의 연관성을 보였다.
• 민감도 분석:
  • 이웃 크기를 변경하거나 Kendall 상관관계를 사용했을 때 약하고 이질적인 환자 효과가 나타났으며 (중앙값 범위 $-0.025$~$-0.032$), 모두 유의하지 않았다.
  • MSigDB Hallmark 저산소증 시그니처를 사용했을 때는 추세가 역전되어 양의 중앙값 환자 상관관계($0.047$)를 생성하였으며, 이는 관찰된 방향성이 "저산소증"의 특정 전사적 정의에 매우 의존적임을 나타낸다.

주요 기여 (Key Contributions)

• 의사 반복 위험의 입증: 본 연구는 단일 공간 전사체 섹션 내의 지점 수준 유의성이 환자 수준의 생물학적 일관성을 보장하지 않는다는 실증적 증거를 제공한다.
• 방법론적 프레임워크: 환자를 반복 단위로 취급하여 공간적으로 의존적인 데이터를 결합함으로써 통계적 검정력을 부풀리는 것을 방지하는, 공간 전사체학을 위한 엄격한 2단계 워크플로우를 확립하였다.
• 이질성 특성화: 공간적 종양-면역 관계(특히 저산소증-세포독성)가 환자 간에 매우 이질적이며, 유전자 시그니처 및 이웃 정의의 선택에 민감하다는 것을 강조한다.

의의 및 주장 (Significance and Claims)
본 논문은 저산소증과 세포독성 배제 사이의 강력한 단일 섹션 연관성이 유방암 환자 전체에 걸친 일반적인 규칙이 아니라고 주장한다. 대신, 증거는 보편적인 저산소증 관련 세포독성 배제 패턴보다는 **공간적 맥락 의존성(spatial context dependence)**을 지지한다.

저자들은 외부 검증 없이 단일 섹션의 지점 수준 공간 연관성을 환자 수준의 생물학적 효과로 해석하는 것에 대해 경고한다. 그들은 환자 간 이질성이 기술적 번거로움이 아니라 생물학적 결과로 취급되어야 한다고 주장한다. 본 연구는 향 far의 공간 전사체 연구가 단순히 단일 섹션 내의 지점 수를 늘리는 것이 아니라, 독립적인 환자를 추가하고 환자 수준의 추정치(estimands)를 사전에 지정하는 것을 우선시해야 한다고 결론짓는다. 연구 결과는 전사체만으로 된 공간적 연관성이 추가적인 검증 없이 대사 경쟁과 세포 풍부도 또는 조직 구조의 차이와 같은 메커니즘을 구분할 수 없음을 시사한다.