최소 적재 클래스 손실과 뮤터레이션 라쳇 억제

초록

본 연구는 비재조합 무성생식 집단에서 유해 돌연변이가 축적되는 뮤터레이션 라쳇 현상을 이론적으로 분석한다. 확산 이론을 이용해 비에피스타시스와 시너지적 에피스타시스가 존재할 때 라쳇 클릭 간 평균 시간(T_J)을 구하고, 비에피스타시스 경우는 최소 적재 클래스 크기와 선택계수(s)에 의해 속도가 결정됨을 보인다. 반면 시너지적 에피스타시스(α>1)에서는 T_J가 급격히 증가해 라쳇 속도가 실질적으로 0에 수렴한다. 분석 결과는 수치 시뮬레이션과 일치한다.

상세 분석

이 논문은 무성생식, 비재조합 haploid 집단이 유해 돌연변이를 되돌릴 수 없이 축적하면서 발생하는 뮤터레이션 라쳇(Muller’s ratchet)의 동역학을 정량적으로 규명한다. 저자는 먼저 개체의 적재 수 k에 대한 적합도 함수를 W(k)= (1‑s)^{k^{α}} 로 정의한다. α=1이면 비에피스타시스(선형 적합도)이며, α>1이면 시너지적 에피스타시스(협동적 악화)를 의미한다. 이 모델은 무한 유전체 가정 하에 백뮤테이션을 무시하고, 변이율 U와 선택계수 s가 주어졌을 때 각 적재 클래스의 평균 빈도 X_J(k) 를 구한다.

정상 상태 해를 구하기 위해 확산 근사와 퀘이시스페시스 방정식을 사용한다. 비에피스타시스 경우, 기존 연구와 일치하게 X_J(k)=e^{-U/s}(U/s)^k/k! 로, 최소 적재 클래스의 평균 크기 n_0 = N e^{-U/s} 가 라쳇 속도를 결정한다는 것을 재확인한다. 특히, 클릭 간 평균 시간 T_J는 n_0와 s의 함수이며, n_0가 클수록 라쳇이 느려진다.

시너지적 에피스타시스(α>1)에서는 적합도 곡선이 급격히 가팔라져, 적재 클래스 J가 증가함에 따라 선택압이 강화된다. 저자는 α=2인 경우에 대해 정확한 해를 구하고, Bessel 함수 I_{2J} 를 이용해 X_J와 X_J(k)를 표현한다. 일반 α>1에 대해서는 대략적인 해를 도출해 X_J ≈ exp