Long-term assessment of relationships between changing environmental conditions and the physiology of southern Beaufort Sea polar bears (Ursus maritimus).

Climate change is influencing polar bear (Ursus maritimus) habitat, diet, and behavior but the effects of these changes on their physiology is not well understood. Blood-based biomarkers are used to assess the physiologic health of individuals but their usefulness for evaluating population health, especially as it relates to changing environmental conditions, has rarely been explored. We describe links between environmental conditions and physiologic functions of southern Beaufort Sea polar bears using data from blood samples collected from 1984 to 2018, a period marked by extensive environmental change. We evaluated associations between 13 physiologic biomarkers and circumpolar (Arctic oscillation index) and regional (wind patterns and ice-free days) environmental metrics and seasonal and demographic co-variates (age, sex, season, and year) known to affect polar bear ecology. We observed signs of dysregulation of water balance in polar bears following years with a lower annual Arctic oscillation index. In addition, liver enzyme values increased over time, which is suggestive of potential hepatocyte damage as the Arctic has warmed. Biomarkers of immune function increased with regional-scale wind patterns and the number of ice-free days over the Beaufort Sea continental shelf and were lower in years with a lower winter Arctic oscillation index, suggesting an increased allocation of energetic resources for immune processes under these conditions. We propose that the variation in polar bear immune and metabolic function is likely indicative of physiologic plasticity, a response that allows polar bears to remain in homeostasis even as they experience changes in nutrition and habitat in response to changing environments.

Estimating abundance of an open population with an N-mixture model using auxiliary data on animal movements.

Accurate assessment of abundance forms a central challenge in population ecology and wildlife management. Many statistical techniques have been developed to estimate population sizes because populations change over time and space and to correct for the bias resulting from animals that are present in a study area but not observed. The mobility of individuals makes it difficult to design sampling procedures that account for movement into and out of areas with fixed jurisdictional boundaries. Aerial surveys are the gold standard used to obtain data of large mobile species in geographic regions with harsh terrain, but these surveys can be prohibitively expensive and dangerous. Estimating abundance with ground-based census methods have practical advantages, but it can be difficult to simultaneously account for temporary emigration and observer error to avoid biased results. Contemporary research in population ecology increasingly relies on telemetry observations of the states and locations of individuals to gain insight on vital rates, animal movements, and population abundance. Analytical models that use observations of movements to improve estimates of abundance have not been developed. Here we build upon existing multi-state mark-recapture methods using a hierarchical N-mixture model with multiple sources of data, including telemetry data on locations of individuals, to improve estimates of population sizes. We used a state-space approach to model animal movements to approximate the number of marked animals present within the study area at any observation period, thereby accounting for a frequently changing number of marked individuals. We illustrate the approach using data on a population of elk (Cervus elaphus nelsoni) in Northern Colorado, USA. We demonstrate substantial improvement compared to existing abundance estimation methods and corroborate our results from the ground based surveys with estimates from aerial surveys during the same seasons. We develop a hierarchical Bayesian N-mixture model using multiple sources of data on abundance, movement and survival to estimate the population size of a mobile species that uses remote conservation areas. The model improves accuracy of inference relative to previous methods for estimating abundance of open populations.

Comparison of WAIC and posterior predictive approaches for N-mixture models.

Hierarchical models are common for ecological analysis, but determining appropriate model selection methods remains an ongoing challenge. To confront this challenge, a suitable method is needed to evaluate and compare available candidate models. We compared performance of conditional WAIC, a joint-likelihood approach to WAIC (WAICj), and posterior-predictive loss for selecting between candidate N-mixture models. We tested these model selection criteria on simulated single-season N-mixture models, simulated multi-season N-mixture models with temporal auto-correlation, and three case studies of single-season N-mixture models based on eBird data. WAICj proved more accurate than the standard conditional formulation or posterior-predictive loss, even when models were temporally correlated, suggesting WAICj is a robust alternative to model selection for N-mixture models.

Serum Virome of Southern Beaufort Sea polar bears (Ursus maritimus) during a period of rapid climate change.

Climate change affects the behavior, physiology and life history of many Arctic wildlife species. It can also influence the distribution and ecology of infectious agents. The southern Beaufort Sea (SB) subpopulation of polar bears (Ursus maritimus) has experienced dramatic behavioral changes due to retreating sea ice and other climate-related factors, but the effects of these changes on physiology and infection remain poorly understood. Using serum from polar bears sampled between 2004 and 2015 and metagenomic DNA sequencing, we identified 48 viruses, all of the family Anelloviridae. Anelloviruses are small, ubiquitous infectious agents with circular single-stranded DNA genomes that are not known to cause disease but, in humans, covary in diversity and load with immunological compromise. We therefore examined the usefulness of anelloviruses as biomarkers of polar bear physiological stress related to climate and habitat use. Polar bear anelloviruses sorted into two distinct clades on a phylogenetic tree, both of which also contained anelloviruses of giant pandas (Ailuropoda melanoleuca), another ursid. Neither anellovirus diversity nor load were associated with any demographic variables, behavioral factors or direct physiological measures. However, pairwise genetic distances between anelloviruses were positively correlated with pairwise differences in sampling date, suggesting that the polar bear "anellome" is evolving over time. These findings suggest that anelloviruses are not a sensitive indicator of polar physiological stress, but they do provide a baseline for evaluating future changes to polar bear viromes.

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.

A hierarchical Bayesian approach for handling missing classification data.

Ecologists use classifications of individuals in categories to understand composition of populations and communities. These categories might be defined by demographics, functional traits, or species. Assignment of categories is often imperfect, but frequently treated as observations without error. When individuals are observed but not classified, these "partial" observations must be modified to include the missing data mechanism to avoid spurious inference.We developed two hierarchical Bayesian models to overcome the assumption of perfect assignment to mutually exclusive categories in the multinomial distribution of categorical counts, when classifications are missing. These models incorporate auxiliary information to adjust the posterior distributions of the proportions of membership in categories. In one model, we use an empirical Bayes approach, where a subset of data from one year serves as a prior for the missing data the next. In the other approach, we use a small random sample of data within a year to inform the distribution of the missing data.We performed a simulation to show the bias that occurs when partial observations were ignored and demonstrated the altered inference for the estimation of demographic ratios. We applied our models to demographic classifications of elk (Cervus elaphus nelsoni) to demonstrate improved inference for the proportions of sex and stage classes.We developed multiple modeling approaches using a generalizable nested multinomial structure to account for partially observed data that were missing not at random for classification counts. Accounting for classification uncertainty is important to accurately understand the composition of populations and communities in ecological studies.