盐城市2013—2023年恙虫病流行病学特征分析

[中图分类号]R513.2 [文献标志码] B [DOI] 10.19767/j.cnki.32-1412.2025.05.025

恙虫病又名丛林斑疹伤寒，是由恙虫病立克次体引起的急性传染病，系一种自然疫源性疾病，啮齿类为主要传染源，恙螨幼虫为传播媒介。该病起病急，临床主要表现为高热、毒血症、皮疹、焦痂和淋巴结肿大等，可导致多脏器损害[2-3]。盐城市所辖东台县1986年首次在长江以北地区报告恙虫病，2006年所辖滨海县发生恙虫病暴发疫情，之后将恙虫病纳入监测报告传染病，在2019年报告病例数最多，采取一系列综合防控措施后，疫情逐渐下降。为更好了解盐城市近几年恙虫病流行趋势，为科学防控提供决策依据，现对盐城市2013—2023年恙虫病流行特征进行分析。

1资料与方法

1.1资料来源数据来源于"中国疾病预防控制信息系统”传染病监测信息系统，统计时间为2013年1月1日—2023年12月31日，按发病日期和现住址统计，内容包括地区、年龄、性别、发病时间等传染病报卡信息及发病数、发病率。

1.2分析方法运用Excel2016软件录入、整理数据，SAS9.4统计学软件进行数据分析。三间分布采用描述性流行病学方法进行分析；计数资料以 n（%） 表示，组间比较采用 χ² 检验;应用 Joinpoint Regres-sionProgram5.3.0.0软件拟合回归模型，计算年度变化百分比（annual percent change，APC）、平均年度变化百分比（average annual percent change，AAPC）并分析变化趋势；采用集中度法和圆形分布法描述季节特征，使用Cochran-Armitage检验有序分类变量的线性趋势。 Plt;0.05 为差异有统计学意义。

利用ArcGis10.8软件进行空间自相关分析。采用全局莫兰指数（MoransI）衡量盐城市发病率的整体分布情况，当Morans Igt;0 且 Z 值 ⩾1.96 时为正相关，表明整体呈现聚集性分布；当Morans Ilt;0 且 Z 值 ≤-1.96 时为负相关，表示空间上呈现离散性分布；当MoransI接近0且 Z 值介于-1.96～1.96，说明整体呈随机性分布。MoransI绝对值越大，表明空间相关性越大。全局MoransI可能会忽视盐城市局部聚集性，因此本研究还采用安瑟伦局部莫兰指数（AnselinLocalMoransI）分析局部地区的聚集现象。聚集性可表现为“高-高”“高-低”、“低-高”“低-低"4种聚集模式，高-高聚集表明高值区域被高值区域包围。

2结果

2.1发病概况盐城市2013—2023年累计报告恙虫病5326例，无死亡病例报告，年平均报告发病率为6.84/10万，2019年报告发病率最高为11.63/10万，2013年报告发病率最低为3.50/10万，各年份报告发病率差异有统计学意义（ χ²=685.334，Plt; 0.001）。发病长期趋势直观呈双峰分布，发病率较高的年度分别为2015年、2016年和2019年。见图1。

图1盐城市2013一2023年恙虫病报告发病数及报告发病率

Joinpoint回归模型对恙虫病发病长期趋势分析显示，2013—2023年 AAPC=2.9767%（95%CI -6.9386%，11.7012% ），发病率呈上升趋势，但差异无统计学意义（ P=0.5251 ）;2013—2015年 APC= 60.8831%（95%CI：2.3454%，159.4649% ，发病率呈上升趋势，差异有统计学意义（ P=0.0356 ;2015-2023年 APC=-7.8922%（95%CI：-31.0452%，-0.6919） ，发病率呈下降趋势，差异有统计学意义（ P=0.0420 。见图2。

2.2季节分布2013—2023年盐城市各年恙虫病发病呈明显单峰分布（图3），集中度法分析显示：Rx=-0.922，Ry=0.238 ，集中度 （M）_=0.953 ，说明盐城市恙虫病发病具有严格季节性；圆形分布法分析显示： ΔX=0.463，Y=-0.832 ，雷氏检验结果为 Z= 4830.120，差异有统计学意义 （Plt;0.05） ，即平均角存在 ，对应高峰日为10月31日，高峰期为10月13日—11月18日。

图2盐城市2013—2023年恙虫病报告发病率变化趋势Joinpoint回归分析

图3盐城市恙虫病发病数据季节分解结果

2.3地区分布盐城市各县区均有恙虫病病例报告，年平均报告发病率居前3位依次为响水县（22.79/10万）东台市（11.48/10万）大丰区（8.38/10万）。大市区（亭湖区、盐都区）年平均报告发病率（1.49/10万）低于县区年平均报告发病率（8.48/10万），差异具有统计学意义 （χ²=1003.184，Plt;0.001） ;沿海地带（响水县、滨海县、射阳县、大丰区、东台市）年平均报告发病率（9.42/10万），显著高于非沿海地带（阜宁县、建湖县、盐都区、亭湖区）年平均报告发病率（3.67/10万） 。各县（市、区）发病分布情况见图4、5。

图4盐城市2013一2023年羔虫病报告发病率热点图

图5盐城市2013—2023年恙虫病报告空间分布图

2.4人群分布盐城市2013—2023年共报告恙虫病男性病例2379例，女性病例2947例，男女性别比为 0.81：1 ;女性年平均报告发病率（7.59/10万）高于男性（6.08/10万），差异有统计学意义（ ⋅x²=65.161 ，Plt;0.0001 ）。盐城市2013—2023年各年龄组均有病例报告，各组报告发病率差异有统计学意义（ χ²= 5765.193， Plt;0.0001 ）。报告病例主要集中在50～80岁 （4132/5326，77.58% ），其中以60～70岁年龄组最多 （1739/5326，32.65%） 。发病人群以农民为主（4791/5326，89.95% ），其次为离退人员 （155/5326，2.91% ）、家务及待业人员 （99/5326，1.86% ）。

2.5空间自相关分析全局空间自相关分析结果显示，盐城市2019一2022年恙虫病发病率呈全局空间正相关，2021年MoransI最大，表明该年盐城市恙虫病的发病率聚集程度最强。见表1。AnselinLocalMoransI显示，建湖县2013—2016年呈现高-低聚集，2017年之后盐城市局部呈现低-低聚集。

表1 2013—2023年盐城市恙虫病发病率全局空间自相关分析结果

3讨论

盐城市2013—2023年累计报告恙虫病5326例，年平均报告发病率6.84/10万，高于省内及周边省份[46。Joinpoint回归分析结果显示，盐城市2013—2015年恙虫病发病率持续增长，年度变化百分比达到 60.99% 。我省自2011年起积极推行秸秆禁烧，由于秸秆机械化还田技术尚未完全普及，大量秸秆堆放为媒介生物恙螨的活跃提供了条件，增加农民等高风险群体的接触机会，2015年后恙虫病发病呈下降趋势，APC为 -7.89% 。盐城市2013—2023年恙虫病报告发病率总体呈上升趋势，但差异无统计学意义，与常州市、连云港市等地恙虫病发病趋势一致[-，这可能与2019年报告发病数大幅增加，作为离群值对趋势拟合产生影响有关。

图6盐城市2013—2023年恙虫病AnselinLocalMoransI分析结果

本研究集中度与圆形分布法结果显示，恙虫病发病具严格季节性，发病高峰期在10月13日至11月18日，属秋季型疫区，符合长江以北地区恙虫病流行特点。10月一11月正值秋收农忙时节，农民于田地、野外等劳作过程中接触恙虫病传染源及传播媒介的机会更多；此外，还可能与宿主动物黑线姬鼠及与媒介生物恙螨的数量、密度、活动性的增加有关，与国内其他地区研究类似[12-13]。

本研究结果显示，盐城市2013—2023年恙虫病女性发病率高于男性，发病职业人群以农民为主，发病年龄段以50～80岁为主。其原因在于农村年轻男性多外出务工，女性及老年人留守务农，增加了暴露机会，故呈现女性、农民、老年高发的特征。

盐城市各县（市、区）均有发病，且县区发病率高于大市区、沿海地带发病率高于非沿海地带。这与各县区较大市区更多致力于农业生产活动有关，由于区域农业种植规模化快速发展，土地利用率提高，人员流动加快，利于恙虫病流行。此前研究发现，盐城市恙虫病发病高峰前1至2个月会出现鼠密度升高现象，且盐城市沿海滩涂地区以黑线姬鼠为绝对优势鼠种，为沿海地带恙虫病发病率高于非沿海地带提供佐证。

全局空间自相关分析结果显示，盐城市2019—2022年恙虫病存在空间聚集性，由于2019年响水县报告发病数较多，使得恙虫病明显辐射周边地区，造成发病率的升高以及全局的空间聚集。2013—2023年建湖县以高-低聚集模式转变为低-低聚集模式，表明该地区防控措施取得成效，发病率水平趋近周边地区。

在今后恙虫病的防控中要重点关注秋收季节后疫情监测，强化中老年农民健康教育和个人防护等工作，不断提高基层医疗机构识别、发现恙虫病患者的能力。

[参考文献]

[1]李贵昌，栗冬梅，李燚，等.2006—2016年我国恙虫病流行特征分析[J].疾病监测，2018，33（2）：139-143.

[2]彭培英，徐蕾，王谷仙，等.1952—1989 年和2006—2017年中国大陆恙虫病流行及时空分布特征[J]中国人兽共患病学报，2022，38（9）：818-823.

[3]岳玉娟，王玉姣，李贵昌，等.2006一2018年中国大陆恙虫病高发区流行病学特征分析[J].疾病监测，2020，35（4）：301-306.

[4] TILAK R，KUNTE R. Scrub typhus strikes back：are weready[J].MedJArmed ForcesIndia，2019，75（1）：8-17.

[5] LUO L，GUO Z N，LEI Z，et al. Epidemiology of tsutsuga-mushi disease and itsrelationship with meteorological factorsin Xiamen city，China[J].PLoSNegl TropDis，202O，14（10）：e0008772.

[6]张育富，刘大鹏，褚宏亮.江苏省2008—2017年恙虫病流行特征和时空聚集性分析[J].中国媒介生物学及控制杂志，2019，30（3）：237-243.

[7]陈安喜，孔卫荣，田野，等.江苏省常州市2012—2021年恙虫病流行特征分析[J].中国媒介生物学及控制杂志，2023，34（3）：372-377.

[8]郑艳泽，陈玉，张嘉陵.2010-2022年连云港市恙虫病流行病学特征分析[J].公共卫生与预防医学2024，35（3）：95-98.

[9]邝仕成，吴彪，陈贻彪，等.海口市某医院330例恙虫病流行和临床特征分析［J].中国感染与化疗杂志，2021，21（5）：512-516.

[10]张嘉溪，谭盛葵.恙虫病流行病学研究新进展[J].中国热带医学，2022，22（3）：274-278.

[11]李孟磊，王文华，李懿，等.河南省 2013—2019 年恙虫病流行病学特征分析[J].实用预防医学杂志，2021，28（10）：1184-1186.

[12]薛莉颖，谢晓菲，王高玉，等.海南省680例不明原因发热患者恙虫病东方体感染调查[J].中国热带医学，2022，22（10）：930-934.

[13]马蒙蒙，王大虎，鲁影，等.广州市2015一2019 年恙虫病流行特征分析[J].医学动物防制，2021，37（8）：744-747.

（本文编辑缪宏建）