- Open Access
Spatiotemporal interactions between wild boar and cattle: implications for cross-species disease transmission
Veterinary Research volume 45, Article number: 122 (2014)
Controlling infectious diseases at the wildlife/livestock interface is often difficult because the ecological processes driving transmission between wildlife reservoirs and sympatric livestock populations are poorly understood. Thus, assessing how animals use their environment and how this affects interspecific interactions is an important factor in determining the local risk for disease transmission and maintenance. We used data from concurrently monitored GPS-collared domestic cattle and wild boar (Sus scrofa) to assess spatiotemporal interactions and associated implications for bovine tuberculosis (TB) transmission in a complex ecological and epidemiological system, Doñana National Park (DNP, South Spain). We found that fine-scale spatial overlap of cattle and wild boar was seasonally high in some habitats. In general, spatial interactions between the two species were highest in the marsh-shrub ecotone and at permanent water sources, whereas shrub-woodlands and seasonal grass-marshlands were areas with lower predicted relative interactions. Wild boar and cattle generally used different resources during winter and spring in DNP. Conversely, limited differences in resource selection during summer and autumn, when food and water availability were limiting, resulted in negligible spatial segregation and thus probably high encounter rates. The spatial gradient in potential overlap between the two species across DNP corresponded well with the spatial variation in the observed incidence of TB in cattle and prevalence of TB in wild boar. We suggest that the marsh-shrub ecotone and permanent water sources act as important points of TB transmission in our system, particularly during summer and autumn. Targeted management actions are suggested to reduce potential interactions between cattle and wild boar in order to prevent disease transmission and design effective control strategies.
Most pathogens of concern to livestock are able to cross-infect multiple host species, including wildlife, and therefore in areas where wildlife and livestock co-occur (i.e. interface areas), pathogens can emerge and establish in these sympatric host populations . For example, foot and mouth disease, rabies, anthrax, brucellosis and bovine tuberculosis (TB) have all been shown to be reciprocally transmissible between livestock and wildlife –. In this context, the demography and behaviour of the hosts’ populations can play an important role in intra- and interspecific pathogen transmission by determining contact rates and environmental exposure . If resources that are commonly used by both domestic and wild species are aggregated, this can result in high spatial and/or temporal overlap between two or more species –, further increasing the probability of disease transmission. How habitat use by hosts affects direct and indirect interactions among hosts is fundamental in understanding multi-host disease transmission , and is critical for designing scientifically-based disease control strategies . Nonetheless, the role that spatial and temporal interactions between livestock and wildlife play in exposure to pathogens and disease transmission remains mostly unknown ,.
Tuberculosis caused by the Mycobacterium tuberculosis complex is an important re-emerging zoonotic disease shared between domestic cattle and wildlife, and the control of this disease is largely limited by the existence of wildlife reservoirs –. In the United Kingdom, for instance, cattle may become infected with TB by using farm buildings (feed stores and cattle sheds) and grazing on grass that has been contaminated with urine, faeces, sputum or wound exudates of badgers (Meles meles) ,. In the United States, white-tailed deer (Odocoileus virginianus) and cattle often share rangeland resources, including water sources and feeders, although temporal segregation between these species is often observed . In the Iberian Peninsula, wild boar (Sus scrofa) are the main wild maintenance host of TB . Recent studies from Spain suggest that TB infection can spread not only by direct contact among individuals but also by indirect transmission , with water sources being high risk areas where pathogen transmission can occur between wildlife and cattle through consumption of short-term infected water . However, epidemiological studies at the interface between livestock and important disease-carrying wildlife, such as wild boar, remain scarce.
Despite compulsory testing and culling of infected cattle, TB infection rates in cattle populations are persistently high in Doñana National Park (DNP), southern Spain . The TB-host community of DNP includes wild boar, red deer (Cervus elaphus) and fallow deer (Dama dama), all of which occur sympatrically in areas used for traditional cattle husbandry. Interestingly, the populations of these three wildlife hosts exhibit common spatial patterns of TB infection across DNP, which may be explained by resource use and behaviour of these species ,. Recent advances in global positioning system (GPS) technology for monitoring wildlife has proven useful for assessing fine-scale spatiotemporal interactions among species , and thus may provide a fundamental understanding of the risk of TB transmission at the wildlife/livestock interface. We deployed GPS technology on cattle and wild boar in DNP to test the hypothesis that patterns of resource selection and spatiotemporal overlap between these two species increase the local risk of interspecific disease transmission. Specifically, we aimed to determine where and when the activity patterns of cattle and wild boar overlapped and whether areas with the greatest potential overlap corresponded with areas with high incidence of TB in cattle.
Material and methods
We conducted the study in DNP (37°0′ N, 6°30′ W), a protected nature reserve located on the Atlantic coast of southern Spain (Figure 1). The region has a Mediterranean climate, classified as dry sub-humid with marked seasons. In the wet season (December–May), marshlands are flooded and ungulates graze in elevated shrublands. Ungulates in DNP are mostly food limited during summer (June–September), when wetlands and natural water bodies dry up causing senescence of herbaceous vegetation. However, a north–south humid ecotone habitat exists year-round between the elevated shrublands and the low dry marshlands; vegetation within this ecotone is dominated by Scirpus maritimus and Galiopalustris sp. with Juncus maritimus associations (see Additional file 1 for further details on habitat types).
The study area has moderate to high densities of red deer, fallow deer and wild boar throughout DNP. A traditional breed of cattle (locally called “marismeña”) is farmed within five cattle management areas in DNP (Figure 1). Coto del Rey (CR) is the northern border of DNP and contains no cattle husbandry. The central area includes three cattle enclosures: SO (n = 350 cattle; density = 5.7 cattle/km2), BR (n = 168 cattle; density = 2.6 cattle/km2), and PU (n = 152 cattle; density =4.0 cattle/km2). Marismillas (MA) is the southern-most area (n = 318 cattle; density = 3.1 cattle/km2). Each cattle management area is surrounded by a cattle-proof fence, which limits the movements of each herd to within their designated management area. However, social groups (overwhelmingly females) showing individual ranging behaviour may be differentiated within each cattle management area . The incidence of TB in cattle is high within DNP (9.23% per year on average), and TB prevalence in wild boar (45–52%), red deer, and fallow deer (14–19%) populations is also high ,. DNP has been proposed as a natural scenario for describing the epidemiology of shared diseases in wild and domestic ungulates ,,.
Animal capture and monitoring
We used data from 18 wild boar and 12 head of cattle from the marismeña breed that were equipped with GPS radio-collars between July 2011 and October 2013. Animal capture followed a protocol approved by the Animal Experiment Committee of Castilla-La Mancha University and by the Spanish Ethics Committee, and designed and developed by scientists (B and C animal experimentation categories) in accordance with EC Directive 86/609/EEC for animal handling and experiments. We captured wild boar using six padded foothold cage traps monitored using camera traps (see  for further details). Captured wild boar were anaesthetized, weighed, ear tagged, radio-collared and assessed for condition, age and sex. The anaesthetic protocol (3 mg/kg of tiletamine-zolazepam and 0.05 mg/kg of medetomidine) followed Barasona et al. . Of the collared wild boar, 11 were males (3 subadults, < 24 months; 8 adults) and seven were females (2 subadults; 5 adults). All collared cattle were adult females. We captured wild boar in different trapping areas across DNP in order to collar a sample of animals from multiple social groups. Cattle from different social groups were radio-collared during routine veterinary inspections of cattle restrained in the farm’s cattle yards. Although cattle were collared in several cattle management areas, the intense trapping efforts were carried out in BR cattle management area, where both species were concurrently monitored (Figure 1).
Radiocollars were programmed to acquire one GPS location per hour and to transmit accumulated packets of 20 locations using GSM (Microsensory System, Spain) . Data collected included date, time, geographic coordinates, and location acquisition time (LAT, which is a measure of the precision of a fix and ranges between 0–160 s). First, we screened GPS locations with LAT ≥ 154 s to detect anomalous fixes (manufacturer’s technical data; Microsensory System, Spain). Using this criterion, 189 and 66 GPS fixes were considered anomalous and thus removed from wild boar and cattle databases, respectively. We also discarded GPS locations obtained during the day of collar deployment and of collar retrieval to avoid possible anomalous behaviour associated with handling procedures, even though differences in behaviour post-handling were not detected elsewhere . Positional error associated with GPS locations averaged 26.6 m (SD = 23.5 m), based on stationary tests from 19 collars (1637 locations in total) carried out in the centre of our study area (i.e., open sky). Fix-rate success averaged 81.2% and 94.0% for wild boar and cattle, respectively. We explored whether the lower fix-rate success obtained for wild boar introduced habitat-induced biases  or dial-induced biases (e.g., due to wild boar using dense vegetation as rest sites during the day ,. However, no significant differences were found in mean LAT values among habitats (Kruskal-Wallis test, z = 48.00, p > 0.05) or between day and night (z = −1.88, p > 0.05). Consequently, we did not correct for habitat-induced fix-rate bias.
Coarse-scale spatial overlap between wild boar and cattle
We estimated annual and seasonal (winter, spring, summer and autumn) home-ranges (HR; 95% Utilization Distribution, UD; ) and core-areas (CA; 50% UD) used by each collared animal using the fixed-kernel function from the ADEHABITAT package  in R version 2.15.2 . Kernels were estimated using the reference bandwidth method  because the least-squares cross-validation method failed to converge for six animals with large sample sizes . Fixed-kernel density estimators allow identification of disjunct areas of activity , which can be particularly important for assessing interspecific patterns of space use in heterogeneous environments like DNP.
Home ranges and CA for each individual animal were used to estimate annual and seasonal spatial overlap  between wild boar and cattle within the BR cattle management area (Figure 1), where both species were concurrently monitored. Spatial overlap was calculated as the area of overlap in HR or CA between wild boar and cattle divided by (1) the total area of HR or CA for wild boar (i.e., overlap for wild boar relative to cattle), or (2) the total area of HR or CA for cattle (i.e., overlap for cattle relative to wild boar).
Fine-scale spatial interaction between wild boar and cattle
The extent of overlap in HR and CA provides only a coarse indicator of the potential interactions between two species because HR and CA estimators represent only the outline of a distribution of locations . To assess annual and seasonal fine-scale interactions and differences in the use of available resources between cattle and wild boar, we estimated latent selection difference functions (LSDs) ,. The GPS locations were transformed into 26 m radius circular buffers (to account for GPS positional error) , and within each buffer we calculated: straight-line distance (km) to nearest artificial water hole (DW); straight-line distance (km) to nearest marsh–shrub ecotone (DE); proportional cover of dense scrub (LT1); proportional cover of low-clear shrubland (LT2); proportional cover of herbaceous grassland (LT3); proportional cover of woodland (LT4); proportional cover of bare land (LT5); and proportional cover of watercourse vegetation (LT6; see Additional file 1). These predictor variables were selected because of their biological relevance for explaining ungulate distribution in DNP (see ). Landcover data was obtained from Andalusia Environmental Information . Collinearity among predictor variables was screened using a Spearman’s pairwise correlation coefficient value of |r| > 0.5 .
We estimated LSDs using logistic regression  and the “RMS” package  in R. For this analysis, we coded locations from cattle as 1 and those from wild boar as 0, i.e. we assessed cattle resource selection or avoidance relative to wild boar. In this analysis, landcover variables with significant positive coefficients indicate those most preferred by cattle relative to wild boar, whereas those with significant negative coefficients indicate those most avoided by cattle relative to wild boar. Distance to-variables, however, should be interpreted the opposite way. Variables with non-significant coefficients represent those habitats with the highest potential for interspecific interactions, because there is no difference in the use or selection of these resources between the two species. The results from LSD analyses can then be used to make inferences about the differences or similarities in fine-scale habitat use and spatial overlap between the two species. The main assumption of LSDs is that all resources should be equally available to both species within the study area. To fulfill this assumption, we only used locations from collared animals located in the central cattle management area (BR; 13 wild boar and 10 cattle; Figure 1) that occurred within annual and seasonal inter-species CA (50% UD) overlap contours, i.e. the area where the greatest inter-species interactions could occur. To account for an unbalanced sampling design and non-independence of observations from the same individual, we estimated robust standard errors using the Huber–White sandwich estimator , grouping data by individual.
We randomly split the annual and seasonal datasets, using 70% of locations to parameterize the models (training datasets) and the remaining 30% of locations for model validation (validation datasets) . The best annual and seasonal models were obtained using a forwards–backwards stepwise procedure on the training datasets based on Akaike Information Criteria (AIC) . We assessed predictive capacity of the best annual and seasonal models using calibration plots. Calibration plots were constructed by testing the annual or seasonal best models on the corresponding validation dataset, and then plotting the observed and predicted frequency of observations in each of 10 equal-size intervals of predicted probabilities (0–1). A model with high predictive capacity should show perfectly aligned points along a 45° line (see ). We also assessed the predictive capacity of each model with the area under the receiver operating characteristic curve (AUC), to rate the ability of the models to correctly discriminate between cattle and wild boar locations. The AUC ranges from 0.5 for models with no discrimination ability to 1 for models with perfect discrimination ability .
Annual and seasonal best LSD models were used to spatially map the relative probability of use (P) by cattle relative to wild boar. Areas with values of P of approx. 0.5 were considered as those where the highest relative probability of spatial interaction between both species could occur . Accordingly, we constructed a spatial interspecific interaction (SII) index for the whole of BR using the rule: SII = (1 – P) if P ≥ 0.5, and SII = P if P < 0.5. Further, because models estimating resource use by a species can be used to predict the species’ distribution in other geographical areas (e.g. ), we used models trained with data from BR to extrapolate P and derived SII index within 1 ha cells across the whole of DNP. Predicted SII values across DNP were correlated with TB epidemiological data (see below).
Sampling and TB diagnosis
Between 2006 and 2013, 570 wild boar were opportunistically shot by park rangers in DNP, and necropsied as part of the DNP health-monitoring programme (see  for details). We recorded the location where animals were shot, the year in which they were sampled, and their TB lesion score, gender and age. Necropsies were performed by qualified wildlife veterinarians that had extensive experience in the diagnosis of macroscopic TB-compatible lesions. Veterinarians performed detailed inspections of the entire animal, including lymph nodes and abdominal and thoracic organs . Cultures using pyruvate-enriched Löwenstein-Jensen medium were performed to confirm TB infection. During the same time period and as part of the TB control programme in DNP, cattle populations within the SO, BR, PU and MA management areas were tested for TB by veterinary authorities using skin tests, and slaughtered if found positive. Prevalence and incidence (as used for cattle because the entire population was tested during annual sanitary campaigns) of TB were estimated for each cattle management area for wild boar and cattle populations, respectively. Finally, we assessed whether there were significant differences in the SII index (annually as well as seasonally) among cattle management areas with high and low TB-incidence in cattle using Mann–Whitney U-tests. All statistical analyses were performed in R version 2.15.2 .
We collected 44 699 locations from wild boar and 47 213 locations from cattle during the study period. Collared wild boar were distributed across all five cattle management areas, whereas collared cattle were only present in BR and MA (Figure 1). The GPS locations were homogeneously collected throughout the study period for all seasons and for both species (see Additional file 2).
There was a stark contrast between estimated annual HR and CA sizes for cattle and wild boar (Figure 2), with cattle using significantly larger areas (average ± SE, HR = 1787.78 ± 826 ha; CA = 346.24 ± 174 ha) than wild boar (HR = 551.33 ± 260 ha; CA = 86 ± 77 ha), (ANOVA, F1, 28 = 16.57 for HR; F1, 28 = 15.21 for CA; both p < 0.001). There were significant seasonal differences in HR sizes for cattle (F3, 8 = 3.69, p = 0.023), but not for wild boar (F3, 14 = 2.47, p > 0.05) (Figure 3). The percent overlap in HR and CA between cattle and wild boar varied among seasons, with percent overlap being highest in autumn and lowest in winter (Table 1). Overall, > 60% of wild boar HR overlapped areas used by cattle, whereas ≤ 40% of the HR of cattle overlapped areas used by wild boar. Wild boar CA showed high overlap with areas used by cattle in spring, summer and autumn (66–78% overlap) but not in winter (only 23%).
Fine-scale assessment of spatial interactions between wild boar and cattle revealed that the environmental variables explaining relative habitat selection by cattle and wild boar differed among seasons (Table 2). During winter and spring, cattle used areas significantly further from water sources (DW) than wild boar; however, use of water sources did not differ significantly between the two species during summer and autumn. Cattle and wild boar selection for themarsh–shrub ecotone (DE) did not differ in any of the seasons analyzed. Conversely, cattle showed consistent avoidance of areas with a higher proportion of dense scrub (LT1) relative to wild boar across all seasons. Relative to wild boar, cattle showed avoidance of areas with higher proportions of low-clear shrubland (LT2), herbaceous grassland (LT3) and watercourse vegetation (LT6) during winter and spring. However, during summer and autumn, cattle and wild boar did not differ in their use of these three habitats. Annually, cattle and wild boar did not differ in their selection for areas close to water sources (DW) or for the marsh–shrub ecotone (DE). Validation of LSD models showed that all had good discriminatory power (all AUC > 0.7) and predictive reliability (Additional file 3), supporting their use in extrapolating spatial patterns of SII across the whole of DNP (Figure 4; see also Additional file 4 for annual SII). In general, areas with a high probability of use by both species (high potential interaction) were mostly associated with the marsh-shrub ecotone and permanent water sources (Table 2 and Additional file 4), especially during summer and autumn (dark areas in Figure 4), whereas shrub-woodlands and temporal grass-marshlands had a low probability of interaction.
TB and spatial overlap
Based on culture confirmed lesions, infection was detected in 55.7% (SE = 4.1%; n = 570) of wild boar tested. The prevalence of TB in wild boar was 45.9% (SE = 3.8%; n = 174) in MA, 64.7% (SE = 5.8%; n = 68) in PU, 46.6% (SE = 3.7%; n = 174) in BR, 73.9% (SE = 4.6%; n = 92) in SO, and 72.6% (SE = 5.7%; n = 62) in CR. Official skin testing of 1,139 cattle in DNP from 2006 to 2013 revealed a mean incidence of 9.0% TB reactors (SE = 4.9%). Mean prevalence of TB in wild boar differed significantly among cattle management areas when these areas were grouped into low (MA = 4.1% and BR = 5.6% TB-incidence in cattle; averagewild boar TB-prevalence = 46.3%) or high (PU = 18.1% and SO = 11.8% TB-incidence in cattle; averagewild boar TB-prevalence = 69.3%) TB-incidence in cattle (F1, 2 = 24.96; p < 0.05). Interestingly, the mean predicted value of annual SII (fine-scale spatial interspecific interaction) was also significantly higher in high TB-incidence areas than in low TB-incidence areas (Z = 88; p < 0.05; Figure 5). These differences were significant in spring, summer and autumn but not in winter (Figure 5).
We assessed fine-scale spatiotemporal interactions between wild and domestic hosts of TB in order to better understand what role resource selection may play in cross-species disease transmission. To our knowledge, this is the first study that has conducted a fine-scale spatial analysis aimed at explaining the patterns of disease transmission at the wild boar/cattle interface. We found that similar use of water resources by cattle and wild boar resulted in high potential interspecific interaction around these landscape features, especially during the dry season. This high spatial overlap at such small spatial extents (e.g. waterholes are only 15 m in diameter) could influence interspecific transmission rates of TB in this Mediterranean system. Our research contributes to an applied understanding of multi-host disease ecology and will help to better target actions and implement control strategies for TB at the wildlife/livestock interface.
We found that cattle had larger HR and CA than wild boar, indicating that these two species have different space use requirements as well as ranging behaviours. Large-scale ranging behaviour in “marismeña” cattle within DNP is mostly determined by human decisions rather than by species-specific traits . Each cattle management area contains a free-ranging cattle herd which is controlled and regulated according to the Cattle Use Plan . The mean HR and CA for wild boar recorded in our study area are somewhat lower than those reported in previous studies in other Mediterranean areas ,. Differences may be related to food availability, population density, activity behaviour and/or composition of social groups ,. Additional research would be required to elucidate which factors regulate ungulate spatial behaviour in DNP. The spatial distribution of cattle (using 95% or 50% UD) overlapped a large proportion of wild boar HR (97%) and CA (63%), and cattle frequently traversed areas used intensively by wild boar. Given that wild boar HR and CA were comparatively smaller than those of cattle, the concentration of important resources and cattle activity within areas commonly used by wild boar created situations that facilitated interaction (or at least spatial overlap) between the two species. Additionally, fine-scale spatial analyses suggest that within areas intensively used by wild boar there was limited spatial interspecific avoidance (also see ,).
Although there were some similarities in the patterns of resource use in all seasons, wild boar and cattle generally used different resources during winter and spring in DNP. Conversely, limited differences in resource selection during summer and autumn resulted in negligible spatial segregation, and thus probably high encounter rates, between the two species during these seasons. This suggests that interspecific contact and subsequent disease transmission between cattle and wild boar is likely to be highest in drier seasons. Interestingly, the spatial distribution of the interspecific interaction index was consistently high at the marsh–shrub ecotone during all seasons. This is probably because this heterogenous habitat offers important resources for ungulates throughout the year, such food and shelter (also see ). Furthermore, the predicted increase in spatial overlap and fine-scale interactions between wild boar and cattle during summer and autumn is likely related to the increased use of areas where forage and water are still available, when seasonal drought severely reduces the availability of resources in Mediterranean areas .
Characterizing and quantifying the potential interactions and the likelihood of disease transmission between domestic and wild hosts is crucial to understanding the complex dynamics of multi-host systems ,. The results from our assessment of the spatial ecology and interactions of wild boar and cattle suggest that environmental and/or interspecific behavioural factors could favour disease transmission at the livestock/wildlife interface. We found that spatial variation in the incidence of TB in cattle in DNP was positively associated with the prevalence of TB in necropsied wild boar, which is consistent with the hypothesis that TB transmission occurs among ungulates, as has previously been argued from both field and molecular epidemiology ,. In the case of wild boar, the high disease prevalence based on culture (up to 50%) observed in DNP is remarkable and indicative of a high risk of disease transmission ,, with about one third of pigs in a random sample expected to be actively excreting mycobacteria by several routes (mainly oro-nasally) .
The epidemiological interaction between the two host species described above was further supported by the fact that areas with high TB-incidence in cattle were also the areas with higher predicted spatial interaction between cattle and wild boar. This suggests that the dynamics of disease transmission in DNP are partly driven by the presence of environmental features that facilitate spatiotemporal overlap between hosts, as indicated by LSD models. The humid marsh–shrub ecotone and the surrounding water holes were the habitats with the highest potential interaction between wild boar and cattle. These landscape features may act as potential sources of M. tuberculosis complex for the host community ,,. For instance, Kukielka et al.  showed that shared water resources in South Central Spain were risky points where TB transmission could occur by indirect contact. Interestingly, we found that the predicted spatial interspecific interaction was highest in areas with high TB-incidence in cattle (and high TB-prevalence in wild boar) during summer and autumn, i.e. the time of year when species are most water-limited. These complex epidemiological scenarios have also been described in dry areas from Africa where cattle share water holes and diseases with wildlife . In South Spain, a recent study reviewed the environmental persistence of M. tuberculosis complex and found that wildlife/livestock interactions occur much more often at water sources than would be expected by chance alone . Aggregation of ungulates is promoted around water points, and this subsequently enhances the opportunities for transmission of diseases ,. This may arise because ungulates come into contact with either a higher proportion of individuals from the same or different species, or with a more heavily contaminated environment (i.e. direct and/or indirect mycobacteria transmission). Spread of TB may occur indirectly from contaminated vegetation, water, mud or fomites . Wild boar activity around these water sources (such as wallowing, brushing, drinking, defecating, urinating, and mating) is likely to result in environmental contamination and TB transmission to other hosts.
We used GPS telemetry data from concurrently monitored domestic cattle and wild boar to describe spatiotemporal interactions by means of new analytical procedures , and from these inferred associated implications for TB transmission. Within this framework, we considered fine-scale spatial overlap in habitats selected by cattle and wild boar as a proxy of interspecific contact. Although we did not measure contact directly, the difficulty of estimating realistic frequency of contact between species, most of which are predicted to be indirect, has been highlighted previously (e.g. ,). However, recent studies have attempted to measure interspecific contact rates in relation to the dynamics of disease transmission. For example, contact rates have been estimated by direct observation of domestic and wild animals in open habitats where they are easily observed, such as alpine meadows (e.g. ). However, this approach was not feasible in our study area because visibility is impeded by closed scrub. Other recent studies using telemetry data have defined critical time and space windows between pairs of GPS locations, and thus only assumed that interspecific contact had occurred within this critical window . Approaches based on proximity loggers potentially have the ability to estimate contact rates between individuals often to within a few meters; however, the performance of these devices is often poor, providing data that is only indicative of contact rates rather than actual contact rates where interactions occur . Further, within an epidemiological context, their utility is constrained to direct rather than indirect disease transmission. The LSD modelling procedure  we used proved a reliable tool to estimate annual and seasonal similarities in the use of shared resources, which is valuable for the study of diseases for which direct as well as indirect interactions among sympatric species are of importance in transmission dynamics, as our case . However, the approach was limited in that we could not demonstrate that the spatial overlap between cattle and wild boar occurred within a sufficiently fine-scale temporal window to be directly related to the transmission of TB. Despite this limitation, the LSD approach can provide spatial predictions which can be extrapolated to a larger area where hypotheses related to the spatial risk of interspecific disease transmission can be tested. Additionally, future research could use a combination of proximity loggers and GPS technology to validate rates of interspecific contact, quantify the potential for indirect disease transmission, and identify habitats where both these events occur most frequently.
Epidemiologists and policy makers need to understand the complex interspecific interactions among potential hosts to identify risk factors for disease transmission and prescribe targeted management actions . Our results highlight aspects of the hosts’ ecology and behaviour that are likely to affect the probability of interspecific disease transmission. Further, our results identify factors that need to be considered in order to prevent interactions between wild and domestic ungulates at key disease reservoir sources, such as permanent water sources in ecosystems with marked dry seasons. Although welfare of wild animals must be considered, it may be possible to segregate livestock and wild ungulates in areas surrounding permanent water sources. For example, farm biosecurity measures, like small-scale fencing, exclusion gates or deterrents , could be implemented at points such as water sources to prevent wild ungulate access to these areas. Recent innovations in South Spain showed that effective segregation strategies of wild ungulates at water points have the potential to reduce interspecific contact and TB transmission at the wildlife/cattle interface . Furthermore, research is being conducted currently into field vaccination of wild boar against M. bovis using oral baits . Ideally, tools from several fields of study should be combined into integrated control plans to minimize pathogen transmission  and to improve the cost-effectiveness of strategies such as host population control through random or selective culling or through habitat management .
Gortázar C, Ferroglio E, Höfle U, Frölich K, Vicente J: Diseases shared between wildlife and livestock: a European perspective. Eur J Wildl Res. 2007, 53: 241-256. 10.1007/s10344-007-0098-y.
Frölich K, Thiede S, Kozikowski T, Jakob W: A review of mutual transmission of important infectious diseases between livestock and wildlife in Europe. Ann N Y Acad Sci. 2002, 969: 4-13. 10.1111/j.1749-6632.2002.tb04343.x.
Artois M: Wildlife infectious disease control in Europe. J Mt Ecol. 2003, 7: 89-97.
Ward AI, Tolhurst BA, Delahay RJ: Farm husbandry and the risks of disease transmission between wild and domestic mammals: a brief review focusing on bovine tuberculosis in badgers and cattle. Anim Sci. 2006, 82: 767-773. 10.1017/ASC2006102.
Cooper SM, Scott HM, de la Garza GR, Deck AL, Cathey JC: Distribution and interspecies contact of feral swine and cattle on rangeland in south Texas: implications for disease transmission. J Wildl Dis. 2010, 46: 152-164. 10.7589/0090-3558-46.1.152.
Proffitt KM, Gude JA, Hamlin KL, Garrott RA, Cunningham JA, Grigg JL: Elk distribution and spatial overlap with livestock during the brucellosis transmission risk period. J Appl Ecol. 2011, 48: 471-478. 10.1111/j.1365-2664.2010.01928.x.
Dobson A: Population dynamics of pathogens with multiple host species. Am Nat. 2004, 164: 64-78. 10.1086/424681.
Nunn CL, Thrall PH, Kappeler PM: Shared resources and disease dynamics in spatially structured populations. Ecol Model. 2014, 272: 198-207. 10.1016/j.ecolmodel.2013.10.004.
Cooper SM, Perotto-Baldivieso HL, Owens MK, Meek MG, Figueroa-Pagan M: Distribution and interaction of white-tailed deer and cattle in a semi-arid grazing system. Agricult Ecosys Environ. 2008, 127: 85-92. 10.1016/j.agee.2008.03.004.
Hudson PJ, Rizzoli AP, Grenfell BT, Heesterbeek JAP, Dobson AP: Ecology of Wildlife Diseases. 2002, Oxford University Press, Oxford
Böhm M, Hutchings M, White P: Contact networks in a wildlife-livestock host community: identifying high-risk individuals in the transmission of bovine TB among badgers and cattle. PLoS One. 2009, 4: e5016-10.1371/journal.pone.0005016.
Martin C, Pastoret PP, Brochier B, Humblet MF, Saegerman C: A survey of the transmission of infectious diseases/infections between wild and domestic ungulates in Europe. Vet Res. 2011, 42: 70-10.1186/1297-9716-42-70.
Phillips C, Foster C, Morris P, Teverson R, Foster C, Morris P, Teverson R: The transmission of Mycobacterium bovis infection to cattle. Res Vet Sci. 2003, 74: 1-15. 10.1016/S0034-5288(02)00145-5.
Gortazar C, Delahay RJ, Mcdonald RA, Boadella M, Wilson GJ, Gavier‐Widen D, Acevedo P: The status of tuberculosis in European wild mammals. Mammal Rev. 2012, 42: 193-206. 10.1111/j.1365-2907.2011.00191.x.
Martínez-López B, Barasona JA, Gortázar C, Rodríguez-Prieto V, Sánchez-Vizcaíno JM, Vicente J: Farm-level risk factors for the occurrence, new infection or persistence of tuberculosis in cattle herds from South-Central Spain. Prev Vet Med. 2014, 116: 268-278. 10.1016/j.prevetmed.2013.11.002.
Garnett BT, Delahay RJ, Roper TJ: Use of cattle farm resources by badgers (Meles meles) and risk of bovine tuberculosis (Mycobacterium bovis) transmission to cattle. Proc Biol Sci. 2002, 269: 1487-1491. 10.1098/rspb.2002.2072.
Delahay RJ, Smith GC, Ward AI, Cheeseman CL: Options for the management of bovine tuberculosis transmission from badgers (Meles meles) to cattle: evidence from a long-term study. Mamm Study. 2005, 30 (Suppl 1): S73-S81. 10.3106/1348-6160(2005)30[S73:OFTMOB]2.0.CO;2.
Naranjo V, Gortazar C, Vicente J, de la Fuente J: Evidence of the role of European wild boar as a reservoir of Mycobacterium tuberculosis complex. Vet Microbiol. 2008, 127: 1-9. 10.1016/j.vetmic.2007.10.002.
Kukielka E, Barasona JA, Cowie CE, Drewe JA, Gortazar C, Cotarelo I, Vicente J: Spatial and temporal interactions between livestock and wildlife in South Central Spain assessed by camera traps. Prev Vet Med. 2013, 112: 213-221. 10.1016/j.prevetmed.2013.08.008.
Barasona JA, VerCauteren KC, Saklou N, Gortazar C, Vicente J: Effectiveness of cattle operated bump gates and exclusion fences in preventing ungulate multi-host sanitary interaction. Prev Vet Med. 2013, 111: 42-50. 10.1016/j.prevetmed.2013.03.009.
Gortazar C, Torres MJ, Acevedo P, Aznar J, Negro JJ, de la Fuente J, Vicente J: Fine-tuning the space, time, and host distribution of mycobacteria in wildlife. BMC Microbiol. 2011, 11: 27-10.1186/1471-2180-11-27.
Gortázar C, Torres MJ, Vicente J, Acevedo P, Reglero M, de la Fuente J, Negro JJ, Aznar-Martín J: Bovine tuberculosis in Doñana biosphere reserve: the role of wild ungulates as disease reservoirs in the last Iberian lynx strongholds. PLoS One. 2008, 3: e2776-10.1371/journal.pone.0002776.
Latham ADM, Latham MC, Anderson DP, Cruz J, Herries D, Hebblewhite M: The GPS craze: six questions to address before deciding to deploy GPS technology on wildlife.New Zeal J Ecol, in press.,
Lazo A: Ranging behaviour of feral cattle (Bos taurus) in Donana National Park, SW Spain. J Zool. 1995, 236: 359-369. 10.1111/j.1469-7998.1995.tb02718.x.
Romero B, Aranaz A, Sandoval Á, Álvarez J, de Juan L, Bezos J, Sánchez C, Galka M, Fernández P, Mateos A, Domínguez L: Persistence and molecular evolution of Mycobacterium bovis population from cattle and wildlife in Doñana National Park revealed by genotype variation. Vet Microbiol. 2008, 132: 87-95. 10.1016/j.vetmic.2008.04.032.
Barasona JA, López-Olvera JR, Beltrán-Beck B, Gortázar C, Vicente J: Trap-effectiveness and response to tiletamine-zolazepam and medetomidine anaesthesia in Eurasian wild boar captured with cage and corral traps. BMC Vet Res. 2013, 9: 107-10.1186/1746-6148-9-107.
Cano-Manuel J, Granados JE, Castillo A, Serrano E, Pérez JM, Soriguer RC, Fandos P, Travesí R: Nuevas tecnologías aplicadas al seguimiento de ungulados silvestres en Sierra Nevada: Collares GPS-GSM. Biodiversidad y Conservación de Fauna y Flora en Ambientes Mediterráneos. Edited by: Barea JM, Ballesteros E, Luzón JM, Moleón M, Tierno JM. 2007, Sociedad Granatense de Historia Natural Sierra Nevada, Granada, Spain, 691-705.
Frair JL, Nielsen SE, Merrill EH, Lele SR, Boyce MS, Munro RH, Stenhouse GB, Beyer HL: Removing GPS collar bias in habitat selection studies. J Appl Ecol. 2004, 41: 201-212. 10.1111/j.0021-8901.2004.00902.x.
Lemel J, Truvé J, Söderberg B: Variation in ranging and activity behaviour of European wild boar Sus scrofa in Sweden. Wildl Biol. 2003, 9 (Suppl 1): 29-36.
DeCesare NJ, Squires JR, Kolbe JA: Effect of forest canopy on GPS-based movement data. Wildl Soc Bull. 2005, 33: 935-941. 10.2193/0091-7648(2005)33[935:EOFCOG]2.0.CO;2.
Sodeikat G, Pohlmeyer K: Escape movements of family groups of wild boar Sus scrofa influenced by drive hunts in Lower Saxony. Germany Wildl Biol. 2003, 9 (Suppl 1): 43-49.
Calenge C: The package adehabitat for the R software: a tool for the analysis of space and habitat use by animals. Ecol Model. 2006, 197: 516-519. 10.1016/j.ecolmodel.2006.03.017.
R Development Core Team R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing; ., [http://www.R-project.org]
Seaman DE, Millspaugh JJ, Kernohan BJ, Brundige GC, Raedeke KJ, Gitzen RA: Effects of sample size on kernel home range estimates. J Wildl Manage. 1999, 63: 739-747. 10.2307/3802664.
Hemson G, Johnson P, South A, Kenward R, Ripley R, Macdonald D: Are kernels the mustard? Data from global positioning system (GPS) collars suggests problems for kernel home‐range analyses with least‐squares cross‐validation. J Anim Ecol. 2005, 74: 455-463. 10.1111/j.1365-2656.2005.00944.x.
Fieberg J, Kochanny CO: Quantifying home-range overlap: the importance of the utilization distribution. J Wildl Manage. 2005, 69: 1346-1359. 10.2193/0022-541X(2005)69[1346:QHOTIO]2.0.CO;2.
Latham ADM, Latham MC, Boyce MS: Habitat selection and spatial relationships of black bears (Ursus americanus) with woodland caribou (Rangifer tarandus caribou) in northeastern Alberta. Can J Zool. 2011, 89: 267-277. 10.1139/z10-115.
Peters W, Hebblewhite M, DeCesare N, Cagnacci F, Musiani M: Resource separation analysis with moose indicates threats to caribou in human altered landscapes. Ecography. 2013, 36: 487-498. 10.1111/j.1600-0587.2012.07733.x.
Recio MR, Mathieu R, Denys P, Sirguey P, Seddon PJ: Lightweight GPS-tags, one giant leap for wildlife tracking?An assessment approach. PLoS One. 2011, 6: e28225-10.1371/journal.pone.0028225.
Braza F, Alvarez F: Habitat use by red deer and fallow deer in Doñana National Park. Miscellània Zoològica. 1987, 11: 363-367.
Consejería de Medio Ambiente y Ordenación del Territorio, Andalucía, España: Red de Información Ambiental de Andalucía, REDIAM. ., [http://www.juntadeandalucia.es/medioambiente/site/rediam]
Tabachnick BG, Fidell LS: Using Multivariate Statistics. 1996, HarperCollins, New York, USA
Hosmer DW, Lemeshow S: Applied Logistic Regression. 2000, John Wiley & Sons, New York, USA
Frank E, Harrell J: rms: Regression Modeling Strategies. R Package Version 4.0-0; 2013. ., [http://CRAN.R-project.org/package=rms]
Freedman DA: On the so-called “Huber sandwich estimator” and “robust standard errors”. Am Stat. 2006, 60: 299-302. 10.1198/000313006X152207.
Boyce MS, Vernier PR, Nielsen SE, Schmiegelow FKA: Evaluating resource selection functions. Ecol Model. 2002, 157: 281-300. 10.1016/S0304-3800(02)00200-4.
Akaike H: A new look at the statistical model identification. IEEE Trans Automat Control. 1974, 19: 716-723. 10.1109/TAC.1974.1100705.
Pearce J, Ferrier S: Evaluating the predictive performance of habitat models developed using logistic regression. Ecol Model. 2000, 133: 225-245. 10.1016/S0304-3800(00)00322-7.
Acevedo P, González-Quirós P, Prieto JM, Etherington TR, Gortázar C, Balseiro A: Generalizing and transferring spatial models: a case study to predict Eurasian badger abundance in Atlantic Spain. Ecol Model. 2014, 275: 1-8. 10.1016/j.ecolmodel.2013.12.011.
Vicente J, Höfle U, Garrido JM, Fernández-De-Mera IG, Juste R, Barral M, Gortazar C: Wild boar and red deer display high prevalences of tuberculosis-like lesions in Spain. Vet Res. 2006, 37: 107-119. 10.1051/vetres:2005044.
Espacio Natural Doñana: Spain: Plan de Aprovechamiento Ganadero Del Parque Nacional de Doñana; 2000. ., [http://www-rbd.ebd.csic.es/gestion/territorioyrecursos/ganaderia/PlanGanadero.pdf]
Boitani L, Mattei L, Nonis D, Corsi F: Spatial and activity patterns of wild boars in Tuscany, Italy. J Mammal. 1994, 75: 600-612. 10.2307/1382507.
Massei G, Genov PV, Staines BW, Gorman ML: Factors influencing home range and activity of wild boar (Sus scrofa) in a Mediterranean coastal area. J Zool. 1997, 242: 411-423. 10.1111/j.1469-7998.1997.tb03845.x.
Bugalho MN, Milne JA: The composition of the diet of red deer (Cervus elaphus) in a Mediterranean environment: a case of summer nutritional constraint?. Forest Ecol Manag. 2003, 181: 23-29. 10.1016/S0378-1127(03)00125-7.
Rodríguez-Prieto V, Martínez-López B, Barasona J, Acevedo P, Romero B, Rodriguez-Campos S, Gortázar C, Sánchez-Vizcaíno J, Vicente J: A Bayesian approach to study the risk variables for tuberculosis occurrence in domestic and wild ungulates in South Central Spain. BMC Vet Res. 2012, 8: 148-10.1186/1746-6148-8-148.
Barasona JA, Torres MJ, Armenteros JA, Diez-Delgado I, Gortázar C, Vicente J: Environmental presence of M. bovis at aggregation points in the wildlife/livestock interface in Mediterranean areas. Proceedings of the VI International M. bovis Conference: 16–19 June 2014; Cardiff, Wales. Edited by: Hewinson Glyn. 2014, British Cattle Veterinary Association, Weybridge, 68-
Fine AE, Bolin CA, Gardiner JC, Kaneene JB: A study of the persistence of Mycobacterium bovis in the environment under natural weather conditions in Michigan. USA Vet Med Int. 2011, 2011: 765430-
Munyeme M, Muma JB, Skjerve E, Nambota AM, Phiri IGK, Samui KL, Dorny P, Tryland M: Risk factors associated with bovine tuberculosis in traditional cattle of the livestock/wildlife interface areas in the Kafue basin of Zambia. Prev Vet Med. 2008, 85: 317-328. 10.1016/j.prevetmed.2008.03.006.
Vicente J, Höfle U, Garrido JM, Acevedo P, Juste R, Barral M, Gortazar C: Risk factors associated with the prevalence of tuberculosis-like lesions in fenced wild boar and red deer in south central Spain. Vet Res. 2007, 38: 451-464. 10.1051/vetres:2007002.
Vicente J, Barasona JA, Acevedo P, Ruiz‐Fons JF, Boadella M, Diez‐Delgado I, Beltran-Beck B, González-Barrio D, Queirós J, Montoro V, de la Fuente J, Gortazar C: Temporal trend of tuberculosis in wild ungulates from Mediterranean Spain. Transbound Emerg Dis. 2013, 60: 92-103. 10.1111/tbed.12167.
Morris RS, Pfeiffer DU, Jackson R: The epidemiology of M. bovis infections. Vet Microbiol. 1994, 40: 153-177. 10.1016/0378-1135(94)90053-1.
Miguel E, Grosbois V, Caron A, Boulinier T, Fritz H, Cornélis D, de Garine-Wichatitsky M: Contacts and foot and mouth disease transmission from wild to domestic bovines in Africa. Ecosphere. 2013, 4: 51-10.1890/ES12-00239.1.
Richomme C, Gauthier D, Fromont E: Contact rates and exposure to inter-species disease transmission in mountain ungulates. Epidemiol Infect. 2006, 134: 21-30. 10.1017/S0950268805004693.
Drewe JA, Weber N, Carter SP, Bearhop S, Harrison XA, Dall SR, McDonald RA, Delahay RJ: Performance of proximity loggers in recording intra-and inter-species interactions: a laboratory and field-based validation study. PLoS One. 2012, 7: e39068-10.1371/journal.pone.0039068.
Nishi JS, Shury T, Elkin BT: Wildlife reservoirs for bovine tuberculosis (Mycobacterium bovis) in Canada: strategies for management and research. Vet Microbiol. 2006, 112: 325-338. 10.1016/j.vetmic.2005.11.013.
Reidy MM, Campbell TA, Hewitt DG: Evaluation of electric fencing to inhibit feral pig movements. J Wildl Manage. 2008, 72: 1012-1018. 10.2193/2007-158.
Beltrán-Beck B, Romero B, Sevilla IA, Barasona JA, Garrido JM, González-Barrio D, Díez-Delgado I, Minguijón E, Casal C, Vicente J, Gortázar C, Aranaz A: Assessment of an oral Mycobacterium bovis BCG vaccine and an inactivated M. bovis preparation for wild boar in terms of adverse reactions, vaccine strain survival, and uptake by nontarget species. Clin Vaccine Immunol. 2014, 21: 12-20. 10.1128/CVI.00488-13.
Boadella M, Vicente J, Ruiz-Fons F, de la Fuente J, Gortazar C: Effects of culling Eurasian wild boar on the prevalence of Mycobacterium bovis and Aujeszky’s disease virus. Prev Vet Med. 2012, 107: 214-221. 10.1016/j.prevetmed.2012.06.001.
Authors would like to acknowledge many students and collaborators who contributed to field work over the study period, also wish to express their gratitude to the Doñana National Park, the monitoring team of EBD-CSIC for their help with the fieldwork and ICTS of Doñana Biological Station. The present work has benefited from the financial aid of research grants JCCM (PEII10-0262-7673), EU (FP7 grant 613779 WildTBVac) and MINECO (AGL2013-48523-C3-1-R). JAB holds an FPU pre-doctoral scholarship. PA is supported by MINECO-UCLM through “Ramón y Cajal” contract (RYC-2012-11970) and partly by EMIDA-ERA-NET grant APHAEA (219235-FP7-ERA-NET-EMIDA). We also thank Graham Nugent (Landcare Research, Lincoln, New Zealand) for his valuable comments on the study.
The authors declare that they have no competing interests.
JAB, JV contributed to the conception, design, data collection, laboratory work, data analysis, drafting and writing of the manuscript. MCL, PA, ADML contributed to design, data analysis and drafting of the manuscript. JAA, CG, FC, RCS participated in the data collection and drafting of the manuscript. All authors have read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Covariates used in the spatial analysis. Environmental predictors, descriptions, mean values (M) and standard deviations (SD) of GPS locations buffers versus total study area grids used in the analysis of resource separation patterns between cattle and wild boar at Doñana National Park. (DOCX 13 KB)
Additional file 3: Calibration plots of the predictive performance of the models. Assessment of the predictive performance of the best seasonal and annual models (see Table 2). Each plot shows the relationship between the predicted probability to be used by cattle in relation to wild boar and the observed proportion of cattle locations on the validation dataset. (DOCX 212 KB)
Additional file 4: Pattern of annual interspecific interaction. Spatial gradient in predicted annual interspecific interaction index (0 = low interaction, 0.5 = maximum interaction) between domestic cattle and wild boar in Doñana National Park, Spain, July 2011–October 2013. Predicted probability of interaction between the two species was derived from an annual Latent Selection Difference model (see Table 2). (DOCX 383 KB)
About this article
Cite this article
Barasona, J.A., Latham, M.C., Acevedo, P. et al. Spatiotemporal interactions between wild boar and cattle: implications for cross-species disease transmission. Vet Res 45, 122 (2014) doi:10.1186/s13567-014-0122-7
- Global Position System
- Wild Boar
- Resource Selection
- Fallow Deer
- Proportional Cover