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

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