JAK2V617F 양성 골수증식성 신생물의 발생과 진행 실생활 발병률 및 수학 모델링 기반 분석

초록

**
본 연구는 JAK2V617F 변이를 가진 전형적 골수증식성 신생물(다혈구증가증·혈소판증가증)의 발병 시점과 종양 성장 기간을 실생활 진단 연령 데이터를 이용해 수학적으로 모델링하였다. 변이 발생률이 연령에 따라 Gompertz 형태로 급증한다는 가설을 검증하고, 첫 번째 종양 세포가 나타나는 평균 연령을 63.1 ± 13 년, 종양이 진단에 이르기까지의 성장 기간을 8.8 년으로 추정하였다. 또한, 다혈구증가증이 혈소판증가증보다 약 1.5년 늦게 변이가 활성화된다는 차이를 확인하였다.

**

상세 분석

**
이 논문은 골수증식성 신생물(MPN)의 발생 메커니즘을 두 단계(T₁, T₂)로 분리하고, 각각을 확률론적 모델로 기술한다. T₁은 배아기부터 변이가 발생하고 증식 가능한 세포가 되는 시간으로, 변이 발생률 τ와 증식 성공 확률 p의 곱인 활성 변이율 δ=τp를 도입한다. 기존 연구에서 가정된 일정한 변이율 대신, 저자들은 δ가 연령 t에 따라 Gompertz 함수 A·exp(k·t)로 증가한다는 가설을 세우고, 이를 실제 연령별 진단 데이터에 적합시켰다. 이 과정에서 사망률을 보정한 연령별 인구 분포를 사용해 진단 연령의 실제 발생 빈도를 추정하였다.

T₂는 최초 종양 세포가 진단 기준인 10¹² 세포(≈10⁶ 세대)까지 증식하는 기간이다. 저자들은 세포 분열 시간의 변동성이 평균화되어 개인 간 차이가 작다고 보고, T₂를 고정값 α(=8.8 년) 혹은 로그정규 분포로 모델링하였다. 두 모델(A.1: 고정 α, A.2: 로그정규) 모두 EM 알고리즘을 이용해 파라미터(A, k, α 또는 μ, σ²)를 추정했으며, χ² 적합도 검정으로 A.1이 더 간결하면서도 데이터에 충분히 부합함을 확인했다.

결과적으로, 활성 변이율이 연령과 함께 급격히 상승함을 보여주며, 이는 노화와 연관된 DNA 복구 저하, 염증성 미세환경 변화, 줄기세포 활성화 등 생물학적 요인과 일치한다. 변이 발생 평균 연령 63 년은 기존 임상 보고와 부합하지만, 변이가 실제 종양으로 진행되기까지는 약 9년이라는 일정한 기간이 존재함을 강조한다. 또한, 다혈구증가증(PV)과 혈소판증가증(ET) 사이의 차이는 변이 활성화 시점에만 존재하고, 성장 속도는 동일하다는 점에서 두 질환의 병리학적 차이를 정량적으로 설명한다.

이 연구는 (1) 연령 의존적 변이율 모델링, (2) EM 기반 파라미터 추정, (3) 외부 코호트 검증이라는 세 가지 방법론적 강점을 가진다. 다만, 변이 발생을 단일 파라미터 δ에 귀속시키는 단순화와, 진단 기준을 고정된 세포 수치로 가정한 점은 실제 임상에서의 진단 다양성을 충분히 반영하지 못한다는 한계가 있다. 향후에는 개별 환자의 클론형성능력, 동시 발생하는 보조 변이(TET2, ASXL1 등) 및 치료 개입 효과를 포함한 다변량 모델링이 필요할 것이다.

**