牛支原体（Mycoplasma bovis）是全球牛群中一种重要的病原体；然而，关于中国国家级的综合流行病学证据仍然有限。为了填补这一知识空白，研究人员首次对2000年至2025年间发表的61项研究进行了全国性综合，涵盖了50,325头牛，旨在描述患病率的变异程度及其时空分布，以及中国大陆牛群中M. bovis感染的潜在流行病学相关性。遵循PRISMA指南，研究人员对61项研究进行了系统评价和荟萃分析（meta-analysis），并利用随机效应模型（random-effects models）汇集了动物水平的患病率估计值，同时应用亚组分析、荟萃回归、空间映射和时间序列预测来量化和解释不同研究、地区和时间之间的异质性。各研究间的患病率估计值差异很大，动物水平的全国汇总估计值为33.0%（95% CI: 29.0–37.0%），表明存在显著但高度异质性的流行病学负担。尽管显著的跨研究异质性和潜在的“小样本效应”需要谨慎解释，但敏感性分析得出的估计值一致（约32–33%）。时间分析表明，虽然存在明显的但适度的年度下降（-1.5个百分点），但在不同的研究时期和数据来源之间存在变异性。患病率在不同的流行病学背景下存在显著差异，其中牦牛和奶牛、犊牛和青年牛以及表现出牛呼吸道疾病（BRD）的动物中观察到的估计值较高。时空映射进一步突显了显著的地理异质性，确定了华中（湖北、河南）、西南（重庆、四川）、西北（新疆、青海）以及华北/东北地区（黑龙江、内蒙古）持续存在的高患病率集群。探索性ARIMA预测显示，2026–2028年全国患病率将大致稳定在约27%，这与明显的时间模式一致，但也伴随着反映数据稀疏和模型不确定性的宽置信区间。区域预测也显示出实质性的差异，华东地区保持高位（约60%），华中地区相对稳定（约32%），而东北地区预计到2028年将出现急剧下降，尽管由于数据有限，这一估计应谨慎解释。始终存在高度的异质性——跨越研究、宿主特征、地区和时间——表明单一的汇总估计值可能会掩盖重要的流行病学差异。这些发现表明，M. bovis感染在中国仍然广泛存在且异质性强。加强协调的监测系统以及对高风险宿主和年龄组的针对性监测，对于改善疾病检测和制定基于证据的、针对特定地区和环境的控制策略至关重要。
研究背景与意义
牛支原体（Mycoplasma bovis）是引起牛呼吸道疾病（BRD）、乳腺炎、关节炎及生殖障碍的全球性重要病原体，其持续感染、易形成亚临床感染以及对常用抗菌药物的低敏感性，使得根除极为困难，给全球养牛业造成了巨大的经济损失。中国作为世界上最大的乳制品和牛肉生产国之一，密集的生产系统和频繁的牲畜流动加剧了该病原体的传播负担。尽管该病原体自2008年在中国首次分离以来已在多个地区被检测到，但由于现有调查受限于诊断方法异质性、地理覆盖不均及时间代表性不足，导致中国大陆牛群中该疾病的国家层面流行病学特征仍不明确。为了填补这一空白，明确其时空分布规律及风险因素，研究人员开展了这项涵盖2000年至2025年的全国性综合研究，该研究成果发表在国际知名期刊《Journal of Dairy Science》上。该研究旨在通过整合多源数据，揭示流行规律，为制定针对性的疾病控制策略和监测体系提供科学依据。
关键技术方法
研究人员采用了PRISMA指南指导的系统评价与荟萃分析框架。首先，通过检索多个中英文数据库，依据严格的纳入排除标准筛选出61项符合条件的研究，共计包含50,325头样本。在数据处理阶段，研究人员应用了随机效应模型进行患病率汇总，并利用逆方差加权法计算总体估计值。为了量化异质性，研究采用了Cochran's Q统计量、τ²及I²指标。在空间分析方面，结合地理信息系统（GIS）技术，将省级数据映射到Albers等面积圆锥投影上，并按时间间隔绘制choropleth地图。在时间序列分析中，研究人员采用自回归积分滑动平均模型（ARIMA）对2026至2028年的流行趋势进行了预测。此外，研究还运用了加权平均最小二乘（WLS）元回归方法，评估了年份和地理区域对患病率的联合影响。
研究结果
研究筛选过程
通过系统检索，最初识别出1,130条记录。经过去重、标题和摘要筛选、全文评估以及质量评价后，最终共有61项研究被纳入定量分析。这些研究主要集中在西北地区（如新疆、宁夏），其次是西南地区，而东北和中部地区的文献相对较少。
研究特征与质量评估
纳入的研究质量普遍较高，79%的研究被评为高质量。时间分布上，大部分研究集中在2015年至2024年间。这表明近年来国内对该领域的关注度显著提升，但地理分布的不均衡也为后续的区域对比分析带来了挑战。
汇总患病率与发表偏倚
动物水平的全国汇总患病率估计为33.0%（95% CI: 29.0–37.0%）。尽管存在极高的异质性（I² = 99.2%）和一定程度的发表偏倚（小样本研究倾向于报告更高患病率），但通过留一法敏感性分析验证了结果的稳健性，排除任何单一研究均未对整体估计值造成实质性改变。
通过亚组分析识别并量化患病率荟萃分析中的异质性来源
通过对不同维度的亚组分析，研究人员明确了导致总体异质性显著的关键来源。
• 地理变异 ：患病率存在显著的地理差异。粗患病率以华中地区最高（47.6%），东北地区最低（15.0%）。经模型调整后，华中地区的汇总患病率高达77.5%，而华南和华北地区相对较低。
• 省份水平变异 ：各省之间存在巨大差异。湖北省报告的粗患病率最高（96.5%），山东省最低（5.9%）。新疆和宁夏虽然粗患病率中等，但由于研究数量较多，其汇总估计值在30%至35%之间。
• 动物类型差异 ：牦牛的粗患病率最高（41.6%），其次为奶牛（33.4%）和肉牛（24.4%）。调整后的汇总估计值分别为牦牛40.0%、奶牛33.2%和肉牛23.0%。
• 年龄组效应 ：犊牛（≤1岁）的感染率显著高于成年牛（≥3岁）。犊牛粗患病率为30.3%，汇总估计达40.1%；而成年牛粗患病率为28.3%，汇总估计降至23.4%。
• 临床症状的影响 ：表现出牛呼吸道疾病的牛只患病率明显高于无特定临床症状的群体。BRD病例的粗患病率为45.0%，汇总估计为48.4%；而无特定症状群体的粗患病率为26.3%，汇总估计为27.5%。
• 样本类型的影响 ：关节液检测出的粗患病率最高（46.2%），其次为鼻拭子（33.5%）、肺组织（32.0%）、血液（28.2%）和牛奶（29.5%）。经模型调整后，各样本类型的汇总估计值与粗患病率趋势基本一致。
• 研究持续时间的影响 ：随着研究持续时间的延长，观察到的患病率呈上升趋势。短期研究（≤1年）粗患病率为23.8%，中期（1-3年）升至36.5%，长期（>3年）则达到47.1%。
• 诊断方法比较 ：分子学诊断方法检测出的粗患病率为31.8%，血清学方法为26.0%。两者在统计学上没有显著差异，表明诊断方法的差异并非导致研究间异质性的唯一因素。
中国牛支原体患病率的时态趋势
从2008年至2025年，全国患病率呈现出明显的年度波动，总体呈现下降趋势。2008年的估计值超过80%，期间经历数次起伏，2021年后稳定在17%至31%之间。WLS元回归分析证实了每年约1.5个百分点的显著下降。然而，区域趋势各异：西北地区呈显著上升态势；华北地区波动剧烈，无明显单调趋势；华东地区极不稳定；西南和华南地区虽有波动但总体处于稳定流行状态；而东北和华中地区则表现出显著的下降趋势。
M. bovis患病率的预测趋势
基于ARIMA模型的预测显示，2026年至2028年全国平均患病率约为27%，但伴随较宽的预测区间。区域预测差异显著：华东地区预计维持在高位（约60%），华中地区相对稳定（约32%），东北地区预计将持续下降至接近零；其他如西南、华北、华南和西北地区则预计维持在20%至30%的中等水平。
中国牛群牛支原体感染的空间分布
高患病率区域呈现出明显的时空聚集性。早期（2008-2010年）主要在黑龙江和青海；随后（2011-2019年）高负担区域转移至新疆、四川、重庆及湖北等地；近期（2023-2025年）热点区域演变为河南、内蒙古和新疆。这种空间转移反映了局部传播放大效应和疾病负担的地理再分布。
讨论与结论
讨论总结
讨论部分指出，这项研究首次全面揭示了中国大陆牛群中M. bovis感染的复杂流行病学图景。研究发现的高异质性不仅源于生物学因素（如宿主类型、年龄、临床症状），也深受方法论差异（如样本基质、诊断窗口期、研究设计）的影响。例如，呼吸道样本和组织样本因病原体载量高而检出率高，血清学抗体滴度则受免疫动态影响。此外，研究强调了时空模式的动态变化，表明国家层面的总体趋势无法替代地方层面的具体驱动机制。尽管预测模型受限于数据的非连续性，但其结果仍为识别高风险区域提供了有价值的线索。
结论翻译
这项首例全国性综合分析表明，M. bovis对中国牛群构成了沉重、持续且极具异质性的负担，各项研究的患病率估计值差异显著，且在不同流行病学背景下存在广泛的变异性（汇总估计：33.0%）。该病原体的流行病学似乎与关键的生产和宿主因素密切相关，在集约化奶牛养殖系统和牦牛中患病率显著较高，且犊牛和青年牛表现出明显的年龄特异性易感性。尽管出现了明显的适度全国性下降，但基于模型的预测表明，由于短期预测指示到2028年仍将维持约27%的持续负担，且存在显著的区域变异性及潜在的分歧轨迹，因此自发消除的可能性不大。观察到的患病率模式异质性，连同M. bovis与BRD的关联，以及在纵向研究和呼吸道组织采样中报告的较高检出率，均进一步证实患病率估计值高度依赖于具体环境，无法直接在不同研究间进行比较。这凸显了统一监测和防控方法的局限性。因此，研究人员强调迫切需要制定针对特定区域的、基于风险的防控策略。优先行动应包括：（1）在高负担地区加强具有代表性的、按年龄分层的纵向监测；（2）在各实验室间标准化多基质诊断方案和数据报告；以及（3）实施有针对性的生物安全措施和流动控制，特别关注奶牛养殖系统、幼龄牛群及BRD管理。总体而言，始终存在的高程度异质性表明，依赖单一的汇总估计值可能会掩盖重要的流行病学差异，防控工作应依据针对特定地区和环境的证据来指导。将这些综合性的流行病学见解转化为实际行动，对于实现减轻M. bovis负担的目标具有重要意义。
打赏
