Risk factors for bovine Tuberculosis at the national level in Great Britain

The control of an infectious disease is often complicated when more than one host species are involved and one of those is a wildlife reservoir. This problem is exacerbated if the reservoir host is poorly understood or if the disease is difficult to manage in the reservoir population [1, 2]. One of the best described examples of this is that of bovine Tuberculosis (bTB) in Great Britain (GB), caused by Mycobacterium bovis. Largely eliminated in many industrialised countries, in GB the disease was eliminated from most cattle herds during the 1960s [3]. However, since then prevalence has been gradually rising and the area affected growing in size [4]. It was soon suspected that this rise was related to the presence of the Eurasian badger (Meles meles), once it was recognised that badgers were a host for M. bovis[5, 6]. However, the nature of the relationship to cattle bTB remained controversial and so the randomised badger culling trial (RBCT) was set up to investigate the effectiveness of badger culls as a control for bTB spread [7]. While the RBCT showed that badger culling operations and cattle herd breakdowns rates are strongly and significantly associated with each other [8–10], relatively little is understood about the role of badgers in herd breakdowns across GB, though a recent paper has analysed risk factors for badger presence [11]. These analyses have been extrapolated to produce a prediction of badger distribution in GB, features that have previously been suggested could be used to predict risk of bTB to cattle [12].

Bovine TB cases are identified using the Single Intradermal Comparative Cervical Test (SICCT) or “skin test” which has relatively low sensitivity [13]. An animal reacting positively to the skin test is called a reactor and is slaughtered. Further, limited ante-mortem testing is carried out using gamma-interferon test, which is believed to have greater sensitivity [14]. Furthermore, carcasses of slaughtered animals are inspected at the abattoir for lesions that could indicate bTB. The case is confirmed by successfully culturing the pathogen from a suspect lesion post-mortem. Routine testing is carried out using the skin test on a one-to-four yearly basis depending upon local disease incidence. Additional testing is carried out following a herd breakdown, on neighbouring herds and on herds with traced contacts with breakdown herds. Cattle moved from high risk areas are also tested before movement to other areas; further, all cattle from England and Wales are tested before and after movement to Scotland.

The number of detected bTB breakdowns in GB has been rising year-on-year for a number of years, punctuated by a decrease in cases in 2001 caused by the epidemic of foot-and-mouth disease (FMD) in GB that restricted the number of tests that could be conducted that year [15]. The culmination of the year-on-year rise in incidence was in 2008 (the peak in incidence), when 2,764 herd breakdowns were confirmed (a prevalence of 5%) with a total of 39,302 animals slaughtered as reactors, and with movement restrictions imposed upon 7,957 (9.2%) herds [4]. Although the pattern of increase and spread in bTB is well documented, the reasons for this increase remain unclear.

While badgers are now recognised as a critical component for bTB persistence, the contribution of badgers to increasing incidences of bTB has not been quantified. The restocking of herds following the 2001 FMD epidemic has been demonstrated to have facilitated the spread of bTB to new areas [16]. However, restocking caused by FMD only explains a small part of the bTB expansion and various studies have looked into risk factors for herd breakdowns to assist understanding of the expansions of these hotspots [15, 17–20]. There are a number of potential routes of transmission and previous studies [15, 18] have demonstrated that cattle movements are important for transmission, with an estimated 15% of herd breakdowns are due to recorded movements of cattle, 75% of breakdowns to local effects within High Risk areas (HRAs), and 9% unexplained [18]. These differences will vary with the local disease prevalence, with movements and the unexplained component being most important in three and four year testing areas. In addition to movements of potentially exposed animals within GB, there are a large number of imports of cattle from Ireland. Despite post-import testing of animals arriving from Ireland, this is likely to represent a risk of introduction of disease from outside GB due to the arrival of falsely negative animals [17].

Specific risk factor studies have identified larger herds as both a risk factor for having a breakdown [20–22] and for the duration of the breakdown [23]. Movements of animals to and from high incidence areas [16, 24], and from farm sales [19, 21] have also been demonstrated to be associated with breakdowns. Other factors associated with herd breakdowns include being a dairy farm and various farming practices such as operating over multiple premises, the spreading of fertilizers and the area of land farmed [7, 16, 19, 20, 25]. However, all of these analyses have been applied to small portions of high incidence areas. Thus far, no studies have looked at risk factors at the national level and in particular none have compared them to those in lower incidence areas.

The extent of the bTB endemic area in GB continues to expand as identified by changes in the parish testing intervals (Figure 1). Under the bTB control policy in operation prior to 2009, all herds are routinely tested on a one, two, three or four year routine herd testing (RHT) interval based upon the perceived risk in the parish of the farm, or the risk in neighbouring parishes. The testing period of the parish is reviewed on a quarterly basis and the extent of the one and two yearly tested parishes are considered high risk areas (HRAs), as they are intended to reflect the current extent of bTB endemicity.

A total of 15,358 bTB breakdowns were confirmed between 2002 and 2008 inclusive, of which 11,599 were not recorded as part of a follow-up and so can be considered a new breakdown for the purposes of these analyses (Table 1). The number of herds in GB is declining, whilst the number of breakdowns is increasing (Figure 2). Furthermore, the number of herds in 1 year testing parishes (high incidence parishes) is increasing at the expense of those in 4 year testing parishes (Figure 2). Analysis of the period of time spent in each parish testing interval shows that the majority of herds in 4 year testing intervals are there for the duration of the study period (Figure 3). However, intervals 2 and 3 are to a much greater extent a transitional testing interval, with the majority of herds spending only one or two years in 2 and 3 year testing intervals (Figure 3).

		Testing interval
	Unit	All herds		1 year		2,3 year		4 year
		Non-TB	TB	Non-TB	TB	Non-TB	TB	Non-TB	TB
Country	England	373849	9068	91768	6607	52630	1491	187359	284
	Wales	102840	2400	25107	1539	20908	539	95854	131
	Scotland	95854	131	0	0	0	0	34143	52
Herd type	Beef	213231	3566	44569	2649	24199	552	124649	132
	Fatt’	58847	1163	10286	764	6133	181	37872	115
	Suckler	167501	2770	31721	1927	22144	515	92457	102
	Dairy	113234	3807	26576	2618	18435	728	51255	96
	Store	19730	293	3723	188	2627	54	11123	22
Herd size	0 – 10	137438	427	27406	320	17158	61	78065	18
	11 – 100	242294	4021	51514	2984	31853	660	130805	90
	>100	192811	7151	37955	4842	24527	1309	108486	359
Irish imports	0	567291	11505	116560	8116	73312	2019	312810	423
	>0	5252	94	315	30	226	11	4546	44
History	No	556467	9744	103245	6430	72053	1923	317053	464
	Yes	16076	1855	13630	1716	1485	107	303	3
High risk moves	0	290232	393	15031	164	9958	39	246715	147
	1-10	102995	923	22602	622	14758	122	45405	62
	>10	179316	10283	79242	7360	48822	1869	25236	258
FMD	No	561173	11296	114834	7958	72219	1970	310725	448
	Yes	11370	303	2041	188	1319	60	6631	19
Badgers	Mean prob'	0.259	0.281	0.293	0.295	0.283	0.286	0.223	0.219
Block Totals		572543	11599	116875	8146	73538	2030	317356	467

The size of the local cattle population was a significant risk factor as a univariate predictor, however once incorporated in the model along with the herd size category, the sign changes and the predictor becomes a significant protective effect. This is due to a correlation between the number of cattle in the 100 km² cells and the herd size covariate, due to a clustering of larger herds. There was an additional correlation between badger numbers and the local cattle density (R² = 0.244, p < 0.001). As the number of animals on a herd usually changes during the year, sensitivity analysis compared animal numbers on 1^st January with those on 1^st June and showed no significant change to risk factor estimates. The estimated probability of badger occurrence was a significant risk factor in the highest bTB incidence areas (1 year testing) and the history of bTB was a significant risk factor in 1, 2 and 3 year testing areas. Neither of these factors were significant in the lowest bTB incidence (four year testing) areas. Sensitivity analysis compared the implementation of history of infection here – infections between 1997 and 2001 versus a history of infections in the previous two years. There was no significant change in the odds ratios between the two implementations. Having had animals culled during the 2001 FMD epidemic was only significant in two and three year testing areas.

As the incidence of bTB continues to rise in GB it is important to develop an improved understanding of the characteristics of the herds that are becoming infected. This can aid control and surveillance of the disease to help to ensure that the disease is contained with minimum surveillance effort. The analyses presented in this paper have demonstrated that it is possible to identify herds that are likely to have a bTB breakdown. Variables associated with the risk of infection appeared to vary depending upon the local incidence of disease. In high incidence parishes (the one year testing parishes) it is features of the locality (badger density and the history of bTB) that define the likelihood of infection. In low-incidence areas (the four year testing parishes) the risk of infection is defined solely by herd demographics and whether the herd brings in high-risk animals from other areas of GB or from Ireland; animals from Ireland predominantly arrive in four year testing parishes (Table 2). The significance of these high risk imports underlines the need to a more sensitive test for bTB as these animals are tested prior to movement.

Whilst demographic factors are important in one year testing parishes, the risk of recording a breakdown is also defined by the probability of badger presence and whether there is a history of bTB on the herd between 1997 and 2001 – the period prior to the analysis in this study. This is the first time that badger-related factors have been shown to be important at a national scale. A sensitivity analysis that defined a history of infection as the two years prior to year k did not result is a significant change from analysis of the period 1997–2001. This indicates that the at-risk farms have been defined prior to the period of study and that the risk of infection is not dynamic but static. This opens up attractive options for risk-based surveillance of the at risk holdings by targeting such herds. In medium incidence (two and three year testing parishes) badgers are not significant as a risk factor, but a history of infection is. Additionally herds that had animals culled during the 2001 FMD epidemic are a significant risk factor in these areas, but not in others, which may be a result of where herds sourced their stock following the FMD epidemic. It may also be a result of disease being seeded in areas that had the properties to sustain infection over a prolonged period, or the effect being overwhelmed by an increase in incidence.

Only fattening herds were significantly different from the other herd types. There are several reasons for this. As fattening herds typically source only older stock and source them from a variety of sources they are likely to act as “sinks” for infection. Furthermore, as fattening herds send relatively large numbers of stock to the slaughterhouse, infection is more likely to be detected on fattening herds via abattoir surveillance and indeed, a greater proportion of breakdowns on fattening herds are identified at the slaughterhouse. Therefore the amount of testing is heterogeneous compounded by many fattening herds being exempt from routine herd testing due to their high stock turnover rate. However, for this reason they are also likely to be less important epidemiologically; the animals go to slaughter and as there are likely to be fewer movements off the fattening herds onto other herds, so fattening herds form the end of the transmission chain. The significance of fattening herds is in contrast with the results of Vial et al. [20] who found greater risk of breakdowns among dairy herds as well as larger herds. However, the study of Vial et al. [20] was restricted to much smaller areas and did not include movements of animals. This may explain the differences between the sets of results.

These results support the work of Green et al. [18] who demonstrated that 75% of cases are due to local effects (the additional local effects, such as badgers and the history of bTB, seen in high, and to a lesser extent, medium bTB incidence areas), whilst 15% are due to recorded cattle movements (dominating the risks seen in lower bTB incidence areas). The remaining 9% of unexplained breakdowns may be the result of heterogeneities in herd risk, reflecting the demographic factors seen here, in particular herd size and fattening herds and their potentially older age distribution.

Despite the disease being managed independently by the respective governments in England, Scotland and Wales, few differences were found between the three countries in terms of relative risk of infection. However, when just the four year-testing regions of GB are compared, there is a significantly greater risk of breakdown in Scotland relative to England and Wales. This is likely to be due to the significant protective effect of the y-coordinate in the model. There are large numbers of non-breakdown herds in four year testing areas of northern England, however, there have been a large cluster of breakdown herds in four year testing areas of southern Scotland, the significant difference of Scotland is therefore offsetting the effect of the y-coordinate protective effect in these areas. If the y-coordinate is dropped from the model, Scotland is no longer significant. Generally, despite the different prevalences and management strategies, there are no national level differences in bTB risks.

A key consideration when considering the role of wildlife reservoirs is whether or not the ‘spillover host’ (in this case, the cattle population) can sustain the disease epidemic on their own or produce transmission chains resulting in locally important problems. While it is known that transmission chains can occur [17], here we show that cattle alone are unlikely to sustain an epidemic under the existing test and slaughter regime; there is only a distinct and relatively stable subset of herds that can be defined as being at elevated risk of breaking down. In areas perceived to be at low risk, these factors are largely demographic rather than epidemiological. Indeed, whilst the factors that define the risk of infection vary depending upon the local prevalence of disease, the factors do not vary depending upon the country of GB. This shows that management factors are of secondary importance and that environmental factors are defining transmission and spread, though this analysis does not preclude within-herd transmission factors being associated with identified risks such as herd history. Therefore, control may be enhanced and streamlined by developing a risk-based strategy for bTB surveillance. In low risk areas the disease can be controlled by monitoring “risky herds” and potential high risk movements of animals and imports. However, this study has also demonstrated that there are no differences between the at-risk herd types in high risk areas compared with low risk areas; it is merely the environmental conditions that determine the risk in those areas. Therefore, should the specific environmental conditions become established in low risk areas then bTB may become established in these areas.

This paper has demonstrated that it is possible to identify cattle herds that are more likely to have a bTB breakdown, but the determinants of breakdown change as the background risk of infection changes. Within perceived higher risk areas it is principally locality factors that determine the risk of breakdown, whilst in lower risk areas the principal determinants are the amount of contacts that herds have with higher risk herds and animals. These determinants of infection could be used to develop a strategy for risk-based surveillance for bTB.