随机Wolbachia动力学：持久性、不变测度与基于概率的释放策略

《Axioms》：Stochastic Wolbachia Dynamics: Persistence, Invariant Measures, and Probability-Based Release Strategies

　　

研究人员在环境变异性下建立了Wolbachia在蚊群中入侵的随机离散时间框架。从具有Allee型阈值的确定性频率依赖差分方程出发，考虑了随机迭代X_n+1=F_Θn(X_n)，其中Θ_n代表适合度代价、细胞质不相容性和母系传播的波动。研究人员证明了存在性、唯一性、正性、[0,1]的前向不变性、连续性以及路径序保持，并将该模型表述为随机动力系统(random dynamical system, RDS)。在独立环境强迫下，该过程具有马尔可夫性(Markov性)；相应的马尔可夫算子(Markov算子)具有Feller性质，状态空间的紧性蕴含不变概率测度(invariant probability measure)的存在性。灭绝附近的局部行为由随机李雅普诺夫指数(Lyapunov指数)λ=E[log((1-M)(1-s_f))]刻画，λ<0时无Wolbachia状态具有局部指数稳定性。研究人员进一步引入了有限时域建立概率(finite-horizon establishment probability)和基于概率的释放阈值(probability-based release threshold)，其单调性由路径比较原理得出。数值模拟展示了环境变异性的大小、来源和依赖结构如何影响入侵概率和释放需求。该框架将非线性差分方程、随机动力系统、马尔可夫算子、不变测度理论和随机稳定性统一起来，对阈值依赖的Wolbachia入侵进行了综合分析。

**随机环境下的Wolbachia入侵：统一随机框架的建立与概率释放策略——解读《Axioms》最新研究**

**一、研究背景与动机**

蚊媒传染病（登革热、寨卡、基孔肯雅热等）仍是全球重大公共卫生挑战。快速城市化、气候变化和杀虫剂抗性演化使传统化学防治效果持续下降，Wolbachia作为自然存在于多种昆虫但缺失于埃及伊蚊(Aedes aegypti)的胞内共生菌，通过细胞质不相容性(cytoplasmic incompatibility, CI)修饰蚊虫生殖并降低其媒介能力，已在多国大规模现场部署中展现显著控病效果。

确定性频率模型表明，Wolbachia入侵受三个生物机制共同主导：感染适合度代价(s_f)、CI强度(s_h)和母系传播保真度(1-M)。在生物学合理参数下，该类模型呈现双稳结构，一个不稳定中间平衡态构成Allee型入侵阈值——低于阈值感染趋于灭绝，高于阈值实现种群替换。然而，确定性模型假设参数在入侵全程恒定。现实中蚊群经历季节温度、降雨、湿度、幼虫资源、释放不完整等环境波动，可同时影响s_f、s_h和M，导致相同初始释放频率在不同环境下产生不同结局。已有随机Wolbachia研究多聚焦特定噪声形式（如随机微分方程或特定扰动），缺乏一个能容纳一般环境过程、统一整合随机动力系统、马尔可夫算子、不变概率测度与概率释放策略的离散时间分析框架。本研究即为此目的而开展。

**二、研究内容与主要贡献**

研究人员将确定性参数向量替换为平稳随机环境过程Θ_n=(s_f, s_h, M)，建立随机离散时间模型X_n+1=F_Θn(X_n)，该表述自然容纳独立同分布、时间相关与遍历平稳环境。主要贡献包括：证明随机模型解的存在唯一性、正性、[0,1]前向不变性、连续性与路径比较原理；发展随机演化算子框架；在独立环境强迫下证明马尔可夫性、Feller性质及不变概率测度存在性；通过随机李雅普诺夫指数刻画灭绝局部稳定性；定义有限时域建立概率与q水平释放阈值并证明其单调性。研究发表于《Axioms》。

**三、关键技术方法**

研究采用随机动力系统理论（随机演化算子与逐路径分析）、Feller-马尔可夫链理论（转移核与马尔可夫算子）、Krylov-Bogolyubov定理（紧状态空间上不变概率测度存在性）、随机李雅普诺夫指数与遍历定理，以及蒙特卡洛数值模拟。数值实验系统考察环境变异性大小、来源及参数相关结构对入侵概率和释放阈值的影响。研究未使用现场样本队列数据，模拟参数基于生物学合理设置而非拟合特定现场试验。

**四、主要研究结果**

**4.1 路径wise动力学性质（第4节）**：定理1证明在生物适定性(H1)与可测性(H2)假设下，随机递归对任意可行初值存在唯一解，且状态空间前向不变。命题2证明对每个环境状态θ，映射F_θ连续且单调不减，故随机演化映射亦连续保序。路径比较原理（推论1）保证相同环境实现下初始频率的序关系保持，零状态为吸收态，完全感染状态在母系传播不完美时非不变。

**4.2 马尔可夫动力学与长期统计行为（第5节）**：在环境过程独立同分布假设下，随机递归定义时齐马尔可夫链；转移核与马尔可夫算子给出单代演化的概率描述。定理3证明马尔可夫算子的Feller性质，定理4经由Krylov-Bogolyubov定理保证不变概率测度存在。由于F_θ(0)=0恒成立，灭绝测度δ₀总是不变测度，故存在性不蕴含唯一性或正支撑。Frobenius-Perron方程给出了绝对连续不变密度的积分刻画。

**4.3 持久性、灭绝与概率释放阈值（第6节）**：定理5证明在平稳遍历环境下，若随机李雅普诺夫指数λ=E[log((1-M)(1-s_f))]<0，则灭绝状态局部指数稳定几乎必然成立。由于λ≤0恒成立，从极低频率的入侵无法自发性增长，成功建立依赖有限规模释放以越过Allee障碍。研究定义有限时域建立概率P_L,η^est(x)（感染频率在规划时域最后L代维持在目标η以上），并证明其对初始频率单调不减；由此定义q水平释放阈值x_q^rel=inf{x: P_L,η^est(x)≥q}。无随机性且传播完美时该阈值退化为确定性Allee阈值。

**4.4 数值模拟（第7节）**：六个数值实验揭示：（1）环境随机性将尖锐确定性阈值转变为概率过渡区，阈值附近释放对波动最敏感；（2）释放阈值随环境变异性增大而升高，高目标概率下更显著；（3）参数敏感性分析显示s_f波动影响最强，s_h与M影响较弱，三者同时波动效应最大；（4）s_f与M的正相关升高释放阈值；（5）长期经验分布呈现阈值依赖的随机双稳现象：低初始频率集中于灭绝状态，高初始频率集中于高感染状态附近。

**五、讨论与结论**

讨论部分强调，该框架与确定性模型的根本区别在于将管理问题从"释放是否超过固定阈值"转变为"需多大释放方能以规定概率实现成功建立"。生物学上，热应激可降低Wolbachia密度、削弱CI并降低母系传播保真度，wMel感染埃及伊蚊的现场实验和实验室实验均证实了温度对上述机制的显著影响，这为将s_f、s_h和M视为环境变化量提供了实证支持。该框架区别于SIR/SIRS水平病原传播模型及生物膜空间模型，其可迁移性在于随机环境方法本身而非Wolbachia特异的递推式。自然扩展方向包括时间相关环境、空间扩散、重复或自适应释放及随机积分差分方程模型。

研究结论指出：环境随机性并非简单地扰动确定性入侵阈值，而是改变了相关管理问题的本质——从判断释放是否超过固定阈值，转变为确定需要多大释放规模才能以规定概率实现成功建立。该研究为阈值依赖的Wolbachia入侵提供了一种将非线性差分方程、随机动力系统、马尔可夫算子、不变测度理论与随机稳定性整合于一体的统一分析框架，其方法论亦可推广至更广泛的离散时间生态与进化系统。