2016—2020年青藏高原東部地區(qū)棘球蚴病空間聚集性分析

施丹丹，蕢嫣，韓帥，伍衛(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.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 結(jié) 果

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。

3 討論

空間流行病學(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)研究。