배경

암의 진행은 유전적으로 다른 암세포 집단이 선택되는 과정인 동시에, 저산소증·면역억제·섬유화가 축적되는 생태계 변화다. 그러나 단일세포 RNA 시퀀싱은 조직을 해리하므로 세포의 원래 위치가 사라지고, 다중영역 시퀀싱은 제한된 지점만 채취해 세포 상태와 계통의 연속적인 공간 구조를 놓치기 쉽다. 암세포 아클론이 어디에서 자라며 주변 세포를 어떻게 바꾸는지, 전이 능력을 얻은 세포가 원발 종양의 어느 구역에서 출현하는지는 제대로 규명되지 않았다.

연구진은 이러한 간극을 메우기 위해 고해상도 공간전사체와 시간에 따라 유전 표지가 누적되는 계보추적을 결합했다. 대상은 Kras 활성화와 Trp53 결손으로 폐선암이 발생하는 KP-Tracer 생쥐다. 8∼12주령 생쥐의 폐포 2형 세포에 Cre 바이러스를 투여해 종양 발생과 CRISPR–Cas9 계보 기록을 동시에 시작했다. 암세포의 현재 전사 상태뿐 아니라 공통 조상과 증식 이력까지 조직 좌표 위에서 복원하려는 설계다.

핵심 발견

연구진은 넓은 조직 면적을 준세포 수준으로 읽는 슬라이드 시퀀싱(Slide-seq)과 단일 핵 수준의 민감도를 제공하는 슬라이드 태그(Slide-tags)를 함께 사용했다. 총 49개 공간전사체 어레이에서 100개가 넘는 종양을 분석했으며, 계보 표지가 누락된 위치에는 주변 암세포의 공간·유전 정보를 이용한 보정 알고리즘을 적용했다. 실제 데이터에서 가려 둔 표지를 복원한 중앙 정확도는 90%였고, 반복 보정으로 데이터셋별 누락 정보의 4∼58%, 평균 31%를 회수했다. 무작위 예측의 중앙 정확도는 67%였다.

종양이 진행되면서 폐포형 상태를 유지한 암세포는 주로 가장자리에, 상피간엽전이(EMT)와 저산소 유전자 프로그램을 띤 세포는 내부에 자리했다. 빠르게 확장한 고적합도 아클론 주변에는 Arg1 양성 면역억제성 종양연관 대식세포와 Postn 양성 근섬유아세포형 암연관 섬유아세포가 모였다. 계통수에서 계산한 적합도는 기존 전사체 기반 적합도 지표와 피어슨 상관계수 0.4를 보였으며, 저산소·EMT 공간 공동체가 가장 높은 적합도와 결부됐다. 이는 아클론 증식, 산소 고갈, 면역·기질 재편, 전이 촉진성 세포 상태가 순차적으로 맞물린다는 설명을 뒷받침한다. 공간 분석은 Piezo2·Robo1·Pecam1을 발현하면서 Nkx2-1을 유지하고 Vim은 발현하지 않는 미보고 암세포 상태도 찾아냈다.

전이가 광범위하게 일어난 생쥐에서는 폐를 200∼500마이크로미터 간격으로 절단해 원발 종양을 입체적으로 추적했다. 종격동 림프절·늑골·횡격막의 전이 병소는 모두 원발 종양 T2 안의 공간적으로 제한된 아클론과 계통적으로 연결됐다. 전이 병소의 TGF-β 프로그램은 해당 원발 아클론과 거의 같았지만(log2 배수변화 −0.14, P=1.0), 콜라겐 관련 발현은 크게 증가했다(log2 배수변화 3.81, P<10⁻⁵). COL3A1 단백질 염색도 전이 부위의 섬유화를 뒷받침했다. Nature Genetics 논문

의미와 전망

이번 연구는 공격적인 암세포 상태가 세포 내부 변이만으로 고정되는 것이 아니라, 아클론 확장이 만든 저산소·면역억제·섬유화 환경과 함께 형성된다는 공간적 근거를 제시했다. 전이의 씨앗이 원발 종양 전체에 고르게 퍼진 것이 아니라 특정 국소 생태계에 집중됐다는 결과도 주목된다. 전이 후에는 기존 TGF-β·EMT 프로그램이 유지되는 동시에 콜라겐 침착이 더해져, 전이 정착에 맞는 새로운 환경이 조성됐다.

다만 결과는 Kras;Trp53 기반 생쥐 폐선암 모델에 주로 의존한다. 얇은 조직 절편은 종양의 전체 3차원 구조를 대표하지 못할 수 있고, 공간 자료의 계보 표지 다양성도 해리형 단일세포 분석보다 낮다. 누락 표지를 이웃 세포로 보정하는 방식은 이동성이 큰 암세포에서 오류가 생길 여지도 있다. 사람 폐선암 코호트와 치료 전후 검체에서 같은 공간 구조가 재현되는지, 저산소증이나 섬유화 신호를 차단했을 때 아클론의 전이 능력이 실제로 약해지는지 검증해야 한다.