고온에 의한 Wolbachia 손실을 고려한 충격 방출 전략

초록

본 논문은 고온 스트레스로 인한 Wolbachia 감염 손실을 포함한 충격 미분 방정식 모델을 제시하고, 방출 규모와 주기, 열 손실 강도 사이의 임계조건을 분석한다. 이론적 임계값과 수치 시뮬레이션을 통해 wMelPop, wMel, wAlbB 균주의 열 내성 및 적합도 비용 차이가 장기 감염 유지에 미치는 영향을 비교한다.

상세 분석

논문은 먼저 야생 모기(S₁)와 Wolbachia 감염 모기(S₂)의 상호작용을 로지스틱 성장 형태와 세포질 불일치(cytoplasmic incompatibility, CI) 파라미터 γ를 포함한 연속 미분식(1)로 기술한다. 여기서 ψᵢ는 출산률, δᵢ는 사망률, K는 공동 환경 수용량을 의미한다. 가정(2)와 (3)에서 ψ₁>δ₁, ψ₂>δ₂, ψ₂<ψ₁, δ₂>δ₁을 두어 감염 모기가 야생 모기에 비해 적합도 비용을 갖는 현실을 반영한다.

정리 1은 초기값이 비음수이면 해가 존재·유일하고, 양의 직교평면이 불변이며, 모든 해가 전역적으로 유계임을 증명한다. 이는 모델이 생물학적으로 의미 있는 해를 제공함을 보장한다.

정리 2에서는 기본 번식 지수 R₁=ψ₁/δ₁, R₂=ψ₂/δ₂를 도입하고, 세 개의 평형점 E₁(야생만 존재), E₂(감염만 존재), E₃(공존)과 그 존재 조건을 제시한다. E₁은 R₁>1일 때 존재하고 R₁>R₂이면 안정적이며, E₂는 R₂>1 및 R₂>(1−γ)R₁일 때 안정적이다. 공존 평형 E₃는 1−γ<R₁R₂<1을 만족할 때만 존재하고, 항상 새들 포인트(saddle)로 불안정함을 보인다. 이는 CI 효과가 충분히 강할 경우 감염 모기가 야생 모기를 완전히 대체할 수 있음을 시사한다.

핵심 기여는 온도에 의한 급격한 감염 손실을 충격 항(term)으로 모델에 삽입한 점이다. 고온 사건은 일정 주기 τₗ에 따라 감염 비율을 α(0<α<1)만큼 감소시키고, 감소된 비율만큼 야생 모기 수가 즉시 증가하도록 설계된다. 동시에 정기적인 방출 충격 τᵣ에 따라 Q개의 감염 모기가 투입된다. 두 충격을 결합한 시스템은 다음과 같은 형태를 가진다:

• 연속 구간: (1)식 그대로 진행
• t = nτᵣ (방출 시점): S₂⁺ = S₂⁻ + Q
• t = mτₗ (고온 손실 시점): S₁⁺ = S₁⁻ + αS₂⁻ , S₂⁺ = (1−α)S₂⁻

이 구조를 이용해 Poincaré 매핑을 구성하고, 주기적 해의 존재와 안정성을 분석한다. 저자는 고온 충격 빈도와 강도(α, τₗ) 및 방출 규모 Q, 방출 주기 τᵣ 사이의 관계를 나타내는 임계식
\