新疆烏蘇市羊棘球蚴病個(gè)體流行率測量

，，，，

（1. 新疆維吾爾自治區(qū)動(dòng)物衛(wèi)生監(jiān)督所，新疆烏魯木齊 830011；2. 中國動(dòng)物衛(wèi)生與流行病學(xué)中心，山東青島 266032）

新疆烏蘇市羊棘球蚴病個(gè)體流行率測量

閆 昊1，張 睿1，沈朝建2，林漢亮1，張繼紅1

（1. 新疆維吾爾自治區(qū)動(dòng)物衛(wèi)生監(jiān)督所，新疆烏魯木齊 830011；2. 中國動(dòng)物衛(wèi)生與流行病學(xué)中心，山東青島 266032）

[目的]了解新疆烏蘇市羊棘球蚴病真實(shí)流行率。[方法]首先用兩種酶聯(lián)免疫吸附試驗(yàn)試劑盒，對(duì)采自烏蘇市屠宰場的172份羊血清同時(shí)進(jìn)行抗體檢測，匯總交叉結(jié)果，通過建立貝葉斯模型，確定先驗(yàn)分布，完成Gibbs抽樣及相關(guān)參數(shù)計(jì)算，評(píng)價(jià)兩種試劑盒（A和B）對(duì)羊棘球蚴病的診斷價(jià)值，選擇診斷價(jià)值較高的試劑盒，對(duì)烏蘇市羊棘球蚴病個(gè)體真實(shí)流行率進(jìn)行測量。其次，采用橫斷面研究，采取兩階段抽樣策略。第一階段為鄉(xiāng)鎮(zhèn)抽樣，采用隨機(jī)抽樣，預(yù)定流行率為85%，CI為90%，可接受誤差為15%；第二階段為個(gè)體抽樣，采用按養(yǎng)殖大小成比例抽樣（PPS），預(yù)期流行率為30%，CI為90%，可接受誤差為5%。采集的1 787份血清，利用診斷價(jià)值取高的試劑盒進(jìn)行檢測，根據(jù)檢測結(jié)果、敏感性和特異性，計(jì)算真實(shí)流行率。[結(jié)果]通過貝葉斯模型分析，試劑盒A和B的敏感性分別為91.8%、76.6%，特異性分別為95.4%、80.7%，試劑盒A的診斷價(jià)值較高。烏蘇市羊棘球蚴病個(gè)體表觀流行率（AP）為48.6%，統(tǒng)計(jì)分析得出個(gè)體真實(shí)流行率（TP）為50.5%，95%置信區(qū)間（CI）為48.2%～52.8%。[結(jié)論]貝葉斯估計(jì)是抽樣信息對(duì)先驗(yàn)分布做出調(diào)整的結(jié)果，其更接近參數(shù)的真實(shí)情況。本研究結(jié)果既可以作為先驗(yàn)信息，也可以作為其他血清學(xué)檢查方法的評(píng)價(jià)和現(xiàn)場應(yīng)用提供參考信息。本研究首次通過兩階段抽樣策略對(duì)羊棘球蚴病的個(gè)體流行率進(jìn)行了測量，相對(duì)以往的監(jiān)測方法更加科學(xué)，更接近研究時(shí)段內(nèi)該地區(qū)的真實(shí)情況。在今后的羊棘球蚴病流行病學(xué)調(diào)查、監(jiān)測和檢測工作中可逐步推廣ELISA診斷方法。

羊棘球蚴?。涣餍新?；酶聯(lián)免疫吸附試驗(yàn)；貝葉斯模型；敏感性；特異性

棘球蚴?。‥chinococcosis）又稱為包蟲?。℉ydatid disease，Hd），是由細(xì)粒棘球絳蟲（Echinococcus granulosus，Eg）的中絳期——棘球蚴（Echinococcus）寄生在哺乳動(dòng)物臟器內(nèi)形成機(jī)械性壓迫，造成不同程度的組織損傷和功能障礙，同時(shí)可引起過敏反應(yīng)，甚至突然死亡，是一種嚴(yán)重的人獸共患寄生蟲病。綿羊是其最主要的中間宿主，犬是終末宿主。

關(guān)于棘球絳蟲屬的種分類，因其在蟲種上變異較大，曾先后報(bào)道有16個(gè)種和13個(gè)亞種[1]。目前在分類學(xué)上公認(rèn)的有5個(gè)種[2]，即細(xì)粒棘球絳蟲（Echinococcus granulosus，Eg）、多房棘球絳蟲（E. multilocularis，Em）、少節(jié)棘球絳蟲（E. oligauthrus，Eo）、福氏棘球絳蟲（E. vogeli，Ev）、石渠棘球絳蟲[3]（Echincococus shiquicus，Es）。目前我國仍以細(xì)粒棘球絳蟲蚴蟲在人體和獸體內(nèi)引起的囊型包蟲?。–ystic Echinococcosis，CE）和多房棘球絳蟲引起的泡型包蟲?。ˋlveolar Echinococcosis，AE）為主。據(jù)獸醫(yī)部門流行病學(xué)調(diào)查數(shù)據(jù)推算，我國每年患病家畜在5 000萬頭以上，因家畜死亡和臟器廢棄造成的直接經(jīng)濟(jì)損失逾30億元[2]。

目前羊棘球蚴病的感染調(diào)查仍以世界動(dòng)物衛(wèi)生組織（OIE）推薦的屠宰場現(xiàn)場檢查為主。該方法存在問題有：內(nèi)臟棘球蚴包囊檢查只能在家畜屠宰時(shí)進(jìn)行，存在很大局限性，同時(shí)家畜大多年齡偏大，可能感染其它動(dòng)物疫病，增加了流行病學(xué)調(diào)查的偏差。棘球蚴病感染初期或家畜年齡較小時(shí)，內(nèi)臟表面未能形成包囊或包囊不明顯，導(dǎo)致檢查出現(xiàn)假陰性。同時(shí)，檢查人員易將棘球蚴包囊與其它寄生蟲病或傳染病形成的結(jié)節(jié)狀病變混淆，導(dǎo)致出現(xiàn)假陽性。

近年來，多種新的免疫學(xué)方法被應(yīng)用于羊棘球蚴病檢測，其是利用抗原、抗體特異性結(jié)合的原理，制備特異性抗原，檢測血清中的抗體。血清抗體檢測方法的敏感性、特異性、實(shí)用性優(yōu)于傳統(tǒng)的檢測方法。

新疆全區(qū)棘球蚴病呈現(xiàn)地方性流行。根據(jù)2014年新疆家畜棘球蚴病流行病學(xué)調(diào)查數(shù)據(jù)，烏蘇市羊棘球蚴病肝臟感染率為30.15%、肺臟感染率為24.72%。烏蘇市位于新疆中部，轄10鎮(zhèn)7鄉(xiāng)，是以農(nóng)業(yè)為主、農(nóng)牧結(jié)合的半農(nóng)半牧區(qū)，有5個(gè)以牧業(yè)為主的牧場，養(yǎng)殖模式以傳統(tǒng)游牧為主。本文通過對(duì)兩種羊棘球蚴病檢測試劑盒進(jìn)行比對(duì)，選取其中一種對(duì)新疆烏蘇市羊棘球蚴病流行率進(jìn)行測量，從而掌握當(dāng)?shù)卣鎸?shí)流行情況。

1 試劑盒比對(duì)的材料和方法

1.1樣品

采自烏蘇市屠宰場的羊血清172份。

1.2主要試劑盒

A：綿羊棘球蚴?。òx?。〦LISA檢測試劑盒，蘭州獸醫(yī)研究所產(chǎn)品；B：包蟲抗體檢測試劑盒（酶聯(lián)免疫抑制法），北京天恩泰生物科技有限公司產(chǎn)品。

1.3檢測方法

綿羊棘球蚴病（包蟲?。〦LISA檢測試劑盒參照《動(dòng)物棘球蚴病診斷技術(shù)》（NY/T 1466—2007）進(jìn)行操作；包蟲抗體檢測試劑盒（酶聯(lián)免疫抑制法）參照說明書操作。

1.4貝葉斯模型建立

采用二試驗(yàn)一種群條件相關(guān)模型[4-5]。將兩種檢測試驗(yàn)結(jié)果與設(shè)立的向量y對(duì)應(yīng)參數(shù)一一對(duì)應(yīng)。檢測結(jié)果陽性用“+”“1”表示，檢測結(jié)果陰性用“-”“0”表示。交叉參數(shù)情況見表1。

表1 兩種檢測結(jié)果的交叉參數(shù)

用 向 量 p=（p11，p10，p01，p00） 表 示向量y對(duì)應(yīng)的發(fā)生概率，n為樣品總數(shù)，則y～Multinomial[n，（p11，p10，p01，p00）]。

假設(shè)感染率為π，兩種檢測方法的敏感性分別為Se1、Se2，特異性分別為Sp1、Sp2，CovD+為真值在陽性條件下兩種檢測方法結(jié)果之間的協(xié)方差，CovD-為真值在陰性條件下兩種檢測方法結(jié)果之間的協(xié)方差，則：

P10=π[Se1(1-Se2)-CovD+]+(1-π)[(1-Sp1) Sp2-CovD-]；

P01=π[(1-Se1)Se2-CovD+]+(1-π)[Sp1(1-Sp2)-CovD-]；

P00=π[(1-Se1)(1-Se2)+CovD+]+(1-π) (Sp1Sp2+CovD-)。

1.5確定先驗(yàn)分布

根據(jù)共軛先驗(yàn)分布原理，選擇β分布Beta（a，b）作為感染率、敏感性和特異性的先驗(yàn)分布[6]。由于未查閱到國內(nèi)對(duì)應(yīng)用這兩種檢測試劑盒的相關(guān)文獻(xiàn)報(bào)道，本研究結(jié)合前期預(yù)試驗(yàn)，以及專家經(jīng)驗(yàn)、意見，給出敏感性和特異性的取值范圍。然后，以區(qū)間中點(diǎn)為先驗(yàn)分布的期望值，以區(qū)間的1/4距離的平方為先驗(yàn)分布的方差[6]，聯(lián)立兩方程式，解之即可得到敏感性和特異性的先驗(yàn)分布參數(shù)（表2）。

表2 兩種檢測方法的先驗(yàn)分布參數(shù)

π由于缺乏先驗(yàn)信息，所以采用無先驗(yàn)信息分布Beta（1，1）。由于CovD+和CovD-缺乏先驗(yàn)信息，且(Se1-1)(1-Se2)≤CovD+≤min(Se1，Se2)-Se1Se2，(Sp1-1)(1-Sp2)≤CovD-≤min(Sp1，Sp2)-Sp1Sp2[7]，故取其先驗(yàn)分布為：

CovD+～Uniform[(Se1-1)(1-Se2)，min(Se1，Se2)-Se1Se2]；

CovD-～Uniform[(Sp1-1)(1-Sp2)，min(Sp1，Sp2)-Sp1Sp2]

(1)情報(bào)信息機(jī)構(gòu)科技在服務(wù)于創(chuàng)新領(lǐng)域的工作過程中，其知識(shí)服務(wù)的理論和實(shí)踐研究在2013年前總體呈上升態(tài)勢(shì)，并在2009年和2013年達(dá)到高位，之后其文獻(xiàn)發(fā)表數(shù)量開始呈下降趨勢(shì)，表明近幾年情報(bào)信息機(jī)構(gòu)在知識(shí)服務(wù)領(lǐng)域的研究活躍性降低。

1.6計(jì)算方法

根據(jù)上述分析，在WinBUGS1.4.2中編寫程序，完成Gibbs抽樣及相關(guān)參數(shù)的計(jì)算[8]。

2 流行率測量的材料和方法

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

為了測量羊棘球蚴病的真實(shí)流行率，本研究選擇采用橫斷面研究。烏蘇市轄區(qū)內(nèi)所有羊只為此次調(diào)查的總體樣本，采樣分為兩階段抽樣：第一階段為鄉(xiāng)鎮(zhèn)抽樣，第二階段為個(gè)體抽樣[9]。

2.2抽樣策略

第一階段，鄉(xiāng)鎮(zhèn)選取采用隨機(jī)抽樣策略。對(duì)轄區(qū)內(nèi)9個(gè)鄉(xiāng)鎮(zhèn)編號(hào)后作為抽樣框，根據(jù)專家意見，預(yù)定鄉(xiāng)鎮(zhèn)流行率為85%（鄉(xiāng)鎮(zhèn)估計(jì)流行率基于個(gè)體流行率＞10%，即確定為流行鄉(xiāng)鎮(zhèn)），置信水平（CL）為90%，可接受絕對(duì)誤差為15%，采用WinEpi軟件，按照估計(jì)流行率抽樣方法，計(jì)算抽取鄉(xiāng)鎮(zhèn)數(shù)，用EXCEL軟件的隨機(jī)數(shù)發(fā)生器抽取鄉(xiāng)鎮(zhèn)。

第二階段，個(gè)體抽樣采用按養(yǎng)殖規(guī)模大小，成比例抽樣（PPS）策略。根據(jù)專家意見以及當(dāng)?shù)貧v年監(jiān)測數(shù)據(jù)，個(gè)體預(yù)期流行率設(shè)定為30%，置信水平（CL）為90%，可接受誤差為5%，采用WinEpi軟件，按照估計(jì)流行率抽樣方法，計(jì)算每個(gè)縣個(gè)體樣本數(shù)，將每個(gè)鄉(xiāng)鎮(zhèn)2015年秋季重大動(dòng)物疫病免疫檔案信息整理編號(hào)后作為抽樣框，用EXCEL軟件的隨機(jī)數(shù)發(fā)生器抽取個(gè)體。

2.3實(shí)驗(yàn)室檢測及實(shí)驗(yàn)方法特性

將所有采集的血清樣品采用敏感性和特異性較高試劑盒進(jìn)行檢測，試驗(yàn)敏感性Se、特異性Sp通過試劑盒比對(duì)獲得。

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

羊棘球蚴病個(gè)體表觀流行率（AP）為檢測的血清學(xué)陽性數(shù)占所有檢測數(shù)的比值；個(gè)體真實(shí)流行率（TP）根據(jù)公式TP=（AP+Sp-1）/（Se+Sp-1）計(jì)算得到；真實(shí)流行病的95%置信區(qū)間（CI）根據(jù)公式CI=P±Z×[p×（1-p）/n]-1計(jì)算得出，其中，P是真實(shí)流行率，Z為標(biāo)準(zhǔn)正態(tài)分布概率p=0.05所對(duì)應(yīng)的值，n為抽樣數(shù)量。

3 結(jié)果

3.1試劑盒比對(duì)結(jié)果

用兩種試劑盒對(duì)172份血清的抗體檢測結(jié)果見表3，匯總交叉結(jié)果見表4。將相關(guān)參數(shù)值代入程序，進(jìn)行100 000次迭代。

表3 兩種試劑盒檢測結(jié)果 （單位：份）

表4 兩種檢測結(jié)果的交叉結(jié)果

通過監(jiān)視歷史路徑圖和自相關(guān)性，用隨后的99 500次迭代進(jìn)行參數(shù)估計(jì)（圖1），貝葉斯模型模擬的具體結(jié)果見表5。試劑盒A的敏感性和特異性估計(jì)值分別為91.8%、95.4%，試劑盒B的敏感性和特異性估計(jì)值分別為76.6%、80.7%，因此選用試劑盒A的診斷價(jià)值較高。

圖1 參數(shù)核密度情況

表5 兩種試劑盒敏感性和特異性的貝葉斯估計(jì)

3.2流行率測量結(jié)果

3.2.1抽樣結(jié)果。第一階段，隨機(jī)選取出6個(gè)鄉(xiāng)鎮(zhèn)（表6）。但由于經(jīng)費(fèi)、人員與實(shí)際情況等因素限制，只選取了4個(gè)鄉(xiāng)鎮(zhèn)（表7）。其中1個(gè)為半農(nóng)半牧區(qū)，3個(gè)為牧區(qū)，牲畜存欄量占總體的53.22%。地理分布情況見圖2。第二階段，計(jì)算每個(gè)鄉(xiāng)鎮(zhèn)需要抽取226份樣品，具體樣品來源見表6。

3.2.2檢測結(jié)果。根據(jù)個(gè)體流行率計(jì)算公式可得到各鄉(xiāng)鎮(zhèn)的個(gè)體表觀流行率、真實(shí)流行率及其置信區(qū)間，結(jié)果見表7。

表6 樣品來源 （單位：個(gè)/只/份）

表7 各鄉(xiāng)鎮(zhèn)檢測結(jié)果

圖2 鄉(xiāng)鎮(zhèn)抽樣地理分布

3.2.3各鄉(xiāng)鎮(zhèn)真實(shí)流行率的描述性統(tǒng)計(jì)。根據(jù)表7中列出的結(jié)果，對(duì)各鄉(xiāng)鎮(zhèn)羊棘球蚴病的真實(shí)流行率進(jìn)行描述性統(tǒng)計(jì)分析，結(jié)果見表8。

表8 各鄉(xiāng)鎮(zhèn)真實(shí)流行率的描述性統(tǒng)計(jì)

3.2.4個(gè)體流行率?？偣矊?duì)1 787份血清進(jìn)行檢測，檢出抗體陽性869份，根據(jù)個(gè)體流行率計(jì)算公式可得到烏蘇市羊棘球蚴病個(gè)體表觀流行率（AP）為48.6%。根據(jù)3.1中結(jié)果，該檢測方法敏感性Se=91.8%、特異性Sp=95.4%，計(jì)算所得個(gè)體真實(shí)流行率（TP）為50.5%，95%置信區(qū)間為48.2%～52.8%，

4 討論

4.1試劑盒比對(duì)

4.1.1貝葉斯模型建立。目前，國內(nèi)外眾多學(xué)者利用貝葉斯統(tǒng)計(jì)評(píng)價(jià)診斷試驗(yàn)診斷價(jià)值研究日益增多[11-16]。因?yàn)樨惾~斯統(tǒng)計(jì)基本原理是根據(jù)抽樣信息對(duì)先驗(yàn)信息做出調(diào)整的結(jié)果，所以它更接近參數(shù)的真實(shí)情況[7]。由于國內(nèi)對(duì)羊棘球蚴病不同ELISA診斷方法間的比對(duì)沒有相關(guān)文獻(xiàn)報(bào)道，本研究中的先驗(yàn)信息來自前期試驗(yàn)數(shù)據(jù)與相關(guān)專家意見，通過加工、分析后，引入到相關(guān)分析中。為了減少先驗(yàn)信息對(duì)模擬結(jié)果的影響，本次研究給出的兩種試劑盒的部分先驗(yàn)分布區(qū)間較寬（0.5～1.0），但本研究的結(jié)果可以作為今后相關(guān)棘球蚴病診斷試驗(yàn)診斷價(jià)值評(píng)價(jià)的先驗(yàn)信息，以供參考。

4.1.2ELISA方法檢測。兩種檢測試劑盒所用包被抗原均用棘球蚴病包囊囊液制備，不同感染時(shí)間下棘球蚴發(fā)育程度有所差異，導(dǎo)致試劑盒的敏感性、特異性產(chǎn)生偏移，所以在監(jiān)測過程中，對(duì)樣品背景信息和試劑盒信息的掌握顯得尤為重要。

4.1.3局限性。本研究選取的采樣地區(qū)屬于羊棘球蚴病高發(fā)區(qū)，同時(shí)樣品來源信息較為復(fù)雜，研究設(shè)計(jì)過程中未能考慮到不同流行地區(qū)，以及與其他疾病的交叉情況，可能對(duì)試驗(yàn)敏感性、特異性造成影響，故在其他地區(qū)采用ELISA方法檢測時(shí)，可能會(huì)出現(xiàn)不確定性。

4.2流行率測量

4.2.1抽樣策略。在第二階段抽樣過程中，將每個(gè)鄉(xiāng)鎮(zhèn)2015年秋季重大動(dòng)物疫病免疫檔案進(jìn)行整理，將養(yǎng)殖戶、飼養(yǎng)量等信息錄入采樣框，按養(yǎng)殖大小成比例抽樣，從而使轄區(qū)內(nèi)的羊被抽中的概率一致。由于當(dāng)?shù)仫曫B(yǎng)模式和飼養(yǎng)情況特殊（每戶養(yǎng)殖量在300～1 200只之間），當(dāng)養(yǎng)殖戶與采樣量確定后，進(jìn)入養(yǎng)殖戶采樣時(shí)，隨機(jī)抽取相應(yīng)數(shù)量，相對(duì)于將該鄉(xiāng)鎮(zhèn)全部羊編號(hào)后簡單隨機(jī)抽樣而言，該方法具有較好的可操作性和實(shí)用性，可在今后大范圍監(jiān)測抽樣時(shí)逐步推廣，有利于減少因便利采樣帶來的誤差。

4.2.2流行率測量。該地區(qū)歷年監(jiān)測數(shù)據(jù)來源于屠宰場現(xiàn)場檢查，數(shù)量為每年200只羊。監(jiān)測結(jié)果顯示該地區(qū)羊棘球蚴病流行率在12%～35%之間。本次研究首次通過兩階段抽樣，對(duì)羊棘球蚴病的個(gè)體流行率進(jìn)行了測量，相對(duì)以往的監(jiān)測方法更加科學(xué)，試驗(yàn)結(jié)果更接近研究時(shí)段內(nèi)該地區(qū)的真實(shí)情況。在牧區(qū)，養(yǎng)殖戶公用同一片草場，所有羊只在同一片草場活動(dòng)、采食，轉(zhuǎn)場時(shí)，經(jīng)過同一條牧道，所以本研究認(rèn)為，在牧區(qū)羊棘球蚴病感染情況同質(zhì)化水平很高，具體情況有待進(jìn)一步研究。

4.2.3防控策略。本次調(diào)查發(fā)現(xiàn)，該地區(qū)羊只調(diào)出量遠(yuǎn)遠(yuǎn)大于調(diào)入量，羊只飼養(yǎng)模式屬于牧區(qū)自繁自養(yǎng)、傳統(tǒng)游牧，而牧區(qū)飼養(yǎng)量占全地區(qū)的50%以上。牧區(qū)淘汰羊只收購至半農(nóng)半牧區(qū)進(jìn)行育肥飼養(yǎng)；在牧區(qū)，牧羊犬隨羊只活動(dòng)，而半農(nóng)半牧區(qū)的犬只以拴養(yǎng)為主，從犬與羊只接觸頻次來說，牧區(qū)要遠(yuǎn)遠(yuǎn)大于半農(nóng)半牧區(qū)。綜上所述，本研究認(rèn)為該地區(qū)羊棘球蚴病防控重點(diǎn)應(yīng)以牧區(qū)為主，半農(nóng)半牧區(qū)為輔，加大牧區(qū)犬只驅(qū)蟲與新生羔羊疫苗免疫力度，從而達(dá)到降低該地區(qū)羊棘球蚴病流行率的目的。

4.2.4局限性。由于受經(jīng)費(fèi)、人員等實(shí)際情況限制，只選取了4個(gè)鄉(xiāng)鎮(zhèn)進(jìn)行了第二階段抽樣，這可能造成選擇誤差。本研究認(rèn)為該地區(qū)同質(zhì)化水平較高，所以在第二階段抽樣過程中，加大了樣本采樣量，從設(shè)計(jì)的226份增加到450份，從而降低隨機(jī)誤差。

4.2.5公共衛(wèi)生評(píng)價(jià)?？紤]到棘球蚴病在我國的流行現(xiàn)狀及其成因[17]，根據(jù)此次研究數(shù)據(jù)和棘球蚴病特性，建議如下：（1）加快建立家/牧犬的登記管理制度，并集中收容、處理流浪犬只。（2）加大以“犬犬投藥、月月驅(qū)蟲”的實(shí)施力度，尤其是牧區(qū)。（3）對(duì)轄區(qū)內(nèi)新生羔羊?qū)嵤┭蚣蝌什∫呙缑庖?，尤其是牧區(qū)，提高其對(duì)棘球蚴病的抵抗力。（4）加強(qiáng)屠宰檢疫管理，杜絕“私屠濫宰”行為，嚴(yán)禁將有包囊的家畜臟器亂棄或投喂犬只。（5）加大棘球蚴病科普知識(shí)的宣傳力度，翻譯為多種民族語言，使農(nóng)牧民了解棘球蚴病的危害及防控知識(shí)。

5 結(jié)論

本研究通過對(duì)兩種羊棘球蚴病檢測試劑盒進(jìn)行比對(duì)，掌握了其各自診斷價(jià)值；選取一種檢測試劑盒，對(duì)新疆烏蘇市羊棘球蚴病流行率進(jìn)行了測量，從而掌握了當(dāng)?shù)卣鎸?shí)流行情況。研究得到以下結(jié)論：

屠宰場現(xiàn)場對(duì)家畜內(nèi)臟進(jìn)行棘球蚴病包囊檢查簡單易行、成本低廉，一直被用于羊棘球蚴病檢疫、監(jiān)測及流行病學(xué)調(diào)查。但是由于該方法的敏感性、特異性依賴于檢查者的專業(yè)技能、羊感染時(shí)間的長短、棘球蚴發(fā)育情況等因素的限制，易出現(xiàn)假陽性和假陰性，增加了流行病學(xué)調(diào)查的偏差。相比之下，血清抗體檢測方法的檢測敏感性、特異性、實(shí)用性要優(yōu)于傳統(tǒng)的檢測方法，為流行病學(xué)調(diào)查與疫情監(jiān)測、疾病診斷提供了新的技術(shù)手段和途徑。

通過建立貝葉斯模型，評(píng)價(jià)兩種羊棘球蚴病酶聯(lián)免疫吸附試驗(yàn)試劑盒對(duì)羊棘球蚴病的診斷價(jià)值，結(jié)果顯示試劑盒A的診斷價(jià)值較高，可在今后的羊棘球蚴病流行病學(xué)調(diào)查、監(jiān)測和檢測工作中逐步推廣該診斷方法。

該地區(qū)歷年監(jiān)測數(shù)據(jù)來源屠宰場現(xiàn)場檢查，監(jiān)測結(jié)果顯示該地區(qū)羊棘球蚴病流行率在12%～35%之間。本研究首次通過兩階段抽樣策略，對(duì)羊棘球蚴病的個(gè)體流行率進(jìn)行了測量，從方法上相對(duì)以往的監(jiān)測方法更加科學(xué)，監(jiān)測結(jié)果更接近研究時(shí)段內(nèi)該地區(qū)的真實(shí)情況。

致謝：

本研究是在導(dǎo)師沈朝建博士的悉心指導(dǎo)下完成的，在此謹(jǐn)向沈朝建導(dǎo)師表示最衷心的感謝和敬意！同時(shí)，也衷心感謝新疆維吾爾自治區(qū)動(dòng)物衛(wèi)生監(jiān)督所（新疆動(dòng)物疾病預(yù)防控制中心）盛卓君所長、林漢亮主任對(duì)本研究的大力支持。最后，特別感謝項(xiàng)目主辦方（聯(lián)合國糧農(nóng)組織中國辦公室）與協(xié)辦方（中國動(dòng)物衛(wèi)生與流行病學(xué)中心）各位導(dǎo)師、工作人員，F(xiàn)ETPV的特聘老師，以及CFETPVⅢ期培訓(xùn)班的學(xué)員對(duì)本研究提供的幫助。

[1]李炳軍. 包蟲病流行病學(xué)調(diào)查現(xiàn)狀[J]. 醫(yī)學(xué)綜述，2009，15（8）：1195-1198.

[2]張壯志，張文寶. 我國包蟲病防控及其面臨的困難[J].獸醫(yī)導(dǎo)刊，2011（6）：27-29.

[3]XIAO N. Biological Features of a New Echinococcus Species（Echinococcus shiqucus）in the East of Qinghai-Tibet Plateau[J]. Chin J Parasitol Parasit Dis，2008，26（4）：307-311.

[4]BRANSCUM A J，GARDNER I A，JOHNSON W O. Estimation of diagnostic-test sensitivity and speci fi city through Bayesian modeling[J]. Prevent Vet Med，2005，68（2/3/4）：145-163.

[5]茆詩松，呂曉玲. 數(shù)理統(tǒng)計(jì)學(xué)[M]. 北京：中國人民大學(xué)出版社，2011：143-159.

[6]顧海燕，陳啟光. 無金標(biāo)準(zhǔn)情況下診斷試驗(yàn)的評(píng)價(jià)方法[J]. 中國衛(wèi)生統(tǒng)計(jì)，1999，16（4）：203-205.

[7]孫世偉，范偉興. 綿羊布魯菌病兩種血清學(xué)診斷方法的貝葉斯評(píng)價(jià)[J]. 動(dòng)物醫(yī)學(xué)進(jìn)展，2014，35（9）：52-54.

[8]韋來生，張偉平. 貝葉斯分析[M]. 合肥：中國科學(xué)技術(shù)大學(xué)出版社，2013：235.

[9]黃保續(xù). 獸醫(yī)流行病學(xué)[M]. 北京：中國農(nóng)業(yè)出版社，2011：189.

[10]茆詩松，程依明，濮曉龍. 概率論與數(shù)理統(tǒng)計(jì)教程[M].北京：高等教育出版社，2011：335.

[11]劉沛，舒華嵩，張言超，等. 混合血清法篩檢丙型肝炎抗體的血清流行病學(xué)研究[J]. 中華流行病學(xué)雜志，1994，15（3）：178-181.

[12]顧海燕，陳啟光. 兩分類獨(dú)立診斷和篩檢試驗(yàn)的評(píng)價(jià)方法[J]. 中國衛(wèi)生統(tǒng)計(jì)，1998，15（2）：57-62.

[13]劉沛，陳啟光. 貝葉斯統(tǒng)計(jì)及其在診斷和篩選試驗(yàn)評(píng)價(jià)中的應(yīng)用[J]. 中國衛(wèi)生統(tǒng)計(jì)，2006，23（4）：361-363.

[14]尉潔，宋覺覺，等. 基于貝葉斯估計(jì)的診斷試驗(yàn)ROC曲線回歸模型[J]. 中國衛(wèi)生統(tǒng)計(jì)，2010，27（2）：152-154.

[15]劉沛，史志旭，等. 無金標(biāo)準(zhǔn)條件下診斷評(píng)價(jià)貝葉斯相關(guān)模型構(gòu)建及應(yīng)用[J]. 中國衛(wèi)生統(tǒng)計(jì)，2012，29（2）：187-190.

[16]楊麗，武海濱，等. 無金標(biāo)準(zhǔn)診斷試驗(yàn)靈敏性和特異度的貝葉斯估計(jì)方法[J]. 中國醫(yī)院統(tǒng)計(jì)，2011，18（2）：111-115.

[17]劉平，李金花，李印，等. 包蟲病在我國的流行現(xiàn)狀及成因分析[J]. 中國動(dòng)物檢疫，2016，33（1）：48-51.

（責(zé)任編輯：朱迪國）

Measurement on Individual Prevalence Rate of Ovine Echinococcosis in Wusu City of Xinjiang

Yan Hao1，Zhang Rui1，Shen Chaojian2，Lin Hanliang1，Zhang Jihong1
（1. Xinjiang Animal Health Supervision Institution，Urumqi，Xinjiang 830011；2. China Animal Health and Epidemiology Center，Qindao，Shandong 266032）

[Objective]In order to recognize the true prevalence of ovine Echinococcosis in Wusu city of Xinjiang. [Methods]Two kinds of ELISA kits（kit A and kit B）were used to detect antibodies existing in 172 collected sheep serum samples from slaughter houses. After a series of steps including summarizing of the cross-reaction results，confirming prior distribution by establishing Bayesian model and completing Gibbs sampling and related parameter calculations，the value of the two kinds of kits in diagnosis was evaluated to choose the one with higher value，so as to measure the true prevalence of individuals. Next，adopting ways of cross-sectional study and two-stage sampling strategy，1 787 serum samples were tested using kit A and the true prevalence rates were calculated on basis of the test results，sensitivity and speci fi city. As to the two-stage sampling， fi rst was random sampling in the counties，with an estimated prevalence of 85%，90% CI，and the acceptable error of 15%. The second stage was individual sampling by PPS strategy（sampling in proportion to the production scale），with an estimated prevalence of 30%，90% CI，and the acceptable error of 5%. [Results]Based on the results of Bayesian model analysis，the sensitivity of the kit A andB were 91.8% and 76.6%，respectively，their speci fi city were 95.4% and 80.7%，respectively. Hence kit A showed higher diagnostic value. The apparent prevalence（AP）of individuals was 48.6%. Through statistical analysis，the true prevalence（TP）of 50.5% was obtained（95% CI：48.2%～52.8%）. [Conclusion]Because Bayesian estimation is the result of adjusting sampling information in view of prior distribution，it is more close to the truth of parameters. The analysis results could be considered as both prior information and an assessment for other serological detection methods，as well as providing information for fi eld application. In this research，two-stage sampling strategy was used for the fi rst time to measure the individual prevalence rate of ovine echinococcosis. Compared with the past methods，it was more scienti fi c and much closer to the real situation of the region during the research phase. As a conclusion，the ELISA method could be gradually applied in fi elds of epidemiological survey，surveillance and detection of ovine echinococcosis.

ovine echincoccusis；prevalence rate；ELISA；eayesian model；sensitivity；speci fi city

S851.3

：A

：1005-944X（2017）02-0021-06

10.3969/j.issn.1005-944X.2017.02.006

聯(lián)合國糧農(nóng)組織-中國獸醫(yī)現(xiàn)場流行病學(xué)培訓(xùn)項(xiàng)目；科技部科技基礎(chǔ)性工作專項(xiàng)（2012FY111000）

張繼紅