Detection and Phylogenetic Analysis of Orf Virus and Muskox Rhadinovirus 1 from Muskoxen (Ovibos moschatus) in the Canadian Arctic

Orf virus (genus Parapoxvirus) has been associated with gross skin lesions on muskoxen (Ovibos moschatus) from Victoria Island, Nunavut, Canada, where muskox populations are experiencing population declines. Orf virus causes painful proliferative and necrotizing dermatitis upon viral replication and shedding, which may lead to animal morbidity or mortality through secondary infections and starvation. Herpesvirus, known to cause gross lesions on skin and mucosa during active viral replication, has also been documented in muskoxen but to date has not been associated with clinical disease. Our objective was to characterize the variation of orf virus and herpesvirus in wild muskoxen of the Canadian Arctic. Tissue samples including gross skin lesions from the nose, lips, and/or legs were opportunistically collected from muskoxen on Victoria Island, Nunavut and Northwest Territories, and mainland Nunavut, Canada, from 2015 to 2017. Sampled muskoxen varied in age, sex, location, hunt type, and body condition. Tissues from 60 muskoxen were tested for genetic evidence of orf virus and herpesvirus infection using PCR targeting key viral genes. Tissues from 38 muskoxen, including 15 with gross lesions, were also examined for histological evidence of orf virus and herpesvirus infection. Eleven muskoxen (10 from Victoria Island and one from mainland Nunavut) with gross lesions had microscopic lesions consistent with orf virus infection. Muskox rhadinovirus 1, a gammaherpesvirus endemic to muskoxen, was detected in 33 (55%) muskoxen including 17 with gross lesions. In all tissues examined, there was no histological evidence of herpesvirus-specific disease. Sequencing and characterization of amplified PCR products using phylogenetic analysis indicated that a strain of orf virus, which appears to be unique, is likely to be endemic in muskoxen from Victoria Island and mainland Nunavut. Many of the muskoxen are also subclinically infected with a known muskox-endemic strain of herpesvirus.

Variation in Muskox (Ovibos moschatus) Guard Hair Growth Rates: Implications for Measuring Chronological Biomarkers

Segmental analyses of hair may be useful for measuring biomarkers over several seasons to years from a single sample. To attribute hair segments to specific time periods, a known chronological marker, or a hair growth rate, is needed. We examined guard hair growth rates of captive muskoxen (Ovibos moschatus) in Fairbanks, Alaska, USA. We sought to determine if a general growth rate could be applied across muskox populations, thus facilitating the use of segmental analyses for various biomarkers. We used archived samples from 16 muskoxen that had guard hairs sampled at six, 14, and 30 wk after shaving. We measured the lengths of 10 guard hairs per sample, calculated weekly and annual growth rates, and then fitted linear mixed-effects models to assess the effect of different covariates on hair growth rate. The period in which hair had been grown had a significant effect (P<0.05) on growth rate. Extrapolated annual hair growth rates were 277±40 mm/yr (weeks 0-6), 248±47 mm/yr (weeks 7-14), and 165±36 mm/yr (weeks 15-30), with an overall average rate of 210±14 mm/yr. These rates were significantly faster than those of free-ranging Greenland muskoxen-78 mm/yr as measured by stable isotope analyses-and varied intra-annually. This suggests that a universal growth rate cannot be generalized across muskox populations and time.

Population genomics of the muskox' resilience in the near absence of genetic variation

Genomic studies of species threatened by extinction are providing crucial information about evolutionary mechanisms and genetic consequences of population declines and bottlenecks. However, to understand how species avoid the extinction vortex, insights can be drawn by studying species that thrive despite past declines. Here, we studied the population genomics of the muskox (Ovibos moschatus), an Ice Age relict that was at the brink of extinction for thousands of years at the end of the Pleistocene yet appears to be thriving today. We analysed 108 whole genomes, including present-day individuals representing the current native range of both muskox subspecies, the white-faced and the barren-ground muskox (O. moschatus wardi and O. moschatus moschatus) and a ~21,000-year-old ancient individual from Siberia. We found that the muskox' demographic history was profoundly shaped by past climate changes and post-glacial re-colonizations. In particular, the white-faced muskox has the lowest genome-wide heterozygosity recorded in an ungulate. Yet, there is no evidence of inbreeding depression in native muskox populations. We hypothesize that this can be explained by the effect of long-term gradual population declines that allowed for purging of strongly deleterious mutations. This study provides insights into how species with a history of population bottlenecks, small population sizes and low genetic diversity survive against all odds.

Accounting for behaviour in fine-scale habitat selection: A case study highlighting methodological intricacies

Animal habitat selection-central in both theoretical and applied ecology-may depend on behavioural motivations such as foraging, predator avoidance, and thermoregulation. Step-selection functions (SSFs) enable assessment of fine-scale habitat selection as a function of an animal's movement capacities and spatiotemporal variation in extrinsic conditions. If animal location data can be associated with behaviour, SSFs are an intuitive approach to quantify behaviour-specific habitat selection. Fitting SSFs separately for distinct behavioural states helped to uncover state-specific selection patterns. However, while the definition of the availability domain has been highlighted as the most critical aspect of SSFs, the influence of accounting for behaviour in the use-availability design has not been quantified yet. Using a predator-free population of high-arctic muskoxen Ovibos moschatus as a case study, we aimed to evaluate how (1) defining behaviour-specific availability domains, and/or (2) fitting separate behaviour-specific models impacts (a) model structure, (b) estimated selection coefficients and (c) model predictive performance as opposed to behaviour-unspecific approaches. To do so, we first applied hidden Markov models to infer different behavioural modes (resting, foraging, relocating) from hourly GPS positions (19 individuals, 153-1062 observation days/animal). Using SSFs, we then compared behaviour-specific versus behaviour-unspecific habitat selection in relation to terrain features, vegetation and snow conditions. Our results show that incorporating behaviour into the definition of the availability domain primarily impacts model structure (i.e. variable selection), whereas fitting separate behaviour-specific models mainly influences selection strength. Behaviour-specific availability domains improved predictive performance for foraging and relocating models (i.e. behaviours with medium to large spatial displacement), but decreased performance for resting models. Thus, even for a predator-free population subject to only negligible interspecific competition and human disturbance we found that accounting for behaviour in SSFs impacted model structure, selection coefficients and predictive performance. Our results indicate that for robust inference, both a behaviour-specific availability domain and behaviour-specific model fitting should be explored, especially for populations where strong spatiotemporal selection trade-offs are expected. This is particularly critical if wildlife habitat preferences are estimated to inform management and conservation initiatives.

Metabarcoding of fecal pellets in wild muskox populations reveals negative relationships between microbiome and diet alpha diversity

Microbiome diversity and diet composition concomitantly influence species health, fitness, immunity, and digestion. In environments where diet varies spatially and temporally, microbiome plasticity may promote rapid host adaptation to available resources. For northern ungulates in particular, metabarcoding of noninvasively collected fecal pellets presents unprecedented insights into their diverse ecological requirements and niches by clarifying the interrelationships of microbiomes, key to deriving nutrients, in context of altered forage availability in changing climates. Muskoxen (Ovibos moschatus) are Arctic-adapted species that experience fluctuating qualities and quantities of vegetation. Geography and seasonality have been noted to influence microbiome composition and diversity in muskoxen, yet it is unclear how their microbiomes intersect with diet. Following observations from other species, we hypothesized increasing diet diversity would result in higher microbiome diversity in muskoxen. We assessed diet composition in muskoxen using three common plant metabarcoding markers and explored correlations with microbiome data. Patterns of dietary diversity and composition were not fully concordant among the markers used, yet all reflected the primary consumption of willows and sedges. Individuals with similar diets had more similar microbiomes, yet in contrast to most literature, yielded negative relationships between microbiome and diet alpha diversity. This negative correlation may reflect the unique capacities of muskoxen to survive solely on high-fiber Arctic forage and provide insight into their resiliency to exploit changing dietary resources in a rapidly warming Arctic altering vegetation diversity.

Grazing and topography control nutrient pools in low Arctic soils of Southwest Greenland

Soil nutrient pools in the dry low Arctic are likely to be released under climatic change and this bioavailability has the potential to increase both terrestrial and aquatic productions. As well as the direct effect of warming, external disturbances such as nutrient deposition and grazing can also drive ecosystem change. This study in the low Arctic Kangerlussuaq area of southwest Greenland compared soil nutrient pools in terms of both topographic position on a catena and by soil depth in two small catchments with contrasting muskox abundance. We tested the hypotheses that there were differences between soil carbon (C), nitrogen (N) and phosphorus (P) across a soil catena (ridge - slope - valley) and by soil depth (litter - 0-5 cm - 25-30 cm) for the two sites (SS17b, muskox present, versus - SS85, no muskox). Total C and N concentrations of soils were on average lower at SS17b compared to SS85. Moreover, the soil N concentration increased downslope in the catena with higher amounts in the valleys compared to the slopes and ridges. Soil P concentration (0.70 g P kg-1) was similar between catchments; however, litter P content was substantially different. The difference in soil nutrients between the two catchments was most likely due to the presence of muskox at SS17b, and hence grazing associated processes (defecation, altered microbiology and nutrient cycling). This study emphasises the heterogeneity of arctic landscapes and need for ecosystem specific research.

Highlights

Soil nutrient pools in two low-arctic catchments in Greenland were compared.Grazing and dung inputs by muskox affect soil nutrient pools in Greenland.Soil P stores in Kangerlussuaq are similar to intensively managed farmland in Europe.The heterogeneity of arctic landscapes and need for ecosystem-specific research are emphasised.

Qiviut cortisol is associated with metrics of health and other intrinsic and extrinsic factors in wild muskoxen (Ovibos moschatus)

Glucocorticoid (GC) levels are increasingly and widely used as biomarkers of hypothalamic-pituitary-adrenal (HPA) axis activity to study the effects of environmental changes and other perturbations on wildlife individuals and populations. However, identifying the intrinsic and extrinsic factors that influence GC levels is a key step in endocrinology studies to ensure accurate interpretation of GC responses. In muskoxen, qiviut (fine woolly undercoat hair) cortisol concentration is an integrative biomarker of HPA axis activity over the course of the hair's growth. We gathered data from 219 wild muskoxen harvested in the Canadian Arctic between October 2015 and May 2019. We examined the relationship between qiviut cortisol and various intrinsic (sex, age, body condition and incisor breakage) and extrinsic biotic factors (lungworm and gastrointestinal parasite infections and exposure to bacteria), as well as broader non-specific landscape and temporal features (geographical location, season and year). A Bayesian approach, which allows for the joint estimation of missing values in the data and model parameters estimates, was applied for the statistical analyses. The main findings include the following: (i) higher qiviut cortisol levels in males than in females; (ii) inter-annual variations; (iii) higher qiviut cortisol levels in a declining population compared to a stable population; (iv) a negative association between qiviut cortisol and marrow fat percentage; (v) a relationship between qiviut cortisol and the infection intensity of the lungworm Umingmakstrongylus pallikuukensis, which varied depending on the geographical location; and (vi) no association between qiviut cortisol and other pathogen exposure/infection intensity metrics. This study confirmed and further identified important sources of variability in qiviut cortisol levels, while providing important insights on the relationship between GC levels and pathogen exposure/infection intensity. Results support the use of qiviut cortisol as a tool to monitor temporal changes in HPA axis activity at a population level and to inform management and conservation actions.

Geography, seasonality, and host-associated population structure influence the fecal microbiome of a genetically depauparate Arctic mammal

The Canadian Arctic is an extreme environment with low floral and faunal diversity characterized by major seasonal shifts in temperature, moisture, and daylight. Muskoxen (Ovibos moschatus) are one of few large herbivores able to survive this harsh environment. Microbiome research of the gastrointestinal tract may hold clues as to how muskoxen exist in the Arctic, but also how this species may respond to rapid environmental changes. In this study, we investigated the effects of season (spring/summer/winter), year (2007-2016), and host genetic structure on population-level microbiome variation in muskoxen from the Canadian Arctic. We utilized 16S rRNA gene sequencing to characterize the fecal microbial communities of 78 male muskoxen encompassing two population genetic clusters. These clusters are defined by Arctic Mainland and Island populations, including the following: (a) two mainland sampling locations of the Northwest Territories and Nunavut and (b) four locations of Victoria Island. Between these geographic populations, we found that differences in the microbiome reflected host-associated genetic cluster with evidence of migration. Within populations, seasonality influenced bacterial diversity with no significant differences between years of sampling. We found evidence of pathogenic bacteria, with significantly higher presence in mainland samples. Our findings demonstrate the effects of seasonality and the role of host population-level structure in driving fecal microbiome differences in a large Arctic mammal.

Severe Hypoxemia in Muskoxen ( Ovibos moschatus ) Immobilized with Etorphine and Xylazine Corrected with Supplemental Nasal Oxygen

Twenty-three muskoxen ( Ovibos moschatus ) housed in a captive facility for rewilding in Sweden were chemically immobilized for annual health evaluations and hoof trimming. The muskoxen were darted in May to September (2012-15) in their holding pen with etorphine (0.015 mg/kg) and xylazine (0.1 mg/kg) intramuscularly. Twenty-two of the 23 animals were immobilized with a single dart injection. The mean (SD) induction time was 4 (2) min. Arterial blood gases were collected from 18 animals. All animals were severely hypoxemic with varying degrees of respiratory acidosis. The hypoxemia resolved in 17 of 18 animals with intranasal oxygen supplementation at 1 L/min per 100 kg. Relative arterial oxygen saturation (SpO₂) measured by pulse oximetry was significantly higher than the arterial oxygen saturation calculated from the partial pressure of arterial oxygen (SaO₂) obtained by a blood gas analyzer. Based on these findings, muskox can be immobilized successfully with etorphine (0.015 mg/kg) and xylazine (0.1 mg/kg) but should receive supplemental oxygen.