분자 발산 시간 추정 모델 비교 연구

초록

분자 진화 속도 변화 모델은 속도 자체, 로그속도, 역속도 세 가지를 기준으로 하고, 시간에 따른 변동은 무시, 선형, 속도·시간 곱 세 가지 형태로 나뉜다. 총 아홉 가지 모델을 포유류 계통수에 적용해 상관관계와 적합도를 비교했으며, 스케일링 팩터를 포함하면 보정점 위치에 관계없이 상대적 시간 추정이 일관됨을 확인했다.

상세 분석

본 논문은 분자 발산 시간 추정에 사용되는 확률적 모델을 체계적으로 분류하고, 실제 데이터에 적용해 그 차이를 정량화한다. 먼저 속도 변화에 대한 페널티 함수를 정의하는 기준을 세 가지(속도, 로그속도, 역속도)로 설정하고, 시간 의존성을 ‘시간에 무관’, ‘시간에 선형 비례’, ‘속도와 시간의 곱’ 세 가지 형태로 확장한다. 이 조합은 3×3=9개의 모델을 만든다. 각 모델은 랜덤 워크와 브라운 운동 형태를 포함하며, 파라미터 추정 시 변동성(variance) 스케일링 팩터를 도입할지 여부도 변수가 된다.

연구자는 포유류 계통수(루트 고정, 내부 노드 하나 고정)를 대상으로 각 모델을 적용해 추정된 분기 시점을 비교하였다. 모델 간 유사성은 추정 시점의 피어슨 상관계수로 측정하고, 이를 최소제곱 트리 형태로 시각화했다. 또한 모델 적합도는 로그우도와 Q‑Q 플롯을 통해 검증했으며, 특히 라우사이아테리아(Laurasiatheria)와 수프라프리마테스(Supraprimates) 내부 클레이드의 연령 차이가 모델에 따라 크게 변동함을 발견했다.

특히 화석 교정점으로 자주 사용되는 말/코뿔소와 고래/하마 노드가 수프라프리마테스 내부 연대와 충돌하는 경우가 빈번했으며, 60백만 년 이전의 설치류 치아 화석이 발견되지 않은 점과도 불일치한다. 이는 화석 교정점 선택이 전체 시간 추정에 미치는 영향을 강조한다.

스케일링 팩터를 포함한 모델은 교정점 위치가 바뀌어도 상대적 연대 비율이 거의 일정하게 유지되는 반면, 스케일링 없이 단순히 페널티만 최소화하는 모델은 비현실적으로 긴 시간 추정을 초래한다. 따라서 변동성 스케일링은 모델의 안정성을 크게 향상시킨다.

결론적으로, 모델 선택에 따라 특정 클레이드의 연대가 크게 달라질 수 있음을 보여주며, 특히 라우사이아테리아 내의 오래된 분기 시점은 모델 의존성이 높다. 연구자는 다중 모델 비교와 스케일링 팩터 도입을 통해 보다 견고한 시간 추정 프레임워크를 제안한다.