基于文献的贵州省羊群布鲁氏菌感染情况分析

刘平，杨娟，彭启凤，盖文燕，张毅，孙向东，王幼明

（1.中国动物卫生与流行病学中心，山东青岛 266032；
2.遵义医科大学公共卫生学院，贵州遵义 563000；
3.山东畜牧兽医职业学院，山东潍坊 261061）

布鲁氏菌病（brucellosis，以下简称布病）是由布鲁氏菌属细菌引起的一种人兽共患传染病，主要侵害动物生殖系统，临床特征表现为母羊流产、空怀和公羊睾丸炎[1]，为我国二类动物疫病。

贵州省位于我国西南部，当地生态环境优越，具有发展山地特色养殖的良好条件。近几年，全省按照“调结构、上规模、创品牌”的生态畜牧业发展思路，全力推进标准化生产、规模化发展，使畜牧业取得很大发展[2]，牛羊饲养量显著增加。贵州省是我国布病防控二类地区，原则上不实施免疫。在当前畜群“调运频繁，流通无界”的大环境下，动物布病防控难度加大。有资料[3]显示：贵州省人感染布鲁氏菌事件时有发生，且以山羊养殖职业人群感染为主。贵州省羊群布鲁氏菌感染情况，目前研究结果参差不齐。

Meta 分析是通过对同一主题多个小样本研究结果进行综合分析，以达到增大样本量的目的，从而提高原结果的统计效能，使结果更具有代表性[4]，目前在国际上[5-6]及国内公共卫生领域[7-8]流行病学研究中使用广泛，在兽医领域[9-10]也逐渐被推广应用。为了解贵州省羊群布鲁氏菌感染情况，本研究采用系统综述和Meta 分析方法，对既往已发表文献的感染数据进行汇总分析，以期掌握贵州全省羊群布鲁氏菌感染情况，为下一步制定防控政策，开展技术研究提供依据。

1.1 文献检索

从中国期刊全文数据库（CNKI）、维普数据库（VIP）、万方数据库和PubMed数据库，分别以“贵州”“布病”“流行”等为检索词，在主题词项目下检索发表的研究贵州省羊群布病流行率的文献，时间范围从建库至2021 年12 月。

1.2 文献纳入标准

将符合下列条件的文献纳入：①研究区域在贵州省内；
②研究群体包含羊群；
③检测的疫病含布病；
④布病检测方法为虎红平板与试管凝集垂直试验；
⑤检测样品数大于50 份，且结果有明确数值，或通过计算可得到；
⑥当相同数据重复发表时，选择信息全面的文献。

1.3 提取信息

提取文献的题目、作者、发表刊物、发表时间、调查时间、调查地区、羊群类别、羊群数量、阳性羊群数、采集样品数、阳性样品数等信息。

1.4 文献质量评价

依据Grading of Recommendations Assessment，Development and Evaluation（GRADE）评估体系[8]，根据文献是否有明确的调查目的、检测方法、调查时间及亚组分类等信息，进行评价打分。分值范围为0～4 分，其中3～4 分为文献质量较高，2 分为质量中等，0～1 分为质量较差。

1.5 统计学方法

采用Excel 软件提取数据进行初步整理，用Stata12.0 软件进行Meta 分析。利用卡方检验和I2值确定研究间异质性，当I2≥50%时，采用随机效应模型合并效应量，否则采用固定效应模型合并。对研究地市、时间段和羊群类别进行亚组分析，分别估计感染率。采用漏斗图、Egger 和Begg 检验评价纳入文献的潜在发表偏倚。所有检验中，P＜0.05 为差异有统计学意义。

2.1 文献筛选

从中英文数据库检索后共收集到181 篇文献，利用endnote 软件删除重复文献32 篇。依据文献纳入标准，通过浏览标题、摘要和阅读全文，最后共纳入24 篇中文文献进行定量分析。文献纳入路径详见图1。

2.2 研究资料基本情况

纳入的24 篇文献，涉及贵阳市、毕节市、黔南州、黔东南州、黔西南州、铜仁市、六盘水市等7 个地市，无遵义市和安顺市信息，共包含506 843 份样品。文献详细信息见表1。

图1 纳入文献路径

2.3 Meta 分析

2.3.1 各地市个体阳性率 24 篇文献的数据分析结果（图2、表2）显示：贵州省羊群布病平均个体阳性率为2.0%（95%CI：1.5%～2.6%）。黔东南州个体阳性率最高，为4.7%（95%CI：2.4%～6.9%）；
其次是毕节市和铜仁市阳性率为2.0%～3.0%，黔西南州和六盘水市为1.0%～2.0%；
黔南州和贵阳市阳性率较低，分别为0.3%和0.2%。值得注意的是，其中3 篇文献在贵州省9个地市开展，研究结果显示，个体阳性率为3.3%（95%CI：2.7%～3.9%），高于本研究合并后的平均结果。

表1 纳入文献的详细信息

2.3.2 时间和群间分析 由表2 可知：以5 年为1 个阶段，2001—2005 年布病阳性率最高，为2.2%（95%CI：0～6.2%）；
2006 年以来，个体阳性率略有下降，但2016—2020 年仍维持在1.6%。8篇文献[19，21，23-25，29-30，34]详细区分种畜场、商品场和散养户的感染情况。汇总分析发现，商品场的个体阳性率最高，为2.5%（95%CI：1.7%～3.3%），其次是散养户，为1.7%（95%CI：1.0%～2.3%），种畜场个体阳性率最低，为1.4%（95%CI：0.6%～2.2%）。

图2 贵州省羊群布病流行情况森林图

2.3.3 群阳性率 有5 篇文献[22，27，29，31，33]提供群体水平羊群布鲁氏菌感染情况，除1 篇2016 年的研究文献外，其他均是近3 年开展的，涉及黔西南州、贵阳市、铜仁市等地市。汇总分析结果（图3）显示，贵州省羊布病群阳性率为6.0%（95%CI：3.0%～9.0%）。

图3 贵州省羊群布病群体水平感染情况森林图

表2 贵州省羊群布病个体水平三间分布情况

2.4 发表偏倚检验

本研究漏斗图（图4）未发现有明显的不对称，进行Begg 秩相关检验（Z=0.94，＜1.96，P=0.346，＞0.05）和Egger 回归分析（t=1.59，P=0.126）均未提示存在发表偏倚。

图4 贵州省羊群布病阳性率漏斗图

近年来随着系统综述和Meta 分析方法在多个领域的推广，应用此法分析动物疫病流行情况的研究逐渐增多，如有研究[7，35-36]聚焦在我国山羊、奶牛、牦牛布鲁氏菌感染情况。本研究只关注贵州省羊群布鲁氏菌感染情况，但方法学与之类似。

有研究[35]对我国2000—2018 年山羊和绵羊的布鲁氏菌感染情况做过Meta 分析，结果发现羊群布病个体阳性率由2000—2009 年的1.0%上升至2010—2018 年的3.2%，其中对贵州省的分析显示，共检测16 639 份样品，阳性率为5.8%，高于本研究结果（2.0%）。本研究包含506 843 份样品，时间跨度为2001—2020 年，相比上个研究，分析结果可能更切近实际。

综合贵州省本地调研数据，当地羊群布病平均群体阳性率为6.0%，平均个体阳性率为2.0%，不同地市略有差异。研究显示，2017—2019 年，云南省昆明市羊布病场群和个体阳性率平均为2.51%和1.31%[37]，2020 年四川省甘洛县羊布病场群和个体阳性率分别为0.047 3%和0.015 8%[38]，可见与西南地区其他省份相比，贵州省羊群布病感染率较高，羊群布病防控形势较严峻。从2016 年起，贵州省已开始实施家畜布病监测和净化工作，对阳性家畜实施“检、消、杀、烧、埋”综合防控措施，到2020 年防控工作取得显著成效[39]，这从本研究个体流行率的时间变化趋势也可体现。

布病等人兽共患病防治工作事关畜牧业高质量发展和人民群众身体健康，事关公共卫生安全和国家生物安全。2022 年9 月份，农业农村部印发《全国畜间人兽共患病防治规划（2022—2030 年）》，其中明确提出布病防控要继续坚持免疫与净化相结合的综合防治措施。在2022 年印发的全国布鲁氏菌病免疫县和免疫奶牛场名单中，未见有贵州省上报的免疫县，推测该省布病防控仍延续之前的监测净化策略。布病净化是艰辛、复杂和漫长的工作，可向走在前列的地区借鉴经验，做好重点畜群（引进畜、种公畜、繁殖母畜）监测，管理好种畜场和阳性场，重点抓住配种期、分娩期和出栏期等，严格落实动物移动监管，严格监测扑杀和处置灭源[40]。也可响应国家政策号召，在确保生物安全的前提下，积极尝试开展阳性动物集中无害化处理和资源化利用试点工作。

本研究纳入24 篇文献，定量分析贵州省羊群布病流行情况，且无明显发表偏倚，结果较为可信。但仍存在一些局限性：一是多数文献对采样动物来源群体没有详细描述，不清楚羊只是来自健康羊群，还是有流产症状羊群，是来自种畜场还是散养户等；
二是文献对羊只没有详细描述，是本地羊还是外调羊，羊只的月龄、性别等也不清楚。这些因素的差异均会影响布病流行率，也是导致本研究数据合并分析异质性较强的主要原因。此外，本研究只收集了上述数据库中的已发表文章，对于在其他数据库发表的文献或未发表文献无法收集相关数据，这对结果可能会造成一定影响。

本研究通过定量分析发现，贵州省羊群布病防控情况不容乐观。在今后的工作中，要强化羊群布病监测和阳性动物扑杀，同时对发现阳性动物进行紧急调查，做好职业人群的宣传和预防，严格做好检疫和调运监管。通过开展“检疫、监测、消毒、扑杀、宣传”的综合性防控，稳步推进布病净化工作。