Metagenomic sequencing sheds light on microbes putatively associated with pneumonia-related fatalities of white-tailed deer (Odocoileus virginianus).

With emerging infectious disease outbreaks in human, domestic and wild animal populations on the rise, improvements in pathogen characterization and surveillance are paramount for the protection of human and animal health, as well as the conservation of ecologically and economically important wildlife. Genomics offers a range of suitable tools to meet these goals, with metagenomic sequencing facilitating the characterization of whole microbial communities associated with emerging and endemic disease outbreaks. Here, we use metagenomic sequencing in a case-control study to identify microbes in lung tissue associated with newly observed pneumonia-related fatalities in 34 white-tailed deer (Odocoileus virginianus) in Wisconsin, USA. We identified 20 bacterial species that occurred in more than a single individual. Of these, only Clostridium novyi was found to substantially differ (in number of detections) between case and control sample groups; however, this difference was not statistically significant. We also detected several bacterial species associated with pneumonia and/or other diseases in ruminants (Mycoplasma ovipneumoniae, Trueperella pyogenes, Pasteurella multocida, Anaplasma phagocytophilum, Fusobacterium necrophorum); however, these species did not substantially differ between case and control sample groups. On average, we detected a larger number of bacterial species in case samples than controls, supporting the potential role of polymicrobial infections in this system. Importantly, we did not detect DNA of viruses or fungi, suggesting that they are not significantly associated with pneumonia in this system. Together, these results highlight the utility of metagenomic sequencing for identifying disease-associated microbes. This preliminary list of microbes will help inform future research on pneumonia-associated fatalities of white-tailed deer.

Season of death, pathogen persistence and wildlife behaviour alter number of anthrax secondary infections from environmental reservoirs.

An important part of infectious disease management is predicting factors that influence disease outbreaks, such as R, the number of secondary infections arising from an infected individual. Estimating R is particularly challenging for environmentally transmitted pathogens given time lags between cases and subsequent infections. Here, we calculated R for Bacillus anthracis infections arising from anthrax carcass sites in Etosha National Park, Namibia. Combining host behavioural data, pathogen concentrations and simulation models, we show that R is spatially and temporally variable, driven by spore concentrations at death, host visitation rates and early preference for foraging at infectious sites. While spores were detected up to a decade after death, most secondary infections occurred within 2 years. Transmission simulations under scenarios combining site infectiousness and host exposure risk under different environmental conditions led to dramatically different outbreak dynamics, from pathogen extinction (R < 1) to explosive outbreaks (R > 10). These transmission heterogeneities may explain variation in anthrax outbreak dynamics observed globally, and more generally, the critical importance of environmental variation underlying host-pathogen interactions. Notably, our approach allowed us to estimate the lethal dose of a highly virulent pathogen non-invasively from observational studies and epidemiological data, useful when experiments on wildlife are undesirable or impractical.

Development of conjugated secondary antibodies for wildlife disease surveillance.

Disease monitoring in free-ranging wildlife is a challenge and often relies on passive surveillance. Alternatively, proactive surveillance that relies on the detection of specific antibodies could give more reliable and timely insight into disease presence and prevalence in a population, especially if the evidence of disease occurs below detection thresholds for passive surveillance. Primary binding assays, like the indirect ELISA for antibody detection in wildlife, are hampered by a lack of species-specific conjugates. In this study, we developed anti-kudu (Tragelaphus strepsiceros) and anti-impala (Aepyceros melampus) immunoglobulin-specific conjugates in chickens and compared them to the binding of commercially available protein-G and protein-AG conjugates, using an ELISA-based avidity index. The conjugates were evaluated for cross-reaction with sera from other wild herbivores to assess future use in ELISAs. The developed conjugates had a high avidity of >70% against kudu and impala sera. The commercial conjugates (protein-G and protein-AG) had significantly low relative avidity (<20%) against these species. Eighteen other wildlife species demonstrated cross-reactivity with a mean relative avidity of >50% with the impala and kudu conjugates and <40% with the commercial conjugates. These results demonstrate that species-specific conjugates are important tools for the development and validation of immunoassays in wildlife and for the surveillance of zoonotic agents along the livestock-wildlife-human interface.

Variation in herbivore space use: comparing two savanna ecosystems with different anthrax outbreak patterns in southern Africa.

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.

The persistence of time: the lifespan of Bacillus anthracis spores in environmental reservoirs.

Anthrax is a lethal bacterial zoonosis primarily affecting herbivorous wildlife and livestock. Upon host death Bacillus anthracis vegetative cells form spores capable of surviving for years in soil. Anthrax transmission requires host exposure to large spore doses. Thus, conditions that facilitate higher spore concentrations or promote spore survival will increase the probability that a pathogen reservoir infects future hosts. We investigated abiotic and pathogen genomic variation in relation to spore concentrations in surface soils (0-1 cm depth) at 40 plains zebra (Equus quagga) anthrax carcass sites in Namibia. Specifically, how initial spore concentrations and spore survival were affected by seasonality associated with the timing of host mortality, local soil characteristics, and pathogen genomic variation. Zebras dying of anthrax in wet seasons-the peak season for anthrax in Etosha National Park-had soil spore concentrations 1.36 orders of magnitude higher than those that died in dry seasons. No other variables considered affected spore concentrations, and spore survival rates did not differ among sites. Surface soils at these pathogen reservoirs remained culture positive for a range of 3.8-10.4 years after host death. Future research could evaluate if seasonal patterns in spore concentrations are driven by differences in sporulation success or levels of terminal bacteremia.

Agricultural land use shapes dispersal in white-tailed deer (Odocoileus virginianus).

BACKGROUND: Dispersal is a fundamental process to animal population dynamics and gene flow. In white-tailed deer (WTD; Odocoileus virginianus), dispersal also presents an increasingly relevant risk for the spread of infectious diseases. Across their wide range, WTD dispersal is believed to be driven by a suite of landscape and host behavioral factors, but these can vary by region, season, and sex. Our objectives were to (1) identify dispersal events in Wisconsin WTD and determine drivers of dispersal rates and distances, and (2) determine how landscape features (e.g., rivers, roads) structure deer dispersal paths. METHODS: We developed an algorithmic approach to detect dispersal events from GPS collar data for 590 juvenile, yearling, and adult WTD. We used statistical models to identify host and landscape drivers of dispersal rates and distances, including the role of agricultural land use, the traversability of the landscape, and potential interactions between deer. We then performed a step selection analysis to determine how landscape features such as agricultural land use, elevation, rivers, and roads affected deer dispersal paths. RESULTS: Dispersal predominantly occurred in juvenile males, of which 64.2% dispersed, with dispersal events uncommon in other sex and age classes. Juvenile male dispersal probability was positively associated with the proportion of the natal range that was classified as agricultural land use, but only during the spring. Dispersal distances were typically short (median 5.77 km, range: 1.3-68.3 km), especially in the fall. Further, dispersal distances were positively associated with agricultural land use in potential dispersal paths but negatively associated with the number of proximate deer in the natal range. Lastly, we found that, during dispersal, juvenile males typically avoided agricultural land use but selected for areas near rivers and streams. CONCLUSION: Land use-particularly agricultural-was a key driver of dispersal rates, distances, and paths in Wisconsin WTD. In addition, our results support the importance of deer social environments in shaping dispersal behavior. Our findings reinforce knowledge of dispersal ecology in WTD and how landscape factors-including major rivers, roads, and land-use patterns-structure host gene flow and potential pathogen transmission.

CAUSE OF DEATH, PATHOLOGY, AND CHRONIC WASTING DISEASE STATUS OF WHITE-TAILED DEER (ODOCOILEUS VIRGINIANUS) MORTALITIES IN WISCONSIN, USA.

White-tailed deer (WTD; Odocoileus virginianus) are a critical species for ecosystem function and wildlife management. As such, studies of cause-specific mortality among WTD have long been used to understand population dynamics. However, detailed pathological information is rarely documented for free-ranging WTD, especially in regions with a high prevalence of chronic wasting disease (CWD). This leaves a significant gap in understanding how CWD is associated with disease processes or comorbidities that may subsequently alter broader population dynamics. We investigated unknown mortalities among collared WTD in southwestern Wisconsin, USA, an area of high CWD prevalence. We tested for associations between CWD and other disease processes and used a network approach to test for co-occurring disease processes. Predation and infectious disease were leading suspected causes of death, with high prevalence of CWD (42.4%; of 245 evaluated) and pneumonia (51.2%; of 168 evaluated) in our sample. CWD prevalence increased with age, before decreasing among older individuals, with more older females than males in our sample. Females were more likely to be CWD positive, and although this was not statistically significant when accounting for age, females were significantly more likely to die with end-stage CWD than males and may consequently be an underrecognized source of CWD transmission. Presence of CWD was associated with emaciation, atrophy of marrow fat and hematopoietic cells, and ectoparasitism (lice and ticks). Occurrences of severe infectious disease processes clustered together (e.g., pneumonia, CWD), as compared to noninfectious or low-severity processes (e.g., sarcocystosis), although pneumonia cases were not fully explained by CWD status. With the prevalence of CWD increasing across North America, our results highlight the critical importance of understanding the potential role of CWD in favoring or maintaining disease processes of importance for deer population health and dynamics.

A framework for integrating inferred movement behavior into disease risk models.

Movement behavior is an important contributor to habitat selection and its incorporation in disease risk models has been somewhat neglected. The habitat preferences of host individuals affect their probability of exposure to pathogens. If preference behavior can be incorporated in ecological niche models (ENMs) when data on pathogen distributions are available, then variation in such behavior may dramatically impact exposure risk. Here we use data from the anthrax endemic system of Etosha National Park, Namibia, to demonstrate how integrating inferred movement behavior alters the construction of disease risk maps. We used a Maximum Entropy (MaxEnt) model that associated soil, bioclimatic, and vegetation variables with the best available pathogen presence data collected at anthrax carcass sites to map areas of most likely Bacillus anthracis (the causative bacterium of anthrax) persistence. We then used a hidden Markov model (HMM) to distinguish foraging and non-foraging behavioral states along the movement tracks of nine zebra (Equus quagga) during the 2009 and 2010 anthrax seasons. The resulting tracks, decomposed on the basis of the inferred behavioral state, formed the basis of step-selection functions (SSFs) that used the MaxEnt output as a potential predictor variable. Our analyses revealed different risks of exposure during different zebra behavioral states, which were obscured when the full movement tracks were analyzed without consideration of the underlying behavioral states of individuals. Pathogen (or vector) distribution models may be misleading with regard to the actual risk faced by host animal populations when specific behavioral states are not explicitly accounted for in selection analyses. To more accurately evaluate exposure risk, especially in the case of environmentally transmitted pathogens, selection functions could be built for each identified behavioral state and then used to assess the comparative exposure risk across relevant states. The scale of data collection and analysis, however, introduces complexities and limitations for consideration when interpreting results.

Immunological Evidence of Variation in Exposure and Immune Response to Bacillus anthracis in Herbivores of Kruger and Etosha National Parks.

Exposure and immunity to generalist pathogens differ among host species and vary across spatial scales. Anthrax, caused by a multi-host bacterial pathogen, Bacillus anthracis, is enzootic in Kruger National Park (KNP), South Africa and Etosha National Park (ENP), Namibia. These parks share many of the same potential host species, yet the main anthrax host in one (greater kudu (Tragelaphus strepsiceros) in KNP and plains zebra (Equus quagga) in ENP) is only a minor host in the other. We investigated species and spatial patterns in anthrax mortalities, B. anthracis exposure, and the ability to neutralize the anthrax lethal toxin to determine if observed host mortality differences between locations could be attributed to population-level variation in pathogen exposure and/or immune response. Using serum collected from zebra and kudu in high and low incidence areas of each park (18- 20 samples/species/area), we estimated pathogen exposure from anti-protective antigen (PA) antibody response using enzyme-linked immunosorbent assay (ELISA) and lethal toxin neutralization with a toxin neutralization assay (TNA). Serological evidence of pathogen exposure followed mortality patterns within each system (kudus: 95% positive in KNP versus 40% in ENP; zebras: 83% positive in ENP versus 63% in KNP). Animals in the high-incidence area of KNP had higher anti-PA responses than those in the low-incidence area, but there were no significant differences in exposure by area within ENP. Toxin neutralizing ability was higher for host populations with lower exposure prevalence, i.e., higher in ENP kudus and KNP zebras than their conspecifics in the other park. These results indicate that host species differ in their exposure to and adaptive immunity against B. anthracis in the two parks. These patterns may be due to environmental differences such as vegetation, rainfall patterns, landscape or forage availability between these systems and their interplay with host behavior (foraging or other risky behaviors), resulting in differences in exposure frequency and dose, and hence immune response.

Disease or drought: environmental fluctuations release zebra from a potential pathogen-triggered ecological trap.

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 spatio-temporal 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.

The roles of environmental variation and parasite survival in virulence-transmission relationships.

Disease outbreaks are a consequence of interactions among the three components of a host-parasite system: the infectious agent, the host and the environment. While virulence and transmission are widely investigated, most studies of parasite life-history trade-offs are conducted with theoretical models or tractable experimental systems where transmission is standardized and the environment controlled. Yet, biotic and abiotic environmental factors can strongly affect disease dynamics, and ultimately, host-parasite coevolution. Here, we review research on how environmental context alters virulence-transmission relationships, focusing on the off-host portion of the parasite life cycle, and how variation in parasite survival affects the evolution of virulence and transmission. We review three inter-related 'approaches' that have dominated the study of the evolution of virulence and transmission for different host-parasite systems: (i) evolutionary trade-off theory, (ii) parasite local adaptation and (iii) parasite phylodynamics. These approaches consider the role of the environment in virulence and transmission evolution from different angles, which entail different advantages and potential biases. We suggest improvements to how to investigate virulence-transmission relationships, through conceptual and methodological developments and taking environmental context into consideration. By combining developments in life-history evolution, phylogenetics, adaptive dynamics and comparative genomics, we can improve our understanding of virulence-transmission relationships across a diversity of host-parasite systems that have eluded experimental study of parasite life history.

The Epidemiological Signature of Pathogen Populations That Vary in the Relationship between Free-Living Parasite Survival and Virulence.

The relationship between parasite virulence and transmission is a pillar of evolutionary theory that has implications for public health. Part of this canon involves the idea that virulence and free-living survival (a key component of transmission) may have different relationships in different host-parasite systems. Most examinations of the evolution of virulence-transmission relationships-Theoretical or empirical in nature-Tend to focus on the evolution of virulence, with transmission being a secondary consideration. Even within transmission studies, the focus on free-living survival is a smaller subset, though recent studies have examined its importance in the ecology of infectious diseases. Few studies have examined the epidemic-scale consequences of variation in survival across different virulence-survival relationships. In this study, we utilize a mathematical model motivated by aspects of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) natural history to investigate how evolutionary changes in survival may influence several aspects of disease dynamics at the epidemiological scale. Across virulence-survival relationships (where these traits are either positively or negatively correlated), we found that small changes (5% above and below the nominal value) in survival can have a meaningful effect on certain outbreak features, including R0, and on the size of the infectious peak in the population. These results highlight the importance of properly understanding the mechanistic relationship between virulence and parasite survival, as the evolution of increased survival across different relationships with virulence may have considerably different epidemiological signatures.

Drivers of community turnover differ between avian hemoparasite genera along a North American latitudinal gradient.

The latitudinal diversity gradient (LDG) is an established macroecological pattern, but is poorly studied in microbial organisms, particularly parasites. In this study, we tested whether latitude, elevation, and host species predicted patterns of prevalence, alpha diversity, and community turnover of hemosporidian parasites. We expected parasite diversity to decrease with latitude, alongside the diversity of their hosts and vectors. Similarly, we expected infection prevalence to decrease with latitude as vector abundances decrease. Lastly, we expected parasite community turnover to increase with latitudinal distance and to be higher between rather than within host species. We tested these hypotheses by screening blood and tissue samples of three closely related avian species in a clade of North American songbirds (Turdidae: Catharus, n = 466) across 17.5° of latitude. We used a nested PCR approach to identify parasites in hemosporidian genera that are transmitted by different dipteran vectors. Then, we implemented linear-mixed effects and generalized dissimilarity models to evaluate the effects of latitude, elevation, and host species on parasite metrics. We found high diversity of hemosporidian parasites in Catharus thrushes (n = 44 lineages) but no evidence of latitudinal gradients in alpha diversity or prevalence. Parasites in the genus Leucocytozoon were most prevalent and lineage rich in this study system; however, there was limited turnover with latitude and host species. Contrastingly, Plasmodium parasites were less prevalent and diverse than Leucocytozoon parasites, yet communities turned over at a higher rate with latitude and host species. Leucocytozoon communities were skewed by the dominance of one or two highly prevalent lineages with broad latitudinal distributions. The few studies that evaluate the hemosporidian LDG do not find consistent patterns of prevalence and diversity, which makes it challenging to predict how they will respond to global climate change.

A classification framework for Bacillus anthracis defined by global genomic structure.

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.

Coalescence modeling of intrainfection Bacillus anthracis populations allows estimation of infection parameters in wild populations.

Bacillus anthracis, the etiological agent of anthrax, is a well-established model organism. For B. anthracis and most other infectious diseases, knowledge regarding transmission and infection parameters in natural systems, in large part, comprises data gathered from closely controlled laboratory experiments. Fatal, natural anthrax infections transmit the bacterium through new host-pathogen contacts at carcass sites, which can occur years after death of the previous host. For the period between contact and death, all of our knowledge is based upon experimental data from domestic livestock and laboratory animals. Here we use a noninvasive method to explore the dynamics of anthrax infections, by evaluating the terminal diversity of B. anthracis in anthrax carcasses. We present an application of population genetics theory, specifically, coalescence modeling, to intrainfection populations of B. anthracis to derive estimates for the duration of the acute phase of the infection and effective population size converted to the number of colony-forming units establishing infection in wild plains zebra (Equus quagga). Founding populations are small, a few colony-forming units, and infections are rapid, lasting roughly between 1 d and 3 d in the wild. Our results closely reflect experimental data, showing that small founding populations progress acutely, killing the host within days. We believe this method is amendable to other bacterial diseases from wild, domestic, and human systems.

The Expectations and Challenges of Wildlife Disease Research in the Era of Genomics: Forecasting with a Horizon Scan-like Exercise.

The outbreak and transmission of disease-causing pathogens are contributing to the unprecedented rate of biodiversity decline. Recent advances in genomics have coalesced into powerful tools to monitor, detect, and reconstruct the role of pathogens impacting wildlife populations. Wildlife researchers are thus uniquely positioned to merge ecological and evolutionary studies with genomic technologies to exploit unprecedented "Big Data" tools in disease research; however, many researchers lack the training and expertise required to use these computationally intensive methodologies. To address this disparity, the inaugural "Genomics of Disease in Wildlife" workshop assembled early to mid-career professionals with expertise across scientific disciplines (e.g., genomics, wildlife biology, veterinary sciences, and conservation management) for training in the application of genomic tools to wildlife disease research. A horizon scanning-like exercise, an activity to identify forthcoming trends and challenges, performed by the workshop participants identified and discussed 5 themes considered to be the most pressing to the application of genomics in wildlife disease research: 1) "Improving communication," 2) "Methodological and analytical advancements," 3) "Translation into practice," 4) "Integrating landscape ecology and genomics," and 5) "Emerging new questions." Wide-ranging solutions from the horizon scan were international in scope, itemized both deficiencies and strengths in wildlife genomic initiatives, promoted the use of genomic technologies to unite wildlife and human disease research, and advocated best practices for optimal use of genomic tools in wildlife disease projects. The results offer a glimpse of the potential revolution in human and wildlife disease research possible through multi-disciplinary collaborations at local, regional, and global scales.

Modeling R0 for Pathogens with Environmental Transmission: Animal Movements, Pathogen Populations, and Local Infectious Zones.

How a disease is transmitted affects our ability to determine R0, the average number of new cases caused by an infectious host at the onset of an epidemic. R0 becomes progressively more difficult to compute as transmission varies from directly transmitted diseases to diseases that are vector-borne to environmentally transmitted diseases. Pathogens responsible for diseases with environmental transmission are typically maintained in environmental reservoirs that exhibit a complex spatial distribution of local infectious zones (LIZs). Understanding host encounters with LIZs and pathogen persistence within LIZs is required for an accurate R0 and modeling these contacts requires an integrated geospatial and dynamical systems approach. Here we review how interactions between host and pathogen populations and environmental reservoirs are driven by landscape-level variables, and synthesize the quantitative framework needed to formulate outbreak response and disease control.

Spores and soil from six sides: interdisciplinarity and the environmental biology of anthrax (Bacillus anthracis).

Environmentally transmitted diseases are comparatively poorly understood and managed, and their ecology is particularly understudied. Here we identify challenges of studying environmental transmission and persistence with a six-sided interdisciplinary review of the biology of anthrax (Bacillus anthracis). Anthrax is a zoonotic disease capable of maintaining infectious spore banks in soil for decades (or even potentially centuries), and the mechanisms of its environmental persistence have been the topic of significant research and controversy. Where anthrax is endemic, it plays an important ecological role, shaping the dynamics of entire herbivore communities. The complex eco-epidemiology of anthrax, and the mysterious biology of Bacillus anthracis during its environmental stage, have necessitated an interdisciplinary approach to pathogen research. Here, we illustrate different disciplinary perspectives through key advances made by researchers working in Etosha National Park, a long-term ecological research site in Namibia that has exemplified the complexities of the enzootic process of anthrax over decades of surveillance. In Etosha, the role of scavengers and alternative routes (waterborne transmission and flies) has proved unimportant relative to the long-term persistence of anthrax spores in soil and their infection of herbivore hosts. Carcass deposition facilitates green-ups of vegetation to attract herbivores, potentially facilitated by the role of anthrax spores in the rhizosphere. The underlying seasonal pattern of vegetation, and herbivores' immune and behavioural responses to anthrax risk, interact to produce regular 'anthrax seasons' that appear to be a stable feature of the Etosha ecosystem. Through the lens of microbiologists, geneticists, immunologists, ecologists, epidemiologists, and clinicians, we discuss how anthrax dynamics are shaped at the smallest scale by population genetics and interactions within the bacterial communities up to the broadest scales of ecosystem structure. We illustrate the benefits and challenges of this interdisciplinary approach to disease ecology, and suggest ways anthrax might offer insights into the biology of other important pathogens. Bacillus anthracis, and the more recently emerged Bacillus cereus biovar anthracis, share key features with other environmentally transmitted pathogens, including several zoonoses and panzootics of special interest for global health and conservation efforts. Understanding the dynamics of anthrax, and developing interdisciplinary research programs that explore environmental persistence, is a critical step forward for understanding these emerging threats.