2016—2020年青藏高原东部地区棘球蚴病空间聚集性分析

施丹丹，蒉嫣，韩帅，伍卫平

中国疾病预防控制中心寄生虫病预防控制所（国家热带病研究中心），国家卫生健康委员会寄生虫病原与媒介生物学重点实验室，世界卫生组织热带病合作中心，国家级热带病国际联合研究中心，上海 200025

棘球蚴病（俗称包虫病）是棘球绦虫的幼虫（棘球蚴）引起的人兽共患寄生虫病，在全球范围内广泛分布[1]。棘球蚴病严重威胁人群健康，是重要的公共卫生问题之一，同时给畜牧业造成了巨大经济损失。我国是囊型和泡型棘球蚴病的流行区，受威胁人口数和患病人数均为全球首位[2-3]。青藏高原东部地区位于青海省和四川省交界处，包括青海省玉树州、果洛州和四川省甘孜州、阿坝州，各县（市）均为棘球蚴病的流行区，且流行程度较重。2012年调查数据显示，青海省人群棘球蚴病患病率为1.10%，其中果洛州为5.20%，玉树州为4.54%[4]；四川省人群棘球蚴病患病率为1.08%，其中甘孜州为1.86%，阿坝州为0.79%[5]。近年来我国加大了棘球蚴病防治力度，并取得了明显成效，但局部地区的患病率仍然较高，存在传播风险[6]。

目前，空间统计分析方法在棘球蚴病流行病学研究中应用广泛，有助于评估患病率差别、识别疾病分布特征和聚集性、分析疾病影响因素[7-8]。了解和掌握棘球蚴病高流行地区的空间分布特征，对棘球蚴病防治策略的制定和资源的合理分配有重要参考意义。本研究通过分析2016—2020 年青藏高原东部地区人群棘球蚴病发病率空间分布特征，探索其空间聚集性，为棘球蚴病的防治提供参考。

1 资料与方法

1.1 研究区域 选取青藏高原东部地区，包括四川省的甘孜州和阿坝州、青海省的玉树州和果洛州共计43个县（市）作为研究区域，见图1。

1.2 数据来源 2016—2020年研究区域的棘球蚴病发病数据来源于2020 年全国传染病监测报告数据库。监测方法：按照各省棘球蚴病流行种类及流行程度的比例，计算抽取各类流行村的数量，采用随机抽样的方法抽取监测点；在监测点选择一个或相邻的几个村的全部居民作为监测人群，每年进行超声检查，并以血清学方法检测抗棘球蚴抗体进行辅助诊断[9]。2016—2018年青藏高原东部地区仅有28 个县（市）开展了监测工作，2019—2020 年43 个县（市）均开展了监测工作。

1.3 空间自相关分析 以县（市）为单位，对2016—2020 年青藏高原东部地区43 个县（市）棘球蚴病发病率进行全局空间自相关分析、局部空间自相关分析和热点分析，2016—2018年未进行监测工作的县（市）数据以缺失值来进行分析。其中空间自相关是指发病率在空间位置上的相关性，可反映发病率在研究县（市）的聚集程度。若发病率在空间上呈现中心与周围分布一致，则称为空间正相关；若呈现中心与周围分布相反，则称为空间负相关；若发病率在空间分布上呈现随机性，则表明空间相关性不显著[10]。局部空间自相关分析选用Local Moran′sI统计量，利用LISA（local indications of spatial association）来进行分析，变量在空间上存在聚集区域的类型，可分为四种：“高-高”聚集（H-H），表示某病发病率高于均值的地区被高于均值的领域包围；“高-低”聚集（H-L），表示某病发病率高于均值的地区被低于均值的领域包围；“低-低”聚集（L-L），表示某病发病率低于均值的地区被低于均值的领域包围；“低-高”聚集（L-H），表示某病发病率低于均值的地区被高于均值的领域包围[11]。采用空间局部Getis-Ord Gi*热点分析发现变量在空间分布上存在热点地区（高值聚集区）和冷点地区（低值聚集区）[12]。

1.4 统计学分析 运用Excel 2007 软件建立数据库。应用SPSS 19.0 对数据进行统计分析，对棘球蚴病发病率进行时间趋势分析。利用ArcGIS 10.1软件进行全局空间自相关分析、局部空间自相关分析和热点分析，并对结果进行可视化展示。检验水准采用α=0.05。

2 结 果

2.1 棘球蚴病发病率 2016—2020年青藏高原东部地区人群棘球蚴病累计新发病例598 例，年平均发病率为0.032%，各年发病率为0.013%～0.553%。其中，2016 年发病率最高，为0.553%；2020 年最低，为0.013%，棘球蚴病发病率随时间推移逐渐下降（χ2=802.991，P＜0.05）。见表1。

表1 2016—2020年青藏高原东部地区棘球蚴病发病率

2.2 棘球蚴病发病率地区分布 新发病例发病率＞0.1%的县（市），2016年为青海省曲麻莱县、久治县、玛沁县、玉树市、达日县、称多县和四川省石渠县、白玉县、甘孜县、德格县、红原县、理塘县；2017年为青海省达日县、玛沁县、曲麻莱县、杂多县、称多县和四川省白玉县；2018 年为青海省称多县、杂多县和四川省德格县、马尔康市、白玉县、道孚县；2019年为青海省治多县、杂多县和四川省阿坝县；2020年仅为青海省久治县和四川省红原县。见图2。

2.3 棘球蚴病发病率空间自相关分析

2.3.1 全局空间自相关分析 2016—2018 年全局Moran′sI值分别为0.005、0.088、-0.046（Z=0.281、1.196、-0.046，P均＞0.05），表明2016—2018 年研究区域棘球蚴病发病率在空间分布上不具有空间自相关性，可能是随机分布。2019、2020 年全局Moran′sI值为0.179、0.166（Z=2.273、2.316，P均＜0.05），表明棘球蚴病发病率呈空间正相关分布，存在高值与高值相邻或低值与低值相邻的聚集现象。

2.3.2 局部空间自相关分析 2019—2020 年青藏高原东部地区棘球蚴病发病率局部空间自相关分析结果显示，2019 年发病率呈现“高-低”和“低-高”聚集，四川省德格县、道孚县、泸定县呈“高-低”聚集，青海省玉树市呈“低-高”聚集；2020年发病率呈现“高-高”“高-低”和“低-高”聚集，四川省阿坝县、红原县、马尔康市和青海省治多县呈“高-高”聚集，四川省白玉县呈“高-低”聚集，四川省若尔盖县、松潘县、理县呈“低-高”聚集。见图3。

2.4 棘球蚴病发病率热点分析 空间局部Getis-Ord Gi*热点分析结果显示，2019年热点99%置信区间地区包含青海省治多县和杂多县，95%置信区间地区包含青海省玉树市；2020年热点99%置信区间地区包含四川省阿坝县、红原县、理县，95%置信区间地区包含四川省若尔盖县、松潘县、黑水县、马尔康市。见图4。

3 讨论

空间流行病学旨在运用地理信息系统等空间信息技术研究疾病分布特征和关键影响因素，为疾病预防控制工作提供策略支持[13]。随着现代空间信息技术的快速发展，空间流行病学正越来越多地被应用于各种疾病研究中，在公共卫生领域发挥着重要作用[14-15]。

2006 年以来，全国逐步开始实施地方棘球蚴病防治项目，通过采取“以控制传染源为主，中间宿主防控与病人查治相结合”的综合防治策略，棘球蚴病的防治工作取得了一定成效[6]。本研究发现，在青藏高原东部地区青海省人群棘球蚴病发病率略高于四川省，这与之前研究结果一致[2]。调查的青藏高原东部地区43 个县（市）2016—2020 年发病率呈下降趋势，2020 年发病率已低至0.013%，说明近年来棘球蚴病防治工作取得了显著成效。棘球蚴病流行趋势趋于平稳，但部分县偶有回升，出现个别新发病例，提示仍需加强对棘球蚴病的动态监测。同时，2016—2020 年棘球蚴病发病率＞0.1%的地区也逐年减少，2016—2018年青藏高原东部地区棘球蚴病发病率较高的地区主要分布在中部交界地区，包括青海省曲麻莱县、称多县、玛多县和达日县，四川省石渠县、德格县、白玉县和阿坝县；2019年发病率较高地区为青海省治多县和杂多县，四川省阿坝县；2020年发病率较高地区为青海省久治县和四川省红原县。可能由于2016年开始实施《全国包虫病监测方案》，部分县刚开展人群监测工作，新发现的病例中可能包含了既往发病但未及时检查出来的病例。2017年青海省达日县新发现121例病例，可能由于当时基层专业人员能力有限，该县2016年未完成全部监测工作，部分监测任务2017年才得以完成。

大量研究显示，棘球蚴病的分布具有一定的空间聚集性。Schmidberger等[16]对1992—2016年德国棘球蚴病数据进行空间自相关性分析，结果显示整个德国的病例分布具有聚集性，局部空间自相关分析结果表明巴登-符腾堡州东南部地区和巴戈利亚州西南部地区为高危区。He 等[17]对四川省2016—2019年人群棘球蚴病流行情况进行调查，全局空间自相关分析显示存在空间正相关，局部空间自相关分析发现“高-高”聚集区在四川省西北部，“低-低”聚集区在四川省东南部。赵月等[18]对四川省甘孜州2014年棘球蚴病的空间分析结果显示棘球蚴病空间分布存在聚集性。本研究发现2016—2018 年青藏高原东部地区棘球蚴病发病率不具有空间聚集趋势，2019—2020年存在空间正相关，具有明显的空间聚集特征，存在“高-高”“高-低”和“低-高”聚集区，可能由于2016—2018年研究区域的筛查地区不完整，有15个县（市）的数据缺失。局部Getis′G热点分析结果与局部空间自相关分析的“高-低”和“高-高”聚集地区较为一致，2019年热点区域位于青海省治多县、杂多县、玉树市，2020年热点区域位于四川省阿坝县、红原县、理县、若尔盖县、松潘县、黑水县、马尔康市。可能由于这些地区位于青南高原，寒冷潮湿的气候环境有利于棘球绦虫虫卵存活，且多房棘球绦虫的中间宿主啮齿类动物数量多、分布广。此外这些地区以畜牧业为主，自然环境和牲畜类型都有利于棘球绦虫完成生活史循环，造成棘球蚴病的流行。甘孜州作为棘球蚴病的重点防控区域，发病率逐渐降低，而阿坝州近年来出现个别新发病例，这可能是热点区域转移的主要原因。

综上，青藏高原东部地区棘球蚴病防治工作取得了一定成效，但仍存在空间分布聚集性。对于“高-低”和“高-高”聚集区域和热点区域，应强化棘球蚴病的防治工作，加强犬只管理，开展健康宣传，同时在既往人群筛查的基础上，加强重点区域的全人群普查。此外，本研究以县级为单位，研究结果仅可反映县级聚集情况，后续可从更小的尺度进行研究；另外，此次研究不涉及影响因素，今后可开展相关研究。