비유전성 진화의 파동 접근법: 뮬러의 톱니와 적응 속도

초록

본 논문은 무성생식 집단에서 뮬러의 톱니(Muller’s ratchet)와 양성 선택에 의한 적응 속도를 설명하기 위해 ‘여행파동(traveling‑wave)’ 이론을 적용한다. 집단의 대다수는 결정론적 방정식으로, 가장 적합한 유전형(스톱스틱 엣지)은 확률론적으로 다루어 전체 적합도 변화의 핵심을 포착한다. 저자는 스톱스틱 엣지의 새로운 모델링 방법과 이산 적합도 계층을 연속체로 근사할 때 발생하는 오차를 보정하는 해석적 교정항을 제시한다. 결과적으로 뮬러의 톱니 속도는 높은 정확도로 예측되고, 적응 속도는 넓은 파라미터 영역에서 로그‑형태로 인구 규모에 따라 증가한다는 결론을 얻는다.

상세 분석

여행파동 이론은 무성생식 집단의 진화를 두 부분으로 분리한다. 첫 번째는 대다수 개체가 차지하는 ‘대량부(bulk)’로, 이 부분은 평균 돌연변이율과 선택압을 이용한 연속적인 적합도 분포 방정식으로 기술된다. 두 번째는 가장 높은 적합도를 가진 ‘최고 적합도 계층(high‑fitness class)’이며, 이 계층은 개체 수가 매우 적어 확률적 요인이 지배한다. 기존 연구에서는 이 계층을 단순히 ‘정규분포의 꼬리’로 근사했지만, 저자는 ‘스톱스틱 엣지(stochastic edge)’라는 개념을 도입해 개별 계층의 사라짐과 재생성을 명시적으로 모델링한다.

뮬러의 톱니 현상은 무성생식 집단에서 무해한 돌연변이가 축적되어 평균 적합도가 감소하는 과정이다. 여행파동 접근법에 따르면, 톱니가 한 번 ‘클릭’될 때는 최고 적합도 계층이 완전히 사라지고, 그 다음 계층이 새로운 엣지가 된다. 저자는 이 과정의 평균 시간 간격을 ‘클릭 시간’이라고 정의하고, 확률적 소멸률과 결정론적 파동 전파 속도를 결합해 정확한 식을 도출한다. 특히, 이산 적합도 계층을 연속적인 파동으로 근사할 때 발생하는 ‘계층 간격 오류’를 보정하기 위해, 파동 전단(gradient)과 전파 속도에 대한 1차 및 2차 교정항을 추가하였다. 이러한 교정은 작은 돌연변이 효과(s)와 큰 집단 규모(N)에서도 오차를 10⁻² 이하로 낮춘다.

양성 선택 하에서의 적응 속도는 최고 적합도 계층이 지속적으로 ‘앞당겨지는’ 형태로 나타난다. 여기서도 스톱스틱 엣지의 확률적 부양이 핵심이다. 저자는 적합도 파동이 일정한 형태를 유지하면서 평균 적합도가 선형이 아닌 로그 형태로 증가한다는 점을 보였다. 구체적으로, 적응 속도 v는
v ≈ (2s ln N / ln ( s/U ))
와 같은 식으로 근사되며, 여기서 s는 선택계수, U는 전체 돌연변이율이다. 이 식은 N이 충분히 클 때는 로그‑성장으로 전환되며, N이 극도로 큰 경우에만 선형 성장으로 변한다는 예측을 제공한다.

시뮬레이션 결과는 이론적 예측과 매우 높은 일치도를 보였다. 뮬러의 톱니에서는 클릭 간격이 이론값과 5 % 이내, 적응 속도에서는 로그‑성장 구간에서 10 % 이내의 차이만을 보였다. 특히, 교정항을 포함하지 않은 기존 모델은 작은 N·U 영역에서 과도하게 빠른 적응을 예측했으나, 본 논문의 보정이 이를 정확히 잡아낸다.

이러한 결과는 무성생식 미생물, 바이러스, 그리고 클론성 암세포 집단 등에서 관찰되는 급격한 적합도 변동을 이해하는 데 중요한 통찰을 제공한다. 특히, ‘스톱스틱 엣지’ 개념은 유전적 다양성이 제한된 상황에서도 집단이 어떻게 새로운 적합도 수준을 탐색하는지를 설명한다는 점에서 진화 이론에 새로운 시각을 제시한다.