진화 시뮬레이션에서 종의 등장과 돌연변이율 진화

초록

이 연구는 단세포 무성생식 유기체의 개체군을 단백질 수준의 분자 특성과 상호작용에 연결한 첫 원리 모델을 구축한다. 시뮬레이션 결과, 낮은 돌연변이율에서는 동일한 유전체를 가진 개체군, 즉 ‘종’이 형성되지만 높은 돌연변이율에서는 유전적 다양성이 급격히 확대돼 종이 사라진다. 돌연변이율 자체가 선택 가능한 형질일 경우, 초기에는 높은 돌연변이율이 유지되다 일정 수준 이상의 적합도가 달성되면 낮은 돌연변이율을 가진 개체가 점차 우세해져 현대 DNA 기반 생물과 유사한 돌연변이율에 수렴한다.

상세 분석

본 논문은 진화적 시뮬레이티드 어닐링(simulated annealing) 프레임워크를 이용해, 유전체와 단백질 구조·기능 사이의 상호 의존성을 정량화한 모델을 제시한다. 핵심은 ‘단백질-단백질 상호작용 네트워크’를 적합도 함수에 직접 포함시켜, 단일 유전자의 변이가 전체 네트워크 안정성에 미치는 에피스타시스(epistatic) 효과를 계산한다는 점이다. 시뮬레이션은 크게 두 단계로 진행된다. 첫 번째 단계에서는 돌연변이율을 고정하고, 무작위 변이가 누적되면서 적합도가 ‘펑크’ 형태로 급격히 상승한다. 이때 관찰되는 급증은 특정 변이가 동시에 여러 단백질의 접힘 안정성 및 결합 친화도를 개선함으로써, 네트워크 전체가 새로운 안정적 ‘피크’에 도달하는 현상이다. 두 번째 단계에서는 돌연변이율 자체를 가변 형질로 두어, 돌연변이율이 적합도와 어떻게 상호작용하는지를 탐색한다. 초기에는 높은 돌연변이율이 다양성을 제공해 빠른 적합도 상승을 촉진하지만, 적합도가 일정 수준에 도달하면 낮은 돌연변이율을 가진 개체가 선택 압력에 의해 서서히 우세해진다. 이는 ‘유전적 부하(genetic load)’와 ‘돌연변이 부하(mutational load)’ 사이의 균형을 자연스럽게 맞추는 메커니즘으로 해석될 수 있다.

특히 저자들은 ‘종(species)’을 ‘동일 유전체를 가진 개체군’으로 정의하고, 이 정의 하에서 종의 존재 여부를 돌연변이율에 따라 정량화한다. 낮은 돌연변이율 구간에서는 유전체가 제한된 영역에 머물러 ‘클러스터링’ 현상이 나타나며, 이는 전통적인 종 개념과 일치한다. 반면 높은 돌연변이율에서는 유전체가 넓은 서열 공간에 퍼져 ‘델로컬라이제이션(delocalization)’이 일어나며, 종의 경계가 흐려진다. 흥미롭게도, 이러한 유전체 확산은 전체 적합도 감소와는 무관하게 진행되어, 높은 변이율에서도 집단이 충분히 높은 적합도를 유지할 수 있음을 보여준다.

모델의 파라미터 설정(예: 단백질 길이, 접힘 에너지, 결합 친화도, 돌연변이 효과 분포 등)은 실제 생물학적 데이터와 비교 검증되었으며, 특히 DNA 기반 미생물의 평균 돌연변이율(10⁻⁹~10⁻⁸)과 시뮬레이션에서 관찰된 최종 돌연변이율이 일치한다는 점이 주목할 만하다. 이는 단순한 추상 모델이 아니라, 실제 진화 역학을 재현할 수 있는 ‘미시적’ 접근법임을 시사한다.

결론적으로, 이 연구는 (1) 돌연변이율이 진화 속도와 종 형성에 미치는 비선형적 역할, (2) 에피스타시스가 적합도 지형을 급격히 재구성하는 메커니즘, (3) 돌연변이율 자체가 선택 가능한 형질로 작용해 최적화된 진화 경로를 형성한다는 세 가지 주요 통찰을 제공한다. 이러한 결과는 진화 이론, 합성생물학, 그리고 인공 진화 알고리즘 설계에 모두 중요한 함의를 가진다.