Space matters: Host spatial structure and the dynamics of plague transmission

The development of models to elucidate the transmission pathways and dynamics of wildlife diseases remains challenging. Sylvatic plague, caused by the bacterium Yersinia pestis (Yp), is an infectious zoonotic disease that primarily affects wild rodents, including prairie dogs (Cynomys spp.) in North America. Proposed transmission pathways for Yp include flea bites, direct contacts between hosts, and environmental reservoirs (e.g. soil, carcasses). We developed a spatially explicit, agent-based model of Yp transmission to explore the effects of alternative transmission pathways, different disease initiation mechanisms (host or fleas), parameter uncertainty, and spatial structure of hosts. A particularly novel aspect of our model was the integration of ecological models with traditional disease models. Specifically, we used estimates from spatial capture-recapture models to generate data-driven spatial distributions, densities, and contact rates to capture the spatial structure of prairie dogs. We simulated ~9 million scenarios across a wide range of parameter values and conducted sensitivity analyses to determine the most influential parameters on the number of flea-days (sum of the mean number of fleas on hosts each day of the simulation), number of newly infected hosts per day, the time to depopulation (<20 prairie dogs remaining), and the proportion of the prairie dog population remaining at the end of the simulation (after 150 days). When including spatial structure, we found the probability of transmission via environmental sources of Yp (i.e. carcasses) had the greatest influence on model results when Yp infection was initiated in prairie dog hosts, rather than in fleas. Conversely, the mechanism of transmission by fleas to prairie dogs had the greatest influence on model results when Yp infection was initiated in fleas (i.e. via introduction by carnivores, a migrant prairie dog, or other mammalian host). Uncertainty in parameter estimates, particularly those related to the transmission pathways of Yp, continue to hamper efforts to realistically model plague dynamics in wild rodents. Our results elucidate the complexity of the flea-plague-prairie dog system and reiterate the importance of research on Yp transmission mechanisms to provide a full understanding of this disease. Our results also emphasize the importance of realistic estimates of spatial structure for exploring transmission dynamics of wildlife diseases and provide a framework for generating a data-driven description of spatial structure.