施丹丹，蕢嫣，韓帥，伍衛(wèi)平

中國疾病預防控制中心寄生蟲病預防控制所（國家熱帶病研究中心），國家衛(wèi)生健康委員會寄生蟲病原與媒介生物學重點實驗室，世界衛(wèi)生組織熱帶病合作中心，國家級熱帶病國際聯(lián)合研究中心，上海 200025

棘球蚴?。ㄋ追Q包蟲?。┦羌蚪{蟲的幼蟲（棘球蚴）引起的人獸共患寄生蟲病，在全球范圍內廣泛分布[1]。棘球蚴病嚴重威脅人群健康，是重要的公共衛(wèi)生問題之一，同時給畜牧業(yè)造成了巨大經(jīng)濟損失。我國是囊型和泡型棘球蚴病的流行區(qū)，受威脅人口數(shù)和患病人數(shù)均為全球首位[2-3]。青藏高原東部地區(qū)位于青海省和四川省交界處，包括青海省玉樹州、果洛州和四川省甘孜州、阿壩州，各縣（市）均為棘球蚴病的流行區(qū)，且流行程度較重。2012年調查數(shù)據(jù)顯示，青海省人群棘球蚴病患病率為1.10%，其中果洛州為5.20%，玉樹州為4.54%[4]；四川省人群棘球蚴病患病率為1.08%，其中甘孜州為1.86%，阿壩州為0.79%[5]。近年來我國加大了棘球蚴病防治力度，并取得了明顯成效，但局部地區(qū)的患病率仍然較高，存在傳播風險[6]。

目前，空間統(tǒng)計分析方法在棘球蚴病流行病學研究中應用廣泛，有助于評估患病率差別、識別疾病分布特征和聚集性、分析疾病影響因素[7-8]。了解和掌握棘球蚴病高流行地區(qū)的空間分布特征，對棘球蚴病防治策略的制定和資源的合理分配有重要參考意義。本研究通過分析2016—2020 年青藏高原東部地區(qū)人群棘球蚴病發(fā)病率空間分布特征，探索其空間聚集性，為棘球蚴病的防治提供參考。

1 資料與方法

1.1 研究區(qū)域 選取青藏高原東部地區(qū)，包括四川省的甘孜州和阿壩州、青海省的玉樹州和果洛州共計43個縣（市）作為研究區(qū)域，見圖1。

1.2 數(shù)據(jù)來源 2016—2020年研究區(qū)域的棘球蚴病發(fā)病數(shù)據(jù)來源于2020 年全國傳染病監(jiān)測報告數(shù)據(jù)庫。監(jiān)測方法：按照各省棘球蚴病流行種類及流行程度的比例，計算抽取各類流行村的數(shù)量，采用隨機抽樣的方法抽取監(jiān)測點；在監(jiān)測點選擇一個或相鄰的幾個村的全部居民作為監(jiān)測人群，每年進行超聲檢查，并以血清學方法檢測抗棘球蚴抗體進行輔助診斷[9]。2016—2018年青藏高原東部地區(qū)僅有28 個縣（市）開展了監(jiān)測工作，2019—2020 年43 個縣（市）均開展了監(jiān)測工作。

1.3 空間自相關分析 以縣（市）為單位，對2016—2020 年青藏高原東部地區(qū)43 個縣（市）棘球蚴病發(fā)病率進行全局空間自相關分析、局部空間自相關分析和熱點分析，2016—2018年未進行監(jiān)測工作的縣（市）數(shù)據(jù)以缺失值來進行分析。其中空間自相關是指發(fā)病率在空間位置上的相關性，可反映發(fā)病率在研究縣（市）的聚集程度。若發(fā)病率在空間上呈現(xiàn)中心與周圍分布一致，則稱為空間正相關；若呈現(xiàn)中心與周圍分布相反，則稱為空間負相關；若發(fā)病率在空間分布上呈現(xiàn)隨機性，則表明空間相關性不顯著[10]。局部空間自相關分析選用Local Moran′sI統(tǒng)計量，利用LISA（local indications of spatial association）來進行分析，變量在空間上存在聚集區(qū)域的類型，可分為四種：“高-高”聚集（H-H），表示某病發(fā)病率高于均值的地區(qū)被高于均值的領域包圍；“高-低”聚集（H-L），表示某病發(fā)病率高于均值的地區(qū)被低于均值的領域包圍；“低-低”聚集（L-L），表示某病發(fā)病率低于均值的地區(qū)被低于均值的領域包圍；“低-高”聚集（L-H），表示某病發(fā)病率低于均值的地區(qū)被高于均值的領域包圍[11]。采用空間局部Getis-Ord Gi*熱點分析發(fā)現(xiàn)變量在空間分布上存在熱點地區(qū)（高值聚集區(qū)）和冷點地區(qū)（低值聚集區(qū)）[12]。

1.4 統(tǒng)計學分析 運用Excel 2007 軟件建立數(shù)據(jù)庫。應用SPSS 19.0 對數(shù)據(jù)進行統(tǒng)計分析，對棘球蚴病發(fā)病率進行時間趨勢分析。利用ArcGIS 10.1軟件進行全局空間自相關分析、局部空間自相關分析和熱點分析，并對結果進行可視化展示。檢驗水準采用α=0.05。

2 結 果

2.1 棘球蚴病發(fā)病率 2016—2020年青藏高原東部地區(qū)人群棘球蚴病累計新發(fā)病例598 例，年平均發(fā)病率為0.032%，各年發(fā)病率為0.013%～0.553%。其中，2016 年發(fā)病率最高，為0.553%；2020 年最低，為0.013%，棘球蚴病發(fā)病率隨時間推移逐漸下降（χ2=802.991，P＜0.05）。見表1。

表1 2016—2020年青藏高原東部地區(qū)棘球蚴病發(fā)病率

2.2 棘球蚴病發(fā)病率地區(qū)分布 新發(fā)病例發(fā)病率＞0.1%的縣（市），2016年為青海省曲麻萊縣、久治縣、瑪沁縣、玉樹市、達日縣、稱多縣和四川省石渠縣、白玉縣、甘孜縣、德格縣、紅原縣、理塘縣；2017年為青海省達日縣、瑪沁縣、曲麻萊縣、雜多縣、稱多縣和四川省白玉縣；2018 年為青海省稱多縣、雜多縣和四川省德格縣、馬爾康市、白玉縣、道孚縣；2019年為青海省治多縣、雜多縣和四川省阿壩縣；2020年僅為青海省久治縣和四川省紅原縣。見圖2。

2.3 棘球蚴病發(fā)病率空間自相關分析

2.3.1 全局空間自相關分析 2016—2018 年全局Moran′sI值分別為0.005、0.088、-0.046（Z=0.281、1.196、-0.046，P均＞0.05），表明2016—2018 年研究區(qū)域棘球蚴病發(fā)病率在空間分布上不具有空間自相關性，可能是隨機分布。2019、2020 年全局Moran′sI值為0.179、0.166（Z=2.273、2.316，P均＜0.05），表明棘球蚴病發(fā)病率呈空間正相關分布，存在高值與高值相鄰或低值與低值相鄰的聚集現(xiàn)象。

2.3.2 局部空間自相關分析 2019—2020 年青藏高原東部地區(qū)棘球蚴病發(fā)病率局部空間自相關分析結果顯示，2019 年發(fā)病率呈現(xiàn)“高-低”和“低-高”聚集，四川省德格縣、道孚縣、瀘定縣呈“高-低”聚集，青海省玉樹市呈“低-高”聚集；2020年發(fā)病率呈現(xiàn)“高-高”“高-低”和“低-高”聚集，四川省阿壩縣、紅原縣、馬爾康市和青海省治多縣呈“高-高”聚集，四川省白玉縣呈“高-低”聚集，四川省若爾蓋縣、松潘縣、理縣呈“低-高”聚集。見圖3。

2.4 棘球蚴病發(fā)病率熱點分析 空間局部Getis-Ord Gi*熱點分析結果顯示，2019年熱點99%置信區(qū)間地區(qū)包含青海省治多縣和雜多縣，95%置信區(qū)間地區(qū)包含青海省玉樹市；2020年熱點99%置信區(qū)間地區(qū)包含四川省阿壩縣、紅原縣、理縣，95%置信區(qū)間地區(qū)包含四川省若爾蓋縣、松潘縣、黑水縣、馬爾康市。見圖4。

3 討論

空間流行病學旨在運用地理信息系統(tǒng)等空間信息技術研究疾病分布特征和關鍵影響因素，為疾病預防控制工作提供策略支持[13]。隨著現(xiàn)代空間信息技術的快速發(fā)展，空間流行病學正越來越多地被應用于各種疾病研究中，在公共衛(wèi)生領域發(fā)揮著重要作用[14-15]。

2006 年以來，全國逐步開始實施地方棘球蚴病防治項目，通過采取“以控制傳染源為主，中間宿主防控與病人查治相結合”的綜合防治策略，棘球蚴病的防治工作取得了一定成效[6]。本研究發(fā)現(xiàn)，在青藏高原東部地區(qū)青海省人群棘球蚴病發(fā)病率略高于四川省，這與之前研究結果一致[2]。調查的青藏高原東部地區(qū)43 個縣（市）2016—2020 年發(fā)病率呈下降趨勢，2020 年發(fā)病率已低至0.013%，說明近年來棘球蚴病防治工作取得了顯著成效。棘球蚴病流行趨勢趨于平穩(wěn)，但部分縣偶有回升，出現(xiàn)個別新發(fā)病例，提示仍需加強對棘球蚴病的動態(tài)監(jiān)測。同時，2016—2020 年棘球蚴病發(fā)病率＞0.1%的地區(qū)也逐年減少，2016—2018年青藏高原東部地區(qū)棘球蚴病發(fā)病率較高的地區(qū)主要分布在中部交界地區(qū)，包括青海省曲麻萊縣、稱多縣、瑪多縣和達日縣，四川省石渠縣、德格縣、白玉縣和阿壩縣；2019年發(fā)病率較高地區(qū)為青海省治多縣和雜多縣，四川省阿壩縣；2020年發(fā)病率較高地區(qū)為青海省久治縣和四川省紅原縣?？赡苡捎?016年開始實施《全國包蟲病監(jiān)測方案》，部分縣剛開展人群監(jiān)測工作，新發(fā)現(xiàn)的病例中可能包含了既往發(fā)病但未及時檢查出來的病例。2017年青海省達日縣新發(fā)現(xiàn)121例病例，可能由于當時基層專業(yè)人員能力有限，該縣2016年未完成全部監(jiān)測工作，部分監(jiān)測任務2017年才得以完成。

大量研究顯示，棘球蚴病的分布具有一定的空間聚集性。Schmidberger等[16]對1992—2016年德國棘球蚴病數(shù)據(jù)進行空間自相關性分析，結果顯示整個德國的病例分布具有聚集性，局部空間自相關分析結果表明巴登-符騰堡州東南部地區(qū)和巴戈利亞州西南部地區(qū)為高危區(qū)。He 等[17]對四川省2016—2019年人群棘球蚴病流行情況進行調查，全局空間自相關分析顯示存在空間正相關，局部空間自相關分析發(fā)現(xiàn)“高-高”聚集區(qū)在四川省西北部，“低-低”聚集區(qū)在四川省東南部。趙月等[18]對四川省甘孜州2014年棘球蚴病的空間分析結果顯示棘球蚴病空間分布存在聚集性。本研究發(fā)現(xiàn)2016—2018 年青藏高原東部地區(qū)棘球蚴病發(fā)病率不具有空間聚集趨勢，2019—2020年存在空間正相關，具有明顯的空間聚集特征，存在“高-高”“高-低”和“低-高”聚集區(qū)，可能由于2016—2018年研究區(qū)域的篩查地區(qū)不完整，有15個縣（市）的數(shù)據(jù)缺失。局部Getis′G熱點分析結果與局部空間自相關分析的“高-低”和“高-高”聚集地區(qū)較為一致，2019年熱點區(qū)域位于青海省治多縣、雜多縣、玉樹市，2020年熱點區(qū)域位于四川省阿壩縣、紅原縣、理縣、若爾蓋縣、松潘縣、黑水縣、馬爾康市?？赡苡捎谶@些地區(qū)位于青南高原，寒冷潮濕的氣候環(huán)境有利于棘球絳蟲蟲卵存活，且多房棘球絳蟲的中間宿主嚙齒類動物數(shù)量多、分布廣。此外這些地區(qū)以畜牧業(yè)為主，自然環(huán)境和牲畜類型都有利于棘球絳蟲完成生活史循環(huán)，造成棘球蚴病的流行。甘孜州作為棘球蚴病的重點防控區(qū)域，發(fā)病率逐漸降低，而阿壩州近年來出現(xiàn)個別新發(fā)病例，這可能是熱點區(qū)域轉移的主要原因。

綜上，青藏高原東部地區(qū)棘球蚴病防治工作取得了一定成效，但仍存在空間分布聚集性。對于“高-低”和“高-高”聚集區(qū)域和熱點區(qū)域，應強化棘球蚴病的防治工作，加強犬只管理，開展健康宣傳，同時在既往人群篩查的基礎上，加強重點區(qū)域的全人群普查。此外，本研究以縣級為單位，研究結果僅可反映縣級聚集情況，后續(xù)可從更小的尺度進行研究；另外，此次研究不涉及影響因素，今后可開展相關研究。