- Research article
- Open Access
Bovine tuberculosis visible lesions in cattle culled during herd breakdowns: the effects of individual characteristics, trade movement and co-infection
BMC Veterinary Research volume 13, Article number: 400 (2017)
Bovine tuberculosis (bTB), caused by Mycobacterium bovis, remains a significant problem for livestock industries in many countries worldwide including Northern Ireland, where a test and slaughter regime has utilised the Single Intradermal Comparative Cervical Tuberculin (SICCT) test since 1959.
We investigated the variation in post-mortem confirmation based on bTB visible lesion (VL) presence during herd breakdowns using two model suites. We investigated animal-level characteristics, while controlling for herd-level factors and clustering. We were interested in potential impacts of concurrent infection, and therefore we assessed whether animals with evidence of liver fluke infection (Fasciola hepatica; post-mortem inspection), M. avium reactors (animals with negative M. bovis-avium (b-a) tuberculin reactions) or Bovine Viral Diarrhoea Virus (BVDV; RT-PCR tested) were associated with bTB confirmation.
The dataset included 6242 animals removed during the 14 month study period (2013–2015). bTB-VL presence was significantly increased in animals with greater b-a reaction size at the disclosing SICCT test (e.g. b-a = 5-9 mm vs. b-a = 0 mm, adjusted Odds ratio (aOR): 14.57; p < 0.001). M. avium reactor animals (b-a < 0) were also significantly more likely to disclose VL than non-reactor animals (b-a = 0; aOR: 2.29; p = 0.023). Animals had a greater probability of exhibiting lesions with the increasing number of herds it had resided within (movement; log-herds: aOR: 2.27–2.42; p < 0.001), if it had an inconclusive penultimate test result (aOR: 2.84–3.89; p < 0.001), and with increasing time between tests (log-time; aOR: 1.23; p = 0.003). Animals were less likely to have VL if they were a dairy breed (aOR: 0.79; p = 0.015) or in an older age-class (e.g. age-quartile 2 vs. 4; aOR: 0.65; p < 0.001). Liver fluke or BVDV variables were not retained in either multivariable model as they were non-significantly associated with bTB-VL status (p > 0.1).
Our results suggest that neither co-infection of liver fluke nor BVDV had a significant effect on the presence of VLs in this high-risk cohort. M. avium tuberculin reactors had a significantly increased risk of disclosing with a bTB lesion, which could be related to the impact of co-infection with M. avium subsp. paratuberculosis (MAP) affecting the performance of the SICCT however further research in this area is required. Movements, test history, breed and age were important factors influencing confirmation in high-risk animals.
Despite decades of research into bovine tuberculosis (bTB), caused by Mycobacterium bovis, there still remains limitations in our ability to explain the variation in bTB detectability and confirmation infection within cattle [1,2,3,4]. Furthermore, animal-level analyses relating to pathological findings from surveillance have not been as frequently investigated, in comparison with epidemiological studies concentrating on herd-level factors [5, 6].
Bovine TB in Northern Ireland remains a significant chronic problem, with continual maintenance of the disease in the national herd despite a coordinated compulsory statutory scheme in place since the 1950s [7, 8]. This national scheme has been predicated on an ante-mortem test and cull regime, and supplemented with additional continual slaughter house surveillance systems . The Single Intradermal Comparative Cervical Tuberculin (SICCT) test is the statutory ante-mortem test, and as dictated by legislation, all animals deemed to be standard reactors are compulsorily culled. Additional animals may be culled if ‘severe’ interpretation of the test is implemented during herd breakdowns, or through the ancillary use of the Interferon-γ test or serological tests, in an attempt to remove exposed animals [9,10,11]. However, all of these tests have limitations in terms of sensitivity (especially SICCT and serological tests) and specificity (especially the interferon-γ test), meaning some infected animals can be missed from breakdown herds or some non-infected animals needlessly culled [9, 12,13,14].
Recently, co-infection of other pathogens has been highlighted as a potential factor affecting the efficacy of the current systems to identify and remove infected animals from endemic bTB countries like the UK and Ireland [15,16,17,18,19,20,21,22]. Bovine Viral Diarrhoea Virus (BVDV) , liver fluke (Fasciola hepatica; ) and Johne’s disease (caused by Mycobacterium avium subsp. paratuberculosis [21, 23]) have all been suggested to be implicated in the epidemiology of bovine TB either by actively exacerbating infection, potentially leading to increase M. bovis transmission, increasing host susceptibility, or by affecting the performance of tests used to detect bTB infection.
Here we investigate the variation in post-mortem confirmation (presence of bTB visible lesions; VL) of at-risk animals culled in Northern Ireland during 2013–2015. All animals removed were ante-mortem positive for bTB, or negative at-risk animals in direct contact with infected animals. Building on the previous work by O’Hagan et al.  in Northern Ireland, we hypothesised that confirmation would be significantly associated with variations in the SICCT test reaction size, the animal’s test history, movement history (trade), and herd-level characteristics where the animal resided before slaughter [1, 3, 4]. Controlling for such risk factors, we then also utilised these data to explore the hypotheses that co-infection (exposure to BVDV, liver fluke or non-bTB Mycobacterium) may be significantly associated (either in positively or negatively) with the probability of bTB confirmation.
The study design was a retrospective observational case-control study. The study population was animals slaughtered as part of bTB herd breakdowns (where bTB has been detected in the herd via one or more positive statutory tests) in Northern Ireland. Cases were exposed animals with bTB visible lesions (VL), a metric of post-mortem confirmation, as used elsewhere (e.g. [1,2,3,4, 24,25,26,27]). Controls were animals without VL detected at slaughter. 97% of animals with lesions were laboratory confirmed (histopathology and/or culture) and a previous risk factor study of bovine reactors from Northern Ireland suggested that little difference in parameter estimates occurred when comparing VL with bacteriologically confirmed samples . Extensive details of the TB program in Northern Ireland are presented in , and specific details of the processing of reactors and the determination of cases is presented in . Briefly, reactor animals are sent to slaughter where post-mortem examination (PME) within the abattoir for evidence of infection through the presence of VLs is undertaken . During the study period, 99% of samples were processed in one specific slaughterhouse, which would reduce the inter-abattoir variability in detection of lesions [1, 2]. The routine PME includes palpation of the lungs, liver and tongue, and incision of the lungs, heart, and liver. In reactor animals, examination of the lymph nodes of the head, chest and mesenterium, lungs, pleura, peritoneum, and the prescapular, popliteal, iliac, precrural lymph nodes was undertaken .
The dataset contained records from animals that were culled from the 27th of November 2013 to the 27th of January 2015. All animals were culled as they were deemed reactors to the Single Intradermal Comparative Cervical Tuberculin (SICCT) test (either standard interpretation or severe interpretation), gamma interferon test positive or animals in-contact with test positive animals.
Metrics of co-infection
The primary goal of this study was to investigate potential interactions between endemic diseases and M. bovis confirmation (VL). We investigated whether there was any evidence of interaction between BVDV positivity or exposure, liver fluke post-mortem status, and Mycobacterium avium complex exposure. All culled animals in the dataset had a tissue sample (ear-tag notch) available, which was stored at −70°c. These samples were used to prospectively test for BVDV virus using a commercial real time RT-PCR test (Virotype BVDV RT-PCR, Qiagen) as described before [22, 28]. All samples were first batch screened by RT-PCR, and individual samples within positive batches were then tested using antigen-capture ELISA test (Herd Check BVDV antigen/serum plus, IDEXX Laboratories) to detect viraemic animals. BVDV-bTB interaction was then assessed using two variables: (i) a direct BVDV RT-PCR positive test assigned to the animal, indicating the animal has active infection with BVDV (ii) indirect herd-level exposure, with all animals within a herd with a positive BVDV animal being considered exposed. Liver fluke status was also assessed with three metrics of infection using surveillance data : (i) Fascioliasis recorded at slaughter, indicating that active infection (the presence of Fasciola hepatica in the liver) was occurring at time of death (1 = active infection; 0 = non-active infection) (ii) Liver fluke damage recorded at slaughter, indicating that the animal had been infected in the past, with resulting scar tissue in the liver, however the animal was not currently infected (1 = historic infection; 0 = non-historic infection) (iii) the animal was either found to have Fascioliasis or liver fluke damage at slaughter (1 = active infection or historic infection; 0 = non-active and non-historic infection). Please note, that category three is related to the other categories and so could not be entered into the same model.
Potential exposure to non-bTB Mycobacterium was assessed using the relative reaction size to the bovine and M. avium Purified Protein Derivative (PPD) tuberculin used during the disclosing test [30,31,32,33]. The comparative skin test (SICTT) used in Northern Ireland assigns test status to animals based on the bovine tuberculin reaction size subtracted by the avian reaction size (b-a), with positive bias values indicating a reaction to the bovine tuberculin greater than to the avian tuberculin (i.e. net bias), and negative reaction sizes indicating a reaction to avian tuberculin greater than to the bovine tuberculin. A bTB “standard reactor”, one deemed positive to the statutory test and required to be culled by law, is one whereby the b-a > 4 mm. The relationship between PPD reaction size and VL presence was non-linear (Additional file 1; and see Additional file 2) – and was modelled by categorising the values into one of six types, three non-standard reactors (1. Avium reactor (b-a < 0 mm), 2. Non-reactor (b-a = 0 mm); 3. Inconclusive bovine reactor (b-a > 0 & < =4 mm)) and three standard reactor categories (4. “Small bovine reactors” (b-a > 4 & b-a < 10 mm); 5. “Large bovine reactors” (b-a > =10 & b-a < 20 mm); 6. “Hyper bovine reactors” (b-a > =20 mm)). We considered M. avium reactors as potentially indicative of exposure to non-bTB Mycobacterium, including Mycobacterium avium subsp. paratuberculosis the agent of Johne’s disease in cattle [30, 32, 33], a notifiable pathogen in Northern Ireland. To further explore the impact of M. avium reactors on lesion presence, the effect of the individual bovine PPD reaction was assessed using descriptive statistics and univariable models. For this, a single intradermal tuberculin (SIT) test results was established as positive if the net bovine PPD reaction was >4 mm.
We concentrated on animal-level variables in this study, while controlling for herd-level clustering. Factors investigated in our models are presented Table 1. These variables were chosen as we hypothesised that they could be influential in the ante-mortem test variability, confirmation, or pathology or because previous research suggested that they may be important (e.g. [3,4,5, 16, 25]).
Response variable and modelling approach
We developed models where the outcome was binary - bTB VL presence or absence [1,2,3,4]. The primary source of data for these analyses was the Northern Ireland Animal and Public Health Information System (APHIS; ). The presence of a lesion at slaughter was modelled using a multivariable binary logit random effects model, with the random effect capturing the potential non-independence of animals coming from the same herds . The significance of this clustering was measured using a likelihood ratio test, compared with a regular logit regression model. Unlike previous studies , we did not include slaughter house as a random effect, as 99% of all animals within the dataset were processed at the same abattoir.
Throughout we used a two-step model building strategy – univariable model screening and multivariable model building. The functional form was assessed using LOWESS graphing, with suitable transforms of predictors implemented to improve model fit (e.g. natural logarithm transform; Additional file 1; Additional file 3; and see Additional file 2). We fitted separate univariable models for each independent variable of interest, and all variables that were significant at p < 0.15 were candidates for a multivariable model building process. The variables tested during the model building for the logistic regression model are presented in Table 1 below. Correlations amongst independent variables were assessed, and where necessary, variables were dropped if highly correlated (i.e. to reduce the possibility of multi-collinearity; [35, 36]). The decision to include a variable amongst a correlated pair was based on which variable improved the model fit more, as assessed by Akaike’s Information Criteria (AIC). Backwards elimination was then used to establish parsimonious models , with final competing models compared using AIC and Bayesian Information Criteria (BIC) values [37, 38]. Models with lower Information Criteria (IC) values were considered preferred models. It should be noted, that ICs could not be used to compare models of differing number of observations . In that case (for example, model with or without the b-a variable), were considered simply as separate competing models, and duly reported. First order interaction terms were fitted for variables retained after backwards model building, and only retained if they improved the fit of the model and were significant at an p < 0.05. Covariate patterns were assessed using visual plots. The discriminatory ability of the final models were assessed by estimating the Area Under the ROC Curve (AUC). Models with an AUC > 0.7 were considered “adequate” . All modelling was undertaken with Stata SE 14 (Stata Corp., 2015).
Descriptive results of the dataset
The dataset included a total of 6242 animals with an ante-mortem and post-mortem record, from 1335 herds; however not all variables in the dataset were complete. Overall, 5198 animals were recorded in the dataset as “positive” skin reactors (5198/6242; 83%). A further 458 animals (458/6242; 7.3%) were designated as gamma positive, with 586 (586/6242; 9.4%) being designated as negative in contacts (NIC). Of these 6242 animals with an associated post-mortem result, 40% (2504) had visible lesions disclosed at slaughter. Of the 5198 standard reactor animals within the dataset, 2444 (47%) disclosed with visible lesions at slaughter.
Overall, there were 5698 animals with a b-a value recorded at the disclosing SICCT test (Table 2). Descriptive statistics relating to the breakdown and VL percentage are presented in Table 2. Not all animals were tested for gamma interferon, as this test was a supplementary test to the SICCT tests. There were 355 of these animals with gamma interferon as their disclosing test (355/5698; 6.2%), 20 of which were VL positive (5.6%). There were 20 animals from 5698 that were negative to all ante-mortem tests (20/5698; 0.35%); 2 were VL positive (10%). Five hundred ninety-six animals were inconclusive reactors to both standard and severe interpretation SICCT tests (596/5698; 10%), 142 of which were VL positive (24%). There were 259 animals that were inconclusive to severe interpretation but negative to both gamma and standard (259/5698; 4.6%; Table 2), 19 of which were VL positive (7.3%).
Factors associated with lesion presence - Univariable models
Multivariable random effect logit model – Lesion presence
After univariable model fitting, there were 14 independent variables offered to the full model (p < 0.15 in univariable models). Correlation analysis amongst independent variables found that the lifetime number of market movements and the number of herds (log(no. associated herds)) the animal resided in during its lifetime were strongly correlated (r = 0.83; p < 0.001). We retained only log(no. associated herds) variable to avoid collinearity in our model, and because this was a far superior predictor of lesion presence in the univariable model (ΔAIC: 208). Purchase also significantly correlated with log(no. associated herds) (r = 0.46; p < 0.001), but was not significantly associated with the outcome at univariable level (p = 0.114) and so was also dropped from the model building process.
Two competing models were run before model building with either the b-a reaction size included or the test result (skin test positive, gamma positive, or NIC animals). This was due to the number of animals without a b-a value (n = 544 missing values); 81% (440/544) of these animals were NIC animals.
The final model (Table 4) including b-a as a predictor significantly explained variation in the dataset (Wald χ2 (df: 9) = 772.86; p < 0.001), and displayed significant variation in animal risk across herds (σ: 1.13; ρ: 0.28; LR test of ρ = 0: p < 0.001). The model exhibited adequate discriminatory ability with an AUC of 0.75 (95%CI: 0.74–0.77). The model was based on 5698 observations, with an average of 4.3 animals per herd. The probability of an animal disclosing a lesion at slaughter varied significantly by b-a reaction size, the number of herds the animal resided in during its lifetime, breed and penultimate test result. Animals with net positive bovine increase (b-a > 0) had a significantly higher risk of having a lesion at slaughter (P < 0.001), than non-reactors. Furthermore, this pattern increased with increasing b-a size; for example, in comparison with non-reactors (b-a = 0), animals with reaction sizes of 1-4 mm were 6.31 times more likely to have a lesion (p < 0.001), whereas animals with reaction sizes 5-9 mm were 14.57 times more likely to have a lesion (p < 0.001). Animals who were M. avium tuberculin reactors (b-a < 0; n = 151) were also significantly more likely to disclose with a lesion relative to non-reactor animals (OR: 2.29; p = 0.023).
There was significant positive relationship between the log-number of herds an animal resided within and lesion risk (log-herds: OR: 2.27; p < 0.001). Dairy breed animals had significantly lower odds of having a lesion at slaughter relative to non-dairy animals (OR: 0.79; P = 0.015). Animals that were positive reactors on their penultimate skin test under severe interpretation were 3.98 times more likely to have a lesion at slaughter than negative animals. There was no significant difference between animals without a penultimate test result and those who tested negative (p = 0.201).
The second model significantly explained variation in the dataset (Wald χ2 (df: 9) = 448.18; p < 0.001; Table 5), and revealed significant variation in the probability of lesion disclosure across herds (significant herd random effect; LR test of ρ = 0: p < 0.001). The model exhibited fair to adequate discriminatory ability with an AUC of 0.69 (95%CI: 0.67–0.70). The final model contained 5823 animals, from 1286 herds, with an average of 4.5 animals per herd. Animals were significantly more likely to disclose with lesions if they were removed due to being SICCT test reactors (OR: 14.62; p < 0.001) in comparison with negative in-contact animals. There was no significant difference in the probability of VL presence between negative in-contacts and animals removed due to gamma interferon results (p = 0.738). Animals with a positive penultimate severe interpretation result were significantly more likely to have a lesion than (severe interpretation) test negative animals (OR = 2.84; p < 0.001). There was no difference between animals without a penultimate test result and test negative animals (p = 0.135). There was significant variation in the probability of a lesion across age quartiles; the highest risk was associated with quartiles 1 and 2. Animals in age quartile 3 (OR: 0.63) and quartile 4 (OR: 0.65) were of significantly lower odds of having a lesion relative to quartile 2 (P < 0.001), respectively. There was no significant difference between age quartile 1 and 2 (p = 0.446). There was a significant higher odds of lesion detection with the increasing (log) number of herds an animal resided within during their life (OR per unit increase: 2.42; p < 0.001). Finally, there was a significant positive relationship between lesion presence and time since last (penultimate) test (OR per log unit increase in days: 1.23; p = 0.003).
Additional analysis of non-reactors and net M. avium PPD reactors
There were 151 animals with net increases in M. avium PPD reactions (b-a < 0), 37 of which would have been SIT reactors (24.5%). 11.92% (18/151) of these animals had VLs. 7/37 SIT positive animals had VLs (18.9%); 11/103 non-SIT animals had VLs (9.7%). This difference in the probability of lesion disclosure was non-significant (OR: 2.18; p = 0.136). The mean M. avium PPD reaction size for avium reactors with lesions was 5.82 mm (n = 11; SD: 3.12; Range: 2-12 mm), while the mean M. avium PPD reaction size for avium reactors without lesions was 4.01 mm (n = 103; SD: 2.05; Range: 1–12 mm), a significant difference (linear regression: p = 0.010).
For non-reactors (b-a = 0), 22/468 (4.7%) of animals had VLs. There was a significantly greater proportion of VL positive animals with positive SIT test results (16.0% (4/25)) than SIT negative animals (4.1%; 18/443; OR: 4.50; p = 0.011). For animals with a b-a = 0 test reading, the mean M. avium PPD reaction size was non-significantly larger in animals with a VL (0.83 mm; range 0-4 mm) than animals without a VL (0.49 mm; range 0-4 mm) at slaughter (linear regression: p = 0.242).
Overall, this study has shown that there is significant variation in the probability of cattle exhibiting a VL at slaughter that were culled as part of herd breakdowns in Northern Ireland. This variation can be attributed to animal-level characteristics, exposure during life histories, and potentially to exposure to other Mycobacteria.
Co-infections have been highlighted as a significant modulator of immune responses to the challenge of bacterial infections, including tuberculosis [40,41,42,43,44]. Recent studies suggest that bovine TB could be affected, either through mechanisms of increased pathology, increased transmission, susceptibility or decreased detection, by pathogens such as BVDV [15, 22, 45], liver fluke (F. hepatica; 18, 19, 44], and Johne’s Disease (JD; [20, 42, 46]; but see ).
During the present study BVDV exposure was assessed in two ways – a positive BVDV RT-PCR test result, or an indirect exposure metric based on the presence of a test positive animal disclosed within the animal’s final herd. The first metric represent Persistently Infected (PI) animals, which were rare in the dataset (0.27%), with the second metric representing potentially ‘exposed’ cattle. We did not find any significant (at P < 0.05) animal-level associations between these BVDV statuses and the probability of an animal having a lesion at slaughter. This could be partially due to the small number of test positive animals found within the cohort, reducing the power to detect an effect. Other research from this cohort suggested that the apparent indirect exposure of animals in herds where BVDV positive cattle were found was negatively associated with the minimum lesion counts found in animals – a finding which would not be generally consistent with the proposed mechanisms by which BVDV might impact bTB epidemiology [15, 22, 45, 48, 49]. Some have suggested that acute BVDV infection could lead to the rapid progression of bTB in co-infected animals . However, experimental studies failed to show significant differences in M. bovis shedding patterns between co-infected (BVDV and M. bovis) and infected (M. bovis only) animals . Previous research suggested that acute BVDV infection could temporarily compromise ante-mortem tests for bTB, by suppressing the immunological response to tuberculin PPDs . Herd-level research from Northern Ireland suggested that, when controlling for confounders (especially herd-size which is a herd-level risk factor for both bTB and BVDV), there was no significant relationship between bTB and BVDV herd risk based on 2827 herds . Results from populations outside Northern Ireland suggest that acute BVDV cases may in some cases exacerbate infection  or impact on ante-mortem tests ability to detect M. bovis in truly infected animals , but that these affects are not readily demonstrated with herd- and animal-level epidemiological population-level datasets in Northern Ireland . However, as pointed out previously , biosecurity practices implemented in the management of BVDV could potentially have benefits for bTB management (but see ).
Liver fluke (F. hepatica) infection has been suggested to negatively impact the control of M. bovis in Britain and Ireland [18, 19, 44, 52]. Experimental work has found that a proportion of co-infected animals with Bacille Calmette Guerin (BCG) can be misdiagnosed using IFN-γ or SICCT tests . Claridge et al.  presented both experimental and epidemiological evidence to suggest that co-infection with liver fluke was associated with reduced SICCT tests disclosure (underascertainment of 27–38%) in England and Wales. Experimental data presented in that paper suggested that co-infected animals exhibited significantly lower b-a reactions to the skin test at 10 and 21 weeks post-innoculation . However, at both time points all co-infected animals would still have been deemed standard reactors (i.e. b-a > 4 mm). Despite these experimental and epidemiological data, during the present study, three metrics of liver fluke infection (Fascioliasis, fluke damage, and either Fascioliasis/fluke damage) were tested and none were significantly associated with the probability of an animal having a visible lesion at slaughter. However, other research from this population has suggested a negative relationship between maximum visible lesion size and liver fluke presence has been found . This appears perhaps contrary to the hypothesis to an immunological suppression hypothesis [54,55,56], but reflects other recent data suggesting the potential for the induction of latency . Further research is on-going into this latter hypothesis (Byrne et al. In prep.).
There was evidence from the present study that reaction to M. avium Purified Protein Derivative (PPD) tuberculin (negative b-a values), was significantly associated with visible lesion presence at slaughter. These animals would not have been classed as standard reactors to bovine tuberculosis, suggesting that animals with bTB visible lesions that strongly react to M. avium tuberculin will have increased risk of going undetected. The M. avium PPD comparator was introduced to increase specificity in areas where non-specific reactions could occur , however at a potential loss of sensitivity relative to the SIT test. For non-reactor animals and M. avium PPD reactors, a proportion of animals strongly reacted to the bovine PPD, and would have been classed as SIT reactors. It is possible that such M. avium reactors have been exposed to non-bovis mycobacteria, including Mycobacterium avium subsp. paratuberculosis (MAP) the etiological agent of Johne’s disease (JD; [30, 32, 33, 58]). Previous research has suggested that the comparative skin test can be used as a preliminary diagnostic for JD in cattle . Gilardoni et al.  suggests animals with JD react positively to both bovine and avian tuberculins, but generally greater intensity to avian PPD. However, because of the non-specific nature of the avium tuberculin reaction, it is not possible to disaggregate exposure to MAP or other members of the M. avium complex. Recently in Ireland, Kennedy et al.  found a significant association between the reactions sizes (most markedly in animals reacting to both avian and bovine PPD) during a SICCT test and an ELISA test for MAP. Increases in tuberculin sizes at day 10 and day 16 post-inoculation was associated with increases in sample/positive ratio for the JD ELISA blood test. However, this study failed to elucidate whether tuberculin reaction sizes can be indicative of exposure to MAP, or if their results were due to cross reactivity. Our findings need to be further investigated with a larger dataset, and from cohorts that include non-breakdown herds, to verify whether this is a general pattern. Further research investigating specifically JD and bTB is on-going in Northern Ireland (Lahuerta-Marin, Byrne et al., personal communication). It should be noted that previous research has found significant effects of JD on bTB in cattle herds ([32, 42, 46, 59, 60], and vice-versa: ), and so co-management should be considered if managing both infections when present in the same herd, and careful attention to the interpretation of comparative tests where co-infection is suspected.
Individual characteristics, herd-level factors, and life histories
Animals classed as of a dairy breed exhibited significant differences to beef breed types in post-mortem confirmation of bTB VL. Dairy breeds were dominated by Friesian and Holstein breeds (91% of total); in comparison, beef breeds were made up of a number of dominant breed types (Limousin, Charolais, Aberdeen Angus, Simmental, and Belgian Blue (86% of total)). We found in this high risk cohort that dairy breed animals tended to be significantly less likely to disclose with a VL, suggesting diagnostics for bTB may be behaving differently amongst herd types . Recent research suggests that some tests may have differing characteristics depending on breed type , with data from chronic herds in Northern Ireland suggesting that the SICCT test may not perform as well in dairy herds as it does in beef herds . Genetic susceptibility may also be a factor , with evidence from an Irish cattle population to suggest that Holstein cattle are the least susceptible breed examined, with Simmental and Charolais beef breeds being most susceptible . In the current study dataset, for breeds with >300 animals represented, the highest proportion of animals confirmed were Charolais (46.08% confirmed bacteriologically) and Simmental (45.05% confirmed), while the least likely breeds to confirm were Friesian (32.22% confirmed) and Holstein (37.97% confirmed) at post-mortem (Additional file 4; also see ). It is currently difficult to ascertain whether the effect is due to increased resistance, non-specific skin test reactions or a reduced pathogenesis and lesion formation in dairy cattle . However, there is evidence to suggest that confirmed breakdown recurrence is more frequent in dairy relative to non-dairy herds in Britain and Northern Ireland [65, 66]. Dairy herds tend to have an older age profile, which also affects risk (see below), however in this cohort there was little difference in the mean age of dairy breed animals (4.0 years) and non-dairy animals (3.8 years).
Age-class was another significant variable affecting the probability of an animal being confirmed (with VL; similar to ). The results of the present study suggested an increased risk for animals within the first and second quartiles, relative to other age-classes. Brooks-Pollock et al.  working in GB found that bTB risk increased with age up to the highest risk age group at 12 and 36 months old before risk declining. Potentially this pattern emerged due to the detection and removal of infected animals, with surviving cohorts of animals could have lower rates of infection. The same authors presented data on VL confirmation which suggested that confirmation risk declined with age. Anergy due to repeated testing could be a factor affecting the immunological response to the tuberculin tests with age [68, 69]. Interestingly, Clegg et al.  when undertaking a case-control study of non-reactor bTB confirmed animals in Ireland found that cases tended to be older than the national average age at slaughter, mirroring previous results  that suggested confirmation risk increased with age in non-reactors.
We found that the movement history of animals was a significant factor associated with increased probability of an animal being confirmed to have lesions found at abattoir. This suggests that movement through multiple premises may increase an animal’s risk significantly, or increase the pathological response if exposed during their lifetime. Lifetime moves were not previously associated with animal-level risk in Northern Ireland  in a study assessing the post-mortem confirmation risk of reactors. However, movement has been highlighted as a bTB risk factor in other studies (e.g. [71, 72]). In the present study, there was a significant positive relationship between increasing age and movement (p < 0.001), and with beef animals relative to dairy breed animals (p < 0.001). Therefore, the compounding effects of potential increased exposure (across multiple premises and animals), older animals, and being more likely a beef animal may have contributed to this finding. Stress due to movement may also be a significant driver of this pattern, as stress is known to negatively affect the immunocompetency of cattle [73, 74]. Lifetime number of herds an animal resided within significantly correlated with lifetime market moves and the number of animal purchases (i.e. buying in) the disclosing herds made in the 3 years before the disclosing test. Therefore, animals which resided in many herds, also made many market moves, and markets are known to be highly connected nodes within cattle trade networks in Northern Ireland  and elsewhere (e.g. GB, Denmark and Italy; [76,77,78]), which potentially increases exposure to infectious hosts and contributes to disease spread [70, 71]. Herds that purchased many animals also tended to buy-in animals with greater market-move histories (e.g. animals in herds in the fourth quartile of purchasing had on average 2.2 more market moves than animals in the first purchasing quartile; linear regression; p < 0.001). Finally, there was also an association between movement history and the time between the penultimate and disclosing tests – for example, for each additional market move, the time between tests increased by 5.5 days (linear regression; p < 0.001). While animals should be tested on average once per year in Northern Ireland (as part of an Annual Herd Test (AHT) or other follow-up tests), this may not always be the case with animals that move frequently across herds which undertake testing at different times of the year. While the average time between penultimate and disclosing test in this high-risk cohort was approximately 6 months (mean: 156 days), the longest time period between tests was almost 2 years (665 days). This pattern could be a factor in the time to detection of infection.
Finally, we provide evidence that the penultimate SICCT test result was significantly associated with post-mortem presence of a lesion at slaughter. Severe interpretation reactors were significantly more likely to confirm with a lesion at death. Similar results with a significantly increased risk of post-mortem confirmation in reactors has been demonstrated in other studies in Northern Ireland  and Ireland . This could indicate that such animals in high-risk herds represent exposed animals in the early stages of infection that may have been missed by the SICCT test and so could pose a future risk [79, 80]. Indeed, Clegg et al. [79, 80] found that inconclusive reactors (animals which had reaction sizes of >2 mm on the bovine tuberculin injection site and between 1 and 4 mm > the avian response) had significantly increased future risk of bTB relative to non-reactor animals within the same herds or in comparison with animals in the national herds in Ireland. Furthermore, Lahuerta-Marin et al.  also showed how SICCT negative, but INF-γ test positive animals left on farm constituted to farms a significant future breakdown risk. In herds with significant recurrent problems, indicators of future risk could be considered when managing herd breakdown risk.
Limitations and future direction
There were some limitations with regards this study. Firstly, post-mortem diagnostics for bovine TB has limitations in terms of variable sensitivity (e.g. see [12, 13]), though the slaughter house surveillance does exhibit high specificity (>99%). Similarly, abattoir carcass assessment for liver fluke infection can exhibit limited sensitivity and specificity at the animal level (e.g. 68%SE, 88% SP in Scotland; ), though there is no assessments of Northern Irish abattoir performance yet nor how transferrable the results are from other countries . Part of the issue relates to inter-abattoir variation in surveillance effectiveness [1, 2], however, we largely controlled for this by having 99% of samples being processed through a single experienced slaughter house. While the sensitivity and specificity of the RT-PCR BVDV test is very high [>99%; 28], the animal level prevalence was very low (lower than expected given routine surveillance data), which may have diminished the power to detect an effect, if present. We did however, endeavour to overcome this, by assessing indirect effects of exposure, as BVDV is highly contagious within herds with actively excreting animals . Future studies will directly assess the impact of MAP on bTB test diagnostics in Northern Ireland, to help refine and assess the validity of the present findings. Liver fluke exposure will be assessed systematically in Northern Ireland using a prospective study design assessing bulk milk antigens using high performance ELISA tests at herd levels, to assess the herd level nexus between liver fluke and bTB risk.
There is significant animal-level variation in the post-mortem presence of VL from animals removed as part of bTB breakdowns in Northern Ireland. In this dataset, we found limited robust evidence of BVDV impacts on metrics of post-mortem confirmation, concurring with previous herd-level analyses. We did not find evidence to support liver fluke infection impacting on VL presence during this study. Exposure to Mycobacteria other than M. bovis may also have an impact on the SICCT test performance, with M. avium PPD reactors exhibiting higher probability of confirmation than non-reactor animals. We hypothesize that this effect could be partially due to exposure to Mycobacterium bovis subsp. paratuberculosis (MAP) the etiological agent of Johne’s disease, but this needs further confirmatory research. There was a strong significant relationship between visible lesion presence and increasing positive bovine tuberculin bias over the avian tuberculin reaction. Visible lesion presence in this high risk cohort was related to breed (decreased in dairy breeds relative to non-dairy), age (older animals are less likely to have lesions), movement history (increased risk with how many herds an animal resided in) and whether the animal was inconclusive to the SICCT at its penultimate test (previous test results can be a predictor of future risk).
Area under the ROC curve
Bovine Viral Diarrhoea Virus
Mycobacterium avium subsp. paratuberculosis
bTB visible lesion
Frankena K, White PW, O'Keeffe J, Costello E, Martin SW, van Grevenhof EM, More SJ. Quantification of the relative efficiency of factory surveillance in the disclosure of tuberculosis lesions in attested Irish cattle. Vet Rec. 2007;161:679–84.
Olea-Popelka F, Freeman Z, White P, Costello E, O'Keeffe J, Frankena K, Martin W, More S. Relative effectiveness of Irish factories in the surveillance of slaughtered cattle for visible lesions of tuberculosis, 2005-2007. Irish Vet J. 2012;65:2.
Shittu A, Clifton-Hadley RS, Ely ER, Upton PU, Downs SH. Factors associated with bovine tuberculosis confirmation rates in suspect lesions found in cattle at routine slaughter in great Britain, 2003–2008. Prev Vet Med. 2013;110(3):395–404.
O’Hagan MJ, Courcier EA, Drewe JA, Gordon AW, McNair J, Abernethy DA. Risk factors for visible lesions or positive laboratory tests in bovine tuberculosis reactor cattle in Northern Ireland. Prev Vet Med. 2015 Jul 1;120(3):283–90.
Humblet MF, Boschiroli ML, Saegerman C. Classification of worldwide bovine tuberculosis risk factors in cattle: a stratified approach. Vet Res. 2009;40(5):1–24.
Skuce RA, Allen AR, McDowell SW. Herd-level risk factors for bovine tuberculosis: a literature review. Vet Med Int. 2012 Jun;28:2012.
Abernethy DA, Denny GO, Menzies FD, McGuckian P, Honhold N, Roberts AR. The Northern Ireland programme for the control and eradication of Mycobacterium bovis. Vet Microbiol. 2006;112(2):231–7.
Robinson PA. Farmers and bovine tuberculosis: Contextualising statutory disease control within everyday farming lives. J. Rural. Stud. 2017;55:168–80.
Lahuerta-Marin A, McNair J, Skuce R, McBride S, Allen M, Strain SA, Menzies FD, McDowell SJ, Byrne AW. Risk factors for failure to detect bovine tuberculosis in cattle from infected herds across Northern Ireland (2004–2010). Res Vet Sci. 2016 Aug 31;107:233–9.
McCallan L, Brooks C, Couzens C, Young F, McNair J, Byrne AW. Assessment of serological tests for diagnosis of bovine tuberculosis. Vet Rec. 2017 Apr 5:vetrec-2016. doi:10.1136/vr.104272
McCallan L, Brooks C, Couzens C, Young F, McNair J, Byrne AW. Performance of serological tests for bovine tuberculosis in cattle herds in Northern Ireland. BioRxiv. 235184. https://doi.org/10.1101/235184.
Lahuerta-Marin, A., Milne, M.G., McNair, J., Skuce, R., McBride, S., Menzies, F., McDowell, S.J.W., Byrne, A.W., Handel, I.G., Bronsvoort, M.B.C. 2017. Bayesian latent class estimation of sensitivity and specificity parameters of diagnostic tests for bovine tuberculosis in chronic herds-Northern Ireland. Vet J, In review.
De la Rua-Domenech R, Goodchild AT, Vordermeier HM, Hewinson RG, Christiansen KH, Clifton-Hadley RS. Ante mortem diagnosis of tuberculosis in cattle: a review of the tuberculin tests, γ-interferon assay and other ancillary diagnostic techniques. Res Vet Sci. 2006;81(2):190–210.
Nuñez-Garcia J, Downs SH, Parry JE, Abernethy DA, Broughan JM, Cameron AR, Cook AJ, de la Rua-Domenech R, Goodchild AV, Gunn J, More SJ. Meta-analyses of the sensitivity and specificity of ante-mortem and post-mortem diagnostic tests for bovine tuberculosis in the UK and Ireland. Prev Vet Med. 2017. doi:10.1016/j.prevetmed.2017.02.017
Kao RR, Gravenor MB, Charleston B, Hope JC, Martin M, Howard CJ. Mycobacterium bovis shedding patterns from experimentally infected calves and the effect of concurrent infection with bovine viral diarrhoea virus. J R Soc Interface. 2007;4(14):545–51.
Broughan JM, Judge J, Ely E, Delahay RJ, Wilson G, Clifton-Hadley RS, Goodchild AV, Bishop H, Parry JE, Downs SH. A review of risk factors for bovine tuberculosis infection in cattle in the UK and Ireland. Epidemiol. Infect. 2016;144:2899–926.
Godfray HC, Donnelly CA, Kao RR, Macdonald DW, McDonald RA, Petrokofsky G, Wood JL, Woodroffe R, Young DB, McLean AR. A restatement of the natural science evidence base relevant to the control of bovine tuberculosis in great Britain. Proc R Soc B Biol Sci. 2013 Oct 7;280(1768):20131634.
Flynn RJ, Mulcahy G, Welsh M, Cassidy JP, Corbett D, Milligan C, Andersen P, Strain S, McNair J. Co-infection of cattle with Fasciola hepatica and Mycobacterium bovis–immunological consequences. Transbound Emerg Dis. 2009;56(6–7):269–74.
Claridge J, Diggle P, McCann CM, Mulcahy G, Flynn R, McNair J, Strain S, Welsh M, Baylis M, Williams DJ. Fasciola hepatica is associated with failure to detect bovine tuberculosis in dairy cattle. Nat Commun. 2012;3:853.
Álvarez J, De Juan L, Bezos J, Romero B, Sáez JL, Gordejo FR, Briones V, Moreno MÁ, Mateos A, Domínguez L, Aranaz A. Interference of paratuberculosis with the diagnosis of tuberculosis in a goat flock with a natural mixed infection. Vet Microbiol. 2008;128(1):72–80.
Álvarez J, De Juan L, Bezos J, Romero B, Sáez JL, Marqués S, Domínguez C, Mínguez O, Fernández-Mardomingo B, Mateos A, Domínguez L. Effect of paratuberculosis on the diagnosis of bovine tuberculosis in a cattle herd with a mixed infection using interferon-gamma detection assay. Vet Microbiol. 2009;135:389–93.
Byrne AW, Guelbenzu-Gonzalo M, Strain SA, McBride S, Graham J, Lahuerta-Marin A, Harwood R, Graham DA, McDowell S. Assessment of concurrent infection with bovine viral diarrhoea virus (BVDV) and Mycobacterium bovis: a herd-level risk factor analysis from Northern Ireland. Prev Vet Med. 2017;141:38–47.
Aagaard C, Govaerts M, Meikle V, Gutiérrez-Pabello JA, McNair J, Andersen P, Suárez-Güemes F, Pollock J, Espitia C, Cataldi A. Detection of bovine tuberculosis in herds with different disease prevalence and influence of paratuberculosis infection on PPDB and ESAT-6/CFP10 specificity. Prev Vet Med. 2010;96(3):161–9.
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(3):451–64.
Downs SH, Broughan JM, Goodchild AV, Upton PA, Durr PA. Responses to diagnostic tests for bovine tuberculosis in dairy and non-dairy cattle naturally exposed to Mycobacterium Bovis in great Britain. Vet J. 2016 Oct 31;216:8–17.
Nalapa DP, Muwonge A, Kankya C, Olea-Popelka F. Prevalence of tuberculous lesion in cattle slaughtered in Mubende district, Uganda. BMC Vet Res. 2017 Mar 21;13(1):73.
Gallagher MJ, Higgins IM, Clegg TA, Williams DH, More SJ. Comparison of bovine tuberculosis recurrence in Irish herds between 1998 and 2008. Prev. Vet. Med. 2013 Sep 1;111(3):237–44.
Guelbenzu-Gonzalo MP. Benchmarking and control of bovine viral diarrhoea in Northern Ireland dairy and suckler herds. Thesis, Queens University Belfast. 2015.
Byrne AW, McBride S, Lahuerta-Marin A, Guelbenzu M, McNair J, Skuce RA, McDowell SW. Liver fluke (Fasciola Hepatica) infection in cattle in Northern Ireland: a large-scale epidemiological investigation utilising surveillance data. Parasit Vectors. 2016;9(1):209.
Ameni G, Aseffa A, Engers H, Young D, Gordon S, Hewinson G, Vordermeier M. High prevalence and increased severity of pathology of bovine tuberculosis in Holsteins compared to zebu breeds under field cattle husbandry in central Ethiopia. Clin Vaccine Immunol. 2007;14(10):1356–61.
Bermingham ML, More SJ, Good M, Cromie AR, Higgins IM, Brotherstone S, Berry DP. Genetics of tuberculosis in Irish Holstein-Friesian dairy herds. J Dairy Sci. 2009;92(7):3447–56.
Gilardoni LR, Paolicchi FA, Mundo SL. Bovine paratuberculosis: a review of the advantages and disadvantages of different diagnostic tests. Rev Argent Microbiol. 2012;44(3):201–15.
Kennedy AE, Byrne N, O'Mahony J, Sayers RG. Investigations and implications of associations between mycobacterial purified protein derivative hypersensitivity and MAP-antibody ELISA in Irish dairy cows. Res. Vet. Sci. 2017;115:13–6.
Houston R. A computerised database system for bovine traceability. Rev Sci Tech. 2001;20(2):652.
Dohoo I, Martin W, Stryhn H. Veterinary epidemiologic research. Canada: University of Prince Edward Island; 2009.
Belsley DA, Kuh E, Welsch RE. Regression diagnostics: identifying influential data and sources of collinearity (Vol. 571). New York: Wiley; 2005.
Burnham KP, Anderson DR. Model selection and multimodel inference a practical information-theoretic approach. New York: Springer; 2002.
Schwarz G. Estimating the dimension of a model. Ann Stat. 1978;6(2):461–4.
Hosmer, D.W., Lemeshow, S. 2000. Applied logistic regression. 2. New York, NY: John Wiley & Sons, USA.
Kamal SM, El Sayed Khalifa K. Immune modulation by helminthic infections: worms and viral infections. Parasite Immunol. 2006 Oct 1;28(10):483–96.
Resende Co T, Hirsch CS, Toossi Z, Dietze R, Ribeiro-Rodrigues R. Intestinal helminth co-infection has a negative impact on both anti-mycobacterium tuberculosis immunity and clinical response to tuberculosis therapy. Clin. Exp. Immunol. 2007;147(1):45–52.
Bezos J, de Juan L, Romero B, Álvarez J, Mazzucchelli F, Mateos A, Domínguez L, Aranaz A. Experimental infection with Mycobacterium Caprae in goats and evaluation of immunological status in tuberculosis and paratuberculosis co-infected animals. Vet Immunol Immunopathol. 2010 Feb 15;133(2):269–75.
Salgame P, Yap GS, Gause WC. Effect of helminth-induced immunity on infections with microbial pathogens. Nat Immunol. 2013;14(11):1118–26.
Lucena AN, Cuartero LG, Mulcahy G, Zintl A. The immunoregulatory effects of co-infection with Fasciola hepatica: From bovine tuberculosis to Johne's disease. Vet J. 2017;222;9–16.
Charleston B, Hope JC, Carr BV, Howard CJ. Masking of two in vitro immunological assays for Mycobacterium Bovis (BCG) in calves acutely infected with non-cytopathic bovine viral diarrhoea virus. Vet Rec. 2001;149(16):481–4.
Aranaz A, Bezos J, Álvarez J, Romero B, Lozano F, Paramio JL, López-Sánchez J, Mateos A, Domínguez L. Assessment of diagnostic tools for eradication of bovine tuberculosis in cattle co-infected with Mycobacterium Bovis and M. Avium subsp. paratuberculosis. Vet Res. 2006;37(4):593–606.
Dunn JR, Kaneene JB, Grooms DL, Bolin SR, Bolin CA, Bruning-Fann CS. Effects of positive results for Mycobacterium Avium subsp paratuberculosis as determined by microbial culture of feces or antibody ELISA on results of caudal fold tuberculin test and interferon-γ assay for tuberculosis in cattle. J Am Vet Med Assoc. 2005;226:429–35.
Menzies FD, Neill SD. Cattle-to-cattle transmission of bovine tuberculosis. Vet J. 2000;160:92–106.
Pollock JM, Neill SD. Mycobacterium bovis infection and tuberculosis in cattle. Vet J. 2002;163:115–27.
Monies RJ, Head JC. Bovine tuberculosis in housed calves. Vet Rec. 1998;145(25):743.
Cardwell JM, Van Winden S, Beauvais W, Mastin A, De Glanville WA, Hardstaff J, Booth RE, Fishwick J, Pfeiffer DU. Assessing the impact of tailored biosecurity advice on farmer behaviour and pathogen presence in beef herds in England and Wales. Prev Vet Med. 2016;135:9–16.
Flynn RJ, Mannion C, Golden O, Hacariz O, Mulcahy G. Experimental Fasciola hepatica infection alters responses to tests used for diagnosis of bovine tuberculosis. Infect Immun. 2007;75:1373–81.
Byrne AW, Graham J, Brown C, Donaghy A, Guelbenzu-Gonzalo M, McNair J, Skuce R, Allen A, McDowell S. Modelling the variation in skin-test tuberculin reactions, post-mortem lesion frequency and pathology in cattle culled during bovine tuberculosis herd breakdowns. Transbound Emerg Dis. 2018. In press.
Ezenwa VO, Etienne RS, Luikart G, Beja-Pereira A, Jolles AE. Hidden consequences of living in a wormy world: nematode-induced immune suppression facilitates tuberculosis invasion in African buffalo. Am Nat. 2010;176:613–24.
Jolles AE, Ezenwa VO, Etienne RS, Turner WC, Olff H. Interactions between macroparasites and microparasites drive infection patterns in free-ranging African buffalo. Ecology. 2008;89:2239–50.
O'Leary S, O'sullivan MP, Keane J. IL-10 blocks phagosome maturation in mycobacterium tuberculosis–infected human macrophages. Am J Respir Cell Mol Biol. 2011;45(1):172–80.
Garza-Cuartero L, O'sullivan J, Blanco A, McNair J, Welsh M, Flynn RJ, Williams D, Diggle P, Cassidy J, Mulcahy G. Fasciola Hepatica infection reduces Mycobacterium Bovis burden and mycobacterial uptake and suppresses the pro-inflammatory response. Parasite Immunol. 2016 Jul 1;38(7):387–402.
Kalis CH, Collins MT, Hesselink JW, Barkema HW. Specificity of two tests for the early diagnosis of bovine paratuberculosis based on cell-mediated immunity: the Johnin skin test and the gamma interferon assay. Vet. Microbiol. 2003;97(1):73–86.
Walravens K, Marché S, Rosseels V, Wellemans V, Boelaert F, Huygen K, Godfroid J. IFN-γ diagnostic tests in the context of bovine mycobacterial infections in Belgium. Vet Immunol Immunopathol. 2002 Sep 10;87(3):401–6.
Seva J, Sanes JM, Ramis G, Mas A, Quereda JJ, Villarreal-Ramos B, Villar D, Pallares FJ. Evaluation of the single cervical skin test and interferon gamma responses to detect Mycobacterium Bovis infected cattle in a herd co-infected with Mycobacterium Avium subsp. paratuberculosis. Vet Microbiol. 2014;171(1):139–46.
Lilenbaum W, Ferreira R, Marassi CD, Ristow P, Oelemann WMR, Fonseca LDS. Interference of tuberculosis on the performance of ELISAs used in the diagnosis of paratuberculosis in cattle. Braz J Microbiol. 2007;38(3):472–7.
Allen AR, Minozzi G, Glass EJ, Skuce RA, McDowell SWJ, Woolliams JA, Bishop, SC. Bovine tuberculosis: the genetic basis of host susceptibility. Proc R Soc B. 2010. p.rspb20100830.
Richardson IW, Bradley DG, Higgins IM, More SJ, McClure J, Berry DP. Variance components for susceptibility to Mycobacterium bovis infection in dairy and beef cattle. Genet. Sel. Evol. 2014;46(1):77.
Wright DM, Allen AR, Mallon TR, McDowell SW, Bishop SC, Glass EJ, Bermingham ML, Woolliams JA, Skuce RA. Field-isolated genotypes of Mycobacterium bovis vary in virulence and influence case pathology but do not affect outbreak size. PloS One. 2013;8(9):e74503.
Doyle LP, Gordon AW, Abernethy DA, Stevens K. Bovine tuberculosis in Northern Ireland: risk factors associated with time from post-outbreak test to subsequent herd breakdown. Prev Vet Med. 2014;116(1):47–55.
Karolemeas, K., McKinley, T.J., Clifton-Hadley, R.S., Goodchild, A.V., Mitchell, A., Johnston,W.T., Conlan, A.J., Donnelly, C.A.,Wood, J.L., 2011. Recurrence of bovine tuberculosis breakdowns in great Britain: risk factors and prediction. Prev Vet Med. 102, 22–29.
Brooks-Pollock E, Conlan AJ, Mitchell AP, Blackwell R, McKinley TJ, Wood JL. Age-dependent patterns of bovine tuberculosis in cattle. Vet Res. 2013;44(1):97.
Radunz BL, Lepper AWD. Suppression of skin reactivity to bovine tuberculin in repeat tests. Australian Vet J. 1985;62(6):191–4.
Coad M, Clifford D, Rhodes SG, Hewinson RG, Vordermeier HM, Whelan AO. Repeat tuberculin skin testing leads to desensitisation in naturally infected tuberculous cattle which is associated with elevated interleukin-10 and decreased interleukin-1 beta responses. Vet Res. 2010;41(2):1–12.
Clegg TA, Good M, More SJ. Risk factors for cattle presenting with a confirmed bTB lesion at slaughter, from herds with no evidence of within-herd transmission. Prev Vet Med. 2016;126:111–20.
Gilbert M, Mitchell A, Bourn D, Mawdsley J. Cattle movements and bovine tuberculosis in great Britain. Nature. 2005;435:491.
Ramirez-Villaescusa AM, Medley GF, Mason S, Green LE. Herd and individual animal risks associated with bovine tuberculosis skin test positivity in cattle in herds in south west England. Prev Vet Med. 2009;92:188–98.
Carroll JA, Forsberg NE. Influence of stress and nutrition on cattle immunity. Vet Clin N Am Food Anim Pract. 2007;23:105–49.
Ryan TJ, Livingstone PG, Ramsey DS, De Lisle GW, Nugent G, Collins DM, Buddle BM. Advances in understanding disease epidemiology and implications for control and eradication of tuberculosis in livestock: the experience from New Zealand. Vet Microbiol. 2006;112:211–9.
Brown E, Marshall AH, Mitchell H, Byrne AW. Analysing cattle movements in Northern Ireland using social network analysis. In: Brennan M, Lindberg A, editors. Society for Veterinary Epidemiology and Preventive Medicine – Proceedings. Tallinn, Estonia. 21–23 March. 2018.
Robinson SE, Christley RM. Exploring the role of auction markets in cattle movements within great Britain. Prev Vet Med. 2007;81:21–37.
Bigras-Poulin M, Thompson RA, Chriél M, Mortensen S, Greiner M. Network analysis of Danish cattle industry trade patterns as an evaluation of risk potential for disease spread. Prev Vet Med. 2006;76:11–39.
Natale F, Giovannini A, Savini L, Palma D, Possenti L, Fiore G, Calistri P. Network analysis of Italian cattle trade patterns and evaluation of risks for potential disease spread. Prev Vet Med. 2009;92:341–50.
Clegg TA, Good M, Duignan A, Doyle R, More SJ. Shorter-term risk of Mycobacterium bovis in Irish cattle following an inconclusive diagnosis to the single intradermal comparative tuberculin test. Prev Vet Med. 2011;102:255–64.
Clegg TA, Good M, Duignan A, Doyle R, Blake M, More SJ. Longer-term risk of Mycobacterium bovis in Irish cattle following an inconclusive diagnosis to the single intradermal comparative tuberculin test. Prev Vet Med. 2011;100(3):147–54.
Lahuerta-Marin A, Gallagher M, McBride S, Skuce R, Menzies F, McNair J, McDowell SW, Byrne AW. Should they stay, or should they go? Relative future risk of bovine tuberculosis for interferon-gamma test-positive cattle left on farms. Vet Res. 2015;46(1):90.
Mazeri S, Sargison N, Kelly RF, Barend M, Handel I. Evaluation of the performance of five diagnostic tests for Fasciola hepatica infection in naturally infected cattle using a Bayesian no gold standard approach. PLoS One. 2016;11(8):e0161621.
We would like to acknowledge the work of TB teams within the veterinary service of Northern Ireland, and helpful discussions with Dr. Angela Lahuerta-Marin at AFBI.
This work was funded as part of a broader study of concurrent infections and their role in bTB epidemiology in Northern Ireland (48005 (122035): 15/3/10 − bTB Endemic Diseases; PI: Dr. Andrew Byrne) funded under the Evidence and Innovation (E&I) Scheme by the Department of Agriculture and Rural Development (DAERA; www.daera-ni.gov.uk).
Availability of data and materials
The data that support the findings of this study are available from DAERA but restrictions apply to the availability of these data, which were used under license for the current study, and so are not publicly available (due to data protection confidentiality). Data are however available from the authors upon reasonable request and with permission of DAERA (www.daera-ni.gov.uk/access-information-0).
Ethics approval and consent to participate
All ante-mortem testing and post-mortem samples were obtained in compliance with official guidelines for the control and eradication of bovine TB in Northern Ireland (Department of Agriculture, Environment, and Rural Affairs (DAERA)), in accordance with legislation (Tuberculosis Control order (NI) 1999 No.263; Tuberculosis (Examination & Testing) Scheme order (NI) 1999 No.264) and under a DSA with the data controllers (DAERA; Data Protection Act 1998).
The authors declare that they have no competing interests.
Consent for publication
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Graphical exploratory assessment of the relationship between the probability of bovine TB visible lesion (VL) presence and the bovine minus avian tuberculin (b-a) reaction sizes (mm) at the disclosing test for cattle culled in Northern Ireland during bTB breakdowns (n = 5698). The LOWESS curve (grey solid line) is a locally weighted regression line (bandwidth: 0.5); orange line represent b-a as a log-transformed predictor. Note, that the log-transformed predictor fails to address the increasing risk at negative b-a values, due to this the b-a variable was modelled as a categorical variable. (PDF 159 kb)
Detailed description of univariable associations between candidate predictors of TB-visible lesion presence. (PDF 333 kb)
Relationship between the probability of an animal having a bTB visible lesion at slaughter and the time between the penultimate and disclosing skin tests in cattle in Northern Ireland. Orange line represents a locally weighted regression line (non-linear regression fit; LOWESS); Grey line represents the predicted fit from the log-transformed predictor variable (log(time in days)). (PDF 157 kb)
The point estimated linear predictions from a random effects logit regression model for the presence of bTB visible lesions (VL) at slaughter amongst recorded cattle breeds for Mycobacterium bovis exposed cattle slaughtered in Northern Ireland. The grey dots represent point estimates; dots farther from zero line represent animals with higher or lower risk of disclosing with a lesion at slaughter relative to the population average. Orange circles represent weighted sample sizes. Outliers include (grey arrows): BGA = Belted Galloway highest risk of VL; SHB = Short horn beef cattle lowest risk of VL. Highly represented breeds are highlighted. (PDF 239 kb)
About this article
Cite this article
Byrne, A.W., Graham, J., Brown, C. et al. Bovine tuberculosis visible lesions in cattle culled during herd breakdowns: the effects of individual characteristics, trade movement and co-infection. BMC Vet Res 13, 400 (2017) doi:10.1186/s12917-017-1321-z
- Bovine TB
- Mycobacterium bovis
- Liver fluke
- Johne’s disease
- Gamma interferon
- Skin tests