Genomic inbreeding coefficient predicts extinction risk

기술 요약: 게놈 근친교배 계수가 멸종 위험을 예측함

문제 제로 (Problem Statement)
동형접합성 연속 구간(runs of homozygosity, $F_{ROH}$) 내의 게놈 비율로 정의되는 게놈 근친교배 계수는 야생 개체군의 근친교배를 평가하는 표준 지표이다. 그러나 $F_{ROH}$ 추정치의 분산 원인은 여전히 명확히 이해되지 않고 있다. 기존 연구들은 $F_{ROH}$가 IUCN 적색 목록(Red List) 상태를 신뢰성 있게 예측할 수 있는지에 대해 상충하는 결과를 내놓았으며, 일부 연구에서는 $F_{ROH}$와 IUCN 적색 목록 상태 사이에 연관성이 없음을 발견하기도 했다. 또한, 관찰된 $F_{ROH}$의 변이가 방법론적 이질성(예: 시퀀싱 깊이, ROH 호출 알고리즘, 길이 임계값)에서 기인하는지, 아니면 종의 정체성, 계통 발생적 역사, 생애사 전략과 같은 생물학적 요인에서 기인하는지도 불분명하다. 이러한 불확실성은 $F_{ROH}$를 보전 평가를 위한 비교 기준선으로 사용하는 데 있어 걸림돌이 된다.

방법론 (Methodology)
저자들은 세 가지 주요 질문을 해결하기 위해 개체군 수준의 $F_{ROH}$ 추정치에 대한 메타 분석적 합성을 수행하였다: (1) 무엇이 $F_{ROH}$의 변이를 가장 크게 설명하는가 (방법론적, 종 수준, 또는 개체군 수준 요인)? (2) 위협받는 종은 위협받지 않는 종보다 높은 $F_{ROH}$를 보이는가? (3) 서로 다른 생애사 전략을 가진 종들은 $F_{ROH}$에서 차이를 보이는가?

• 데이터 수집: 본 연구는 129개의 일차 실증 연구와 4개의 기존 방법론적 표준화 데이터셋을 집계하였다. 최종 데이터셋은 포유류, 조류, 파충류, 어류 및 기타 분류군을 아우르는 454개 종(14,000명 이상의 개체)의 1,088개 야생 개체군을 나타내는 1,491개의 기록로 구성되었다.
• 데이터 추출: 평균 및 표준 편차의 $F_{ROH}$, 표본 크기, 시퀀싱 깊이, ROH 호출기(예: PLINK, BCFtools), 그리고 최소 ROH 길이 임계값에 관한 정보가 추출되었다.
• 공변량: 보전 상태는 IUCN 적색 목록으로부터 얻었다. 생애사 특성(수명, 성숙 연령, 체중, 번식력)은 "느린 생애사 지수(slow life index)"를 구축하기 위해 데이터베이스(예: AnAge, Amniote Life History Database)로부터 수집되었다.
• 통계 모델링: 저자들은 MCMCglmm 패키지를 사용하여 베이지안 계통 발생 일반 선형 혼합 모델(PGLMM)을 채택하였다.
  • 고정 효과 (Fixed Effects): 기술적 변수(ROH 호출기, 길이 임계값, 시퀀싱 깊이) 및 생애사 특성.
  • 변량 효과 (Random Effects): 개체군 정체성, 연구 정체성, 그리고 계통 발생적으로 구조화된 종 항(TimeTree에서 유도된 계통 발생 공분산 행렬을 통해 모델링됨).
  • 결측치 처리: 결측된 표준 편차와 시퀀싱 깊이를 처리하기 위해 다중 대치법(Multiple imputation)이 사용되었다.
  • 민감도 분석: 데이터의 견고성을 테스트하기 위해 특정 분류군, 특정 ROH 길이 임계값(100 kbp 대 1 Mbp), 그리고 합성 데이터셋을 제외한 하위 집합을 사용하여 모델을 재적합하였다.

주요 기여 (Key Contributions)

1. 분산 분할: 본 연구는 $F_{ROH}$ 분산의 포괄적인 분해를 제공하며, 계통 발생적으로 구조화된 종의 정체성이 방법론적 및 개체군 수준의 요인을 압도하는 지배적인 변이 원인임을 입증하였다.
2. 멸종 위험과의 연관성: 방법론적 이질성과 계통 발생적 비독립성을 통제한 후, $F_{ROH}$와 IUCN 멸종 위험 사이의 강력한 양의 상관관계를 확립하였다.
3. 생애사 연결: 생애사 전략과 게놈 근친교배 사이의 관계를 정량화하여, "느린" 생애사(긴 수명, 큰 체구, 늦은 성숙, 낮은 번식력)를 가진 종들이 더 높은 기저 $F_{ROH}$를 갖는 경향이 있음을 보여주었다.
4. 방법론적 통찰: 기술적 선택의 구체적인 영향을 명확히 하여, 높은 ROH 길이 임계값과 특정 호출기(예: PLINK 대 BCFtools)가 $F_{ROH}$ 추정치를 체계적으로 변화시킨다는 점을 확인하였다.

결과 (Results)

• 변이의 원인: 전체 데이터셋에서 계통 발생적으로 구조화된 종 성분이 가장 큰 비율의 분산을 설명하였다(사후 평균 ~61.8%). 방법론적 변수(기술적 공변량)는 변산의 약 7.7%만을 설명하였다. 개체군 및 연구 정체성은 각각 5.8%와 16.7%를 설명하였다.
• 멸종 위험: 멸종 위험과 $F_{ROH}$ 사이에는 일관된 양의 상관관계가 존재한다. IUCN 상태를 순서형 변수(관심 대상에서 위기 단계까지)로 취급했을 때, 위험도가 한 단계 증가할 때마다 $F_{ROH}$는 0.16 로짓 단위만큼 증가하였다. 위협받는 종(취약, 위기, 위급)은 일반적으로 관심 대상 종보다 높은 $F_{ROH}$를 나타냈다.
• 생애사 특성: "느린 생애사 지수"는 $F_{ROH}$와 양의 상관관계를 보였다. 구체적으로 체중은 $F_{ROH}$와 양의 상관관계를, 번식력은 음의 상관관계를 보였다. 이 패턴은 포유류와 조류 전반에서 유지되었으나, 어류의 경우 제한된 표본 크기로 인해 결과가 결정적이지 않았다.
• 계통 발생 신호: $F_{ROH}$의 계통 발생 신호는 0.62(중앙값 0.72)로 추정되었으며, 이는 가까운 친척 관계의 종들이 유사한 기저 수준의 게놈 근친교 de를 공유함을 나타낸다.
• 방법론적 효과: 최소 ROH 길이 임계값이 강력한 영향을 미쳤다. 임계값을 10배 높이면 추정된 $F_{ROH}$가 약 0.82 로짓 단위 감소한다. PLINK는 BCFtools보다 높은 추정치를 생성하였다.

의의 및 주장 (Significance and Claims)
본 논문은 $F_{ROH}$가 멸종 위험의 유효한 비교 지표라고 주장하지만, 그 효용성은 절대적인 값이 보편적인 임계치가 아님을 인식하는 데 달려 있다고 주장한다. 대신, $F_{ROH}$는 진화적 역사와 생애사 전략에 의해 형성된 계통 특이적 기저치와 관련하여 해석되어야 한다.

저자들은 이전 연구들이 $F_{ROH}$와 멸종 위험 사이의 연관성을 감지하지 못한 이유로 두 가지 요인을 제시한다: (1) 최근의 인구학적 감소가 아직 게놈 근친교배로 완전히 나타나지 않은 시간적 불일치("표류 부채(drift-debt)" 가설), 특히 장수하는 종에서 그러하다. (2) 높은 $F_{ROH}$가 최근의 인위적 압력이 아닌 조상의 인구학적 역사나 보존된 생애사 특성을 반영하는 계통 발생적 불일치이다.

따라서 본 연구는 보전 평가가 단일한 높은 $F_{ROH}$ 값을 멸종 위험의 진단 도구로 의존해서는 안 된다고 주장한다. 대신, 평가는 여러 개체군을 샘플링하고, 비교 기준을 설정하기 위해 밀접하게 관련된 종을 활용하며, 방법론적 변이와 계통 발생적 맥락을 고려해야 한다. 이 프레임워크는 기존 문헌의 불일치를 해소하고 생물 다양성 보전에서 게놈 근친교배 추정치의 적용을 개선하는 것을 목표로 한다.