基于标记-释放-重捕数据的蚊虫动态机制-统计推断

《Population Ecology》：Mechanistic-Statistical Inference of Mosquito Dynamics From Mark-Release-Recapture Data

　　

针对蚊媒传染病的生物防治策略需要可靠的野外估计，以评估释放雄蚊的扩散和生存。研究人员提出了一种用于标记-释放-重捕（MRR）数据的机制-统计框架，该框架将（i）具有随机死亡和捕获的基于个体的二维伊藤（It?）扩散模型，与（ii）其反应-扩散近似联系起来，后者产生期望密度和特定诱捕器的捕获量。通过求解反应-扩散系统并将其嵌入泊松（Poisson）观测模型中，对每日诱捕计数进行推断，并通过参数自举（parametric bootstrap）评估不确定性。研究人员首先使用与实验数据高度相似的模拟数据验证参数可识别性。然后，研究人员分析了在埃尔卡诺（El Cano，古巴哈瓦那）进行的一项城市MRR试验，该试验包括四次每周释放的不同标记的雄性埃及伊蚊（Aedes aegypti）以及一个室内BG-Sentinel诱捕器网络。对于当前数据集，齐次迁移率模型和城市-非居住区异质替代模型给出了非常相似的结果。估计结果表明，释放后预期寿命约为5天（总成体年龄8天），5天后的典型位移为180米，诱捕效率与迁移和生存同时被识别。与总结诱捕回报（例如，每日生存的指数拟合或平均重捕距离）的经验分析不同，该机制公式明确耦合了迁移、死亡和捕获，并减少了诱捕器布局带来的偏差——从而实现联合估计和生物可解释参数。该方法提供了一种经过校准、计算高效的流程，用于从稀疏的MRR数据中提取生物可解释参数，并为城市环境中不育昆虫技术（SIT）类干预措施的设计和评估提供了一种有原则的经验总结替代方案。

**论文解读**

**研究背景与问题**
蚊虫是多种人类疾病的主要媒介，包括登革热、寨卡病毒等。由于传统杀虫剂面临耐药性和环境问题，生物防治策略如不育昆虫技术（sterile insect technique, SIT）受到关注。SIT需要释放大量不育雄蚊，但野外释放后雄蚊的扩散能力和生存时间是关键参数，直接影响释放策略的设计与优化。然而，现有方法多依赖经验性总结（如指数拟合或平均重捕距离），未能明确耦合运动、死亡和捕获过程，且容易受诱捕器布局偏差影响。此外，蚊子寿命短、重捕率低导致数据稀疏，传统种群模型难以准确推断。因此，亟需一种能联合估计扩散、生存和诱捕效率的机制性方法。

**研究内容与结论**
研究人员设计了一项在古巴哈瓦那埃尔卡诺（El Cano）进行的标记-释放-重捕（mark-release-recapture, MRR）实验：在固定点连续4周每周释放10,000只经不同颜色荧光粉标记的绝育雄性埃及伊蚊（Aedes aegypti），并在释放点周围布置21个室内BG-Sentinel诱捕器（间距约50–400米），每日检查并记录重捕数据。研究提出一个机制-统计框架（mechanistic-statistical framework）：微观层面采用基于个体的二维伊藤（It?）扩散模型，整合随机死亡和捕获；宏观层面推导出反应-扩散偏微分方程（reaction-diffusion partial differential equation, PDE）作为其确定性极限，计算期望密度和各诱捕器捕获量。通过将PDE解嵌入泊松（Poisson）观测模型，利用最大似然估计（MLE）联合推断扩散系数（μ）、死亡率（λ）和基准诱捕率（α），并采用参数自举（parametric bootstrap）量化不确定性。模拟验证表明参数可识别性良好（齐次模型R2均>0.95，异质模型除一个参数外均>0.94）。应用于实际数据后，齐次迁移模型与异质（城市-非居住区）迁移模型经AIC比较（589.33 vs. 591.34）无显著差异，因此采用更简洁的齐次模型。估计结果显示：释放后预期寿命约5天（加上释放前3天，总成体寿命约8天）；5天后典型位移约180米；每小时在诱捕器位置的捕获概率约为0.14。该方法避免了经验分析中因捕获截断和诱捕器布局产生的偏差，实现了运动、生存与诱捕效率的联合估计。

**技术方法**
研究采用的关键技术包括：① 基于个体的二维伊藤扩散模型（individual-based It? diffusion model），模拟每只蚊子的随机运动、死亡和捕获事件；② 反应-扩散偏微分方程（reaction-diffusion PDE），作为微观模型的确定性极限，计算期望密度与各诱捕器累计捕获数；③ 泊松观测模型（Poisson observation model），将每日诱捕计数视为泊松分布，期望值由PDE解给出；④ 最大似然估计（MLE）结合BFGS算法优化参数；⑤ 参数自举（parametric bootstrap）评估估计标准差，通过IBM生成合成数据集重新估计。样本队列来源于古巴哈瓦那埃尔卡诺的现场MRR实验，包括四次连续释放和21个室内诱捕器网络。

**研究结果**
**3.1 估计程序在模拟数据集上的验证**
通过100个随机参数向量生成合成数据集（模拟实际实验设计），分别测试齐次模型（μ, λ, α）和异质模型（μ_urb, μ_non, λ, α）。估计值与真实值高度一致：齐次情况下μ、λ、α的决定系数（R2）分别为0.99、0.98、0.95；异质情况下分别为0.99、0.87（μ_non）、0.97、0.94。异质模型中非居住区扩散系数（μ_non）在某些参数配置下可识别性稍弱，但整体无系统性偏差。这表明推断方法在整个允许参数空间内稳健。

**3.2 使用实验数据进行参数估计**
将齐次和异质模型拟合到实际MRR数据，通过AIC比较（齐次AIC=589.33，异质AIC=591.34），异质模型因增加一个参数并未带来显著改进，故采用齐次模型。齐次模型最大似然估计（MLE）为：μ=64.00 m2/min，λ=0.212 day^-1，α=293.66 day^-1。据此计算：释放后预期寿命约1/λ≈4.7天（原文表述为约5天），加上释放前3天总寿命约8天；5天后期望位移约180米；每小时在诱捕器位置的捕获概率约为14%。异质模型给出μ_urb=64.24 m2/min，μ_non=76.29 m2/min，但标准差较大（8.74），无法确认需要两个不同扩散系数。模型预测的累计捕获量与观测数据吻合良好，表明该方法能再现早期快速积累和后续饱和趋势。

**讨论总结与结论翻译**
讨论部分强调该机制的-统计框架相比传统经验分析的优势：明确耦合运动、死亡和捕获过程，避免经验性总结的偏差（如捕获截断和诱捕器布局影响）；通过AIC检验表明在当前数据尺度下齐次模型足够解释观测；估计的扩散参数（μ）可用于计算真实种群水平位移，而非受诱捕器布局影响的均值；与Gato et al.（2021）在古巴另一试验中估计的寿命（约3.76天释放后）和总寿命（约7天）基本一致，而扩散位移（180米/5天）高于其报告的平均重捕距离（77.3米），符合机制框架可纠正偏差的预期。计算效率上，PDE基推断比纯IBM蒙特卡洛方法快得多。局限性包括未考虑环境协变量（天气、宿主密度）以及使用固定高斯诱捕核。未来工作可扩展至时变参数和多源释放设计。

**研究结论翻译**：该方法提供了一种经过校准、计算高效的流程，用于从稀疏的MRR数据中提取生物可解释参数（扩散、生存和诱捕效率），并为城市环境中SIT类干预措施的设计和评估提供了一种有原则的经验总结替代方案。