[中圖分類號]R513.2 [文獻標志碼] B [DOI] 10.19767/j.cnki.32-1412.2025.05.025
恙蟲病又名叢林斑疹傷寒,是由恙蟲病立克次體引起的急性傳染病,系一種自然疫源性疾病,嚙齒類為主要傳染源,恙螨幼蟲為傳播媒介。該病起病急,臨床主要表現為高熱、毒血癥、皮疹、焦痂和淋巴結腫大等,可導致多臟器損害[2-3]。鹽城市所轄東臺縣1986年首次在長江以北地區報告恙蟲病,2006年所轄濱??h發生恙蟲病暴發疫情,之后將恙蟲病納入監測報告傳染病,在2019年報告病例數最多,采取一系列綜合防控措施后,疫情逐漸下降。為更好了解鹽城市近幾年恙蟲病流行趨勢,為科學防控提供決策依據,現對鹽城市2013—2023年恙蟲病流行特征進行分析。
1資料與方法
1.1資料來源數據來源于\"中國疾病預防控制信息系統”傳染病監測信息系統,統計時間為2013年1月1日—2023年12月31日,按發病日期和現住址統計,內容包括地區、年齡、性別、發病時間等傳染病報卡信息及發病數、發病率。
1.2分析方法運用Excel2016軟件錄入、整理數據,SAS9.4統計學軟件進行數據分析。三間分布采用描述性流行病學方法進行分析;計數資料以 n(%) 表示,組間比較采用 χ2 檢驗;應用 Joinpoint Regres-sionProgram5.3.0.0軟件擬合回歸模型,計算年度變化百分比(annual percent change,APC)、平均年度變化百分比(average annual percent change,AAPC)并分析變化趨勢;采用集中度法和圓形分布法描述季節特征,使用Cochran-Armitage檢驗有序分類變量的線性趨勢。 Plt;0.05 為差異有統計學意義。
利用ArcGis10.8軟件進行空間自相關分析。采用全局莫蘭指數(Moran'sI)衡量鹽城市發病率的整體分布情況,當Moran's Igt;0 且 Z 值 ?1.96 時為正相關,表明整體呈現聚集性分布;當Moran's Ilt;0 且 Z 值 ≤-1.96 時為負相關,表示空間上呈現離散性分布;當Moran'sI接近0且 Z 值介于-1.96~1.96,說明整體呈隨機性分布。Moran'sI絕對值越大,表明空間相關性越大。全局Moran'sI可能會忽視鹽城市局部聚集性,因此本研究還采用安瑟倫局部莫蘭指數(AnselinLocalMoran'sI)分析局部地區的聚集現象。聚集性可表現為“高-高”“高-低”、“低-高”“低-低\"4種聚集模式,高-高聚集表明高值區域被高值區域包圍。
2結果
2.1發病概況鹽城市2013—2023年累計報告恙蟲病5326例,無死亡病例報告,年平均報告發病率為6.84/10萬,2019年報告發病率最高為11.63/10萬,2013年報告發病率最低為3.50/10萬,各年份報告發病率差異有統計學意義( χ2=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萬),差異具有統計學意義 (χ2=1003.184,Plt;0.001) ;沿海地帶(響水縣、濱??h、射陽縣、大豐區、東臺市)年平均報告發病率(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萬),差異有統計學意義( ?x2=65.161 ,Plt;0.0001 )。鹽城市2013—2023年各年齡組均有病例報告,各組報告發病率差異有統計學意義( χ2= 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年Moran'sI最大,表明該年鹽城市恙蟲病的發病率聚集程度最強。見表1。AnselinLocalMoran'sI顯示,建湖縣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年恙蟲病AnselinLocalMoran'sI分析結果

本研究集中度與圓形分布法結果顯示,恙蟲病發病具嚴格季節性,發病高峰期在10月13日至11月18日,屬秋季型疫區,符合長江以北地區恙蟲病流行特點。10月一11月正值秋收農忙時節,農民于田地、野外等勞作過程中接觸恙蟲病傳染源及傳播媒介的機會更多;此外,還可能與宿主動物黑線姬鼠及與媒介生物恙螨的數量、密度、活動性的增加有關,與國內其他地區研究類似[12-13]。
本研究結果顯示,鹽城市2013—2023年恙蟲病女性發病率高于男性,發病職業人群以農民為主,發病年齡段以50~80歲為主。其原因在于農村年輕男性多外出務工,女性及老年人留守務農,增加了暴露機會,故呈現女性、農民、老年高發的特征。
鹽城市各縣(市、區)均有發病,且縣區發病率高于大市區、沿海地帶發病率高于非沿海地帶。這與各縣區較大市區更多致力于農業生產活動有關,由于區域農業種植規?;焖侔l展,土地利用率提高,人員流動加快,利于恙蟲病流行。此前研究發現,鹽城市恙蟲病發病高峰前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.
(本文編輯繆宏建)