[논문 리뷰] A Simple Model of Epidemics with Pathogen Mutation

Michelle Girvan, Duncan S. Callaway외 2명|RePEc: Research Papers in Economics|2001. 05. 18.

Mathematical and Theoretical Epidemiology and Ecology Models인용 수 6

한 줄 요약

이 논문은 병원체 변이와 숙주 면역 기억을 통합한 전파 역학의 계산 모델을 제안하며, 변이율과 면역 지속 기간 간의 상호작용으로 인해 진동하는 유행과 동적 멸종이 발생함을 보여준다. 주요 발견은 고도로 연결된 인구에서 병원체가 지속되기 위해 빠르게 변이해야 한다는 점으로, 밀도 높은 숙주 네트워크에서 면역 회피를 위한 강력한 진화적 압력이 존재함을 시사한다.

ABSTRACT

We study how the interplay between the memory immune response and pathogen mutation affects epidemic dynamics in two related models. The first explicitly models pathogen mutation and individual memory immune responses, with contacted individuals becoming infected only if they are exposed to strains that are significantly different from other strains in their memory repertoire. The second model is a reduction of the first to a system of difference equations. In this case, individuals spend a fixed amount of time in a generalized immune class. In both models, we observe four fundamentally different types of behavior, depending on parameters: (1) pathogen extinction due to lack of contact between individuals, (2) endemic infection (3) periodic epidemic outbreaks, and (4) one or more outbreaks followed by extinction of the epidemic due to extremely low minima in the oscillations. We analyze both models to determine the location of each transition. Our main result is that pathogens in highly connected populations must mutate rapidly in order to remain viable.

연구 동기 및 목표

• 병원체 변이와 숙주 면역 기억이 유행 역학에 미치는 상호작용을 이해하기 위해.
• 유행 시스템에서 만성 감염, 주기적 유행, 동적 멸종으로의 전이를 모델링하기 위해.
• 진동이 발생하고 병원체가 낮은 감염 최소 수준으로 멸종하는 조건을 규명하기 위해.
• 접촉률과 면역 지속 기간이 유행의 지속 또는 종료에 미치는 영향을 분석하기 위해.

제안 방법

• 비트스트링 모델은 고정 길이 ℓ의 이진 수열로 병원체 유전자형을 표현하며, 면역 반응은 이전에 만난 유전자형과의 하밍 거리 기반으로 한다.
• 개개인은 메모리 레퍼토리에 있는 모든 유전자형과 하밍 거리가 임계값을 초과하는 경우에만 감염 가능하다.
• 면역을 단순화하기 위해 차분 방정식 모델을 사용하여 면역 지속 기간 τ를 고정하고, τ 시간 단위 동안 면역 상태를 유지한 후 다시 감염 가능 상태로 복귀한다.
• 모델은 인구 수준의 역학을 시뮬레이션하기 위해 맵 기반 접근법을 사용하며, 시간에 따라 감염 가능자, 감염자, 면역자 비율을 추적한다.
• 홉프 분기 분석을 통해 면역 지속 기간 τ와 접촉률 z에 따라 진동 행동가 시작되는 조건를 규명한다.
• 최소 감염 비율이 1/N 이하로 떨어질 때의 동적 멸종 조건을 유도하기 위해 점근적 근사법을 사용한다.

실험 결과

연구 질문

• RQ1 병원체가 변이하고 면역 기억이 존재하는 인구에서 어떤 조건에서 진동하는 유행이 발생하는가?
• RQ2 면역 지속 기간이 접촉률과 어떻게 상호작용하여 유행이 지속되거나 멸종하는가?
• RQ3 진동 역학에서 감염의 최소 수준은 무엇에 의해 결정되며, 이는 언제 동적 멸종으로 이어지는가?
• RQ4 인구 연결성(접촉률)은 병원체 변이에 대한 진화적 압력에 어떻게 영향을 미치는가?
• RQ5 지속적인 진동과 멸종을 구분하는 면역 지속 기간과 접촉률 간의 임계 관계는 무엇인가?

주요 결과

• 면역 지속 기간 τ와 접촉률 z가 중간 범위에 있을 경우, 진동 역학이 호프 분기 현상으로 발생한다.
• 진동의 시작은 면역 지속 기간과 접촉률의 곱이 임계값을 초과할 때 발생하며, 이는 감염 가능자 재충전과 병원체 지속성 간의 균형에 의해 결정된다.
• 진동에서 최소 감염 비율이 1/N 이하로 떨어질 때 동적 멸종이 발생하며, 이 임계값은 점근적 공식 p_min = z^(-1/2(τ²−τ+2)) (τ log z)^(τ−1) 로 예측된다.
• 대규모 인구에서 점근적 근사로 예측한 멸종 경계는 수치 시뮬레이션과 매우 유사하며, 특히 높은 접촉률에서 잘 일치한다.
• 높은 접촉률을 가진 고도로 연결된 인구에서는 병원체가 멸종을 피하기 위해 빠르게 변이해야 하며, 이는 면역 회피를 위한 강력한 선택 압력을 시사한다.
• 모델은 매우 짧거나 매우 긴 면역 지속 기간이 모두 진동을 억제함을 보여주며, 짧은 기간은 감염 가능자 풀의 축적을 방해하고, 긴 기간은 낮은 최소 감염 수준으로 인해 병원체 멸종을 초래한다.