1950-2015年中国布鲁菌病空间流行病学分析

丁言稳,刘 婷,蔡 奔,贾贤杰

(蚌埠医科大学 公共卫生学院 流行病与卫生统计学教研室,安徽 蚌埠 233030)

布鲁菌病(以下简称布病)是一种由布鲁菌群侵入机体后通过淋巴途径传播到血液,再通过血液导致全身感染引起的人畜共患的传染变态反应性疾病[1]。该病病人负重骨关节最易受累,脊柱最常见,其次是髋、膝、骶髂、肩和踝关节等,同时可累及神经系统[2]。急性期以发热、多汗、乏力和关节痛为主要表现,可重复感染或转为慢性,严重者可致残疾。同时布病发病症状缺乏特异性,临床诊断鉴别难度大,易发生误诊误治[2]。21世纪以来,我国布病发病率持续上升,自2006年始,布病稳居甲乙类传染病报告发病数前10位,布病疫情呈现由北向南从传统牧区、半牧区向农区扩散,成为危及居民健康的重要公共卫生问题之一[3]。本文在空间流行病学基础上回顾性分析1950-2015年布病趋势,通过总结布病趋势分析,以期加深临床医生对布病认知,减少误诊率,做到早诊断、早治疗,改善预后,降低布病慢性化比例[4-5],减轻群众健康的损害。现作报道。

1 资料与方法

1.1资料来源

数据资料来自国家人口健康共享杯服务平台公共卫生科学数据中心提供的《1950-2015年布病数据》。使用资料内容包括1950-2009年全国各年、各月布病发病人数、死亡人数、发病率和死亡率;2010-2015年全国各省(自治区、直辖市)布病月发病率。

1.2方法

采用Excel2010构建1950-2015年布病数据数据库,并绘制60年间布病发病率、死亡率和月平均发病例数趋势图作直观分析。通过Arcgis10.2软件空间统计工具中的聚类分析模块,对2000-2015年布病发病率作空间自相关分析,其中,空间自相关分析包括全局空间自相关分析和局部空间自相关分析。

1.2.1 全局空间自相关分析 该分析是对整个区域属性值空间特征的描述,表示同一个变量在不同空间位置上的相关性[6-7]。本文通过全局自相关分析描述全国各省(区、市)布病发病的总体分布情况,反映指标为全局空间自相关系数Moran′sI,其取值一般在-1～1,>0表示正相关,<0表示负相关,=0表示不相关;Moran′sI值越趋近于1,总体空间差异越小,值越趋近于-1,总体空间差异越大[8]。Moran′sI计算公式如下:

1.2.2 局部空间自相关性分析 反映指标为LocalMoran′sI,能够反映区域单元属于高值集聚或低值集聚的空间分布,空间关联模式可分为4种类型:高高关联、高低关联、低高关联和低低关联[8]。LocalMoran′sI计算公式如下:

式中n是参与分析的空间单元数,即全国各省市数;xi和xj分别表示在i省(区、市)和j省(区、市)的发病率,i≠j;wij是空间权重矩阵,表示i省(区、市)和j省(区、市)的邻近关系[7]。

2 结果

2.11950-2009年布病发病及死亡分布

1950-2009年60年间,我国共报告布病发病927 987例,死亡547例,年平均发病率为0.201 2(1/100 000),年平均死亡率为0.000 3(1/100 000)。其中1950-1979年间布病发病率波动幅度较大;1980-1998年间布病发病率较为平稳,无明显的升高或降低;1999-2009年间布病发病率持续上升(见图1)。死亡率方面,1950-1999年间我国布病死亡率相对较高;2000-2009年间虽发病率升高,但死亡率较低。其中布病发病致死较严重年份包括1957年(0.002)、1960年(0.002)、1964年(0.001)、1967年(0.004)、1973年(0.001)、1976年(0.001)(见图2)。

2.22000-2015年布病发病率地区分布

通过2000-2015年全国各省(市)布病发病率标准化处理结果可见,我国布病发病率地区分布总体趋势较为平稳,其中发病率较高的为内蒙古、黑龙江、山西、新疆等省、自治区,其次为吉林、河北、河南、辽宁等省,发病率较低的为安徽、江西、湖南、广东、广西、贵州、云南、四川等省、自治区。

2.3布病的空间分布

2.3.1 布病全局空间自相关分析 基于Arcgis10.2统计软件空间分析功能,计算2000-2015年全国布病年发病率全局Moran′sI值,结果显示,2000-2015年我国布病发病率呈现空间负相关性,其中2008-2013年全局空间自相关指数差异有统计学意义(P<0.05),表明其发病情况的全局空间分布关系是离散模式;2000-2007年、2014-2015年我国布病发病情况的全局空间分布关系则为随机模式(P>0.05)(见表1)。

表1 2000-2015年布病的全局空间自相关分析

2.3.2 布病局部空间自相关分析 对2000-2015年我国各省(区、市)布病发病率进行聚类和异常值分析,结果显示,高高空间关联模式存在于山西(2000年)、新疆(2014-2015年)和宁夏(2014-2015年);2002-2003年、2005-2013年内蒙古均存在高低空间关联模式;2001-2003年、2013-2014年山西存在高低空间关联模式;2004年西藏存在高低空间关联模式(见表2)。

表2 2000-2015年全国各省(区、市)布病发病率聚类和异常值分析

3 讨论

布病是由布鲁菌侵入机体引起的一种人畜共患的传染变态反应性疾病[9],是《中华人民共和国传染病防治法》规定的乙类传染病,近年我国布病流行呈逐年上升趋势,已成为严重的公共卫生问题[10]。本研究通过时间、地区和空间三个维度探究我国布病流行趋势,着重采用空间自相关分析法描述我国布病发病的空间分布状况。结果显示,我国1950-2009年60年间,共报告布病发病927 987例,死亡547例,年平均发病率为0.2012(1/100 000),年平均死亡率为0.000 3(1/100 000),1999-2009年间我国布病发病率持续上升,年平均增长率为0.232 1(1/100 000),2009年发病率高达2.696 9(1/100 000)。这与国内相关研究结论一致,LI等[11]明确从1999年开始我国布病发病急剧上升;白永飞等[12]研究指出,国家相关干预政策与布病疫情发展有关,影响布病的危险因素虽持续存在,但居民就诊意识和医疗机构报告敏感性提高,导致布病报告发病率迅速上升。我国布病报告发病例数平均每月约370例,1-6月逐月上升,6-11月逐月减少,布病发病高峰期是3-8月。这与杨丽等[13-14]研究一致,布病病例报告主要集中在春夏季节,可能与季节放牧、羔羊接生等季节性活动有关。

本研究显示,2000-2015年我国布病发病率地区分布总体较为稳定,其中布病发病率较高的是内蒙古、黑龙江、山西、新疆等省和自治区。这与国内相关研究一致,如LIN等[15]在文章中指出我国北方布病发病率较南方高;同时有研究[16]基于2009年全国布病检测数据指出内蒙古、山西省位于我国布病发病率前2位。本研究结果显示,我国布病发病呈现全局空间负相关性,高低空间关联模式多存在于内蒙古,高高空间关联模式多存在于山西、新疆和宁夏。ZHANG等[17]对我国布病2004-2010布病发病率进行空间分析发现,高风险地区主要集中在内蒙古地区及其周围以北的一些省份,这与本文的研究结果一致。

本文采用的数据为全国布病(粗)发病率,无发病人群信息并且有部分数据缺失,这会导致进行空间自相关分析时存在一定偏倚,使计算结果与实际情况存在一定偏差。此外,本研究采用的Moran′sI统计量容易受空间聚集区域大小的影响[18],可能会对真实结果的计算带来一定影响。今后的研究中有待进一步完善布病病例数据,如性别、年龄、职业等指标;同时考虑布病患病相关因素,如温度、湿度等生态指标对布病发病的影响。