鹽城市2013—2023年恙蟲病流行病學(xué)特征分析

[中圖分類號(hào)]R513.2 [文獻(xiàn)標(biāo)志碼] B [DOI] 10.19767/j.cnki.32-1412.2025.05.025

恙蟲病又名叢林斑疹傷寒，是由恙蟲病立克次體引起的急性傳染病，系一種自然疫源性疾病，嚙齒類為主要傳染源，恙螨幼蟲為傳播媒介。該病起病急，臨床主要表現(xiàn)為高熱、毒血癥、皮疹、焦痂和淋巴結(jié)腫大等，可導(dǎo)致多臟器損害[2-3]。鹽城市所轄東臺(tái)縣1986年首次在長(zhǎng)江以北地區(qū)報(bào)告恙蟲病，2006年所轄濱海縣發(fā)生恙蟲病暴發(fā)疫情，之后將恙蟲病納入監(jiān)測(cè)報(bào)告?zhèn)魅静。?019年報(bào)告病例數(shù)最多，采取一系列綜合防控措施后，疫情逐漸下降。為更好了解鹽城市近幾年恙蟲病流行趨勢(shì)，為科學(xué)防控提供決策依據(jù)，現(xiàn)對(duì)鹽城市2013—2023年恙蟲病流行特征進(jìn)行分析。

1資料與方法

1.1資料來(lái)源數(shù)據(jù)來(lái)源于\"中國(guó)疾病預(yù)防控制信息系統(tǒng)”傳染病監(jiān)測(cè)信息系統(tǒng)，統(tǒng)計(jì)時(shí)間為2013年1月1日—2023年12月31日，按發(fā)病日期和現(xiàn)住址統(tǒng)計(jì)，內(nèi)容包括地區(qū)、年齡、性別、發(fā)病時(shí)間等傳染病報(bào)卡信息及發(fā)病數(shù)、發(fā)病率。

1.2分析方法運(yùn)用Excel2016軟件錄入、整理數(shù)據(jù)，SAS9.4統(tǒng)計(jì)學(xué)軟件進(jìn)行數(shù)據(jù)分析。三間分布采用描述性流行病學(xué)方法進(jìn)行分析；計(jì)數(shù)資料以 n（%） 表示，組間比較采用 χ² 檢驗(yàn);應(yīng)用 Joinpoint Regres-sionProgram5.3.0.0軟件擬合回歸模型，計(jì)算年度變化百分比（annual percent change，APC）、平均年度變化百分比（average annual percent change，AAPC）并分析變化趨勢(shì)；采用集中度法和圓形分布法描述季節(jié)特征，使用Cochran-Armitage檢驗(yàn)有序分類變量的線性趨勢(shì)。 Plt;0.05 為差異有統(tǒng)計(jì)學(xué)意義。

利用ArcGis10.8軟件進(jìn)行空間自相關(guān)分析。采用全局莫蘭指數(shù)（Moran'sI）衡量鹽城市發(fā)病率的整體分布情況，當(dāng)Moran's Igt;0 且 Z 值 ?1.96 時(shí)為正相關(guān)，表明整體呈現(xiàn)聚集性分布；當(dāng)Moran's Ilt;0 且 Z 值 ≤-1.96 時(shí)為負(fù)相關(guān)，表示空間上呈現(xiàn)離散性分布；當(dāng)Moran'sI接近0且 Z 值介于-1.96～1.96，說(shuō)明整體呈隨機(jī)性分布。Moran'sI絕對(duì)值越大，表明空間相關(guān)性越大。全局Moran'sI可能會(huì)忽視鹽城市局部聚集性，因此本研究還采用安瑟倫局部莫蘭指數(shù)（AnselinLocalMoran'sI）分析局部地區(qū)的聚集現(xiàn)象。聚集性可表現(xiàn)為“高-高”“高-低”、“低-高”“低-低\"4種聚集模式，高-高聚集表明高值區(qū)域被高值區(qū)域包圍。

2結(jié)果

2.1發(fā)病概況鹽城市2013—2023年累計(jì)報(bào)告恙蟲病5326例，無(wú)死亡病例報(bào)告，年平均報(bào)告發(fā)病率為6.84/10萬(wàn)，2019年報(bào)告發(fā)病率最高為11.63/10萬(wàn)，2013年報(bào)告發(fā)病率最低為3.50/10萬(wàn)，各年份報(bào)告發(fā)病率差異有統(tǒng)計(jì)學(xué)意義（ χ²=685.334，Plt; 0.001）。發(fā)病長(zhǎng)期趨勢(shì)直觀呈雙峰分布，發(fā)病率較高的年度分別為2015年、2016年和2019年。見(jiàn)圖1。

圖1鹽城市2013一2023年恙蟲病報(bào)告發(fā)病數(shù)及報(bào)告發(fā)病率

Joinpoint回歸模型對(duì)恙蟲病發(fā)病長(zhǎng)期趨勢(shì)分析顯示，2013—2023年 AAPC=2.9767%（95%CI -6.9386%，11.7012% ），發(fā)病率呈上升趨勢(shì)，但差異無(wú)統(tǒng)計(jì)學(xué)意義（ P=0.5251 ）;2013—2015年 APC= 60.8831%（95%CI：2.3454%，159.4649% ，發(fā)病率呈上升趨勢(shì)，差異有統(tǒng)計(jì)學(xué)意義（ P=0.0356 ;2015-2023年 APC=-7.8922%（95%CI：-31.0452%，-0.6919） ，發(fā)病率呈下降趨勢(shì)，差異有統(tǒng)計(jì)學(xué)意義（ P=0.0420 。見(jiàn)圖2。

2.2季節(jié)分布2013—2023年鹽城市各年恙蟲病發(fā)病呈明顯單峰分布（圖3），集中度法分析顯示：Rx=-0.922，Ry=0.238 ，集中度 （M）_=0.953 ，說(shuō)明鹽城市恙蟲病發(fā)病具有嚴(yán)格季節(jié)性；圓形分布法分析顯示： ΔX=0.463，Y=-0.832 ，雷氏檢驗(yàn)結(jié)果為 Z= 4830.120，差異有統(tǒng)計(jì)學(xué)意義 （Plt;0.05） ，即平均角存在 ，對(duì)應(yīng)高峰日為10月31日，高峰期為10月13日—11月18日。

圖2鹽城市2013—2023年恙蟲病報(bào)告發(fā)病率變化趨勢(shì)Joinpoint回歸分析

圖3鹽城市恙蟲病發(fā)病數(shù)據(jù)季節(jié)分解結(jié)果

2.3地區(qū)分布鹽城市各縣區(qū)均有恙蟲病病例報(bào)告，年平均報(bào)告發(fā)病率居前3位依次為響水縣（22.79/10萬(wàn)）東臺(tái)市（11.48/10萬(wàn)）大豐區(qū)（8.38/10萬(wàn)）。大市區(qū)（亭湖區(qū)、鹽都區(qū)）年平均報(bào)告發(fā)病率（1.49/10萬(wàn)）低于縣區(qū)年平均報(bào)告發(fā)病率（8.48/10萬(wàn)），差異具有統(tǒng)計(jì)學(xué)意義 （χ²=1003.184，Plt;0.001） ;沿海地帶（響水縣、濱海縣、射陽(yáng)縣、大豐區(qū)、東臺(tái)市）年平均報(bào)告發(fā)病率（9.42/10萬(wàn)），顯著高于非沿海地帶（阜寧縣、建湖縣、鹽都區(qū)、亭湖區(qū)）年平均報(bào)告發(fā)病率（3.67/10萬(wàn)） 。各縣（市、區(qū)）發(fā)病分布情況見(jiàn)圖4、5。

圖4鹽城市2013一2023年羔蟲病報(bào)告發(fā)病率熱點(diǎn)圖

圖5鹽城市2013—2023年恙蟲病報(bào)告空間分布圖

2.4人群分布鹽城市2013—2023年共報(bào)告恙蟲病男性病例2379例，女性病例2947例，男女性別比為 0.81：1 ;女性年平均報(bào)告發(fā)病率（7.59/10萬(wàn)）高于男性（6.08/10萬(wàn)），差異有統(tǒng)計(jì)學(xué)意義（ ?x²=65.161 ，Plt;0.0001 ）。鹽城市2013—2023年各年齡組均有病例報(bào)告，各組報(bào)告發(fā)病率差異有統(tǒng)計(jì)學(xué)意義（ χ²= 5765.193， Plt;0.0001 ）。報(bào)告病例主要集中在50～80歲 （4132/5326，77.58% ），其中以60～70歲年齡組最多 （1739/5326，32.65%） 。發(fā)病人群以農(nóng)民為主（4791/5326，89.95% ），其次為離退人員 （155/5326，2.91% ）、家務(wù)及待業(yè)人員 （99/5326，1.86% ）。

2.5空間自相關(guān)分析全局空間自相關(guān)分析結(jié)果顯示，鹽城市2019一2022年恙蟲病發(fā)病率呈全局空間正相關(guān)，2021年Moran'sI最大，表明該年鹽城市恙蟲病的發(fā)病率聚集程度最強(qiáng)。見(jiàn)表1。AnselinLocalMoran'sI顯示，建湖縣2013—2016年呈現(xiàn)高-低聚集，2017年之后鹽城市局部呈現(xiàn)低-低聚集。

表1 2013—2023年鹽城市恙蟲病發(fā)病率全局空間自相關(guān)分析結(jié)果

3討論

鹽城市2013—2023年累計(jì)報(bào)告恙蟲病5326例，年平均報(bào)告發(fā)病率6.84/10萬(wàn)，高于省內(nèi)及周邊省份[46。Joinpoint回歸分析結(jié)果顯示，鹽城市2013—2015年恙蟲病發(fā)病率持續(xù)增長(zhǎng)，年度變化百分比達(dá)到 60.99% 。我省自2011年起積極推行秸稈禁燒，由于秸稈機(jī)械化還田技術(shù)尚未完全普及，大量秸稈堆放為媒介生物恙螨的活躍提供了條件，增加農(nóng)民等高風(fēng)險(xiǎn)群體的接觸機(jī)會(huì)，2015年后恙蟲病發(fā)病呈下降趨勢(shì)，APC為 -7.89% 。鹽城市2013—2023年恙蟲病報(bào)告發(fā)病率總體呈上升趨勢(shì)，但差異無(wú)統(tǒng)計(jì)學(xué)意義，與常州市、連云港市等地恙蟲病發(fā)病趨勢(shì)一致[-，這可能與2019年報(bào)告發(fā)病數(shù)大幅增加，作為離群值對(duì)趨勢(shì)擬合產(chǎn)生影響有關(guān)。

圖6鹽城市2013—2023年恙蟲病AnselinLocalMoran'sI分析結(jié)果

本研究集中度與圓形分布法結(jié)果顯示，恙蟲病發(fā)病具嚴(yán)格季節(jié)性，發(fā)病高峰期在10月13日至11月18日，屬秋季型疫區(qū)，符合長(zhǎng)江以北地區(qū)恙蟲病流行特點(diǎn)。10月一11月正值秋收農(nóng)忙時(shí)節(jié)，農(nóng)民于田地、野外等勞作過(guò)程中接觸恙蟲病傳染源及傳播媒介的機(jī)會(huì)更多；此外，還可能與宿主動(dòng)物黑線姬鼠及與媒介生物恙螨的數(shù)量、密度、活動(dòng)性的增加有關(guān)，與國(guó)內(nèi)其他地區(qū)研究類似[12-13]。

本研究結(jié)果顯示，鹽城市2013—2023年恙蟲病女性發(fā)病率高于男性，發(fā)病職業(yè)人群以農(nóng)民為主，發(fā)病年齡段以50～80歲為主。其原因在于農(nóng)村年輕男性多外出務(wù)工，女性及老年人留守務(wù)農(nóng)，增加了暴露機(jī)會(huì)，故呈現(xiàn)女性、農(nóng)民、老年高發(fā)的特征。

鹽城市各縣（市、區(qū)）均有發(fā)病，且縣區(qū)發(fā)病率高于大市區(qū)、沿海地帶發(fā)病率高于非沿海地帶。這與各縣區(qū)較大市區(qū)更多致力于農(nóng)業(yè)生產(chǎn)活動(dòng)有關(guān)，由于區(qū)域農(nóng)業(yè)種植規(guī)模化快速發(fā)展，土地利用率提高，人員流動(dòng)加快，利于恙蟲病流行。此前研究發(fā)現(xiàn)，鹽城市恙蟲病發(fā)病高峰前1至2個(gè)月會(huì)出現(xiàn)鼠密度升高現(xiàn)象，且鹽城市沿海灘涂地區(qū)以黑線姬鼠為絕對(duì)優(yōu)勢(shì)鼠種，為沿海地帶恙蟲病發(fā)病率高于非沿海地帶提供佐證。

全局空間自相關(guān)分析結(jié)果顯示，鹽城市2019—2022年恙蟲病存在空間聚集性，由于2019年響水縣報(bào)告發(fā)病數(shù)較多，使得恙蟲病明顯輻射周邊地區(qū)，造成發(fā)病率的升高以及全局的空間聚集。2013—2023年建湖縣以高-低聚集模式轉(zhuǎn)變?yōu)榈?低聚集模式，表明該地區(qū)防控措施取得成效，發(fā)病率水平趨近周邊地區(qū)。

在今后恙蟲病的防控中要重點(diǎn)關(guān)注秋收季節(jié)后疫情監(jiān)測(cè)，強(qiáng)化中老年農(nóng)民健康教育和個(gè)人防護(hù)等工作，不斷提高基層醫(yī)療機(jī)構(gòu)識(shí)別、發(fā)現(xiàn)恙蟲病患者的能力。

[參考文獻(xiàn)]

[1]李貴昌，栗冬梅，李燚，等.2006—2016年我國(guó)恙蟲病流行特征分析[J].疾病監(jiān)測(cè)，2018，33（2）：139-143.

[2]彭培英，徐蕾，王谷仙，等.1952—1989 年和2006—2017年中國(guó)大陸恙蟲病流行及時(shí)空分布特征[J]中國(guó)人獸共患病學(xué)報(bào)，2022，38（9）：818-823.

[3]岳玉娟，王玉姣，李貴昌，等.2006一2018年中國(guó)大陸恙蟲病高發(fā)區(qū)流行病學(xué)特征分析[J].疾病監(jiān)測(cè)，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年恙蟲病流行特征和時(shí)空聚集性分析[J].中國(guó)媒介生物學(xué)及控制雜志，2019，30（3）：237-243.

[7]陳安喜，孔衛(wèi)榮，田野，等.江蘇省常州市2012—2021年恙蟲病流行特征分析[J].中國(guó)媒介生物學(xué)及控制雜志，2023，34（3）：372-377.

[8]鄭艷澤，陳玉，張嘉陵.2010-2022年連云港市恙蟲病流行病學(xué)特征分析[J].公共衛(wèi)生與預(yù)防醫(yī)學(xué)2024，35（3）：95-98.

[9]鄺仕成，吳彪，陳貽彪，等.海口市某醫(yī)院330例恙蟲病流行和臨床特征分析［J].中國(guó)感染與化療雜志，2021，21（5）：512-516.

[10]張嘉溪，譚盛葵.恙蟲病流行病學(xué)研究新進(jìn)展[J].中國(guó)熱帶醫(yī)學(xué)，2022，22（3）：274-278.

[11]李孟磊，王文華，李懿，等.河南省 2013—2019 年恙蟲病流行病學(xué)特征分析[J].實(shí)用預(yù)防醫(yī)學(xué)雜志，2021，28（10）：1184-1186.

[12]薛莉穎，謝曉菲，王高玉，等.海南省680例不明原因發(fā)熱患者恙蟲病東方體感染調(diào)查[J].中國(guó)熱帶醫(yī)學(xué)，2022，22（10）：930-934.

[13]馬蒙蒙，王大虎，魯影，等.廣州市2015一2019 年恙蟲病流行特征分析[J].醫(yī)學(xué)動(dòng)物防制，2021，37（8）：744-747.

（本文編輯繆宏建）