Risk of tuberculosis cattle herd breakdowns in Ireland: effects of badger culling effort, density and historic large-scale interventions

Inferring bovine tuberculosis transmission between cattle and badgers via the environment and risk mapping

Bovine tuberculosis (bTB), caused by Mycobacterium bovis, is one of the most challenging and persistent health issues in many countries worldwide. In several countries, bTB control is complicated due to the presence of wildlife reservoirs of infection, i.e. European badger (Meles meles) in Ireland and the UK, which can transmit infection to cattle. However, a quantitative understanding of the role of cattle and badgers in bTB transmission is elusive, especially where there is spatial variation in relative density between badgers and cattle. Moreover, as these two species have infrequent direct contact, environmental transmission is likely to play a role, but the quantitative importance of the environment has not been assessed. Therefore, the objective of this study is to better understand bTB transmission between cattle and badgers via the environment in a spatially explicit context and to identify high-risk areas. We developed an environmental transmission model that incorporates both within-herd/territory transmission and between-species transmission, with the latter facilitated by badger territories overlapping with herd areas. Model parameters such as transmission rate parameters and the decay rate parameter of M. bovis were estimated by maximum likelihood estimation using infection data from badgers and cattle collected during a 4-year badger vaccination trial. Our estimation showed that the environment can play an important role in the transmission of bTB, with a half-life of M. bovis in the environment of around 177 days. Based on the estimated transmission rate parameters, we calculate the basic reproduction ratio (R) within a herd, which reveals how relative badger density dictates transmission. In addition, we simulated transmission in each small local area to generate a first between-herd R map that identifies high-risk areas.

Impact of test, vaccinate or remove protocol on home ranges and nightly movements of badgers a medium density population

In the British Isles, the European badger (Meles meles) is thought to be the primary wildlife reservoir of bovine tuberculosis (bTB), an endemic disease in cattle. Test, vaccinate or remove ('TVR') of bTB test-positive badgers, has been suggested to be a potentially useful protocol to reduce bTB incidence in cattle. However, the practice of removing or culling badgers is controversial both for ethical reasons and because there is no consistent observed effect on bTB levels in cattle. While removing badgers reduces population density, it may also result in disruption of their social behaviour, increase their ranging, and lead to greater intra- and inter-species bTB transmission. This effect has been recorded in high badger density areas, such as in southwest England. However, little is known about how TVR affects the behaviour and movement of badgers within a medium density population, such as those that occur in Northern Ireland (NI), which the current study aimed to examine. During 2014-2017, badger ranging behaviours were examined prior to and during a TVR protocol in NI. Nightly distances travelled by 38 individuals were determined using Global Positioning System (GPS) measurements of animal tracks and GPS-enhanced dead-reckoned tracks. The latter was calculated using GPS, tri-axial accelerometer and tri-axial magnetometer data loggers attached to animals. Home range and core home range size were measured using 95% and 50% autocorrelated kernel density estimates, respectively, based on location fixes. TVR was not associated with measured increases in either distances travelled per night (mean = 3.31 ± 2.64 km) or home range size (95% mean = 1.56 ± 0.62 km², 50% mean = 0.39 ± 0.62 km²) over the four years of study. However, following trapping, mean distances travelled per night increased by up to 44% eight days post capture. Findings differ from those observed in higher density badger populations in England, in which badger ranging increased following culling. Whilst we did not assess behaviours of individual badgers, possible reasons why no differences in home range size were observed include higher inherent 'social fluidity' in Irish populations whereby movements are less restricted by habitat saturation and/or that the numbers removed did not reach a threshold that might induce increases in ranging behaviour. Nevertheless, short-term behavioural disruption from trapping was observed, which led to significant increases in the movements of individual animals within their home range. Whether or not TVR may alter badger behaviours remains to be seen, but it would be better to utilise solutions such as oral vaccination of badgers and/or cattle as well as increased biosecurity to limit bTB transmission, which may be less likely to cause interference and thereby reduce the likelihood of bTB transmission.

Can Ingoing Contact Chains and other cattle movement network metrics help predict herd-level bovine tuberculosis in Irish cattle herds?

We used logistic regression to investigate whether the risk of an Irish cattle herd undergoing a bovine tuberculosis (bTB) breakdown increased with the size of the Ingoing Contact Chain (ICC) of previous herd to herd cattle movements, in a sequence up to eight moves back from the most recent, direct, movement into the herd. We further examined whether taking into account the bTB test history of each herd in the chain would improve model fit. We found that measures of cattle movements directly into the herd were risk factors for subsequent bTB restrictions, and the number of herds that animals were coming from was the most important of these. However, in contrast to a previous study in Great Britain, the ICC herd count at steps more remote than direct movements into the herd did not result in better fitting models than restricting the count to direct movements. Restricting the ICC counts to herds which had previously or would in the future test positive for bTB resulted in improved model fits, but this was not the case if only the previous test status was considered. This suggests that in many cases bTB infected animals are moving out of herds before being identified through testing, and that risk-based trading approaches should not rely solely on the previous test history of source herds as a proxy for future risk. Model fit was also improved by the inclusion of variables measuring bTB history of the herd, bTB in neighbouring herds, herd size, herd type, the movement network measures "in strength" and "betweenness", altitude, modelled badger abundance and county. Rainfall was not a good predictor. The most influential measures of bTB in nearby herds (a proxy for neighbourhood infection) were the proportion of herds with a history of bTB whose centroids were within 6 km, or whose boundaries were within 4 km, of the index herd. As well as informing national control and surveillance measures, our models can be used to identify areas where bTB rates are anomalously high, to prompt further investigation in these areas.

Disturbance Ecology Meets Bovine Tuberculosis (bTB) Epidemiology: A Before-and-After Study on the Association between Forest Clearfelling and bTB Herd Risk in Cattle Herds

Disturbance ecology refers to the study of discrete processes that disrupt the structure or dynamics of an ecosystem. Such processes can, therefore, affect wildlife species ecology, including those that are important pathogen hosts. We report on an observational before-and-after study on the association between forest clearfelling and bovine tuberculosis (bTB) herd risk in cattle herds, an episystem where badgers (Meles meles) are the primary wildlife spillover host. The study design compared herd bTB breakdown risk for a period of 1 year prior to and after exposure to clearfelling across Ireland at sites cut in 2015–2017. The percent of herds positive rose from 3.47% prior to clearfelling to 4.08% after exposure. After controlling for confounders (e.g., herd size, herd type), we found that cattle herds significantly increased their odds of experiencing a bTB breakdown by 1.2-times (95%CIs: 1.07–1.36) up to 1 year after a clearfell risk period. Disturbance ecology of wildlife reservoirs is an understudied area with regards to shared endemic pathogens. Epidemiological observational studies are the first step in building an evidence base to assess the impact of such disturbance events; however, such studies are limited in inferring the mechanism for any changes in risk observed. The current cohort study suggested an association between clearfelling and bTB risk, which we speculate could relate to wildlife disturbance affecting pathogen spillback to cattle, though the study design precludes causal inference. Further studies are required. However, ultimately, integration of epidemiology with wildlife ecology will be important for understanding the underlying mechanisms involved, and to derive suitable effective management proposals, if required.

Preserving identity in capture–mark–recapture studies: increasing the accuracy of minimum number alive (MNA) estimates by incorporating inter-census trapping efficiency variation

Quantifying abundance is often key to understanding ecological and evolutionary processes in wild populations. Despite shortcomings in producing accurate abundance estimates, minimum number alive (MNA) remains a widely used tool, due to its intuitive computation, reliable performance as an abundance indicator, and linkage to individual life-histories. Here, we propose a novel “efficiency-modified” MNA (eMNA) metric, which aims to preserve MNA’s favourable aspects while remedying its flaws, by incorporating (a) growth correlates to back-age individuals first captured as adults, and (b) estimates of undetected persistence beyond last capture based on time-varying capture efficiency. We evaluate eMNA through samplings of a simulated baseline population parameterised using data from a long-term demographic study of European badgers (Meles meles), under three different levels of capture efficiency (low; intermediate/“real” based on badger field data; high). We differentiate between eMNA’s performance as an abundance estimator—how well it approximates true abundance (accuracy)—and as an abundance indicator—how tightly it correlates with population abundance and changes thereof (precision). eMNA abundance estimates were negatively biased at all capture efficiencies. However, this bias was negligible at intermediate-to-high capture efficiency, particularly once low-information terminal sampling years (the first year and final three years of simulated studies) were removed. Excluding these years, eMNA under-estimated abundance by only 3.5 badgers (1.5% of population) at intermediate (real) capture efficiencies, and performed as a precise abundance indicator, with half the standard deviation of Cormack–Jolly–Seber probabilistic estimates and proving robust to inter-sampling variation in capture efficiency. Using undetected persistence probabilities to parameterise survival regression, we recreated baseline age-based survival relationships, albeit with some negative bias for under-represented ages. We offer considerations on the continued limitations of using eMNA for abundance estimates, minimum study duration for reliability, the metric’s benefits when individual identity is required, and potential for further improvement.

Is there an association between road building and bovine tuberculosis herd risk? A three time-point study in Ireland, 2011-2019

Bovine tuberculosis (bTB) is a chronic infectious disease caused by Mycobacterium bovis which results in a significant economic cost to cattle industries and governments where it is endemic. In Ireland, the European badger is the main wildlife reservoir of infection. In this study, we investigated whether (motorway) road construction was associated with an increased risk of bTB in associated cattle herds. For this study, we considered three observation periods: pre-construction (2011-2014), construction (2015-2017) and post-construction (2018-2019). We selected 1543 herds situated, based on proximity, between >50 m and <5 km of the roadworks, and extracted information about their herd-size, herd-type, inward animal movements, bTB history, and distance to the roadworks. Generalized linear mixed models were performed, whose outcome were whether a herd experienced a bTB breakdown with ≥1 or ≥3 standard reactor/s, respectively. Herds located at a distance of >3 km from the roadworks were found to be at reduced risk of a bTB breakdown over the construction period compared with those situated within 1 km of the roadworks for ≥1 reactor/s (>3 km and construction vs. <1 km: OR: 0.595, 95 % Confidence Interval (CI): 0.354-0.999) or ≥3 reactors (>3 km and construction vs. <1 km: OR: 0.431, 95 % CI: 0.174-1.067). Other previously reported risk factors such as inward movements, herd-size and herd-type were also associated with bTB risk in the final models (≥1 reactor/s and ≥3 reactors). These findings appear to be consistent with bTB breakdowns being a consequence as opposed to coincident to road construction, given the temporal and spatial consistency of the evidence. The potential for badger social group disturbance leading to the spatial spread of infection to cattle herds, as previously described in the United Kingdom, could be a hypothetical mechanism to explain these findings. However, our findings are not consistent with previous Irish studies, including recent work from another road construction project, albeit running alongside and cross over an existing road rather than construction of a new road as in this case, or experiences from national targeted badger removal. Further research is warranted to verify this pattern occurs elsewhere, and the underlying biological mechanism. Until further data are available, we recommend that badgers are vaccinated, as a precautionary measure, in advance of the commencement of major roadworks.

Spatio-temporal models of bovine tuberculosis in the Irish cattle population, 2012-2019

INTRODUCTION

Bovine tuberculosis (bTB) is an important zoonotic disease which has serious and sometimes fatal effects on both human and non-human animals. In many countries it is endemic in the cattle population and has a considerable economic impact through losses in productivity and impacts on trade. The incidence rate in Ireland varies by herd and location and it is hoped that statistical disease-mapping models accounting for both spatio-temporal correlation and covariates might contribute towards explaining this variation.

METHODS

Ireland was divided into equally sized hexagons for computational efficiency (n = 997). Different spatio-temporal random-effects models (e.g. negative binomial Besag-York-Mollié) were explored, using comprehensive data from the national bTB eradication programme to examine the association between covariates and the number of bTB cattle. Leveraging a Bayesian framework, model parameter estimates were obtained using the integrated nested Laplace approximation (INLA) approach. Exceedance probabilities were calculated to identify spatial clusters of cases.

RESULTS

Models accounting for spatial correlation significantly improved model fit in comparison to non-spatial versions where independence between regions was assumed. In our final model at hexagon level, the number of cattle (IR = 1.142, CrI: 1.108 - 1.177 per 1000), the capture of badgers (IR = 5.951, CrI: 4.482 - 7.912), percentage of forest cover (IR = 1.031, CrI: 1.020 - 1.042) and number of farm fragments (IR = 1.012, CrI: 1.009 - 1.015 per 10 fragments) were all associated with an increased incidence of bTB. Habitat suitability for badgers, percentage of dairy herds and the number of cattle movements into the herd were not. As an epidemiological tool and to suggest future work, an interactive online dashboard was developed to monitor disease progression and disseminate results to the general public.

CONCLUSION

Accounting for spatial correlation is an important consideration in disease mapping applications and is often ignored in statistical models examining bTB risk factors. Over time, the same regions in Ireland generally show highest incidences of bTB and allocation of more resources to these areas may be needed to combat the disease. This study highlights national bTB incidence rates. Shifting from national level analysis to smaller geographical regions may help identify localised high-risk areas.

The challenge of estimating wildlife populations at scale: the case of the European badger (Meles meles) in Ireland

Estimating population size in space and time is essential for applied ecology and wildlife management purposes; however, making accurate and precise estimates at large scales is highly challenging. An example is the European badger (Meles meles), a widespread and abundant mammal in Ireland. Due to their role in the epidemiology of bovine tuberculosis, the species has been culled in agriculturally dominant landscapes with the intention of reducing spillback infection to local cattle populations. Despite several studies using different approaches having estimated badger populations at different time points and scales, there remains considerable uncertainty regarding the current population and its future trajectory. To explore this uncertainty, we use published data and expert opinion to estimate a snapshot of probable badger population size using a Monte Carlo approach, incorporating variation in three key components: social group numbers, group size, and culling efficacy. Using this approach, we estimate what the badger population in Ireland would be with/without culling, assuming a steady-state population at carrying capacity, and discuss the limitations of our current understanding. The mean estimate for the badger population size was 63,188 (5–95th percentile, 48,037–79,315). Population estimates were sensitive to the assumption of mean group size across landscape type. Assuming a cessation of culling (in favour of vaccination, for example) in agricultural areas, the mean estimated population size was 92,096 (5–95th percentile, 67,188–118,881). Despite significant research being conducted on badgers, estimates on population size at a national level in Ireland are only approximate, which is reflected in the large uncertainty in the estimates from this study and inconsistencies between recording of data parameters in previous studies. Focusing on carefully estimating group size, factors impacting its variation, in addition to understanding the dynamics of repopulation post-culling, could be a fruitful component to concentrate on to improve the precision of future estimates.

Protective immunity against tuberculosis in a free-living badger population vaccinated orally with Mycobacterium bovis Bacille Calmette-Guérin

Vaccination of badgers with Mycobacterium bovis Bacille Calmette-Guérin (BCG) has been shown to protect badgers against tuberculosis in experimental trials. During the 3-year County Kilkenny BCG vaccine field study, badgers were treated orally with placebo (100% in Zone A), BCG (100% in Zone C) or randomly assigned 50%: 50% treatment with BCG or placebo (Zone B). At the end of the study, 275 badgers were removed from the trial area and subjected to detailed post-mortem examination followed by histology and culture for M. bovis. Among these badgers, 83 (30.2%) were captured for the first time across the three zones, representing a non-treated proportion of the population. Analysis of the data based on the infection status of treated animals showed a prevalence of 52% (95% CI: 40%-63%) infection in Zone A (placebo), 39% (95% CI: 17%-64%) in Zone B (placebo) and 44% (95% CI: 20%-70%) in Zone B (BCG vaccinated) and 24% (95% CI: 14%-36%) in Zone C (BCG vaccinated). There were no statistically significant differences in the proportion of animals with infection involving the lung and thoracic lymph nodes, extra-thoracic infection or in the distribution and severity scores of histological lesions. Among the 83 non-treated badgers removed at the end of the study, the infection prevalence of animals in Zone A (prevalence = 46%, 95% CI: 32%-61%) and Zone B (prevalence = 44%, 95% CI: 23%-67%) was similar to the treated animals in these zones. However, in Zone C, no evidence of infection was found in any of the untreated badgers (prevalence = 0%, 95% CI: 0%-14%). This is consistent with an indirect protective effect in the non-vaccinated badgers leading to a high level of population immunity. The results suggest that BCG vaccination of badgers could be a highly effective means of reducing the incidence of tuberculosis in badger populations.

Future Risk of Bovine Tuberculosis (Mycobacterium bovis) Breakdown in Cattle Herds 2013-2018: A Dominance Analysis Approach

Bovine tuberculosis (bTB) remains a significant endemic pathogen of cattle herds, despite multi-decadal control programmes being in place in several countries. Understanding the risks of future bTB breakdown (BD) and the associated characteristics of herds and index breakdowns could help inform risk categorisation. Such risk categories could then contribute to tailored management and policies. Here, we estimated the future risk of herd BD for the cohort of herds that were derestricted during 2013 in Ireland using multivariable logit regression models, with a dominance analysis approach. One third of herds that were derestricted in 2013 experienced a breakdown during the follow-up five year period (1469/4459; 33%). BD length was a significant predictor of future risk, primarily driven by long BDs > 230 days relative to short BDs < 130 days (OR 95%CI: 1.157–1.851), as was having had a previous BD (OR 95%CI: 1.012–1.366). Herd-size was the dominant predictor of future risk (accounted for 46% of predicted variance), suggesting significant increase in risk of future breakdown with increasing (log) herd-size (OR 95%CI: 1.378–1.609). There was significant spatial variation in future risk across counties, and it was the second most dominant predictor of future risk (25% of predicted variance). The size of index breakdowns was not a strong predictor of future risk over a 5-year period. These findings can inform a risk-based policy development.

Serological test performance for bovine tuberculosis in cattle from herds with evidence of on-going infection in Northern Ireland

The ability to accurately identify infected hosts is the cornerstone of effective disease control and eradication programs. In the case of bovine tuberculosis, accurately identifying infected individual animals has been challenging as all available tests exhibit limited discriminatory ability. Here we assess the utility of two serological tests (IDEXX Mycobacterium bovis Ab test and Enfer multiplex antibody assay) and assess their performance relative to skin test (Single Intradermal Comparative Cervical Tuberculin; SICCT), gamma-interferon (IFNγ) and post-mortem results in a Northern Ireland setting. Furthermore, we describe a case-study where one test was used in conjunction with statutory testing. Serological tests using samples taken prior to SICCT disclosed low proportions of animals as test positive (mean 3% positive), despite the cohort having high proportions with positive SICCT test under standard interpretation (121/921; 13%) or IFNγ (365/922; 40%) results. Furthermore, for animals with a post-mortem record (n = 286), there was a high proportion with TB visible lesions (27%) or with laboratory confirmed infection (25%). As a result, apparent sensitivities within this cohort was very low (≤15%), however the tests succeeded in achieving very high specificities (96–100%). During the case-study, 7/670 (1.04%) samples from SICCT negative animals from a large chronically infected herd were serology positive, with a further 17 animals being borderline positive (17/670; 2.54%). Nine of the borderline animals were voluntarily removed, none of which were found to be infected post-mortem (no lesions/bacteriology negative). One serology test negative animal was subsequently found to have lesions at slaughter with M. bovis confirmed in the laboratory.

Mycobacterium bovis Population Structure in Cattle and Local Badgers: Co-Localisation and Variation by Farm Type

Bovine tuberculosis surveillance in Northern Ireland includes Multiple-Locus Variable number tandem repeat Analysis (MLVA) to determine the Mycobacterium bovis genetic type present in both cattle and the predominant wildlife host, the European badger (Meles meles). These data are useful for investigating clusters of infection and understanding the scale at which interspecific transmission may occur. We utilised a comprehensive dataset of routinely sampled isolates from infected cattle and from badgers killed in road-traffic accidents to investigate the spatial co-location of MLVA types in, and between, the badger and cattle populations. Furthermore, we investigated the hypothesis that the type of farming enterprise might explain some variation in this relationship. MLVA types were spatially co-localised in cattle and road-traffic accident (RTA) badger hosts, indicative of a shared epidemic. Dairy herds were more likely to have at least one MLVA type in common with nearby RTA badgers, compared to non-dairy herd types. Marginally more MLVA spatial clustering was observed in non-dairy herds, which may be a consequence of relatively more between-herd movements. For the cattle population, local transmission mechanisms such as infection from contiguous herds, infectious wildlife and short-range between-herd cattle movements appear primarily to drive the epidemic: there appears to be a more limited role for long-range movements. Animal management practices are likely to be the driving force behind this observation, as beef rearing is associated with elevated numbers of animal movements compared to dairy herds.

Is moving from targeted culling to BCG-vaccination of badgers (Meles meles) associated with an unacceptable increased incidence of cattle herd tuberculosis in the Republic of Ireland? A practical non-inferiority wildlife intervention study in the Republic of Ireland (2011-2017)

Bovine tuberculosis (BTB) remains as a costly disease of cattle-herds in the Republic of Ireland (ROI). This persistence is partially attributable to the presence of M. bovis infection in a wildlife reservoir, the European badger (Meles meles). Thus, both area-wide and limited-area targeted-badger-culling have been part of the ROI-BTB control/eradication program to help reduce the future incidence of a cattle-herd BTB breakdown (i.e. a "new herd-level occurrence of BTB"). However, neither badger-culling practice can be sustained as a major component in the ongoing BTB eradication program in the ROI. Vaccination of badgers with Bacille Calmette-Guérin (BCG) has been proposed as an alternative to badger culling. Thus, in 2011, a five-year non-inferiority study was implemented in seven counties in the ROI. This study was designed to compare and contrast the cattle-herd-BTB-incidence in areas where intramuscular badger vaccination would be implemented versus the cattle-herd-BTB-incidence in the remaining area of the same county where targeted-badger-culling was maintained as the standard treatment response to probable badger-sourced BTB breakdowns. Our outcome of interest was a new cattle-herd-BTB-episode (breakdown) with a total of >2 standard skin-test (SICTT) reactors detected during the episode. Treatments (badger vaccination or targeted badger culling) were cluster allocated based on where the majority of the herd owner's land was located. To assess the impact of the two treatments, we compared the incidence-risk, of our defined outcome, for cattle herds in the area under vaccination to the outcome incidence-risk for cattle herds in the remainder of the same county after 4 and 5 years of having implemented badger vaccination. A random-effects logit model with adjustment for clustering by treatment, and statistical control of herd-type, herd-size and five-year prior-BTB-episode history was used for our analyses. Although not included in the logistic model, a relative badger density metric based on the annual number of badgers captured-per-sett-night of capturing effort was developed for each treatment area; this metric indicated that relative badger density was approximately 40 % higher in vaccination areas than in the targeted badger-culling areas during our study. Overall, our study results indicated that vaccination was not inferior to targeted badger-culling in four counties and badger vaccination was deemed to produce ambivalent results in one (County Cork North) of the seven study sites in the ROI. A post-study investigation, in County Galway, where vaccination was deemed inferior to target culling, revealed that widespread purchases of cattle from a nearby cattle mart, by herd owners in the vaccination-area, was associated with the increased herd and vaccination-area risk of BTB. No single "biasing hypothesis" was evident for the apparent vaccine inferiority in the second study site (County Monaghan) where vaccination was deemed inferior to targeted culling; hence no further investigations were conducted.

BoLA-DRB3 exon2 polymorphisms among tuberculous cattle: Nucleotide and functional variability and their association with bovine tuberculosis pathology

Bovine tuberculosis (bTB) is caused by Mycobacterium bovis and disseminated worldwide. In Argentina, the highest prevalence occurs in dairy areas. BoLA DRB3.2 is related to the adaptive immunity in mycobacterial infections. Genetic polymorphisms of this marker have been associated with resistance or susceptibility to bovine diseases. We evaluated the association between BoLA DRB3.2 polymorphisms and bTB pathology scores in dairy and beef cattle breeds of Argentina. Most bovines exhibited visible lesions compatible with tuberculosis and, furthermore, 150 (85.7%) were also positive by bacteriology. A pathology index showed a variable degree of disease, from 3 to 76 (median pathology score = 9 (IQR: 7-15)). Thirty-five BoLA DRB3.2 alleles were identified with an associated frequency from 16% to 0.3%, distributed 73% (n = 128) in heterozygosis and 27% (n = 47) in homozygosis, with 12 BoLA DRB3.2 alleles (*0101, *1101, *1501, *0201, *2707 *1001, *1002, *1201, *14011, *0501 *0902 and *0701) representing the 74.7% of the population variability. A functional analysis grouped them in 4 out of 5 clusters (A-D), suggesting a functional overlapping. Among the 90 identified genotypes, *1101/*1101, *1101/*1501 and *0101/*0101 were the most frequent (10%, 8.9% and 8.9%, respectively). No association was detected between the pathology scores and a specific DRB3.2 allele (p > .05). Animals infected with M. bovis spoligotype SB0153 showed a significantly higher pathology score than those affected by the spoligotype SB0145 (p = .018). Furthermore, the Aberdeen Angus breed exhibited highest pathological scores (p < .0001), which were associated with disseminated lesion, thus suggesting that the host component could be important to the disease progression.

Bovine tuberculosis breakdown duration in cattle herds: an investigation of herd, host, pathogen and wildlife risk factors

Background

Despite rigorous controls placed on herds which disclose ante-mortem test positive cattle to bovine tuberculosis, caused by the infection of Mycobacterium bovis, many herds in Northern Ireland (NI) experience prolonged breakdowns. These herds represent a considerable administrative and financial burden to the State and farming community.

Methods

A retrospective observational study was conducted to better understand the factors associated with breakdown duration, which was modelled using both negative binomial and ordinal regression approaches.

Results

Six explanatory variables were important predictors of breakdown length in both models; herd size, the number of reactors testing positive in the initial SICCT test, the presence of a lesioned animal at routine slaughter (LRS), the count of M. bovis genotypes during the breakdown (MLVA richness), the local herd-level bTB prevalence, and the presence of herds linked via management factors (associated herds). We report that between 2008 and 2014, mean breakdown duration in NI was 226 days (approx. seven months; median: 188 days). In the same period, however, more than 6% of herds in the region remained under movement restriction for more than 420 days (13 months); almost twice as long as the mean. The MLVA richness variable was a particularly important predictor of breakdown duration. We contend that this variable primarily represents a proxy for beef fattening herds, which can operate by purchasing cattle and selling animals straight to slaughter, despite prolonged trading restrictions. For other herd types, the model supports the hypothesis that prolonged breakdowns are a function of both residual infection within the herd, and infection from the environment (e.g. infected wildlife, contiguous herds and/or a contaminated environment). The impact of badger density on breakdown duration was assessed by including data on main sett (burrow) density. Whilst a positive association was observed in the univariate analysis, confounding with other variables means that the contribution of badgers to prolonged breakdowns was not clear from our study. We do not fully reject the hypothesis that badgers are implicated in prolonging bTB breakdowns via spillback infection, but given our results, we posit that increased disease risk from badgers is unlikely to simply be a function of increasing badger density measured using sett metrics.

Interspecific visitation of cattle and badgers to fomites: A transmission risk for bovine tuberculosis?

In Great Britain and Ireland, badgers (Meles meles) are a wildlife reservoir of Mycobacterium bovis and implicated in bovine tuberculosis transmission to domestic cattle. The route of disease transmission is unknown with direct, so-called "nose-to-nose," contact between hosts being extremely rare. Camera traps were deployed for 64,464 hr on 34 farms to quantify cattle and badger visitation rates in space and time at six farm locations. Badger presence never coincided with cattle presence at the same time, with badger and cattle detection at the same location but at different times being negatively correlated. Badgers were never recorded within farmyards during the present study. Badgers utilized cattle water troughs in fields, but detections were infrequent (equivalent to one badger observed drinking every 87 days). Cattle presence at badger-associated locations, for example, setts and latrines, were three times more frequent than badger presence at cattle-associated locations, for example, water troughs. Preventing cattle access to badger setts and latrines and restricting badger access to cattle water troughs may potentially reduce interspecific bTB transmission through reduced indirect contact.

Variation in Mycobacterium bovis genetic richness suggests that inwards cattle movements are a more important source of infection in beef herds than in dairy herds

Background

We used genetic Multi-Locus VNTR Analysis (MLVA) data gathered from surveillance efforts to better understand the ongoing bovine tuberculosis (bTB) epidemic in Northern Irish cattle herds. We modelled the factors associated with Mycobacterium bovis MLVA genotype richness at three analytical scales; breakdown level, herd level, and patch level, and compared the results between dairy and non-dairy production types.

Results

In 83% of breakdowns and in 63% of herds, a single MLVA genotype was isolated. Five or more MLVA genotypes were found in less than 3 % of herds. Herd size and the total number of reactors were important explanatory variables, suggesting that increasing MLVA genotype richness was positively related to increases in the number of host animals. Despite their smaller relative size, however, the highest MLVA genotype richness values were observed in non-dairy herds. Increasing inwards cattle movements were important positive predictors of MLVA genotype richness, but mainly in non-dairy settings.

Conclusions

The principal finding is that low MLVA genotype richness indicates that small-scale epidemics, e.g. wildlife, contiguous farms, and within-herd recrudescence, are important routes of M. bovis infection in cattle herds. We hypothesise that these mechanisms will maintain, but may not explicitly increase, MLVA genotype richness. The presence of elevated MLVA richness is relatively rare and likely indicates beef fattening enterprises, which purchase cattle from over long distances. Cattle movements were furthermore an important predictor of MLVA genotype richness in non-dairy herds, but not in dairy herds; this may represent reduced cattle purchasing levels in dairy enterprises, compared to beef. These observations allude to the relative contribution of different routes of bTB infection between production types; we posit that infection associated with local factors may be more evident in dairy herds than beef herds, however in beef herds, inwards movements offer additional opportunities for introducing M. bovis into the herd.

Electronic supplementary material

The online version of this article (10.1186/s12866-019-1530-7) contains supplementary material, which is available to authorized users.

Characteristics of Northern Irish cattle herds without bovine tuberculosis infection

BACKGROUND

Despite ongoing eradication efforts, bovine tuberculosis (bTB) is endemic in cattle herds in Northern Ireland (NI). This disease has serious implications for the economy, farming and animal welfare. Previous research identified a population of herds which have remained free from bTB infection for 10 years (2004-2014). Understanding the characteristics of these herds may have important implications for eradication efforts, such as spatially targeted interventions.

METHODS

A cluster analysis and a retrospective case-control analysis was conducted to compare bTB- free herds with herds which experienced prolonged infection (ie, bTB breakdowns lasting more than ≥ 365 days).

RESULTS

Only small, localised clusters of herds which have remained free from bTB were revealed, thus limiting the potential for spatially targeted interventions. The results illustrated the importance of herd size to disease status; over 27 per cent of the bTB-free herds had up to 10 animals. However, the data also showed that there were no inward movements in the year before the bTB skin test in those herds which remained free from bTB.

CONCLUSIONS

Attention should therefore be given to the cattle movement network in NI to better understand the risk associated with cattle purchasing.

Liver fluke (Fasciola hepatica) co-infection with bovine tuberculosis in cattle: A prospective herd-level assessment of herd bTB risk in dairy enterprises

Co-infection of tuberculosis (TB) and helminths is recognized as a significant problem in regions where such pathogens are endemic and chronic cases exist. Co-infection can modulate the immune system leading to interference with diagnostic tests, increased pathological impacts and pathogen persistence. However, research has found that such interactions between pathogens can be context and species specific. Recent studies have suggested that liver fluke, Fasciola hepatica, infection may impact on immunological responses and diagnostics for bovine tuberculosis (bTB; caused by Mycobacterium bovis) in cattle. Where evidence of such interaction exists, there would be an onus on policy makers to adjust eradication programs to minimize impacts. We assessed the association between herd-level bTB breakdown risk and seasonal variation in liver fluke exposure based on 5,753 bulk tank milk (BTM) samples from 1,494 dairy herds across Northern Ireland. BTM was tested by an IDEXX antibody specific enzyme-linked immunosorbent assay (ELISA) using the 'f2' antigen as a detection agent. The ELISA determined the result based on a sample to (known) positive ratio (S/P%) from which binary status and categories of exposure were derived. Associations were tested using multivariable random effects models. Models predicting bTB risk were not improved with the inclusion of liver fluke exposure levels. Variations in modelling liver fluke exposure (S/P%, binary, categories of exposure) and bTB risk (skin test breakdowns, post-mortem confirmed breakdowns, breakdown size and lag effects) also failed to support associations (neither positive nor negative) between the pathogens at herd-level. These results, along with previously published animal-level data from Northern Ireland, suggest that the nexus between bTB and F. hepatica may have small size effects at the population-level. However, our results also highlight the high prevalence of F. hepatica in cattle in our study population, and therefore we cannot fully discount the potential hypothesis of population-level depression of immune response to M. bovis due to co-infection.

Characteristics of Mycobacterium bovis infected herds tested with the interferon-gamma assay

The IFN-γ (interferon gamma) assay is used in Ireland as an ancillary diagnostic test to the single intradermal comparative tuberculin test (SICTT) to maximise the detection of Mycobacterium bovis infected animals (bTB) in cattle herds. Understanding the relationships between herd and animal risk factors and IFN-γ test results is critical to enable the development and evaluation of policy measures on how best to use the test. In this study, we set out to characterise Irish herds with IFN-γ test positive animals in terms of herd size, number of SICTT reactors and number of IFN-γ positive tests, and to evaluate the IFN-γ test in terms of the test cut-off values. The results showed that larger herds with more SICTT reactors were likely to have more IFN-γ positives in the herd, and herds with an IFN-γ test positive animal that was also positive for bTB lesions at post-mortem had higher numbers of IFN-γ positive animals in the herd. Raising the cut-off values for the IFN-γ test only marginally decreased the combined sensitivity of the IFN-γ and the SICTT for diagnosis of bTB lesioned animals. The analysis has provided valuable information on the performance of the IFN-γ test as it is used under current bTB infection levels in Ireland.

Liver fluke (Fasciola hepatica) co-infection with bovine tuberculosis (bTB) in cattle: A retrospective animal-level assessment of bTB risk in dairy and beef cattle

Bovine tuberculosis (bTB), caused by Mycobacterium bovis, remains a persistent problem for cattle industries in endemic countries. The frequency, quality, and performance of tests, and the presence of wildlife reservoirs, have been identified as impediments to eradication. Recently, exposure to helminth infection (Fasciola hepatica) has been associated negatively with the disclosure of bTB. Here, for the first time, we assess impact of concurrent infections of Fasciola hepatica and the disclosure of bTB at the animal-level using large surveillance datasets. We utilized a dataset of 138,566 animal records from an abattoir from Northern Ireland (2011-2013). The presence of F. hepatica infection was assessed from macroscopic tissue inspection at abattoir. Multivariable models were developed to assess co-infection associations with bTB status based on: Single Intradermal Comparative Tuberculin Test (SICTT), lesion, bacteriological confirmation, including either all animals, or only skin-test negative animals (lesions at routine slaughter; LRS; confirmed nonreactors at routine slaughter; cNRs) or positive (reactors) animals alone, respectively. The relationship between skin tuberculin reaction sizes and fluke status was also explored for a subset of animals with field recordings (n = 24,680). Controlling for known risk factors (e.g., climatic, herd, and individual level characteristics), we did not find significant associations between the SICTT (standard or severe interpretation), lesion, nor confirmation status of animals and their liver fluke status. The only exception was a negative association between liver fluke positivity, and LRS or cNRs, respectively; though effect-sizes were small (e.g., LRS Odds-Ratio: 0.87; 95% CI: 0.76-1.00). There was limited evidence of a relationship between tuberculin reaction sizes during SICTT testing and liver fluke infection status. These data do not support the contention that the detection of bTB using skin-tests or reactor postmortem follow-up may be compromised by co-infection at a population level, but the relationship with lesion formation (pathogenesis) may indicate an impact for postmortem surveillance.

Modeling as a Decision Support Tool for Bovine TB Control Programs in Wildlife

Computer modeling has a long history of association with epidemiology, and has improved our understanding of the theory of disease dynamics and provided insights into wildlife disease management. A summary of badger bovine TB models and their role in decision making is presented, from a simple initial SEI model, to SEIR (inclusion of a recovered category) and SEI₁I₂ (inclusion of two stages of disease progression) variants, and subsequent spatially-explicit individual-based models used to assess historical badger management strategies. The integration of cattle into TB models allowed comparison of the predicted impacts of different badger management strategies on cattle herd breakdown rates, and provided an economic dimension to the outputs. Estimates of R₀ for bovine TB in cattle and badgers are little higher than unity implying that the disease should be relatively easy to control, which is at odds with practical experience. A cohort of recent models have suggested that combined strategies, involving management of both host species and including vaccination may be most effective. Future models of badger vaccination will need to accommodate the partial protection from infection and likely duration of immunity conferred by the currently available vaccine (BCG). Descriptions of how models could better represent the ecological and epidemiological complexities of the badger-cattle TB system are presented, along with a wider discussion of the utility of modeling for bovine TB management interventions. This includes consideration of the information required to maximize the utility of the next generation of models.

Bovine Tuberculosis in Britain and Ireland - A Perfect Storm? the Confluence of Potential Ecological and Epidemiological Impediments to Controlling a Chronic Infectious Disease

Successful eradication schemes for bovine tuberculosis (bTB) have been implemented in a number of European and other countries over the last 50 years. However, the islands of Britain and Ireland remain a significant aberration to this trend, with the recent exception of Scotland. Why have eradication schemes failed within these countries, while apparently similar programs have been successful elsewhere? While significant socio-economic and political factors have been discussed elsewhere as key determinants of disease eradication, here we review some of the potential ecological and epidemiological constraints that are present in these islands relative to other parts of Europe. We argue that the convergence of these potential factors may interact additively to diminish the potential of the present control programs to achieve eradication. Issues identified include heterogeneity of diagnostic testing approaches, the presence of an abundant wildlife reservoir of infection and the challenge of sustainably managing this risk effectively; the nature, size, density and network structure of cattle farming; potential effects of Mycobacterium bovis strain heterogeneity on disease transmission dynamics; possible impacts of concurrent endemic infections on the disclosure of truly infected animals; climatological differences and change coupled with environmental contamination. We further argue that control and eradication of this complex disease may benefit from an ecosystem level approach to management. We hope that this perspective can stimulate a new conversation about the many factors potentially impacting bTB eradication schemes in Britain and Ireland and possibly stimulate new research in the areas identified.

Assessment of the safety of Bacillus Calmette-Guérin vaccine administered orally to badgers (Meles meles)

European badgers (Meles meles) are a wildlife reservoir for Mycobacterium bovis (M. bovis) in parts of England, Wales and Ireland, constituting a potential source of tuberculosis (TB) infection for cattle. Vaccination of badgers against TB is one of the tools available for helping reduce the prevalence of bovine TB in badgers, made possible by the licensing in 2010 of Bacillus Calmette-Guérin (BCG) vaccine for intramuscular administration to badgers (BadgerBCG). However, practical limitations associated with administering an injected vaccine to wild animals make an oral, bait-delivered form of the vaccine highly desirable. Evaluation of the safety of oral BCG to badgers and the environment is a mandatory step on the road to licensing an oral vaccine. This study had the following objectives: (a) to determine whether adverse effects followed the oral administration of BCG vaccine to badgers; (b) to measure the quantity and frequency of BCG excreted in the faeces of vaccinated badgers; and (c) to assess whether there was evidence of the vaccine spreading to unvaccinated, 'sentinel' badgers sharing the same environment as vaccinated animals. We report here that the oral administration per badger of ≥6.4 × 10⁹ cfu BCG, followed 14 days later by a single oral dose of ≥6.4 × 10⁷ cfu BCG caused no adverse physical effects and did not affect the social behaviour and feeding habits of the vaccinated animals. BCG was cultured from the faeces of two of nine vaccinated animals (372 cfu/g and 996 cfu/g, respectively) approximately 48 h after the higher dose of BCG was administered and by one of the nine vaccinated animal (80 cfu/g) approximately 24 h after receiving the lower dose of BCG. We found no evidence for the transmission of BCG to unvaccinated, sentinel, badgers housed with the vaccinated animals despite the occasional excretion of BCG in faeces.

Mycobacterium bovis hosted by free-living-amoebae permits their long-term persistence survival outside of host mammalian cells and remain capable of transmitting disease to mice

Bovine tuberculosis (TB) is a zoonotic disease caused by Mycobacterium bovis. Despite intensive TB control campaigns, there are sporadic outbreaks of bovine TB in regions declared TB free. It is unclear how M. bovis is able to survive in the environment for long periods of time. We hypothesized that Free-living amoebae (FLA), as ubiquitous inhabitants of soil and water, may act as long-term reservoirs of M. bovis in the environment. In our model, M. bovis would be taken up by amoebal trophozoites, which are the actively feeding, replicating and mobile form of FLA. Upon exposure to hostile environmental conditions, infected FLA will encyst and provide an intracellular niche allowing their M. bovis cargo to persist for extended periods of time. Here, we show that five FLA species (Acanthamoeba polyphaga, Acanthamoeba castellanii, Acanthamoeba lenticulata, Vermamoeba vermiformis and Dictyostellium discoideum) are permissive to M. bovis infection and that the M. bovis bacilli may survive within the cysts of four of these species for over 60 days. We further show that exposure of M. bovis-infected trophozoites and cysts to Balb/c mice leads to pulmonary TB. This work describes for the first time that FLA carrying M. bovis can transmit TB.

Assessment of concurrent infection with bovine viral diarrhoea virus (BVDV) and Mycobacterium bovis: A herd-level risk factor analysis from Northern Ireland

Bovine viral diarrhoea virus (BVDV) is a significant pathogen of cattle, leading to severe economic and animal-welfare impacts. Furthermore, the pathogen has been associated with impacting the progression or spread of other pathogens (e.g. Mycobacterium bovis, the causative agent of bovine tuberculosis (bTB)). During this study we investigated (i) risk factors for BVDV at a herd-level and (ii) whether there was any association between BVDV and herd-level bTB risk. The data for this study were gathered from a voluntary BVDV control programme in Northern Ireland (2013-2015) based on the identification of virus positive animals through tissue tag testing of calves. We assigned a herd-level BVDV status to 2827 participating herds, where a herd was assumed "infected" if one or more animals tested positive for BVDV. Two model suites were developed. Firstly, we assessed risk factors for BVDV herd status using multivariable logit random-effects modelling, aggregating to the calendar year level (2013-2015; n=4828; model 1). Secondly, we aggregated data across the three years of the study to give an overall status for the whole study period (n=2827; logistic model 2). Risk factors included year, herd-type, herd size, number of births, inward trade moves, calf mortality, and region. Furthermore, the herd-level bovine tuberculosis status (based on the single intradermal comparative cervical tuberculin (SICCT) test outcomes, or confirmation at post-mortem), or the size of bTB breakdowns (number of SICCT test positive animals), of herds was also investigated to assess whether there was an association (co-infection) with herd BVDV status. The final models suggested that BVDV herd status was positively associated with increased levels of calf mortality, herd size, number of births, the number of BVDV tests undertaken and the number of animals introduced to the herd. There was a significant univariable positive association between BVDV status, and SICCT breakdown risk, breakdown size and confirmed bTB status in model 2. However, there was no evidence of significant associations between bTB status (using SICTT status, confirmed status or herd breakdown size) and BVDV status in final multivariable models when controlling for other significant confounders. These results provide information for action for the future control and eradication of BVDV in Northern Ireland, though these data provide little support for the hypothesised association between BVDV and bTB status at herd-level. Further animal-level analyses are necessary to investigate whether there is support for a BVD-bTB co-infection association, including the impact of co-infection on the severity of infection.

Control and eradication of tuberculosis in cattle: a systematic review of economic evidence

Bovine tuberculosis (TB) is a disease of zoonotic importance for which control and eradication programmes have been carried out in many countries for decades. While the impact of these programmes on public health is still uncertain, the impact on trade is significant because of movement restrictions for animals, costs of testing and culling. The objective of this systematic review was to provide a contribution to the general debate over costs against benefits for the control and eradication of bovine TB in cattle. The search strategy was performed on four electronic databases following the PRISMA (Preferred Reporting Items for Systematic Reviews and Meta-Analyses) guidelines. The selection process, data abstraction and quality appraisal were carried out according to the Cochrane Collaboration guidelines. The search identified 66 articles out of which eight fulfilled the inclusion criteria. The evidence gathered in this review by combining the conclusions of the most methodologically sound articles supports the idea that, when multiple cost and benefit components are taken into account, efforts to control or eradicate bovine TB may be effective in reducing disease prevalence, economically viable and worth doing.

The Region of Difference Four is a Robust Genetic Marker for Subtyping Mycobacterium caprae Isolates and is Linked to Spatial Distribution of Three Subtypes

Alpine Mycobacterium caprae isolates found in cattle and red deer display at least three genetic variations in the region of difference four (RD4) that can be used for further differentiation of the isolates into the subtypes 'Allgäu', 'Karwendel' and 'Lechtal'. Each genomic subtype is thereby characterized by a specific nucleotide deletion pattern in the 12.7-kb RD4 region. Even though M. caprae infections are frequently documented in cattle and red deer, little is known about the transmission routes. Hence, robust markers for M. caprae subtyping are needed to gain insight into the molecular epidemiology. For this reason, a rapid and robust multiplex PCR was developed for the simultaneous detection of three M. caprae RD4 subtypes and was used to subtype a total number of 241 M. caprae isolates from animals (145 cattle, 95 red deer and one fox) from Bavaria and Austria. All three subtypes occur spatially distributed and are found in cattle and in red deer suggesting transmission between the two species. As subtypes are genetically stable in both species it is hypothesized that the described genetic variations developed within the host due to 'within-host replication'. The results of this study recommend the genomic RD4 region as a reliable diagnostic marker for M. caprae subtype differentiation.

In situ adaptive response to climate and habitat quality variation: spatial and temporal variation in European badger (Meles meles) body weight

Variation in climatic and habitat conditions can affect populations through a variety of mechanisms, and these relationships can act at different temporal and spatial scales. Using post-mortem badger body weight records from 15 878 individuals captured across the Republic of Ireland (7224 setts across ca. 15 000 km(2) ; 2009-2012), we employed a hierarchical multilevel mixed model to evaluate the effects of climate (rainfall and temperature) and habitat quality (landscape suitability), while controlling for local abundance (unique badgers caught/sett/year). Body weight was affected strongly by temperature across a number of temporal scales (preceding month or season), with badgers being heavier if preceding temperatures (particularly during winter/spring) were warmer than the long-term seasonal mean. There was less support for rainfall across different temporal scales, although badgers did exhibit heavier weights when greater rainfall occurred one or 2 months prior to capture. Badgers were also heavier in areas with higher landscape habitat quality, modulated by the number of individuals captured per sett, consistent with density-dependent effects reducing weights. Overall, the mean badger body weight of culled individuals rose during the study period (2009-2012), more so for males than for females. With predicted increases in temperature, and rainfall, augmented by ongoing agricultural land conversion in this region, we project heavier individual badger body weights in the future. Increased body weight has been associated with higher fecundity, recruitment and survival rates in badgers, due to improved food availability and energetic budgets. We thus predict that climate change could increase the badger population across the Republic of Ireland. Nevertheless, we emphasize that, locally, populations could still be vulnerable to extreme weather variability coupled with detrimental agricultural practice, including population management.

Understanding and managing bTB risk: perspectives from Ireland

There is substantial variation in herd risk for bovine tuberculosis (bTB) in Ireland, with most herds playing little to no role in the ongoing endemic. In infected areas, bTB persistence (affecting one or a group of herds) is a key feature of the infection. In this paper, we present our current understanding and management of bTB risk in Ireland, based on a detailed review of research and policy. There is close interaction between science and policy in Ireland, seeking both to understand and effectively manage bTB risk. Detailed research on bTB persistence is presented, including current understanding of the relative importance of different infection sources, which can include residual infection in cattle and/or re-infection, either from local sources or following cattle introduction. In recent years, there have been three primary drivers for policy change, including scientific advances, ongoing improvements to programme supports, and ongoing programme review. In this review, three key future programme challenges are identified. Although good progress is being made, eradication has not yet been achieved. Firstly, a key question concerns the additional effort that will be required, to move towards final eradication. Secondly, a percentage of non-infected animals are falsely positive to current testing methods. This is an ongoing challenge, given the imperfect specificity of test methods but will become more so, as the positive predictive value falls with reducing bTB prevalence. Finally, there is a need to re-engage with the farming community, so that they play a much greater role in programme ownership.