계통도 유전체 상위상 복원: 오류·결측·재조합을 SAT로 해결

초록

본 논문은 재조합, 유전체 오류 및 결측치를 동시에 고려한 계통도 기반 하플로타입 추론 문제(HCRE)를 정의하고, 이를 SAT(부울 만족도) 문제로 효율적으로 변환한 뒤 최신 SAT 솔버를 이용해 정확하고 실용적인 해법을 제시한다. 실험에서는 시뮬레이션 및 실제 젖소 집단 데이터를 통해 높은 정확도와 합리적인 실행 시간을 입증한다.

상세 분석

이 연구는 기존 MRHC(Minimum‑Recombinant Haplotype Configuration) 모델이 실제 데이터에서 흔히 나타나는 두 가지 핵심 결함—유전체 오류와 결측값—을 전혀 반영하지 못한다는 점을 지적한다. 이를 보완하기 위해 저자들은 HCRE(Haplotype Configuration with Recombinations and Errors)라는 새로운 문제 정의를 제시한다. HCRE는 (r, e)‑HC 형태로, 허용 가능한 재조합 수 r과 허용 가능한 오류 수 e를 명시적으로 제한함으로써 연구자가 실제 실험 조건(재조합율, 오류율 등)에 맞는 파라미터를 직접 설정할 수 있게 한다.

문제의 계산 복잡도는 APX‑hard임을 증명함으로써 다항시간 정확 알고리즘이 존재하지 않을 가능성을 명확히 한다. 그럼에도 불구하고, 저자들은 SAT(부울 만족도) 문제로의 정밀한 변환을 통해 실용적인 해법을 구현한다. 변환 과정은 크게 세 가지 제약군으로 구성된다. 첫 번째는 멘델리안 유전법칙을 만족하도록 부모‑자식 간의 소스 벡터와 하플로타입 비트를 연결하는 논리식이다. 두 번째는 관측된 유전체와 하플로타입 사이의 일관성을 검증하며, 오류가 발생한 경우 e_i