Correlation of OCT-Based Radiomic Signatures With Dose-Associated Radiation Response in Tumor Spheroids

기술 요약: 종양 구형체(Tumor Spheroids)의 선량 관련 방사선 반응과 OCT 기반 라디오믹스 시그니처의 상관관계

문제 제기
3차원 종양 구형체는 2D 세포 배양과 in vivo 조직 사이의 중요한 가교 역할을 하며, 방사선 생물학 연구에서 세포 간 상호작용 및 확산 기반 구배(예: 산소, 영양분)를 정확하게 재현한다. 그러나 선량 의존적 방사선 반응을 정량화하기 위한 확장 가능하고 재현 가능한 판독법에는 상당한 공백이 존재한다. 전통적인 명시야(Brightfield, BF) 이미징은 2D 투영만을 제공하므로 내부의 3D 구조나 깊이 의존적 변화를 해결하지 못한다. 형광 현미경은 특이성은 높지만, 외부 프로브가 살아있는 샘플을 교란할 수 있으며 주로 종단적 측정에 국한되는 경우가 많다. 공초점 현미경은 광학 단층 촬영을 제공하지만, 시간 소모가 크고 두꺼운 샘플에서는 신호 손실이 발생한다. 결과적으로, 조사(irradiation)로 유도된 복잡한 3차원적 및 동적 변화를 포착할 수 있는 라벨이 필요 없는(label-free), 비파괴적이며 고처리량(high-throughput)인 방법이 필요하다.

방법론
본 in vitro 단면 연구에서는 두 가지 밀도(5,000개 및 10,000개 세포)로 배양된 SAS 인간 설암 세포 구형체를 사용하였으며, 0에서 12 Gy 범위의 선량으로 조사하였다. 이미징은 조사 후 1, 4, 6, 8, 11일에 수행되었다.

• 이미징 양식:
  • 명시야(Brightfield): 기준점 역할을 하기 위해 ImageJ를 사용하여 2D 형태학적 지표(직경 및 투영 면적)를 수동으로 측정하였다.
  • 광간섭 단층 촬영(Optical Coherence Tomography, OCT): 스펙트럼 도메인 OCT 시스템(중심 파장 880 nm)을 사용하여 깊이 분해 볼륨 데이터를 획득하였다. 각 획득 시에는 두 개의 공동 등록된 볼륨, 즉 구조적 강도(intensity) 볼륨과 스펙클 분산(speckle-variance, SV) 볼륨이 생성되었다. SV 볼륨은 세포 역학 및 미세환경 활동과 관련된 시간적 강도 변동을 포착하기 위해 20회의 반복된 B-스캔 프레임으로부터 계산되었다.
• 라디오믹스 파이프라인:
  • 세그멘테이션(Segmentation): OCT 볼륨은 Otsu 임계값 처리와 형태학적 연산을 사용하여 크롭 및 마스킹되어 구형체 조직을 분리하였다.
  • 특징 추출(Feature Extraction): Radiomics.jl을 사용하여 강도 및 SV 볼륨으로부터 형태(강도에서만 추출), 1차 통계량, 그리고 텍스처 기술자(GLCM, GLRLM, GLSZM, GLDM)를 포함한 특징들을 추출하였다.
  • 견고성 필터링(Robustness Filtering): 단일 구형체에 대한 6회의 반복 스캔 실험을 통해 사분위수 범위 대 중앙값 비율(IQR/median)이 35%를 초과하는 특징을 필터링하여, 초기 198개 특징 중 24개를 제거하였다.
  • 특징 선택(Feature Selection): 특징들은 상관관계 제거(Spearman $|\rho| \ge 0.95$) 및 순열 중요도(permutation importance)와 단변량 점수에 기반한 앙상블 랭킹을 거쳤다. 공정한 비교를 위해 "OCT 전용" 및 "결합된 OCT+BF" 특징 세트에 대해 독립적으로 선택을 수행하였다.
• 모델링: 5가지 회귀 모델(Random Forest, Extra Trees, SVR, LightGBM, XGBoost)을 사용하여 조사 선량을 예측하도록 학습시켰으며, 반복 5-겹 교차 검증을 통해 BF 전용, OCT 전용, 결합된 특징 세트를 비교하여 성능을 평가하였다.

주요 기여

1. OCT 라디오믹스의 검증: 본 연구는 OCT 유래 라디오믹스 특징(구조적 강도 및 스펙클 분산)이 기존의 명시야 형태학보다 종양 구형체의 방사선 선량을 더 정확하게 예측할 수 있음을 입증한다.
2. 라벨이 없는 고처리량 판독법: 외부 프로브나 파괴적 샘플링 없이 종단적 평가를 가능하게 하는 파이프라인을 구축하였다.
3. 스펙클 분산의 유용성: 본 연구는 OCT 신호의 시간적 변동이 구조적 이미징에 비해 치료 유도 이질성에 관한 보완적인 정보를 제공함을 보여줌으로써 스펙클 분산 대비(contrast)의 가치를 강조한다.
4. 재현성 분석: 선택된 라디오믹스 특징들이 반복된 획득 시 낮은 변동성을 보임을 보여줌으로써, 관찰된 코호트 수준의 변동성이 측정 노이즈가 아닌 생물학적 차이를 반영함을 입증하였다.

결과

• 성능 비교: OCT 전용 모델은 모든 회귀 모델에서 기존의 명시야 기반 형태학 베이스라인보다 일관되게 우수한 성능을 보였다. OCT 전용 모델의 중앙 $R^2$ 값은 0.77에서 0.85 사이였던 반면, BF 전용 모델은 0.61에서 0.69 사이였다 ($p < 0.001$로 모든 비교에서 유의함).
• 다중 모달 융합의 영향: 명시야 특징을 OCT에 추가하는 것이 일관된 성능 이점을 제공하지는 않았다. 결합된 OCT+BF 모델의 중앙 $R^2$ 값은 0.80에서 0.84 사이였다. Random Forest 모델에서만 OCT 전용 모델 대비 통계적으로 유의미한 개선($p = 0.026$, 낮은 검정력 0.62)이 나타났으며, 다른 모든 모델에서는 유의미하지 않았다. 특히, LightGBM의 경우 OCT 전용 모델이 결합 모델(0.84)보다 약간 높은 중앙 $R^2$(0.85)를 달성하였다.
• 특징 견고성: 4개의 핵심 특징(Brightfield Area, Intensity GLRLM GrayLevelNonUniformity, Shape Sphericity, Variance GLDM GrayLevelNonUniformity)이 일관되게 높은 순위를 차지했다. 반복성 테스트 결과, 6회의 반복 스캔에 대해 낮은 상대적 분산(IQR/median = 0.032–0.047)을 보였으나, 코호트 수준의 변동성은 이보다 몇 자릿수 더 높게 나타나 특징의 변별력을 확인하였다.
• 시각적 상관관계: OCT en face 투영의 질적 분석 결과, 높은 선량(8–12 Gy)은 경계 불규칙성을 증가시키고 매끄러운 동심원 형태의 SV 패턴에서 이질적이고 패치 형태의 분포로의 변화를 유도하였으며, 이는 선택된 텍스처 특징들에 의해 포착되었다.

의의 및 주장
본 논문은 OCT 라디오믹스가 종양 구형체의 방사선 선량 반응을 평가하기 위한 민감하고 재현 가능하며 라벨이 필요 없는 고처리량 판독법을 제공한다고 주장한다. 주요 결론은 OCT 전용 분석이 현재의 명시야 기반 형태학보다 우수하다는 것이며, 이는 2D 투영이 놓치는 심도 분해 구조 및 동적 정보가 상당한 선량 관련 신호를 인코딩하고 있음을 시사한다.

저자들은 명시야의 크기 역학(size dynamics)과 OCT의 텍스처 및 형태 특징을 결합한 작은 규모의 공유 예측자가 나타나기는 했으나, 명시야 데이터의 추가가 일반적으로 OCT 단독 사용보다 실질적인 이점을 일관되게 제공하지 않는다고 완곡하게 주장한다. 이는 OCT 단독 파이프라인의 타당성을 뒷받침하며, 동시적인 명시야 획득이나 수동 주석(annotation)의 필요성을 제거한다. 연구의 한계로는 단일 세포주 사용, 단면 설계(진정한 의미의 구형체 내 종단적 추적 불가), 그리고 선량 상관관계가 생물학적 종단 지표의 간접적 대리물이라는 점을 인정하였다. 저자들은 향후 연구가 진정한 종단적 OCT 획득을 가능하게 하고, 다양한 세포주에 걸쳐 검증하며, 라디오믹스 시그니처를 직교하는 생물학적 종단 지표와 연결하는 데 집중해야 한다고 제안한다.