摘 要:利用計算機(jī)模擬模型——北美動物疾病傳播模型(the North American Animal Disease Spread Model,NAADSM),對多個泰國清邁府和南奔府的豬偽狂犬?。ˋujeszky’s Disease)病毒根除策略進(jìn)行了比較。在模擬暴發(fā)期間,暴發(fā)持續(xù)時間、受感染豬群數(shù)和撲殺豬群數(shù)進(jìn)行了比較。同時對清群、為限制豬的活動范圍而進(jìn)行的劃分區(qū)域、改進(jìn)的檢測和免疫接種策略也進(jìn)行了評估。根據(jù)本研究的結(jié)果,根除豬偽狂犬病最有效的策略是采用清群策略和以受感染豬群為中心半徑為3 km、8 km和16 km劃出一個限制豬活動范圍的劃分區(qū)域策略,同時加強(qiáng)在感染豬群周圍半徑16 km范圍內(nèi)的根除和免疫接種措施。
關(guān)鍵詞:豬;偽狂犬?。桓?;撲殺;泰國
中圖分類號:S858 文獻(xiàn)標(biāo)識碼:C 文章編號:1001-0769(2015)08-0004-08
奧耶斯基?。ˋujeszky’s Disease,AD),又稱偽狂犬病,對全球養(yǎng)豬生產(chǎn)者構(gòu)成嚴(yán)重的經(jīng)濟(jì)威脅。這是一種病毒性疾病,臨床癥狀因感染時的年齡不同而各異。偽狂犬病毒(Pseudorabies virus,PRv)屬于皰疹病毒科下的皰疹病毒亞科。仔豬感染偽狂犬病后會導(dǎo)致中樞神經(jīng)系統(tǒng)紊亂。斷奶豬和育肥豬感染時首先侵染呼吸系統(tǒng),但神經(jīng)系統(tǒng)可能也會受波及。如果公豬和母豬受感染,AD可能會導(dǎo)致呼吸系統(tǒng)紊亂。據(jù)悉豬是偽狂犬病毒的儲庫,是除了靈長類動物和人類以外的大多數(shù)哺乳動物感染此疾病的源頭。
美國對偽狂犬病的首次介紹早至1813年,據(jù)報道牛感染后會發(fā)生嚴(yán)重的騷癢;這使這種疾病被冠以狂犬病的名稱。1902年,Aladar Aujeszky從狗、牛、貓身上分離出這種病毒,并證實它能使豬患上同樣的疾病。在亞洲,首次報道PRv的暴發(fā)在20世紀(jì)50年代的中國。后來,該病傳播到亞洲其他國家,包括1977年傳播到泰國。
PRv遍布全球各地,特別是養(yǎng)豬密集地區(qū),如南美、亞洲和歐洲。挪威、芬蘭或馬耳他未見有PRv感染的報道。德國、奧地利、瑞士、丹麥、英國、加拿大和新西蘭已將該病從國內(nèi)的豬群中根除。2004年,美國國內(nèi)的豬群成功根除PRv。在那些被認(rèn)為沒有PRv的國家,接種疫苗不被允許。然而,PRv仍在美國、德國、波蘭、法國、意大利和其他地方的野公豬或野豬群中流行。
正如前文提到的,PRv于1977年首次在泰國得到診斷,并仍在泰國的豬群中繼續(xù)流行。泰國有意愿消滅PRv。泰國政府通過畜牧發(fā)展部正在考慮執(zhí)行一項計劃,來重審有效的控制措施,并比較和評估哪些措施適用于泰國;這項計劃將使用計算機(jī)仿真模型——北美動物疫病傳播模型(North American Animal Disease Spread Model,NAADSM)來完成。
創(chuàng)建NAADSM的目的是建立用于模仿動物疫病在易感畜群中的傳播和控制,其特點是在畜群水平上得以體現(xiàn),而不是在動物個體上。這種疫病傳播模型是一種從易感到受感染和免疫的狀態(tài)傳播模型(圖1)。隨機(jī)過程被應(yīng)用于每一次模擬的疫病暴發(fā),結(jié)果包括由每一次模擬產(chǎn)生的一系列可能的后果。在某一次疫病發(fā)展中,每一個畜群被分配一個特定的緯度和經(jīng)度,同時疫病的發(fā)展以每天的步長表示。另外,NAADSM還包括一個成本會計部分,這對估算與疫病模擬暴發(fā)相關(guān)的費(fèi)用非常有幫助。因此,NAADSM是一種可評估商業(yè)性生產(chǎn)豬場中PRv根除計劃合適性的模型;可以評估泰國北部地區(qū)清邁府和南奔府PRv根除計劃的有效性;估計清邁府和南奔府在暴發(fā)PRv時所需疫苗的數(shù)量;估算政府執(zhí)行此類根除措施所用的成本(后續(xù)研究)。
1 材料和方法
1.1 地區(qū)和研究豬群
本研究在泰國北部地區(qū)的清邁府和南奔府(圖2)進(jìn)行。選中這個地區(qū)是因為可獲得詳細(xì)的數(shù)據(jù)。本研究包括育肥豬場、產(chǎn)仔-育肥豬場和父母代種豬場。
1.2 數(shù)據(jù)的來源
商業(yè)化豬場的地理坐標(biāo)由一個在清邁-南奔府地區(qū)、Nan府進(jìn)行無口蹄疫區(qū)域原型考察的研究團(tuán)隊收集。豬群用1984年全球大地坐標(biāo)系中的全球定位系統(tǒng)(GPS)進(jìn)行地理編碼。該研究團(tuán)隊還提供了有關(guān)豬群類型(育肥豬場、產(chǎn)仔-育肥豬場和父母代種豬場)和普查數(shù)據(jù)。任何重復(fù)數(shù)據(jù)均被刪除,如果數(shù)據(jù)有相同坐標(biāo)但屬于不同的所有者,這些數(shù)據(jù)會進(jìn)行整合。豬場規(guī)模根據(jù)豬的數(shù)量予以分為三類:小型(<500頭)、中型(500~5 000頭)和大型(≥5 000頭)。
1.3 分析中所包含的因素
本研究是在暗示性假設(shè)下進(jìn)行的:(1)所有的易感豬都是易感的;(2)所有受感染的豬同樣受感染,并且病毒在豬群中傳播;(3)所有感染PRv的豬最終都會表現(xiàn)出臨床癥狀。由于泰國以前沒有對PRv特點的研究,疾病參數(shù)可能與其他國家所見參數(shù)不同。
1.4 模型結(jié)構(gòu)
在本研究中,我們使用了NAADSM 3.1.23版。這個模型側(cè)重于豬群間的傳播,豬群被用作建模單位。為了控制豬群間的變量,豬群細(xì)分為3種不同類型:父母代種豬、產(chǎn)仔-育肥豬和育肥豬。
2 模型參數(shù)
2.1 疾病參數(shù)
NAADSM是一種狀態(tài)空間模型,也是一種群體型模型。我們給每一種模擬的豬群指定一種狀態(tài):易感型、潛伏型、亞臨床感染型、臨床感染型、自然免疫型、疫苗免疫型或死亡。在模擬這些模型的時候,豬群的感染狀態(tài)在這些狀態(tài)中轉(zhuǎn)化。表1列出了疾病過渡狀態(tài)的定義。由潛伏型豬群和其他類型豬群開始的每一種狀態(tài)的模型都是易感的。潛在感染的豬群從中部和高密度飼養(yǎng)地區(qū)選擇。這個中部地區(qū)被選中為以各個方向為二次傳播PRv提供易感豬群。
PRv過渡階段的疾病特點和時間周期根據(jù)暴發(fā)歷史、文獻(xiàn)描述和專家觀點進(jìn)行建模。概率密度函數(shù)用來描述每個豬群水平上的每種狀態(tài)的持續(xù)時間(表2)。我們使用Excel(@RISK Ver4.5)風(fēng)險分析軟件,根據(jù)專家意見使疾病特征符合研究豬群,從而選擇最佳概率密度函數(shù)。對每一種豬群類型,潛伏期被假定為遵循對數(shù)-邏輯分布。α、β、γ的值因豬群類型各異。產(chǎn)仔-育肥生產(chǎn)的值分別為-0.36、2.29和2.00。育肥豬生產(chǎn)的值為-0.42、2.32和2.15,父母代種群生產(chǎn)的值分別為-0.32、2.20和2.02。產(chǎn)仔-育肥和育肥生產(chǎn)類型的亞臨床感染期被認(rèn)為遵循高斯(正常)分布,平均值和標(biāo)準(zhǔn)差分別為5.41 d和0.88 d、14.43 d和1.19 d。父母代種豬被認(rèn)為遵循三角分布,最小、最可能和最大值為2.96 d、5.00 d和8.04 d。臨床感染期被認(rèn)為遵循對數(shù)-邏輯分布,所有的生產(chǎn)類型具有α、β、γ的各種值,產(chǎn)仔-育肥生產(chǎn)類型的值分別為24.85、22.49、2.25,而育肥生產(chǎn)類型的值分別為27.07、104.19和2.18,父母代種群類型的值分別為13.59、42.72和4.75。我們認(rèn)為每個接種疫苗的豬群全年能保持免疫力,覆蓋根據(jù)對數(shù)正態(tài)分布的時間(平均值和標(biāo)準(zhǔn)差分別為 300 d和60 d)。
2.2 傳播選項——接觸傳播和空氣傳播
2.2.1 接觸傳播
這種模型的傳播選項模擬3種同時傳播機(jī)制:直接接觸、間接接觸和空氣傳播。經(jīng)直接接觸的傳播基于模擬接觸或豬在受感染和易感豬群之間的流動。間接接觸基于模擬接觸或人、設(shè)備和車輛的流動,流動方向隨意。如果受感染豬群處于亞臨床感染或臨床感染,經(jīng)直接接觸、間接接觸和空氣傳播的傳染可能發(fā)生。疾病可在不同的生產(chǎn)類型間傳播。
接觸率是指每天來自豬舍的運(yùn)送或者外出接觸的平均數(shù)。正常情況下,對每個配對的生產(chǎn)類型而言,接觸率是單獨(dú)指定的;然而,在這個模型中,我們假設(shè)每個配對的生產(chǎn)類型的接觸率相同,并指定接觸率為0.03(直接)和2.1(間接)。因此,0.03是指豬群之間流動或接觸的平均數(shù)為100 d中3次,間接接觸率運(yùn)貨次數(shù)是每天2次(2.1)。我們使用感染傳播函數(shù)0.2(直接)和0.01(間接);如果接觸發(fā)生,這就是概率;它對轉(zhuǎn)移感染將足夠。
流動方向是隨意的,流動距離建立在兩個概率密度函數(shù)基礎(chǔ)上:一個是直接接觸,另一個是間接接觸。直接接觸和間接接觸的距離分布被認(rèn)為每一種接觸都是一個三角密度函數(shù),直接接觸和間接接觸的最小、最可能和最大值分別是0 d、40 d、190 d和0.5 d、25 d、60 d。
2.2.2 空氣傳播
空氣傳播或本地傳播按照與受感染豬場接近程度來模擬。空氣傳播選項包含的輸入?yún)?shù)有風(fēng)向、傳播速度、距源頭1 km的感染概率和最大傳播距離。在本研究中,風(fēng)傳播被認(rèn)為是隨意的(0°~360°)。來自源頭的疾病傳播速度(按照專家意見)以指數(shù)形式下降,距源頭1 km處每天的傳播概率設(shè)定為0.5。
2.3 控制策略——檢測、追蹤、劃分區(qū)域
2.3.1 檢測
這種模型只能檢測臨床感染豬群。感染豬群的檢測建立在作為輸入?yún)?shù)的兩個概率基礎(chǔ)上——觀察到癥狀的概率和報道癥狀的概率。檢測的總概率等于這兩個概率的乘積。描述檢測的概率是作為一個時間依賴性函數(shù)給出的。它被認(rèn)為具有100 %的特定性,無假陽性情況,同時檢測的這一天被檢測豬場被無意地檢疫。
觀察到臨床癥狀的概率表示假設(shè)感染已經(jīng)在豬群中出現(xiàn)了一些日子了,豬場主或獸醫(yī)向管理機(jī)構(gòu)報告PRv可疑癥狀的概率。觀察到癥狀的概率可根據(jù)生產(chǎn)類型單獨(dú)設(shè)定。對于產(chǎn)仔-育肥豬和父母代種豬這兩種生產(chǎn)類型,可以認(rèn)為超過60 d癥狀將惡化,概率呈線性關(guān)系,即在0 d、14 d和60 d時分別是1 %、50 %和90 %。對于育肥豬的生產(chǎn)類型來說,可以認(rèn)為癥狀在超過60 d時惡化,這種生產(chǎn)類型觀察到和報告癥狀的概率也呈線性關(guān)系,即在第0天、第7天和第60天時分別是0 %、15 %和85 %。
報告觀察到臨床癥狀的概率表示,基于豬場主和獸醫(yī)對最近一次暴發(fā)PRv的了解,豬群被報告給動物衛(wèi)生部門的概率。我們認(rèn)為所有生產(chǎn)類型的這種概率是相同的。時間依賴性函數(shù)最初在第0天為88 %,第8~14天則為899 %。
2.3.2 追溯
在這個模型中,我們對所有生產(chǎn)類型進(jìn)行了追蹤研究(trace-out)。如果接觸的接受者被成功追溯到,它將自動受到檢疫,還有可能提前被撲殺。假設(shè)所有生產(chǎn)類型的追溯都相同,追溯直接接觸會在檢測前60 d進(jìn)行模擬,追溯成功的概率達(dá)98 %。追溯間接接觸會在檢測前14 d進(jìn)行模擬,追溯成功的概率達(dá)80 %。
2.3.3 區(qū)域劃分(限制豬的流動)
區(qū)域劃分(限制豬的流動)是指創(chuàng)建一個圓形區(qū)域,用來將豬的流動限制在監(jiān)控區(qū)內(nèi)并加強(qiáng)檢測。我們對各種分區(qū)策略進(jìn)行了研究,并在撲殺策略后進(jìn)行補(bǔ)充。第一個策略是不劃分區(qū)域;在這種情況下既沒有限制豬的流動,也沒有加強(qiáng)檢測。其他策略包括劃分一個半徑3 km、 8 km和16 km的范圍。
半徑3 km區(qū)域?qū)λ猩a(chǎn)類型的作用是改變直接流動速度(即第0天為100 %,第2天急降為0 %,直到第14天仍保持在0 %)、間接流動速度(即第0天為100 %,第3~4天降為20 %,第14天直接上升到25 %)和檢測概率(即觀察到臨床癥狀的概率乘以2)。
半徑8 km區(qū)域?qū)λ猩a(chǎn)類型的作用是改變直接流動速度(即第0天為100 %,第3~4天降為25 %,第7天降為1 %)、間接流動速度(即第0天為100 %,第3~4天降為50 %,第7天降為25 %)和檢測概率(即觀察到臨床癥狀的概率乘以1)。
半徑16 km區(qū)域?qū)λ猩a(chǎn)類型的影響是改變檢測概率(即觀察到臨床癥狀的概率乘以1)。
3 根除策略——撲殺和接種疫苗
3.1 撲殺
本模型模擬了所有生產(chǎn)類型檢測到PRv后豬群的撲殺。假設(shè)延遲兩天執(zhí)行撲殺方案。據(jù)認(rèn)為撲殺能力在第14天時每天可多達(dá)5個豬場。撲殺的優(yōu)先權(quán)根據(jù)檢測結(jié)果、感染天數(shù)(感染天數(shù)越大,優(yōu)先權(quán)越高)和生產(chǎn)類型(父母代種豬、產(chǎn)仔-育肥豬,然后是育肥豬)來決定。所有檢測的豬群被認(rèn)定為將被撲殺。
3.2 接種疫苗
此模型中的所有生產(chǎn)類型都被模擬進(jìn)行接種疫苗。假設(shè)任何生產(chǎn)類型的兩個發(fā)病豬場在接種疫苗前必須檢測,在接種疫苗之后豬舍的免疫力會有一個14 d的延后期。給豬群接種疫苗的能力被認(rèn)為在第14天時可高達(dá)每天100個豬群,接種疫苗范圍的半徑長達(dá)16 km。接種疫苗優(yōu)先權(quán)基于接種理由(接種范圍大?。⑸a(chǎn)類型(父母代種豬、產(chǎn)仔-育肥和僅有育肥)以及感染天數(shù)。出于NAADSM 3.1.23版的局限性,疫苗在給整個免疫豬群帶來完全的免疫力上須100 %有效。
4 方案和結(jié)果
不同的方案在相同的豬群統(tǒng)計數(shù)、疾病傳播速度、疾病檢測下進(jìn)行研究,但除接種疫苗策略和分區(qū)域策略外。每個方案經(jīng)1 000次重復(fù)后得以建模。
本研究比較了撲殺和不撲殺兩種措施,然后增加了劃分區(qū)域(通過比較無區(qū)域、半徑 3 km和8 km,半徑3 km、8 km和16 km),最后再增加接種疫苗(通過比較不接種疫苗、半徑8 km范圍內(nèi)接種和半徑16 km范圍內(nèi)接種疫苗)。
結(jié)果因關(guān)注的主要項目不同而各異,包括豬群數(shù)量和受感染的豬,被檢測的豬群數(shù)量、接種疫苗的豬群數(shù)和豬數(shù)量、疫情暴發(fā)的持續(xù)天數(shù)。當(dāng)撲殺和接種疫苗完成時,疫情暴發(fā)被認(rèn)為結(jié)束,不再有潛伏或受感染的豬群。
5 結(jié)果
下面的結(jié)果是基于每一種方案1 000次重復(fù)模擬(或暴發(fā))后得到的結(jié)果。由模擬模型產(chǎn)生的所有可能性,有95 %用于概括輸出參數(shù),除非特別說明。最初比較的策略是豬群撲殺策略(表3)。當(dāng)這個模型中包含撲殺時,暴發(fā)的持續(xù)時間是165 d。相比較,如果這個模型不包括撲殺策略時,疫病就變成流行?。ū┌l(fā)超過30 000 d)。如果模型包含豬群撲殺策略,受感染的豬群和動物總數(shù)下降99 %。
表4總結(jié)了劃區(qū)策略的比較結(jié)果。如果考慮到材料和方法這部分介紹的劃分區(qū)域,疾病暴發(fā)持續(xù)時間從176 d~181 d不等(約6個月),具體取決于方案。這些方案中的每一種,受感染豬群檢測的首日是第20天。疾病檢測在所有方案中都是百分之百進(jìn)行的,即使這個模型認(rèn)定檢測具有百分之百的明確性。
對包含半徑3 km和半徑8 km區(qū)域的方案和無劃區(qū)(無流動限制)方案進(jìn)行比較,疾病暴發(fā)的持續(xù)時間多出1 d。大約減少了1 300頭豬(約占受感染豬總數(shù)的1 %)和7個豬群(約占受感染豬群總數(shù)的3 %)的感染。如果加以劃分區(qū)域,受到檢測并撲殺的豬和豬群總數(shù)分別下降1 %和3 %。
另一方面,如果增加半徑3 km、8 km和16 km的分區(qū)域方案,基于 95 %的無分區(qū)域方案暴發(fā)持續(xù)時間將多出5 d。如果增加劃分區(qū)域策略,受感染的豬減少了1 500頭(約占受感染豬總數(shù)的1 %)、豬群減少8個(約占受感染豬群總數(shù)的4 %)。有了劃分區(qū)域后,進(jìn)行檢測并撲殺的豬總數(shù)以及豬群的總數(shù)分別比不劃分區(qū)域時的撲殺豬數(shù)少1 %和4 %。如果接種疫苗,將半徑8 km和16 km的接種疫苗圈與無疫苗接種(表5)比較,半徑8 km和16 km接種疫苗圈內(nèi)的暴發(fā)持續(xù)時間均減少16 d(6 %)。如果接種圈是半徑8 km和16 km,受感染的豬群約減少2 %。半徑8 km接種圈的受撲殺豬群總數(shù)減少3 %,半徑16 km接種圈的受撲殺豬群總數(shù)減少4 %。半徑8 km接種圈,接種疫苗的豬群總數(shù)是420個。當(dāng)接種圈大小增加到半徑為16 km時,接種疫苗豬群的總數(shù)增加6 %。
6 討論
模擬模型的局限性在于它們無法預(yù)測未來,也不代表實時暴發(fā)情況;相反,它們用于幫助決策、規(guī)劃、明確潛在結(jié)果,評估建立在現(xiàn)有數(shù)據(jù)基礎(chǔ)之上的策略。本研究的結(jié)果只是清邁和南奔府商業(yè)化豬場的一次疫病暴發(fā),并不是全國或地區(qū)范圍的疫病。數(shù)據(jù)的準(zhǔn)確性會影響疾病傳播模型的有效性。我們獲得了自2004年以來的數(shù)據(jù),是本研究可用的最新和最準(zhǔn)確的數(shù)據(jù)。如果可以獲得時間更近的信息,豬群數(shù)據(jù)肯定需要修改。
本研究使用的模型——NAADSM 3.1.23版,有其局限性。該模型是為口蹄疫研究開發(fā)的,可能無法準(zhǔn)確地反映PRv的特征。潛伏期是皰疹病毒的一個主要特征,在這個模型中難以限制。此外,在本研究中,此模型假設(shè)豬的流動是隨意的,這一點在養(yǎng)豬行業(yè)中可能并不適用,特別是當(dāng)豬在一種生產(chǎn)系統(tǒng)中流動時。
在泰國,此疾病參數(shù)對PRv毒株并不是特有的,但相反是基于科學(xué)文獻(xiàn)報道的病毒特性。合理的是有不同特性的新病毒株可能出現(xiàn),結(jié)果可能偏離本研究。通常,疫苗并非100 %有效;然而,此模型假設(shè)免疫接種100 %有效。接種疫苗只能防止臨床癥狀出現(xiàn),但無法阻止病毒的感染。接種疫苗的有效性是指它能避免豬出現(xiàn)臨床癥狀的能力。由于NAADSM 3.1.23版不允許調(diào)整疫苗的有效性,暴發(fā)的實際持續(xù)時間(天數(shù))以及受感染和撲殺的豬或豬群的總數(shù)可能不同于本研究。準(zhǔn)備可替代的控制措施可能是有必要的。此模型中商業(yè)化生產(chǎn)豬被認(rèn)為無免疫力,這可能也會影響結(jié)果。
撲殺感染豬群(撲殺策略)是很多國家在試圖消滅PRv時采取的一種常用措施。我們將撲殺和不撲殺兩種策略進(jìn)行了比較。撲殺似乎是消滅PRv感染方法中最有益的一種;當(dāng)施以撲殺策略時,暴發(fā)持續(xù)時間(天數(shù))從超過30 000 d(7年)急劇減少到165 d(不滿1年),受感染的豬群和豬總數(shù)下降了99 %。
劃分區(qū)域(限制豬的流動范圍)策略繼撲殺策略之后實施,用以限制豬在監(jiān)控區(qū)域內(nèi)的流動,并加強(qiáng)檢測。本研究所檢查的3種劃分區(qū)域方案之間似乎并沒有有意義的差異,因為在采用劃分區(qū)域策略時,暴發(fā)持續(xù)時間增加了1 d~5 d。然而,如果我們應(yīng)用半徑3 km、8 km、16 km的劃分區(qū)域,感染似乎會減輕,并且必須被撲殺的豬群數(shù)量減少。
接種圈也作為一種控制和消滅PRv的手段得到研究。本模型在執(zhí)行撲殺和劃分區(qū)域(限制豬的流動范圍)策略之后應(yīng)用了接種疫苗策略。半徑8 km和16 km的接種圈使暴發(fā)持續(xù)時間減少了95 %,從181 d縮短到165 d。半徑16 km的接種圈在減輕暴發(fā)上似乎更有效,并可比半徑8 km的接種圈更多地減少感染率和撲殺率。
7 小結(jié)
顯然,根據(jù)本研究結(jié)果,撲殺策略在消滅PRv方面最有效。從暴發(fā)持續(xù)時間(天數(shù))、受感染、檢測和撲殺的豬和豬群總數(shù)來看,繼撲殺策略之后的劃分區(qū)域和圈內(nèi)接種疫苗在本模型中并沒有表現(xiàn)出重大的影響。然而,當(dāng)劃分區(qū)域策略和接種疫苗策略得到執(zhí)行后,暴發(fā)持續(xù)時間、感染數(shù)量和撲殺豬群均下降。因此,在后續(xù)研究完成之前,本項研究的建議是采用撲殺策略、3種分區(qū)(半徑3 km、8 km和 16 km)限制流動以及在半徑16 km接種圈加強(qiáng)檢測,以在清邁和南奔府消滅PRv。
8 后續(xù)研究
對未來研究的建議是在劃分區(qū)域和接種圈中做出最具成本效益的選擇,還建議評估消滅PRv的其他策略。
原題名:Evaluation of strategies for eradication of Aujeszky’s disease (Pseudorabies) in commercial swine farms in Chiang-Mai and Lampoon provinces, Thailand(英文)
原作者:N. Ketusing等