蚊虫扩散的行为态景观资源背景分析（Mosquito dispersal in context）

《PLOS Complex Systems》：Mosquito dispersal in context

　　

摘要：蚊虫扩散（mosquito dispersal）在蚊虫生态学和蚊媒病原体传播中起重要作用。传统的反应-扩散（reaction-diffusion）模型和基于斑块（patch-based）模型通常采用简单的迁出假设来描述蚊虫迁出，在模拟蚊虫扩散时普遍忽略了蚊虫行为生态学——尤其是搜寻资源（resource searching）的过程。本研究采用高拟真（highly mimetic）模型分析蚊虫运动，以探讨不同建模思路所能揭示的信息。研究人员构建了两类行为态微模拟（behavioral state microsimulation）模型，其中蚊虫处于不同行为态（behavioral state）并在景观中搜寻所需资源，资源以点集形式分布。此类模型构建与分析较为繁琐，因此研究人员开发了 ramp.micro——一个用于在 R 语言中构建、求解、分析及可视化蚊虫行为态微模拟模型的程序包。蚊虫种群扩散是该类复杂系统的涌现（emergent）特征，研究人员据此发展了描述与理解它的新方法。研究结果表明：若不对处于行为态、为完成特定目标而进行连续飞行（sequences of flights）的蚊虫加以考量，仅对观测到的蚊虫运动模式做粗略分析无法揭示种群结构；即便资源呈随机均匀分布，蚊虫种群仍倾向于形成高度空间结构化群落（spatially structured communities）；蚊虫种群密度的部分异质性可归因于由搜寻行为和资源空间分布所定义的网络特征。上述模型凸显了理解蚊虫行为态和资源可获得性对认识蚊虫种群扩散的重要性，并指出局地背景（local context）的关键作用。受此动力学启发，考虑资源可获得性作为影响蚊虫运动和扩散的因素，将有助于改进蚊虫生态学及蚊媒病原体传播的空间模型。

论文解读——《Mosquito dispersal in context》（发表于《PLOS Complex Systems》）

一、研究背景与立项依据

传统上，蚊虫种群扩散多采用反应-扩散（reaction-diffusion）方程或基于斑块的元种群（metapopulation）模型描述，这类唯象（phenomenological）扩散模型忽略蚊虫生物学细节，特别是蚊虫按生理需求在吸血（blood feeding）、产卵（egg laying）和取糖（sugar feeding）间切换的行为态（behavioral state），以及为获取资源而进行的意图性飞行（intentive flight bouts）。标记—释放—重捕（mark-release-recapture, MRR）研究虽可提供经验数据，但难以解析资源共分布对种群空间结构的塑造作用。由于蚊虫实际运动是依行为态在异质资源点集间做序列搜寻，扩散假设可能误导对蚊媒疾病传播范围的估计。为此，研究人员提出基于行为态微模拟（behavioral state microsimulation）的框架，将扩散视为资源搜寻过程的涌现属性，并结合网络分析方法量化其结构。

二、主要关键技术方法

研究人员建立两类离散时间成虫行为态微模拟模型——BQ模型（仅含吸血态B与产卵态Q）和BQS模型（增加取糖态S），将血餐宿主栖息点（haunts，集合记为?）、产卵水生生境（aquatic habitats，集合记为??）和糖源（sugar sources，集合记为??）定义为平面上随机均匀分布的独立点集。蚊虫每日以一定存活概率（p）尝试搜寻当前行为态对应资源，成功则转入下一行为态并按距离加权矩阵（dispersal matrix M，元素由指数衰减权重归一化得到）迁往目标资源点；失败可经参数（φ_BS、φ_QS、φ_QB等）转入取糖态再返吸血态。幼期蚊（immature mosquito population）采用密度依赖的孵化—成熟—羽化模块，雌性成蚊羽化后加入成虫B或S态。研究人员开发R语言程序包 ramp.micro（https://github.com/dd-harp/ramp.micro）实现模型构建、迭代求解、净扩散核（net dispersal kernel，如单次任务完成矩阵 D^B、D^Q，完整生殖营养周期矩阵 K²=D^QD^B或 K²=D^BD^Q）计算、终生卵扩散核G与潜在病原扩散核V（考虑外潜伏期EIP）推导，并将单资源点集上的扩散矩阵转为有向加权图，采用iGraph包walktrap算法做群落检测（community detection），对比原始扩散矩阵与周期闭合矩阵所揭示的空间结构。

三、研究结果

Population dynamics（种群动态）

在256个吸血点、289个产卵点（BQ模型）及额外225个取糖点（BQS模型）随机均匀分布的景观中，各水生生境幼期参数一致，但稳态（steady state）下各资源点处成虫密度呈显著异质分布，证明仅资源点的空间共分布与搜寻规则即可引致种群密度异质性。

Dispersal（扩散）

定义了单飞行段资源间转移矩阵M_?→??、M_??→?、M_?→??等，列和为1描述存活蚊虫目的地概率分布，死亡归入逐日存活参数p。通过迭代可推演完成单次吸血或产卵任务（含多次失败重试）的净转移矩阵D^B（从??出发最终成功抵?吸血）与D^Q（从?出发最终成功抵??产卵），以及完整生殖营养周期在两资源点集内部闭合的运动矩阵K²_?=D^QD^B（血餐后返血餐点）和K²_??=D^BD^Q（产卵后返产卵点）。

Movement to complete a task（完成任务的运动）

计算表明D^B与D^Q构成二分图（bipartite graph），经一次完整吸血—产卵循环后二分图解耦为各自点集上的独立图（K²矩阵），该结构更适于揭示真实种群连通性，而非直接用单步混合矩阵分析。

Movement through a feeding cycle（通过一个取食周期的运动）

块矩阵形式给出完整周期转移：[[0,D^B];[D^Q,0]]，周期结束后得K²_?与K²_??，反映蚊虫在同一资源类型点间的"回返"扩散模式。

Egg laying（产卵扩散）

终生卵扩散核 G = Σ_io·p^g·i·(K²_??)^i-1D^Q（o为每批产卵数，p_g为日成蚊存活率），种群水平净卵扩散 ē = Σ_jN_?,j·G 显示卵并非均匀散布，而依网络拓扑向特定产卵点聚集。

Pathogen dispersal（病原体扩散）

模拟初始吸血队列B₀，累加EIP（extrinsic incubation period，外潜伏期）日后所有传染叮咬得潜在传播矩阵V，加权稳态密度得潜在病原扩散（potential pathogen dispersion）V? = diag(N_?)·V，某些吸血点因拓扑位置接收更多传染叮咬。

Dispersal distances（扩散距离）

基于资源点间距与转移概率质量函数（PMF）可定义扩散距离的概率分布。

Population structure（种群结构）

对K²_?与K²_??做walktrap群落检测发现明显社区（communities）；而对原始单步矩阵（式13）做同样分析所得群落与生物学意义不符，说明忽视行为态序列会掩盖真实结构——粗析观测运动模式若不计入行为态将无法揭示社区结构。

Communities（社区）

即便资源随机均匀分布，因近距离多资源共存引致正反馈，蚊虫自组织形成社区，社区内连接紧密、社区间弱连接；卵扩散与病原扩散的社区在空间上重叠。这些社区可作元种群斑块，且较任意平面剖分更能近似底层过程。

Flows（流）

资源相对邻近度造成转移矩阵不对称，产生由缺资源区向外、全资源区向内的净流（flow）；不对称流强化而非破坏社区凝聚力，高密度点位于社区核心、低密度点位于边缘。

Context and pattern（背景与格局）

限制某类资源（如仅留13个人血宿主点模拟农村焦点、仅13个水生生境模拟城市焦点、仅13个糖源）时，蚊虫聚集于最稀缺资源邻近区域，最稀缺资源主导种群结构；糖为可选项时限制效应较弱。

四、讨论与结论翻译

研究人员通过模拟基本蚊虫行为阐明蚊虫运动如何被资源空间共分布（context）塑造；种群结构由必需资源共分布与蚊虫搜寻行为共同涌现，导致成虫分布异质性与斑块化。显著量的细粒度种群异质性与结构源于在简单异质资源景观上施加简单搜寻与扩散规则。斑块围绕具全部必需资源的位点自组织形成；反之缺任一必需资源的景观位置因蚊虫离寻而低密度。资源被利用的可能性取决于其与其他必需资源的邻近度。蚊虫群落（mosquito communities）——定义为描述资源间蚊虫运动的网络中具连通性的节点组——是类元种群斑块的自组织空间现象，社区内有资源间大量往返运动，社区间为远距离扩散。与元种群模型不同，此处斑块是涌现特征而非人为划定。行为态微模拟虽复杂但 ramp.micro 可降低分析成本并提高可复现性。此框架足够真实以定量界定蚊虫基础生态位（fundamental niche）——依赖行为算法、水生生境质量稳定性、其他资源可得性及搜寻代价。景观异质资源分布（尤其一种或多种必需资源相对匮乏）引致复杂结构化异质蚊虫种群。该高拟真模型提示需重新审视MRR估算蚊虫种群规模的公式偏差，并可用于遗传改造蚊株释放扩散、监测诱蚊器布设及传病热点（focus）识别等应用。综上，研究人员得出结论：蚊虫种群及潜在传播倾向于聚集于具备全部必需资源之处；异质资源分布——特别是至少一种必需资源的相对稀缺——造就复杂结构化异质蚊虫种群，考虑资源可得性的空间模型有助于更精准靶向控蚊与阻断疾病传播。

（注：全文术语首次出现均附英文原词或缩写，上标如^B、下标如_?按原文保留，未引入原文未述之推测内容。）