parasitic infections which emerge in wildlife may have a significant effect on individuals and populations ( mller , 2005 ) .
this is especially a critical issue in the conservation of threatened species ( thompson et  al . , 2010 ) .
therefore , understanding the role of infectious and parasitic agents in wildlife populations , such as how they can influence fitness and in severe cases lead to population decline or extinction , requires accurate data on the diversity and load of potential pathogens in natural ecosystems , especially at the local scale ( smith et  al . , 2009 ) .
the appearance and spread of new parasites in populations of wild animals is not well documented .
most co - evolution studies are conducted with microbial organisms under laboratory conditions , taking advantage of their short generation times ( levin and lenski , 1983 , buckling and rainey , 2002 , gandon et  al . , 2008 ) .
however , such simplified artificial models do not reflect the situation in wild populations where infected animals are additionally affected by a number of environmental factors .
the difficulty in conducting such studies in wildlife results from both the difficulty in capturing the start point of the introduction of a new pathogen to a population and problems associated with long - term continuous acquisition of material for research .
the european bison ( bison bonasus l. , 1758 ) is one of the last mega - herbivores to survive in europe .
extinct in the wild at the beginning of 20th century , it was restored back to the wild from captive survivors .
nowadays there are over 3500 european bison in the wild distributed in 35 isolated populations ( raczyski , 2015 ) .
bison numbers are growing ; however , only five populations number more than 200 individuals . apart from small herd size , isolation and low genetic variation due to a severe genetic bottleneck after extinction in the wild ( tokarska et  al . , 2011 ) , diseases and parasites ( pucek et  al . , 2004 ) are the main threats to bison .
the oldest and largest population of bison is in the biaowiea primeval forest ( bpf ) ( 984 individuals , in both the polish and belarusian parts of bpf ) , constituting the core of the global population of the species . since reintroduction to the wild in 1952 , this population has been parasitologically investigated and monitored at differing intensities . a total of 88 species of parasites have been discovered in european bison , with an increasing trend in species richness as well as in the prevalence and intensity of infections ( drd , 1995 , karbowiak et  al .
this may be a result of transmission from other wild ruminants and cattle , or the increased contact rate between bison in winter due to large aggregations arising from supplementary feeding in fixed locations ( radwan et  al .
one of the most pathogenic parasites discovered , the blood - sucking nematode ashworthius sidemi schulz , 1933 was first found in bison in bpf in 2000 . the other blood - sucking nematode 
haemonchus contortus  was described in bison in bpf in the 1960s , but only in captivity ( drd , 1961 , drd , 1967 ) .
a.  sidemi was also found in the neighboring knyszyn forest in 2009 ( demiaszkiewicz et  al . , 2009b ) .
the species is characteristic of asiatic deer , especially sika deer cervus nippon , and was probably transmitted to poland from red deer cervus elaphus from ukraine and slovak republic , where sika deer has been previously introduced ( kotrl and kotrl , 1973 , kotrl and kotrl , 1973 , 1977 , drd et  al . , 1998 , demiaszkiewicz et  al . ,
the parasite was found in ne and se poland in species such as red deer , roe deer capreolus capreolus , moose alces alces , and european bison ( drd et  al . , 2003 ,
a.  sidemi was also found in deer in france ( ferte et  al . ,
, the presence of a.  sidemi was genetically confirmed in cattle , which indicates possible transmission from wildlife to livestock ( moskwa et  al . , 2015 ) .
previous studies from bpf have shown the influence of winter bison aggregation on a.  sidemi and coccidia infection ( radwan et  al . , 2010 , pyziel et  al . ,
2011 ) and on the seasonal pattern of parasite egg excretion ( koodziej - sobociska et  al . ,
the aim of our study was the long - term analysis of the invasive a.  sidemi spread since its appearance in the free ranging population of european bison in bpf and to determine the factors affecting its intensity .
we hypothesized that after the new , invasive parasite was introduced to the bison population , a phase of rapid growth in prevalence and intensity of a.  sidemi was due to the susceptibility of the animals to infection and a higher likelihood of disease transmission . in the bison population , we presume that a.  sidemi spread is significantly favoured by the high bison densities in large herds which form at fixed locations where supplementary winter feeding is provided , as well as biological factors such as age and sex .
knowledge about the spread of this parasite in wild hosts and particularly about factors which influence its prevalence and intensity may help to model the spread of emerging diseases in wildlife and guide the conservation management of wild endangered fauna .
bison were culled by rifle by biaowiea national park staff , as approved by the management plan .
the bison population in bpf has been regulated , on average , by approximately 11% annually since the 1970s through culling or translocation ( hayward et  al . , 2011 ,
biaowiea national park staff culled bison under permissions issued each year by the ministry of environment and the general directorate for environmental protection ( warsaw , poland ) .
tissue samples were collected from culled animals under the permission of the provincial wildlife conservator in the years 19992005 , and from 2006 under permissions issued each year by the ministry of environment , general directorate for environmental protection ( warsaw , poland ) and regional directorate for environmental protection ( bialystok , poland ) . for this research the guidelines of the polish nature conservation act adopted on october 16 , 1991 ( dz.u .
nr 99 , poz.1079 ) and on april 16 , 2004 ( dz.u . 2004 nr 92 poz .
the study was conducted in the biaowiea primeval forest ( bpf ; 5229-5237n , 2331-2421e ) located on the polish - belarusian border .
it is one of the best preserved lowland forests in europe inhabited by well - preserved communities of animals .
the polish part of the forest ( 600  km ) is covered mainly by deciduous and mixed tree stands ( 94% ) ; open habitats constitute the remaining 6% ( sokoowski , 2004 ) .
the climate of bpf is transitional between atlantic and continental type with clearly marked cold and warm seasons .
mean annual temperature between 2000 and 2015 was  + 8.0  c ; mean temperature of the coldest month , january , was  3.5  c , and the warmest , july , was  + 19.9  c .
the vegetative season lasts , on average , 215 days ( range : 198238 ) and was calculated according to tylkowski ( 2013 ) .
snow cover persisted from 41 to 120 days per year with a maximum recorded depth of 55  cm .
the studied bison population is widespread in the polish part of bpf across an area of 800  km ( kowalczyk et  al . , 2013 ) .
the population of 522 individuals , i.e. 15% of all wild living bison ( raczyski , 2015 ) , is isolated from the population in the belarussian part of the forest by a fence constructed on the border in 1981(kowalczyk et  al . , 2012 ) .
during snow - free periods bison males live solitarily ( 62% of males ) or in small bull groups ( up to 8 ind . ) , while females with calves and sub - adults roam in mixed groups numbering , on  average , 1115 individuals ( krasiska and krasiski , 1995 , krasiska et  al . , 2000 ) .
herd size increases during the rutting season ( august  october ) . in winter , bison are supplementary fed with different intensities and aggregate around seven main feeding sites .
supplementary fed bison stay in a limited area for up to five months creating herds of up to 100 individuals . during the first part of the winter they consume hay stored in roofed haystacks in fixed locations
; later in winter , the supplementary food ( usually hay ) is delivered to the forest floor in these locations .
some bison roam in winter out of the feeding sites , usually utilizing hay left by farmers on mowed meadows out of the forest or grazing on winter crops of rape and cereals ( hofman - kamiska and kowalczyk , 2012 , kowalczyk et  al . , 2011 , kowalczyk et  al . , 2013 ) .
there is relatively low number of cattle grazing in neighboring to the forest pastures , however around 15% of bison , occasionally or seasonally and mainly from autumn until early spring , forage on pastures utilized also by cattle ( kowalczyk et  al . , 2013 ) .
this indicates possibility of interspecies parasite transmission ( moskwa et  al . , 2015 ) .
no infected animals were recorded ( drd , 1961 , drd et  al . ,
to 2000 , 41 bison were investigated and no a.  sidemi was found ( drd et  al . , 1989 ,
in 20002001 , three of the 19 investigated bison were infected with a.  sidemi , but due to incomplete data on bison sex , age , and body weight or herd size , they were not included in the model , but were used for prevalence calculations . from 2002 to 2015 , 165 culled bison were investigated ( from 4 to 22 individuals per winter season ) ( see table  1 ) .
five individuals which had been released from captivity into the wild were excluded from the analysis due to their infection by a high number of nematodes ( max .
this may have resulted from their different susceptibility to infection , related to their lack of contact with the parasite in captivity as a result of regular deworming .
the contents of the whole of the abomasa and a 1  m long section of the duodena of culled bison were examined by scraping and sedimentation method ( drd et  al . , 1998 ,
the sediment preserved in 23% formalin was diluted with water up to 2  l and thoroughly mixed .
a 200  ml aliquot was examined in small portions under a dissecting microscope in order to collect all helminths .
after evaporation of the alcohol , non - permanent preparations were made from all selected nematodes for species identification .
, 1994 , jacquiet et  al . , 1997 and drd et  al .
( 1998 ) and the number of a.  sidemi were counted . for the analysis of infection intensity we used 160 individuals with a full set of data from 2002 to 2015 . to model which factors affected a.  sidemi infection intensity ( actual number of nematodes ) , we used the following set of independent variables : bison sex , age , body weight , herd size , feeding intensity , and the number of years since parasite introduction .
herd size was determined during the annual inventory of bison in bpf carried out by biaowiea national park staff .
examined bison were divided into three groups depending on the intensity of supplementary feeding provided in winter : intensively fed bison  food delivered to the feeding sites three to five times a week ; less intensively fed bison 
food delivered once a week , and non - fed bison  bison not utilizing feeding sites .
bison age was expressed in three age classes : juveniles ( up to 1 year old ) , sub - adults ( 23 years old ) and adults ( over 3 years old ) .
we checked for multi - collinearity among explanatory variables , and the following non - correlated ( r  <  0.5 ) most informative covariates were included in the final model : bison sex , age , body weight , herd size and number of years since introduction .
then , we fitted negative binomial generalized linear models ( negative binomial glm ) to count data , which dealt with observed model over - dispersion ( zuur et  al . , 2009 ) .
the number of a.  sidemi per animal was set as a dependent variable and as main effects , we applied bison age , sex , body mass , and as interactive effects : number of years since introduction with herd size and body mass .
the akaike information criterion ( aic ) with second - order correction for small sample size ( aicc ; burnham and anderson , 2002 ) was used for model ranking .
beforehand , out of the set of all possible combinations of models , based on model nesting , we excluded more complex versions of models which had lower aicc scores ( richards , 2008 , richards et  al . , 2010 ) .
we checked the normality and homoscedasticity in the distribution of the final model residuals by inspecting the quantile - quantile distribution plot and model residuals against plots of fitted values ( estimated responses ) .
all statistical analyses were performed in the r program ( version 3.1.2 ; r development core team , 2012 ) .
prevalence of a.  sidemi infection increased rapidly from 0% in 1999 , by 9% in 2000 , 25% in 2001 , reaching 100% in 2004 ( fig .  1 ) .
sidemi invasion , the median infection intensity was 2890 and varied from 0 to 44,310 nematodes per animal  ( fig .  2 ) .
after its introduction in 2000 , infection intensity increased , reaching the highest median infection intensity of 8200 a.  sidemi in the winter of 2008/2009 .
afterwards , a significant decline in the median infection intensity was observed to the minimum value of 410 nematodes in the winter of 2014/2015 ( nbglm : slope  =  0.22 + -0.09 , z  =  2.45 , p  =  0.01 ) ( fig .  2 ) .
in 20022006 , bison herd size , from which studied individuals originated ranged from 4 to 104 individuals and did not differ significantly across years .
then , from 2007 , herd size significantly decreased ( 2007  max herd size  =  75 individuals , 2015  =  31 individuals ; p  < 
the size of the bison herds was strongly correlated with the intensity of supplementary winter feeding ( fig . 
the size of the bison herds utilizing provided forage grew as the frequency and amount of food delivered to feeding sites increased .
model selection ( based on the aicc criteria ) for the considered negative binomial generalized linear model ( glm ) revealed that changes in a.  sidemi infection intensity were statistically significantly associated with bison age , and the interactive effects between herd size and the number of years since introduction .
a.  sidemi infection intensity was the highest in sub - adults ( median : 10,835 nematodes ) , medium in adults ( median : 2955 ) and the lowest in juveniles ( median : 1700 ) ( fig .  5 ) .
this was also confirmed by a loess model ( local polynomial regression )  infection intensity grew up to the age of approximately 3842 months , then it dropped and after 75 months remained at a fairly stable level ( fig .  6 ) .
the positive effect of bison herd size on a.  sidemi infection intensity increased significantly with the number of years since nematode introduction ( table  3 , fig .  7 ) .
though the final model indicated that bison sex did not have a significant effect on a.  sidemi infection intensity , we found significant differences between sexes when only adult individuals were considered ( males  1785 nematodes , females  3982 ) ( fig .  8 , mann - whitney test , w  =  954 , p  =  0.027 ) .
in this paper we analysed the spread of a pathogenic parasite in a wild ungulate from the initial phase of its progression .
this was possible due to the long - term parasitological studies conducted on european bison in bpf ; this enabled the capture of the initial phase of the invasive blood - sucking nematode introduction , and the investigation of its spread over the following years .
such studies are of great importance because they allow for modelling the spread of a new parasite in a population , and its impact on the host population ( anderson and may , 1978 , hudson et  al . , 1998 ) and ecological processes ( thomas et  al . , 1999 ,
long term field studies of host and parasite populations represent the ideal approach to analyse the dynamics and consequences of co - evolution because they allow the direct detection of host and parasite reciprocal influences ( decaestecker et  al . , 2007 ) .
however , the full course of new parasite introduction and spread in wild vertebrates is difficult to accurately trace .
recently , host - parasite interactions in the extreme environment of the arctic have been reviewed by kutz et  al .
the authors , based on studies of parasites of holarctic ruminants such as moose ( alces alces ) , muskoxen ( ovibos moschatus ) and caribou ( rangifer tarandus ) , concluded that insights gained from studying the history and ecology of host  parasite systems in the arctic will be central to understanding the role that climate change is playing in these complex systems ( hoberg et  al .
, 2007 , kutz et  al . , 2012 , hoberg et  al . , 2012 , steele et  al . , 2013 ,
it can also guide wildlife management and conservation throughout the arctic , and be generalized to provide insights into host 
parasites influence host populations causing deleterious effects on their hosts ( anderson and may , 1979 , prado et  al . , 2009 ) ,
and so they are constantly forced to adapt to one another ( schmid - hempel , 2011 ) .
short phylogenetic parasite - host relationships are known to be more pathogenic for the host than longer co - existing relationships ( gandon et  al . , 2008 ) .
our study revealed that ten years after a.  sidemi introduction , its infection intensity dropped significantly and prevalence varied between 89 and 100% .
the question is whether several years of interaction with the pathogen adaptive mechanisms to limit the number of parasites was developed ?
most co - evolution studies are conducted with microbial organisms under laboratory conditions taking advantage of their short generation times ( levin and lenski , 1983 , buckling and rainey , 2002 , gandon et  al . , 2008 ) , and these studies are n't comparable with mega - herbivore species such as european bison .
long - term co - evolution sometimes leads to a relatively stable relationship tending towards commensalism or mutualism and the establishment of parasite - host balance , as it is in the evolutionary interest of the parasite that its host thrives ( rook , 2007 ) .
although we have not investigated the immunological parameters , it could be expected , based on other host - parasite interaction studies , that co - evolution could have appeared , or will appear in a.  sidemi - bison interactions and affect the course of invasion . due to the lack of data on this topic
european bison are characterized by very low genetic variation ( the lowland line of the european bison originates from only seven founders ) ; this raises the questions about the relationship of genetic diversity to the infection of the species by numerous new parasites and the increasing prevalence and intensity of these infections ( tokarska et  al . , 2009 , tokarska et  al . , 2011 ) . there is considerable evidence that at least part of the natural variation in resistance to nematode infection is under genetic control ( wakelin , 1985 , barger , 1989 , stear et  al . , 1999 ) .
however , analysis of the impact of genetic factors ( the major histocompatibility complex ) on the intensity of a.  sidemi infection in bison has shown that despite high sequence divergence , neither the alleles nor drb heterozygosity were significantly associated with infection intensity ( radwan et  al . , 2007 ,
, the risk associated with the occurrence of any new pathogen in a population of a species with such a limited gene pool should be taken into account .
blood - sucking parasites , such as a.  sidemi , are strongly pathogenic ( osiska et  al . , 2010 ) , especially as a cause of anaemia ( stefaski , 1968 , gatongi et  al . , 1998 ) .
a study by boughton and hardy ( 1935 ) revealed that 2000 h.  contortus females can suck 30  ml of blood per day ; thus the high infection intensity of a.  sidemi ( reaching 44,300 nematodes ) observed in bison in bpf may cause significant blood loss resulting in a serious impact on host condition ( koodziej - sobociska et  al .
moreover , the histopathological examination of tissue from infected bison showed infiltrations of inflammatory cells in the walls of the abomasa and duodena ( mainly lymphoid cells and eosinophils , as well as hyperaemia , oedema and lesions of mucosa and proliferation of lymphatic follicles ) at various levels of intensity ( osiska et  al . , 2010 ) .
such serious histopathological lesions connected with a deterioration in blood parameters in infected animals ( koodziej - sobociska et  al .
, 2016b ) confirm the parasite potentially has a strong impact on the bison population .
there are many mechanisms and factors involved in parasitism : the immune system , toxins , behavioural traits , population density and others ( solomon et  al . , 2015 ) .
one of the major host defence mechanisms is immunity , which is influenced by factors such as age , sex and nutritional status ( bush et  al .
, 2001 , pea et  al . , 2004 , cornell et  al . , 2008 ) .
the higher a.  sidemi infection intensity in sub - adult bison compared to juveniles could be due to longer exposure to the parasite .
difference between results obtained from young animals to these from adults is probably connected to an immature immune system .
young animals , with immature immune systems , are the most susceptible to infection and disease and they spread a high number of infective eggs into the environment in their faeces ( woolhouse , 1998 , cornell et  al . , 2008 ,
neonates and juvenile ruminants are very susceptible to paratuberculosis infection , and that this is related to a high degree of exposure from their mothers and an immature immune system ( thakur et  al . , 2013 ) . also in captive bison , it was found that young bison are usually characterized by the heaviest infection levels and suffer the most severe consequences ( treboganova , 2010 ) . in bpf , the youngest bison ( up to one year old ) have a lower parasitic load than sub - adults .
we hypothesize that this could be a result of both the short exposure time of the host to the parasite , and maternal antibodies ingested by bison calves which provide immunological protection ( jackson and nazar , 2006 , hurley and theil , 2011 ) .
the suckling period in bison lasts up to one year ( krasiska and krasiski , 2013 ) . however , most calves are born in may 
july , and by winter are already a few months old and eat hay provided at feeding sites ( krasiska and krasiski , 2013 ) which exposes them to infectious agents in the contaminated environment ( radwan et  al . , 2010 ,
our detailed analysis of infection dynamics in relation to age showed that infection intensity grew with age up to approximately 3842 months , then dropped , and after 75 months , when bison are fully mature , remained at a fairly stable level . by adulthood animals
have developed a stronger immunity and usually harbour lower infection levels ( brass and stevens , 1982 , kethineni et  al . ,
the lower infection intensity that we observed in bison bulls stands in opposition to the widely observed patterns of the higher susceptibility of males to parasitosis ( addis , 1946 , solomon , 1966 , alexander and stimson , 1988 , zuk and mckean , 1996 ) .
this phenomenon is explained by the association between testosterone and the immune system ; thus , sexually mature male vertebrates are often more susceptible to infection and carry higher parasite burdens in the wild ( zuk and mckean , 1996 ) .
recent studies have revealed that male - biased parasitism is not universal and that there are many other factors that can influence parasite infection such as the age of individuals , sexual size dimorphism , hormone levels , individual host variability ( e.g. behavioural , physiological ) , and immunocompetence ( kiffner et  al . , 2013 ) .
we suggest that the behavioural traits may be responsible for the lower parasitic load in bulls .
bison males live solitarily ( 62% of males ) or in small bull groups ( up to 8 individuals ) , while females with calves and sub - adults roam in groups numbering , on average , 1115 individuals ( krasiska and krasiski , 1995 , krasiska et  al . , 2000 ) .
this leads to sexual segregation and much lower contact rates between males and other bison , which decreases the risk of parasite transmission .
additionally , pregnant and lactating females are immunosuppressed and therefore more susceptible to infection ( lloyd , 1983 , krishnan et  al . , 1996 ) .
european bison have an average gestation period of 264 ( range 254270 ) days ( krasiska and krasiski , 2013 ) and similarly to closely related species such as american bison ( bison bison ) , can give birth each year ( wilson et  al . , 2002 ) . also , during a four year study 60% of cows gave birth to at least 3 calves .
moreover , the infection intensity is not an effect of just one season , but may be influenced by pregnancies in previous years .
these data support the higher a.  sidemi infection intensity observed in bison females in bpf .
we found that herd size , as a factor which significantly shapes a.  sidemi infection intensity in bison , became increasingly important the longer the parasite was present in the host population .
this means that the highest infection intensities were recorded in the biggest bison herds , where the winter supplementary feeding of bison was intensive .
the role of herd size increased , the longer the time since the first detection of the parasite .
this indicates that during the initial phase , the new parasite spread more successfully even in small herds , probably due to the high susceptibility of the host to the new pathogen .
after the host population became saturated with the parasite , other factors such as herd size increasingly shaped the dynamics of the disease .
in large herbivore management , supplementary feeding is mainly intended to reduce herbivore impact on agriculture , to enhance body condition and reproductive performance ( kozak et  al . , 1995 ) , or to supply endangered animals with food or water in crucial periods of their annual life cycle ( loarie et  al . , 2009 ) .
most of the free - ranging bison populations inhabiting forests are supplementary fed in winter to mitigate migrations and reduce farm crop depredation ( kerley et  al .
however , research conducted in the last few years has showed that supplementary feeding may have negative long - term effects on bison and increase parasitic load ( radwan et  al . , 2010 ,
this is a consequence of bison aggregation in fixed locations for several months of the year and may lead to environmental contamination by parasites and their more efficient transmission ( radwan et  al . , 2010 ) . increasing intensity of supplementary
feeding leads to higher aggregation and larger herds and reduces bison mobility and ranging ( schneider , 2008 ) .
this may be more severe during cold winters when bison occupy much smaller winter ranges ( krasiska et  al . , 2000 ) .
an annual survey of parasite excretion in bison dung showed an increase in their prevalence and number during winter months ( koodziej - sobociska et  al . , 2016a ) .
increased winter densities of bison at feeding sites influences the parasitic load of this herbivore , especially in intensively fed herds ( radwan et  al . , 2010 ,
our results have shown that the level of a.  sidemi infection is dependent on the herd size and that the herd size is strongly correlated with the intensity of supplementary feeding .
this in turn may shape the demographic composition and the overall fitness of the european bison populations ( hayward et  al .
management practices also influence the quick spread and growth of infection intensity in bison , especially at the end of winter when bison are in a poorer condition ( hayward et  al . , 2015 , koodziej - sobociska et  al .
the results of this study indicate a strong need to modify supplementary feeding in order to scatter the herds , increase bison ranges and as a consequence reduce parasite transmission .
actions implemented in bpf during the bison conservation program in 20062010 aimed to split large winter aggregations of bison ( numbering initially up to 104 individuals ) , which resulted in an increase in the number of winter herds and a reduction in their size .
, we would predict increased parasitic load strongly related to the herd size and supplementary feeding intensity .
our study , for the first time , has allowed for the long - term analysis of parasite invasion in a large herbivore in the wild : from the appearance of the parasite , through the period of rapid increase in prevalence and intensity of the infection , to the stabilization phase in the host - parasite system .
the spread of animals into new areas , due to climate change and environmental modification , can affect the occurrence and spread of parasitic infections to new areas and new species .
knowledge about the course of biological invasions and the factors influencing them in wildlife can help in the modelling , prediction , prevention , or even reduction of the emergent diseases .
this is especially important in the conservation management of rare and endangered species . in recent years ,
bison have been introduced to new places in europe ; therefore , knowledge about the processes of emergent diseases and the determination of factors influencing the severity of pathogen and parasite infections can help to design conservation management that reduces the risk of infection spread and therefore the consequences for populations of this unique and other wild species .
 the full course of new parasite introductions in wild animals is difficult to accurately trace . 
 we documented and analysed the invasive blood - sucking nematode ashworthius sidemi ( trichostrongylidae ) introduction and spread in european bison ( bison bonasus ) from the initial phase of its progression . in the polish part of the biaowiea primeval forest ( bpf ) the parasite was first found in 2000 . from 2002 to 2015 , 165 culled bison were investigated . the prevalence and intensity of a.  sidemi schulz , 
 1933 infection increased over the following years , reaching 100% of investigated bison four years after introduction and a maximal median intensity of 8200 nematodes per animal in the winter of 2008/2009 . 
 afterwards , a significant decline of median infection intensity was observed to the minimum value of 410 nematodes per animal . between 2011 and 2014 
 prevalence varied from 89 to 100% . among the factors analysed , the number of years since introduction , herd size , age and sex proved to significantly influence infection intensity . 
 a higher infection intensity was recorded in sub - adults compared to juveniles and adults . 
 males had significantly lower infection intensity than females , but this was the case for adults only . 
 the highest infection intensities were recorded in the biggest bison herds , where the winter supplementary feeding of bison is intense . moreover , 
 the longer the parasite was present in the host population , the more important herd size became as a factor . 
 our study indicates that it is not solely biological factors that determine the spread of a newly detected parasite in wildlife , but that management practices can also have a strong influence . 
 this is especially important in endangered species under intensive human care as the management practices may pose a threat to the species . 
Introduction
Materials and methods
Results
Discussion
Conflict of interest