에티오피아 고지대 적응의 유전적 구조

초록

고지대에 사는 에티오피아인(아마라·오로모) 집단을 저지대와 비교해 전장 SNP와 CpG 메틸레이션을 분석하였다. 티베트인에서 발견된 EPAS1·EGLN1 등 혈색소 관련 변이는 에티오피아인에게서는 영향을 주지 않았으며, 아마라에서는 rs10803083이 혈색소와 전장 수준으로 연관됨을 확인했다. 산소포화도와의 연관은 없었다. 대립유전자 빈도 차이 분석에서는 면역·병원체 방어 유전자와 세포주기·DNA 손상·복구 경로가 고지대 적응에 관여할 가능성을 시사한다. 오로모에서는 고·저지대 간 메틸레이션 차이가 몇몇 유전자에서 유의하게 나타났다.

상세 분석

본 연구는 고지대 적응을 이해하기 위해 에티오피아의 두 주요 민족인 아마라와 오로모를 대상으로 고도(고지대 ≈ 3500 m)와 저지대(≈ 1500 m) 거주자를 각각 채집한 뒤, 전장 SNP 유전체와 CpG 메틸레이션 데이터를 통합 분석한 점이 특징이다. 표본 규모는 각각 고·저지대에서 100~150명 수준으로, 인구 구조를 보정하기 위해 PCA와 ADMIXTURE 분석을 수행하였다. GWAS는 혈색소(Hb)와 산소포화도(SpO₂)를 주요 표현형으로 설정하고, 선형 회귀 모델에 연령·성별·고도·주요 인구 구조 주성분을 공변량으로 포함시켰다.

티베트인에서 EPAS1·EGLN1 등 저산소 반응 유전자에 대한 강력한 선택 신호와 Hb와의 연관이 보고된 바 있으나, 에티오피아인에서는 동일 변이가 통계적으로 유의미한 연관을 보이지 않았다. 이는 고지대 적응이 인구마다 다른 유전적 경로를 택할 수 있음을 시사한다. 흥미롭게도 아마라 집단에서는 rs10803083 (염색체 6q21 인근)이 Hb와 전장 수준(p < 5 × 10⁻⁸)으로 연관되었으며, 이 변이는 기존 후보 유전자와는 전혀 겹치지 않는다. 기능적 해석에 따르면 해당 영역은 아직 잘 규명되지 않은 조절 요소일 가능성이 높으며, 향후 기능 실험이 필요하다.

산소포화도와 관련된 GWAS에서는 어느 집단에서도 유의미한 신호가 탐지되지 않았다. 이는 SpO₂가 복합적인 환경·생리적 요인에 크게 좌우되며, 단일 유전자가 큰 효과를 발휘하기 어려운 특성을 반영한다는 해석이 가능하다.

대립유전자 빈도 차이(Fst) 기반 스캔에서는 전통적인 저산소 반응 유전자는 outlier로 나타나지 않았다. 대신 고·저지대 간 차이가 가장 큰 SNP들은 병원체 방어(예: HLA, KIR)와 연관된 유전자에 집중되었으며, 이는 고지대 환경에서 감염 압력이 중요한 선택 요인일 수 있음을 암시한다. 또한, 세포주기 조절·DNA 손상·복구 경로(예: ATM, BRCA1, CHEK2) 유전자군에서 일관된 Fst 상승이 관찰되었다. 고지대의 저산소와 자외선 노출이 DNA 손상률을 높이고, 이를 효율적으로 복구하는 메커니즘이 적응에 기여했을 가능성을 제시한다.

메틸레이션 분석은 오로모 고·저지대 집단을 대상으로 Illumina EPIC 배열을 이용해 수행되었으며, 차등 메틸화된 CpG 사이트는 주로 면역·대사 관련 유전자(예: IL6, SLC2A1)와 연관되었다. 메틸레이션 차이는 전사 조절 수준에서 고지대 적응이 진행될 수 있음을 보여준다.

전체적으로 이 연구는 고지대 적응이 인구별로 서로 다른 유전·표현형 경로를 따를 수 있음을 뒷받침한다. 기존 티베트 연구와는 달리 에티오피아인에서는 저산소 반응 유전자가 아닌 면역·DNA 복구 경로가 선택 압력의 주요 표적이 되었으며, 새로운 GWAS 신호(rs10803083)와 메틸레이션 변이가 고지대 적응 메커니즘을 다각적으로 이해하는 데 중요한 단서를 제공한다. 향후 대규모 메타분석과 기능 검증을 통해 이러한 후보 유전자의 실제 생리학적 역할을 규명할 필요가 있다.