劉愛(ài)玲,馬清霞,韓 雪,沈朝建,康京麗,張 毅,王幼明
(1. 中國(guó)動(dòng)物衛(wèi)生與流行病學(xué)中心,山東青島 266032;2. 青島市動(dòng)物疫病預(yù)防與控制中心,山東青島 266071;3. 北京維通利華實(shí)驗(yàn)動(dòng)物技術(shù)有限公司,北京 100107)
牛結(jié)核?。╞ovine tuberculosis,bTB)是由牛結(jié)核分枝桿菌引起的一種人獸共患慢性傳染病,可侵害多種動(dòng)物,家畜中奶牛最易感染[1]。結(jié)核病病原為分支桿菌屬細(xì)菌,在自然界分布廣泛,對(duì)動(dòng)物有致病性的主要是結(jié)核分支桿菌、牛分支桿菌、禽分支桿菌和副結(jié)核分支桿菌。可引起牛結(jié)核病的主要是牛分支桿菌、結(jié)核分支桿菌、禽分支桿菌,其中牛分支桿菌致病性最強(qiáng)。本屬病原菌均為一種纖細(xì)、平直或稍彎曲的桿菌,革蘭氏染色陽(yáng)性,無(wú)莢膜,不產(chǎn)生內(nèi)外毒素,不形成芽孢,無(wú)運(yùn)動(dòng)性,對(duì)外界環(huán)境抵抗力較強(qiáng),耐干燥和濕冷,但不耐熱[2]。
bTB的潛伏期長(zhǎng)短不一,自然感染病例的潛伏期一般為16~45 d,但有的可在幾個(gè)月以上,甚至數(shù)年。Morris等[3]研究表明,雖然牛結(jié)核分枝桿菌宿主廣泛,但近些年主要感染牛和人。bTB可通過(guò)空氣傳播[4],不容易通過(guò)污染的飼料和牧草傳播,因?yàn)樵诂F(xiàn)實(shí)條件下污染物中感染細(xì)菌的存活時(shí)間相對(duì)較短,動(dòng)物通過(guò)消化道途徑不容易暴露至感染劑量,但可以通過(guò)口咽部黏膜感染[3]。有研究表明,將患牛引進(jìn)到無(wú)疫牛群是bTB傳播的主要風(fēng)險(xiǎn)因素。在動(dòng)物水平方面,bTB傳播的風(fēng)險(xiǎn)因素包括易感動(dòng)物年齡、性別、品種、健康狀況、免疫狀況,以及對(duì)bTB的遺傳抗性和易感性,還有垂直傳播等[5]。
近年來(lái),奶產(chǎn)品銷量迅速上升,奶牛飼養(yǎng)量也明顯提高,因此對(duì)奶牛疾病的預(yù)防不可忽視。尤其是bTB,其不但影響奶牛生產(chǎn)力,而且威脅人類身體健康。因此,加強(qiáng)bTB防控具有重要意義。萊西市是山東省主要奶牛生產(chǎn)基地,擁有“雀巢、迎春樂(lè)、新希望”等多家龍頭奶業(yè)企業(yè)[6]。為評(píng)估該市無(wú)bTB奶牛場(chǎng)(戶)通過(guò)引進(jìn)奶牛感染該病的風(fēng)險(xiǎn),本研究繪制了奶牛引入傳播bTB的風(fēng)險(xiǎn)路徑,建立了定量評(píng)估模型,以期為該市bTB防控策略制定提供技術(shù)支持。
奶牛來(lái)源地區(qū)bTB場(chǎng)群流行率和個(gè)體流行率無(wú)差異。
引進(jìn)奶牛到達(dá)萊西市后,采用單純皮內(nèi)變態(tài)反應(yīng)(SICT)方法,對(duì)所有奶牛進(jìn)行bTB檢測(cè);采用全進(jìn)全出方式隔離30 d,每天對(duì)所有奶牛進(jìn)行臨床癥狀和體溫變化監(jiān)測(cè),隔離后進(jìn)行消毒。
bTB通過(guò)奶牛引進(jìn)傳播的風(fēng)險(xiǎn)路徑見(jiàn)圖1。從存在bTB的養(yǎng)殖場(chǎng)引入奶牛,采用SICT檢測(cè)方法,不能夠發(fā)現(xiàn)所有感染奶牛,引入后存在bTB釋放風(fēng)險(xiǎn);引入的bTB感染奶牛即使未出現(xiàn)癥狀,在隔離階段也存在bTB釋放風(fēng)險(xiǎn);在隔離結(jié)束后混群,動(dòng)物暴露于bTB,存在感染風(fēng)險(xiǎn)。
該風(fēng)險(xiǎn)路徑假設(shè)奶牛能夠持續(xù)感染并傳播疫病,并在引入后釋放感染。本研究對(duì)隔離期間,bTB通過(guò)奶牛直接接觸以及糞便或空氣傳播感染的概率不作探討。
將參數(shù)整合到模型中進(jìn)行概率分布分析,在Microsoft Excel中使用@Risk 軟件,采用Monte Carlo方法將批量為50的模型運(yùn)行10 000次。參數(shù)描述和賦值見(jiàn)表1。參數(shù)估計(jì)來(lái)自萊西市動(dòng)物疫病預(yù)防控制中心調(diào)查數(shù)據(jù)和專家觀點(diǎn)。
使用@Risk軟件,對(duì)影響模型的重要參數(shù)進(jìn)行敏感性分析,確定各參數(shù)的影響程度及相關(guān)性。
根據(jù)下列公式計(jì)算通過(guò)從感染地區(qū)引進(jìn)1頭奶牛傳播bTB的概率為2.93×10-3(95%CI:1.14×10-3~8.25×10-3),概率分布結(jié)果見(jiàn)圖 2。
圖1 奶牛引進(jìn)傳播bTB的風(fēng)險(xiǎn)路徑
根據(jù)下列公式,計(jì)算通過(guò)從感染地區(qū)引進(jìn)50頭奶牛至少有1頭傳播bTB的概率為0.137(95%CI:0.056~0.339),概率分布結(jié)果見(jiàn)圖3。
該模型中的參數(shù)計(jì)算結(jié)果及不確定分析見(jiàn)表2。
模型主要影響參數(shù)的相關(guān)系數(shù)見(jiàn)圖4。結(jié)果顯示,奶牛是否暴露于bTB感染奶牛的概率和群流行率相關(guān)系數(shù)分別為0.63和0.37,對(duì)奶牛引進(jìn)時(shí)傳播bTB的概率影響較大,是主要風(fēng)險(xiǎn)因素。隔離措施相關(guān)系數(shù)分別為-0.52和-0.28,是主要保護(hù)因素。
萊西市是奶牛養(yǎng)殖示范地區(qū),養(yǎng)殖模式分為規(guī)模場(chǎng)、養(yǎng)殖小區(qū)和奶牛集中村3種,規(guī)?;潭缺容^高,生物安全措施較好。本研究數(shù)據(jù)來(lái)自2016年萊西市監(jiān)測(cè)數(shù)據(jù)。該地區(qū)場(chǎng)群流行率為7.3%,感染場(chǎng)內(nèi)個(gè)體流行率為1.9%[8]。引進(jìn)地區(qū)的bTB流行率可能存在差異,因此模型存在一定的不確定性。
在奶牛引進(jìn)時(shí),萊西市采用SICT方法進(jìn)行檢測(cè)。雖然目前SICT仍是bTB檢測(cè)的國(guó)際通用方法,但其敏感性較低并不是最佳方法[9]。IFN-γ檢測(cè)方法也被廣泛應(yīng)于bTB檢測(cè),但是試驗(yàn)的特異性和敏感性也不是最好。診斷試驗(yàn)方法的不足也是造成bTB傳播的一個(gè)重要因素[10]。本研究采用SCIT方法進(jìn)行場(chǎng)群流行率和個(gè)體流行率估算,由于試驗(yàn)敏感性較低,引進(jìn)感染bTB奶牛的可能性較大。
表1 通過(guò)奶牛引進(jìn)傳播bTB的定量風(fēng)險(xiǎn)模型輸入?yún)?shù)
圖2 引進(jìn)1頭奶牛感染的概率模擬結(jié)果
圖3 至少有1頭奶牛感染的概率模擬結(jié)果
表2 模型中重要參數(shù)模擬結(jié)果和不確定性分析
圖4 至少有1頭奶牛感染的概率敏感性分析結(jié)果
通過(guò)分析發(fā)現(xiàn),隔離措施是奶牛引進(jìn)時(shí)降低bTB傳播風(fēng)險(xiǎn)的主要保護(hù)因素。因bTB的潛伏期一般為16~45 d,因此在采取嚴(yán)格隔離措施的同時(shí),要延長(zhǎng)隔離時(shí)間,并再次進(jìn)行抽樣檢測(cè),這樣可以有效降低bTB傳播風(fēng)險(xiǎn)。研究發(fā)現(xiàn),是否暴露于感染動(dòng)物是引進(jìn)時(shí)感染bTB的主要風(fēng)險(xiǎn)因素,因此在混群后,要提高生物安全措施,進(jìn)行嚴(yán)格消毒,同時(shí)加強(qiáng)奶牛臨床癥狀和體溫監(jiān)測(cè)。
風(fēng)險(xiǎn)評(píng)估顯示,萊西市奶牛場(chǎng)戶從感染地區(qū)引進(jìn)1頭奶牛傳播bTB的概率為2.93×10-3(95%CI:1.14×10-3~8.25×10-3); 引 進(jìn) 50頭奶牛至少有1頭傳播bTB的概率為0.137(95%CI:0.056~0.339)。敏感性分析結(jié)果顯示,暴露于感染動(dòng)物、群流行率和隔離措施,對(duì)引進(jìn)奶牛時(shí)傳播bTB的風(fēng)險(xiǎn)影響較大。建議在引進(jìn)奶牛時(shí)進(jìn)行嚴(yán)格隔離,延長(zhǎng)隔離時(shí)間,并在混群前再進(jìn)行1次bTB檢測(cè),混群后進(jìn)行嚴(yán)格消毒并加強(qiáng)監(jiān)測(cè),這樣可以降低引進(jìn)奶牛時(shí)傳播bTB的風(fēng)險(xiǎn)。