【论著】2010—2024年中国人布鲁氏菌病的流行动态及时空扩散趋势分析
作者:张云飞 李新楼 徐强 牟笛 师悦 陈曦 周海健 秦天 阚飙 郑灿军 方立群文章来源:中华预防医学杂志, 2025, 59(6): 884-891. DOI: 10.3760/cma.j.cn112150-20250307-00181.摘要目的分析2010—2024年中国人布鲁氏菌病的流行动态及时空扩散趋势。方法通过“传染病报告信息管理系统”,收集中国31个省份2010年1月1日至2024年12月31日上报的人感染布鲁氏菌病病例资料,包含发病日期、性别、年龄、职业和地址等。采用Joinpoint回归和空间插值方法,从牧区/半牧区、其他区域,城市与农村分层分析人布鲁氏菌病疫情的时空动态与人群分布特征,探究疾病的流行趋势。结果2010—2024年中国牧区/半牧区共报告252 094例布鲁氏菌病病例,年均发病率为36.57/100 000±7.28/100 000;其他区域累计报告病例519 748例,年均发病率均数为2.54/100 000±0.74/100 000。牧区/半牧区人布鲁氏菌病发病率呈现下降-反弹-下降的趋势,2010—2017年为发病率下降阶段(APC=-7.20; PP=0.015),2021—2024年发病率再次呈现下降趋势(APC=-7.53; P=0.027);而在其他区域则呈现波动上升的趋势:2010—2015年为发病率上升阶段(APC=20.37; PPPP=0.735),仍保持较高的发病水平。全国农村发病率高于城市( P<0.001)。2010—2024年中国人布鲁氏菌病呈现从北方疫区向周边邻近区域逐渐扩散和南方地区跳跃性扩散的总体趋势。牧区/半牧区发病人群的年龄结构和其他区域呈现显著不同的分布特征,牧区/半牧区35~49岁人群发病率较高,其他区域55~64岁人群发病率较高。结论我国牧区/半牧区和其他区域的人布鲁氏菌病流行情况具有显著差异;并呈现从北方疫区向周边临近区域逐渐扩散和南方地区跳跃性扩散的总体趋势。布鲁氏菌病是由布鲁氏菌感染引起的一种人畜共患传染病 [ 1 ] 。人布鲁氏菌病又称地中海热、波状热或马耳他热,主要是人通过接触染疫动物或者其流产物、奶或奶制品等引发的感染。人布鲁氏菌病的潜伏期通常为5~60 d,大多数病例在暴露后2~4周出现症状,临床常见以发热、多汗、乏力、关节疼痛或全身肌肉疼痛为主的症状 [ 2 ] 。该病主要在中东、地中海、中美洲和南美洲等区域流行,呈全球性分布,最新研究表明,每年有1.6万~210万新发感染病例 [ 3 , 4 ] ,给人群健康和畜牧业发展带来了巨大影响和损失。自2000年以来,中国人布鲁氏菌病年度报告病例数快速上升,以内蒙古、甘肃和新疆等北方牧区为主,南方等非牧区发病率也在不断上升,并出现食源性感染 [ 5 ] 。张夏男等 [ 6 ] 通过分析2004—2017年我国南北方布鲁氏菌的发病趋势得出,北方省份发病率较高,南方省份出现持续增长的态势;杨向东等 [ 7 ] 通过分析2006—2018年云南省布鲁氏菌病报告病例发现,疾病流行不断上升,老疫点持续存在,新疫点不断增加,发病人群集中在20~69岁年龄组,职业以农民为主。宁夏人布鲁氏菌病疫情呈现上升-下降的波动趋势,空间分布上存在聚集性特征 [ 8 ] ;山西太原病例发病主要集中在2—8月,职业以农牧民为主 [ 9 ] 。可见,既往研究较多侧重于从全国、部分省、市等区域开展布鲁氏菌病的流行特征分析,缺乏该病空间扩散趋势的量化研究。此外,由于布鲁氏菌病的传播、扩散和畜牧养殖、畜牧产品消费和加工等因素高度相关,牧区/半牧区与其他地域、农村与城市等不同地域暴露感染方式的不同可能导致疾病流行的人群与时空差异,分析研究人布鲁氏菌病在不同类型区域的流行特征与传播规律对于该病的针对性防控具有重要价值。鉴于此,本研究通过系统分析2010—2024年中国人布鲁氏菌病的发病数据,从牧区/半牧区和其他区域、城市与农村等不同地域,开展疾病的时空动态分布及扩散特征、人群分布动态特征的分析研究,以期为制定针对性的防控策略提供数据支撑与科学依据。资料与方法一、资料来源(一)人感染布鲁氏菌病发病数据数据来源于中国“传染病报告信息管理系统”,收集我国31个省份(不包括中国香港、澳门和台湾地区)2010年1月1日至2024年12月31日上报的人感染布鲁氏菌病病例资料,包含发病日期、性别、年龄、职业和地址(区县行政代码,不包含个人隐私信息)。(二)人口数据及中国牧区/半牧区数据2010—2024年人口数据来源于国家统计局第七次全国人口普查数据。中国牧区/半牧区数据来源于中国牧区名录(2021版) [ 10 ] ,主要依据自然地理条件、降水量、地形地貌和当地的生产方式等因素,将内蒙古、青海、新疆、黑龙江、吉林、辽宁、河北、山西、宁夏、西藏、云南、四川、甘肃13个省份归为牧区/半牧区,其他省份为其他区域(非牧区)。二、质量控制依据《中华人民共和国传染病防治法》《国家传染病报告管理办法》等法规制度,传染病报告实行首诊负责制。接诊医生发现布鲁氏菌病感染病例或疑似病例后,根据诊断结果,按照规定及时填写传染病报告卡进行报告。医疗机构疫情管理人员及时审核传染病报告卡,进行错项、漏项、逻辑错误等检查,进行必要的核实、补充或订正,将审核的传染病报告卡及时录入“中国疾病预防控制信息系统”;县级疾病预防控制机构疫情管理人员每日对辖区内报告或数据交换的传染病信息进行审核,对有疑问的报告信息及时反馈报告单位或向报告人核实。对误报、重报信息应及时删除,确认无误后通过网络完成确认审核。医疗卫生机构发生报告病例诊断变更、已报告病例因该病死亡或填卡错误时,应由该医疗卫生机构及时进行订正报告,并重新填写传染病报告卡或抽取电子传染病报告卡,卡片类别选择订正项,并注明原报告病名。对报告的疑似病例,应及时进行排除或确诊。另外,责任报告单位发现漏报的病例后,应及时进行补报。三、统计学分析1.描述性分析:采用Stata 17.0软件进行统计学分析。2010—2024年各县区的年发病率(/10万)=年新发病例数/人口数×100 000;年均发病率(/10万)=2010—2024年新发病例之和/15年人口数之和×100 000。采用χ 2检验分析不同性别、年龄、职业等人群组间的发病率差异,职业主要包括农民、牧民,畜产品收购、销售和屠宰,家务及待业、工人、干部职员、离退人员、民工和其他类型,其中农民、牧民、畜产品收购、销售和屠宰人群为风险人群。2.时间变化趋势分析:采用Joinpoint 5.3.0.0软件,基于研究期间的年发病率拟合Joinpoint模型,根据贝叶斯信息准则(bayesian information criterion,BIC)选择最佳模型,计算年度变化百分比(annual percent change,APC),评价2010—2024年度牧区/半牧区和其他区域的发病趋势,如果APC>0,说明发病率逐年增加,如果APC[ 6 ] 。3.空间扩散特征分析:本研究采用ArcGIS 10.7进行空间插值分析,在区县尺度上以年为度量单位,将首次出现病例报告的县区定义为布鲁氏菌病新扩散的区域,评估研究期间布鲁氏菌病的空间扩散趋势分布 [ 11 , 12 ] ,绘制形成空间扩散趋势地图;基于空间扩散趋势地图,采用等值线的方法,计算上一年扩散区域到下一年扩散区域的最小距离,形成研究期间逐年的扩散距离分布矩阵,并绘制扩散速度图;根据实际扩散速度分布,从含置信区间的线性拟合和含置信区间的分数多项式的估计方法中,选择了 R 2较大的估计方法作为最终的估计方法,计算拟合后的扩散速度图,并展示了 M( Q 1, Q 3)和95% CI值。以上分析均为双侧检验,检验水准α=0.05。结果一、布鲁氏菌病发病情况及变化趋势2010—2024年中国牧区/半牧区共报告252 094例布鲁氏菌病病例,年均发病率为(36.57/100 000±7.28/100 000);其他区域累计报告病例519 748例,年均发病率均数为(2.54/100 000±0.74/100 000),见 图1A 、 1B 。Joinpoint回归模型分析结果显示,牧区/半牧区流行情况呈现下降-反弹-下降的趋势:2010—2017年为发病率下降阶段(APC=-7.20;95% CI:-10.84,-4.80; PCI:13.30,28.56; P=0.015);2021—2024年发病率再次呈现下降趋势(APC=-7.53;95% CI:-13.51,-0.91; P=0.027)( 图1C )。其他区域发病率的流行趋势与牧区/半牧区存在差异,呈现波动上升的趋势:2010—2015年为上升阶段(APC=20.37;95% CI:15.61,27.39; PCI:-27.13,-14.51; PCI:18.70,36.24; PCI:-7.45,5.73; P=0.735),仍保持较高的发病水平( 图1D )。全国农村发病率高于城市(χ 2=318 674, P<0.001);均呈现疑似季节分布趋势,即每年的3—8月为流行高发季节,5月份达到发病最高点( 图1E 、 1F )。图1 2010—2024年中国牧区/半牧区和其他区域人布鲁氏菌病的时间分布特征 图A:牧区/半牧区月发病人数和年发病率分布;图B:其他区域月发病人数和年发病率分布;图C:牧区/半牧区流行变化特征;图D:其他区域流行变化特征;图E:牧区/半牧区城市、农村发病的季节分布;图F:其他区域城市、农村发病的季节分布二、布鲁氏菌病的空间分布特征及扩散趋势2010—2024年牧区/半牧区所有区县年均发病率 M( Q 1, Q 3)为14.51/100 000(0.26/100 000,50.52/100 000),其他区域所有区县年均发病率的 M( Q 1, Q 3)为0.48/100 000(0.12/100 000,3.03/100 000)。年均发病率地图显示,发病率较高的区域主要集中在中国北部省份的牧区/半牧区(内蒙古、新疆、甘肃、青海、宁夏)及其邻近省份,如山西、河北、黑龙江、辽宁、吉林西部、山东、河南北部和中部、陕西北部和中部;此外,云南东部、广西东北部存在零星发病率较高的区县。新发地域的空间扩散趋势分析结果显示,2010—2024年布鲁氏菌病呈现从中国北方疫区向周边邻近区域逐渐扩散和南方地区跳跃性扩散的总体趋势。疫情从牧区/半牧区持续向其他区域扩散,特别是向南方邻近的非牧区扩散,且扩散速度较快,南方地域的跳跃性扩散速度低于北方的扩散速度。其中,2010—2012年疫区主要集中在内蒙古的东部及其中部区域、黑龙江的东部、新疆的北部地区,以及陕西、山西北部临接内蒙古的地域;2012年后,疫情快速地从上述疫区不断向周边扩散,特别是到2018年左右,我国南方绝大部分地区均受到布鲁氏菌病的影响。拟合的逐年扩散速度趋势图显示,2010—2024年疫情从早期的40 km/年的速度持续扩散;2011—2020年扩散速度有所下降;2021年以后,扩散的速度降至最低,平均扩散速度在4.59 km/年;2022年以后,扩散速度逐步增大,2024年平均扩散速度为14.38 km/年。见 图2 。图2 2010—2024年中国人布鲁氏菌病的扩散速度三、性别和年龄分布特征牧区/半牧区、其他区域的城市和农村性别、年龄组累计发病专率金字塔图显示,不同年龄组的发病专率呈现中间高、两头低的总体发病趋势,城市与农村各年龄组的发病专率分布趋势较为一致,且各年龄组男性发病率均高于女性(均 P表1 。牧区/半牧区中,城市和农村地区的发病率最高年龄组均为40~44岁( 图3A 、 3B ),而35~49岁人群(包括35~39、40~44、45~49岁组)的发病专率高于各其他年龄组(χ²=29 836.72, P图3C 、 3D ),而55~64岁人群(包括55~59、60~64岁组)的发病专率高于各其他年龄组(χ²=112 461.02, P图3 2010—2024年中国不同特征成年人人布鲁氏菌病发病专率(1/100 000)图A:牧区/半牧区的城市;图B:牧区/半牧区的农村;图C:其他区域的城市;图D:其他区域的农村四、职业分布特征牧区/半牧区中,剔除职业分布信息“不详”的1 196例病例后,城市报告病例主要集中在农民(74.3%)、牧民(9.5%)和家务及待业(7.3%);农村病例主要集中在农民(78.8%)、牧民(12.9%)和其他(3.7%)。其他区域中,剔除职业分布信息“不详”的4 213例病例后,城市报告病例主要集中在农民(75.3%)、家务及待业(7.5%)和其他(6.2%);农村病例主要集中在农民(85.4%)、其他(4.6%)和牧民(4.0%)。畜产品收购、销售和屠宰人群病例在牧区/半牧区城市、农村和其他区域的城市、农村分别报告病例分别为112、973、3 297和2 466例。讨论20世纪80、90年代,我国人布鲁氏菌病保持在每年1 000~2 000发病病例的较低水平,2000年后疫情逐步波动上升,2023年发病高达7万多例,达到有报告历史以来的最高峰 [ 13 , 14 ] 。内蒙古、新疆、甘肃、宁夏、山西等北部省份是我国人布鲁氏菌病疫情本地传播流行的主要区域,也是我国牛羊等家畜养殖业的重要区域 [ 15 ] 。本文依据中国牧区名录将全国分为牧区/半牧区和其他区域,其中,牧区/半牧区是我国人布鲁氏菌病疫情主要疫区,该区域人间布鲁氏菌病疫情主要与当地家畜养殖业发展密切相关 [ 16 , 17 ] ,其他区域的疫情可能有当地相对分散、少量的家畜养殖业的影响,更主要可能受从北方牧区/半牧区向其他区域转运交易家畜的影响 [ 13 , 18 ] ,其中除云南等少部分地区出现明显的本土传播疫情外 [ 7 ] ,其他大部分地区尚未见本土传播疫情报道,多数为家畜输入引起的疫情 [ 19 , 20 ] 。本文分析2010—2024年两类区域人布鲁氏菌病疫情的变化情况,结果显示,2010年以来,我国人布鲁氏菌病疫情总体上呈现波动上升趋势,可能主要与我国家畜养殖业的发展、肉类需求增长以及牲畜转运、交易增加有密切关系,有多个研究也均表明疫情上升与该原因有关 [ 13 , 18 , 21 ] 。同时,本文研究结果也显示,牧区/半牧区与其他区域疫情变化态势存在显著的差异。2010—2017年在牧区/半牧区疫情下降的情况下,其他区域仍发生一波高峰疫情;2010—2017年北方牧区/半牧区疫情下降可能主要是因为在《国家布鲁氏菌病防治计划(2016—2020)》颁布后 [ 22 ] ,北方高流行区采取强化家畜免疫措施,畜间疫情得到有效控制,人间疫情随之下降,而南方疫情之所以在北方高流行区下降后仍出现一波疫情高峰,其原因尚需更深入的调查研究。因为从一般规律来看,作为疫情输入源头的牧区/半牧区疫情下降时,其他区域疫情可能也将随之下降,除非另外两个因素发生变化,一是该期间从牧区/半牧区向其他区域转运交易牲畜大量增加,二是其他区域出现一波不依赖于牧区/半牧区的本地传播疫情。尽管目前尚不能确定是哪种因素造成两类区域疫情趋势的差异,但这一现象仍提示需要高度关注其他区域多年的输入疫情是否已经出现本土传播的情况,同时高度重视并加强转运交易家畜时的检疫措施。另外,2018年后牧区/半牧区疫情再次上升并在2021年达到峰值,可能与新型冠状病毒感染疫情期间布鲁氏菌病的防控措施力度受到一定影响有关。有研究表明,近年来我国布鲁氏菌病疫情呈现出空间上不断扩散的趋势 [ 13 , 21 ] 。本文通过计算空间扩散速度分析我国布鲁氏菌病疫情扩散情况,结果显示,近年来我国布鲁氏菌病疫情表现为从中国北方疫区向周边邻近区域逐渐扩散和南方地区跳跃性扩散的总体趋势。有研究也显示出类似结果,施玉静等 [ 23 ] 研究发现我国南方绵羊、山羊和牛的布鲁氏菌病血清学阳性率不断上升,尤其在重庆、贵州和四川等地。这一趋势可能主要与我国牛羊等家畜从北方主要家畜养殖区域向南方区域转运和交易不断增加有关 [ 18 ] ,也可能与新报告发病区域内从事饲养、屠宰、售卖、毛皮加工等职业暴露人群增加,该区域动物检疫能力、家畜预防接种不足、屠宰不规范,职业人员发病意识不足等多种因素有关 [ 24 ] 。另外,气候因素,如温度等也可能影响布鲁氏菌病的传播,而医疗资源的配置,如每千人床位比、医护人员的数量,经济水平等也可能影响该病的时空扩散 [ 25 ] 。另外,本研究发现,2022年以后,扩散速率小幅上升,可能和牛羊消费量的增多有关;通过收集2018—2024年牛羊产量的数据,和传播速度呈正相关( Rs=0.75, P=0.052)。另外,也可能与防控力度有关,2022年以来扩散的96个区县,95个不属于列入《国家布鲁氏菌病防治计划(2016—2020)》明确的“河北、山西、内蒙古、辽宁、吉林、黑龙江、陕西、甘肃、青海、宁夏、新疆等11个省份和新疆生产建设兵团达到并维持控制标准”的区域 [ 22 ] 。另外,布鲁氏菌病疫情较快的扩散速度也与布鲁氏菌病的流行病学特点有关,感染布鲁氏杆菌且能够传播的家畜往往症状较轻,转运和交易时难以发现 [ 26 ] ,使疫情快速扩散成为可能。本文研究结果发现,2010年以来,我国布鲁氏菌病疫情扩散速度逐渐下降,从早期的40 km/年下降至2021年以后的4.59 km/年,可能主要因为疫情经过二三十年的传播扩散,已传播至全国各地大部分地区,覆盖面较广,且已具有了一定的规模,继续快速传播的空间已变小。多数研究认为我国布鲁氏菌病患者年龄主要为40~69岁人群 [ 18 , 27 , 28 ] 。本研究结果显示,牧区/半牧区发病人群的年龄结构和其他区域呈现不同的分布特征,在牧区/半牧区,35~49岁人群发病率较高,而在其他区域,55~64岁人群发病率较高。可能与不同区域的职业暴露人群以及感染方式不同有关 [ 5 , 21 , 23 ] ,但仍需要进一步的调查来分析产生这一差异的真实原因。综上所述,近年来我国布鲁氏菌病疫情已扩散至全国范围,防控形势严峻,北方牧区/半牧区仍是防控重点,但南方非牧区疫情仍然呈现持续空间传播扩散的趋势,且该区域疫情变化趋势与北方有显著差异,需要引起高度关注。建议针对我国南北方不同疫情特点,制订针对性防控目标和策略,同时,卫生、畜牧等部门加强联防联控,强化家畜转运和交易的检疫措施;另外,针对重点人群开展感染风险的主动监测,为疫情的风险评估提供更多的数据支撑。