关键词:
白纹伊蚊
指标
登革热风险
动力学模型
常微分方程
宏病毒组学
摘要:
1 研究背景和目的白纹伊蚊可以传播多种蚊媒传染病,如登革热、寨卡病毒病、基孔肯雅热等等。对白纹伊蚊进行密度消长监测,建立白纹伊蚊种群预测模型和携带病原体检测,可以指导灭蚊工作、评价灭蚊效果、判断疾病暴发风险等,具有重要的公共卫生学意义。2研究方法2017年8月至2018年7月,每月采用布雷图指数(Breteau Index,BI)和诱蚊诱卵指数(Mosquitoand oviposition positive index,MOI)监测广州市白云区的XW村和NF村,从2017年12月开始增加成蚊密度指数(Adult-mosquito Density Index,ADI)加以监测。2014年7月至2015年6月在XW村亦开展了为期1年,每月2次的白纹伊蚊密度调查,该数据也用于模型的验证。收集同期气象数据和白云区各月登革热病例数。基于参考文献和气象数据建立白纹伊蚊8个生活阶段的常微分方程组,现场监测数据用于该模型的验证。另外,选择2017年9月至2018年8月广州市白云区、越秀区和天河区的2个农村和2个城市社区采集的白纹伊蚊样本作为研究对象,按照不同城市和农村、不同季节合样,用于宏病毒组学研究。3研究结果3.1 白纹伊蚊现场监测结果幼虫和成蚊密度消长均受气温影响。两监测点各指数的消长趋势基本一致,三个监测指标之间均存在正相关关系,但大部分月份标化BI都大于标化MOI与标化ADI。2017年8月至9月,BI和MOI指示的登革热风险均为3级,受10月份灭蚊影响,10-11月两者指示的登革热风险并不一致。12月至次年2月,BI指示的登革热风险等级为1-3级,MOI和ADI这三个月反映的登革热风险等级均为0级。3月份后,BI已超过20,指示风险为3级,而XW村和NF村的MOI分别为2.13、3.77,风险为0级,ADI分别为4只/人·小时和6只/人·小时,风险分别为1级和2级,从4月份开始各监测指标反映的登革热风险趋于一致。3.2 白纹伊蚊种群动力学模型结果模拟的白纹伊蚊成蚊主要有两个密度高峰,分别是在6-7月和9-10月,在普遍高温的8月份成蚊密度反而有所下降。模拟相对密度与标化监测数据的平均绝对残差范围是7.40%-11.22%,Pearson相关系数为0.902-0.989,该种群动力学模型总体拟合效果较好。基于模型建立布雷图指数和诱蚊诱卵指数对应关系的参考值范围,与现场数据相比总体符合率为83.33%(40/48)。3.3 白纹伊蚊宏病毒组分析结果4份宏病毒样本共注释到499种病毒,其中农村冬春季样本含有更多的农村冬春季样本含有更多的细长病毒属、正痘病毒属、Negevirus、甲型冠状病毒属、绿藻病毒属、埃可病毒属和青绿藻病毒属,农村夏秋季含有更多的黄病毒属、Phikzlikevirus和Lambdalikevirus,而城市冬春季则含有较多的Yualikevirus、T4likevirus和Phikmvlikevirus。样本层级聚类和主成分分析表明城市和农村的病毒群落结构相似,而冬春季和夏秋季的病毒群落结构之间存在差异。4结论4.1 布雷图指数指示的登革热风险总是高于或等于诱蚊诱卵指数和成蚊密度指数指示的风险。诱蚊诱卵指数和成蚊密度指数能较敏感地反映白纹伊蚊的活动情况,指示登革热风险的准确性高。4.2基于白纹伊蚊生态习性、温度、日照时长和降雨量建立了白纹伊蚊种群动力学模型,模型拟合效果较好。4.3 初步发现城市和农村的病毒群落组成相似,而冬春季和夏秋季的病毒群落组成则存在差异。