2001—2018年西藏自治区包虫病流行趋势分析_赵哲.pdf

中国国境卫生检疫杂志2023年2月第46卷第1期Chinese Frontier Health Quarantine Feb.2023，Vol 46，No.1Analysis of epidemic trend of echinococcosis in Tibet AutonomousRegion，2001-2018ZHAO Zhe*，ZHANG Shuyuan，ZHAO Chenxi，ZHANG Wenkai，CIREN Panduo，LI Li，LIU Kun，LONG Yong*Department of Public Health，Gansu University of Chinese Medicine，Lanzhou，Gansu 730000，ChinaAbstract：ObjectiveTo analyze the epidemic trend of echinococcosis in Tibet Autonomous Region，and to providebasic data for the surveillance，prevention and control of echinococciosis.MethodsThe data of echinococcosisincidence of each province，autonomous region and municipality in China from 2001 to 2018 was collected from thePublic Health Science Data Sharing Center，and the annual incidence trend in Tibet was estimated by Joinpointregression analysis model.The data was analyzed by two time periods，2001-2009 years，2010-2018 years.Thespatial trend of echinococcosis was analyzed by Global Moranlys I statistic and hotspot analysis was carried outusing Geographic Information System（GIS）.ResultsFrom 2001 to 2018，the incidence of echinococcosis in Tibetshowed an increasing trend，with an average annual change percentage（AAPC）of 36.6.During the same period，thenational data showed a trend of increasing first and then slowly，with an average annual change percentage of 75.31and 5.45.The average incidence in Tibet in the second period was higher（4.46/100000）than that in the first period（0.14/100000）.During the both periods，Global Morans I analysis of the average incidence in China was 0.108（Z=2.090，P0.05）and 0.154（Z=2.948，P0.05），there was significant positive spatial autocorrelation.Local Morans Iidentified that Tibet appeared a high-high correlation model in the second period.Hotspot analysis revealed thatTibet was identified as a hotspot region in both periods.ConclusionThe incidence of echinococcosis wasincreased from 2001 to 2018 in Tibet.It was suggested to strengthen the scientific monitoring of echinococcosis inTibet，strictly implement prevention and control measures，and formulate relevant emergency response prescheme.Keywords：Echinococcosis；Incidence；Regression analysis；Spatial analysis疾病监测20012018 年西藏自治区包虫病流行趋势分析赵哲1，2，张淑媛2，3，赵晨锡2，4，张文凯2，3，次仁潘多5，李丽6，刘昆2，龙泳21.甘肃中医药大学公共卫生学院，甘肃 兰州730000；2.空军军医大学军事预防医学系军队防疫与流行病学教研室，特殊作业环境危害评估与防治教育部重点实验室；3.陕西中医药大学公共卫生学院；4.包头医学院公共卫生学院；5.亚东县卫生服务中心；6.西藏大学附属阜康医院摘要：目的分析西藏自治区包虫病的流行趋势，为包虫病的监测和防治提供数据资料。方法收集公共卫生科学数据共享中心20012018年全国各省、自治区、直辖市的包虫病发病率数据，利用Joinpoint回归分析模型估计西藏年发病率变化趋势；以20012009年为第一时间段，20102018年为第二时间段，用地理信息系统分析空间变化趋势。结果20012018年西藏包虫病发病率呈递增趋势，平均年度变化百分比为36.6；同期全国呈先快后慢递增趋势，平均年度变化百分比分别为75.31和5.45。西藏在第二时间段平均发病率（4.46/10万）于第一时间段（0.14/10万），第一时间段和第二时间段全国各省、自治区、直辖市平均发病率Morans I指数分别为0.108（Z=2.090，P0.05）和0.154（Z=2.948，P0.05），均存在空间正相关性；局部空间自相关表明西藏在第二时间段呈“高-高”关联模式；GetisG结果显示西藏在两个时间段均为热点地区。结论西藏的包虫病发病率20012018年间呈递增趋势，建议加强科学监测，严格落实个性化防控措施，制定相关应急处置预案。关键词：包虫病；发病率；回归分析；空间分析中图分类号：R181.3文献标识码：ADOI：10.16408/j.1004-9770.2023.01.021基金项目：军队后勤科研计划重点项目（BWS20J020）；军队医学科技青年培育计划拔尖项目（21QNPY087）通信作者：龙泳，E-mail：92中国国境卫生检疫杂志2023年2月第46卷第1期Chinese Frontier Health Quarantine Feb.2023，Vol 46，No.1包虫病即棘球蚴病（echinococcosis），是棘球绦虫的幼虫寄生于人体组织引起的人兽共患病1，目前发现有5种绦虫可致人和动物患病：细粒棘球绦虫、泡状棘球绦虫、石渠棘球绦虫、伏氏棘球绦虫和少节棘球绦虫2。包虫病在全球分布广泛，南美洲国家阿根廷、非洲的南非共和国以及亚洲的伊朗和土耳其都有包虫病的病例报告3-6。我国是世界上受包虫病危害较大的国家之一，主要分布于新疆、青海、宁夏、甘肃、内蒙古等西部省或自治区的农牧区域，对居民健康和畜牧业的发展有很大影响。我国的包虫病主要有两种，一种是由细粒棘球绦虫引起的囊型包虫病，另一种是由多房棘球绦虫引起的泡型包虫病，其中囊型包虫病占90%以上2。西藏自治区位于我国西南边陲，牧区广阔，牛羊成群，有包虫病流行的传播条件。近年来，西藏农牧业快速发展，包虫病发病人数逐渐增多，但相关文献报道较少。本研究通过分析20012018年间西藏包虫病发病率等数据，了解西藏包虫病流行趋势，为包虫病科学防控提供基础依据。1资料与方法1.1资料来源1.1.1疫情资料20012018年全国各省、自治区、直辖市的包虫病数据资料来源于公共卫生科学数据共享中心（https：/），数据库标识符为8a47b989-40a6-4492-83e8-8eb7c9d11abb。各省、自治区、直辖市人口资料来源于中国统计年鉴。1.1.2地图资料以中国电子地图为基础地图，审图号GS（2019）1822号。地图数据来源于国家基础地理信息数据库。1.2方法1.2.1Joinpoint回归模型分析采用美国国家癌症研 究 院（National Cancer Institute，NCI）开 发 的Joinpoint Regression Program 4.9.0.1软件7，利用发病数和人口数计算发病率，再近似估计标准误。计算20012018年全国和西藏包虫病发病率的年均变化百分比（average annual percent change，AAPC）和年度变化百分比（annual percent change，APC），由此分析西藏和全国20012018年间包虫病的发病趋势。该软件可以模拟包虫病的发病趋势变化，转折点表示上升趋势和下降趋势的连接点。APC0时，表示呈上升趋势；APC0时，表示呈下降趋势；若无转折点，表示APC=AAPC。用95%CI值判断转折点变化是否有意义8。P0表示整体区域在空间上存在正相关性，值越大证明空间聚集性越强；I=0表示整体区域呈随机分布，无空间聚集性；I0表示整体呈现负相关，值越小证明空间差异性越大9。局部Morans I自相关可得出高-高（H-H）、低-低（L-L）、低-高（L-H）高-低（H-L）4种聚集模式，以显示空间正相关性和负相关性。空间局部GetisG热点分析用来表示空间自相关程度，以发现热点区（高发病聚集区）和冷点区（低发病聚集区）10。2结果2.1包虫病发病率 Joinpoint 趋势分析20012018年西藏包虫病发病率Joinpoint趋势分析显示，APC=36.60，95%CI=28.445.3，P0.001，呈递增趋势；同时期全国在2007年出现转折点，20012007年递增较快（APC=75.31，95%CI=15.1167.0，P=0.013）；20072018年递增缓慢（APC=5.45，95%CI=2.98.1，P0.001）。见表1。表120012018年包虫病发病率趋势分析Tab.1Analysis of the incidence trend of echinococcosis，2001-2018区域起止年份APC（%）95%CI值AAPC95%CI值P值西藏2001201836.6028.445.336.628.445.30.001全国2001200775.3115.1167.026.210.244.50.013200720185.452.98.10.00193中国国境卫生检疫杂志2023年2月第46卷第1期Chinese Frontier Health Quarantine Feb.2023，Vol 46，No.1图3包虫病热点分布Fig.3Hot spot distribution of echinococcosis图2包虫病局部空间自相关Fig.2Local spatial autocorrelation of echinococcosis图1全国包虫病发病率Fig.1Incidence of echinococcosis in China2.2包虫病发病率的空间分析2.2.1发病率分析西藏的包虫病平均发病率在第二时间段为4.46/10万，高于第一时间段（0.14/10万），在全国范围内处于中下水平。见图1。2.2.2空间自相关分析全局空间自相关分析结果显示，第一时间段全国各省、自治区、直辖市包虫病发病率的Morans I值为0.108（Z=2.090，P0.05），存在正向空间自相关性；第二时间段的Morans I值为0.154（Z=2.948，P2.58）地区均有西藏，说明西藏是热点地区，是包虫病高发聚集区。见图3。3讨论包虫病是我国法定报告的丙类传染病，对民众生活的影响虽不像甲类传染病那么严重，但我国近年来仍有新发病例发生，仍需加强包虫病的监测和防护。中国疾病预防控制中心调查了20162019年全国包虫病的发病情况，发现全国包虫病流行基本稳定，但发病数逐渐增多11-12。西藏自治区拉萨市13和昌都市14对当地的包虫病调查表明，包虫病在所报道地区流行形势严峻。本研究分析了全国包虫病的流行特征，发现西藏自治区的包虫病发病率逐步增高，继而对西藏自治区的流行情况做了调查，并在空间层面对两个时间段的发病情况做出空间自相关分析、热点区域分析以及趋势分析，将数据模型分析和空间分析相结合，更全面地分析包虫病在西藏的流行趋势。这是本研究与之前包虫病的调查的不同之处。3.120012018 年西藏包虫病发病情况分析本研究显示西藏包虫病年发病率逐渐上升，而包虫病的发生和传播与农牧业息息相关。包虫病的宿主以牛、羊及犬类为主，西藏的牧区和农业区地域广阔，牛、羊、犬较多，所以当地居民接触到包虫虫卵的机会较多，包虫病感染概率增加。另外，藏族人民平时收集动物粪便作为生活燃料，可能也是感染概率增加的一个原因。3.220012018 年西藏包虫病空间趋势分析用地理信息系统（geographic information system，GIS）对全国两个时间段包虫病发病情况作空间自相关分析，显示第一时间段包虫病发病率虽存在正向空间自相关性（P0.05），但西藏不是聚集地区。局部空间自相关分析表明第二时间段西藏呈“高-高”聚集模94中国国境卫生检疫杂志2023年2月第46卷第1期Chinese Frontier Health Quarantine Feb.2023，Vol 46，No.1式。GetisG分析显示，两个时间段热点地区均有西藏，表明西藏较易发生和流行包虫病。此外，西藏之所以是热点地区，可能与其区域面积广而医疗条件及了解防控知识能力相对短缺有一定的关联；居民牧业活动频繁，与牲畜接触机会增多，包虫病感染率随之增加。针对西藏包虫病近年发病率和发病人数日趋递增的情况，应加强个人卫生与环境卫生，规范犬类管理，加强包虫病防治的宣传教育。同时加大包虫病的监测力度，出现病例及时治疗，控制包虫病的流行。参考文献1孟庆贺.我国包虫病的诊断和治疗进展J.中国城乡企业卫生，2014，29（5）：16-18.Meng QH.Advances in diagnosis and treatment of echinococ-cosis in China J.Chinese Journal of Urban and Rural Enter-prise Hygiene，2014，29（5）：16-18.（in Chinese）2金珍.包虫病的研究进展J.畜禽业，2020，31（4）：118.Jin Z.Research progress of echinococcosis J.Livestock andPoultry Industry，2020，31（4）：118.（in Chinese）3Uchiumi L，Mujica G，Araya D，et al.Prevalence of human cys-tic echinococcosis in the towns of orquinco and Ramos Mexiain Rio Negro Province，Argentina，and direct risk factors for infec-tionJ.Parasit Vectors，2021，14（1）：262.4Miambo RD，Afonso SMS，Noormahomed EV，et al.Echinococ-cosis in humans and animals in Southern Africa developmentcommunity countries：a systematic review J.Food WaterborneParasitology，2020，20：e00087.5Mahmoudi S，Mamishi S，Banar M，et al.Epidemiology ofechinococcosis in Iran：a systematic review and meta-analysisJ.BMC Infectious Diseases，2019，19（1）：929.6Eroglu F，Dokur M，Ulu Y.The molecular epidemiology of cys-tic and alveolar echinococcosisin Southeast TurkeyJ.IranianJournal of Parasitology，2021，16（2）：327-335.7National Cancer Institute.Joinpoint trend analysis soft-ware，version 4.4EB/OL.（2017-01-04）2018-03-05.https：/surveil-lance.cancer.gov/joinpoint/.8Kim HJ，Fay MP，Feuer EJ，et al.Permutation tests for join-point regression with applications to cancer rates J.Statisticsin Medicine，2000，19（3）：335-3519姜晓峰，冯克民，郝慧霞，等.2013-2018年内蒙古自治区棘球蚴病报告病例空间分布特征J.中国血吸虫病防治杂志，2020，32（3）：268-272.Jiang XF，Feng KM，Hao HX，et al.Spatial distribution of re-ported echinococcosis cases in Inner Mongolia Autonomous Re-gion from 2013 to 2018 J.Chinese Journal of SchistosomiasisControl，2020，32（3）：268-272.（in Chinese）10张见增，泽仁拉措，洛桑扎西，等.2017-2019年西藏丁青县包虫病流行特征及空间聚集性分析J.中华卫生杀虫药械，2021，27（3）：249-253.ZhangJZ，ZeRLaCo，LuosangTS，etal.EpidemiologicalcharacteristicsandspatialclusteringofechinococcosisinDingqing County，Tibet from 2017 to 2019 J.Chinese Journalof Hygienic Insecticides&Equipments，2021，27（3）：249-253.（inChinese）11魏思慧，伍卫平，韩帅，等.2016-2017年全国棘球蚴病监测结果分析J.中国病原生物学杂志，2020，15（8）：924-928.WeiSH，WuWP，HanS，etal.Analysisofechinococcosissurveillance results in China from 2016 to 2017 J.Journal ofPathogen Biology，2020，15（8）：924-928.（in Chinese）12蒉嫣，伍卫平，韩帅，等.2018-2019年全国棘球蚴病监测分析J.中国病原生物学杂志，2021，16（9）：1025-1029.LeiY，WuWP，HanS，etal.Surveillanceandanalysisofechinococcosis in China from 2018 to 2019 J.Journal ofPathogen Biology，201，16（9）：1025-1029.（in Chinese）13格桑卓嘎，汤晓佳，王政和，等.2017年拉萨地区人群棘球蚴病流行现状与特征J.中国寄生虫学与寄生虫病杂志，2020，38（1）：36-40.Gesang ZG，Tang XJ，Wang ZH，et al.Prevalence and charac-teristics of echinococcosis in Lhasa in 2017J.Chinese Journalof Parasitology and Parasitic Diseases，20，38（1）：36-40.（in Chi-nese）14贡桑曲珍，李斌，陈伟奇，等.昌都市棘球蚴病流行现状分析J.中国寄生虫学与寄生虫病杂志，2018，36（1）：68-74.Gong SQZ，LI B，Chen WQ，et al.Epidemiology of echinococco-sis in Changdu City J.Chin J Parasitol Parasit，2018，36（1）：68-74.（in Chinese）收稿日期：2022-06-0695