[논문 리뷰] Interbreeding conditions for explaining Neandertal DNA in living humans: the nonneutral case

A. Neves|arXiv (Cornell University)|2011. 05. 30.

Pleistocene-Era Hominins and Archaeology참고 문헌 17인용 수 3

한 줄 요약

이 논문은 이전의 인간-네안데르탈 교잡에 대한 중립 모델을 보완하여 적합도 차이를 도입함으로써, 해부학적으로 현대 인간의 1%의 적합도 우위조차도 중동 지역에서의 장기적 공존을 설명할 수 없을 정도로 네안데르탈의 멸종을 너무 빨리 일으킨다는 점을 보여준다. 비중립 모델은 비록 관찰된 비아프리카인들 내의 1–4% 네안데르탈 DNA를 달성하기 위해 더 높은 교잡 비율이 필요하다고 예측하지만, 이와 같은 빠른 멸종은 이스라엘 동굴에서 관찰된 네안데르탈과 현대 인간 간의 장기적 교대와 모순된다. 따라서 중립 모델은 유전적 자료와 고고학적 자료를 모두 설명하는 데 더 타당하다는 결론을 내린다.

ABSTRACT

We consider here an extension of a previous work by Neves and Serva, still unpublished, which estimates the amount of interbreeding between anatomically modern Africans and Neandertals necessary for explaining the experimental fact that 1 to 4% of the DNA in non-African living humans is of Neandertal origin. In that work we considered that Africans and Neandertals had the same fitness (neutral hypothesis) and Neandertal extinction was thus an event of fortune. In this work we consider that Africans had larger fitnesses. We show results for four values for the fitness difference: 1%, 5%, 10% and 20% and compare them with the corresponding neutral results. Some technical differences with respect to the neutral case appear. We conclude that even with 1% fitness difference Neandertals extinction comes up in too small a time, so the neutral model looks more suitable for explaining the known data on occupation of some caves in Israel for a very long time, alternately by Africans and Neandertals.

연구 동기 및 목표

• 비중립 모델(현대 인간이 네안데르탈보다 높은 적합도를 가짐)이 비아프리카 인구 내 1–4%의 네안데르탈 DNA를 설명할 수 있는지 조사하기.
• 현대 인간과 네안데르탈 간의 적합도 차이가 중동 지역의 장기적 공존 고고학적 기록과 호환되는지 평가하기.
• 이전에 연구된 중립 모델과의 비교를 통해 비중립 모델의 예측(특히 필요한 교잡 비율과 멸종 시점)을 분석하기.
• 다양한 인구 규모와 적합도 차이에서 비중립 모델이 사용한 결정론적 근사의 타당성을 평가하기.

제안 방법

• 연구는 상수 인구 크기 $ N $ 를 가진 두 개의 하위집단(현대 인간과 네안데르탈)을 대상으로 하는 결정론적 평균장 모델을 사용한다. 세대 간에 겹치지 않는 모델이다.
• 현대 인간이 네안데르탈보다 상대적 적합도 우위를 가지는 적합도 차이 $ h $ 를 도입하고, 시간에 따른 인간 집단 내 네안데르탈 유전자형 비율을 모델링한다.
• 결과의 $ N $-독립성을 확보하기 위해, 교잡 비율을 스케일링한 파라미터 $ \beta = \alpha / N $ 를 사용한다. 여기서 $ \alpha $ 는 세대당 교잡 비율이다.
• 비중립 가정 하에 네안데르탈 DNA 기여의 확률 분포를 유도하고, 수치 시뮬레이션과 해석적 근사법을 사용하여 중립 사례($ h = 0 $)와 비교한다.
• 멸종 시기는 초기 인간 비율과 적합도 차이 $ h $ 의 함수로 결정론적으로 계산되며, 세대 길이를 20년으로 가정한다.
• 모델의 타당성은 조건 $ \alpha \ll 1 / \log_2 N $ 를 점검하여 평가되며, 이는 평균장 근사가 유지됨을 보장한다.

실험 결과

연구 질문

• RQ1 현대 인간에게 1%의 적합도 우위를 도입할 경우, 비아프리카인들 내 1–4%의 네안데르탈 DNA를 달성하기 위해 필요한 교잡 비율은 어떻게 변하는가?
• RQ2 비중립 모델은 중동 지역 동굴에서 관찰된 네안데르탈과 현대 인간 간의 장기적 교대(예: 스키울, 카프자 등)를 설명할 수 있는가?
• RQ3 적합도 차이에 따라 네안데르탈 멸종 시점은 어떻게 변하며, 이는 스키울과 카프자와 같은 유적지에서 관찰된 13만 년 간의 인구 교대와 일치하는가?
• RQ4 실제 인구 규모와 교잡 비율을 고려할 때, 비중립 모델에서 사용된 결정론적 근사는 어느 정도 유효한가?
• RQ5 유전적 자료와 공존 고고학적 증거를 동시에 설명하는 데 있어 중립 모델이 비중립 모델보다 더 적합한가?

주요 결과

• 현대 인간에게 1%의 적합도 우위($ h = 0.01 $)가 있더라도, 네안데르탈 멸종은 약 600~1,000년 내에 발생하며, 이는 이스라엘 동굴에서 관찰된 13만 년 간의 인구 교대와는 너무 빠르게 일치하지 않는다.
• $ h = 0.01 $ 일 때 필요한 교잡 파라미터 $ \beta $ 는 약 $ 9.3 \times 10^{-4} $ 이며, 확률 밀도의 최대값은 $ \beta_{\mathrm{max}} \approx 5.0 \times 10^{-4} $ 에서 발생하여, 매우 높은 교잡 비율이 필요함을 시사한다.
• 적합도 차이 $ h $ 가 증가할수록 필요한 $ \beta $ 는 증가하고 멸종 시간은 감소한다. $ h = 0.2 $ 일 경우 멸종은 600~1,000년 내에 발생하며, $ h = 0.01 $ 일 경우 11,000~20,000년 내에 발생한다.
• 모델은 $ N = 1000 $ 일 때 $ \beta \gg 10^{-4} $ 인 경우 결정론적 근사가 붕괴됨을 보여주며, 이는 더 큰 $ \beta $ 에서의 결과를 주의 깊이 해석해야 함을 시사한다. 다만 더 작은 $ N $ 은 정확도 향상에 기여할 수 있다.
• 비록 비중립 모델이 수학적으로 단순하지만, 중립 모델은 여전히 장기적 공존 고고학적 기록과 더 잘 일치하므로, 유전적 자료와 고고학적 자료를 동시에 설명하는 데 더 적합하다.