A Bayesian spatio-temporal model for forecasting the prevalence of antibodies to Ehrlichia species in domestic dogs within the contiguous United States

Yan Liu, Robert B. Lund, Shila K. Nordone, Michael J. Yabsley, Christopher S. McMahan

Dogs in the United States are hosts to a diverse range of vector-borne pathogens, several of which are important zoonoses. This paper describes factors deemed to be significantly related to the prevalence of antibodies to Ehrlichia spp. in domestic dogs, including climatic conditions, geographical factors, and societal factors. These factors are used in concert with a spatio-temporal model to construct an annual seroprevalence forecast. The proposed method of forecasting and an assessment of its fidelity are described. Approximately twelve million serological test results for canine exposure to Ehrlichia spp. were used in the development of a Bayesian approach to forecast canine infection. Data used were collected on the county level across the contiguous United States from routine veterinary diagnostic tests between 2011-2015. Maps depicting the spatial baseline Ehrlichia spp. prevalence were constructed using Kriging and head-banging smoothing methods. Data were statistically analyzed to identify factors related to antibody prevalence via a Bayesian spatio-temporal conditional autoregressive (CAR) model. Finally, a forecast of future Ehrlichia seroprevalence was constructed based on the proposed model using county-level data on five predictive factors identified at a workshop hosted by the Companion Animal Parasite Council and published in 2014: annual temperature, percentage forest coverage, percentage surface water coverage, population density and median household income. Data were statistically analyzed to identify factors related to disease prevalence via a Bayesian spatio-temporal model. The fitted model and factor extrapolations were then used to forecast the regional seroprevalence for 2016. The correlation between the observed and model-estimated county-by-county Ehrlichia seroprevalence for the five-year period 2011-2015 is 0.842, demonstrating reasonable model accuracy. The weighted correlation (acknowledging unequal sample sizes) between 2015 observed and forecasted county-by-county Ehrlichia seroprevalence is 0.970, demonstrating that Ehrlichia seroprevalence can be forecasted accurately. The forecast presented herein can be an a priori alert to veterinarians regarding areas expected to see expansion of Ehrlichia beyond the accepted endemic range, or in some regions a dynamic change from historical average prevalence. Moreover, this forecast could potentially serve as a surveillance tool for human health and prove useful for forecasting other vector-borne diseases.