施丹丹,蕢嫣,韓帥,伍衛(wèi)平
中國(guó)疾病預(yù)防控制中心寄生蟲病預(yù)防控制所(國(guó)家熱帶病研究中心),國(guó)家衛(wèi)生健康委員會(huì)寄生蟲病原與媒介生物學(xué)重點(diǎn)實(shí)驗(yàn)室,世界衛(wèi)生組織熱帶病合作中心,國(guó)家級(jí)熱帶病國(guó)際聯(lián)合研究中心,上海 200025
棘球蚴?。ㄋ追Q包蟲病)是棘球絳蟲的幼蟲(棘球蚴)引起的人獸共患寄生蟲病,在全球范圍內(nèi)廣泛分布[1]。棘球蚴病嚴(yán)重威脅人群健康,是重要的公共衛(wèi)生問題之一,同時(shí)給畜牧業(yè)造成了巨大經(jīng)濟(jì)損失。我國(guó)是囊型和泡型棘球蚴病的流行區(qū),受威脅人口數(shù)和患病人數(shù)均為全球首位[2-3]。青藏高原東部地區(qū)位于青海省和四川省交界處,包括青海省玉樹州、果洛州和四川省甘孜州、阿壩州,各縣(市)均為棘球蚴病的流行區(qū),且流行程度較重。2012年調(diào)查數(shù)據(jù)顯示,青海省人群棘球蚴病患病率為1.10%,其中果洛州為5.20%,玉樹州為4.54%[4];四川省人群棘球蚴病患病率為1.08%,其中甘孜州為1.86%,阿壩州為0.79%[5]。近年來我國(guó)加大了棘球蚴病防治力度,并取得了明顯成效,但局部地區(qū)的患病率仍然較高,存在傳播風(fēng)險(xiǎn)[6]。
目前,空間統(tǒng)計(jì)分析方法在棘球蚴病流行病學(xué)研究中應(yīng)用廣泛,有助于評(píng)估患病率差別、識(shí)別疾病分布特征和聚集性、分析疾病影響因素[7-8]。了解和掌握棘球蚴病高流行地區(qū)的空間分布特征,對(duì)棘球蚴病防治策略的制定和資源的合理分配有重要參考意義。本研究通過分析2016—2020 年青藏高原東部地區(qū)人群棘球蚴病發(fā)病率空間分布特征,探索其空間聚集性,為棘球蚴病的防治提供參考。
1.1 研究區(qū)域 選取青藏高原東部地區(qū),包括四川省的甘孜州和阿壩州、青海省的玉樹州和果洛州共計(jì)43個(gè)縣(市)作為研究區(qū)域,見圖1。
1.2 數(shù)據(jù)來源 2016—2020年研究區(qū)域的棘球蚴病發(fā)病數(shù)據(jù)來源于2020 年全國(guó)傳染病監(jiān)測(cè)報(bào)告數(shù)據(jù)庫(kù)。監(jiān)測(cè)方法:按照各省棘球蚴病流行種類及流行程度的比例,計(jì)算抽取各類流行村的數(shù)量,采用隨機(jī)抽樣的方法抽取監(jiān)測(cè)點(diǎn);在監(jiān)測(cè)點(diǎn)選擇一個(gè)或相鄰的幾個(gè)村的全部居民作為監(jiān)測(cè)人群,每年進(jìn)行超聲檢查,并以血清學(xué)方法檢測(cè)抗棘球蚴抗體進(jìn)行輔助診斷[9]。2016—2018年青藏高原東部地區(qū)僅有28 個(gè)縣(市)開展了監(jiān)測(cè)工作,2019—2020 年43 個(gè)縣(市)均開展了監(jiān)測(cè)工作。
1.3 空間自相關(guān)分析 以縣(市)為單位,對(duì)2016—2020 年青藏高原東部地區(qū)43 個(gè)縣(市)棘球蚴病發(fā)病率進(jìn)行全局空間自相關(guān)分析、局部空間自相關(guān)分析和熱點(diǎn)分析,2016—2018年未進(jìn)行監(jiān)測(cè)工作的縣(市)數(shù)據(jù)以缺失值來進(jìn)行分析。其中空間自相關(guān)是指發(fā)病率在空間位置上的相關(guān)性,可反映發(fā)病率在研究縣(市)的聚集程度。若發(fā)病率在空間上呈現(xiàn)中心與周圍分布一致,則稱為空間正相關(guān);若呈現(xiàn)中心與周圍分布相反,則稱為空間負(fù)相關(guān);若發(fā)病率在空間分布上呈現(xiàn)隨機(jī)性,則表明空間相關(guān)性不顯著[10]。局部空間自相關(guān)分析選用Local Moran′sI統(tǒng)計(jì)量,利用LISA(local indications of spatial association)來進(jìn)行分析,變量在空間上存在聚集區(qū)域的類型,可分為四種:“高-高”聚集(H-H),表示某病發(fā)病率高于均值的地區(qū)被高于均值的領(lǐng)域包圍;“高-低”聚集(H-L),表示某病發(fā)病率高于均值的地區(qū)被低于均值的領(lǐng)域包圍;“低-低”聚集(L-L),表示某病發(fā)病率低于均值的地區(qū)被低于均值的領(lǐng)域包圍;“低-高”聚集(L-H),表示某病發(fā)病率低于均值的地區(qū)被高于均值的領(lǐng)域包圍[11]。采用空間局部Getis-Ord Gi*熱點(diǎn)分析發(fā)現(xiàn)變量在空間分布上存在熱點(diǎn)地區(qū)(高值聚集區(qū))和冷點(diǎn)地區(qū)(低值聚集區(qū))[12]。
1.4 統(tǒng)計(jì)學(xué)分析 運(yùn)用Excel 2007 軟件建立數(shù)據(jù)庫(kù)。應(yīng)用SPSS 19.0 對(duì)數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析,對(duì)棘球蚴病發(fā)病率進(jìn)行時(shí)間趨勢(shì)分析。利用ArcGIS 10.1軟件進(jìn)行全局空間自相關(guān)分析、局部空間自相關(guān)分析和熱點(diǎn)分析,并對(duì)結(jié)果進(jìn)行可視化展示。檢驗(yàn)水準(zhǔn)采用α=0.05。
2.1 棘球蚴病發(fā)病率 2016—2020年青藏高原東部地區(qū)人群棘球蚴病累計(jì)新發(fā)病例598 例,年平均發(fā)病率為0.032%,各年發(fā)病率為0.013%~0.553%。其中,2016 年發(fā)病率最高,為0.553%;2020 年最低,為0.013%,棘球蚴病發(fā)病率隨時(shí)間推移逐漸下降(χ2=802.991,P<0.05)。見表1。
表1 2016—2020年青藏高原東部地區(qū)棘球蚴病發(fā)病率
2.2 棘球蚴病發(fā)病率地區(qū)分布 新發(fā)病例發(fā)病率>0.1%的縣(市),2016年為青海省曲麻萊縣、久治縣、瑪沁縣、玉樹市、達(dá)日縣、稱多縣和四川省石渠縣、白玉縣、甘孜縣、德格縣、紅原縣、理塘縣;2017年為青海省達(dá)日縣、瑪沁縣、曲麻萊縣、雜多縣、稱多縣和四川省白玉縣;2018 年為青海省稱多縣、雜多縣和四川省德格縣、馬爾康市、白玉縣、道孚縣;2019年為青海省治多縣、雜多縣和四川省阿壩縣;2020年僅為青海省久治縣和四川省紅原縣。見圖2。
2.3 棘球蚴病發(fā)病率空間自相關(guān)分析
2.3.1 全局空間自相關(guān)分析 2016—2018 年全局Moran′sI值分別為0.005、0.088、-0.046(Z=0.281、1.196、-0.046,P均>0.05),表明2016—2018 年研究區(qū)域棘球蚴病發(fā)病率在空間分布上不具有空間自相關(guān)性,可能是隨機(jī)分布。2019、2020 年全局Moran′sI值為0.179、0.166(Z=2.273、2.316,P均<0.05),表明棘球蚴病發(fā)病率呈空間正相關(guān)分布,存在高值與高值相鄰或低值與低值相鄰的聚集現(xiàn)象。
2.3.2 局部空間自相關(guān)分析 2019—2020 年青藏高原東部地區(qū)棘球蚴病發(fā)病率局部空間自相關(guān)分析結(jié)果顯示,2019 年發(fā)病率呈現(xiàn)“高-低”和“低-高”聚集,四川省德格縣、道孚縣、瀘定縣呈“高-低”聚集,青海省玉樹市呈“低-高”聚集;2020年發(fā)病率呈現(xiàn)“高-高”“高-低”和“低-高”聚集,四川省阿壩縣、紅原縣、馬爾康市和青海省治多縣呈“高-高”聚集,四川省白玉縣呈“高-低”聚集,四川省若爾蓋縣、松潘縣、理縣呈“低-高”聚集。見圖3。
2.4 棘球蚴病發(fā)病率熱點(diǎn)分析 空間局部Getis-Ord Gi*熱點(diǎn)分析結(jié)果顯示,2019年熱點(diǎn)99%置信區(qū)間地區(qū)包含青海省治多縣和雜多縣,95%置信區(qū)間地區(qū)包含青海省玉樹市;2020年熱點(diǎn)99%置信區(qū)間地區(qū)包含四川省阿壩縣、紅原縣、理縣,95%置信區(qū)間地區(qū)包含四川省若爾蓋縣、松潘縣、黑水縣、馬爾康市。見圖4。
空間流行病學(xué)旨在運(yùn)用地理信息系統(tǒng)等空間信息技術(shù)研究疾病分布特征和關(guān)鍵影響因素,為疾病預(yù)防控制工作提供策略支持[13]。隨著現(xiàn)代空間信息技術(shù)的快速發(fā)展,空間流行病學(xué)正越來越多地被應(yīng)用于各種疾病研究中,在公共衛(wèi)生領(lǐng)域發(fā)揮著重要作用[14-15]。
2006 年以來,全國(guó)逐步開始實(shí)施地方棘球蚴病防治項(xiàng)目,通過采取“以控制傳染源為主,中間宿主防控與病人查治相結(jié)合”的綜合防治策略,棘球蚴病的防治工作取得了一定成效[6]。本研究發(fā)現(xiàn),在青藏高原東部地區(qū)青海省人群棘球蚴病發(fā)病率略高于四川省,這與之前研究結(jié)果一致[2]。調(diào)查的青藏高原東部地區(qū)43 個(gè)縣(市)2016—2020 年發(fā)病率呈下降趨勢(shì),2020 年發(fā)病率已低至0.013%,說明近年來棘球蚴病防治工作取得了顯著成效。棘球蚴病流行趨勢(shì)趨于平穩(wěn),但部分縣偶有回升,出現(xiàn)個(gè)別新發(fā)病例,提示仍需加強(qiáng)對(duì)棘球蚴病的動(dòng)態(tài)監(jiān)測(cè)。同時(shí),2016—2020 年棘球蚴病發(fā)病率>0.1%的地區(qū)也逐年減少,2016—2018年青藏高原東部地區(qū)棘球蚴病發(fā)病率較高的地區(qū)主要分布在中部交界地區(qū),包括青海省曲麻萊縣、稱多縣、瑪多縣和達(dá)日縣,四川省石渠縣、德格縣、白玉縣和阿壩縣;2019年發(fā)病率較高地區(qū)為青海省治多縣和雜多縣,四川省阿壩縣;2020年發(fā)病率較高地區(qū)為青海省久治縣和四川省紅原縣??赡苡捎?016年開始實(shí)施《全國(guó)包蟲病監(jiān)測(cè)方案》,部分縣剛開展人群監(jiān)測(cè)工作,新發(fā)現(xiàn)的病例中可能包含了既往發(fā)病但未及時(shí)檢查出來的病例。2017年青海省達(dá)日縣新發(fā)現(xiàn)121例病例,可能由于當(dāng)時(shí)基層專業(yè)人員能力有限,該縣2016年未完成全部監(jiān)測(cè)工作,部分監(jiān)測(cè)任務(wù)2017年才得以完成。
大量研究顯示,棘球蚴病的分布具有一定的空間聚集性。Schmidberger等[16]對(duì)1992—2016年德國(guó)棘球蚴病數(shù)據(jù)進(jìn)行空間自相關(guān)性分析,結(jié)果顯示整個(gè)德國(guó)的病例分布具有聚集性,局部空間自相關(guān)分析結(jié)果表明巴登-符騰堡州東南部地區(qū)和巴戈利亞州西南部地區(qū)為高危區(qū)。He 等[17]對(duì)四川省2016—2019年人群棘球蚴病流行情況進(jìn)行調(diào)查,全局空間自相關(guān)分析顯示存在空間正相關(guān),局部空間自相關(guān)分析發(fā)現(xiàn)“高-高”聚集區(qū)在四川省西北部,“低-低”聚集區(qū)在四川省東南部。趙月等[18]對(duì)四川省甘孜州2014年棘球蚴病的空間分析結(jié)果顯示棘球蚴病空間分布存在聚集性。本研究發(fā)現(xiàn)2016—2018 年青藏高原東部地區(qū)棘球蚴病發(fā)病率不具有空間聚集趨勢(shì),2019—2020年存在空間正相關(guān),具有明顯的空間聚集特征,存在“高-高”“高-低”和“低-高”聚集區(qū),可能由于2016—2018年研究區(qū)域的篩查地區(qū)不完整,有15個(gè)縣(市)的數(shù)據(jù)缺失。局部Getis′G熱點(diǎn)分析結(jié)果與局部空間自相關(guān)分析的“高-低”和“高-高”聚集地區(qū)較為一致,2019年熱點(diǎn)區(qū)域位于青海省治多縣、雜多縣、玉樹市,2020年熱點(diǎn)區(qū)域位于四川省阿壩縣、紅原縣、理縣、若爾蓋縣、松潘縣、黑水縣、馬爾康市??赡苡捎谶@些地區(qū)位于青南高原,寒冷潮濕的氣候環(huán)境有利于棘球絳蟲蟲卵存活,且多房棘球絳蟲的中間宿主嚙齒類動(dòng)物數(shù)量多、分布廣。此外這些地區(qū)以畜牧業(yè)為主,自然環(huán)境和牲畜類型都有利于棘球絳蟲完成生活史循環(huán),造成棘球蚴病的流行。甘孜州作為棘球蚴病的重點(diǎn)防控區(qū)域,發(fā)病率逐漸降低,而阿壩州近年來出現(xiàn)個(gè)別新發(fā)病例,這可能是熱點(diǎn)區(qū)域轉(zhuǎn)移的主要原因。
綜上,青藏高原東部地區(qū)棘球蚴病防治工作取得了一定成效,但仍存在空間分布聚集性。對(duì)于“高-低”和“高-高”聚集區(qū)域和熱點(diǎn)區(qū)域,應(yīng)強(qiáng)化棘球蚴病的防治工作,加強(qiáng)犬只管理,開展健康宣傳,同時(shí)在既往人群篩查的基礎(chǔ)上,加強(qiáng)重點(diǎn)區(qū)域的全人群普查。此外,本研究以縣級(jí)為單位,研究結(jié)果僅可反映縣級(jí)聚集情況,后續(xù)可從更小的尺度進(jìn)行研究;另外,此次研究不涉及影響因素,今后可開展相關(guān)研究。