2006-2017年全國(guó)布魯氏菌病流行特征變化分析

張夏男，張洪偉，紀(jì) 穎，王秋芳，謝立璟，葉 萌

布魯氏菌病(簡(jiǎn)稱布病)是由布魯氏菌引起的人獸共患的一種慢性傳染病[1]，是我國(guó)的乙類傳染病。布病的主要臨床表現(xiàn)是發(fā)熱、多汗、乏力，肌肉和關(guān)節(jié)疼痛等[2]，布病在急性期如治療不徹底，極易轉(zhuǎn)為慢性，病程長(zhǎng)達(dá)數(shù)年甚至幾十年，甚至喪失勞動(dòng)力[3]。我國(guó)大部分省市報(bào)告人間布魯氏菌病的病例數(shù)有升高趨勢(shì)[4-6]，本研究擬采用GeoDa-1.14.0.0軟件、集中度法及圓形分布法、構(gòu)成比對(duì)全國(guó)31個(gè)省市布病流行特征的變化趨勢(shì)進(jìn)行分析。

1 資料與方法

1.1資料來(lái)源 2006-2017年布病疫情資料來(lái)自公共衛(wèi)生科學(xué)數(shù)據(jù)共享中心網(wǎng)站(http://www.phsciencedata.cn)，該網(wǎng)站提供的數(shù)據(jù)庫(kù)涵蓋了39種傳染病、慢性非傳染性疾病、健康危險(xiǎn)因素、生命登記等內(nèi)容，其中傳染病收集的主要是自2004年傳染病網(wǎng)絡(luò)直報(bào)以來(lái)報(bào)告的全部傳染病數(shù)據(jù)，主要包括分地區(qū)、年齡組的發(fā)病人數(shù)、發(fā)病率、死亡人數(shù)和死亡率等。為了分析2006-2017年布病病例的空間和時(shí)間屬性，將每3年的新發(fā)病例數(shù)進(jìn)行合并形成4個(gè)時(shí)間段的數(shù)據(jù)，第1時(shí)間段為2006-2008年，第2時(shí)間段為2009-2011年，第3時(shí)間段為2012-2014年，第4時(shí)間段為2015-2017年。每個(gè)時(shí)間段各省市的累積發(fā)病率以3年新發(fā)病例數(shù)的總和作為分子，各省3年平均人口數(shù)為分母計(jì)算。

1.2 分析方法

1.2.1地區(qū)分布及時(shí)空分布 將4個(gè)時(shí)間段的累積發(fā)病率(1/10萬(wàn))數(shù)據(jù)導(dǎo)入GeoDa-1.14.0.0，以全國(guó)基礎(chǔ)地理數(shù)據(jù)為地圖，建立Rook權(quán)重矩陣文件，呈現(xiàn)全國(guó)31個(gè)省市累積發(fā)病率的箱地圖，同時(shí)進(jìn)行全局空間自相關(guān)分析和局部空間自相關(guān)分析。全局空間自相關(guān)分析使用全局Moran’sⅠ指數(shù)反映空間要素屬性值的聚集程度；局部空間自相關(guān)分析反映空間單元與相鄰單元疫情的相關(guān)程度，可用LISA集聚圖識(shí)別四類聚集區(qū)域：高-高聚集、高-低聚集、低-高聚集和高-低聚集。將局部Moran’sⅠ指數(shù)通過(guò)10%顯著性水平檢驗(yàn)的區(qū)域分類顯示于LISA集聚圖中。

1.2.2季節(jié)性分析 采用集中度法的M值和圓形分布法的R值對(duì)4個(gè)時(shí)間段布病的季節(jié)性特點(diǎn)進(jìn)行分析。M值和R值的計(jì)算公式及評(píng)定標(biāo)準(zhǔn)參見(jiàn)文獻(xiàn)[7-8]。評(píng)定標(biāo)準(zhǔn)如下：越接近1季節(jié)性越強(qiáng)，越接近0季節(jié)性越弱，0.9以上代表嚴(yán)格季節(jié)性，0.7～0.9代表很強(qiáng)的季節(jié)性，0.5～0.7代表明顯的季節(jié)性，0.3～0.5代表一定的季節(jié)性，0.3以下代表季節(jié)性較弱。

1.2.3年齡分布 利用構(gòu)成比分析各年齡段分布的變化情況，并用趨勢(shì)性χ2檢驗(yàn)分析老少組(14歲及以下和65歲及以上)和非老少組(15～64歲組)在4個(gè)時(shí)間段的構(gòu)成比有無(wú)統(tǒng)計(jì)學(xué)差異，檢驗(yàn)水準(zhǔn)α=0.05。

2 結(jié) 果

2.1全國(guó)布病疫情現(xiàn)狀 2006-2017年全國(guó)布病疫情總體呈上升趨勢(shì)，年平均發(fā)病率為2.836/10萬(wàn)，發(fā)病率最低為2006年(1.454/10萬(wàn))，最高為2014年(4.223/10萬(wàn))，4個(gè)時(shí)間段的年平均發(fā)病率分別是1.686/10萬(wàn)、2.691/10萬(wàn)、3.457/10萬(wàn)和3.469/10萬(wàn)。如圖1所示，布病病例數(shù)與發(fā)病率曲線完全吻合。

圖1 2006-2017年全國(guó)布病發(fā)病曲線Fig.1 Brucellosis incidence curve in China from 2006 to 2017

2.2各省市布病累積發(fā)病率 北方省份和南方省份的界定參照文獻(xiàn)[9]。16個(gè)北方省份和15個(gè)南方省份中，第1時(shí)間段江西省、海南省、貴州省均沒(méi)有報(bào)告病例，之后的3個(gè)時(shí)間段31省市均有報(bào)告病例。高發(fā)省份(內(nèi)蒙古自治區(qū)、山西省、黑龍江省、吉林省、河北省)布病的累積發(fā)病率在第4個(gè)時(shí)間段有所下降，但寧夏回族自治區(qū)和新疆維吾爾自治區(qū)的累積發(fā)病率出現(xiàn)明顯提高，如圖2所示。

圖2 各省市布病在4個(gè)時(shí)間段的累積發(fā)病率Fig.2 Cumulative incidence rate of brucellosis in four time periods in various provinces

2.3布病累積發(fā)病率的空間分布 4個(gè)時(shí)間段的箱式地圖顯示全國(guó)布病累積發(fā)病率的空間分布。第1時(shí)間段的累積發(fā)病率在0至113.41/10萬(wàn)之間，離群省為內(nèi)蒙古、黑龍江、吉林、河北、山西和陜西。第2時(shí)間段的累積發(fā)病率在0.002 3/10萬(wàn)至205.69/10萬(wàn)之間，離群省為內(nèi)蒙古、黑龍江、吉林、山西、河北和寧夏。第3時(shí)間段的累積發(fā)病率在0.029 4/10萬(wàn)至124.76/10萬(wàn)之間，離群省為內(nèi)蒙古、黑龍江、山西、寧夏和新疆。第4時(shí)間段的累積發(fā)病率在0.08/10萬(wàn)至100.94/10萬(wàn)之間，離群省為內(nèi)蒙古、黑龍江、山西、寧夏和新疆。

2.4全局空間自相關(guān)分析 以4個(gè)時(shí)間段31個(gè)省市布病累積發(fā)病率為單位進(jìn)行全局Moran’s Ⅰ空間自相關(guān)分析，Moran’s Ⅰ指數(shù)分別為0.212、0.190、0.377和0.264(Z值分別為3.800、4.368、4.108和2.861，均有P<0.05)，表明4個(gè)時(shí)間段31個(gè)省市布病累積發(fā)病率存在空間正相關(guān)，即相鄰省市之間有相似的布病發(fā)病風(fēng)險(xiǎn)，隨著時(shí)間的推移，發(fā)病風(fēng)險(xiǎn)有上升趨勢(shì)。

2.5局部空間自相關(guān)分析 4個(gè)時(shí)間段的LISA聚集圖主要表現(xiàn)為“高高”“低低”聚集，且只有“高高”“低低”“低高”3種局域相關(guān)類型，無(wú)“高低”聚集區(qū)。前3個(gè)時(shí)間段“高高”聚集發(fā)生在中國(guó)的北方區(qū)域主要包括內(nèi)蒙古、黑龍江、吉林、河北、山西等5省市；第4個(gè)時(shí)間段“高高”聚集發(fā)生在中國(guó)的北方區(qū)域內(nèi)蒙古、寧夏和甘肅，并且隨著時(shí)間的推移寧夏從“低高”聚集轉(zhuǎn)變?yōu)椤案吒摺本奂?，提示寧夏布病發(fā)病率有顯著的升高。甘肅先從“非顯著性區(qū)域”轉(zhuǎn)變?yōu)椤暗透摺本奂?，最后轉(zhuǎn)變?yōu)椤案吒摺本奂８拭C聚集類型的轉(zhuǎn)變除了與甘肅本身布病發(fā)病率提高有關(guān)外，還與寧夏和新疆發(fā)病率的提高有關(guān)?！暗透摺本奂霈F(xiàn)在“高高”聚集周圍，陜西省在4個(gè)時(shí)間段均為“低高”聚集?！暗偷汀本奂l(fā)生在中國(guó)的南方區(qū)域主要包括云南、貴州、廣西、廣東、湖南、江西和浙江。這種空間聚集現(xiàn)象顯示南北方布病發(fā)病率的差異依然存在。

2.6時(shí)間分布 隨著時(shí)間的推移，M值和R值均有減小的趨勢(shì)，逐漸低于0.3，見(jiàn)表1，提示布病從具有一定季節(jié)性逐漸轉(zhuǎn)變?yōu)榧竟?jié)性不明顯。

表1 全國(guó)布病在4個(gè)時(shí)間段的集中度和圓形分布法分析Tab.1 Analysis of concentration and circular distribution of brucellosis in four time periods in China

2.7年齡分布 因未獲得2017年分年齡組的病例數(shù)，所以第4個(gè)時(shí)間段只統(tǒng)計(jì)2015和2016年的分年齡組的病例數(shù)。4個(gè)時(shí)間段內(nèi)，均以35～54歲組的發(fā)病人數(shù)居多，老少組的構(gòu)成比分別是7.3%、7.4%、9.9%和12.7%。趨勢(shì)性χ2檢驗(yàn)的結(jié)果表明，隨時(shí)間變化，老少組和非老少組構(gòu)成比的趨勢(shì)變化差異有統(tǒng)計(jì)學(xué)意義(χ2=1 946.95，P<0.001)。

3 討 論

本研究以每3年作為一個(gè)時(shí)間段研究全國(guó)布病在2006-2017年流行特征的變化情況，既反映了布病的長(zhǎng)期變化趨勢(shì)，也避免了布病發(fā)病的逐年波動(dòng)性，所以具有較強(qiáng)的實(shí)際意義。結(jié)果表明2006-2017年全國(guó)布病疫情總體呈上升趨勢(shì)，但自2015年來(lái)布病的年發(fā)病率出現(xiàn)下降趨勢(shì)，尤其是布病的高發(fā)省份，如內(nèi)蒙古、黑龍江、山西。《國(guó)家中長(zhǎng)期動(dòng)物疫病防治規(guī)劃(2012-2020年)》[10]和《國(guó)家布魯氏菌病防治計(jì)劃(2016-2020年)》[11]的頒布和實(shí)施對(duì)遏制布病疫情上升態(tài)勢(shì)起到了積極的作用。而南方省份布病累積發(fā)病率明顯低于北方省份，但處于持續(xù)的上升態(tài)勢(shì)。這可能與南方省份作為新發(fā)疫區(qū)，缺乏經(jīng)驗(yàn)有一定的關(guān)系。

空間分析方法可用來(lái)分析肺結(jié)核、病毒性肝炎、手足口病等多種傳染病的流行病學(xué)特征[12-13]。不同研究者在不同省份根據(jù)該方法均發(fā)現(xiàn)布病分布存在空間自相關(guān)[5-6, 14]，本研究顯示，2006-2017年全國(guó)布病在4個(gè)時(shí)間段的Moran’s Ⅰ指數(shù)為0.190～0.377，說(shuō)明存在明顯聚集性，主要集中在內(nèi)蒙古附近的北方省份，與2006-2012年的研究結(jié)果一致[15]。LISA集聚圖顯示除新疆維吾爾自治區(qū)外，高高聚集區(qū)與高發(fā)病率區(qū)基本一致；另外寧夏回族自治區(qū)和甘肅省隨著時(shí)間的變化逐漸演變?yōu)楦吒呔奂瘏^(qū)，陜西省在4個(gè)時(shí)間段均為低-高聚集區(qū)。低-高聚集區(qū)要采取積極的措施控制牲畜及其產(chǎn)品從高發(fā)區(qū)域流入，從而預(yù)防疫情蔓延。我國(guó)南方省份雖屬于低-低聚集區(qū)，但病例數(shù)卻逐年遞增，建議這些區(qū)域也不能放松警惕，加強(qiáng)與北方省份的交流，積極探索布病的防控措施。

集中度法的M值和圓形分布的R值在4個(gè)時(shí)間段呈現(xiàn)遞減趨勢(shì)，說(shuō)明布病的季節(jié)性越來(lái)越弱。本研究中，2012-2014年的R值為0.273 5，與2012-2016年的R值接近[16]。2004-2010年濟(jì)南市的M值和R值為0.44左右[7]，2004-2012年山東省的M值為0.375 3[17]，2004-2015年北京市M值和R值為0.3左右[18]，提示不同省份和不同時(shí)間段的M值和R值存在差異。

4個(gè)時(shí)間段內(nèi)，均以35～54歲組的發(fā)病人數(shù)居多，但隨著時(shí)間的推移，老少組的構(gòu)成比升高。2011-2016年福建省的病例中有11.91%是散居兒童、學(xué)生、教師、餐飲食品業(yè)等，14起聚集性疫情[6]。2013-2018年甘肅省學(xué)生占比為2.42%[14]。這些證據(jù)提示布病從職業(yè)人群向一般人群蔓延，可能與感染方式的多樣化有關(guān)，所以在布病防控過(guò)程中既要關(guān)注職業(yè)人群也要關(guān)注一般人群。

綜上所述，全國(guó)布病的南北差異依然存在，季節(jié)性逐漸趨于不明顯，老少人數(shù)有增多趨勢(shì)。布病防控形勢(shì)依然比較嚴(yán)峻，既要關(guān)注職業(yè)人群也要關(guān)注一般人群，高發(fā)省份依然要認(rèn)真落實(shí)各項(xiàng)防控措施，低發(fā)省份加強(qiáng)監(jiān)測(cè)預(yù)防。

利益沖突：無(wú)