Unveiling the Microbial Realm with VEBA 2.0: A modular bioinformatics suite for end-to-end genome-resolved prokaryotic, (micro)eukaryotic, and viral multi-omics from either short- or long-read sequencing.

The microbiome is a complex community of microorganisms, encompassing prokaryotic (bacterial and archaeal), eukaryotic, and viral entities. This microbial ensemble plays a pivotal role in influencing the health and productivity of diverse ecosystems while shaping the web of life. However, many software suites developed to study microbiomes analyze only the prokaryotic community and provide limited to no support for viruses and microeukaryotes. Previously, we introduced the Viral Eukaryotic Bacterial Archaeal (VEBA) open-source software suite to address this critical gap in microbiome research by extending genome-resolved analysis beyond prokaryotes to encompass the understudied realms of eukaryotes and viruses. Here we present VEBA 2.0 with key updates including a comprehensive clustered microeukaryotic protein database, rapid genome/protein-level clustering, bioprospecting, non-coding/organelle gene modeling, genome-resolved taxonomic/pathway profiling, long-read support, and containerization. We demonstrate VEBA's versatile application through the analysis of diverse case studies including marine water, Siberian permafrost, and white-tailed deer lung tissues with the latter showcasing how to identify integrated viruses. VEBA represents a crucial advancement in microbiome research, offering a powerful and accessible platform that bridges the gap between genomics and biotechnological solutions.

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.

Cranes soar on thermal updrafts behind cold fronts as they migrate across the sea.

Thermal soaring conditions above the sea have long been assumed absent or too weak for terrestrial migrating birds, forcing obligate soarers to take long detours and avoid sea-crossing, and facultative soarers to cross exclusively by costly flapping flight. Thus, while atmospheric convection does develop at sea and is used by some seabirds, it has been largely ignored in avian migration research. Here, we provide direct evidence for routine thermal soaring over open sea in the common crane, the heaviest facultative soarer known among terrestrial migrating birds. Using high-resolution biologging from 44 cranes tracked across their transcontinental migration over 4 years, we show that soaring performance was no different over sea than over land in mid-latitudes. Sea-soaring occurred predominantly in autumn when large water-air temperature difference followed mid-latitude cyclones. Our findings challenge a fundamental migration research paradigm and suggest that obligate soarers avoid sea-crossing not due to the absence or weakness of thermals but due to their low frequency, for which they cannot compensate with prolonged flapping. Conversely, facultative soarers other than cranes should also be able to use thermals over the sea. Marine cold air outbreaks, imperative to global energy budget and climate, may also be important for bird migration.

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.

Dynamic winter weather moderates movement and resource selection of wild turkeys at high-latitude range limits.

For wide-ranging species in temperate environments, populations at high-latitude range limits are subject to more extreme conditions, colder temperatures, and greater snow accumulation compared with their core range. As climate change progresses, these bounding pressures may become more moderate on average, while extreme weather occurs more frequently. Individuals can mitigate temporarily extreme conditions by changing daily activity budgets and exhibiting plasticity in resource selection, both of which facilitate existence at and expansion of high-latitude range boundaries. However, relatively little work has explored how animals moderate movement and vary resource selection with changing weather, and a general framework for such investigations is lacking. We applied hidden Markov models and step selection functions to GPS data from wintering wild turkeys (Meleagris gallopavo) near their northern range limit to identify how weather influenced transition among discrete movement states, as well as state-specific resource selection. We found that turkeys were more likely to spend time in a stationary state as wind chill temperatures decreased and snow depth increased. Both stationary and roosting turkeys selected conifer forests and avoided land covers associated with foraging, such as agriculture and residential areas, while shifting their strength of selection for these features during poor weather. In contrast, mobile turkeys showed relatively weak resource selection, with less response in selection coefficients during poor weather. Our findings illustrate that behavioral plasticity in response to weather was context dependent, but movement behaviors most associated with poor weather were also those in which resource selection was most plastic. Given our results, the potential for wild turkey range expansion will partly be determined by the availability of habitat that allows them to withstand periodic inclement weather. Combining hidden Markov models with step selection functions is broadly applicable for evaluating plasticity in animal behavior and dynamic resource selection in response to changing weather. We studied turkeys at northern range limits, but this approach is applicable for any system expected to experience significant changes in the coming decade, and may be particularly relevant to populations existing at range peripheries.

Comparing invasive and noninvasive faecal sampling in wildlife microbiome studies: A case study on wild common cranes.

In ecological and conservation studies, responsible researchers strive to obtain rich data while minimizing disturbance to wildlife and ecosystems. We assessed if samples collected noninvasively can be used for faecal microbiome research, comparing microbiota of noninvasively collected faecal samples to those collected from trapped common cranes at the same sites over the same periods. We found significant differences in faecal microbial composition (alpha and beta diversity), which likely did not result from noninvasive sample exposure to soil contaminants, as assessed by comparing bacterial oxygen use profiles. Differences might result from trapped birds' exposure to sedatives or stress. We conclude that if all samples are collected in the same manner, comparative analyses are valid, and noninvasive sampling may better represent host faecal microbiota because there are no trapping effects. Experiments with fresh and delayed sample collection can elucidate effects of environmental exposures on microbiota. Further, controlled tests of stressing or sedation may unravel how trapping affects wildlife microbiota.

Linking migration and microbiota at a major stopover site in a long-distance avian migrant.

Migration is one of the most physical and energetically demanding periods in an individual bird's life. The composition of the bird's gut or cloacal microbiota can temporarily change during migration, likely due to differences in diets, habitats and other environmental conditions experienced en route. However, how physiological condition, migratory patterns, and other drivers interact to affect microbiota composition of migratory birds is still unclear. We sampled the cloacal bacterial microbiota of a long-distance migrant, the steppe buzzard (Buteo buteo vulpinus), at an important spring stopover bottleneck in Eilat, Israel, after crossing the ca. 1800 km Sahara Desert. We examined whether diversity and composition of the cloacal microbiota varied with body condition, sex, movement patterns (i.e., arrival time and migration distance), and survival. Early arrival to Eilat was associated with better body condition, longer post-Eilat spring migration distance, higher microbial α-diversity, and differences in microbiota composition. Specifically, early arrivals had higher abundance of the phylum Synergistota and five genera, including Jonquetella and Peptococcus, whereas the phylum Proteobacteria and genus Escherichia-Shigella (as well as three other genera) were more abundant in later arrivals. While the differences in α-diversity and Escherichia-Shigella seem to be mainly driven by body condition, other compositional differences associated with arrival date could be indicators of longer migratory journeys (e.g., pre-fueling at wintering grounds or stopover habitats along the way) or migratory performance. No significant differences were found between the microbiota of surviving and non-surviving individuals. Overall, our results indicate that variation in steppe buzzard microbiota is linked to variation in migratory patterns (i.e., capture/arrival date) and body condition, highlighting the importance of sampling the microbiota of GPS-tracked individuals on multiple occasions along their migration routes to gain a more detailed understanding of the links between migration, microbiota, and health in birds.

Including a spatial predictive process in band recovery models improves inference for Lincoln estimates of animal abundance.

Abundance estimation is a critical component of conservation planning, particularly for exploited species where managers set regulations to restrict harvest based on current population size. An increasingly common approach for abundance estimation is through integrated population modeling (IPM), which uses multiple data sources in a joint likelihood to estimate abundance and additional demographic parameters. Lincoln estimators are one commonly used IPM component for harvested species, which combine information on the rate and total number of individuals harvested within an integrated band-recovery framework to estimate abundance at large scales. A major assumption of the Lincoln estimator is that banding and recoveries are representative of the whole population, which may be violated if major sources of spatial heterogeneity in survival or harvest rates are not incorporated into the model. We developed an approach to account for spatial variation in harvest rates using a spatial predictive process, which we incorporated into a Lincoln estimator IPM. We simulated data under different configurations of sample sizes, harvest rates, and sources of spatial heterogeneity in harvest rate to assess potential model bias in parameter estimates. We then applied the model to data collected from a field study of wild turkeys (Meleagris gallapavo) to estimate local and statewide abundance in Maine, USA. We found that the band recovery model that incorporated a spatial predictive process consistently provided estimates of adult and juvenile abundance with low bias across a variety of spatial configurations of harvest rate and sampling intensities. When applied to data collected on wild turkeys, a model that did not incorporate spatial heterogeneity underestimated the harvest rate in some subregions. Consistent with simulation results, this led to overestimation of both local and statewide abundance. Our work demonstrates that a spatial predictive process is a viable mechanism to account for spatial variation in harvest rates and limit bias in abundance estimates. This approach could be extended to large-scale band recovery data sets and has applicability for the estimation of population parameters in other ecological models as well.

First detection of avian influenza subtype H4N6 in Israel in a wild mallard (Anas platyrhynchos).

Avian influenza viruses (AIV) are a worldwide threat to animal and human health. As wild waterfowl circulate and spread these viruses around the world, investigations of AIV prevalence in wild populations are critical for understanding pathogen transmission, as well as predicting disease outbreaks in domestic animals and humans. Surveillance efforts in this study have isolated H4N6 for the first time in Israel from a faecal sample of a wild mallard (Anas platyrhynchos). Phylogenetic analyses of the HA and NA genes revealed that this strain is closely related to isolates from Europe and Asia. This Eurasian origin, together with Israel serving as an important migratory bottleneck of the mid Palearctic-African flyway, suggests a potential introduction of this strain by migratory birds. Additional phylogenetic analysis of the isolate's internal genes (PB1, PB2, PA, NP, M and NS) revealed high levels of phylogenetic relatedness with other AIV subtypes, indicating previous reassortment events. High reassortment rates are characteristic for H4N6 viruses, which, together with this subtype's ability to infect pigs and adaptability to the human receptor binding domain, raises the concern that it would potentially become zoonotic in the future. These results emphasize the importance of continuous AIV monitoring in migratory birds.

PATHOGEN SURVEY AND PREDICTORS OF LYMPHOPROLIFERATIVE DISEASE VIRUS INFECTION IN WILD TURKEYS (MELEAGRIS GALLOPAVO).

Growing populations of Wild Turkeys (Meleagris gallopavo) may result in increased disease transmission among wildlife and spillover to poultry. Lymphoproliferative disease virus (LPDV) is an avian retrovirus that is widespread in Wild Turkeys of eastern North America, and infections may influence mortality and parasite co-infections. We aimed to identify individual and spatial risk factors of LPDV in Maine's Wild Turkeys. We also surveyed for co-infections between LPDV and reticuloendotheliosis virus (REV), Mycoplasma gallisepticum, and Salmonella pullorum to estimate trends in prevalence and examine covariance with LPDV. From 2017 to 2020, we sampled tissues from hunter-harvested (n=72) and live-captured (n=627) Wild Turkeys, in spring and winter, respectively, for molecular detection of LPDV and REV. In a subset of captured individuals (n=235), we estimated seroprevalence of the bacteria M. gallisepticum and S. pullorum using a plate agglutination test. Infection rates for LPDV and REV were 59% and 16% respectively, with a co-infection rate of 10%. Seroprevalence for M. gallisepticum and S. pullorum were 74% and 3.4%, with LPDV co-infection rates of 51% and 2.6%, respectively. Infection with LPDV and seroprevalence of M. gallisepticum and S. pullorum decreased, whereas REV infection increased, between 2018 and 2020. Females (64%), adults (72%), and individuals sampled in spring (76%) had higher risks of LPDV infection than males (47%), juveniles (39%), and individuals sampled in winter (57%). Furthermore, LPDV infection increased with percent forested cover (ß=0.014±0.007) and decreased with percent agriculture cover for juveniles (ß=-0.061±0.018) sampled in winter. These data enhance our understanding of individual and spatial predictors of LPDV infection in Wild Turkeys and aid in assessing the associated risk to Wild Turkey populations and poultry operations.

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.

PREVALENCE AND RISK FACTORS OF ANAPLASMA INFECTIONS IN EASTERN MOOSE (ALCES ALCES AMERICANA) AND WINTER TICKS (DERMACENTOR ALBIPICTUS) IN MAINE, USA.

Eastern moose (Alces alces americana) are heavily parasitized by winter ticks (Dermacentor albipictus), the dominant cause of increased calf mortality in the northeastern US. Although much work has focused on the direct negative effects of winter tick on moose, it remains unknown whether diseases transmitted by ticks may also affect moose health or pose a risk to other species. We explored the role that moose and winter ticks play in transmission of the tick-borne bacterial pathogens, Anaplasma spp., which cause mild to severe illness in humans and domestic animals. Our objectives were to 1) estimate the prevalence of Anaplasma spp. in moose and winter ticks; 2) determine the phylogenetic placement of these strains with respect to those found in other hosts and vectors; and 3) explore risk factors of Anaplasma infection in moose. A total of 157 moose (142 calves, 15 adults) were captured in western (n=83) and northern (n=74) Maine in 2017 and 2018. We screened for Anaplasma spp. in moose whole blood samples using a genus-specific PCR assay targeting the 16S rRNA gene. Over half (54%) of the moose were infected with Anaplasma bacteria, with a greater proportion of moose harboring Anaplasma-infections in the western (67%) versus northern study areas (38%). Male moose exhibited a higher prevalence than did females (63% vs. 47%). In contrast, Anaplasma spp. prevalence in winter ticks was low (<1%). Sequencing and phylogenetic analysis revealed that the single Anaplasma strain in moose was highly divergent from the strain in winter ticks and most closely related to an uncharacterized North American cervid strain. We conclude that winter ticks are unlikely to play a significant role in Anaplasma transmission to moose; however, high infection prevalence warrants further investigation into the impacts of Anaplasma spp. infection on moose health.

Drivers of change and stability in the gut microbiota of an omnivorous avian migrant exposed to artificial food supplementation.

Human activities shape resources available to wild animals, impacting diet and probably 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 analysing the effect of artificially supplemented food provided as part of a crane-agriculture management programme. We sampled faecal droppings in Russia (nonsupplemented, premigration) and in Israel in late autumn (nonsupplemented, postmigration) and winter (supplemented and nonsupplemented, 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 nonsupplemented 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 nonsupplemented groups spent most of their time in agricultural fields, probably feeding on residual grain when available, while food-supplemented cranes spent most of their time at the feeding station. Thus, nonsupplemented and food-supplemented diets probably diverge only in winter, when crop rotation and depletion of anthropogenic resources may lead to a more variable diet in nonsupplemented 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 behaviour towards understanding how human-induced changes shape the gut microbiota in wild animals.

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.

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.