Age-dependent patterns of bovine tuberculosis in cattle
© Brooks-Pollock et al.; licensee BioMed Central Ltd. 2013
Received: 25 February 2013
Accepted: 30 September 2013
Published: 16 October 2013
Bovine tuberculosis (BTB) is an important livestock disease, seriously impacting cattle industries in both industrialised and pre-industrialised countries. Like TB in other mammals, infection is life long and, if undiagnosed, may progress to disease years after exposure. The risk of disease in humans is highly age-dependent, however in cattle, age-dependent risks have yet to be quantified, largely due to insufficient data and limited diagnostics. Here, we estimate age-specific reactor rates in Great Britain by combining herd-level testing data with spatial movement data from the Cattle Tracing System (CTS). Using a catalytic model, we find strong age dependencies in infection risk and that the probability of detecting infection increases with age. Between 2004 and 2009, infection incidence in cattle fluctuated around 1%. Age-specific incidence increased monotonically until 24–36 months, with cattle aged between 12 and 36 months experiencing the highest rates of infection. Beef and dairy cattle under 24 months experienced similar infection risks, however major differences occurred in older ages. The average reproductive number in cattle was greater than 1 for the years 2004–2009. These methods reveal a consistent pattern of BTB rates with age, across different population structures and testing patterns. The results provide practical insights into BTB epidemiology and control, suggesting that targeting a mass control programme at cattle between 12 and 36 months could be beneficial.
Bovine tuberculosis (BTB) is an infectious disease of cattle caused by the TB species Mycobacterium bovis. In Great Britain, the number of cattle slaughtered due to BTB increased from ~200 in 1986 to almost 30 000 in 2010, despite extensive and expensive control measures . Like TB in other mammals, disease is characterised by granulomas that form in the respiratory system from which infectious mycobacteria are excreted [3, 4]. However, the focus for control is on the removal of infected cattle at an earlier stage of infection, detected by a positive sensitivity reaction to the single intradermal comparative cervical tuberculin (SICCT) test.
Cattle herds in Great Britain are subject to compulsory regular testing for the purposes of both surveillance and control . Herds in the highest incidence areas have all cattle over 6 weeks tested annually, whereas herds in low incidence areas must test all cows and male breeders every four years . Reactor cattle (cattle that react to the SICCT test) are slaughtered and the breakdown herd is subject to movement restrictions until a number of follow-up tests are passed. Although the surveillance system generates a wealth of data, it has not been possible (until very recently) to estimate age-specific BTB risks because negative results were systematically recorded at the herd level only.
Quantifying the age-specific risk is a key tool for estimating critical vaccination thresholds and the basic reproductive ratio of a disease, R0[7, 8]. Once an epidemic has reached a steady state, R0 is given by the ratio of the distribution of life times to the distribution of ages at first infection . The theory behind disease control is to reduce R0, which results in an older average age of infection. An age-varying risk of infection can be exploited for control purposes to reduce the effort required for eradication . For instance, if infection rates in calves were higher than in older cattle, then increasing the average age of infection would have a bigger impact than expected under the assumption of a risk that is independent of age.
Age-stratified serological surveys have been used to estimate the age-specific risk of infection for many childhood diseases . Age-stratified case reports can also be used, provided that two conditions hold: that all individuals experience infection during their lifetimes and that there are no age-related differences in detection . For BTB in Great Britain, only a small percentage of animals are infected during their lifetimes and the age-related tested patterns are unknown, so that without additional data, reactor reports alone cannot be used to quantify the risk of infection with age. In this paper, we combine herd-level test data and spatial movement data from the Cattle Tracing System (CTS) to estimate background testing patterns in order to provide estimates of age- and time-dependent BTB risks in cattle. We use the age-specific reactor rates to parameterise a general age-structured catalytic model of BTB infection in order to explore the patterns with age and impact of testing.
Materials and methods
Great Britain has a rich BTB dataset dating back to the 1950s when the first test-and-slaughter scheme was introduced to control disease . Herd-level test results for test-negative herds and animal-level test results for reactor cattle, inconclusive reactor cattle and tests resulting from contact tracing are contained in the database VetNet, collated and managed by the Animal Health and Veterinary Laboratories Agency (AHVLA), which is part of the UK department for Food, the Environment and Rural Affairs (Defra).
Cattle demographic and movement data are contained within Cattle Tracing System (CTS). Introduced in 1996, the CTS contains the births, movements and deaths of all registered cattle in Great Britain, cattle data (sex, breed) and location data . The CTS is run by the British Cattle Movement Service (BCMS), also part of Defra. We used an extract from VetNet covering BTB tests between 2004 and 2009 and a CTS extract for the same period. The data were provided by Defra via the AHVLA and RADAR (Rapid Analysis and Detection of Animal-related Risks). We collated the data using PostgreSQL  and carried out the modelling in R .
Herd test types
All herds in GB are subject to regular SICCT testing at a frequency determined by the local incidence of infection [6, 15]. Herds in annual testing areas have a whole herd test (WHT) once a year of all cattle over 6 weeks that are present on the date of the test. Herds in four-yearly testing areas are subject to a routine herd test (RHT) once every 4 years of female animals that have calved, male breeders and male calves intended for breeding that have been purchased since the last test. Herds with confirmed infection at a surveillance test are subject to two short interval (SI) tests at 60 days and 120 days, then additional follow-up tests after 6 months (6M) and twelve months (12M). For a more detailed description of the testing programme, see [6, 15].
The number of reactors was collated directly from the “animal” table in VetNet. Reactor cattle were identified by the test result “R”. Inconclusive reactors or other results were not included. Reactors were classified by age in whole years (test date minus birth date), sex (M or F), calendar year of test, purpose (beef, dairy or mixed) based on breed and test type (routine herd tests (RHT), whole herd tests (WHT), short-interval tests following a breakdown (SI) and follow-up tests at 6 months (6M) and 12 months (12 M) after passing breakdown release).
Inferring negative tests
In general, there is no historic information in VetNet about cattle that tested negative to the SICCT test. In order to estimate underlying testing patterns, we reconstructed the herd on the test date by combining herd-level testing data from VetNet with animal-level data from the CTS. From the list of cattle that were present during a test, we used Defra’s eligibility guidelines to determined which of the cattle had been eligible for the test. We assumed that all eligible present animals were tested. In detail, we extracted the date on which herds were tested and the type of test that was carried out from the VetNet testing table. If testing took place over several days and was recorded multiple times in VetNet, we took the date of the first day of the test. Then the CPH (county-parish-holding) number was matched to the CTS locations table to obtain a CTS id number. Using the CTS id number, we extracted a list of cattle id numbers that had moved onto the premises before the test date and left after the test date (or had not left at all) from the CTS movement table. The cattle id numbers were used to determine the sex and age at the time of the test. We determined eligibility according to the type of test that was carried out:
Whole herd tests (WHT) and follow-up tests (SI, 6M, 12M): all animals aged over 6 weeks.
Routine herd tests (RHT): female animals that have calved, male breeders, male calves intended for breeding that have been purchased since the last test.
Female animals that had calved and male breeders were identified using the relationships data in the standard CTS extract and, as the relationships data were not complete, by the age of the animal. It was difficult to identify male calves intended for breeding, so instead we included male animals that had yet to be tested in that herd. Between 2004 and 2009, 700 979 out of 701 577 VetNet records matched a CTS record. In terms of premises, 87 177 distinct CPH numbers were recorded as having a test between 2004 and 2009. We were able to match 92,251 (99.6%) in the CTS location table. Of the cattle records, 1.4% were removed due to missing birth dates.
Once a list of cattle had been obtained, we classified the tested cattle (as for reactors) by age, sex, purpose and tests by the test type and calendar year.
Because the denominator is the number of cattle tested not the total cattle population, the reactor rates are effectively weighted towards areas with higher levels of testing. Given the imperfect sensitivity and specificity of the SICCT test, the reactor rates only approximate true prevalence. However, as reactor rates are relative low (≤ 2.5%), the discrepancy is negligible and does not impact on the overall results .
Other demographic data
Collated mortality data were extracted from the CTS livestock table by calendar year, age, sex and breed purpose (beef or dairy).
Full derivations of equations (6) and (7) are in the Additional file 1 and further details of similar derivations can be found in  (Appendix D, pages 677–678), with the difference that we explicitly include mortality as the total population size decreases due to loss of reactors. Initial conditions are S0 = 1, I0 = 0, and for the other age groups the initial conditions are computed from the equations.
We investigated two alternative formulations of λ i : one where λ i is independent of age but dependent on calendar year such that λ iy = β y , and one in where the effects of age and year are separable and distinguishable, such that λ iy = α i β y .
From these age-specific equations, we calculate the average age at infection as and the average age at detection as . A lower bound for the average number of secondary cases per infected animal (Rƒ) is given by the reciprocal at the proportion of susceptible cattle, calculated as .
where p iy = γI iy /(S iy + I iy ) is calculated using the model (equations 6 and 7) as a function of λ iy . To improve convergence for λ, we re-parameterised the model to estimate log (λ) and sampled from a flat prior log([1e− 10, 1e− 1]). We used the optim function to obtain maximum likelihood estimates for the parameters that were used as starting points for the MCMC chain. Following visual inspection of the MCMC trace plots, we generated 10 000 samples of the posterior distributions from a chain of 40 000 samples (i.e. sampling every 4th sample) and with a burn-in of 10 000 samples. We checked convergence and mixing with a visual inspection of the MCMC trace plots and model fit and by comparing multiple MCMC chains.
Where n is the number of parameters and where k is the number of observations and log (L) is the mean log likelihood of the MCMC chain. We identified the model with the largest explanatory power by the lowest BIC score.
Impact of removal rate
Infected cattle are removed in the model via testing and death only (there is no recovery or death due to infection). In fitting the model, we found that the removal rate was highly correlated with the force-of-infection experienced by the first cohort and therefore could not be estimated independently. Using an unconstrained prior allowed solutions in which the removal rate was close to zero, which is not believed to be the case [19, 20]. Therefore, we used an informative prior for the removal rate with distribution ~ N (0.6, 0.02)which led to posterior estimates ranging from 50% to 70%. Furthermore, we tested the sensitivity of the results to the removal rate. As is discussed in the results, the removal rate affects the magnitude of the λ i s for young cattle but not the patterns with age. For model simplicity, we did not include variable sensitivity and specificity associated with different test interpretations [19, 21], age-specific test sensitivity or any desensitisation associated with repeated testing.
Demography and patterns of cattle testing
Between 2004 and 2009, the number of cattle tests increased by 28%. Combining VetNet and the CTS, we find that the age distribution of tests remained largely consistent over that period, with on average 66% of tests conducted on cattle 3 years or younger (Additional file 2). In 2005, the number of tests of beef cattle exceeded the number of tests on dairy cattle and remained higher for the following years. Beef cattle are on average tested at younger ages compared to dairy cattle (37 months versus 48 months). We found that January born animals of all breeds were tested at older ages than animals born in other months (50 months versus 41 months).
Age-specific reactor rates
Although the number of tests increased year on year between 2004 and 2009, we find elevated reactor rates in 2005, 2008 and 2009. The lowest reactor rates were in 2006 when we estimate that 0.73% of tests were failed tests. Although there were over 7000 more reactors in 2009 compared to 2005, we estimate that the reactor rates were equal up to two significant figures (1.12% in 2005 compared to 1.10% in 2009).
Age-specific catalytic model
Fitting the λ jy s independently, we estimate a force-of-infection ranging between 0.003 years-1 (beef cattle aged 3–4 years in 2006) and 0.036 years-1 (dairy cattle aged 2–3 years during 2008). For both dairy and beef cattle and across all years, the highest rates of infection were in cattle aged 2–4 years.
Infection rates for cohort 0 (0 to 1 years) were similar for dairy and beef cattle across all years, suggesting that these youngest cattle experience similar risks despite being located in different parts of the country and in different herd structures. Infection rates for cohort 0 increased from 0.0083 years-1 in 2004 to a maximum of 0.013 years-1 in 2008. The infection risk in these youngest cattle may reflect true infection annual risk, as they will not have experienced infection in previous years.
In cohort 4, at 3 years, the number of cattle removed as reactors exceeds the number of new infections generated. In ages up to 3 years, the imperfect removal rate results in infected cattle left behind in the population. After 3 years, the number of new infections is less than the number of reactors and residual infection is cleared. The differential impact of removal rate with age is reflected in the 2d posterior density plots (Additional file 6). In cohort 0, there is a strong co-dependence between infection rate and removal rate, where higher removal rates are balanced by lower infections rates, and vice versa, higher infection rates can be balanced by a lower removal rate to produce equally good model fits. However, after 3 years, the trade-off between infection and removal rate is no longer observed and the infection rates are well defined.
The age-dependent risks were different for dairy and beef cattle. In dairy cattle (Figure 3b), all ages had an increased risk of reacting relative to cohort 0, with the greatest relative risk experienced by cattle aged 1–2 years. In beef cattle (Figure 3e), cattle aged 1 to 2 years experienced a greater risk, but then cattle over 3 years experienced a lower risk than cattle under 1 year.
These estimates for force-of-infection can be interpreted in terms of cattle ages. We estimate that the mean age of infection in dairy cattle was 4.4 (3.9, 4.6) years. The mean age at detection over the same period was 8 months later at 5.1 years. The mean percentage of susceptible dairy cattle was 97.9% (calculated using equation 6), which suggests that the reproductive number in the population was not less than 1.03. In beef cattle, the mean age at detection was 3.8 years, but we estimate that the mean age at infection was 3.3 years. In beef cattle, the mean proportion of susceptible animals was 98.6%, again suggesting a reproductive number of not less than 1.015.
To our knowledge, this work is the first attempt to characterize the age-specific dynamics of TB in cattle. We demonstrated that the age-specific reactor rates cannot be explained by a cumulative constant hazard over an animal’s lifetime. Knowledge of age-specific risks is essential for accurately interpreting trends in pooled epidemiological data and for the design and assessment of control strategies such as vaccination.
The method we describe here provides a way of estimating age-specific reactor rates using herd-level testing data combined with cattle movement records and can be applied at a local or national scale. Combining the detailed spatial demography contained in the Cattle Tracing System with the national epidemiological data in VetNet has provided insight into BTB epidemiology and control in Great Britain, although there is an important distinction between test results and true infection rates. As the majority of tests are conducted in high incidence areas in Britain, the results presented here reflect disease dynamics in high incidence regions. In low incidence settings, beef cattle intended for slaughter are not routinely tested, and therefore are not fully captured here. Mitchell et al. demonstrated that many cattle, particularly in low incidence areas, are never tested during their lifetimes because of turnover and movement between herds . Indeed, this is reflected by the low numbers of young animals tested during routine herd tests. We believe that the increase in incidence with age reflects an intrinsic property of bTB transmission. However, to estimate the magnitude of the force-of-infection experienced by a different population, it would be necessary to measure both reactor numbers and underlying testing patterns with age.
This analysis is based on routinely collected surveillance data does not address the underlying immunology of BTB infection, nor attempt to capture the complex BTB transmission process. Our motivation was to obtain a parsimonious interpretation of the available data without assumptions about the relationship between test status and infectiousness or the interaction between infected cattle, infected wildlife and the local environment. Previous approaches to modelling BTB have included multiple stages of infection [23–25], including a latent or occult period immediately following infection when infected animals may be unresponsive to the skin test. Including an occult period in our model formulation would increase the estimated infection risks as a higher proportion of infected cattle would die before being detected. The natural next step for this work is to include an explicit transmission process to capture the contact process between susceptible and infectious cattle in a “who acquires infection from whom?” framework. We did not include an explicit transmission process in this model because of the complex and fiercely debated role of badgers in cattle transmission dynamics.
Nevertheless, the substantial variation of BTB rates with age has implications for both the understanding of BTB infection and the design and implementation of control programmes. With TB in humans, the tuberculin skin test has a lower sensitivity for disease in young children than in other age groups  and it has been suggested that vaccinating neonates may be suboptimal due to differential response to infection. Comparison can also be drawn with Mycobacterium avium subspecies paratuberculosis (MAP) infection in cattle. Response to MAP has been shown to vary with age with young calves at increased risk of developing high bacterial loads [27, 28]. Extrapolating from TB and MAP to BTB, it could be that young calves have a lower chance of detection but have the potential to be high shedders. Further analysis of VetNet indicates that 50-60% of reactors under 2 years are found with visible lesions, compared to around 30% in older reactors (Figure 4). Further work is needed to quantify the role of calves and older cattle in the transmission process.
Many studies have found differences in BTB risk between dairy and beef cattle [29–31], and indeed the unadjusted age profiles of beef and dairy reactor numbers are hugely different (Figure 1b). However, we found that when reactor numbers are normalised by test patterns the age dependencies are remarkably comparable. In fact, we found that beef calves under 1 year experience similar risks of reacting to dairy cattle of the same age. This highlights interesting questions about the drivers of the differences between older beef and dairy cattle and the interaction between breed, herd structure and age. Could the differential risks in age and breed be exploited for improved control? A detailed understanding of age-dependent patterns in different countries with different infection burdens will help to elucidate the true characteristics of BTB infection in cattle.
- SICCT test:
Single intradermal comparative cervical tuberculin
Cattle Tracing System.
EBP is funded by an Engineering and Physical Sciences Research Council (EPSRC) fellowship. JLNW is supported by the Alborada Trust and the RAPIDD program of the Science & Technology Directorate, U.S. Department of Homeland Security, and the Fogarty International Center, U.S. National Institutes of Health. AJKC is supported by Defra grant no. SE-3127. TJM is supported by the BBSRC. We thank Steve Holdship and Rose Nicholson at Defra the AHVLA team for providing access to the CTS and VetNet.
- Smith NH, Hewinson RG, Kremer K, Brosch R, Gordon SV: Myths and misconceptions: the origin and evolution of Mycobacterium tuberculosis. Nat Rev Microbiol. 2009, 7: 537-544. 10.1038/nrmicro2165.View ArticlePubMedGoogle Scholar
- Animal Health and Veterinary Laboratories Agency (AHVLA): Bovine tuberculosis: Infection status in cattle in GB. 2012, London: Statistical ReportGoogle Scholar
- Bloom B, Murray C: Tuberculosis: commentary on a reemergent killer. Science. 1992, 257: 1055-1064. 10.1126/science.257.5073.1055.View ArticlePubMedGoogle Scholar
- Menzies FD, Neill SD: Cattle-to-cattle transmission of bovine tuberculosis. Vet J. 2000, 160: 92-106.View ArticlePubMedGoogle Scholar
- Abernethy DA, Upton P, Higgins IM, McGrath G, Goodchild AV, Rolfe SJ, Broughan JM, Downs SH, Clifton-Hadley R, Menzies FD, De la Rua-Domenech R, Blissitt MJ, Duignan A, More SJ: Bovine tuberculosis trends in the UK and the Republic of Ireland, 1995–2010. Vet Rec. 2013, 172: 312-10.1136/vr.100969.View ArticlePubMedGoogle Scholar
- Krebs J, The Independent Scientific Review Group: Bovine TB in Cattle and Badgers. 1997, London: Department for the Environment Food and Rural Affairs (Defra), Report PB3423Google Scholar
- Anderson RM, May RM: Age-related changes in the rate of disease transmission: implications for the design of vaccination programmes. J Hyg (Lond). 1985, 94: 365-436. 10.1017/S002217240006160X.View ArticleGoogle Scholar
- Anderson RM, May RM: Infectious Diseases of Humans. Dynamics and Control. 1992, Oxford: Oxford University PressGoogle Scholar
- Dietz K: Transmission and control of arbovirus diseases. Epidemiology. Edited by: Lugwig D, Cooke KL. 1975, Philadelphia: Society for Industrial and Applied Mathematics, 104-121.Google Scholar
- Farrington CP, Kanaan MN, Gay NJ: Estimation of the basic reproduction number for infectious diseases from age-stratified serological survey data. Appl Stat. 2001, 50: 251-292.Google Scholar
- Reynolds D: A review of tuberculosis science and policy in Great Britain. Vet Microbiol. 2006, 112: 119-126. 10.1016/j.vetmic.2005.11.042.View ArticlePubMedGoogle Scholar
- Mitchell A, Bourn D, Mawdsley J, Wint W, Clifton-Hadley R, Gilbert M: Characteristics of cattle movements in Britain–an analysis of records from the cattle tracing system. Anim Sci. 2005, 80: 265-273.View ArticleGoogle Scholar
- The PostgreSQL Global Development Group: PostgreSQL Database Management System. version 9.0
- R Core Team: R: A language and environment for statistical computing. 2012, Vienna, Austria: R Foundation for Statistical ComputingGoogle Scholar
- Green LE, Cornell SJ: Investigations of cattle herd breakdowns with bovine tuberculosis in four counties of England and Wales using VETNET data. Prev Vet Med. 2005, 70: 293-311. 10.1016/j.prevetmed.2005.05.005.View ArticlePubMedGoogle Scholar
- Rogan WJ, Gladen B: Estimating prevalence from the results of a screening test. Am J Epidemiol. 1978, 107: 71-76.PubMedGoogle Scholar
- Vynnycky E, White RG: An Introduction to Infectious Disease Modelling. Oxford University Press. 2010, Oxford: Oxford University PressGoogle Scholar
- Martin AD, Quinn KM, Park JH: MCMCpack: Markov Chain Monte Carlo in R. J Stat Softw. 2011, 42: 22-View ArticleGoogle Scholar
- Conlan AJK, McKinley TJ, Karolemeas K, Brooks-Pollock E, Goodchild AV, Mitchell AP, Birch CP, Clifton-Hadley RS, Wood JL: Estimating the hidden burden of bovine tuberculosis in Great Britain. PLoS Comput Biol. 2012, 8: e1002730-10.1371/journal.pcbi.1002730.PubMed CentralView ArticlePubMedGoogle Scholar
- Monaghan ML, Doherty ML, Collins JD, Kazda JF, Quinn PJ: The tuberculin test. Vet Microbiol. 1994, 40: 111-124. 10.1016/0378-1135(94)90050-7.View ArticlePubMedGoogle Scholar
- De la Rua-Domenech R, Goodchild AV, Vordermeier HM, Hewinson RG, Christiansen KH, Clifton-Hadley RS: Ante mortem diagnosis of tuberculosis in cattle: a review of the tuberculin tests, gamma-interferon assay and other ancillary diagnostic techniques. Res Vet Sci. 2006, 81: 190-210. 10.1016/j.rvsc.2005.11.005.View ArticlePubMedGoogle Scholar
- Mitchell AP, Green LE, Clifton-Hadley R, Mawdsley J, Sayers R, Medley GF: An analysis of single intradermal comparative cervical test (SICCT) coverage in the GB cattle population. 2006, Exeter: Proceedings of the Society for Veterinary Epidemiology and Preventive Medicine (SVEPM), 70-86.Google Scholar
- Barlow ND, Kean JM, Hickling G, Livingstone PG, Robson AB: A simulation model for the spread of bovine tuberculosis within New Zealand cattle herds. Prev Vet Med. 1997, 32: 57-75. 10.1016/S0167-5877(97)00002-0.View ArticlePubMedGoogle Scholar
- Kao RR, Roberts MG, Ryan TJ: A model of bovine tuberculosis control in domesticated cattle herds. Proc Biol Sci. 1997, 264: 1069-1076. 10.1098/rspb.1997.0148.PubMed CentralView ArticlePubMedGoogle Scholar
- Fischer EA, Van Roermund HJ, Hemerik L, Van Asseldonk MA, De Jong MC: Evaluation of surveillance strategies for bovine tuberculosis (Mycobacterium bovis) using an individual based epidemiological model. Prev Vet Med. 2005, 67: 283-301. 10.1016/j.prevetmed.2004.12.002.View ArticlePubMedGoogle Scholar
- Moyo S, Isaacs F, Gelderbloem S, Verver S, Hawkridge AJ, Hatherill M, Tameris M, Geldenhuys H, Workman L, Pai M, Hussey G, Hanekom WA, Mahomed H: Tuberculin skin test and QuantiFERON® assay in young children investigated for tuberculosis in South Africa. Int J Tuberc Lung Dis. 2011, 15: 1176-1181. 10.5588/ijtld.10.0770.View ArticlePubMedGoogle Scholar
- Nielsen SS, Ersbøll AK: Age at occurrence of Mycobacterium avium subspecies paratuberculosis in naturally infected dairy cows. J Dairy Sci. 2006, 89: 4557-4566. 10.3168/jds.S0022-0302(06)72505-X.View ArticlePubMedGoogle Scholar
- Mitchell RM, Medley GF, Collins MT, Schukken YH: A meta-analysis of the effect of dose and age at exposure on shedding of Mycobacterium avium subspecies paratuberculosis (MAP) in experimentally infected calves and cows. Epidemiol Infect. 2012, 140: 231-246. 10.1017/S0950268811000689.View ArticlePubMedGoogle Scholar
- Ramírez-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-198. 10.1016/j.prevetmed.2009.08.011.View ArticlePubMedGoogle Scholar
- Skuce RA, Allen AR, McDowell SWJ: Herd-level risk factors for bovine tuberculosis: a literature review. Vet Med Int. 2012, 2012: 621210-PubMed CentralView ArticlePubMedGoogle Scholar
- Karolemeas K, McKinley TJ, Clifton-Hadley RS, Goodchild AV, Mitchell A, Johnston WT, Conlan AJK, Donnelly CA, Wood JLN: Predicting prolonged bovine tuberculosis breakdowns in Great Britain as an aid to control. Prev Vet Med. 2010, 97: 183-190. 10.1016/j.prevetmed.2010.09.007.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.