유전체 데이터가 밝힌 로마인 역사의 퍼즐

초록

본 연구는 6개 로마 집단을 대상으로 수십만 개 SNP를 분석하여, 로마인의 유전적 조상이 약 80% 서유럽·남아시아 혼합이며, 남아시아와 유럽 조상의 혼합 시기는 약 850년 전이라고 추정한다. 동유럽이 주요 유럽 기원이며, 인도 북서부가 남아시아 기원의 핵심 지역임을 확인했다. 또한, 인도 출발 직후 강력한 창시자 효과가 있었고, 유럽 도착 후 급격한 인구 확장이 일어났음을 제시한다.

상세 분석

이 논문은 로마인(Roma)의 기원과 이동 경로를 고해상도 전장 유전체(genome‑wide) 데이터로 재구성한 최초의 시도 중 하나이다. 연구팀은 인도, 동유럽, 서유럽, 중동 등 다양한 인구집단과 비교하기 위해 6개의 로마인 집단(동유럽·발칸반도·서유럽 등)에서 600,000여 개 SNP를 정밀히 genotyping하였다. 주요 분석 흐름은 다음과 같다.

1. 조상 구성비 추정: ADMIXTURE와 f4‑statistics를 활용해 로마인의 유전적 구성요소를 분해하였다. 결과는 서유럽·동유럽·중동 등 서유라시아 계열과 인도 남아시아 계열이 각각 약 60%와 20%를 차지하고, 나머지 20%는 추가적인 소수 기원을 포함한다는 점을 보여준다. 특히, 인도 내에서도 북서부(펀자브·라자스탄 인근) 지역이 가장 높은 기여도를 보였으며, 이는 언어학적·역사적 증거와 일치한다.

2. 혼합 시기 추정: ALDER와 ROLLOFF 같은 LD‑decay 기반 방법을 적용해 남아시아와 유럽 조상의 혼합 시점을 추정하였다. 평균 850년 전(±≈100년)이라는 결과는 로마인이 중세 초기에 유럽에 정착하기 시작한 시점과 일치한다. 지역별로는 동유럽(폴란드·루마니아)에서 약 900년, 발칸반도에서 약 800년 정도로 미세한 차이가 관찰되었다.

3. 창시자 효과와 인구 규모 변동: ROH(run of homozygos이션)와 IBD(identical‑by‑descent) 분석을 통해 로마인 집단 내 높은 동형접합 구간이 존재함을 확인했다. 특히, 평균 ROH 길이가 2–4 Mb에 달하며, 이는 강력한 창시자 효과를 의미한다. 이어서 MSMC와 SMC++를 이용한 인구 규모 추정에서는 인도 출발 직후 급격한 인구 감소(Ne 약 1,000~2,000) 후, 유럽 도착 후 10세기 내에 Ne가 30,000 이상으로 급증한 패턴을 보였다. 이는 ‘병목 현상 + 급격한 인구 팽창’ 모델을 강력히 지지한다.

4. 지리적 기원 탐색: qpAdm와 qpGraph를 활용해 다양한 후보 모델을 비교했으며, 최적 모델은 ‘인도 북서부(펀자브) + 동유럽(폴란드·우크라이나) + 소규모 남부 유럽(이탈리아)’ 조합으로 설명되었다. 특히, 동유럽이 유럽 조원의 주된 공급원이라는 점은 이전 연구에서 제시된 ‘동부·중부 유럽 경유’ 가설을 유전학적으로 뒷받침한다.

5. 문화·역사와의 연계: 유전적 결과를 로마인의 언어학적 분류와 중세 문헌 기록과 비교했을 때, 약 9세기~12세기 사이에 ‘인도 → 중동 → 동유럽 → 서유럽’ 순으로 이동했으며, 이 과정에서 현지 인구와의 혼합이 지속적으로 일어났음을 시사한다. 또한, 창시자 효과는 초기 소수 집단이 장거리 이동을 감당하기 위해 강력한 혈연 결속을 유지했을 가능성을 제시한다.

전반적으로 이 연구는 로마인의 복합적인 유전적 역사를 고해상도 SNP 데이터와 최신 인구유전학 모델을 결합해 정량화함으로써, 기존의 Y‑염색체·미토콘드리아 기반 연구가 놓쳤던 ‘전체 유전체 수준의 혼합·병목·팽창’ 과정을 명확히 밝혀냈다.