2018—2022年內(nèi)蒙古自治區(qū)人間布魯氏菌病時(shí)空聚集性分析

宋 健, 韓 松, 王虎虎, 高永明, 初 迪, 尉瑞平

布魯氏菌病(bruollosis,簡(jiǎn)稱(chēng)布病)又稱(chēng)地中海熱、馬爾他熱等,是由布氏菌通過(guò)多種途徑侵入機(jī)體而引起的一種人獸共患變態(tài)反應(yīng)性傳染病。臨床表現(xiàn)以發(fā)熱、乏力、多汗、關(guān)節(jié)疼痛、肝脾及淋巴結(jié)腫大為主要特征,常造成多系統(tǒng)多臟器損傷[1]。布病在全世界廣泛流行,主要在地中海地區(qū)、非洲、亞洲等低收入和中等收入國(guó)家[2],我國(guó)布病疫情也較為嚴(yán)重。內(nèi)蒙古自治區(qū)為畜牧業(yè)發(fā)達(dá)省份,是人間布病的主要流行區(qū),全區(qū)布病疫情范圍不斷擴(kuò)大且檢出布魯氏菌為侵襲力和致病力最強(qiáng)的羊生物種[1,3-4],提示全區(qū)布病疫情處于嚴(yán)峻態(tài)勢(shì)。目前,全區(qū)關(guān)于布病的研究主要以監(jiān)測(cè)分析和血清學(xué)監(jiān)測(cè)為主,未見(jiàn)數(shù)據(jù)的空間屬性。本文通過(guò)利用空間自相關(guān)和回顧性時(shí)空掃描分析方法,探討人間布病空間聚集性,為布病的預(yù)防和控制工作提供參考。

1 資料與方法

1.1 數(shù)據(jù)來(lái)源 收集2018—2022年“中國(guó)傳染病報(bào)告信息管理系統(tǒng)”中內(nèi)蒙古地區(qū)報(bào)告布病臨床診斷和實(shí)驗(yàn)室確診病例信息(年齡、性別、職業(yè)、發(fā)病時(shí)間等),人口資料來(lái)源于內(nèi)蒙古自治區(qū)統(tǒng)計(jì)年鑒。

1.2 研究方法

1.2.1 空間自相關(guān)分析 空間自相關(guān)分析已經(jīng)廣泛用于傳染病的時(shí)空特征研究,本文利用相鄰權(quán)重法來(lái)構(gòu)建空間權(quán)重矩陣,以縣為單位通過(guò)Moran′sI和局部空間自相關(guān)系數(shù)(LISA)進(jìn)行全局空間自相關(guān)分析和局域空間自相關(guān)分析。本文利用Moran′sI進(jìn)行全局空間自相關(guān)和局域空間自相關(guān)分析。全局空間自相關(guān)用于整個(gè)區(qū)域內(nèi)的空間分布情況,判斷是否有空間聚集性,但無(wú)法在地圖上顯示出聚集情況具體發(fā)生的區(qū)域。全局Moran′sI統(tǒng)計(jì)量取值范圍在[-1,1],I小于0時(shí)代表空間負(fù)相關(guān),越接近-1代表空間單元分布差異越大;I=0時(shí)代表空間區(qū)域內(nèi)發(fā)病無(wú)相關(guān);大于0時(shí)代表空間正相關(guān),越接近1代表空間單元關(guān)系越密切。

局域空間自相關(guān)是評(píng)估每個(gè)空間單元與鄰近單元的相關(guān)程度。局域自相關(guān)的聚集模式可結(jié)合Moran散點(diǎn)圖和局域Moran′sI檢驗(yàn)進(jìn)行判斷,并以L(fǎng)ISA聚集性地圖進(jìn)行直觀展示?？臻g關(guān)聯(lián)模式分為4種:高-高聚集、低-低聚集、高-低聚集、低-高聚集。其中高-高聚集為發(fā)病的“熱點(diǎn)地區(qū)”、低-低聚集為發(fā)病的“冷點(diǎn)”區(qū)域[5-6],2種情況具有較強(qiáng)的空間正相關(guān)性;高-低模式和低-高模式表示其空間屬性與相鄰地區(qū)的空間屬性不同,呈現(xiàn)空間負(fù)相關(guān)性。

1.2.2 時(shí)空掃描聚類(lèi)分析 時(shí)空掃描聚類(lèi)分析可探測(cè)疾病在時(shí)空范圍內(nèi)發(fā)生數(shù)量的異常改變,檢驗(yàn)改變是否由于隨機(jī)變異造成的,是通過(guò)動(dòng)態(tài)變化的圓柱形掃描窗口對(duì)不同的時(shí)間和地理區(qū)域進(jìn)行掃描,探測(cè)所有病例可能發(fā)生聚集的時(shí)間和地理位置,圓柱的高代表掃描時(shí)間的范圍,圓柱的底部代表掃描范圍,通過(guò)計(jì)算動(dòng)態(tài)變化的窗口下的對(duì)數(shù)似然比(LLR)來(lái)判斷可能的聚集區(qū),最終選擇LLR值最大的窗口為1類(lèi)聚集區(qū),其他LLR經(jīng)過(guò)檢驗(yàn)具有統(tǒng)計(jì)學(xué)意義的窗口定義為2類(lèi)及以上聚集區(qū),LLR值越大的窗口,表明越有可能是聚集區(qū)域。本文采用Satscan9.5對(duì)2018—2022年全區(qū)布病發(fā)病數(shù)據(jù)進(jìn)行時(shí)空掃描分析,將模型最大空間掃描區(qū)域設(shè)定為全區(qū)總?cè)丝跀?shù)的30%,掃描時(shí)間以月為單位,縣(旗、區(qū))為空間掃描單位,時(shí)空分析選擇Poisson分布模型進(jìn)行。各年發(fā)病率差異分析采用χ2檢驗(yàn),檢驗(yàn)水準(zhǔn)α=0.05。

2 結(jié) 果

圖1 2018—2022年內(nèi)蒙古自治區(qū)布病各年發(fā)病情況Fig.1 Incidence of brucellosis in the Inner Mongolia Autonomous Region from 2018 to 2022

2.2 全局空間自相關(guān)分析 2018—2022年各年全局空間自相關(guān)系數(shù)Moran′sI在0.231～0.285,經(jīng)蒙特卡羅檢驗(yàn)后,2018—2022年全區(qū)布病發(fā)病率全局空間自相關(guān)系數(shù)Moran′sI均有統(tǒng)計(jì)學(xué)差異,說(shuō)明全區(qū)布病病例呈聚集性分布,提示全區(qū)布病發(fā)病率呈現(xiàn)明顯的空間正自相關(guān),存在一定的空間聚集區(qū)(表1)。

表1 2018-2022年內(nèi)蒙古自治區(qū)布病發(fā)病率全局空間自相關(guān)分析Tab.1 Spatial autocorrelation analysis of the incidence rate of brucellosis in the Inner Mongolia Autonomous Region from 2018 to 2022

2.3 局域空間自相關(guān)分析 對(duì)各年發(fā)病率進(jìn)行局域空間自相關(guān)分析顯示。2018年“高-高”聚集區(qū)共有10個(gè),主要集中在通遼市和錫林郭勒盟。2019年“高-高”聚集區(qū)共有8個(gè),主要分布在通遼市和巴彥淖爾市。2020年“高-高”聚集區(qū)共有10個(gè),主要集中在通遼市、巴彥淖爾市和包頭市。2021年“高-高”聚集區(qū)共有6個(gè),主要分布在巴彥淖爾市和包頭市。2022年“高-高”聚集區(qū)共有6個(gè),主要集中在巴彥淖爾市、包頭市和鄂爾多斯市。近5年,“高-高聚集區(qū)”存在由蒙東地區(qū)向蒙西地區(qū)轉(zhuǎn)移、蒙西“高-高”聚集區(qū)范圍擴(kuò)大的現(xiàn)象(圖2)。

(審圖號(hào):(2020)4619號(hào))圖2 2018-2022年內(nèi)蒙古自治區(qū)布病發(fā)病率聚集性地圖Fig.2 Cluster map of the brucellosis incidence rate in the Inner Mongolia Autonomous Region from 2018 to 2022

2.4 時(shí)空掃描聚類(lèi)分析 對(duì)2018—2022年全區(qū)布病報(bào)告發(fā)病進(jìn)行時(shí)空掃描分析,結(jié)果顯示全區(qū)布病發(fā)病存在時(shí)空聚集性分布,共掃描到7個(gè)聚集區(qū),其中1類(lèi)聚集區(qū)(LLR 值最大,為 3 954.31)主要分布以鄂爾多斯市鄂托克旗為中心,覆蓋半徑為103.01 km,涉及3個(gè)地區(qū),分別為鄂爾多斯市烏審旗、杭錦旗以及烏海市的海南區(qū),聚集時(shí)間為2020年6月至2022年11月,該聚集區(qū)發(fā)病風(fēng)險(xiǎn)為周?chē)貐^(qū)的5.83倍。2類(lèi)聚集區(qū)以呼倫貝爾市根河市為中心(LLR為2 772.54),覆蓋半徑為156.49 km,涉及2個(gè)地區(qū),分別為呼倫貝爾市鄂倫春自治旗和額爾古納市,聚集時(shí)間為 2020年6月至2022年11月,該聚集區(qū)發(fā)病風(fēng)險(xiǎn)為周?chē)貐^(qū)的5.81倍。3類(lèi)聚集區(qū)以錫林郭勒盟阿巴嘎旗為中心(LLR為2 190.98),覆蓋半徑為234.55 km,涉及旗縣最多為9個(gè)地區(qū),分別為錫林郭勒盟的錫林浩特市、蘇尼特左旗、蘇尼特右旗、正藍(lán)旗、正鑲白旗、鑲黃旗、東烏珠穆沁旗、西烏珠穆沁旗和赤峰市的克什克騰旗,聚集時(shí)間為 2019年5月至2021年10月,該聚集區(qū)發(fā)病風(fēng)險(xiǎn)為周?chē)貐^(qū)的3.22倍。見(jiàn)表2和圖3。

表2 2018-2022年內(nèi)蒙古自治區(qū)布病時(shí)空聚集性分析Tab.2 Spatiotemporal clustering analysis of brucellosis in the Inner Mongolia Autonomous Region from 2018 to 2022

(審圖號(hào):(2020)4619號(hào))圖3 2018-2022年內(nèi)蒙古自治區(qū)布病時(shí)空聚集性分析示意圖Fig.3 Spatiotemporal clustering analysis of brucellosis in the Inner Mongolia Autonomous Region from 2018 to 2022

3 討 論

內(nèi)蒙古自治區(qū)是人間布病發(fā)病最嚴(yán)重的省份之一,近年因畜牧養(yǎng)殖業(yè)發(fā)展迅速,牲畜交易頻繁,人間布病疫情大幅度回升,發(fā)病數(shù)居高不下,近5年布病年發(fā)病率從38.36/10萬(wàn)上升到88.56/10萬(wàn),疫情呈現(xiàn)高發(fā)態(tài)勢(shì)。布病流行特征與其他省份相同[7-10],發(fā)病數(shù)呈現(xiàn)“三高”現(xiàn)象,即男性高于女性、農(nóng)民高于牧民以及青壯年高發(fā),與該群體為從事畜牧業(yè)的主要?jiǎng)趧?dòng)力、高頻次接觸傳染源和個(gè)人防護(hù)意識(shí)薄弱有關(guān),應(yīng)作為健康教育和疫情防控的重點(diǎn)人群。

傳染病監(jiān)測(cè)數(shù)據(jù)中既包含時(shí)間屬性又有空間屬性,通過(guò)空間聚集性分析方法對(duì)布病開(kāi)展時(shí)空研究來(lái)探索在時(shí)間和空間的聚集情況[11]。通過(guò)自相關(guān)分析和時(shí)空掃描聚類(lèi)分析相結(jié)合,更加系統(tǒng)全面的分析內(nèi)蒙古布病的時(shí)空分布特征。分析發(fā)現(xiàn),歷年布病發(fā)病率在空間分布上均呈聚集性,存在空間正相關(guān),2018年“高-高”聚集區(qū)主要位于內(nèi)蒙古蒙東和蒙中地區(qū)的通遼市和錫林郭勒盟周邊區(qū)域,2022年“高-高”聚集區(qū)主要位于內(nèi)蒙古蒙西地區(qū)的鄂爾多斯市和巴彥淖爾市周邊區(qū)域,五年布病的“熱點(diǎn)”區(qū)域由蒙東地區(qū)轉(zhuǎn)移到蒙西地區(qū),“熱點(diǎn)”地區(qū)的數(shù)量由10個(gè)降到6個(gè),提示聚集區(qū)域呈現(xiàn)縮小的趨勢(shì)。通遼市退出“熱點(diǎn)”區(qū)域,從“高-高”聚集分布轉(zhuǎn)變?yōu)殡S機(jī)分布,但其報(bào)告發(fā)病率依然處于高水平。2017年通遼市出臺(tái)《關(guān)于進(jìn)一步加強(qiáng)動(dòng)物防疫工作的通知》《通遼市職業(yè)病防治規(guī)劃(2017—2020年)》等相關(guān)文件用于控制人間和畜間布病疫情,分析結(jié)果說(shuō)明該地區(qū)近年針對(duì)布病的防控措施起到一定作用,但巴彥淖爾市和鄂爾多斯市布病疫情發(fā)展較快,在今后的工作中將其列入重點(diǎn)防控區(qū)域。

由于空間自相關(guān)分析未考慮時(shí)間因素,本研究通過(guò)時(shí)空掃描分析方法進(jìn)一步探索布病時(shí)空聚集性。結(jié)果顯示,各年均存在聚集時(shí)間,聚集時(shí)間1類(lèi)聚集區(qū)集中在內(nèi)蒙古蒙西地區(qū)(鄂爾多斯市和巴彥淖爾市部分區(qū)域),聚集時(shí)間主要集中在2020—2022年,結(jié)果與局域空間自相關(guān)分析結(jié)果相似。鄂爾多斯市鄂托克旗為1類(lèi)聚集中心,羊絨產(chǎn)業(yè)是鄂托克旗最具特色的傳統(tǒng)產(chǎn)業(yè),隨著鄂托克旗羊肉已正式成為首批中歐地理標(biāo)志協(xié)定的地理標(biāo)志產(chǎn)品,養(yǎng)殖規(guī)模呈現(xiàn)擴(kuò)大趨勢(shì),這可能導(dǎo)致該區(qū)域存在聚集性發(fā)病。進(jìn)一步分析結(jié)果顯示,除第6類(lèi)聚集區(qū)外,其余聚集區(qū)高發(fā)聚集時(shí)間均持續(xù)1年以上,提示這些地區(qū)可能存在疾病長(zhǎng)期的流行趨勢(shì)。

綜上所述,內(nèi)蒙古地區(qū)布病發(fā)病的重點(diǎn)人群以青壯年農(nóng)牧民為主,主要發(fā)病時(shí)間集中在3-8月,因此在每年發(fā)病高峰前,加強(qiáng)對(duì)高危人群的主動(dòng)干預(yù),通過(guò)發(fā)放防護(hù)用品、提高健康教育頻次等方式,降低人群感染風(fēng)險(xiǎn)。由于布病是一種人獸共患病,應(yīng)將“One Health”理念融合到布病的防治過(guò)程中,加強(qiáng)跨地區(qū)、跨部門(mén)、跨學(xué)科的交流合作[12-13],如畜牧部門(mén)將規(guī)?；B(yǎng)殖和散養(yǎng)戶(hù)納入監(jiān)管,加大對(duì)牲畜的檢疫,提高畜間免疫接種;市場(chǎng)監(jiān)管部門(mén)應(yīng)采取定期和突擊相結(jié)合的方式對(duì)市場(chǎng)開(kāi)展檢查,對(duì)銷(xiāo)售使用未檢疫牲畜產(chǎn)品的行為依法嚴(yán)處,消除食品安全隱患;交通運(yùn)輸部門(mén)、公安部門(mén)要積極配合畜牧部門(mén)加大對(duì)未附有動(dòng)物檢疫證明調(diào)運(yùn)牲畜等違法違規(guī)行為的查處力度,堅(jiān)決阻斷違法販運(yùn)鏈條。通過(guò)以上方式建立健康-畜牧-食品衛(wèi)生-交通運(yùn)輸?shù)牟疾》乐尉C合體,以有效遏制全區(qū)布病的傳播和蔓延。

利益沖突：無(wú)