[论文解读] A Simple Model of Epidemics with Pathogen Mutation

Michelle Girvan, Duncan S. Callaway+2 more|RePEc: Research Papers in Economics|May 18, 2001

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) 预测。
• 对于大种群规模，渐近近似预测的灭绝边界与数值模拟结果高度吻合，尤其在高接触率下。
• 在高度连通的种群（高 z）中，病原体必须快速突变才能避免灭绝，表明存在强烈的免疫逃逸进化压力。
• 模型表明，极短和极长的免疫力持续时间均会抑制振荡：前者因无法积累易感者池，后者则因最低感染水平过低导致病原体灭绝。