소규모 집단의 성 결정 확률 변동에 대한 멸종 회피 메커니즘

초록

본 논문은 고정된 개체수 N을 가진 작은 성재생 집단이 무작위 성 결정에 의해 전 세대가 남성 또는 여성만으로 구성될 위험을 분석한다. Moran 모델을 이용해 매 시간 단계마다 한 개체가 선택돼 이성 교배 후 자손으로 교체되는 과정을 모사하고, 한 세대( N 단계) 동안 여성 개체 수 k가 무작위 보행을 한다는 점을 강조한다. 초기 성비와 무관하게 빠르게 이항분포에 수렴하지만, 이 분포는 불안정해 결국 k=0 또는 k=N 에 도달해 멸종한다. 평균 멸종 시간은 약 2^N/N 세대로 추정되며, 자손 성 비에 편향이 존재할 경우에도 비슷한 규모의 회복력을 보인다.

상세 분석

이 연구는 작은 개체군이 겪는 ‘성별 고정’ 현상을 확률론적 관점에서 정량화한다. 기본 가정은 개체수가 고정된(N) Moran 과정이며, 매 시간 단계에서 무작위로 한 개체가 선택돼 이성 파트너와 교배한다. 자손은 남·여가 될 확률이 ½(또는 편향된 p)이며, 부모는 사망해 교체된다. 이렇게 정의된 과정은 여성 개체 수 k 를 0과 N 사이의 상태공간을 갖는 일차원 마코프 체인으로 만든다. 전이 확률은
P(k→k+1)= (N−k)/N·k/N·p,
P(k→k−1)= (N−k)/N·k/N·(1−p),
P(k→k)=1−P(k→k+1)−P(k→k−1)
와 같이 표현된다. 여기서 (N−k)/N는 남성 선택 확률, k/N는 여성 선택 확률을 의미한다. 따라서 k 는 평균적으로 ½ N을 중심으로 무작위 보행을 하며, 편향 p≠½일 경우 평균이 p·N으로 이동한다.

분석의 핵심은 두 가지 시간 척도이다. 첫 번째는 ‘혼합 시간’으로, 초기 성비가 어떻든 1세대(N 단계) 안에 k의 분포가 정규화된 이항분포 B(N,½) (또는 B(N,p))에 수렴한다는 점이다. 이는 마코프 체인의 고유값 스펙트럼을 이용해 전이 행렬의 두 번째 고유값 λ₂≈1−2/N이므로, 수렴 속도가 O(N)임을 보여준다. 두 번째는 ‘멸종 시간’이다. 이항분포는 흡수 상태 k=0, N이 아닌 ‘반사’ 상태처럼 보이지만, 실제로는 흡수 경계가 존재한다. 첫 회귀 분석과 전이 행렬의 역행렬을 이용하면 평균 흡수 시간 τ≈2^N/N가 도출된다. 이는 ‘중심극한정리’와 ‘대수적 성장’이 결합된 결과로, 개체수가 20~30 정도면 멸종까지 수천 세대가 걸리지만, N이 50을 넘으면 τ는 천문학적 규모가 된다.

편향된 성 결정(p≠½)도 동일한 스케일을 유지한다. p가 0.6 정도로 크게 치우쳐도 평균 흡수 시간은 여전히 O(2^N/N)이며, 단지 평균 성비가 p·N으로 이동한다. 이는 ‘강인성(robustness)’이라는 용어가 의미하는 바와 일치한다. 실제 생물학적 사례로는 Wolbachia 감염에 의해 수컷을 사망시키는 ‘세균성 성전환’이나, 파충류의 온도 의존성 성 결정이 있다. 이러한 메커니즘이 존재해도, 개체군이 충분히 작다면 확률적 변동에 의해 전 세대가 단일 성으로 전락할 위험은 여전히 존재한다.

수학적 결과는 시뮬레이션과도 일치한다. 저자들은 Monte Carlo 실험을 통해 N=10~30 범위에서 평균 멸종 시간이 2^N/N와 근접함을 확인했으며, 편향 p를 바꾸어도 차이는 미미했다. 또한, 초기 조건(전부 남성, 전부 여성, 혹은 임의의 비율)과 무관하게 1세대 내에 이항분포에 수렴한다는 점을 실증했다. 이러한 결과는 작은 개체군이 ‘성비 균형’을 유지하기 위해 외부 이주나 인공 교배가 필요함을 시사한다.

결론적으로, 이 논문은 작은 집단이 겪는 성별 고정 현상이 통계적으로 거의 불가피함을 수학적으로 증명하고, 그 시간 척도가 2^N/N라는 간단한 식으로 요약될 수 있음을 보여준다. 이는 보전생물학, 특히 멸종 위기 곤충이나 파충류 집단 관리에 직접적인 함의를 제공한다.