吳秀娟,李 印,楊顯超,李凱航,周錦萍
(1. 上海市動物疫病預(yù)防控制中心,上海 201103;2. 中國動物衛(wèi)生與流行病學中心,山東青島 266032)
上海市滬淮批發(fā)市場從某省引進感染H7N9流感病毒活雞的定量風險分析
吳秀娟1,李 印2,楊顯超1,李凱航1,周錦萍1
(1. 上海市動物疫病預(yù)防控制中心,上海 201103;2. 中國動物衛(wèi)生與流行病學中心,山東青島 266032)
2013年3月底,上海市報道了全國首例人H7N9流感確診病例。隨后全國其他地區(qū)又有多例病例被確診,而且多被證實有活禽接觸史。根據(jù)調(diào)查,上海市市場售賣的活禽大多來自外?。ㄊ校渲袦磁l(fā)市場的活禽全部來源于某省大型養(yǎng)殖公司,主要是黃羽肉雞。本研究利用定量風險分析的方法,評估滬淮批發(fā)市場在2015年12月從該省引進至少1羽感染H7N9流感病毒活雞的可能性。根據(jù)情景樹,建立了定量風險分析模型,并利用@RISK軟件,模擬計算出2015年12月至少有1羽活雞攜帶H7N9流感病毒進入批發(fā)市場混群的概率為11.9%(95%CI:6.22%~19.82%)。該結(jié)果的不確定性為中等。模型敏感性結(jié)果顯示,檢測技術(shù)敏感性、H7N9流感個體流行率和H7N9流感場群流行率是重要參數(shù)。評估結(jié)果提示,養(yǎng)殖場應(yīng)提高生物安全防護意識,有效降低養(yǎng)殖場中病原存在的概率。建議相關(guān)部門進一步優(yōu)化檢測技術(shù),提升產(chǎn)地檢疫效能,通過多種方式減少活禽交易量,推廣規(guī)范化交易模式,提高交易各環(huán)節(jié)的生物安全水平。
H7N9流感;批發(fā)市場;活禽調(diào)運;定量風險分析
高致病性禽流感作為重要的人獸共患病,不僅給家禽養(yǎng)殖業(yè)造成巨大經(jīng)濟損失,也給人類健康帶來重大威脅。過去的十多年,在亞洲、中東、東歐以及非洲等地,數(shù)以億萬的家禽因高致病性禽流感死亡或被撲殺[1-2]。
2013年3月31日以來,上海、安徽、江蘇、浙江、河南、北京、山東、福建、江西等省(市)陸續(xù)發(fā)生多起H7N9流感病毒感染人事件,引發(fā)民眾恐慌,引起全球關(guān)注。
自2013年4月在活禽市場樣品中首次檢測到H7N9流感病毒陽性后,上海市立即關(guān)閉了全市所有活禽交易市場,暫停一切活禽交易,撲殺活禽交易市場內(nèi)所有禽類,并對市場進行徹底清洗和消毒。這些措施有效切斷了H7N9流感病毒的傳播途徑,阻止了病毒向人間進一步擴散[3]。據(jù)流行病學調(diào)查發(fā)現(xiàn),上海市是典型的活禽輸入?yún)^(qū)域,市場上售賣的活禽大多來自外省市。因此,H7N9流感病毒由外部傳入上海市的可能性較大。另外,黃羽肉雞是上海市活禽交易市場的主要交易品種,交易量約占禽類總交易量的70%,因此病毒的傳播很可能與肉雞的跨省調(diào)運有關(guān)[4]。
上海市活禽市場交易網(wǎng)絡(luò)調(diào)查結(jié)果顯示,滬淮批發(fā)市場的活禽全部來源于某省大型養(yǎng)殖公司。本研究利用定量風險分析的方法,評估滬淮批發(fā)市場在2015年12月從該省引進至少1羽肉雞攜帶H7N9流感病毒的可能性,為今后H7N9流感的預(yù)防與控制奠定基礎(chǔ)。
1.1 數(shù)據(jù)來源和分析方法
通過前期問卷調(diào)查和現(xiàn)場調(diào)查獲得的活禽交易網(wǎng)絡(luò)數(shù)據(jù);通過文獻檢索獲得的H7N9流感流行率。對所有數(shù)據(jù)通過EXCEL進行統(tǒng)計分析,采用@Risk軟件,模擬測算至少有1羽感染H7N9流感病毒活雞進入滬淮批發(fā)市場的概率。
1.2 風險評估模型建立的假設(shè)
1.2.1 家禽來源的假設(shè)。滬淮批發(fā)市場的家禽均來自某省2個大型家禽公司的交易平臺。此交易平臺的家禽來源復雜,可能來自多個簽約農(nóng)戶,但本研究的風險評估模型假定每批雞均來自同一個養(yǎng)殖場戶。
1.2.2 H7N9流感場群流行率的假設(shè)。因數(shù)據(jù)缺乏,本研究中某省的H7N9流感場群流行率采用全國的場群流行率,并假設(shè)不同地區(qū)的H7N9流感場群流行率差異不顯著。
1.2.3 運輸過程中的假設(shè)。由于每批雞都有產(chǎn)地檢疫證,因此假設(shè)每羽活雞都經(jīng)過產(chǎn)地檢疫并實施嚴格的陽性淘汰,并且運輸過程中H7N9流感群內(nèi)個體流行率不變。
1.3 風險評估模型中的參數(shù)設(shè)定
1.3.1 H7N9流感場群流行率(HP)。依據(jù)2013年全國家禽H7N9流感監(jiān)測數(shù)據(jù)(共檢測12 792個場群,陽性場群35個),用Beta分布模擬場群流行率。
1.3.2 H7N9流感群內(nèi)個體流行率(P)。依據(jù)2015年下半年禽流感定點調(diào)查數(shù)據(jù),2014年10—12月,江蘇、上海、福建3個省(市)活禽市場中的活雞H7N9流感個體檢出率為1.01%,用Beta分布模擬場群內(nèi)個體流行率。
1.3.3 活禽批發(fā)市場調(diào)運數(shù)據(jù)。上海市活禽網(wǎng)絡(luò)調(diào)查數(shù)據(jù)顯示,滬淮批發(fā)市場2015年12月調(diào)入活雞28批次,總共38 385羽,平均每批次1 371羽。
1.3.4 產(chǎn)地檢疫檢測方法的敏感性和特異性。采用血凝抑制試驗,檢測活雞血清中的H7N9流感抗體。參考專家意見,設(shè)定該試驗的敏感性(Se)為85%~95%,特異性(Sp)為80%~90%。采用Pert分布,模擬該試驗的敏感性和特異性。
2.1 情景樹
繪制感染H7N9流感病毒活雞進入滬淮批發(fā)市場的情景樹(圖1)。
圖1 感染H7N9流感病毒活雞進入滬淮批發(fā)市場的情景樹
2.2 至少有1羽活雞攜帶H7N9流感病毒的概率
根據(jù)情景樹,按照下列公式計算至少有1羽活雞在混群時攜帶H7N9流感病毒的概率。
HP為場群流行率,P為個體流行率,Se和Sp分別為產(chǎn)地檢疫檢測方法的敏感性和特異性。
經(jīng)過計算并模型仿真10 000次,得到至少1羽活雞在混群時攜帶H7N9流感病毒的概率為3.33×10-6(95% CI:1.67×10-6~5.77×10-6)。 具體結(jié)果見圖2。
圖2 至少1羽活雞在混群時攜帶H7N9病毒的概率分布
2.3 每批次至少有1羽活雞攜帶H7N9病毒進入市場的概率
每批次活雞的平均數(shù)量為1 371羽。按照下列公式計算每批次至少有1羽活雞攜帶H7N9流感病毒進入市場混群的概率。
經(jīng)過模型仿真5 000次,得到每批次至少有1羽活雞攜帶H7N9流感病毒進入市場混群的概率為0.46%(95%CI:0.23%~0.79%)。具體結(jié)果見圖3。
圖3 每批次至少有1羽活雞攜帶H7N9流感病毒的概率分布
2.4 至少有1羽活雞攜帶H7N9流感病毒進入市場的概率
2015年12月,滬淮批發(fā)市場共采購28批次活雞。按照下列公式計算該月至少有1羽活雞攜帶H7N9流感病毒進入市場的概率。
經(jīng)過模型仿真10 000次,得到至少有1羽活雞攜帶H7N9流感病毒進入市場混群的概率為11.9%(95%CI:6.22%~19.82%)。具體結(jié)果見圖4。
圖4 至少有1羽活雞攜帶H7N9流感病毒進入市場的概率分布
2.5 模型的敏感性分析
分析結(jié)果顯示,在定量風險分析模型中,檢測技術(shù)敏感性、H7N9流感個體流行率和H7N9流感場群流行率是最重要的參數(shù)(圖5)。
圖5 模型中的參數(shù)敏感性分析結(jié)果
2.6 模型的不確定性分析
對定量風險分析模型中應(yīng)用到的參數(shù),進行定性的不確定性分析,發(fā)現(xiàn)定量分析模型的不確定性為中等(表1)。
表1 模型中參數(shù)不確定分析結(jié)果
通過定量分析模型分析,得到2015年12月至少有1羽活雞攜帶H7N9流感病毒進入滬淮批發(fā)市場混群的概率為11.9%(95%CI:6.22%~19.82%)。該結(jié)果的不確定性為中等。此結(jié)果因建立在諸多假設(shè)之上,具有一定的局限性。
首先,采用的場群流行率和個體流行率都是監(jiān)測結(jié)果中的病毒陽性檢出率??紤]到病原學試驗的敏感性和特異性以及是否處于排毒期等因素,H7N9流感的實際場群流行率、個體流行率和文中引用的數(shù)據(jù)可能存在偏倚。其次,2015年的H7N9流感流行情況較2013年會有所變化,其流行可能更為普遍。該家禽來源省作為養(yǎng)殖大省,其H7N9流感場群流行率可能高于全國流行率。再次,本文假設(shè)每羽活雞都經(jīng)過產(chǎn)地檢驗并實施陽性淘汰,但交易網(wǎng)絡(luò)調(diào)查結(jié)果顯示,并非每批進入滬淮批發(fā)市場的活雞都經(jīng)過嚴格地產(chǎn)地檢疫和陽性淘汰。最后,由于滬淮批發(fā)市場的電子數(shù)據(jù)經(jīng)人工輸入,無法避免出現(xiàn)漏輸、錯輸?shù)痊F(xiàn)象。這些因素可能導致本模型的最終計算結(jié)果存在偏倚,因此活雞攜帶H7N9流感病毒進入滬淮批發(fā)市場的實際概率可能比本文計算的結(jié)果更高。另外,批發(fā)市場內(nèi)存在不同批次、種類活禽混合飼養(yǎng)現(xiàn)象,也存在當天售賣不出去的活禽在市場內(nèi)滯留的現(xiàn)象(約2~5天)。這些都為H7N9流感病毒在批發(fā)市場內(nèi)增殖和傳播提供了有利條件。
本研究只考慮了活雞攜帶H7N9流感病毒進入批發(fā)市場的風險,并且假設(shè)運輸途中活雞不可能通過污染的籠具感染,群體內(nèi)個體流行率不變。但現(xiàn)實中,攜帶H7N9流感病毒進入批發(fā)市場的不僅僅是活雞,還有其他種類活禽,也存在有關(guān)介質(zhì)因素,如籠具、車輛、人員等。
依據(jù)模型敏感性結(jié)果分析發(fā)現(xiàn),如果滬淮批發(fā)市場引進活雞數(shù)量減少,感染H7N9流感病毒的活雞進入市場的風險就會降低。該分析結(jié)果為今后上海市制定H7N9流感防控策略奠定了基礎(chǔ)。
基于本定量風險分析結(jié)果,給予以下建議:(1)在養(yǎng)殖生產(chǎn)環(huán)節(jié),加強宣傳干預(yù),提高生物安全防護意識,有效降低養(yǎng)殖場中病原存在的概率。(2)在產(chǎn)地檢疫環(huán)節(jié),進一步優(yōu)化檢測技術(shù),以養(yǎng)殖場為單元,正確處理陽性活禽,提升產(chǎn)地檢疫效能。(3)在活禽調(diào)運環(huán)節(jié),一要通過多種方式減少活禽交易量,如引導冰鮮家禽消費,引入信息技術(shù)平臺,推行終端到終端的交易模式,減少流通環(huán)節(jié),進而降低活禽調(diào)運導致病毒傳播和擴散的風險;二要建立活禽交易模式或平臺示范,推廣規(guī)范化活禽交易模式,提升交易各環(huán)節(jié)的生物安全水平。
[1] MAGALH?ES R J S,ORTIZ-PELAEZ A,THI K L L,et al. Associations between attributes of live poultry trade and HPAI H5N1 outbreaks:a descriptive and network analysis study in northern Vietnam. [J]. BMC Vet Res,2010,6(1):10.
[2] KERKHOVE M D V,VONG S,GUITIAN J,et al. Poultry movement networks in Cambodia: implications for surveillance and control of highly pathogenic avian influenza(HPAI/H5N1)[J]. Vaccine,2009,27(45):6345-6352.
[3] HE Y,LIU P,TANG S,et al. Live poultry market closure and control of avian influenza A (H7N9),Shanghai,China[J]. Emerg Infect Dis,2014,20(9):1565-1566.
[4] ZHOU X,LI Y,WANG Y,et al. The role of live poultry movement and live bird market biosecurity in the epidemiology of in fl uenza A(H7N9):A cross-sectional observational study in four eastern China provinces[J]. J Infect,2015,71(4):470-479.
(責任編輯:朱迪國)
Quantitative Risk Analysis on the Entry of Chickens Carrying H7N9 Virus into Huhuai Wholesale Market of Shanghai City
Wu Xiujuan1,Li Yin2,Yang Xianchao1,Li Kaihang1,Zhou Jinping1
(1. Shanghai Animal Disease Prevention and Control Center,Shanghai 201103;2. China Animal Health and Epidemiology Center,Qingdao,Shandong 266032)
The first confirmed case of human infection with influenza A(H7N9)virus was reported in Shanghai city in late March 2013. Subsequently,many other infection cases were also con fi rmed in other areas of China,and contact history with live poultry in majority of the con fi rmed cases was found. According to investigation,the live poultry sold in markets mostly came from other provinces and cities outside Shanghai city. Among the live bird markets,poultry(mainly yellow broiler)in the Huhuai wholesale market were all derived from a large-scale poultry farm of one province. Using quantitative risk analysis,the possibility of introducing at least one yellow broiler carrying H7N9 in fl uenza virus into the market from that province in December 2015 was evaluated. According to the scenario tree,a quantitative risk analysis model was established. By @RISK analog calculation,the rate of the incident was 11.9%(95%CI :6.22%~19.82%),and the uncertainty level was medium. Results of sensitivity analysis of the model indicated three important parameters of the sensitivity of the detection technology,individual and herd prevalence of H7N9. The assessment results indicated poultry farmers should raise awareness of bio-safety containment and effectively reduce the existence of pathogens in poultry farms. In addition,some suggestions were also given:further optimizing detection techniques to improve ef fi ciency in inspection and quarantine at the place of origin,using multiple ways toreduce live poultry trading volume,popularizing the standardized trading model,and promoting the bio-safety level in each link of the live poultry trade networks.
H7N9 influenza;wholesale markets;allocation and transportation of live poultry;quantitative risk analysis
S851.3
B
1005-944X(2017)03-0005-04
10.3969/j.issn.1005-944X.2017.03.002
上?;钋菔袌鼋灰拙W(wǎng)絡(luò)與禽流感病毒陽性率的相關(guān)性研究(滬農(nóng)青字〔2016〕第2-4號)
周錦萍