[논문 리뷰] Identifying expanding TCR clonotypes with a longitudinal Bayesian mixture model and their associations with cancer patient prognosis, metastasis-directed therapy, and VJ gene enrichment

David Swanson, Alexander D. Sherry외 3명|arXiv (Cornell University)|2026. 01. 08.

T-cell and B-cell Immunology인용 수 0

한 줄 요약

이 논문은 시간에 따라 확장하는, 수축하는, 및 정적 TCR 클론타입을 식별하기 위한 종단 Bayesian 혼합 모델을 도입하고, 동적 클론타입을 MDT 및 예후와 연결하며, 페널라이즈드 로그-선형 모델링을 사용하여 VJ 유전자 강화를 연구합니다.

ABSTRACT

Examination of T-cell receptor (TCR) clonality has become a way of understanding immunologic response to cancer and its interventions in recent years. An aspect of these analyses is determining which receptors expand or contract statistically significantly as a function of an exogenous perturbation such as therapeutic intervention. We characterize the commonly used Fisher's exact test approach for such analyses and propose an alternative formulation that does not necessitate pairwise, within-patient comparisons. We develop this flexible Bayesian longitudinal mixture model that accommodates variable length patient followup and handles missingness where present, not omitting data in estimation because of structural practicalities. Once clones are partitioned by the model into dynamic (expanding or contracting) and static categories, one can associate their counts or other characteristics with disease state, interventions, baseline biomarkers, and patient prognosis. We apply these developments to a cohort of prostate cancer patients who underwent randomized metastasis-directed therapy or not. Our analyses reveal a significant increase in clonal expansions among MDT patients and their association with later progressions both independent and within strata of MDT. Analysis of receptor motifs and VJ gene enrichment combinations using a high-dimensional penalized log-linear model we develop also suggests distinct biological characteristics of expanding clones, with and without inducement by MDT.

연구 동기 및 목표

• 암에서 TCR 클로나타입 다이나믹스(확장/수축 vs 정적) 분석을 촉진한다.
• 가변 추적 및 누락 데이터를 수용하기 위한 2-구성요소 종단 Bayesian 혼합 모델을 개발한다.
• 동적 클로나타입을 임상 결과, 기저 바이오마커, 치료(MDT)와 연결한다.
• 페널라이즈드 완전 포화 로그 선형 프레임워크를 사용하여 의미 있는 상호작용을 식별하기 위해 VJ 유전자 강화 패턴을 평가한다.

제안 방법

• 클론 카운트의 동적(lambda_ijk) 및 정적(lambda_ij) 비율을 갖는 2-구성요소 포아송 기반 혼합 모델을 제안하고, 감마 사전분포와 환자 간 하이퍼파라미터의 계층적 공유를 도입한다.
• Stan에서 Hamiltonian Monte Carlo로 모델을 적합시켜 각 클론의 동적 소속에 대한 후방 확률을 얻는다.
• 베이지안 혼합 접근법과 Fisher의 정확 검정/ Beta-이항 방법을 비교하고 MAR 처리 및 유연한 추적 추적의 이점을 강조한다.
• L1 정규화가 적용된 페널라이즈드 완전 포화 로그-선형 모델을 적용하여 Clonotype 확장 상태와 VJ 유전자 강화 상호작용을 테스트하고 교차 검증으로 패널티를 선택한다.

실험 결과

연구 질문

• RQ1 비정형 추적 시간에서 TCR 클로나타입을 확장, 수축 또는 정적으로 신뢰성 있게 분류하는 종단 Bayesian 혼합 모델이 가능한가?
• RQ2 확장 클로나타입이 메타스시스-지향 치료 및 전립선암에서 무진행 생존과 연관이 있는가?
• RQ3 클론 다이나믹스와 기저 바이오마커, MDT와의 상호작용을 포함한 관계는 무엇인가?
• RQ4 클론 다이나믹스를 고려할 때 확장 클로나타입은 특정 VJ 유전자 가족 강화 패턴을 보이는가?
• RQ5 MDT 유도와 비-MDT 유도 확장 클론 간 수용체 모티프나 시퀀스 수준 특성이 차이가 있는가?

주요 결과

• MDT는 기저선-추적 분석에서 확장 클론의 수를 유의하게 증가시킨다.
• 확장 클론의 수는 MDT 구간 및 모델 사양에 따라 변하는 무진행 생존과의 예후 연관성을 보인다.
• 패널라이즈드 로그-선형 모델링은 특정 V 및 J 유전자 가족 조합의 유의미한 강화 및 확장 관련 특정 상호작용(TCRBJ02*TCRBV02..BV29 등)을 나타내는 비무작위 강화 패턴을 시사한다.
• 기저 바이오마커(upr, il15, vegfa 등)가 클론 확장 또는 수축과 연관되며, penalized 회귀 하에서 gmcsf*mip1b, il15*vegfa와 같은 특정 상호작용이 나타난다.
• Lorenz-곡선 분석은 클론 확장이 빈도 변화로 확장되기보다는 번들링되는 경향이 있음을 시사한다.
• MDT 상태는 동적 클론의 수와 예후 연관성 모두에 영향을 미쳐 확장-무진행 생존( PFS ) 연관 해석에서의 교란 요인을 강조한다.