埃及伊蚊卵滞育对沃尔巴克氏体感染蚊释放防控虫媒病毒效能的影响机制研究

《Theoretical Ecology》：The impact of egg quiescence on the efficacy of Wolbachia-infected mosquito releases for arbovirus control

　　

本研究提出一个常微分方程模型，用以解析卵滞育对释放沃尔巴克氏体（Wolbachia）感染蚊防控虫媒病毒传播效能的作用机制。该模型最多可呈现五个平衡点及四种具有生物学意义的场景：两种群均灭绝；未感染种群持续存在而感染种群灭绝；感染种群持续存在而未感染种群灭绝；以及两种群共存。共存场景的出现源于系统存在双稳态（bistability）特性。敏感性分析结果表明，滞育对蚊虫适合度的影响取决于不同环境条件下存活与发育进程之间的权衡。由于携带沃尔巴克氏体的蚊卵无法在滞育中存活，或滞育后羽化的成蚊不育，滞育对感染蚊的适合度产生负面影响，从而降低蚊群中感染的流行率。这增加了遭遇未感染蚊的概率，而相较于沃尔巴克氏体感染蚊，未感染蚊更可能携带登革病毒，因此在触发滞育的环境胁迫期间或之后，登革传播风险可能上升。更为重要的是，滞育可通过增加释放期间所需感染蚊的数量，削弱感染在蚊群中的建立。因此，在滞育高发地区采用沃尔巴克氏体感染蚊防控虫媒病毒传播的策略可能面临失效风险。

论文解读：卵滞育对沃尔巴克氏体感染蚊释放防控虫媒病毒效能的影响

一、研究背景与问题提出

登革热、寨卡、基孔肯雅热及黄热病等虫媒病毒病在全球热带与亚热带地区持续流行，其主要传播媒介为埃及伊蚊（Aedes aegypti）和白纹伊蚊（Aedes albopictus）。在缺乏有效且可负担疫苗的情况下，控制媒介蚊虫种群仍是减缓病原传播的核心策略。近年来，释放沃尔巴克氏体（Wolbachia）感染蚊已成为一种新兴的替代性防控手段，通过种群替换或种群抑制策略降低或消除媒介种群，从而阻断虫媒病毒向人群传播。沃尔巴克氏体诱导的细胞质不相容（cytoplasmic incompatibility, CI）和母系遗传（maternal inheritance）特性，以及其对宿主病毒复制的抑制作用，是实现该策略的生物学基础。

然而，环境胁迫因素对沃尔巴克氏体感染蚊释放效能的潜在影响尚未得到充分研究。滞育（diapause）与静止（quiescence）是蚊虫应对环境胁迫的重要休眠机制。其中，静止是一种非季节性的、直接的、快速且可逆的发育停滞状态，常见于热带蚊虫的卵期。静止卵在不利环境条件下可保持活力长达120天，并通过经卵传播（transovarial transmission）维持病原体传播潜力。已有实验研究表明，沃尔巴克氏体感染蚊卵在长期储存（即静止状态）后活力显著下降，且由此羽化的雌成蚊常表现不育。尽管温度对细胞质不相容和母系遗传效率的直接影响已被数学模型所探讨，但触发休眠机制的环境胁迫条件对释放策略效能的影响此前尚未被建模分析。

二、研究内容与核心结论

针对上述空白，研究人员构建了一个包含八个仓室的常微分方程模型，将未感染与沃尔巴克氏体感染的蚊群分别划分为卵、静止卵、幼虫加蛹（未成熟期）及成蚊四个阶段，系统解析了卵静止对两种蚊群适合度、感染流行率及虫媒病毒传播风险的影响。该模型最多可呈现五个平衡点，对应四种生物学场景：两种群灭绝、仅未感染种群持续、仅感染种群持续以及两种群共存。共存场景的存在源于系统在特定参数条件下呈现双稳态特性，即系统的长期动态对初始条件敏感。

核心结论包括：第一，静止对未感染蚊群适合度的影响取决于环境条件——在不利环境（如极端温度或低湿度）下，静止卵的存活-孵化潜力高于活跃卵，增加静止率可提升未感染蚊群适合度；而在有利环境下则相反。第二，由于沃尔巴克氏体感染蚊卵在静止后存活率低或羽化成蚊不育，静止对感染蚊群适合度产生负面影响，降低感染在蚊群中的流行率。第三，感染流行率下降增加了人类遭遇未感染蚊（更可能携带登革病毒）的概率，从而在触发静止的环境胁迫期间或之后可能轻微增加登革传播风险。第四，更为关键的是，静止可削弱感染在蚊群中的建立，显著增加释放期间所需感染蚊的数量阈值，在静止高发地区可能危及沃尔巴克氏体蚊释放策略的有效性。

三、关键技术方法

研究人员采用的理论框架为非线性常微分方程（ODE）动力学模型，将蚊群按感染状态与生活史阶段划分为八个仓室，并引入密度制约的种内竞争机制。模型解析方面，通过求解平衡点并运用Routh-Hurwitz判据等稳定性分析技术，确立了各平衡点的存在条件与局部渐近稳定性条件。数值模拟采用Python的SciPy库进行长时间（10?年）仿真以确保系统达到平衡。全局敏感性分析采用拉丁超立方采样（Latin Hypercube Sampling, LHS）结合偏秩相关系数（Partial Rank Correlation Coefficient, PRCC）方法，量化各参数对关键阈值参数（R_u、R_w、R_uw）及感染流行率、传播风险输出的影响。模型参数主要来源于已发表的蚊虫生物学与沃尔巴克氏体感染研究文献。

四、研究结果

平衡点的存在性与稳定性分析：通过解析求解与数值模拟相结合，研究人员证明了系统最多存在五个平衡点。当R_u＞max{1, R_w}时，感染-free平衡点P_(u,0)局部渐近稳定；当R_w＞max{1, νR_u}且母系遗传效率ζ=1时，感染种群独占平衡点P_(0,w)稳定。这一结果表明，仅当母系传递完全时，感染蚊才能实现对未感染蚊群的完全替换。共存平衡点的存在性与数量由二次多项式A?ō2+B?ō+C?=0的根决定，在特定参数条件下可出现两个共存平衡点（P_(u,w)^?稳定、P_(u,w)⁺为鞍点），呈现双稳态特性。

静止对蚊群适合度的影响：通过归一化敏感性指数分析，研究人员发现静止率η_o对未感染蚊群适合度R_u的影响方向取决于σ_q/μ_q与σ_o/μ_o的比值关系：在不利环境下（σ_q/μ_q＞σ_o/μ_o），增加静止率提升适合度；在有利环境下则相反。对于感染蚊群，由于μ_q^w显著高于μ_q（反映静止卵存活后羽化成蚊不育导致的繁殖潜力损失），σ_q^w/μ_q^w＜σ_o^w/μ_o^w的条件更易满足，因此增加静止率η_o^w降低感染蚊群适合度R_w，可能危及其长期持续存在。

双稳态与释放阈值：相位图分析表明，在共存场景下系统呈现双稳态，长期动态对初始条件高度敏感。稳定性图显示，释放时感染蚊与未感染蚊成蚊数量比（A^w/A)(0)存在一个临界阈值，低于该阈值感染无法建立，高于该阈值则可实现共存或感染蚊占优。随着静止率η_o增加，达到共存所需的感染蚊释放数量显著增加，说明静止削弱了感染建立的效率。

虫媒病毒传播风险评估：在有利与不利两种环境场景下模拟发现，不利环境导致总体蚊群数量下降，未感染蚊比例在环境胁迫期间及之后有所上升。当静止率η_o与η_o^w从0.022增至0.2时，遭遇未感染蚊的比值比（odds ratio, OR_A）随之增加，表明在触发静止的环境胁迫期间或之后，登革病毒向人群传播的风险可能轻微上升。

全局敏感性分析：PRCC分析显示，R_u、R_w、R_uw三个阈值参数对产卵率、发育率、死亡率及静止率等参数具有不同程度的敏感性。R_uw（衡量感染蚊因母系传递不完全而产生未感染后代的贡献）对沃尔巴克氏体流行率和传播风险的影响最为显著：R_uw增加与感染流行率下降、传播风险上升正相关，凸显了母系传递效率的流行病学重要性。传播风险对流行率的敏感性最高，且与R_w负相关、与R_u和R_uw正相关。

五、讨论与结论

研究讨论指出，静止对蚊群适合度的影响取决于存活与发育进程之间的权衡，未感染蚊群会根据环境线索优化其适合度，在孵化率与静止率之间进行动态调节。对于沃尔巴克氏体感染蚊群，静止的负面影响源于感染卵在静止后存活率低或羽化成蚊不育，这实质性地降低了感染蚊的适合度，削弱了感染在蚊群中的流行率。更重要的是，静止增加了感染建立所需的释放蚊数量，在静止高发地区可能危及沃尔巴克氏体蚊释放策略的有效性。研究结论强调，在制定沃尔巴克氏体蚊释放计划时，必须充分考虑目标地区的卵静止发生频率与环境胁迫条件，否则防控效果可能大打折扣，甚至在环境胁迫期间增加登革传播风险。