Immunogenomic Intertumor Heterogeneity Across Primary and Metastatic Sites in a Patient With Lung Adenocarcinoma
한 줄 요약
이 연구는 치료 전 상태(treatment-naïve)인 32세 비흡연 여성 lung adenocarcinoma (LUAD) 환자의 사후 검안 표본을 통해, primary tumor와 19개의 synchronous metastases 간 immunogenomic intertumor heterogeneity를 분석했다. Whole-exome sequencing (WES), DNA methylation profiling, 그리고 T-cell receptor (TCR) sequencing을 병행하여 genomic, epigenomic, immune landscape의 공간적 이질성을 규명했다. 주요 발견으로, 모든 종양에서 공통적으로 존재하는 clonal TP53 mutation은 초기 사건임을 시사하지만, 다른 somatic mutations과 neoantigens는 metastatic sites 간에 현저한 이질성을 보였다. 특히 left-side thoracic tumors에서는 특정 prevalent T cell clonotypes (예: CACRPGNEAFF)이 완전히 배제되는 spatial restriction 현상이 관찰되었다. 또한, methylated neoantigen coding mutations의 비율은 T-cell density, richness, clonality와 음의 상관관계를 보여, neoantigen methylation이 immunosuppression을 유도하는 기전일 가능성을 제시했다. 이 연구는 치료 압력이 없는 상태에서 T cell repertoire가 synchronous metastases의 진화와 genomic landscape를 어떻게 형성하는지를 보여주는 중요한 사례이다.
방법
Study design과 대상 환자: 본 연구는 사후 검안(autopsy)을 통해 수집된 32세 여성 never-smoker LUAD 환자의 표본을 대상으로 한 단일 사례 연구(single-case study)이다. 환자는 급성 venous thromboembolism으로 입원 후, CT와 MRI를 통해 광범위한 metastasis가 확인되었으며, liver biopsy로 poorly differentiated LUAD로 진단되었다. 환자는 palliative radiation therapy (C5-C7 spine metastases, 10×30 cGy)와 carboplatin/paclitaxel chemotherapy 1회만 시행한 후 13일 만에 사망했으며, 사후 검안 시 primary tumor (left lung)를 포함하여 thyroid gland, pleural cavity, lymph nodes, heart, GI tract, liver, vertebrae, adrenal gland, kidneys, ovaries 등 총 20개의 tumor samples과 1개의 normal GI tissue sample을 수집했다.
Genomic 및 Epigenomic 분석 프로토콜: DNA는 FFPE tissues에서 AllPrep DNA/RNA FFPE Kit (Qiagen)를 사용하여 추출했으며, Nano Drop 1000 Spectrophotometer로 정량하고 Agilent 2200 Tape Station으로 fragment size를 평가했다. Whole-exome sequencing (WES)은 Agilent SureSelect Human All Exon V4 kit을 이용해 library prep 후 Illumina HiSeq 2000 플랫폼에서 paired-end multiplex sequencing 수행했으며, 평균 coverage depth는 178x (range: 63x-225x)였다. Somatic single nucleotide variants (SNVs)는 MuTect v1.1.4, VarScan 2, Strelka2를 사용하여 호출(call)한 후, sequencing depth ≥ 20× (tumor), ≥ 10× (germline), variant allele frequency (VAF) ≥ 0.02 (tumor), < 0.01 (germline), supporting reads ≥ 4, 그리고 population databases (ESP6500, 1000 genome, EXAC)에서 frequency < 0.01인 변이만 필터링하여 최종 목록을 작성했다. Phylogenetic tree는 R package phangorn v2.0.2를 사용하여 binary presence/absence matrix로 생성했다. Neoantigen 예측은 NetMHC3.4 알고리즘을 통해 patient-specific HLA-A, -B, -C alleles에 대한 9-mer 및 10-mer peptides의 binding affinity (IC50 < 500 nM for binders, IC50 < 50 nM for high affinity)를 평가했다. DNA methylation profiling은 EZ DNA Methylation Kit (Zymo Research)로 bisulfite conversion 후 Infinium Human Methylation 450k arrays (Illumina)에 hybridization 수행했으며, SWAN 알고리즘으로 normalization하고 MethylCIBERSORT를 사용하여 cellular deconvolution을 수행했다.
Immune profiling 및 통계 분석: T-cell repertoire 분석은 immunoSEQ (Adaptive Biotechnologies, hsTCR β Kit) 프로토콜을 따라 TCR β chain의 CDR3 regions을 amplification 후 Illumina MiSeq 플랫폼에서 150 cycle paired-end sequencing 수행했다. T-cell density는 TCR-β template counts를 housekeeping genes의 PCR amplification으로 결정된 usable DNA 총량으로 정규화하여 계산했으며, richness는 unique T-cell rearrangements 수로, clonality는 $1 - Pielou's evenness$ ($Clonality = 1 - \frac{-\sum p_i \ln p_i}{\ln N}$)로 정의했다. 두 샘플 간 T cell repertoire 유사도는 Jaccard index (JI)로 계산했으며, 연속 변수 간 연관성은 Spearman's rank correlation, 범주형 변수는 Kruskal-Wallis test를 사용했고, p-value < 0.05를 통계적 유의성으로 판단했다.
주요 결과
Genomic heterogeneity와 Phylogenetic analysis: 총 228개의 non-silent mutations이 검출되었으며, 샘플당 평균 76개 (range: 57-98)였다. 20개 샘플 모두에서 공유되는 mutation은 단 10개였으며, 이 중 TP53 mutation (chromosome 17_7578382, stop gain, G->C spectrum)만이 모든 종양에 존재하여 early genomic event임을 시사했다. 나머지 cancer gene mutations (CDKN2A, ASXL1, MET 등)은 subclonal diversification 이후의 late events로 나타났다. Phylogenetic tree 분석 결과, 종양 간 진화적 경로가 현저히 달랐으며, Jaccard index (JI)는 0.14에서 0.82 사이 (average=0.49)로 분포했다. Primary tumor와 metastases 간 homology는 0.14~0.73 (average=0.33)였으며, left pleural metastasis, left hilar lymph node, right pleural metastasis 등 thoracic lesions가 primary tumor와 가장 높은 유전적 유사성을 보였다.
T cell repertoire의 공간적 이질성과 Neoantigen 연관성: T-cell density는 3%~38% (average=13%), richness는 4,168~23,487 unique rearrangements (average=14,344), clonality는 0.02~0.05 (average=0.04)로 측정되었다. Non-thoracic tumors는 thoracic tumors보다 T-cell density (p<0.01), richness (p<0.0001), clonality (p<0.01)가 유의하게 높았다. 특히 left-side thoracic tumors (left lung, left pleura, left hilar LN 등)에서는 "other" tumors에서 가장 흔한 clonotype인 CACRPGNEAFF가 완전히 검출되지 않는 현상이 관찰되었다. T cell repertoire 간 JI는 평균 0.35 (range: 0.12-0.61)였으며, 이는 이전 보고된 multi-region NSCLC cohort (average JI=0.17)보다 유의하게 높았다 (p<0.0001). Thoracic tumors는 non-thoracic tumors보다 shared T cell clones의 비율과 빈도가 더 높았다 (p<0.01).
Neoantigen methylation과 Immunosuppression 기전: 평균 39개의 predicted neoantigens (IC50 < 500 nmol/L)가 검출되었으며, high binding affinity neoantigens (IC50 < 50 nmol/L)는 평균 11개였다. T cell repertoire와 neoantigen homology 간에는 약하지만 통계적으로 유의한 양의 상관관계 (r=0.12, p=0.0162)가 있었다. 핵심적으로, methylated neoantigen coding mutations의 비율은 T-cell density (r=-0.46, p=0.0549), richness (r=-0.55, p=0.0152), clonality (r=-0.61, p=0.0055)와 음의 상관관계를 보였다. MethylCIBERSORT deconvolution 결과, CD8+ T cell fraction은 richness (r=0.75, p=0.0002) 및 clonality (r=0.43, p=0.063)와 양의 상관관계를 보였으며, CD8 to Treg ratio도 richness (r=0.68, p=0.0012) 및 clonality (r=0.64, p=0.032)와 양의 상관관계를 보였다. 또한 CD8+ T cell fraction은 methylated neoantigen coding mutations와 음의 상관관계 (r=-0.48, p=0.0391)를 나타냈다.
강점
이 연구는 치료 전(treatment-naïve) 상태에서 primary tumor와 다발성 synchronous metastases를 동시에 분석할 수 있는 희귀한 autopsy case를 활용하여, 치료 요법에 의한 선택 압력(selective pressure) 없이 intrinsic intertumor heterogeneity를 규명했다는 점에서 방법론적으로 강점이 있다. WES, DNA methylation, TCR sequencing을 통합하여 genomic, epigenomic, immune landscape의 상호작용을 다각도로 분석했으며, 특히 neoantigen methylation과 T cell infiltration 간의 정량적 상관관계를 제시함으로써 immunosuppression의 새로운 기전을 제안했다. 또한, 기존 연구들이 주로 longitudinal changes나 단일 종양 내 intratumor heterogeneity에 집중했던 것과 달리, spatial distribution (thoracic vs non-thoracic)에 따른 T cell repertoire의 차이를 명확히 구분하여 보고했다.
한계
본 연구는 단일 환자(single patient)를 대상으로 한 사례 연구이므로, 결과를 일반적인 LUAD population으로 일반화(generalize)하기에는 표본 규모가 부족하다. 사후 검안 표본은 FFPE 처리 과정과 necrosis 등으로 인해 DNA quality에 영향을 받을 수 있으며, 이는 sequencing depth와 mutation calling의 정확도에 잠재적 편향을 초래할 수 있다. 또한, T cell repertoire 분석이 TCR β chain에만 국한되어 전체 immune response를 완전히 반영하지 못할 가능성이 있으며, neoantigen 예측은 in silico 알고리즘(NetMHC3.4)에 의존하여 실제 antigen presentation과 기능적 반응과의 일치성을 검증하지 못했다. 마지막으로, methylated neoantigen coding mutations와 T cell metrics 간의 상관관계는 통계적으로 유의미하지만 인과관계를 입증하기에는 observational data의 한계가 존재한다.
해석
이 연구는 lung adenocarcinoma에서 synchronous metastases 간에 genomic 및 immunogenomic heterogeneity가 광범위하게 존재하며, 이는 T cell repertoire의 공간적 분포와 밀접하게 연관되어 있음을 보여준다. 특히 left-side thoracic tumors에서의 특정 T cell clonotype 배제 현상은, 종양의 anatomical location이 immune surveillance escape에 중요한 역할을 할 수 있음을 시사한다. Neoantigen methylation이 T cell density 및 clonality를 감소시키는 기전으로 작용할 가능성은, epigenetic modification을 표적으로 한 immunotherapy combination 전략의 이론적 근거를 제공한다. LLM Wiki의 oncology 문헌들과 비교할 때, 이 연구는 treatment-naïve 상태에서의 baseline heterogeneity를 제공하여, 기존 치료 후 변화된 landscape 연구들과 대비되는 중요한 참조점이 된다. 향후 multi-omics integration을 통한 precision oncology 접근에서, metastatic site-specific immune profiling의 중요성을 강조하는 근거로 활용될 수 있다.