[논문 리뷰] Shrinkage-based random local clocks with scalable inference

Alexander A. Fisher, Xiang Ji외 3명|arXiv (Cornell University)|2021. 05. 15.

Genomics and Phylogenetic Studies인용 수 4

한 줄 요약

이 논문은 로그 분지별 시계 속도의 증분에 무거운 尾를 가진 베이지안 브릿지 사전분포를 적용하여, 시계 위치에 대한 사전 지식 없이도 지역 시계의 효율적 추론을 가능하게 하는 확장 가능한 베이지안 수축-시계 모델을 제안한다. 닫힌 형식의 기울기와 선형 시간 순환 알고리즘을 활용한 하미льт로니안 몬테카를로(Hamiltonian Monte Carlo)를 통해, 무작위 지역 시계(RLC) 대비 3배 이상의 속도 향상을 달성하였으며, 쥐와 인플루엔자 바이러스 계통수에서 알려진 시계를 정확히 복원하였다. 이는 대규모 인플루엔자 표면 당단백질 분석에서도 성공적으로 적용되었다.

ABSTRACT

Local clock models propose that the rate of molecular evolution is constant within phylogenetic sub-trees. Current local clock inference procedures scale poorly to large taxa problems, impose model misspecification, or require a priori knowledge of the existence and location of clocks. To overcome these challenges, we present an autocorrelated, Bayesian model of heritable clock rate evolution that leverages heavy-tailed priors with mean zero to shrink increments of change between branch-specific clocks. We further develop an efficient Hamiltonian Monte Carlo sampler that exploits closed form gradient computations to scale our model to large trees. Inference under our shrinkage-clock exhibits an over 3-fold speed increase compared to the popular random local clock when estimating branch-specific clock rates on a simulated dataset. We further show our shrinkage-clock recovers known local clocks within a rodent and mammalian phylogeny. Finally, in a problem that once appeared computationally impractical, we investigate the heritable clock structure of various surface glycoproteins of influenza A virus in the absence of prior knowledge about clock placement.

연구 동기 및 목표

• 큰 나무에서 무작위 지역 시계 모델의 계산 불가능성을 해결하기 위해 확장 가능한 추론 프레임워크를 개발하는 것.
• 자기상관 시계 모델의 과도한 수축 문제를 해결하기 위해, 중대한 갑작스러운 속도 변화를 허용하는 무거운 꼬리 사전분포를 도입하는 것.
• 분자 시계 분석에서 시계 위치나 수의 사전 지정이 필요 없도록 하는 것.
• 고급 MCMC 샘플링을 통해 대규모 계통수에서 유전적 시계 속도 변화의 효율적이고 정밀한 추론을 가능하게 하는 것.
• 실제 데이터, 특히 인플루엔자 A 바이러스 당단백질에 모델을 적용하여 숙주 특이적 시계 패턴을 밝혀내는 것.

제안 방법

• 부모-자식 노드 간 로그 시계 속도 증분에 평균이 0이고 꼬리가 두꺼운 베이지안 브릿지 사전분포를 적용한 베이지안 수축-시계 모델을 제안한다.
• 베이지안 브릿지의 축약된 스파이크-슬래브 표현을 활용하여, 속도 변화의 흐름을 촉진하면서도 필요한 경우 큰 변화를 허용한다.
• 브릿지 사전분포의 가우시안 척도 혼합 형태의 미분 가능성에 기반해 증분 공간에서 하미틀로니안 몬테카를로(HMC) 샘플링을 수행한다.
• 로그 사후확률의 공동 기울기를 O(N) 시간에 계산하기 위해 후위 순환 알고리즘을 개발하여 나무 크기에 비례하는 선형 확장성을 확보한다.
• HMC 질량 행렬을 사전분포의 헤시안에 기반해 조건화하여 혼합성과 수렴 속도를 향상시킨다.
• BEAST에 모델을 구현하여 시간 측정 가능 나무에 대한 기존 계통수 추론 파이프라인과 통합한다.

실험 결과

연구 질문

• RQ1 무거운 꼬리 사전분포를 가진 베이지안 수축 모델이 시계 위치에 대한 사전 지식 없이도 대규모 계통수에서 지역 시계를 효율적으로 추론할 수 있는가?
• RQ2 제안된 HMC 기반 추론 프레임워크는 큰 나무에 대해 효과적으로 확장되며 고성능 샘플링을 유지하는가?
• RQ3 시뮬레이션 및 실제 데이터 세트에서 제안된 모델이 무작위 지역 시계(RLC) 모델보다 정확도와 속도에서 얼마나 우수한가?
• RQ4 인플루엔자 A 바이러스 표면 당단백질의 유전적 시계 구조는 조류 및 말 숙주에서 어떻게 나타나며, 숙주 특이적 속도 변화는 감지 가능한가?
• RQ5 추론된 시계 위치는 브릿지 지수 α의 선택에 얼마나 민감한가? 커버리지와 계산 비용 사이의 상충 관계는 무엇인가?

주요 결과

• 20개의 시뮬레이션 데이터셋(40개 종)에서 무작위 지역 시계(RLC) 대비 효과적인 표본 크기(ESS)를 1초당 3배 이상 빠르게 확보하였다.
• 모델은 쥐와 척추동물의 적응 방사에서 알려진 4개의 지역 시계를 정확히 복원하였으며, RLC 추정치 및 이전 연구 결과와 일치하였다.
• 인플루엔자 A N7 및 H7 당단백질 계통수에서 각각 1개와 7개의 지역 시계를 확인하였으며, N7 계통수에서는 동서 아메리카 조류 계통 간에 유의미한 속도 변화를 보였다.
• 후행 평균 루트 높이 추정치는 NA N7의 경우 1798(95% HPD: 1733–1855)이며, HA H7의 경우 1853(1808–1897)였으며, 이는 이전 추정치와 5년 이내로, 나무 사전분포가 다름에도 불구하고 일치하였다.
• 모델은 말 선진의 속도 감소를 확인하여 이전 결과를 재확인하였고, NA N7의 동부 조류 클래드에서의 가능성이 있는 시계를 식별하였다.
• 사용자는 브릿지 지수 α를 조정하여 시계 탐지 민감도를 조절할 수 있으며, 더 작은 α 값은 더 높은 시계 해상도를 제공하지만 계산 시간이 증가한다.