정적 환경에서 수평 유전자 전달이 단세포 집단 평균 적합도에 미치는 영향

초록

**
본 논문은 접합에 의한 수평 유전자 전달(HGT)이 무성생식 단세포 원핵생물 집단의 돌연변이‑선택 평형에 미치는 영향을 수학적으로 모델링한다. 집단은 고정된 농도의 항생제 존재 하에서 저항성을 획득해야만 생존할 수 있다고 가정한다. 항생제에 의해 유도되는 1차 사멸 속도 상수가 낮은 경우와 높은 경우 두 극한에서 모델의 거동을 분석했으며, 평균 적합도는 낮은 접합률에서 가장 높게 나타난다. 접합률이 임계값을 초과하면 평균 적합도는 최소값으로 감소한 뒤, 접합률이 무한히 커질 때까지 점진적으로 상승하지만, 그 한계값은 낮은 접합률에서 얻은 평균 적합도보다 낮다. 고려한 파라미터 범위에서는 평균 적합도의 접합률 의존성이 비교적 작으며, 항생제 사멸 상수가 0에 접근할수록 이 의존성은 사라진다. 항생제 사멸 상수가 크게 될 경우, 평균 적합도의 거동에 대한 해석적 해를 도출했으며, 이는 시뮬레이션 결과와 잘 일치한다. 연구 결과는 접합 매개 HGT가 돌연변이‑선택 평형 상태에서 평균 적합도에 약간의 해로운 영향을 미친다는 것을 시사한다. 따라서 HGT는 새로운 혹은 변하는 환경에 대한 빠른 적응을 가능하게 함으로써 선택적 이점을 제공한다는 결론을 내린다. 이러한 결과는 환경 스트레스가 HGT를 촉진한다는 관측과 일치한다.

**

상세 요약

**
이 연구는 수평 유전자 전달(HGT)이 미생물 집단의 진화 역학에 미치는 영향을 정량적으로 평가하기 위해, 접합(conjugation) 과정을 명시적으로 포함한 확률론적 모델을 구축하였다. 모델은 세 가지 주요 요소—(1) 무성생식에 의한 복제, (2) 돌연변이에 의한 유전적 부하, (3) 항생제에 의한 1차 사멸—를 결합하고, 접합을 통해 플라스미드나 저항성 유전자를 다른 세포로 전달하는 과정을 첫 번째 차수 반응으로 기술한다. 이러한 설계는 실제 미생물 군집이 겪는 ‘정적’ 환경, 즉 항생제 농도가 시간에 따라 변하지 않는 상황을 가정함으로써, 복잡한 환경 변동성을 배제하고 HGT 자체의 효과를 명확히 드러내는 데 목적이 있다.

모델 해석은 두 극단, 즉 항생제 사멸 상수 ( \kappa ) 가 거의 0에 가까운 경우와 매우 큰 경우로 나뉜다. ( \kappa \to 0 ) 일 때는 항생제 압력이 미미하므로, 평균 적합도는 주로 돌연변이 부하와 접합에 의해 발생하는 유전적 교환에 의해 결정된다. 이 경우 접합률 ( \gamma ) 가 낮을수록 평균 적합도가 최댓값에 가깝게 유지되며, 접합이 거의 없을 때와 동일한 수준을 보인다. 이는 HGT가 별다른 선택적 이점을 제공하지 못하고, 오히려 불필요한 유전적 교환으로 인해 적합도가 소폭 감소할 가능성을 시사한다.

반면 ( \kappa ) 가 크게 증가하면 항생제에 대한 저항성 유전자를 보유한 세포만이 살아남을 수 있다. 여기서 접합은 저항성 유전자를 빠르게 퍼뜨리는 메커니즘으로 작용한다. 저접합률 영역에서는 저항성 유전자가 충분히 전파되지 않아 평균 적합도가 낮게 유지되지만, 접합률이 임계값을 초과하면 저항성 유전자의 급격한 확산으로 인해 평균 적합도가 일시적으로 최소값에 도달한다. 이는 과도한 접합이 비저항성 유전자를 동시에 섞어 넣어 전체 적합도를 낮추는 ‘유전적 부하’를 초래하기 때문이다. 접합률이 무한대로 향하면 시스템은 ‘완전 혼합’ 상태에 가까워지며, 평균 적합도는 새로운 평형값에 수렴한다. 그러나 이 평형값은 저접합률에서 얻은 최적 적합도보다 낮으며, 이는 HGT가 장기적으로는 약간의 비용을 동반한다는 결론을 뒷받침한다.

수치 시뮬레이션과 해석적 해의 일치성은 모델이 실제 미생물 집단의 동태를 충분히 포착하고 있음을 보여준다. 특히, 높은 ( \kappa ) 영역에서 도출된 폐쇄형 해는 접합률과 사멸 상수 사이의 비선형 상호작용을 명확히 설명한다. 그러나 모델은 몇 가지 제한점을 가진다. 첫째, 접합을 1차 반응으로 단순화함으로써 플라스미드 복제 비용, 수용체 세포의 상태, 접합 효율성의 변동성을 무시한다. 둘째, 항생제 농도가 일정하다는 가정은 실제 임상 혹은 환경 상황에서 흔히 관찰되는 농도 변동을 반영하지 못한다. 셋째, 돌연변이율을 일정하게 유지했지만, 스트레스 상황에서는 돌연변이율이 상승할 수 있다는 점을 간과한다.

향후 연구에서는 (1) 플라스미드 유지 비용과 전이 효율을 파라미터화한 확장 모델, (2) 동적 항생제 농도 프로파일을 포함한 비정상 환경 시나리오, (3) 스트레스‑유도 돌연변이와 HGT 간의 상호작용을 고려한 다중 스케일 모델링을 제안한다. 이러한 접근은 HGT가 ‘빠른 적응’이라는 선택적 이점을 제공하면서도, 장기적인 평균 적합도에는 미미하거나 약간 부정적인 영향을 미칠 수 있다는 현재의 결론을 보다 정교하게 검증할 수 있을 것이다.

**