다단계 암 발생을 위한 돌연변이 대기시간 모델

논문은 암 발생을 설명하기 위해 세포 집단 내에서 m개의 특정 돌연변이가 순차적으로 축적되는 시간을 분석한다. 모델은 고정된 크기 N의 세포 집단을 대상으로 하며, Moran 과정에 기반한다. 각 세포는 평균 수명 1을 갖고, 사망 후 무작위로 선택된 다른 세포가 복제된다. 돌연변이는 세포가 살아 있는 동안에 발생하며, 유형 j−1 세포가 유형 j 세포로 변하는 속도는 u_j이다. 먼저 m=2인 경우를 다룬다. 유형 1 돌연변이는 전체 집단에서 발생률 N u₁을 갖고, 발생한 유형 1 세포는 분기 과정을 통해 자손을 만든다. 분기 과정의 총 생산량 Z는 안정법칙 지수 ½에 속하며, Z의 꼬리 확률은 P(Z>n)≈(πn)^{-1/2}이다. 따라서 최초의 유형 2 돌연변이가 나타날 때까지 필요한 유형 1 가계 수 M은 O(1/√u₂)이며, 전체 대기시간 τ₂는 (N u₁ √u₂)⁻¹ 스케일의 지수분포에 근접한다. 이 결과는 기존 문헌의 비공식적 추정(식 1.1)을 엄밀히 증명한 것으로, Theorem 1에서 위험 함수 h(t)=1−e^{−2t/λ} (λ>0) 형태를 제시한다. 다음으로 일반 m≥2에 대해 확장한다. 저자들은 복합 비율 r_{j,m}=u₁^{1/2^{j+1}} u₁^{1/4^{j+2}} … u₁^{1/2^{m−j}} 를 정의하고, q_{j,m}≈r_{j,m}가 유형 j 세포가 최종적으로 유형 m 세포를 낳을 확률임을 보인다. 주요 가정은 네 가지이다: (i) N u₁→0 (돌연변이율이 충분히 작다), (ii) 연속적인 돌연변이율 비율 u_{j+1}/u_j이 일정 상수 b_j>0보다 크다, (iii) N^{a} u_m→0 (가장 늦은 단계에서도 돌연변이율이 작다), (iv) N r_{1,m}→∞ (분기 과정 근사가 유효). 이러한 가정 하에 τ_m·N r_{0,m}는 Exp(1)으로 수렴한다 (Theorem 2). 즉, 대기시간은 N과 각 단계의 돌연변이율이 결합된 복합 스케일 N r_{0,m}에 의해 결정된다. 경계 상황인 N r_{1,m}→γ>0에서는 Theorem 3이 적용된다. 여기서는 유형 1 세포가 고정될 확률을 무시할 수 없으므로, 두 단계 모델을 변형해 u₁ τ_m의 극한 분포가 Exp(α) 형태가 된다. α는 γ에 대한 무한급수 α=∑_{k≥1}γ^{k}k!/(k−1)!² 로 정의되며, 이는 Wright–Fisher 확산 과정에서 첫 번째 유형 m 돌연변이가 발생하기 전까지 과정이 0 또는 1에 도달할 확률을 해석함으로써 도출된다. 식 (1.5)와 (1.6)에서 보듯, 확산 과정의 경계값 문제를 풀어 α를 구한다. 이 경우 대기시간은 u₁ τ_m≈Exp(α)이며, 시뮬레이션 결과와 매우 잘 맞는다. 논문의 나머지 부분은 정리와 증명이다. 섹션 2에서는 Theorem 1의 엄밀한 증명을 제공하고, 섹션 3에서는 두 유형 모델에 대한 기본 결과(분기 과정, Wright–Fisher 확산, 위험 함수)를 정리한다. 섹션 4에서는 귀납적 방법을 사용해 r_{j,m}과 q_{j,m}의 관계를 증명하고, 전체 m단계 결과를 도출한다. 마지막으로, N u₁→∞인 경우는 별도 연구(