Integrating genomics and evolution into the conservation and breeding management of a Threatened Australian frog

기술 요약: 위기 종인 호주 개구리의 보전 및 번식 관리에 유전체학 및 진화론 통합

문제 제술
보전 번식 프로그램(CBP)은 멸종 위기 종을 보호하는 데 매우 중요하지만, 가계 기록이 불완전하거나 없는 경우 유전적 다양성을 효과적으로 관리하는 데 필요한 유전체 데이터가 부족한 경우가 많습니다. 호주 뉴사우스웨일스(NSW) 동부 온대림에 서식하는 고유종이자 대형 지상성 양서류인 남부 스터터링 개구리(Mixophyes australis)는 서식지 상실과 항암균(Chytridiomycosis, Batrachochytrium dendrobatidis에 의해 발생)으로부터 심각한 위협을 받고 있습니다. 이 종은 상당한 범위 축소를 경험했으며, 시드니 남쪽의 개체군들은 대부분 절멸했습니다. 2022년에 보험 개체군을 조성하고 재도입을 촉진하기 위해 새로운 CBP가 설립되었습니다. 그러나 인구 구조, 유전자 흐름 패턴, 그리고 창시자(founders)들의 유전적 연관성에 대한 명확한 이해 없이는 근친 교배 저하, 적응 잠재력 상실, 또는 서로 다른 진화 계통이 혼합됨에 따른 잡종 저하(outbreeding depression)의 위험이 있습니다. 본 연구는 새로 기술된 이 위기 종의 관리를 위한 의사 결정을 알리기 위해 진화 생물학과 유전체학을 통합하는 과제를 다룹니다.

연구 방법론
본 연구는 두 가지 별개의 데이터 세트를 분석하기 위해 다양성 어레이 기술(DArTseq)을 통해 생성된 전 유전체 단일 염기 다형성(SNP) 데이터를 활용하였습니다:

1. 보전 데이터 세트: Symbio 야생 동물 공원(두 곳의 야생 계곡 시스템, 개체군 A와 B에서 채집)의 CBP 창시자 73개 샘플과 재도입된 개체군(R1) 5개 샘플.
2. 야생 데이터 세트: Mahony 등(2023)으로부터 확보한 59개 샘플로, 현존 범위(포트 맥쿼리에서 시드니 사이)의 Mixophyes australis 샘플 55개(개체군 S1–S4)와 멸종된 범위(시드니에서 노우라 사이)의 샘플 4개를 포함합니다.

데이터 처리는 dartRverse 패키지를 사용하여 RStudio에서 엄격한 필터링 과정을 거쳤으며, 이를 통해 재현성이 낮은 로커스(loci)를 제거하고, 읽기 깊이(read depth)를 조정하며, 미세 대립유전자 빈도(minor allele frequency)를 필터링했습니다. 서로 다른 시기에 생성된 데이터 세트 간의 검증 편향(ascertainment bias)을 해결하기 위해, 분석은 종종 독립적으로 수행되었으며, 교차 비교가 필요한 경우 지표들이 표준화되었습니다.

주요 분석 접근 방식은 다음과 같습니다:

• 개체군 구조: 주좌표 분석(PCoA), 쌍별 고정 지수($F_{ST}$) 계산, Mantel 테스트를 통한 격리 거리(IBD) 테스트, 그리고 베이지안(STRUCTURE) 및 희소 비음 행렬 분해(sNMF) 방법을 이용한 클러스터링 분석.
• 유전적 다양성: 관찰된($H_O$) 및 기대($H_E$) 헤테로접합성, 근친 교배 계수($F_{IS}$), 대립유전자 풍부도(AR), 사적 대립유전자(PA), 그리고 유효 개체군 크기($N_e$) 계산.
• 연관성 추정: 창시자들 사이의 친족 계수($\Phi$)를 결정하기 위해 세 가지 서로 다른 추정치인 유전적 관계 행렬(GRM), EMIBD9(우도 기반 기대값-최대화 방법), 그리고 COANCESTRY에 구현된 TrioML을 비교하였습니다. 여기에는 가까운 친족을 식별하기 위한 연관성 네트워크 시각화가 포함되었습니다.

주요 결과

• 개체군 구조: 야생 개체군은 상당한 유전적 구조화를 보입니다. 헌터 강(Hunter River) 지역을 경계로 북부 개체군(S1–S3)과 남부 개체군(S4) 사이에 명확한 분리가 존재합니다. 남부 개체군(S4)은 유전적으로 독특하며 북부 그룹들과 매우 높은 차별성($F_{ST} > 0.25$)을 보입니다. 북부 범위 내에서는 격리 거리가 뚜렷하게 나타나지만, S1과 S2는 지리적 분리에도 불구하고 더 밀접하게 군집을 이룹니다. 반면, CBP 창시자(개체군 A와 B) 및 재도입된 개체군(R1)은 유전적으로 유사한 소스 풀에서 유래했음을 나타내는 최소한의 유전적 차별성($F_{ST} < 0.1$)을 보입니다.
• 유전적 다양성: 남부 개체군(S4)은 낮은 헤테로접합성, 높은 근친 교배($F_{IS} = 0.599$), 낮은 대립유전자 풍부도, 그리고 매우 낮은 유효 개체군 크기($N_e \approx 1.7$)를 특징으로 하는 심각하게 침식된 유전적 다양성을 나타냅니다. 북부 개체군들은 중간 정도의 헤테로접합성을 보이지만 근친 교배 수준은 다양합니다. S3는 인구학적으로 가장 건강해 보이며, S1과 S2는 더 높은 고유 다양성(사적 대립유전자)을 보유하고 있습니다. 보전 개체군(A+B)은 중간 정도의 헤테로접합성과 낮은 근친 교배를 보이며, 이는 건강한 야생 개체군과 유사하여 창시자들이 유전적 변이를 대표성 있게 포착했음을 시사합니다.
• 연관성 및 방법론: 세 가지 연관성 추정치(GRM, EMIBD9, TrioML) 모두 쌍별 연관성의 상대적 순위를 일관되게 도출하였으나, 절대값은 달랐습니다. GRM과 TrioML은 번식 개체군에 대해 낮은 배경 친족 관계($\Phi \approx 0.007–0.015$)를 추정하였는데, 이는 3촌 관계 또는 무관한 개인 수준입니다. EMIBD9은 매우 근친 교배된 작은 개체군을 위해 설계된 특성상 더 높은 추정치($\Phi \approx 0.09–0.10$)를 냈습니다. 26개의 활성 번식 개체군 하위 집단 내에서, GRM과 TrioML은 약 8개의 부모-자식/형제 관계, 1개의 이복 형제, 그리고 2개의 사촌 관계를 식별했습니다. 반면 EMIBD9은 더 적은 수의 가까운 관계를 식별했습니다.

의의 및 주장
본 논문은 가계 기록이 없는 상황에서 유전체 데이터가 어떻게 직접적으로 보전 번식 결정을 안내할 수 있는지를 입증함으로써, M. australis의 장기적인 유전적 관리를 위한 기초적인 프레임워크를 제공한다고 주장합니다. 본 연구는 다음과 같은 기여를 강조합니다:

1. 관리 단위 정의: 북부와 남부 개체군 사이의 유의미한 유전적 단절을 식별하였으며, 이는 남부 혈통(S4)에 별도의 관리가 필요함을 시사합니다. 따라서 잔존 남부 개체군에 대한 추가 샘플링이 이루어질 때까지 북부와 남부 창시자를 혼합하는 것에 주의가 필요합니다.
2. 창시자 평가: 현재의 CBP 창시자들이 유전적으로 다양하고 유의미한 구조화를 보이지 않음을 확인하였으며, 이는 즉각적인 잡종 저하의 위험 없이 프로그램 내에서 유연한 배우자 선택 전략을 사용할 수 있음을 의미합니다.
3. 방법론적 지침: 연관성 추정치를 비교함으로써, 저자들은 번식 쌍을 선택할 때 절대적인 값보다는 연관성의 상대적 순위에 집중할 것을 권고하며 관리자들에게 실질적인 프레임워크를 제공합니다. 저자들은 속도를 위해 GRM을, 견고함을 위해 TrioML을 추천하며, 더 크고 근친 교배가 적은 데이터 세트에 대한 EMIBD9 사용은 경계할 것을 조언합니다.
4. 보전 전략: 본 연구 결과는 국지적 적응을 보존하면서 근친 교배를 최소화하기 위해 통제된 교환을 동반한 여러 개의 반격리(semi-isolated) 하위 개체군을 유지하는 전략을 뒷받침합니다. 연구는 지속적인 유전적 모니터링, 잔존 남부 개체군에 대한 추가 샘플링, 그리고 재도입 전 개체군 혼합의 생존력을 검증하기 위한 실험적 교배 시험의 필요성을 강조합니다.

궁극적으로, 본 논문은 진화 과정과 유전체학을 보전 평가에 통합하는 것이 항암균과 같은 신종 질병에 직면한 위기 종의 적응 잠재력을 유지하는 데 필수적이라고 주장합니다.