Persistence of the coral Porites harrisoni in the Persian/Arabian Gulf is underpinned by strong selection on few genomic loci

*기술 요약: 페르시아만/아라비아해 내 Porites harrisoni의 지속성*

문제 제기
산호초를 형성하는 산호들은 상한 열적 한계치에 도달하고 있으며, 가속화되는 해양 온난화에 따라 파괴적인 손실이 발생할 것이라는 전망이 나오고 있다. 페르시아만/아라비아해(PAG)는 여름철 해수 온도가 34°C를 초월하여 대부분의 열대 산호초 열적 한계를 2°C 이상 상회하는 극한 환경을 나타낸다. 이러한 가혹한 조건에도 불구하고, 경성 산호인 Porites harrisoni는 PAG에서 생존하고 있다. 그러나 이러한 열 내성의 유전적 근거, PAG 개체군과 인접한 보다 온화한 환경(예: 오만만, GO)을 연결하는 집단 유전 구조, 그리고 이 맥락에서의 숙주-조류 공생의 특이성은 여전히 잘 알려져 있지 않다. 이러한 메커니즘을 이해하는 것은 보조 유전자 흐름(assisted gene flow)과 같이 다른 지역의 산호 회복력을 높이기 위한 능동적 개입 전략을 평가하는 데 매우 중요하다.

연구 방법론
본 연구는 세 곳의 지점(열적으로 극단적인 남부 PAG의 Al Saada Reef와 Saadiyat Reef, 그리고 환경적으로 온화한 오만만 Al Aqah)에서 P. harrisoni를 조사하기 위해 집단 유전체 접근법을 사용하였다.

• 샘플링 및 시퀀싱: 2022년에 총 119개의 산호 군체를 샘질하였다. Illumina NovaSeq6000을 이용한 전유전체 샷건(WGS) 시퀀싱을 수행하였으며, 평균 커버리지는 45X를 기록하였다.
• 변이 호출 및 필터링: GATK Best Practices를 따라 리드(reads)를 처리하였다. 두 가지 별도의 SNP 필터링 파이프라인을 적용하였는데, 하나는 집단 구조 분석용(저빈도 변이 및 클론 제거)이었고, 다른 하나는 선택적 스윕(selective sweep) 탐지용(고정 신호를 포착하기 위해 희귀 변이 유지)으로 더 허용적인 방식이었다.
• 집단 구조: 분석에는 주성분 분석(PCA), 계통 추정(ADMIXTURE), 쌍별 $F_{ST}$ 계산, 그리고 인구 통계학적 역사 재구성을 위한 Stairway Plot v2가 포함되었다.
• 숙주-공생체 연관성: 호스트 유전형을 Mantel 테스트를 사용하여 조류 공생체 군집(ITS2 유형 프로파일)과 상관 분석함으로써 유전형-유전형 연관성을 평가하였다.
• 선택의 신호: 양의 선택(positive selection) 하에 있는 후보 영역은 Weir & Cockerham의 $F_{ST}$(유전적 차별화)와 교차 집단 확장 하플로타입 동형성(XP-EHH) 통계량을 통합하여 식별하였다. 아웃라이어(outlier) SNP는 SnpEff을 사용하여 주석을 달았으며, 뉴클레오타이드 다양성($\pi$)과 연관 불평형(LD) 붕괴를 평가하여 선택적 스윕의 성격을 규명하였다.

주요 결과

• 집단 구조: 세 지점 간의 전반적인 유전적 차별화는 낮았으며($F_{ST}$ 범위 0.021 ~ 0.041), 이는 호르무즈 해협을 통한 유전자 흐름이 존재함을 나타낸다. 그러나 ADMIXTURE 분석($K=2$) 결과, PAG 개체군은 GO와 비교하여 뚜렷한 조상 계통을 공유하고 있었으며, 가장 남쪽의 PAG 지점(SA)은 혼합(admixture)이 최소화되어 있었다. 인구 통계 재구성 결과, 홀로세(Holocene) 동안의 PAG 형성(~6,000년 전)과 일치하는 유효 개체군 크기($N_e$)의 유의미한 감소가 나타났으며, 이는 창시자 효과(founder effect) 또는 병목 현상을 시사한다.
• 숙주-공생체 특이성: 숙주의 유전적 유사성이 공생체의 정체성을 엄격하게 결정하지는 않았으나, 숙주와 Cladocopium 간의 유전적 거리 사이에서 중간 정도의 유의미한 연관성이 발견되었다. PAG 군체는 주로 열에 강한 Cladocopium thermophilum과 연관되어 있었던 반면, GO 군체는 Durusdinium trenchii 및 Symbiodinium microadriaticum을 보유하고 있었다.
• 선택의 유전적 신호: 본 연구는 GO 산호에서는 관찰되지 않는, PAG 산호에서 강력한 양의 선택을 받는 두 가지 특정 유전 영역을 식별하였다:
  1. Contig 29: 가장 높은 $F_{ST}$ 값(최대 0.93)과 확장된 하플로타입을 보였다. 이 영역은 세포외 기질(ECM) 리모델링 관련 유전자(예: COL4A1, COL2A1), TGF-$\beta$ 신호 전달(ACVR2B, BMPR2), 그리고 세포 주기 조절(CDK2, ORC4) 관련 유전자를 포함한다.
  2. Contig 62: 매우 긴 하플로타입(정규화된 XP-EHH 최대 10.47)과 높은 차별화($F_{ST}$ 최대 0.69)를 보였다. 이 영역은 면역 반응, 스트레스 반응(예: KDELC1), 세포 사멸(PDCD6), 그리고 RNA 가공에 관여하는 유전자가 풍부하다.
• 선택의 성격: 대립유전자 빈도 도표는 GO 산호가 아닌 PAG 산호에서 특정 대립유전자의 거의 완전한 고정(near-fixation)을 보여주었다. 확장된 하플로타입과 함께 뚜렷한 뉴클레오타이드 다양성($\pi$) 감소가 나타나지 않는 것은, 선택이 단일한 새로운 돌연변이가 아닌 기존의 유전적 변이(standing genetic variation)에 작용했음을 시사하는 "연성(soft)" 선택적 스윕임을 의미한다.

의의 및 주장
본 논문은 PAG의 극한 환경에서 P. harrisoni의 지속성이 소수의 유전 좌위에 작용하는 강력한 방향성 선택에 의해 뒷받침된다고 주장한다. "협상 불가능한 SNP(non-negotiable SNPs)"—즉, PAG 개체군에서 거의 고정에 이른 대립유전자—의 식별은 극한의 열 및 염도 스트레스에서 생존하는 데 있어 이러한 특정 유전 변이의 필수적인 역할을 강조한다.

본 연구는 선택적 육종을 통한 보조 유전자 흐름과 같은 향후 능동적 개입 전략을 위한 유전적 표적을 식별하는 첫 단계라고 주장한다. 그러나 저자들은 중요한 제약 조건을 강조한다: PAG 산호의 적응은 (높은 염도를 포함한) 걸프 지역의 특정 환경 조건과 밀접하게 연결되어 있다. 따라서 이러한 유전형을 더 차갑고 염도가 낮은 환경(예: GO)으로 이전하는 것은 장기적으로 실행 불가능할 수 있다. 본 연구의 결과는 지역적 유전 다양성이 귀중한 자원이기는 하나, 극한 환경의 좁은 선택적 체제(selective regime)가 미래의 적응 잠재력을 제한할 수 있음을 강조하며, 국지적 개입만으로는 부족한 전 지구적 기후 완화 전략이 더욱 시급함을 시사한다.