MSI-stratified analysis of TREM2 and HAVCR2 expression in mucinous colorectal adenocarcinoma: a structured analytical workflow with signature de-circularization

기술 요약: 점액성 대장암(Mucinous Colorectal Adenocarcinoma) 내 TREM2 및 HAVCR2의 MSI-층화 분석

문제 제기
TREM2(Triggering Receptor Expressed on Myeloid cells 2)와 HAVCR2(TIM-3 인코딩)는 대장암(CRC)에서 부상하는 면역 항암 표적이다. 대규모 전사체 연구를 통해 이들의 벌크 전사체 발현 사이의 강력한 양의 상관관계가 관찰되었으나, 이러한 연관성이 점액성 선암종(MC)의 특이적인 생물학적 특징을 반영하는 것인지, 아니면 주로 MC 아형 내에 높은 미세부수체 불안정성-높음(MSI-H) 상태의 유병률에 의해 주도되는 것인지는 여전히 해결되지 않은 과제로 남아 있다. MC 종양은 비점액성 선암종(NMC)에 비해 MSI-H가 풍부하며(20–40% vs. 10–15%), 이는 혼란을 야기하는 시나리오를 만든다. 즉, 점액성 조직학적 특성에 기인한 것으로 간주되는 분자적 차이가 실제로는 MSI 조성에서 비롯될 수 있다는 것이다. 또한, MC와 NMC를 비교하는 기존의 벌크 전사체 연구들은 흔히 MSI 상태를 조정하지 못하여 가짜 연관성(spurious associations)을 초래할 위험이 있다. 더불어, 이 분야의 유전자 시그니처 분석은 관심 유전자가 해당 시그니처 측정에 의도치 않게 포함되는 "순환성(circularity)" 문제로 인해 자주 어려움을 겪는다.

방법론
본 연구는 TCGA-COAD 및 TCGA-READ 코호트에서 추출한 251명의 고유 환자 데이터(MC 32명, NMC 219명; MSI-H 35명, non-MSI-H 216명)를 사용하여 구조화된 분석 워크플로우를 적용하였다. 분석은 다음과 같은 엄격한 단계로 진행되었다:

1. MSI-층화 상관관계 분석: 전체 코호트와 조직학적 아형(MC vs. NMC), MSI 상태(MSI-H vs. MSS), 그리고 해부학적 위치(COAD vs. READ)별로 층화하여 TREM2와 HAVCR2 사이의 피어슨(Pearson) 및 스피어만(Spearman) 상관관계를 계산하였다. 하위 그룹 간의 상관계수를 비교하기 위해 Fisher's z-변환을 사용하였다.
2. 다변량 모델링: MC × MSI 상호작용 항을 테스트하기 위해 선형 모델을 적합시켜 병리적 단계(pathological stage)와 해부학적 위치를 조정하였으며, 이를 통해 MSI와 유전자 발현 사이의 관계를 점액성 조직학이 수정하는지 확인하였다. MSI, 위치, 단계를 순차적으로 조정하여 편상관관계(partial correlations)를 산출하였다.
3. 유전자 시그니처 탈-순환화 감사(De-circularization Audit): 7개의 큐레이션된 유전자 시그니처(TREM2+ 대식세포, Goblet/점액성, MUC1 크로스토크 포함)에 대해 체계적인 감사를 수행하였다. 자기 포함 편향(예: TREM2가 결과값일 때 TREM2+ 대식세포 시그니처에서 TREM2 제외)과 교차 세트 공선성(예: 중복 세트에서 ERBB2 및 GRB2 제거)을 방지하기 위해 중복되는 유전자를 제거하였다.
4. 혼란 요인 입증: 연구는 EREG를 대조 사례로 활용하여, MSI 층화를 통해 어떻게 MSI 구성이 연관성의 방향을 역전시킬 수 있는지 보여주었다.
5. 외부 검증 및 MR: GSE119409 직장암 코호트(n=56, 반응 데이터 포함)를 사용하여 화학방사선요법(CRT) 반응과의 연관성을 맥락화하였고, 혈액 유래 eQTL 도구를 사용한 탐색적 멘델 무작위 배정(Mendelian Randomization, MR)을 통해 인과적 추론을 평가하였다.

주요 결과

• 상관관계 강도 및 특이성: TREM2와 HAVCR2는 전체 코호트에서 강한 양의 상관관계를 보였다(Pearson r = 0.70, p = 2.8e-38). 이 상관관계는 non-MSI-H 종양(r = 0.66)보다 MSI-H 종양(r = 0.83)에서 더 강하게 나타났다(p = 0.041). 결정적으로, 이 상관관계는 MC(r = 0.80)와 NMC(r = 0.67) 사이에 유의미한 차이를 보이지 않았다(p = 0.17).
• 상호작용 테스트: 공식적인 상호작용 모델을 통해 MC 상태가 TREM2 또는 HAVCR2 발현에 미치는 MSI의 효과를 수정하지 않음을 확인하였다(TREM2에 대한 상호작용 p = 0.76; HAVCR2에 대한 상호작용 p = 0.50). MSI-H의 주효과는 두 유전자 모두에서 유의미했으나, MC의 주효과는 유의미하지 않았다.
• 혼란 증거: EREG 유전자는 MSI 층화에 따라 연관성 방향이 완전히 역전됨을 보여주었다. 미조정 분석에서 EREG는 MC에서 낮게 나타났다(Cohen d = -0.51). 그러나 MSI-H 층 내에서 EREG는 MC에서 높게 나타났으며(d = +0.87), 이는 MC 내의 높은 MSI-H 유병률이 어떻게 오도된 한계 연관성(marginal associations)을 생성할 수 있는지를 보여준다.
• 시그니처 분석: 탈-순환화 및 Benjamini-Hochberg 교정을 거친 후, Goblet/점액성 시그니처(d = +0.88)와 TREM2+ 대식세포 시그니처(d = +0.63)는 MC와 NMC 사이에 유의미한 차이를 보였다. 그러나 TREM2+ 대식세포의 상승은 MSI에 대해 조정되지 않았다.
• 외부 검증: GSE119409 코호트에서 TREM2+ 대식세포 시그니처는 CRT 민감도와 유의미한 연관성을 보였으나(P_BH < 0.034), 해당 데이터셋에는 MC 조직학적 주석이 부족했다. 단일 유전자 TREM2 및 HAVCR2는 완만한 예측 가치를 보였다(AUC ~0.68–0.69).
• 멘델 무작위 배정(MR): MR 분석에서 Bonferroni 유의성을 충족하는 유전자는 없었으며, 공통 국소화(colocalization) 분석에서도 혈액 유래 eQTL과 CRC GWAS 신호 사이의 공유된 인과적 변리에 대한 증거를 찾지 못했다.

의의 및 주장
본 논문은 TREM2와 HAVCR2 사이의 강력한 벌크 전사체 상관관계가 점액성 조직학에 특이적인 것이 아니라, 주로 MSI-H 상태와 관련이 있다고 주장한다. 본 연구는 MSI를 적절히 층화했을 때 "MC 특이적 TREM2-TIM-3 면역 축"이라는 가설이 가용한 벌크 전사체 데이터에 의해 지지되지 않음을 논한다.

저자는 세 가지 방법론적 기여를 강조한다:

1. 분석 프레임십: 공통적인 분석 오류를 피하기 위해 체계적인 유전자 세트 탈-순환화를 포함하는, MSI-층화 면역 전사체 분석을 위한 재현 가능한 워크플로우를 제공한다.
2. 혼란 정량화: EREG 사례 연구를 통해, 조직학적 아형 간의 불균등한 MSI-H 분포가 어떻게 가짜 연관성을 생성할 수 있는지 정량적으로 입증하고, 한계 비교(marginal comparisons)와 함께 층별 보고(stratum-specific reporting)를 권고한다.
3. 부정적 결과 보고: 연구는 MC 중심의 TREM2/TIM-3 실험에 대한 중복적인 전임상 투자를 방지하기 위해 명시적으로 부정적인 일차 결과를 보고하며, 자원을 MSI-층화 조사나 다른 측면의 MC 생물학(예: 점액 생산 경로)으로 돌리는 것이 더 나을 수 있음을 시사한다.

저자는 작은 MC 샘플 크기(n=32), 벌크 RNA-seq이 세포 유형별 패턴을 해결할 수 없는 점, 그리고 MR 분석에서의 조직 특이적 도구 부족 등의 한계를 언급하며 자신의 결과에 대해 겸허한 태도를 유지한다. 결론적으로, TREM2-HAVCR2 축은 CRC에서 생물학적으로 유의미하지만, 그 동인은 점액성 분화 자체보다는 MSI-H의 면역원성 있는 성격인 것으로 보인다.