초고속 시퀀싱과 타깃 캡처로 밝힌 조류 진화계통도

초록

본 연구는 조류의 주요 분기인 Neoaves의 계통 관계를 해소하기 위해 초보존 유전자(ultraconserved elements, UCE) 1,541개를 타깃 캡처와 고속 시퀀싱으로 수집하였다. 32종의 Neoaves와 닭을 포함한 33종을 분석했으며, 결측 데이터를 허용한 87% 완전성 매트릭스와 100% 완전성 매트릭스(416 loci)를 각각 Bayesian 및 최대우도(ML) 방법으로 추정하였다. 결과는 1,541 loci를 이용한 경우 높은 해상도와 일관된 트리를 제공했으며, 기존 연구와 일치하는 물새·육새 군집의 단일계통성, 조류와 앵무새의 친연 관계, 맹금류의 비단일계통성 등을 재확인했다. 특히, 형태적으로 크게 차이나는 황새와 해오라미, 뿔새와 투라코스 간의 예상치 못한 근연 관계가 새롭게 제시되었다. 데이터 양을 늘릴수록 트리의 안정성이 향상된다는 점을 강조하며, 아직 해결되지 않은 깊은 분기들에 대한 추가 연구 필요성을 제시한다.

상세 분석

이 논문은 Neoaves라는 조류 대분류 내에서 급격하고 오래된 방사형 분화가 초래한 계통학적 난제를 최신 계통유전체학 기법으로 접근한다. 핵심 방법론은 초보존 유전자(ultraconserved elements, UCE) 타깃 캡처와 Illumina 기반 고속 시퀀싱이다. UCE는 진핵생물 전반에 걸쳐 보존된 비코딩 영역으로, 변이율이 낮아 깊은 진화적 시점에서도 정렬이 용이하고, 주변 가변 부위가 충분히 존재해 계통 신호를 제공한다는 장점이 있다. 연구팀은 기존에 설계된 조류 전용 UCE 프로브 세트를 활용해 32종 Neoaves와 닭(Outgroup)에서 총 1,541개의 UCE를 목표로 하였으며, 각 종당 평균 30배 이상의 커버리지를 확보했다.

데이터 전처리 단계에서는 Trimmomatic으로 저품질 리드와 어댑터를 제거하고, SPAdes를 이용해 각 종의 어셈블리를 수행했다. 이후 PHYLUCE 파이프라인을 통해 UCE 리드 매핑, 컨센서스 서열 생성, 다중 정렬(MAFFT) 및 트리 구축을 진행했다. 두 종류의 데이터 매트릭스를 구성했는데, 하나는 결측 데이터를 허용해 87% 완전성을 유지한 1,541 loci(총 1.2 Mbp)이며, 다른 하나는 100% 완전성을 강제해 416 loci(약 300 kbp)만을 남겼다.

계통 추정은 Bayesian (MrBayes)과 최대우도 (RAxML) 두 방법을 병행했으며, 각각 10⁷ 세대의 MCMC와 1000 부트스트랩 반복을 수행했다. 흥미로운 점은 1,541 loci 매트릭스에서는 Bayesian과 ML 트리가 전반적으로 일치했으며, 주요 분기(예: 물새군, 육새군, 조류-앵무새 군)의 지지값이 95% 이상으로 매우 강했다. 반면 416 loci 매트릭스에서는 두 방법 간 위상 차이가 몇몇 깊은 분기에서 나타났고, 부트스트랩 지지도도 낮아 데이터 양이 트리 안정성에 미치는 영향을 명확히 보여준다.

결과적으로, 기존 연구에서 제시된 몇몇 논란의 관계가 재확인되었다. 특히, Passeriformes(참새목)와 Psittaciformes(앵무목)의 친연 관계는 높은 확률로 지지되었으며, Falconidae와 Accipitridae가 각각 독립적인 군집을 형성한다는 비단일계통성도 확인되었다. 새롭게 제시된 관계로는 Phaethontidae(황새)와 Eurypygidae(해오라미)가 근연이며, Otididae(뿔새)와 Musophagidae(투라코스)가 서로 가까운 친연 관계를 가진다. 이러한 결과는 형태적·생태적 차이에도 불구하고 유전적 신호가 깊은 공통 조상을 반영한다는 점을 시사한다.

연구는 또한 데이터 결측이 트리 구축에 미치는 영향을 평가했는데, 87% 완전성 매트릭스가 충분히 높은 지지를 제공함을 보여준다. 이는 대규모 유전체 데이터에서 완전성을 강제하는 것이 오히려 정보 손실을 초래할 수 있음을 의미한다. 마지막으로, 저자들은 현재 남아 있는 불확실한 깊은 분기(특히 물새·육새 외부의 몇몇 작은 군집)를 해결하기 위해서는 더 많은 UCE loci와 전장 게놈 시퀀싱이 필요하다고 강조한다.