棘球蚴病的流行病学分布特征

曾祥嫚，官亚宜，伍卫平

棘球蚴病，俗称包虫病，是一种由棘球绦虫的幼虫寄生于动物或人体内而引起的人兽共患寄生虫病。该病有四种形式，即细粒棘球蚴引起的囊型棘球蚴病，多房棘球蚴引起的泡型棘球蚴病、伏氏棘球蚴引起的多囊型棘球蚴病及少节棘球蚴引起的单囊型棘球蚴病。囊型棘球蚴病和泡型棘球蚴病是人类感染此病的两种最主要形式。

棘球蚴病危害十分严重。同一时段，全球可能有超过100万人患有此病，其中很多人会患有严重的致命性临床综合征，即使经过治疗，也常常面临生活质量的下降。据估计，囊型棘球蚴病每年至少损失100万伤残调整生命年(DALYs)；泡型棘球蚴病每年损失约65万伤残调整生命年(DALYs)。棘球蚴病造成的疾病负担超过盘尾丝虫病，并与非洲锥虫病接近，每年与该病相关的治疗费用及牲畜业损失预计为30亿美元[1]，俨然成为全球重大公共卫生问题。

流行病学是研究人群中疾病与健康状况的分布及其影响因素，并研究防治疾病及促进健康的策略及措施的科学。现已广泛应用于传染性疾病、慢性非传染性疾病、突发公共卫生事件、伤害等各种与健康有关的状态及卫生事件的研究，并使之取得很大进展。本文对棘球蚴病的流行病学分布特征进行综述，旨在为我国对该病的预防控制提供帮助。

1 空间分布特征

1.1囊型棘球蚴病(CE)的空间分布 细粒棘球绦虫能适应多种中间宿主，遍及世界各大陆，涉及极地、温带、亚热带和热带地区。因此，全球除南极洲外，均有CE流行。

亚洲：西亚的塞浦路斯、土耳其、黎巴嫩、约旦、巴勒斯坦、以色列、叙利亚、科威特、沙特阿拉伯、伊朗，中亚的哈萨克斯坦、乌兹别克斯坦、吉尔吉斯斯坦、塔吉克斯坦和土库曼斯坦以及东亚的中国和蒙古有较严重的CE流行[2-5]。

欧洲：西班牙、意大利、法国、塞尔维亚、阿尔巴尼亚、希腊、波斯尼亚和黑塞哥维那等地中海国家正受威胁[2]。

南美洲：主要分布在秘鲁、智利、阿根廷、乌拉圭和巴西南部地区[6]。

北美洲：加拿大、阿拉加斯州等地细粒棘球绦虫循环于狼、犬与驼鹿、驯鹿和其它鹿科动物之间。另外，美国(除阿拉加斯州)、墨西哥也有病例报道(美国以输入性病例为主，墨西哥为本地感染病例)[6-7]。

非洲：尼日利亚、肯尼亚、埃塞俄比亚、苏丹、坦桑尼亚、摩洛哥、突尼斯、利比亚、阿尔及利亚、埃及是CE的重要流行区[8-9]。

1.2CE的空间分布特征 CE的流行是个复杂的过程，受到多种因素的影响。如生物因素——细粒棘球绦虫有10个株G1～G10，其中G1，G2主要在犬-羊中循环，G3主要在犬-水牛中循环，G4主要在犬-马中循环，G5主要在犬-牛中循环，G6主要在犬-骆驼中循环，G7和G9主要在犬-猪中循环，G8和G10主要在狼-鹿中循环，不同株细粒棘球绦虫致病力不同，导致CE流行亦有差异[9]；自然因素——温度、湿度、降雨量通过影响棘球绦虫虫卵在环境中的寿命而影响棘球蚴病的流行[10]；社会因素——生产方式、文化习俗、社会变革也能影响CE流行[10]。因此，在不同国家、同一国家不同地区之间CE流行不同。

1.2.1不同国家之间流行不同 总的来说，亚非等第三世界国家流行严重，部分国家或地区牛、绵羊、山羊、骆驼、猪的感染率分别达40%、60%、60%、60%和20%以上；人群患病率达5%、发病率达30/10万以上；犬细粒棘球绦虫感染率达30%以上。如王虎等[5]于1995-2005年在中国青海柴达木盆地、青南高原、祁连山地和河湟谷地等地区调查显示绵羊感染率分别为50.70%、57.31%和52.99%，牦牛感染率分别为52.90%、57.21%和53.03%，山羊感染率分别为47.91%、3.5%和17.51%，猪感染率分别为21.91%、4.89%和3.67%；Koskei等[11]在埃塞俄比亚阿萨拉、默克莱和德布雷地区分别检查14、15和15只犬，细粒棘球绦虫感染率为50%、20%和33.3%；Ernest等[12]在坦桑尼亚恩戈罗恩戈罗地区5个屠宰场调查357头牛、105只绵羊和619只山羊，感染率分别为48.7%、63.8%和34.7%；Nelson等[13]研究中显示，在图尔卡纳与苏丹、埃塞俄比亚、乌干达交接的东北部地区，人群患病率5%～10%，外科CE发病率达198/10万。

1.2.2同一国家不同地区流行亦有差异 Omer等[14]于2001-2003年在苏丹中部、西部、南部地区调查的4 893头牛中感染率为6%，中部地区最高为20%(50/250)，其次为南部7.1%(23/325)、西部地区5.2%(226/4 318);调查的10 422只绵羊中感染率为11%，西部地区最高为11.9%，其次为南部2.7%(23/325)、中部地区2.5%(10/400)；调查的5565只山羊中感染率为2%，南部地区7.1%(3/42)，较西部地区1.9%(103/5 523)高；调查的779头骆驼中，西部地区61.4%(347/565)，较中部地区55.6%(119/214)高。Azlaf等[15]在摩洛哥里夫、卢科斯、拉巴特、卡萨布兰卡、阿特拉斯山脉、南部地区调查618头牛、2 948只绵羊、2 337只山羊、482头骆驼(均来自南部地区)和455匹马(均来自中部地区)，感染率分别为22.98%、10.58%、1.88%、12.03%和17.80%。其中阿特拉斯山脉、里夫和南部地区牛较其它家畜感染率高；卢科斯牛、羊均有较高的感染率，分别为37.61%和31.65%；中部地区马感染率最高17.80%，其次为牛10.29%、绵羊2.66%，未发现山羊感染。在中国，CE在四川、青海、新疆、宁夏、甘肃、内蒙古等地流行。感染率：牦牛10.47%～57.21%，绵羊9.03%～57.31%，山羊3.5%～60%，猪3.67%～21.91%，人群患病率0.22%～8.7%，犬细粒棘球绦虫感染率13.2%～33%[5,16-19]。另外，西藏也有较为严重的CE流行。龚学红等[20]收集1996-2000年那曲地区人民医院、第一人民医院、第二人民医院CE病例709例，其中第一人民医院、第二人民医院收集的162例，那曲占46.9%，拉萨占21.6%，山南占13.6%，昌都占6.2%，林芝占4.4%，日喀则占2.5%，阿里占1.2%。

1.3泡型棘球蚴病(AE)的空间分布特征 AE的流行同样受生物因素，如宿主的种类及密度；自然因素，如海拔、气温、地理景观(地形地貌、地理景观构成、植被覆盖类型、草场覆盖比例)；社会因素，如经济条件、生产生活方式、文化习俗、生活习惯的影响[10,21]。因此，其发生具有一定的地域性，多发生在植被覆盖较高的山区，这些地区主要分布着针叶松森林、草原和高比例的自然植被，具有较低的年平均大气气温、较高的年平均降水量和较低的年平均土壤表面温度[21]。如中国青藏高原地区，欧洲，北美及中亚一些国家、地区。

中国：四川、青海高度流行。在四川，对石渠、甘孜两县233只黑唇鼠兔、14只灰尾兔、12只松田鼠和70只小家鼠调查，阳性率分别为5.6%、7.1%、25%和0%；对16只犬、171只藏狐检查，多房棘球绦虫感染率分别为16.98%和44.4%；对7 138位居民进行B超检查，AE患病率3.1%[22-23]。在青海，柴达木盆地灰仓鼠感染率1.80%；青南高原田鼠、灰尾兔、小家鼠感染率分别为0.45%、12.50%和0%；祁连山地和河湟谷地高原鼠兔感染率1.8%；称多、甘德、班玛、久治、泽库人群AE调查表明患病率分别为0.8%、0.71%、5.09%、2.58%和0.29%[5]。另外，宁夏、甘肃、新疆等地也有流行[17-18]。

欧洲：主要分布在奥地利、德国南部、瑞士和法国东部，原来认为的非流行区，如波兰、斯洛伐克、立陶宛、斯洛文尼亚、比利时、匈牙利等也有病例报告[24-26]。

北美：以前认为仅限于加拿大和阿拉加斯州北部苔原地带，后来发现其它地区也有流行。野生动物感染多见[27]。

中亚：主要分布在草原、山区，如吉尔吉斯坦、哈萨克斯坦部分地区[3,28]。

2 棘球蚴病的时间分布特征

2.1CE的时间分布特征 一些国家、地区经过防治，流行有所下降。如塞浦路斯，19世纪70年代以前CE严重流行，外科CE年发病率12.9/10万。1971年，“国家畜牧局”启动防治规划。之后，认为CE已在犬和家畜之间消灭，先后于1974年、1985年在北部、南部地区终止活动。然而，1989年后在家畜中陆续出现CE暴发流行，遂于1993年重新启动防治规划。2000-2003年，犬细粒棘球绦虫感染率为0.012%，2004年牛和羊感染率分别为6.61%和1.53%，山羊感染率低于0.5%，均较以前有所下降[2]；火地岛1976-1977年犬感染率80%，绵羊感染率55%，1979年后，在部分地区实行防治，1997-1999年对防治地区20个农场5 916只绵羊调查，感染率降为1.1%[29]。希腊1967-1971年人群发病率14.8/10万，之后，国家采取检验肉类及管理流浪狗等措施，2008年人群发病率降为0.3/10万[30]。

然而，也有一些国家、地区(特别是中亚部分国家)流行趋于严重。哈萨克斯坦1974-1994年外科CE发病率0.9-1.4/10万，1994-2003年上升至6.4/10万。患儿的比例也明显上升(1993-2000年患儿比例从19.9%上升至32.2%)[3]。乌兹别克斯坦1990年牛、绵羊、山羊、猪、骆驼、驴的感染率分别为24.5%、45.1%、8.0%、6.0%、25.4%和32.0%，2002年分别上升至45.5%、62.2%、11.1%、8.0%、35.0%和38.5%。2002-2010年全国新发病例8 014例，其中2002年2 421例，2010年上升至2 966例[3,31]。吉尔吉斯斯坦1991年人群年发病率5.4/10万，2000年上升至18/10万。比什凯克某医院1990年约有CE病人21例，1998年、1999年分别上升至127例和124例。儿童CE也有增加：1990年仅2例，2000年上升至82例。另外，在北部农村对1 486人进行超声检查，患病率约1.35%，根据医院报告的肝肺囊肿患者所占比例，推算农村地区有高达3.4%的人群处于CE亚临床期[32]。土库曼斯坦2000年人群平均发病率17/10万，比1988年增加3倍，其中发病率最高的阿哈尔达30/10万[3]。

2.2AE的时间分布特征 欧洲和中亚一些国家地区发病呈上升趋势，且在欧洲，原来认为的非流行区(波兰、斯洛伐克、立陶宛、斯洛文尼亚、比利时、匈牙利等)也有病例报告并呈上升趋势。奥地利1991-2000年，人群AE年发病率2.4/10万，2001年-2010年上升至2.8/10万[33]。瑞士1993-2000年人群年均发病率0.10/10万，2001-2005年上升至0.26/10万，其中2003年最高达0.38%，约28例新发病例。赤狐密度1980年约1万左右，1995-2004年达4万/年[24]。波兰1990-2004年AE新发病例不到50例，2005-2009年上升至55例，2010-2011年仅23个月时间23例。另外，赤狐密度由1995年67 000上升至2006年220 000，平均感染率2.6%，最高地区达62.9%[25]。立陶宛某医院1997-2001年新发病例0-4例/年，2002-2006年上升至10-16例/年[26]。吉尔吉斯坦1996年发现第1例AE病人，以后每年约有0-3例新病例。然而，自2004年开始，新发病例呈逐年上升趋势，至2011年约有新病例61例。另外，调查的151只赤狐，多房棘球绦虫的感染那率为64%，平均感染虫数8 669[28]。

3 棘球蚴病的人群分布特征

3.1性别特征 从世界各国的资料分析来看，棘球蚴病的患病率大多为女性高于男性[5,18,34]。男女性别比在0.40～0.86∶1之间。但是约旦的调查资料显示，仅15岁以上年龄组女性高于男性，男女性别比0.40～0.80∶1；15岁以下年龄组女性低于男性，男女性别比1.6∶1[35]。

3.2年龄特征 几乎所有的研究都发现，棘球蚴病的患病率有随着年龄增加逐渐增高的趋势。Shambesh等[34]研究显示0～4岁年龄组人群患病率最低为0.06%，15～30岁年龄组达2%，是前者的33倍，而50岁以上年龄组人群患病率更高，达4%。研究者最初认为出现这种情况的原因是因为包囊生长缓慢，早期感染者多数为无症状感染者。无法被发现。随着年龄的增长包囊逐渐增大，被发现的可能性增加，使得患病率增加。但实际上，各年龄组的小包囊是没有差异的。即各年龄组被感染的几率是一样的。王虎等[5]研究中1～10岁年龄组人群患病率0.66%，11～20岁、20～30岁年龄组明显上升，分别为2.12%、4.89%。61～70岁年龄组达到最高12.88%。且各年龄组的分布差异存在着统计学差异。这可能是在长期的生产过程中，年龄越大，参与生产、生活活动越多，接触虫卵机会越多，使感染具有累积性和重复性。Yang等[18]研究中CE在0～10岁年龄组无患者，21～30岁年龄组明显上升，60岁以上组达到最高为6.0%；而AE也是0～10岁组无患者，30～40岁组明显上升为2.2%，60岁以上组达到最高为6.8%。

3.3其它特征 不同职业人群棘球蚴病患病率不同，如在中国，大量研究表明牧民的患病率最高[5,16-18,22]。农村人口发病率高于城市人口。Acosta-Jamett等[36]在智利科金博地区收集1995-2006年外科CE病例，发现农村人口居多的Limarí和Choapa地区外科CE发病率分别为6.2～10.8/10万和1.0～4.8/10万，高于城镇人口居多的Elqui地区(1.5～3.1/10万)，差异具有统计学意义。另外，不同民族、宗教信仰，棘球蚴病的流行也会不同[5,22]。

4 结 语

鉴于以上方面，在今后研究中，应根据棘球蚴病流行病学分布特征，进行相关流行因素的研究，从而为制定科学合理的综合防控策略和措施提供依据。

参考文献：

[1]Budke CM, Deplazes P, Torgerson PR. Global socioeconomic impact of cystic echinococcosis[J]. Emerg Infect Dis, 2006, 2(2): 296-303.DOI: 10.3201/eid1202.050499

[2]Dakkak A. Echinococcosis/hydatidosis: a severe threat in Mediterranean countries[J]. Vet Parasitol, 2010, 174(1-2): 2-11. DOI: 10.1016/j.vetpar.2010.08.009

[3]Torgerson PR, Oguljahan B, Muminov AE, et al. Present situation of cystic echinococcosis in Central Asia[J]. Parasitol Int, 2006, 55(Suppl): S207-212. DOI: 10.1016/j.parint.2005.11.032

[4]Jabbar A, Narankhajid M, Nolan MJ, et al. A first insight into the genotypes of Echinococcus granulosus from humans in Mongolia[J]. Mol Cell Probes, 2011, 25(1): 49-54. DOI: 10.1016/j.mcp.2010.11.001

[5]Wang H, Zhang JX, Schantz PM, et al. Epidemiologic survey and analysis on echinococcosis in humans and animals from 1995 to 2005 in Qinghai province[J]. Chin J Zoonoses, 2006, 22(12): 1129-1133. DOI:10.3969/j.issn.1002-2694.2006.12.011 (in Chinese )

王虎, 张静宵, Schantz PM, 等. 1995-2005年青海省棘球蚴病流行病学调查分析[J]. 中国人兽共患病学报, 2006, 22(12): 1129-1133.

[6]Moro P, Schantz PM. Cystic echinococcosis in the Americas[J]. Parasitol Int, 2006, 55(Suppl): S181-S186. DOI: 10.1016/j.parint.2005.11.048

[7]Steta J, Torre A. Mexican-native human echinococcosis:case report of an underestimated disease[J]. Ann Hepatol, 2009, 8(3): 251-254.

[8]Wahlers K, Menezes CN, Wong ML, et al. Cystic echinococcosis in sub-Saharan Africa[J]. Lancet Infect Dis, 2012, 12(11): 871-880. DOI: 10.1016/S1473-3099(12)70155-X

[9]Sadjjadi SM. Present situation of echinococcosis in the Middle East and Arabic North Africa[J]. Parasitol Int, 2006, 55(Suppl): S197-202. DOI: 10.1016/j.parint.2005.11.030

[10]Otero-Abad B, Torgerson PR. A systematic review of the epidemiology of echinococcosis in domestic and wild animals[J]. PLoS Negl Trop Dis, 2013, 7(6): e2249. DOI: 10.1371/journal.pntd.0002249

[11]Koskei P, Janitschke K, Feseha G. Prevalence ofEchinococcusgranulosusin some selected sites of Ethiopia[J]. East Afr J Public Health, 2011, 8(3): 170-175.

[12]Ernest E, Nonga HE, Kassuku AA, et al. Hydatidosis of slaughtered animals in Ngorongoro district of Arusha region, Tanzania[J]. Trop Anim Health Prod, 2009, 41(7): 1179-1185. DOI: 10.1007/s11250-008-9298-z

[13]Nelson GS. Hydatid disease: research and control in Turkana, Kenya. 1. Epidemiological observations[J]. Trans R Soc Trop Med Hyg, 1986, 80(2): 177-182.

[14]Omer RA, Dinkel A, Romig T, et al. A molecular survey of cystic echinococcosis in Sudan[J]. Vet Parasitol, 2010, 169(3-4): 340-346. DOI: 10.1016/j.vetpar.2010.01.004

[15]Azlaf R, Dakkak A. Epidemiological study of the cystic echinococcosis in Morocco[J]. Vet Parasitol, 2006, 137(1-2): 83-93. DOI: 10.1016/j.vetpar.2006.01.003

[16]Wang GZ, He XH, Chu XD, et al. Epidemiological study on human echinococcosis in Hobukesar Mongolian Autonomous County of Xinjiang[J]. Chin J Endemiol, 2009, 28(2): 214-217. DOI: 10.3760/cma.j.issn.1000-4955.2009.02.027 (in Chinese)

王桂芝, 冯晓辉, 初向东, 等. 2007年新疆和布克赛尔蒙古自治县棘球蚴病现况调查[J]. 中国地方病学杂志, 2009, 28(2): 214-217.

[17]Zhao YM, Zhong SG, Cai XP, et al. Epidemiological survey on hydatidosis in Tibetan region of southern Gansu[J]. Chin J Epidemiol, 2008, 29(11): 1100. (in Chinese)

赵玉敏, 种世桂, 才学鹏, 等. 甘南藏区包虫病流行病学调查[J]. 中华流行病学杂志, 2008, 29(11): 1100.

[18]Yang YR, Sun T, Li Z, et al. Community surveys and risk factor analysis of human alveolar and cystic echinococcosis in Ningxia Hui Autonomous Region, China[J]. Bull World Health Organ, 2006, 84(9): 714-721.

[19]Tao B, Zhang B, Du BB. Epidemiological survey and effect of prevention of echinococcosis in Inner Mongolian Autonomous Region[J]. Chin J Zoonoses, 2011, 27(7): 677-678. DOI: 10.3969/j.issn.1002-2694.2011.07.025 (in Chinese)

涛波, 张斌, 杜宝彪. 内蒙古自治区包虫病流行病学调查与防治效果[J]. 中国人兽共患病学报, 2011, 27(7): 677-678.

[20]Gong XH, Ci WY, Ze YG, et al. A clinical analysis of 709 cases with cystic echinococcosis in Tibet Autonomous Region[J]. Chin J Parasital Parasit Dis, 2001, 19(2): 128. DOI: 10.3969/j.issn.1000-7423.2001.02.025 (in Chinese)

龚学红, 次旺仁增, 泽永革, 等. 西藏自治区709例囊型棘球蚴病的临床分析[J]. 中国寄生虫学与寄生虫病杂志, 2001, 19(2): 128.

[21]Miterpáková M, Dubinsky P, Reiterová K, et al. Climate and environmental factors influencingEchinococcusmultilocularisoccurrence in the Slovak Republic[J]. Ann Agric Environ Med, 2006, 13(2): 235-242.

[22]Wang Q, Qiu J, Yang W, et al. Socioeconomic and behavior risk factors of human alveolar echinococcosis in Tibetan communities in Sichuan, People's Republic of China[J]. Am J Trop Med Hyg, 2006, 74(5): 856-862.

[23]He JG, Qiu JM, Liu FJ, et al. Epidemiological survey on hydatidosis in Tibetan region of western Sichuan Ⅱ. Infection situation among domestic and wild animals[J]. Chin J Zoonoses, 2000, 16(5): 62-65. DOI: 10.3969/j.issn.1002-2694.2000.05.019 (in Chinese)

何金戈, 邱加闽, 刘凤洁, 等.四川西部藏区包虫病流行病学研究Ⅱ.牲畜及野生动物两型包虫病感染状况调查[J]. 中国人兽共患病杂志, 2000, 16(5): 62-65.

[24]Schweiger A, Ammann RW, Candinas D, et al. Human alveolar echinococcosis after fox population increase, Switzerland[J]. Emerg Infect Dis, 2007, 13(6): 878-882. DOI: 10.3201/eid1306.061074

[25]Nahorski WL, Knap JP, Pawlowski ZS, et al. Human Alveolar Echinococcosis in Poland: 1990-2011[J]. PLoS Negl Trop Dis, 2013, 7(1): e1986. DOI: 10.1371/journal.pntd.0001986

[26]Bruzinskaite R, Marcinkute A, Strupas K, et al. Alveolar echinococcosis, Lithuania[J]. Emerg Infect Dis, 2007, 13(10): 1618-1619. DOI: 10.3201/eid1310.061161

[27]Storandt ST, Kazacos KR.Echinococcusmultilocularisidentified in Michigan with additional records from Ohio[J]. J Parasitol, 2012, 98(4): 891-893. DOI: 10.1645/GE-3057.1

[28]Ziadinov I, Deplazes P, Mathis A, et al. Human alveolar echinococcosis in Kyrgyzstan. Frequency distribution of Echinococcus multilocularis and other helminths of foxes in Kyrgyzstan[J]. Vet Parasitol, 2010, 171(3-4): 286-292. DOI: 10.1016/j.vetpar.2010.04.006

[29]Zanini F, Gonzalo R, Pérez H, et al. Epidemiological surveillance of ovine hydatidosis in Tierra del Fuego, Patagonia Argentina, 1997-1999[J]. Vet Parasitol, 2006, 138(3/4): 377-381. DOI: 10.1016/j.vetpar.2006.01.051

[30]Sotiraki S, Chaligiannis I. Cystic echinococcosis in Greece. Past and present[J]. Parasite, 2010, 17(3): 205-210.

[31]Hong ST, Jin Y, Anvarov K, et al. Infection status of hydatid cysts in humans and sheep in Uzbekistan[J]. Korean J Parasitol, 2013, 51(3): 383-385. DOI: 10.3347/kjp.2013.51.3.383

[32]Torgerson PR, Karaeva RR, Corkeri N, et al. Human cystic echinococcosis in Kyrgystan: an epidemiological study[J]. Acta Trop, 2003, 85(1): 51-61.

[33]Schneider R, Aspock H, Auer H. Unexpected increase of alveolar echinococcosis, Austria, 2011[J]. Emerg Infect Dis, 2013, 19(3): 475-477. DOI: 10.3201/eid1903.120595

[34]Shambesh MA, Craig PS, Macpherson CN, et al. An extensive ultrasound and serologic study to investigate the prevalence of human cystic echinococcosis in northern Libya[J]. Am J Trop Med Hyg, 1999, 60(3): 462-468.

[35]Al-Qaoud KM, Craig PS, Abdel-Hafez SK. Retrospective surgical incidence and case distribution of cystic echinococcosis in Jordan between 1994 and 2000[J]. Acta Trop, 2003, 87(2): 207-214. DOI: 10.1016/S0001-706X(03)00022-6

[36]Acosta-Jamett G, Cleaveland S, Cunningham AA, et al.Echinococcusgranulosusinfection in humans and livestock in the Coquimbo region, north-central Chile[J]. Vet Parasitol, 2010, 169(1-2): 102-110. DOI: 10.1016/j.vetpar.2009.12.009