다중인구 SNP 주파수 스펙트럼을 이용한 공동 인구역사 추정

초록

본 논문은 1‑locus 2‑allele Wright‑Fisher 모델의 확산 근사를 이용해 최대 3개 인구의 다차원 SNP 주파수 스펙트럼을 계산하고, 이를 복합 가능도 방식으로 인구역사 모델을 추정한다. 부트스트랩을 통해 연관성을 보정하고, 선택 효과까지 포함한 확장 모델을 제시한다. 인간의 아프리카 외부 확산과 신대륙 정착을 4개 인구(YRI, CHB, CEU, MXL)의 5 Mb 비코딩 데이터에 적용했으며, 비동의성 변이의 빈도 분포도 성공적으로 재현한다.

상세 분석

이 연구는 기존 인구역사 추정 방법이 주로 2인구 간의 1차원 주파수 스펙트럼에 의존하던 점을 넘어, 3인구까지 확장 가능한 다차원 스펙트럼을 수치적으로 계산한다는 점에서 혁신적이다. 핵심은 Wright‑Fisher 이산시간 마코프 체인을 연속시간 확산 방정식으로 근사하고, 경계조건을 통해 인구 규모 변화, 이주, 혼합 등을 모델링한다는 점이다. 저자들은 3차원 편미분 방정식을 유한 차분법으로 풀어 기대 스펙트럼을 얻으며, 이는 복합 가능도(likelihood) 함수에 직접 삽입된다. 복합 가능도는 각 SNP가 독립이라고 가정하지만, 실제 데이터에서는 연관(LD) 때문에 분산이 과소평가될 수 있다. 이를 보완하기 위해 저자들은 블록 부트스트랩을 사용해 실제 연관 구조를 보존한 재표본을 생성하고, 파라미터 추정치의 표준 오차와 가설 검정의 p값을 계산한다.

선택을 포함한 모델 확장은 특히 단일 부위에 대한 선택계수(s)와 새로운 변이의 발생률을 별도로 추정한 후, 해당 선택 강도가 인구 간 주파수 분포에 미치는 영향을 확산 방정식에 추가함으로써 구현된다. 이때 선택은 시간에 따라 변할 수 있으며, 인구 규모 변동과 이주율과도 상호작용한다. 결과적으로, 선택이 포함된 모델은 비동의성(비동의성) 변이의 스펙트럼을 정확히 재현하며, 이는 기존의 중립 모델이 설명하지 못했던 고빈도 비동의성 변이의 과잉을 해소한다.

실제 데이터 적용에서는 68명의 개인(4개 인구)으로부터 5 Mb 비코딩 영역을 재시퀀싱한 데이터를 사용했다. 모델 파라미터는 아프리카에서의 초기 인구 규모, 비아프리카 인구의 급격한 확장, 동아시아와 유럽 인구 간의 이주율, 그리고 멕시코계 라틴아메리카(MXL)의 혼합 비율 등을 포함한다. 부트스트랩 결과는 각 파라미터에 대한 95% 신뢰구간을 제공하고, 가설 검정은 ‘아프리카 외부 확산이 단일 단계가 아니라 연속적인 확장 과정이다’라는 기존 가설을 통계적으로 기각한다.

또한, 저자들은 이전 연구에서 추정된 아미노산 변이의 선택 효과 분포와 현재 추정된 인구역사 모델을 결합해, YRI, CHB, CEU 세 인구의 비동의성 변이 주파수 스펙트럼을 예측한다. 예측값은 실제 관측값과 높은 상관관계를 보이며, 특히 저빈도 변이와 중간 빈도 변이 구간에서 오차가 최소화된다. 이는 인구역사와 선택이 동시에 작용하는 복합 모델이 실제 유전체 데이터 해석에 매우 유용함을 시사한다.

전반적으로 이 논문은 다인구 다차원 스펙트럼을 수치적으로 계산하고, 복합 가능도와 부트스트랩을 결합함으로써 인구역사와 선택을 동시에 추정할 수 있는 강력한 프레임워크를 제공한다. 향후 더 많은 인구와 전장 유전체 데이터를 다룰 때도 확장 가능성이 높으며, 특히 인간 진화사 연구와 선택 스캔의 널 모델 구축에 큰 영향을 미칠 것으로 기대된다.