Management boundaries are often delimited by political and social factors, whereas animal movements are affected by ecological and geophysical constraints. Thus, understanding connectivity among distinct management units is of considerable importance, particularly for harvested species, where quotas set in ignorance of connectivity may fail to meet management goals. Our goal was to construct an individual-based model (IBM) to better understand Wild Turkey movements at large scales, benefiting from multiple data sources that are often available for harvested species. We built an IBM describing spring seasonal movements of Wild Turkeys, using data from ringed, radio-, and GPS-marked turkeys captured in Maine, USA. Our IBM accommodated variation in individual turkey response to landscape connectivity metrics and identified emergent migratory connectivity dynamics among harvest management regions. We calculated a low degree of connectivity among wildlife management districts (WMD) which, in combination with the substantial number of boundary crossings observed, indicated a more diffuse distribution of turkeys among WMDs. Nevertheless, the number of turkeys moving between districts provided a clear delineation of where immigration was strongest, identifying which WMDs should be managed as singular population units. This approach has widespread utility for any species or system where harvest management decisions are made at finer spatial scales than the movement dynamics affecting population processes.

Background The distribution of resources can affect animal range sizes, which in turn may alter infectious disease dynamics in heterogenous environments. The risk of pathogen exposure or the spatial extent of outbreaks may vary with host range size. This study examined the range sizes of herbivorous anthrax host species in two ecosystems and relationships between spatial movement behavior and patterns of disease outbreaks for a multi-host environmentally transmitted pathogen. Methods We examined range sizes for seven host species and the spatial extent of anthrax outbreaks in Etosha National Park, Namibia and Kruger National Park, South Africa, where the main host species and outbreak sizes differ. We evaluated host range sizes using the local convex hull method at different temporal scales, within-individual temporal range overlap, and relationships between ranging behavior and species contributions to anthrax cases in each park. We estimated the spatial extent of annual anthrax mortalities and evaluated whether the extent was correlated with case numbers of a given host species. Results Range size differences among species were not linearly related to anthrax case numbers. In Kruger the main host species had small range sizes and high range overlap, which may heighten exposure when outbreaks occur within their ranges. However, different patterns were observed in Etosha, where the main host species had large range sizes and relatively little overlap. The spatial extent of anthrax mortalities was similar between parks but less variable in Etosha than Kruger. In Kruger outbreaks varied from small local clusters to large areas and the spatial extent correlated with case numbers and species affected. Secondary host species contributed relatively few cases to outbreaks; however, for these species with large range sizes, case numbers positively correlated with outbreak extent. Conclusions Our results provide new information on the spatiotemporal structuring of ranging movements of anthrax host species in two ecosystems. The results linking anthrax dynamics to host space use are correlative, yet suggest that, though partial and proximate, host range size and overlap may be contributing factors in outbreak characteristics for environmentally transmitted pathogens.

Environmental factors are common forces driving infectious disease dynamics. We compared inter-annual and seasonal patterns of anthrax infections in two multi-host systems in southern Africa: Etosha National Park, Namibia, and Kruger National Park, South Africa. Using several decades of mortality data from each system, we assessed possible transmission mechanisms behind anthrax dynamics, examining 1) within- and between-species case correlations, and 2) associations between anthrax mortalities and environmental factors, specifically rainfall and the Normalized Difference Vegetation Index (NDVI). Anthrax cases in Kruger had wide inter-annual variation in case numbers, and large outbreaks seemed to follow roughly a decadal cycle. In contrast, outbreaks in Etosha were smaller in magnitude and occurred annually. In Etosha, the host species commonly affected remained consistent over several decades, although plains zebra (Equus quagga) became relatively more dominant. In Kruger, turnover of the main host species occurred after the 1990s, where the previously dominant host species, greater kudu (Tragelaphus strepsiceros), was replaced by impala (Aepyceros melampus). In both parks, anthrax infections showed two seasonal peaks, with each species having only one peak in a year. Zebra, springbok (Antidorcas marsupialis), wildebeest (Connochaetes taurinus) and impala cases peaked in wet seasons, while elephant (Loxodonta africana), kudu and buffalo (Syncerus caffer) cases peaked in dry seasons. For common host species shared between the two parks, anthrax mortalities peaked in the same season in both systems. Among host species with cases peaking in the same season, anthrax mortalities were mostly synchronized, which may imply similar transmission mechanisms or shared sources of exposure. Between seasons, outbreaks in one species may contribute to more cases in another species in the following season. Higher vegetation greenness was associated with more zebra and springbok anthrax mortalities in Etosha, but fewer elephant cases in Kruger. These results suggest that host behavioral responses to changing environmental conditions may affect anthrax transmission risk, with differences in transmission mechanisms leading to multi-host biseasonal outbreaks. This study reveals the dynamics and potential environmental drivers of anthrax in two savanna systems, providing a better understanding of factors driving biseasonal dynamics and outbreak variation among locations.

Human activities shape resources available to wild animals, impacting diet and likely altering their microbiota and overall health. We examined drivers shaping microbiota profiles of common cranes (Grus grus) in agricultural habitats by comparing gut microbiota and crane movement patterns (GPS-tracking) over three periods of their migratory cycle, and by analyzing the effect of artificially-supplemented food provided as part of a crane-agriculture management program. We sampled fecal droppings in Russia (non-supplemented, pre-migration) and in Israel in late fall (non-supplemented, post-migration) and winter (supplemented and non-supplemented, wintering). As supplemented food is typically homogenous, we predicted lower microbiota diversity and different composition in birds relying on supplementary feeding. We did not observe changes in microbial diversity with food supplementation, as diversity differed only in samples from non-supplemented wintering sites. However, both food supplementation and season affected bacterial community composition and led to increased abundance of specific genera (mostly Firmicutes). Cranes from the non-supplemented groups spent most of their time in agricultural fields, likely feeding on residual grain when available, while food-supplemented cranes spent most of their time at the feeding station. Thus, non-supplemented and food-supplemented diets likely diverge only in winter, when crop rotation and depletion of anthropogenic resources may lead to a more variable diet in non-supplemented sites. Our results support the role of diet in structuring bacterial communities and show that they undergo both seasonal and human-induced shifts. Movement analyses provide important clues regarding host diet and behavior towards understanding how human-induced changes shape the gut microbiota in wild animals.

When a transmission hotspot for an environmentally persistent pathogen establishes in otherwise high-quality habitat, the disease may exert a strong impact on a host population. However, fluctuating environmental conditions lead to heterogeneity in habitat quality and animal habitat preference, which may interrupt the overlap between selected and risky habitats. We evaluated spatiotemporal patterns in anthrax mortalities in a plains zebra (Equus quagga) population in Etosha National Park, Namibia, incorporating remote-sensing and host telemetry data. A higher proportion of anthrax mortalities of herbivores was detected in open habitats than in other habitat types. Resource selection functions showed that the zebra population shifted habitat selection in response to changes in rainfall and vegetation productivity. Average to high rainfall years supported larger anthrax outbreaks, with animals congregating in preferred open habitats, while a severe drought forced animals into otherwise less preferred habitats, leading to few anthrax mortalities. Thus, the timing of anthrax outbreaks was congruent with preference for open plains habitats and a corresponding increase in pathogen exposure. Given shifts in habitat preference, the overlap in high-quality habitat and high-risk habitat is intermittent, reducing the adverse consequences for the population.

Bacillus anthracis, the causative agent of anthrax, is a considerable global health threat affecting wildlife, livestock, and the general public. In this study whole-genome sequence analysis of over 350 B. anthracis isolates was used to establish a new high-resolution global genotyping framework that is both biogeographically informative, and compatible with multiple genomic assays. The data presented in this study shed new light on the diverse global dissemination of this species and indicate that many lineages may be uniquely suited to the geographic regions in which they are found. In addition, we demonstrate that plasmid genomic structure for this species is largely consistent with chromosomal population structure, suggesting vertical inheritance in this bacterium has contributed to its evolutionary persistence. This classification methodology is the first based on population genomic structure for this species and has potential use for local and broader institutions seeking to understand both disease outbreak origins and recent introductions. In addition, we provide access to a newly developed genotyping script as well as the full whole genome sequence analyses output for this study, allowing future studies to rapidly employ and append their data in the context of this global collection. This framework may act as a powerful tool for public health agencies, wildlife disease laboratories, and researchers seeking to utilize and expand this classification scheme for further investigations into B. anthracis evolution.