Geographic distance rather than climate shapes genetic differentiation in Acacia senegal in Niger

*기술 요약: 니제르의 Acacia senegal에서 유전적 차별화를 형성하는 것은 기후가 아닌 지리적 거리임*

문제 제기
자연 집단 내 유전적 변이의 공간적 분포는 유전자 흐름, 유전적 부동(genetic drift), 인구통계학적 역사, 그리고 자연 선택의 상호작용에 의해 결정된다. 두 가지 주요 메커니즘이 유전적 차별화를 생성한다: 지리적 격리에 따라 유전적 차이가 축적되는 격리 거리 모델(Isolation by Distance, IBD), 그리고 생태적 이질성이 유전자 흐름을 제한하거나 분화된 선택을 유도하는 환경 격리 모델(Isolation by Environment, IBE)이다. 지형적 위치와 기후 조건이 강하게 상관되어 있는 경우—많은 자연계에서 흔히 나타나는 현상—이 두 효과를 분리해내는 것은 경관 유전학(landscape genetics)의 중요한 과제이다. 니제르의 Acacia senegal 집단은 서쪽에서 동쪽으로 이어지는 가파른 건조도 및 온도 구배를 따라 분포하며, 이는 지리적 위치와 기후 조건이 본질적으로 정렬되어 있는 시나리오를 형성한다. 이전의 게놈 연구들은 이 집단들에서 뚜렷한 유전적 클러스터를 식별하고 온도와 관련된 특정 로커스(loci)를 탐지했으나, 게놈 전반의 집단 차별화가 주로 지리적 격리에 의해 발생하는지, 아니면 거리와 독립적인 기후 적응에 의해 발생하는지는 여전히 불분명했다.

연구 방법론
본 연구는 니제르의 세 가지 아라비아 검(gum arabic) 생산 지역(서부, 중부, 동부)에 분포하는 8개의 A. senegal 자연 집단(총 81개 개체)을 분석하였다. 연구진은 이전의 집단 게놈 연구(Assoumane et al., 2026)에서 사용된 쌍별 $F_{ST}$ 추정치와 지리적 좌표를 활용하였으며, 이를 CHELSA V2.1 및 CHELSA-BIOCLIM+ (1981–2010)로부터 얻은 독립적인 기후 데이터와 통합하였다.

주요 방법론적 단계는 다음과 같다:

1. 변수 선택: 열 극한, 열 범위, 수분 공급, 수분 균형, 증발 요구량을 나타내는 5가지 생물 기후 변수(BIO5, BIO7, BIO12, CMI, VPD)를 유지하였다. 중복되는 변수들은 제외되었다.
2. 거리 행렬:
   • 유전적 거리: 선형화된 $F_{ST}$ ($F_{ST}/(1-F_{ST})$).
   • 지리적 거리: 대원 거리(Great-circle distance, Haversine 공식).
   • 기후적 거리: 잔류 변수의 주성분 분석(PCA)으로부터 도출된 유클리드 거리(변동의 80% 이상을 포착함), 그리고 민감도 테스트를 위해 표준화된 변수로부터 직접 도출함.
3. 통계적 프레임워크:
   • 무작위화를 이용한 다중 행렬 회귀(MMR): 반응 행렬의 공동 순열을 통해 유의성을 테스트하며, 지리적 거리와 기후적 거리가 유전적 차별화와 동시에 연관되는 정도를 추정하기 위한 주요 추론 도구로 사용되었다.
   • 맨텔(Mantel) 및 부분 맨텔(Partial Mantel) 검정: 이변량 연관성을 조사하고 혼란 변수를 통제하기 위해 사용되었으나, 공간적으로 구조화된 데이터에서의 한계점을 인정하였다.
   • 거리 기반 분산 분할(db-RDA): 유전적 변이를 고유한 지리적, 고유한 기후적, 공유된, 그리고 설명되지 않는 부분으로 나누기 위한 탐색적 차원 축소 접근법이다.
4. 강건성 검증: 지리적 변환 및 기후 거리 정의를 변경하는 세 가지 민감도 분석과, 개별 집단 제거에 따른 결과의 안정성을 평가하기 위한 "하나의 집단 제외(leave-one-population-out)" 분석을 수행하였다.

주요 결과

1. 지리적 거리의 지배성: 쌍별 유전적 차별화는 지리적 거리와 매우 강한 양의 상관관계를 보였다 (Mantel $r = 0.986$, $p = 0.0001$). 주요 MMR 모델에서 지리적 계수는 크고 매우 유의미했으며 ($\beta_{geo} = 0.980$, $p = 0.0001$), 유전적 거리 행렬 변동의 97.3%를 설명하였다.
2. 기후 연관성의 취약성: 단순 맨텔 검정에서는 유전적 차별화가 기후적 거리와 양의 상관관계($r = 0.604$)를 보였으나, 지리적 거리를 통제했을 때 이 관계는 사라졌다. 동시 MMR 모델에서 기후 계수는 0에 가깝고 유의하지 않았다 ($\beta_{clim} = 0.011$, $p = 0.772$).
3. 강건성: 지리학이 지배적인 동인이라는 결론은 모든 민감도 분석(로그 변환 거리 또는 대안적 기후 정의 사용)과 8회의 "하나의 집단 제외" 반복 분석에서 일관되게 유지되었다. 지리적 계수는 모든 반복에서 유의미하게 나타난 반-면, 기후 계수는 유의미하지 않았다.
4. 분산 분할: 탐색적 db-RDA 결과, 조정된 유전적 변동의 56.0%가 지리학에 고유하게 귀속되었고, 12.1%가 기후에 고유하게 귀속되었다. 고유 기후 분율은 통계적으로 지지되었으나 ($p = 0.0135$), 지리적 분율보다 실질적으로 작았다. 저자들은 이 방법이 MMR과 다른 통계적 가정을 기반으로 하므로 탐색적 보완책으로 간주되어야 한다고 언급하였다.
5. 공간적 vs 기후적 디커플링: 지리적으로는 동부에 위치하지만 기후적으로는 중부 집단과 더 가까운 하나의 집단(Malan Mainari)이 제거되었을 때 모델 계수의 가장 큰 편차를 보였으나, 지리학적 지배성이라는 전체적인 패턴은 지속되었다.

의의 및 주장
본 논문은 니제르의 A. senegal 집단에서 유전적 차별화의 구조를 형성하는 주된 요인이 지리적 거리라고 결론짓는다. 관찰된 패턴은 강한 격리 거리(IBD) 신호와 일치한다.

저자들은 연구의 의의에 대해 다음과 같이 몇 가지 겸허하면서도 구체적인 주장을 펼친다:

• 이전의 로커스 특이적 발견과의 차별성: 본 연구는 SNP와 최저 기온(BIO6) 사이의 로커스 특이적 연관성에 대한 이전의 발견을 부정하는 것이 아니다. 대신, 특정 로커스에서 국지적 적응이 일어날 수는 있지만, 게놈 전반의 집단 차별화는 광범위한 기후 적응보다는 주로 지리적 거리에 의해 구동된다는 점을 명확히 한다.
• 방법론적 필요성: 지리적 구배와 기후 구배가 정렬된 시스템에서는 지리학과 유전 사이의 단순 상관관계만으로는 IBE의 증거로 불충분하다는 것을 본 연구는 입증한다. 두 효과를 분리하기 위해서는 동시 모델링(MMR)과 분산 분할법의 사용이 필수적이다.
• 보전적 함의: 지리적 신호의 안정성은 보전 샘플링 및 생식질 수집 시 니제르 범위 내의 지리적 피복을 우선시하여 유전적 차별화의 주요 축을 포착해야 함을 시사한다. 그러나 저자들은 지리적 차별화 그 자체를 국지적 기후 적응의 증거로 해석해서는 안 된다고 경고하며, 종자 이동 구역(seed-transfer zones)을 정의하기 위해서는 더 밀도 높은 샘플링과 프로비넌스 시험(provenance trials)이 필요하다고 덧붙였다.
• 한계점: 저자들은 적은 수의 집단(n=8)으로 인한 통계적 검정력의 한계와, 미세 기후 또는 역사적 정착 데이터가 아닌 장기 격자형 기후 데이터를 사용했다는 점을 명시적으로 인정하였다. 이들은 거의 완벽한 모델 적합도가 표본 크기가 크기 때문이 아니라 공간적 구배의 강도를 반영하는 것이라고 설명하였다.

요약하자면, 본 논문은 니제르의 A. senegal에 있어 유전적 변이의 공간적 구조가 주로 지리적 거리에 의해 결정되며, 기후의 기여는 부차적이고 더 세밀한 데이터 없이는 지리적 구배와 구별하기 어렵다는 강력한 증거를 제공한다.