遼寧省某市奶牛布病感染狀況及群間傳播風(fēng)險因素調(diào)查

遼寧省某市奶牛布病感染狀況及群間傳播風(fēng)險因素調(diào)查

崔基賢1，王靖飛2，谷志大1

（1.遼寧省動物疫病預(yù)防控制中心/遼寧省動物醫(yī)學(xué)研究院，遼寧沈陽 110164；2.中國農(nóng)業(yè)科學(xué)院哈爾濱獸醫(yī)研究所獸醫(yī)生物技術(shù)國家重點(diǎn)實(shí)驗(yàn)室/疫病診斷中心，黑龍江哈爾濱 150001）

為了解遼寧省某市的奶牛布魯氏菌?。ê喎Q奶牛布病）流行現(xiàn)狀及該市奶牛布病群間傳播風(fēng)險因素，開展本研究。對該市轄區(qū)內(nèi)所有奶牛進(jìn)行布魯氏菌病的普查，獲知該市奶牛布病的表觀個體血清流行率和表觀群血清流行率；采用問卷調(diào)查方式進(jìn)一步開展病例-對照研究，初步篩選出14項(xiàng)該市奶牛布病群間傳播可能風(fēng)險因素，經(jīng)相關(guān)性檢驗(yàn)結(jié)合生物學(xué)知識，選取了5個變量進(jìn)入多因素分析，最終通過非條件logistic回歸方法初步建立了該市奶牛布病群間傳播風(fēng)險因素回歸模型，并建立ROC曲線計(jì)算了模型的預(yù)測概率。研究結(jié)果顯示，該市奶牛個體流行率為0.99%，群流行率12.44%。非條件logistic回歸分析結(jié)果顯示，該市奶牛布病群間傳播的風(fēng)險因素包括：“無專用購銷車輛”（OR=3.57，95%CI: 1.04-12.5，p =0.04），“進(jìn)圈時不穿專用工作服和膠靴”（OR=5.0，95%CI: 0.97-25.00 ，p =0.06），“與犬混養(yǎng)”（OR=12.39，95%CI: 0.97-157.56，p =0.05）。該結(jié)果表明與犬混養(yǎng)可能是該市奶牛布病群間傳播重要的風(fēng)險因素，而奶牛銷售時自有購銷專用車輛和場內(nèi)生物安全措施的采取如進(jìn)圈穿專用工作服可使奶牛布病在場群間傳播的風(fēng)險降低。針對回歸模型建立的ROC曲線下面積為71.6%，表明此模型解釋該市奶牛布病群間傳播風(fēng)險的能力為中等。

奶牛；布魯氏菌?。涣餍新?；風(fēng)險因素

奶牛布病是由布魯氏菌引起的以奶牛流產(chǎn)和繁殖障礙為主要特征的一種慢性人畜共患傳染病，嚴(yán)重危害家畜和人類健康，在世界許多國家和地區(qū)都曾造成巨大的經(jīng)濟(jì)損失，以美國為例，1934—1937年根除布病費(fèi)用達(dá)35 億美元，然而僅在1952年1年由布病引起奶牛產(chǎn)奶量下降和流產(chǎn)而導(dǎo)致的損失就達(dá)400億美元[1]。并對人類公共衛(wèi)生安全產(chǎn)生巨大的威脅。該病在全球范圍內(nèi)廣泛流行，目前僅有北歐、中歐、加拿大、日本、澳大利亞和新西蘭被承認(rèn)已根除布魯氏菌病[2]。在我國，《國家中長期動物疫病防控規(guī)劃（2012—2020）》規(guī)定布魯氏菌病為重點(diǎn)防控的16種動物疫病之一。并且要求到2015年，遼寧省奶牛布病達(dá)到穩(wěn)定控制標(biāo)準(zhǔn)。該市是我省奶牛養(yǎng)殖主產(chǎn)區(qū)之一，以現(xiàn)代流行病學(xué)方法對該市奶牛布病感染狀況及風(fēng)險因素進(jìn)行調(diào)查，對制定奶牛布病針對性防控策略，有效推動我省奶牛布病防控工作，具有一定的示范意義。

為了解該市奶牛場布魯氏菌病流行現(xiàn)狀，并提出防控該病的有效策略，開展了針對該市奶牛場的布病感染狀況和風(fēng)險因素調(diào)查。

1 材料和方法

1.1 研究設(shè)計(jì)

1.1.1 場群流行率和個體流行率。遼寧省對奶牛布病的防控策略為“全群檢疫加撲殺”。結(jié)合省奶牛布病凈化工作，2013年對該市轄區(qū)內(nèi)所有奶牛進(jìn)行布魯氏菌病普查，包括187個場群，37888頭奶牛。獲知該市奶牛布病的個體表觀流行率和群表觀流行率。查閱文獻(xiàn)確定診斷試驗(yàn)敏感性和特異性，通過公式 TP=（AP+Se-1）/（Se+Sp-1）計(jì)算真實(shí)個體流行率。

1.1.2 檢測方法。虎紅平板試驗(yàn)和試管凝集試驗(yàn)的垂直試驗(yàn)。查閱文獻(xiàn)得知該垂直試驗(yàn)敏感性81%，特異性99.9%[3]。

1.1.3 風(fēng)險因素調(diào)查。以場群為單位采用問卷調(diào)查方式開展病例-對照研究。陽性奶牛定義為虎紅平板試驗(yàn)和試管凝集試驗(yàn)的垂直試驗(yàn)結(jié)果為陽性的奶牛；陽性群定義為群內(nèi)至少有一頭陽性奶牛的場或散養(yǎng)戶。陰性群定義為場內(nèi)未檢測到奶牛布病的場或散養(yǎng)戶。所有陽性群均列為調(diào)查對象，并在該市轄區(qū)內(nèi)遵循同一區(qū)域且養(yǎng)殖條件相似的匹配原則選擇3倍數(shù)量于陽性群的陰性群作為對照群進(jìn)行調(diào)查。

1.2 數(shù)據(jù)收集和整理

設(shè)計(jì)了《奶牛布病風(fēng)險因素調(diào)查問卷》， 分別針對飼養(yǎng)管理、人員來訪、與其他動物混養(yǎng)情況、產(chǎn)后胎盤和流產(chǎn)胎兒處理情況、場內(nèi)生物安全措施、奶牛銷售、擠奶情況等7個方面的23項(xiàng)風(fēng)險因素展開調(diào)查，變量均設(shè)計(jì)為二分類變量。收集的數(shù)據(jù)錄入Microsoft excel 2010。

1.3 數(shù)據(jù)分析

所有變量均為二分類變量，以場群布病狀態(tài)為因變量，采用Fisher確切概率檢驗(yàn)對23個協(xié)變量（風(fēng)險因素）進(jìn)行單變量分析。篩選出P<0.2的變量，進(jìn)行多變量logistic回歸模擬，建立非條件logistic回歸分析模型。對于篩選出的變量首先進(jìn)行雙變量相關(guān)性檢驗(yàn)，若有變量間Pearson相關(guān)系數(shù)>0.4，在其中選一個最符合生物學(xué)常識的進(jìn)入logistic回歸模型[3]，將另一個剔除。執(zhí)行回歸操作時采取“全部進(jìn)入逐個排除”的方式，以p<0.05作為檢驗(yàn)水準(zhǔn)。并對最終保留的變量進(jìn)行交互作用檢驗(yàn)[4]。采用Hosmer-Lemeshow法[5]對建立的模型進(jìn)行擬合優(yōu)度檢驗(yàn)。P>0.05說明預(yù)測概率獲得的期望頻數(shù)與預(yù)期頻數(shù)無統(tǒng)計(jì)學(xué)差異，模型擬合較好。對最終建立的logistic模型制作ROC曲線，檢驗(yàn)?zāi)Ｐ皖A(yù)測概率的擬合優(yōu)度。以上統(tǒng)計(jì)學(xué)分析均應(yīng)用SPSS 20.0和EpiInfoTM 7完成。

2 結(jié)果

2.1 群流行率和個體流行率

對該市轄區(qū)所有奶牛進(jìn)行普查，包括187個場群，37888頭奶牛，共檢測到19個陽性場，67頭陽性奶牛，群表觀流行率10.16%，個體表觀流行率0.18%。根據(jù)公式TP=（AP+Sp-1）/（Se+Sp-1），計(jì)算得出真實(shí)個體流行率為0.99%，真實(shí)群流行率12.44%。

2.2 風(fēng)險因素調(diào)查結(jié)果

下發(fā)至19個陽性場和57個陰性場共計(jì)76份調(diào)查問卷全部回收，回收率100%。根據(jù)調(diào)查問卷結(jié)果，對單因素在感染場和陰性場間的分布做了分析（表2），篩選出單因素分析中p<0.2的14項(xiàng)因素，對這14項(xiàng)變量進(jìn)行相關(guān)性檢驗(yàn)，結(jié)合生物學(xué)知識，篩選出“沒有購銷專用車輛”“與犬混養(yǎng)”“進(jìn)圈時不穿工作服及膠靴”等3個變量進(jìn)入最終的logistic回歸模型，分析結(jié)果發(fā)現(xiàn) “無專用購銷車輛”（OR=3.57，95%CI: 1.04-12.5，p =0.04），“進(jìn)圈時不穿專用工作服和膠靴”（OR=5.0，95%CI: 0.97-25.00 ，p =0.06），“與犬混養(yǎng)”（OR=12.39，95%CI: 0.97-157.56，p =0.05）是該市牛場感染布魯氏菌的風(fēng)險因素。針對logistic回歸模型建立了ROC曲線，曲線下面積為71.6%，表明此模型預(yù)測該市奶牛布病群間傳播風(fēng)險能力為中等[6]。

表1 X市奶牛布病群間傳播風(fēng)險因素單變量分析結(jié)果

3 討論

本研究通過對該市奶牛布病的全群普查，掌握了該市奶牛布病的場群流行率和個體流行率，并以問卷調(diào)查的方式開展了一次病例—對照研究。

由于本研究樣本量較小，故通過Fisher確切概率檢驗(yàn)對調(diào)查獲得的風(fēng)險因素?cái)?shù)據(jù)逐一進(jìn)行單因素分析，確定了對該市奶牛布病群間傳播有統(tǒng)計(jì)學(xué)意義的6項(xiàng)因素。但在建立logistic回歸模型時，放寬了進(jìn)入模型變量的選擇條件，準(zhǔn)入條件設(shè)為Fisher 確切概率檢驗(yàn)p<0.2。對符合條件的14項(xiàng)變量進(jìn)行了相關(guān)性檢驗(yàn)，結(jié)合生物學(xué)知識，最終選擇了“自繁自養(yǎng)”“收拾棚圈時穿專用工作服和膠靴”“與犬混養(yǎng)”“上門收購”“無專用購銷車輛”等5個變量建立logistic回歸模型，逐個剔除變量，觀察其他因素OR值變化，同時檢驗(yàn)?zāi)Ｐ妥兓闆r，最終的logistic回歸模型包括“進(jìn)圈時不穿專用工作服和膠靴”“與犬混養(yǎng)”“無專用購銷車輛”等3個變量。針對該模型制作的ROC曲線下面積為71.6%，表明此模型預(yù)測該市奶牛布病群間傳播的能力為中等[5]，該市存在其他影響奶牛群間布病傳播風(fēng)險因素，有待進(jìn)一步研究探索。

表2 X市奶牛布病群間傳播風(fēng)險因素的多變量logistic回歸分析結(jié)果

圖1 非條件logistic回歸模型的ROC曲線

由于本次病例-對照研究是以問卷調(diào)查方式進(jìn)行，而布魯氏菌病相關(guān)防控知識在奶牛飼養(yǎng)場戶中已廣泛普及，在調(diào)查時養(yǎng)殖場戶對某些問題的回答可能存在偏倚，例如流產(chǎn)胎兒的處理方式，回答隨意丟棄和烹飪食用的比例可能比實(shí)際偏低。在以后的研究中還需進(jìn)一步探討如何設(shè)計(jì)更特異的問題，盡可能降低偏倚。

4 結(jié)論

4.1 采用虎紅平板試驗(yàn)和試管凝集試驗(yàn)的垂直試驗(yàn)對該市所有奶牛場中的全部奶牛進(jìn)行布魯氏菌病感染狀況普查，結(jié)果顯示該市2013年奶牛布病個體流行率為0.99%，場群流行率為12.44%。

4.2 進(jìn)行了針對該市奶牛布病群間傳播風(fēng)險因素的病例-對照研究，結(jié)果表明與犬混養(yǎng)可能是該市奶牛布病群間傳播重要風(fēng)險因素，而奶牛銷售時自有購銷專用車輛和提高場內(nèi)生物安全措施如進(jìn)圈穿專用工作服可降低奶牛場群間布病傳播風(fēng)險。

[1]CDC.Brucellosis.2009，http：//www.cdc.gov/brucellosis/ clinicians/brucella-species.html.

[2]http：//www.oie.int/wahis_2/public/wahid.php/Diseaseinformation/Diseasedistributionmap.

[3] Hegazy Y M，Moawad A，Osman S，et al. Ruminant Brucellosis in the Kafr El Sheikh Governorate of the Nile Delta，Egypt：Prevalence of a Neglected Zoonosis[J]. Public Library of Science，2011，5（1）：1-9.

[4]Dohoo IR，Martin SW，Stryhn H. Veterinary Epidemiologic Research[M]. Charlottetown，Canada：AVC Inc，2003.

[5]Hosmer D W，Lemeshow S. Applied Logistic Regression[M]. New York：Wiley J and Sons，2000.

[6]宇傳華. SPSS與統(tǒng)計(jì)分析[M].北京：電子工業(yè)出版社，2007.

Investigation of Prevalence and Risk Factors of Brucellosis in Dairy Cows in X City of Liaoning Province

Cui Jixian1，Wang Jingfei2，Gu Zhida1

（1. Liaoning Center for Animal Disease Prevention and Control/Liaoning Veterinary Research Institute；2. Center for Animal Disease Diagnosis/State Key Laboratory of Biotechnology，Harbin Veterinary Research Institute，Chinese Academy of Agricultural Sciences，Harbin 150001）

To estimate the prevalence and identify risk factors of brucellosis in dairy cows in X city of Liaoning province，all the cows in X city were tested using both rose bengal plate test and tube agglutination test，and the risk factors for Brucella infection on farms were identifed through a case-control study. The results showed that the apparent individual sero-prevalence of brucellosis in cows in X city was 0.18%，while the apparent herd sero-prevalence of brucellosis was 10.16%. Three factors，including presence of dogs in the farm（OR=12.39，95%CI: 0.97-157.56，p =0.05），no specialized vehicles（OR=3.57，95%CI: 1.04-12.5，p =0.04）and no work suits and boots during working（OR=5.0，95%CI: 0.97-25.00，p =0.06）were identifed as the risk factors for brucellosis transmission between farms in X city. It was suggested that brucellosis control measures should be designed according to the results of the investigation in X city.

dairy cow；brucellosis；prevalence；risk factors

S852.23

：B

：1005-944X（2014）11-0051-04