• Nem Talált Eredményt

5. CHAPTERS

5.3. CHAPTER 3

Benzimidazole resistance within red deer, roe deer and sheep populations within a joint habitat in Hungary

Click to BUY NOW!

.tracker-software.c Click to BUY NOW!

.tracker-software.c

ContentslistsavailableatScienceDirect

SmallRuminantResearch

j o u r n al ho me p ag e :w w w . e l s e v i e r . c o m / l o c a t e / s m a l l r u m r e s

Benzimidazoleresistancewithinred deer,roedeerand sheep populationswithinajointhabitatinHungary

GáborNagya,∗,ÁgnesCsivincsikb,LászlóSugárc,AttilaZsolnaid

aKaposvárUniversity,DepartmentofAnimalNutrition,GubaS.u.40.,H-7400Kaposvár,P.O.Box16.,Hungary

bKaposvárUniversity,InstituteofDiagnosticImagingandRadiationOncology,H-7400Kaposvár,GubaS.u.40.,Hungary

cKaposvárUniversity,DepartmentofPhysiologyandAnimalHygiene,H-7400Kaposvár,GubaS.u.40.,Hungary

dKaposvárUniversity,DepartmentofGameBiologyandEthology,H-7400Kaposvár,GubaS.u.40.,Hungary

a r t i c l e i n f o

Articlehistory:

Received22December2016

Receivedinrevisedform8February2017 Accepted9February2017

Theanthelminticresistanceofgastrointestinalnematodesisoneofthemostimportant,economicrisk factorsingrazingruminantsystems,allovertheworld.Wehaveinfinitesimalinformationaboutthe resistancestatusofnematodesindeerspecies.OuraimwastodeterminethepresenceofBZresistance inthegeneralistworm,Haemonchuscontortusinpasturedsheepandfreerangingredandroedeerby RFLP-PCRmethodbasedonthedetectionPhe200Tyrsinglenucleotidepolymorphism.Byinvestigation of70wormsfromeachhostspecies,thehomozygoussusceptiblegenotypewasthemostrepresentative inthereddeer(100%),thehomozygousresistantgenotypewasmostprevalentinthesheep(68.6%)and moderateintheroedeer(17.1%),whiletheheterozygousgenotypewasobservedinequalproportionin thesheepandroedeer(28.6%).Ourresultssuggestthatoverlappinghabitatsofsheepflocksandroedeer couldcontributetotheoccurrenceandspreadofresistantallelewithinwildlife.

©2017ElsevierB.V.Allrightsreserved.

1. Introduction

Theanthelminticresistance(AR)ofgastrointestinalnematodes isoneofthemostimportanteconomicriskfactorsingrazing rumi-nantsystemsallovertheworld(SutherlandandLeathwick,2011;

Papadopoulosetal.,2012;Roseetal.,2015).Thefirstappearingof resistanceagainstbroadspectrumdrugsdatesbacktothe1960’s andits spreading couldbeobservedin everyclass (FAO, 2004;

Kaplan,2004).

Although,thephenomenonofARiswellknowninmajor domes-ticatedruminantspecies,asinthecaseofsheep(Barrereetal., 2013;Pe ˜na-Espinozaetal.,2014;Chaudhryetal.,2015);within theminorfood-producingspecies(eg.deerspecies),wehavemuch less information about the AR status of theirnematodes. In a study(Mackintoshetal.,2014),theresistancetomoxidectinand abamectinwasdemonstratedbyfaecaleggcountreduction test (FECRT)inreddeer.TheauthorssuggestedthatARof Ostertagia-type species was in the background of drug ineffectiveness. A similarresultwasshowedin thecase of albendazoleand iver-mectin (Mylrea et al., 1991) in fallow deer. Both studies were carriedoutunderfarm circumstances.It washypothesizedthat

Correspondingauthor.

E-mailaddress:nagy.gabor@ke.hu(G.Nagy).

failureoftreatmentswasaffectedbyseveralfactorsincludingthe ARofgastrointestinalnematodes.Chintoan-Utaetal.(2014)tested andconfirmedthecrosstransmissionofHaemonchuscontortusand its benzimidazole(BZ) resistancebetween thefallow deer, red deer,roedeer,cattle,and sheep.Bothoftheirresistance detec-tionmethods(egghatchtestandmoleculartest)demonstratedthe presenceofARinthewormsderivedfromroedeer.Threesingle nucleotidepolymorphisms(SNPs)inthe-tubulingeneisotype1 havebeenassociatedwithbenzimidazoleresistanceinH. contor-tus.ThePhe200Tyrmutationhasbeenfoundineverycountryata relativelyhighfrequency.TheSNPsatcodons167and198have alsobeenreportedin multiplecountriesbuthave a more vari-ableoccurrenceand aregenerallypresent ata lowerfrequency thanthePhe200Tyrmutation(Chaudhryetal.,2016).Sincethe mostcommonmolecularmechanismthat confersBZresistance intrichostrongyles insmallruminantsinvolvesa phenylalanine totyrosinemutationatcodon200of -tubulingene isotype1;

therefore,moleculardiagnosticmethodsbasedonthedetection ofthisSNPareusedmostfrequentlytoidentifyBZresistanceinH.

contortus(Colesetal.,2006).

A Hungarian study (Nagy et al., 2016)demonstrated a high homozygousresistantgenotype(70%)andresistanceallele(85%) proportioninfarmedreddeer.Theresultsofamolecular diagnos-ticmethodshowedthatroutine,long-termadministrationofthe

http://dx.doi.org/10.1016/j.smallrumres.2017.02.012 0921-4488/©2017ElsevierB.V.Allrightsreserved.

Click to BUY NOW!

w.tr

acker-software.c w Click to BUY NOW!

.tracker-software.c

samedrugcouldcontributetotheappearanceofARinH.contortus inthereddeer.

Recently,anewapproach,refugiatheory,aroseduringthe inves-tigationof ARin natural environment. It hypothesizes thatthe presenceofwildruminantsonpasturesoflivestockprovidesa sig-nificantpercentageofhomozygoussensitiveonesamongthetotal amountofinfectiouslarvae.Regularinfectionbythismixed pop-ulationimpedestheunlimitedspreadofresistancewithinaworm populationinfectingalivestockflockconcerned.(VanWyk,2001).

In this study, we aimed to determine the presence of BZ resistance in the generalist worm, Haemonchus contortus in a microregion,wheretheresourcespartlysharedbysheepandfree rangingredandroedeer.

2. Materialsandmethods 2.1. Studysite

Our investigationwascarried out in southwesternHungary, withinacontiguousareawithoutanynaturalorman-built isolat-ingelementsthatimpedemovementsofanimals(Fig.1).Thestudy site(approximately30km2)characterizedby145–276maltitude abovesealevel,sub-Mediterraneanclimate,withsome submon-taneoushabitatpatches,with10Cannualmeantemperatureand 630–800mmannualprecipitation.Theproportionofforestsand agriculturalareasonstudysiteare about55%and 45%, respec-tively.Inthecoreofthesite,a6000-hamonoblockforestexists withthedominationofoaks(Quercusspp.),limes(Tiliaspp.), horn-beam(Carpinusbetulus)andEuropeanbeech(Fagussylvatica).The agriculturalareas(viz.grasslands,pastures,arablelands,andold orchards)arethemosttypicalontheperipheryofthestudysiteand theyprovideaheavilyfragmentedlandscapestructure.The aver-agedensityofreddeer,roedeer,andsheepare1.71,0.84and0.17 animal/km2,respectively.Ourdatawerebasedonhunting statis-tics(hunteddeerspecimen/km2);whileinthecaseofsheep,itwas calculatedbytheofficialregistryoftheHungarianSheepandGoat BreedersAssociation.

2.2. Wormcollectionandmoleculardiagnosticprocedure

DuringtheintervalfromSeptemberof2013tillAugustof2015, wecollected38abomasifromreddeer(N=14),roedeer(N=14) andsheep(N=10from5differentflocks).Theorgansofdeerwere collectedfromhuntingbagsinregularhuntingseasons,whilesheep weresampledataregionalabattoir.Afterevisceration,each abo-masumwasplacedseparatelyintoaplasticbagimmediately,and eachwasstoredat 18Cuntilexamination.Aftermelting,the organswereopenedalongsidethebigcurvature,andthecontent wasplacedintoaplasticjar,whilethemucosawaswashed thor-oughly.Inthecaseofdeerspecies,5Haemonchuscontortusmales werepickedrandomlyfromeachorgan,while7specimenswere collectedfromeverysinglesheepabomasum.Formorphological identification,weusedLichtenfels’setal.(1994)work.Thatmethod isbasedonthelengthoftheleftandrightspiculebarbs.

Inordertogenotypethecollectedworms,weusedthe detec-tionofPhe200Tyrsingle-nucleotidepolymorphismsoncodon200 of-tubulingeneisotype1,whichisthemostcommonmolecular mechanismconferringBZresistanceintrichostrongylesofsmall ruminants(Colesetal.,2006).DNAlysatesweremadeseparately from210 adultmale worms. The appliedRestriction Fragment Length Polymorphism-PCR method and the primer sequences (AvikaF:5-CTACCCTTTCCGTCCATCAA-3 andAvikaR:5- TGAA-GACGAGGGAATGGAAC -3) were previously detailed by Tiwari etal.(2006).PCRreactionswereperformedinatotalvolumeof 10l,containing200MofeachdNTP,0.2Mprimers,10×PCR

buffer,0.5unitDynazymeDNApolymerase(FinnzymesOy,Espoo, Finland)and 100ng genomicDNA.ThePCR cyclingprofile con-sistedofdenaturationat94Cfor3min,45cyclesofdenaturationat 94C(for30s),annealingat56C(for30s),andextensionat72C (for30s),followed byafinalextensionat 72C for5min.After amplification,1lTaaIrestrictionendonuclease(5U/l;Thermo FisherScientificInc.,Waltham,Massachusetts,USA),and1.22l 10×digestionbufferwereaddedtothetotalPCRvolume. Diges-tionwascarriedoutat65C,overnight.Digestedfragmentswere resolvedon4%agarosegelstainedwithGRGreenNucleicAcidGel StainandvisualizedunderUVlight.Genotypedeterminationwas basedonthefragmentlengthssuchas305bpforSalleleand257bp forRallele.

2.3. Statisticalmethods

In this survey, the nomenclature of genotype and allele frequency,wefollowed,weredetailedbyPierce(2012).The pro-portionsweredeterminedwith95%confidenceinterval (CI95%) (Reiczigel et al., 2010). For comparison of genotype and allele frequencyintheruminantpopulations,chi-squaretestwith sig-nificancelevelof0.05wasperformedwithBonferronicorrection usingRstatisticalsoftwarei3863.3.0version.

3. Results

Inthisstudy,weexaminedaltogether210maleworms,which wereidentifiedasHaemonchuscontortus,todeterminethe propor-tionof Phe200TyrSNPoncodon200of-tubulingene isotype 1 within three different ruminant species. Among the worms, 110weregenotypedashomozygoussusceptible(SS),40as het-erozygous (RS) and 60 as homozygousresistant (RR), thus the totalfrequencyofSS,RSandRRwas52.4%(CI95%=45.5–59.1%), 19%(CI95%=14.1–24.9%)and28.6%(CI95%=22.7–35.2%), respec-tively. Considering thealleles,the proportionof susceptible (S) andresistant(R)allelewas61.9%(CI95%=57.2–66.5%)and38.1%

(CI95%=33.5–42.8%), respectively. Distribution of the different genotypewasshowedawidevarietyinhosts.TheSSwasthemost representativeinthereddeer,whereallofthewormsbelongedto thisgenotypegroup.TheRRwasmostprevalentinthesheepandit wasmoderateintheroedeer,whiletheRSwasobservedinequal proportioninthesetwohosts(Table1).Differenceofallele fre-quenciesbetweenthehostpopulationswasconfirmedsignificant bychi-squaretest.

4. Discussion

AlthoughseveralstudieswereconcernedwithARinsmall rumi-nants(Papadopoulosetal.,2012;Roseetal.,2015),ourknowledge isverydeficientinconnectionwithwildruminants(Chintoan-Uta etal.,2014;Mackintoshetal.,2014)WeinvestigatedtheBZ resis-tancestatusofH.contortus,collectedfromsheepandtwosympatric deerspeciesand ourresultsconfirmedconsiderabledivergence betweenthehosts. InHungarianveterinary practice,oneof the mostpreferreddruggroupsisBZs;therefore,theirusageinparasite managementisconsideredverygeneral.Ourresultssuggestthat thehighlevelofRalleleproportioninsheepinthisregionisderived fromthehightreatmentfrequencyandthelong-termunrotated drugapplication;asourpreviousworkconfirmedit(under publi-cation).Ontheotherhand,inHungary,theusageofanthelmintics andotherprescriptionmedicinesforhuntingpartiesarerestricted bythefoodsafetylaw,becausetheconsumptionofthe therapeu-ticdosecannotbeensuredandthewithdrawalperiodcannotbe controlledinwildlife.Inconsequence,thehuntingpartiesofthe studiedareahaveneverappliedanyanthelmintics.Basedonthese

Click to BUY NOW!

w.tr

acker-software.c w Click to BUY NOW!

.tracker-software.c

Fig.1.LocalizationofcollectedHaemonchuscontorusfromdifferenthosts(circleswithx=reddeer;triangles=roedeer;blackdots=sheep).(Forinterpretationofthe referencestocolourinthisfigurelegend,thereaderisreferredtothewebversionofthisarticle.)

Table1

Genotypeandallelefrequencies(CI95%)inHaemonchuscontortus(S=susceptible,R=resistant)collectedfromdifferentruminanthosts.

Host Genotypefrequency(%) Allelefrequency(%)

homozygous

facts,wesuggested,theresistantH.contortusrecoveredfromroe deermightderiveratherfromcrossinfectionfromsheep,thanfrom thefailureofhuntingmanagement.

OurobservationcouldnotexplainexactlytheabsenceofRallele inthereddeer.Wehypothesizedthatthepresenceofsheepon pasturesindicatedavoidancebehaviourinthereddeer.Somefree rangestudies(Osborne1984;DeGabrieletal.,2011)showedthe absenceof sheepin a given habitat induced a higher reddeer density.AccordingtoCuartasetal.(2000),thisphenomenonwas attachedtothedisturbanceeffectofsheephusbandryratherthan anydirectorindirectcompetitionbetweenthetwospecies.During ourinvestigation,manysheepflocksusedregularlythepastures ofthestudyarea,thustheirpresencemightbeaffectedthered deeraggregation;however,thisstudydidnotcoverdetailedfield observationofallhostspecies.TheresultofaNorwegianstudy (Mysterudetal.,1999)suggestedthatroedeerreservedits ter-ritorybya trade-offmechanismbetweenshelterandforage. In conclusions,theauthorsfoundthattherewasnoanyevidencefor avoidancebetweenthetworuminants.

Thehomerangeofroedeercanvarywidelyinscale(Jeppesen, 1990; Maillard et al., 2002) and it can be affected by habitat resources.Somestudiesconfirmthat theroedeeradjusthome rangesizeinresponsetolandscapestructure.Fragmented land-scapescanprovideabundant,high-qualityfood,whichincreases deerhabitatcarryingcapacity(Reimoser,2003;SaïdandServanty 2005).Thesevaluable,fragmentedlandscapescanalsobe utiliz-ableforsheepgrazing;asitwasobservedinthestudysite.We

supposedthattheoverlapinhabitatusebetweensheepflocksand roedeercouldcontributetotheoccurrenceandspreadofresistant allelewithinwildlife.Ourresultsraisedanissueabouttheroleof roedeerinAR;whetherroedeerisapotentialtransmitterofthe problemorthepotentialsolutionforit,asasourceofSSgenotyped wormssupportingrefugiaformationonpastures(VanWyk,2001).

Toanswerthisquestion,weneedtoaccomplishmoreregulated andcomplexstudiestoknowtheexacteffectsofthishostspecies.

Ourstudywasbasedonthemoleculardiagnosticdetectionof Phe200TyrSNPoncodon200of-tubulingeneisotype1;while SNPs on codon 167and 198were not investigated. Therefore;

thetrueprevalenceofBZresistancewithinthestudiedH. contor-tuspopulationmightbeunderestimated;notwithstanding,inall threehostspecies;thereby thosewerecomparable.Thesample collectionperiodwasratherlongbecausedeerforpost-mortem examinationcouldbereachedduringtheirhuntingseasonsand inalimitedcount.Thisfactcouldaffectourresultsthroughthe annualcircleofwormburdenincertainhosts.Sheepsamplescould betakenallyearround;therefore,thesedatarepresentmostlythe averageprevalenceofresistance.Roedeersampleswerecollected mainlyin spring and summerwhen theintensity of this para-siteisgenerallyratherhighinruminants.Ontheotherhand,red deerwereharvestedduringautumnandwinterwhenthelevelof infectionislargelymoderateorlow.Consideringthesefacts,our deerdatacannotdescribethestatementofresistanceaccurately.

Theremarkabledifferencebetweenthetwodeerspeciescanbe explainedevenbydifferentsamplingperiods.Roedeerwere

sam-Click to BUY NOW!

w.tr

acker-software.c w Click to BUY NOW!

.tracker-software.c

pledinspringandsummerwhenthesheepwithlambsattheirfeet areonpasturesandshedalotofeggs.Underthesecircumstances, wormsthatinfecttheroedeermightoriginatepartlyfromsheep.

Ontheotherhand,thesamplingofreddeerwascarriedoutafter themainwormburdenofsheep,whenfewlarvaecanbeingested bygrazing.During this period,theresearchercanmeeta natu-rallyselectedpopulationofworms;whichwereaffectedbyalot ofenvironmentalfactors,exceptanthelmintics,inawildruminant stomach.Consideringtheabove,ourresultsalsosupporttherefugia hypothesisifweassumethatabomasalcircumstancesofwild rumi-nantsevolveasortofselectionpressuretoparasites.Amongthese conditions,theselective-advantageattributedtoARmightprove tobeuseless;andthoseworms,whichcarryRallele,willperishin largeramountundertransferredselectionpressure.AsLeathwick (2013)confirmedinhisstudy,wormscarryinganthelmintic resis-tancemechanismarelikelytobelessfit.Forscientificjustification ofthishypothesis,furtherinvestigationswillbeneeded.

Conflictofinterest

Theauthorsdeclarenoconflictofinterest.

Acknowledgements

Theauthorsarevery thankfultotheprofessionalhunters of SEFAGPlc.fordeersamplecollection,totheseniorlecturerofthe KaposvárUniversity,GáborCsizmadiaforhisprofessional linguis-ticproposals;andtotheanonymousreviewersfortheirhelpful criticism.

References

Barrere,V.,Falzon,L.C.,Shakya,K.P.,Menzies,P.I.,Peregrine,A.S.,Prichard,R.K., 2013.AssessmentofbenzimidazoleresistanceinHaemonchuscontortusin sheepflocksinOntario,Canada:comparisonofdetectionmethodsfordrug resistance.Vet.Parasitol.198(1–2),159–165.

Chaudhry,U.,Redman,E.M.,Raman,M.,Gilleard,J.S.,2015.Geneticevidencefor thespreadofabenzimidazoleresistancemutationacrosssouthernIndiafrom asingleoriginintheparasiticnematodeHaemonchuscontortus.Int.J.Parasitol.

45(11),721–728.

Chaudhry,U.,Redman,E.M.,Ashraf,K.,Shabbir,M.Z.,Rashid,M.I.,Ashraf,S., Gilleard,J.S.,2016.MicrosatellitemarkeranalysisofHaemonchuscontortus populationsfromPakistansuggeststhatfrequentbenzimidazoledrug treatmentdoesnotresultinareductionofoverallgeneticdiversity.Parasit.

Vectors9,349,http://dx.doi.org/10.1186/s13071-016-1624-0.

Chintoan-Uta,C.,Morgan,E.R.,Skuce,P.J.,Coles,G.C.,2014.Wilddeeraspotential vectorsofanthelmintic-resistantabomasalnematodesbetweencattleand sheepfarms.Proc.Biol.Sci.281(1780),20132985.

Coles,G.C.,Jackson,F.,Pomroy,W.E.,Prichard,R.K.,vonSamson-Himmelstjerna, G.,Silvestre,A.,Taylor,M.A.,Vercruysse,J.,2006.Thedetectionofanthelmintic resistanceinnematodesofveterinaryimportance.Vet.Parasitol.136(3–4), 167–185.

Cuartas,P.,Gordon,I.J.,Hester,A.J.,Pérez-Barbería,F.J.,Hulbert,I.A.R.,2000.The effectofheatherfragmentationandmixedgrazingonthedietofsheepOvis ariesandreddeerCervuselaphus.ActaTheriol.45(3),309–320.

DeGabriel,J.L.,Albon,S.D.,Fielding,D.A.,Riach,D.J.,Westaway,S.,Irvine,R.J.,2011.

Thepresenceofsheepleadstoincreasesinplantdiversityandreductionsin theimpactofdeeronheather.J.Appl.Ecol.48,1269–1277.

FAO,2004.GuidelinesResistanceManagementandIntegratedParasiteControlin Ruminants.FoodandAgricultureOrganizationoftheUnitedNations(FAO), Romeftp://ftp.fao.org/docrep/fao/010/ag014e/ag014e06.pdf.

Jeppesen,J.L.,1990.(Homerangeandmovementsoffree-rangingroedeer (Capreoluscapreouls)atKalø.DanishRev.GameBiol.14(1),1–18.

Kaplan,R.M.,2004.Drugresistanceinnematodesofveterinaryimportance:a statusreport.TrendsParasitol.20(10),477–481.

Leathwick,D.M.,2013.Managinganthelminticresistance-Parasitefitness,druguse strategyandpotentialforreversiontowardssusceptibility.Vet.Parasitol.198, 145–153.

Lichtenfels,J.R.,Pilitt,P.A.,Hoberg,E.P.,1994.Newmorphologicalcharactersfor identifyingindividualspecimensofHaemonchusspp:(Nematoda:

Trichostrongyloidea)andakeytospeciesinruminantsinNorthAmerica.J.

Parasitol.80(1),107–119.

Mackintosh,C.G.,Cowie,C.,Fraser,K.,Johnstone,P.,Mason,P.C.,2014.Reduced efficacyofmoxidectinandabamectininyoungreddeer(Cervuselaphus)after 20yearsofmoxidectinpour-onuseonaNewZealanddeerfarm.Vet.Parasitol.

199(1–2),81–92.

Maillard,D.,Calenge,C.,Invernia,N.,Gaudin,J.C.,2002.Homerangesizeoffemale roedeerreintroducedintoaMediterraneanhabitat.Z.Jagdwiss.48,194–200.

Mylrea,G.E.,Mulley,R.C.,English,A.W.,1991.Gastrointestinalhelminthosisin fallowdeer(Damadama)andtheirresponsetotreatmentwithanthelmintics.

Aust.Vet.J.68(2),74–75.

Mysterud,A.,Larsen,P.K.,Ims,R.A.,Østbye,E.,1999.Habitatselectionbyroedeer andsheep:doeshabitatrankingreflectresourceavailability?Can.J.Zool.77, 776–783.

Nagy,G.,Csivincsik,Á.,Zsolnai,Á.,Sugár,L.,2016.Benzimidazoleresistancein Haemonchuscontortusrecoveredfromfarmedreddeer.Parasitol.Res.115(9), 3643–3647.

Osborne,B.C.,1984.Habitatuseofreddeer(CervuselaphusL.)andhillsheepin thewesthighlands.J.Appl.Ecol.21,497–506.

Papadopoulos,E.,Gallidis,E.,Ptochos,S.,2012.Anthelminticresistanceinsheepin Europe:aselectedreview.Vet.Parasitol.189,85–89.

Pe ˜na-Espinoza,M.,Thamsbort,Stig.M.,Demler,J.,Enemark,H.L.,2014.Field efficacyoffouranthelminticsandconfirmationofdrug-resistantnematodesby controlledefficacytestandpyrosequencingonasheepandgoatfarmin Denmark.Vet.Parasitol.206(3–4),208–215.

Pierce,B.A.,2012.Genetics.AConceptualApproach,fourthed.W.HFreemanand Company,NewYork.

Reiczigel,J.,Földi,J.,Ózsvári,L.,2010.Exactconfidencelimitsforprevalenceofa diseasewithanimperfectdiagnostictest.Epidemiol.Infect.138,1674–1678.

Reimoser,F.,2003.Steeringtheimpactsofungulatesontemperateforests.J.Natl.

Conserv.10,243–252.

Rose,H.,Rinaldi,L.,Mavrot,F.,deWaal,T.,Skuce,P.,Charlier,J.,Torgerson,P.R., Hertzberg,H.,Hendrickx,G.,Vercruysse,J.,Morgan,E.R.,2015.Widespread anthelminticresistanceinEuropeanfarmedruminants:asystematicreview.

Rose,H.,Rinaldi,L.,Mavrot,F.,deWaal,T.,Skuce,P.,Charlier,J.,Torgerson,P.R., Hertzberg,H.,Hendrickx,G.,Vercruysse,J.,Morgan,E.R.,2015.Widespread anthelminticresistanceinEuropeanfarmedruminants:asystematicreview.