Continental scale genetic diversity supports geo-clustered ecoregions connectivity

기술 요약: 대륙 규모의 유전적 다양성은 지리적으로 클러스터링된 생태 지역 간의 연결성을 뒷받침한다

문제 제기
전 세계적으로 담수 생물 다양성이 감소하고 있으나, 현재 미국의 보존 전략은 이러한 손실을 완화하기에 불충분하다. 효과적인 생태 모니터링은 종종 미세한 수준의 분류학적 정보 부족으로 인해 방해를 받으며, 이는 대규모 스케일에서의 생태계 상태에 대한 지식 격차로 이어진다. 생태적 연결성은 생태계 회복과 탄력성의 근간이지만, 이를 측정하는 것은 어렵다. 구조적 연결성(공간적 배치)이 대리 지표로 자주 사용되지만, 이는 생물의 이동이나 유전자 흐름(기능적 연결성)을 고려하지 않는다. 이동을 직접 관찰하는 것은 까다로우며, 종 분포 데이터($\beta$-다양성)는 개체군이 동일한 종을 공유함에도 불구하고 고립되어 있을 경우 기능적 연결성을 온전히 나타내지 못할 수 있다. 따라서 대륙 규모에서 유전적 구조를 관찰하여 연결성을 이해하고 취약한 지역을 식별하기 위해 유전 데이터를 활용할 필요가 있다.

연구 방법론
본 연구는 미국 국립과학재단(NSF)의 국립 생태 관측 네트워크(NEON)에서 수집된 저서 수생 절지동물의 대규모 DNA 메타바코딩 데이터를 활용하였다. 데이터셋은 2018년부터 2022년 사이에 호수, 강, 하천에서 수집된 19개 생태 지역(생태 기후 영역), 34개 지점, 470개 샘플을 포함한다.

• 데이터 처리: 원시 시퀀스 파일은 MetaWorks 파이프라인을 사용하여 처리되었다. 이 과정에는 페어드 엔드(paired-end) 리드 병합, 프라이머 트리밍, 디노이징(denoising), 키메라 제거를 통한 정확한 서열 변이(Exact Sequence Variants, ESVs) 생성 단계가 포함된다. 분류학적 할당은 RDP 분류기를 사용하였으며, 가성 의사유전자(putted pseudogenes)는 프로파일 은닉 마르코프 모델(HMM)을 사용하여 필터링하였다.
• 클러스터링: ESVs는 종 수준의 그룹을 나타내기 위해 SWIMA를 사용하여 종 경계 클러스터(Species Bound Clusters, SBCs)로 묶였다.
• 연결성 분석:
  • 구조적 연결성: 경계를 공유하는 생태 지역들이 연결된 것으로 정의되는 인접 행렬에 의해 정의되었다.
  • 기능적 연결성: 생태 지역 간 ESV 및 SBC의 존재/부재로부터 계산된 $\beta$-다양성(Sørensen 비유사도)에 의해 정의되었다.
  • 통계적 비교: 켄달 순위 상관관계(Mantel test)를 사용하여 구조적 연결성과 기능적 연결성을 비교하였다. 비유사도 행렬로부터 그래프가 생성되었으며, 엣지(edge)를 필터링하고 지배적인 클러스터를 식별하기 위해 침투 임계값(percolation thresholds)이 적용되었다.
  • 유전적 분화: 생태 지역 간 종 내 평균 쌍 유사도(mean pairwise Sørensen dissimilarity)를 사용하여 종 내 유전적 분화를 계산하였다. 유전적 클러스터링을 시각화하기 위해 근린 결합 트리(Neighbor-joining trees)가 구축되었다.
• 결과 해석:
  • 공간적 자기상관: 구조적 연결성(인접성)과 기능적 연결성($\beta$-다양성)은 유의미하게 상관관계가 있었다. 인접한 생태 지역은 인접하지 않은 지역보다 더 유사한 생물 다양성을 공유했다.
  • 동서 분단: 전반적인 상관관계에도 불구하고, 로키 산맥과 대략 일치하는 뚜렷한 동서 분단이 나타났다. $\beta$-다양성 그래프에서 생태 지역은 두 개의 약하게 연결된 그룹(동부와 서부)으로 나뉘었다.
  • 서부 지역의 차등적 상관관계: 서부 생태 지역 내에서는 종 수준의 $\beta$-다양성이 인접성과 매우 약한 상관관계를 보인 반면, ESV 수준의 $\beta$-다양성은 강하고 유의미한 상관관계를 보였다. 반면, 동부 생태 지역은 두 수준 모두에서 유의미한 상관관계를 보였다.
• 종 내 분화: 서부 생태 지역에만 존재하는 종들은 동부의 종들에 비해 유의미하게 높은 종 내 유전적 분화를 나타냈다. 클러스터 형성을 위한 최소 임계값($d_{min}$) 또한 서부 종들에게서 더 높게 나타났다.
• 고립: 특정 생태 지역(예: 동남부의 D04 및 D18)은 구조적 고립(예: 해양에 의한 분리 또는 제한된 인접성)과 일치하는 기능적 고립 상태인 것으로 확인되었다.

주요 결과

• 공간적 자기상관: 구조적 연결성(인접성)과 기능적 연결성($\beta$-다양성)은 유의미하게 상관관계가 있었다. 인접한 생태 지역은 인접하지 않은 곳보다 더 유사한 생물 다양성을 공유했다.
• 동서 분단: 전반적인 상관관계에도 불구하고, 로키 산맥과 대략 일치하는 뚜렷한 동서 분단이 나타났다. $\beta$-다양성 그래프에서 생태 지역은 두 개의 약하게 연결된 그룹(동부와 서부)으로 나뉘었다.
• 서부 지역의 차등적 상관관계: 서부 생태 지역 내에서 종 수준의 $\beta$-다양성은 인접성과 매우 약한 상관관계를 보였으나, ESV 수준의 $\beta$-다양성은 강하고 유의미한 상관관계를 보였다. 반면, 동부 생태 지역은 두 수준 모두에서 유의미한 상관관계를 보였다.
• 종 내 분화: 서부 생태 지역에만 존재하는 종들은 동부의 종들에 비해 유의미하게 높은 종 내 유전적 분화를 나타냈다. 클러스터 형성을 위한 최소 임계값($d_{min}$) 또한 서부 종들에게서 더 높게 나타났다.
• 고립: 특정 생태 지역(예: 동남부의 D04 및 D18)은 구조적 고립(예: 해양에 의한 분리 또는 제한된 인접성)과 일치하는 기능적 고립 상태인 것으로 확인되었다.

핵-기여
본 논문은 DNA 메타바코딩 데이터가 단순한 종의 존재 여부를 넘어 유전적 구조와 기능적 연결성을 대륙 규모에서 관찰하는 데 어떻게 사용될 수 있는지를 입증한다. ESV(더 미세한 유전적 해상도를 나타냄)와 종을 모두 분석함으로써, 본 연구는 광범위한 종 분포가 겉으로는 연결되어 보일지라도, 기저의 유전적 구조는 유전자 흐름에 대한 상당한 장벽이 존재함을 보여준다. 또한, 본 연구는 로키 산맥을 유전자 흐로를 제한하는 주요 장벽으로 성공적으로 식별하였으며, 이로 인해 동부와 서부 개체군 사이에 뚜렷한 분리가 발생함을 밝혀냈다.

의의
본 연구는 유전자 흐름이 억제되는 지점을 이해하는 것이 보존 계획 수립에 필수적이라고 주장한다. 연구 결과는 미국 서부의 개체군이 더 고립되어 있고 더 높은 종 내 유전적 분화를 보인다는 것을 시사하며, 이는 로키 산맥의 "스카이 아일랜드(sky island)" 효과와 역사적 기후 변화 때문인 것으로 보인다. 결과적으로, 이들 서부 개체군은 인위적 스트레스와 기후 변화에 더 취약할 수 있으며, 제한된 연결성으로 인해 회복 속도가 느릴 수 있다. 이러한 연결성의 취약 지점을 식별함으로써, 보존 노력은 이처럼 취약하고 고립된 지역을 보호하는 데 집중될 수 있다. 저자들은 메타바코딩 기반의 유전적 다양성 분석을 통합하는 것이 기능적 연결성을 매핑하는 데 도움이 되어, 환경 의사 결정 및 탄력성 계획을 개선할 수 있다고 주장한다.