Event Abstract

To be or not to be (confirmed through bacteriology): spatial epidemiology of bovine tuberculosis in a low prevalence region in Spain

  • 1 VISAVET Health Surveillance Centre (UCM), Spain
  • 2 Other, Spain
  • 3 Ministry of Agriculture and Livestock, Junta of Castile and León, Spain
  • 4 Department of Animal Health, Faculty of Veterinary Medicine, Complutense University of Madrid, Spain

Bovine tuberculosis (bTB) is an infectious chronic disease caused by members of the Mycobacterium tuberculosis complex, mainly M. bovis and M. caprae. It is heterogeneously distributed worldwide, and despite the efforts invested on its eradication, it is still prevalent in several parts of the world (OIE, World Organisation for Animal Health, 2019). In Europe, bTB is still endemic in areas of England, Wales, Ireland and Spain (EFSA, 2018). Bovine tuberculosis causes a high economic impact due to movement restrictions, diagnostic costs and insurances payments. Wildlife species (i.e., wild boar, badger, and red deer) may act as reservoirs of bTB for livestock, adding up additional complexity to the epidemiology of the disease (Nugent et al., 2015; Phillips et al., 2003; de la Rua-Domenech et al., 2006). Eradication programs are based on skin test and passive surveillance at the abattoir. All reactors are subjected to post-mortem analysis to confirm the presence of bTB through bacteriological culture. An extensively assessed feature of bTB is its persistence in certain cattle herds, in terms of either herd recurrence or prolonged periods of restriction (Doyle et al., 2015; White et al., 2013). Recurrence and infection persistence in a herd may indicate failure of testing, thereby allowing infected animals not only to remain in the population, but potentially to act as an ongoing source of infection to other herds and wildlife (Broughan et. al., 2016). Accurate diagnosis of bTB in live animals is often difficult. Several factors have been identified that influence test performance and hence can lead to failures in detecting infection, such as limited sensitivity and specificity, the prevalence of the disease or disease stage in the animal, among others (de la Rua-Domenech et al., 2006). In Spain, challenges to the final eradication of bTB include limitations in the accuracy of currently available diagnostic tests, the presence of wildlife reservoirs (mostly in wild boar and red deer), and undetected bTB in imported cattle (Ciaravino et al., 2018). These factors are associated to persistence of the disease at the farm level and in the environment, respectively. For these reasons, the key to controlling TB in Spain is an integrated approach at all levels in the ecosystem. Herd prevalence in Spain has remained relatively constant during the last 15 years (2.2% in 2002, 2.3% in 2017) (MAPA, 2018). However, the high diversity in the epidemiological situations in different areas of Spain adds even more complexity to the problem. The Castilla y Leon autonomous community, which holds 20% of the Spain’s total cattle population, is classified as a high prevalence region (>1%, 1.43% in 2018), but includes provinces with both high and low bTB herd prevalence (MAPA, 2018). Soria, a low bTB prevalence province (Figure 1) with a median number of 285 herds tested every year for bTB, has however a higher apparent herd prevalence compared with other neighboring provinces (10-15% positive herds detected every year) that are however not confirmed through bacteriology (with only 20-40% of skin positive herds being confirmed every year). These results may be due to a lack of specificity of in-vivo diagnostic tests (i.e., too many false positive animal/herds in the skin test) or a lack of sensitivity of the post-mortem tests (i.e., too many false negative animal/herds in bacteriology). For this reason, we performed a study to describe the spatial patterns in the distribution of positive (confirmed and non-confirmed) bTB cattle/herds in Soria, in an attempt to identify factors that could help establish if there is a lack of sensitivity or specificity in the diagnostic tests used in the bTB national eradication program in Soria. For this reason, we collected information on individual (age, years in the herd) and herd (location, history of bTB, herd size, prevalence at a county level) risk factors from animals and herds with a positive reaction in the skin test (SIT)/IFN-Ƴ assay in Soria during 2015-2018. In addition, we analyzed data of reactors detected in Burgos, a neighboring province in Castilla y Leon, to assess whether patterns identified in Soria were similar in other regions. In each of the regions, positive animals to the in-vivo tests were classified as bTB-confirmed or true positive (SIT or IFN- ƴ+/Culture+) and bTB-unconfirmed or false positive (SIT or IFN-ƴ+/Culture-). Parametric and non-parametric tests were used to detect significant differences between individual/herd variables and bacteriology results. We then explored the spatial variation in the risk of being confirmed through bacteriology for herds in Soria and Burgos: clustering of confirmed/unconfirmed animals in positive herds was assessed using the Poisson model of the spatial scan statistic, implemented using the SaTScan software (Kulldorff, 2009). We then assessed the association between being confirmed through culture (outcome) and the available predictor variables (age, years in the herd, confirmed bTB infection in the previous year to the bacteriology result, herd size and county level prevalence in the previous year) using univariable models with a liberal p-value (p <0.3) for Soria and Burgos. Subsequently, all variables potentially associated were evaluated using a mixed multivariate logistic regression model. To account for the lack of independence between the observations, all models (uni- and multivariable) included herd as a random effect. A backward stepwise selection to obtain the adequate model and comparison of the fit of different models based on likelihood-ratio tests were performed. Residuals of the final model were analyzed for spatial autocorrelation using Moran’s test to determine whether any spatial autocorrelation detected had been explained by the final model (Bivand et al., 2013). Our source population consisted of 30,242 animals subjected to bTB testing and located in 137 and 70 bTB positive farms in Burgos (n=17,845, 59% animals) and Soria (n=12,397, 41% animals), respectively during 2015-2018. Recorded information included herd type (beef, dairy, mixed), herd size, herd location (latitude and longitude), date of testing, bTB diagnostic test results at every round of testing for each animal, and age. All herds in Soria were beef herds, while in Burgos there were 108 (78.8%) beef, 19 (13.9%) dairy, and 10 (7.3%) mixed herds. Post-mortem bacteriology (Culture, C) results were available for 2,124 animals, of which 1,538 (66 C+ animals) and 586 (22 C+ animals) of them were located in farms in Burgos and Soria, respectively. Annual number of animals subjected to bacteriological culture ranged between 165-679 (yearly median C+ animals=15) and 73-238 (yearly median C+ animals=6) in Burgos and Soria, respectively. bTB-confirmed animals located in Burgos were significantly (Mann Whitney test, p-value <0.001) older (median= 7.5 years, IQR [3.9-10.1]) than unconfirmed animals (median= 4.7 years, IQR [2.5-8.3]), while no significant differences were found between bTB-confirmed and unconfirmed animals in Soria (median= 4.3 years for both categories). Poisson spatial models identified the presence of clusters of bTB-confirmed animals in the north-east part of Soria both when analyzing data of the cumulative 2015-2018 time-span or separately for years 2015 and 2016 (Figure 2). In addition, non-confirmed reactors were also clustered in the cumulative 2015-2018 dataset and in each of the years of study in different areas in the north and center part of the province. Same pattern was observed for both bTB-confirmed and unconfirmed herds in Burgos (Figure 3). Among the variables analyzed in the univariable models (Table 1), years in the herd and bTB-confirmation in the herd in the previous year were selected for inclusion in the multivariable analysis in Soria (Table 2). County level prevalence in the previous year was also considered in Burgos (Supplementary material 1, Supplementary material 2). The model with bTB-confirmation in the herd in the previous year and the model with county level prevalence in the previous year had the best fit in the final multivariable models in Soria and Burgos, respectively (Supplementary material 2). Moran’s I test performed on the standardized residuals of the final models for Soria and Burgos revealed no spatial autocorrelation once the effect of the selected variables had been considered in the model. According to the final models, animals in herds in which bTB infection had been confirmed the previous year had a 90% reduction in the odds of being a confirmed in the following year in Soria, while in Burgos this variable was not associated with the outcome. The risk of finding positive culture results in reactors animals increased in animals located in counties in Burgos with a bTB-prevalence higher than 0.5%, while in Soria this variable produced no effect. High (in Soria and Burgos) and low (in Soria) risk bTB-confirmation clusters were found, but when accounting for the effect of significantly associated variables to bTB-confirmation, the odds of being a confirmed animal was similar regardless of the location of the herd. In general, we conclude that the previous history of bTB and prevalence in the area seem to be associated with the success in bTB-confirmation in Soria and Burgos, respectively. Notwithstanding, other variables (i.e., animals in herds located within/outside a cluster of bTB-confirmation, reactors to either SIT or IFN-ƴ) should be considered to further assess whether the risk of being a bTB-confirmed animal is due to the lack of performance of diagnostic tests in Soria.

Figure 1
Figure 2
Figure 3

Acknowledgements

We thank Guillermo Rodriguez Diez and Moises Gonzalez Sanchez for help in data collection.

References

Bivand, R.S., Pebesma, E., & Gomez-Rubio, V. (2013). Applied spatial data analysis with R, Second edition. Springer, NY. http://www.asdar-book.org/ Broughan, J. M., Judge, J., Ely, E., Delahay, R. J., Wilson, G., Clifton-Hadley, R. S., Goodchild, A. V., Bishop, H., Parry, J. E., & Downs, S. H. (2016) A review of risk factors for bovine tuberculosis infection in cattle in the UK and Ireland. Epidemiol. Infect. 144, 2899–2926. doi:10.1017/S095026881600131X Ciaravino, G, Garcia-Saenz, A., Cabras, S., Allepuz, A., Casal, J., Gacia-Bocanegra, I., De Koejier, A., Gubbins, S., Saez, J.L., Cano Terriza, D., & Napp, S. (2018). Epidemics. 23:110-120. doi: 10.1016/j.epidem.2018.01.003 Doyle, L.P., Courcier, E.A., Gordon, A.W., O´Hagan M.J.H., & Menzies, F.D. (2016). Bovine tuberculosis in Northern Ireland: Risk factors associated with duration and recurrence of chronic herd breakdowns. Preventive Veterinary Medicine 131 1–7. http://dx.doi.org/10.1016/j.prevetmed.2016.06.016 EFSA (2016). The European Union summary report on trends and sources of zoonoses, zoonotic agents and food‐borne outbreaks in 2017. EFSA Journal, 16(12):5500, 262 pp. doi: https://doi.org/10.2903/j.efsa.2018.5500 Kulldorff M, Information Management Services I (2009) SaTScan TM v8.0: Software for the spatial and space-time scan statistics. Available: http://www.satscan.org/ Ministerio de Agricultura y Pesca, Alimentación (2018). Programa de Erradicación de Tuberculosis Bovina presentado por España para el año 2019. Retrieved from: https://www.mapa.gob.es/es/ganaderia/temas/sanidad-animal-higiene-ganadera/programatb2019verdefinitiva_tcm30-500265.pdf Nugent, G., Buddle, B. M., & Knowles, G. (2015). Epidemiology and control of Mycobacterium bovis infection in brushtail possums (Trichosurus vulpecula), the primary wildlife host of bovine tuberculosis in New Zealand. New Zealand Veterinary Journal, 63(Suppl 1), 28–41. https://doi.org/10.1080/00480169.2014.963791 OIE, Wold Organisation for Animal Health. (2016). OIE-Listed diseases, infections and infestations in force in 2019. Retrieved from http://www.oie.int/en/animal-health-in-the-world/oie-listed-diseases-2019/ Phillips, C. J., Foster, C. R., Morris, P. A., & Teverson, R. (2003). The transmission of Mycobacterium bovis infection to cattle. Research in Veterinary Science, 74(1), 1–15. https://doi.org/10.1016/S0034-5288(02)00145-5 de la Rua-Domenech, R. (2006). Human Mycobacterium bovis infection in the United Kingdom: Incidence, risks, control measures and review of the zoonotic aspects of bovine tuberculosis. Tuberculosis (Edinb), 86 (2), 77–109. https://doi.org/10.1016/j.tube.2005.05.002 White P. W., Matin S. W., De Jong, M.C.M., O´Keeffe J.J., More, S.J., & Frankena, K. (2013) The importance of ‘neighbourhood’ in the persistence of bovine tuberculosis in Irish cattle herds. Preventive Veterinary Medicine 110 (2013) 346– 355. http://dx.doi.org/10.1016/j.prevetmed.2013.02.012

Keywords: Bovine tuberculosis (bTB), Risk factors, Cluster detection, Spatial Epidemiology, Disease mapping, diagnosis, Cattle, Spain

Conference: GeoVet 2019. Novel spatio-temporal approaches in the era of Big Data, Davis, United States, 8 Oct - 10 Oct, 2019.

Presentation Type: Student oral presentation

Topic: Spatio-temporal surveillance and modeling approaches

Citation: Pozo P, Minguez O, Grau A, Bezos J, Nacar J, Romero B, De Juan L and Alvarez J (2019). To be or not to be (confirmed through bacteriology): spatial epidemiology of bovine tuberculosis in a low prevalence region in Spain. Front. Vet. Sci. Conference Abstract: GeoVet 2019. Novel spatio-temporal approaches in the era of Big Data. doi: 10.3389/conf.fvets.2019.05.00072

Copyright: The abstracts in this collection have not been subject to any Frontiers peer review or checks, and are not endorsed by Frontiers. They are made available through the Frontiers publishing platform as a service to conference organizers and presenters.

The copyright in the individual abstracts is owned by the author of each abstract or his/her employer unless otherwise stated.

Each abstract, as well as the collection of abstracts, are published under a Creative Commons CC-BY 4.0 (attribution) licence (https://creativecommons.org/licenses/by/4.0/) and may thus be reproduced, translated, adapted and be the subject of derivative works provided the authors and Frontiers are attributed.

For Frontiers’ terms and conditions please see https://www.frontiersin.org/legal/terms-and-conditions.

Received: 18 Jun 2019; Published Online: 27 Sep 2019.

* Correspondence: Ms. Pilar Pozo, VISAVET Health Surveillance Centre (UCM), Madrid, Madrid, 28040, Spain, ppozo@ucm.es