Penna 모델을 이용한 혼합생식 집단의 게놈 진화

초록

본 연구는 두 개의 개구리 종(Rana lessonae, Rana ridibunda)에서 발생하는 잡종(Rana esculenta)의 혼합생식 현상을 Penna 모델을 통해 시뮬레이션하였다. 잡종은 한 부모의 유전체를 배란 전 제거하고 남은 유전체를 무성으로 전달함으로써 클론과 성적 번식이 결합된 형태를 보인다. 시뮬레이션 결과, 장기간 고립된 잡종 집단은 상보적인 대립유전자를 축적해 새로운 종으로 분화하지만, 외부 잡종 유입이 소량이라도 존재하면 정화 선택이 우세해 종분화가 억제된다.

상세 분석

Penna 모델은 개체의 수명과 발현되는 유전자를 연령에 따라 순차적으로 활성화시키는 비선형 이산 모델로, 각 유전자는 0(정상) 또는 1(결함)으로 표시된다. 본 연구에서는 두 부모 종의 유전체를 각각 32비트(또는 64비트) 길이의 이진 문자열로 설정하고, 잡종은 한 쪽 부모의 전체 염색체를 배란 전 제거한 뒤 남은 염색체를 그대로 전달한다. 이렇게 하면 잡종의 생식세포는 ‘클론’처럼 동일한 haplotype을 유지하면서도, 체세포에서는 두 부모의 대립유전자가 서로 보완되어 이형접합 상태를 유지한다.

시뮬레이션 파라미터는 다음과 같다: 초기 인구 10,000명, 최대 수명 32세, 매 세대마다 2개의 자손을 생산, 돌연변이율 μ=0.001, 그리고 ‘돌연변이 누적 임계치’를 3개 결함이 누적되면 사망하도록 설정하였다. 고립된 잡종 집단에서는 무성 전달된 haplotype이 세대 간 거의 변하지 않기 때문에, 결함이 누적된 대립유전자는 이형접합으로 보완된다. 결과적으로 전체 결함 비율은 증가하지만, 표현형상으로는 정상 개체가 유지된다. 이는 ‘보완적 대립유전자 집합(complementary haplotype set)’이 형성되어 종 내부에 고정되는 메커니즘을 보여준다.

반면, 외부에서 매 세대 1% 수준의 새로운 잡종이 투입되면, 이들 개체는 서로 다른 haplotype을 가지고 있어 재조합이 발생한다. 재조합은 기존 보완적 집합을 파괴하고, 결함이 동일한 염색체에 겹치게 하여 이형접합 보완 효과를 감소시킨다. 따라서 자연 선택은 결함을 가진 개체를 제거하는 ‘정화 선택(purifying selection)’으로 전환된다. 이 과정에서 잡종 집단은 기존의 보완적 구조를 유지하지 못하고, 결국 부모 종과 유사한 유전적 구성을 회복한다.

핵심 인사이트는 다음과 같다. 첫째, 혼합생식은 유전적 ‘고정점(fixed point)’을 만들 수 있는 강력한 메커니즘이며, 이는 장기간 고립 시 새로운 종으로의 분화 가능성을 제공한다. 둘째, 소량의 유전자 흐름(gene flow)이 존재하면 이러한 고정점을 깨뜨리고, 전통적인 다윈식 정화 선택이 우세하게 된다. 셋째, Penna 모델은 연령 의존적 발현과 결함 누적을 동시에 고려함으로써, 실제 개구리 집단에서 관찰되는 ‘연령에 따른 사망률 증가’를 정량적으로 재현한다. 마지막으로, 모델은 잡종의 ‘클론‑성적 혼합’ 전략이 환경 변화와 유전자 흐름에 얼마나 민감한지를 보여주어, 복합 종(complex species) 이론에 실험적 근거를 제공한다.