Muskox status, recent variation, and uncertain future

The geographic distribution, and the biotic and abiotic predictors of select zoonotic pathogen detections in Canadian polar bears

Increasing Arctic temperatures are facilitating the northward expansion of more southerly hosts, vectors, and pathogens, exposing naïve populations to pathogens not typical at northern latitudes. To understand such rapidly changing host-pathogen dynamics, we need sensitive and robust surveillance tools. Here, we use a novel multiplexed magnetic-capture and droplet digital PCR (ddPCR) tool to assess a sentinel Arctic species, the polar bear (Ursus maritimus; n = 68), for the presence of five zoonotic pathogens (Erysipelothrix rhusiopathiae, Francisella tularensis, Mycobacterium tuberculosis complex, Toxoplasma gondii and Trichinella spp.), and observe associations between pathogen presence and biotic and abiotic predictors. We made two novel detections: the first detection of a Mycobacterium tuberculosis complex member in Arctic wildlife and the first of E. rhusiopathiae in a polar bear. We found a prevalence of 37% for E. rhusiopathiae, 16% for F. tularensis, 29% for Mycobacterium tuberculosis complex, 18% for T. gondii, and 75% for Trichinella spp. We also identify associations with bear age (Trichinella spp.), harvest season (F. tularensis and MTBC), and human settlements (E. rhusiopathiae, F. tularensis, MTBC, and Trichinella spp.). We demonstrate that monitoring a sentinel species, the polar bear, could be a powerful tool in disease surveillance and highlight the need to better characterize pathogen distributions and diversity in the Arctic.

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.

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.

The Effect of Species and Sex on the Element Content of Muskox (Ovibos moschatus) and Caribou (Rangifer tarandus groenlandicus) Tissues

Muskox (Ovibos moschatus) and caribou (Rangifer tarandus groenlandicus) are wild ruminants that inhabit the Greenland tundra. They are part of the diet of many Greenlanders, being important sources of protein and micronutrients such as iron. The objective of this study is to analyse the element profiles of three tissues from these species: skeletal muscle, liver and adipose tissue, and to determine if they are affected by species and sex (male vs. female). Samples were obtained from annual hunts in two different regions of West Greenland. Element profiles were analysed using inductively-coupled plasma-optical emission spectrometry. The interaction between species and sex was only detected in Na (sodium) in the muscle and adipose tissue, where male and female caribou had the highest concentrations, respectively. The effect of sex was not statistically significant in the liver samples and only occasionally in the other tissues. Species was the most relevant factor in element profiles found in this study. Caribou had higher concentrations of K (potassium) and S (sulphur) in the muscle and liver. Fewer differences were detected between species in the adipose tissue, compared to the other tissues. These differences may reflect the feeding behaviour and the geographical location of both species. This study contributes to evaluate the element composition of the edible tissue of these wild ungulate species, as well as evaluating the factors of sex and species that could differentiate their composition.

Molecular characterization of Blastocystis and Entamoeba of muskoxen and sheep in Greenland

Molecular characterisation of endobionts that are shared among human and non-human hosts can help shed light on the epidemiology and inform studies that aim to unravel the role of these organisms in health and disease. Two of the most common of shared endobionts include the single-celled intestinal protists Blastocystis and Entamoeba. Here, we present the first known data on genetic diversity and host specificity of these two genera in Greenland. Faecal DNA samples from 243 muskoxen and 44 sheep were submitted to metabarcoding of nuclear small subunit ribosomal DNA. Entamoeba- and Blastocystis-specific sequences were clustered, and consensus sequences were subjected to taxonomic query. Using MinION-based sequencing, near-complete nuclear small subunit ribosomal DNA sequences were obtained from four faecal samples. Of the 243 muskox samples, 180 (74%) and 19 (8%) were positive for Blastocystis and Entamoeba, respectively. Forty (91%) and six (14%) of the 44 sheep samples were positive for Blastocystis and Entamoeba, respectively. Blastocystis subtypes (ST) 10, 14, 21, 24-26, and a novel subtype (ST40) were identified. Colonisation by more than one subtype was common. ST40 was common in muskoxen but limited to Northeast Greenland. Entamoeba bovis and the E. bovis-associated ribosomal lineages (RL) 1 and 8 were found, and three conditional lineages (CL) 3, 4, and 10 were confirmed; CL10 was promoted to RL12. Several novel lineages were identified, all of which were linked to the E. bovis complex. In conclusion, Blastocystis was far more common than Entamoeba and found in approximately three of every four animals; both can be considered common colonisers of large herbivorous mammals in Greenland. Multiple subtypes/lineages of both genera were commonly observed, some of which were novel, but most of which are seen in many other parts of the world.

Next-generation sequencing approach to investigate genome variability of Parapoxvirus in Canadian muskoxen (Ovibos moschatus)

In 2016, the first orf virus, a double-stranded DNA (dsDNA) virus of the genus parapoxvirus, from a muskox was isolated on Victoria Island, Nunavut (NU), Canada. We used deep sequencing on DNA extracted from orf virus-positive tissues from wild muskoxen from locations on Victoria Island and the adjacent mainland. Orf virus sequence reads derived from four samples were nearly identical. The consensus sequences generated from pooled reads of MxOV comprises of a large contiguous sequence (contig) of 131,759 bp and a smaller right terminal contig of 3552 bp, containing all coding sequences identified as Parapoxvirus. Individual gene comparisons reveal that MxOV shares genetic characteristics with reference strains from both sheep and goat origin. Recombination analysis using Bootscan, MAXCHI, GENECONV, CHIMAERA, SISCAN, and RDP algorithms within the RDP4 software predicted recombination events in two virulence factors, and a large 3000 bp segment of the MxOV genome. Partial B2L nucleotide sequences from strains around the world and other North American isolates were compared to MxOV using MUSCLE alignments and RAxML phylogenetic trees. MxOV was identical to our previously characterized isolate, and shared similarity with orf virus isolated from sheep and goats. The phylogenetic grouping of partial B2L nucleotide sequences did not follow the sample geographic distribution. More full genomes of orf virus, or at least full B2L gene squences, in wildlife are needed especially in North America to better understand the epidemiology of the disease in muskoxen.

Spatiotemporal influences of climate and humans on muskox range dynamics over multiple millennia

Processes leading to range contractions and population declines of Arctic megafauna during the late Pleistocene and early Holocene are uncertain, with intense debate on the roles of human hunting, climatic change, and their synergy. Obstacles to a resolution have included an overreliance on correlative rather than process-explicit approaches for inferring drivers of distributional and demographic change. Here, we disentangle the ecological mechanisms and threats that were integral in the decline and extinction of the muskox (Ovibos moschatus) in Eurasia and in its expansion in North America using process-explicit macroecological models. The approach integrates modern and fossil occurrence records, ancient DNA, spatiotemporal reconstructions of past climatic change, species-specific population ecology, and the growth and spread of anatomically modern humans. We show that accurately reconstructing inferences of past demographic changes for muskox over the last 21,000 years require high dispersal abilities, large maximum densities, and a small Allee effect. Analyses of validated process-explicit projections indicate that climatic change was the primary driver of muskox distribution shifts and demographic changes across its previously extensive (circumpolar) range, with populations responding negatively to rapid warming events. Regional analyses show that the range collapse and extinction of the muskox in Europe (~13,000 years ago) was likely caused by humans operating in synergy with climatic warming. In Canada and Greenland, climatic change and human activities probably combined to drive recent population sizes. The impact of past climatic change on the range and extinction dynamics of muskox during the Pleistocene-Holocene transition signals a vulnerability of this species to future increased warming. By better establishing the ecological processes that shaped the distribution of the muskox through space and time, we show that process-explicit macroecological models have important applications for the future conservation and management of this iconic species in a warming Arctic.

Herbivores in Arctic ecosystems: Effects of climate change and implications for carbon and nutrient cycling

Arctic terrestrial herbivores influence tundra carbon and nutrient dynamics through their consumption of resources, waste production, and habitat-modifying behaviors. The strength of these effects is likely to change spatially and temporally as climate change drives shifts in herbivore abundance, distribution, and activity timing. Here, we review how herbivores influence tundra carbon and nutrient dynamics through their consumptive and nonconsumptive effects. We also present evidence for herbivore responses to climate change and discuss how these responses may alter the spatial and temporal distribution of herbivore impacts. Several current knowledge gaps limit our understanding of the changing functional roles of herbivores; these include limited characterization of the spatial and temporal variability in herbivore impacts and of how herbivore activities influence the cycling of elements beyond carbon. We conclude by highlighting approaches that will promote better understanding of herbivore effects on tundra ecosystems, including their integration into existing biogeochemical models, new applications of remote sensing techniques, and the continued use of distributed experiments.

Wolverines (Gulo gulo) in the Arctic: Revisiting distribution and identifying research and conservation priorities amid rapid environmental change

Wolverines (Gulo gulo) occupy most of the globe’s Arctic tundra. Given the rapidly warming climate and expanding human activity in this biome, understanding wolverine ecology, and therefore the species’ vulnerability to such changes, is increasingly important for developing research priorities and effective management strategies. Here, we review and synthesize knowledge of wolverines in the Arctic using both Western science sources and available Indigenous Knowledge (IK) to improve our understanding of wolverine ecology in the Arctic and better predict the species’ susceptibility to change. To accomplish this, we update the pan-Arctic distribution map of wolverines to account for recent observations and then discuss resulting inference and uncertainties. We use these patterns to contextualize and discuss potential underlying drivers of distribution and population dynamics, drawing upon knowledge of food habits, habitat associations, and harvest, as well as studies of wolverine ecology elsewhere. We then identify four broad areas to prioritize conservation and research efforts: (1) Monitoring trends in population abundance, demographics, and distribution and the drivers thereof, (2) Evaluating and predicting wolverines’ responses to ongoing climate change, particularly the consequences of reduced snow and sea ice, and shifts in prey availability, (3) Understanding wolverines’ response to human development, including the possible impact of wintertime over-snow travel and seismic testing to reproductive denning, as well as vulnerability to hunting and trapping associated with increased human access, and (4) Ensuring that current and future harvest are sustainable.

Prehistoric perspectives can help interpret the present: 14,000 years of moose (Alces alces L) in the Western Arctic

Rapidly changing climate at high latitudes has triggered a search for bellwethers of ecological change there. If the initial signs of change can be identified, perhaps we can predict where these changes will lead. Large-bodied, terrestrial, herbivores are potential candidates for bellwether taxa because of the key roles they play in some ecological communities. Here we assemble historical, archaeological and paleontological records of moose (<i>Alces alces</i> Linnaeus, 1758.) from the western Arctic and Subarctic. Results show that rather than having recently invaded tundra regions in response to post-Little Ice Age warming, moose have inhabited river corridors several hundred kilometers north of the closed, boreal forest since they first colonized North America across the Bering Land Bridge ca. 14,000 years ago. The combination of high mobility, fluctuation-prone metapopulations, and reliance on early successional vegetation makes changes in the northern range limits of moose undependable bellwethers for other biotic responses to changing climate. The history of moose at high latitudes illustrates how understanding what happened in prehistory is useful for correctly assigning significance and cause to present-day ecological changes.

Projected bioclimatic distributions in Nearctic Bovidae signal the potential for reduced overlap with protected areas

Assumptions about factors such as climate in shaping species' realized and potential distributions underlie much of conservation planning and wildlife management. Climate and climatic change lead to shifts in species distributions through both direct and indirect ecological pressures. Distributional shifts may be particularly important if range overlap is altered between interacting species, or between species and protected areas. The cattle family (Bovidae) represents a culturally, economically, and ecologically important taxon that occupies many of the world's rangelands. In contemporary North America, five wild bovid species inhabit deserts, prairies, mountains, and tundra from Mexico to Greenland. Here, we aim to understand how future climate change will modify environmental characteristics associated with North American bovid species relative to the distribution of extant protected areas. We fit species distribution models for each species to climate, topography, and land cover data using observations from a citizen science dataset. We then projected modeled distributions to the end of the 21st century for each bovid species under two scenarios of anticipated climate change. Modeling results suggest that suitable habitat will shift inconsistently across species and that such shifts will lead to species-specific variation in overlap between potential habitat and existing protected areas. Furthermore, projected overlap with protected areas was sensitive to the warming scenario under consideration, with diminished realized protected area under greater warming. Conservation priorities and designation of new protected areas should account for ecological consequences of climate change.

Muskoxen homogenise soil microbial communities and affect the abundance of methanogens and methanotrophs

Grazing herbivores may affect soil microbial communities indirectly by impacting soil structure and vegetation composition. In high arctic environments, this impact is poorly elucidated, while having potentially wide-reaching effects on the ecosystem. This study examines how a key arctic herbivore, the muskox Ovibos moschatus, affects the soil microbial community in a high arctic fen. Environmental DNA was extracted from soil samples taken from grazed control plots and from muskox exclosures established 5 years prior. We sequenced amplicons of the 16S rRNA gene to provide insight into the microbial communities. We found that in the grazed control plots, microbial communities exhibited high evenness and displayed highly similar overall diversity. In plots where muskoxen had been excluded, microbial diversity was significantly reduced, and had more uneven intra-sample populations and overall lower ecological richness and evenness. We observed that the composition of microbial communities in grazed soils were significantly affected by the presence of muskoxen, as seen by elevated relative abundances of Bacteroides and Firmicutes, two major phyla found in muskox faeces. Furthermore, an increase in relative abundance of bacteria involved in degradation of recalcitrant carbohydrates and cycling of nitrogen was observed in grazed soil. Ungrazed soils displayed increased abundances of bacteria potentially involved in anaerobic oxidation of methane, whereas some methanogens were more abundant in grazed soils. This corroborates a previous finding that methane emissions are higher in arctic fens under muskox grazing. Our results show that the presence of large herbivores stimulates soil microbial diversity and has a homogenizing influence on the inter-species dynamics in soil microbial communities. The findings of this study, thus, improve our understanding of the effect of herbivore grazing on arctic ecosystems and the derived methane cycling.

Fecal microbiota and diets of muskox female adults and calves

In mammals, the gut microbiome is vertically transmitted during maternal lactation at birth. In this study, we investigated the gut microbiome and diets of muskox, a large herbivore inhabiting in the high Arctic. We compared the microbiota composition using bacterial 16S rRNA gene sequencing and diets using stable isotope analysis of muskox feces of six female adults and four calves on Ella Island, East Greenland. Firmicutes were the most abundant bacterial phylum in both the adults and calves, comprising 94.36% and 94.03%, respectively. Significant differences were observed in the relative abundance of the two Firmicutes families. The adults were primarily dominated by Ruminococcaceae (73.90%), and the calves were dominated by both Ruminococcaceae (56.25%) and Lachnospiraceae (24.00%). Stable isotope analysis of the feces in the study area revealed that both adults and calves had similar ranges of ¹³C and ¹⁵N, likely derived from the dominant diet plants. Despite their similar diets, the different gut microbiome compositions in muskox adults and calves indicate that the gut microbiome of the calves may not be fully colonized to the extent of that of the adults.

Draft Genome Assembly of an Iconic Arctic Species: Muskox (Ovibos moschatus)

Muskoxen (Ovibos moschatus) are Arctic species within the Caprinae subfamily that are economically and culturally significant to northern Indigenous communities. Low genetic diversity from repeated genetic bottlenecks, coupled with the effects of Arctic warming (e.g., heat stress, changing forage, pathogen range expansions), present conservation concerns for this species. Reference genome assemblies enhance our ecological and evolutionary understanding of species (which in turn aid conservation efforts). Herein, we provide a full draft reference genome of muskox using Illumina Hiseq data and cross-species scaffolding. The final reference assembly yielded a genome of 2,621,890,883 bp in length, a scaffold N50 of ~13.2 million, and an annotation identifying ~19.3 k genes. The muskox genome assembly and annotation were then used to reconstruct a phylogenetic tree which estimated muskoxen diverged from other ungulate species~12 Mya. To gain insight into the demographic history of muskoxen we also performed pairwise sequentially Markovian coalescent (PSMC) that identified two population bottlenecks coinciding with major glaciation events contributing to the notoriously low genetic variation observed in muskoxen. Overall, this genome assembly provides a foundation for future population genomic studies, such as latitudinal analyses, to explore the capacity of muskoxen to adapt to rapidly changing environments.

Large herbivores facilitate the persistence of rare taxa under tundra warming

Ecological rarity, characterized by low abundance or limited distribution, is typical of most species, yet our understanding of what factors contribute to the persistence of rare species remains limited. Consequently, little is also known about whether rare species might respond differently than common species to direct (e.g., abiotic) and indirect (e.g., biotic) effects of climate change. We investigated the effects of warming and exclusion of large herbivores on 14 tundra taxa, three of which were common and 11 of which were rare, at an inland, low-arctic study site near Kangerlussuaq, Greenland. Across all taxa, pooled commonness was reduced by experimental warming, and more strongly under herbivore exclusion than under herbivory. However, taxon-specific analyses revealed that although warming elicited variable effects on commonness, herbivore exclusion disproportionately reduced the commonness of rare taxa. Over the 15-year duration of the experiment, we also observed trends in commonness and rarity under all treatments through time. Sitewide commonness increased for two common taxa, the deciduous shrubs Betula nana and Salix glauca, and declined in six other taxa, all of which were rare. Rates of increase or decline in commonness (i.e., temporal trends over the duration of the experiment) were strongly related to baseline commonness of taxa early in the experiment under all treatments except warming with grazing. Hence, commonness itself may be a strong predictor of species' responses to climate change in the arctic tundra biome, but large herbivores may mediate such responses in rare taxa, perhaps facilitating their persistence.

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.

Congenital cataract and spherophakia leading to starvation in a free-ranging muskox neonate from the Northwest Territories, Canada

A muskox neonate (Ovibos moschatus) that died of starvation was diagnosed with congenital lenticular anomalies that included spherophakia and hypermature cataract associated with probable lens-induced lymphocytic uveitis and neutrophilic keratitis. Impaired sight as a result of cataract and associated inflammation likely contributed to abandonment and starvation, although maternal death cannot be excluded definitively. Ocular lesions, such as congenital cataracts and spherophakia in neonates, may be important factors affecting survival in free-ranging animals.

Large and small herbivores have strong effects on tundra vegetation in Scandinavia and Alaska

Large and small mammalian herbivores are present in most vegetated areas in the Arctic and often have large impacts on plant community composition and ecosystem functioning. The relative importance of different herbivores and especially how their specific impact on the vegetation varies across the Arctic is however poorly understood.Here, we investigate how large and small herbivores influence vegetation density and plant community composition in four arctic vegetation types in Scandinavia and Alaska. We used a unique set of exclosures, excluding only large (reindeer and muskoxen) or all mammalian herbivores (also voles and lemmings) for at least 20 years.We found that mammalian herbivores in general decreased leaf area index, NDVI, and abundance of vascular plants in all four locations, even though the strength of the effect and which herbivore type caused these effects differed across locations. In three locations, herbivore presence caused contrasting plant communities, but not in the location with lowest productivity. Large herbivores had a negative effect on plant height, whereas small mammalian herbivores increased species diversity by decreasing dominance of the initially dominating plant species. Above- or belowground disturbances caused by herbivores were found to play an important role in shaping the vegetation in all locations.Synthesis: Based on these results, we conclude that both small and large mammalian herbivores influence vegetation in Scandinavia and Alaska in a similar way, some of which can mitigate effects of climate change. We also see important differences across locations, but these depend rather on local herbivore and plant community composition than large biogeographical differences among continents.

Zoonotic pathogens in wild muskoxen (Ovibos moschatus) and domestic sheep (Ovis aries) from Greenland

The present study aimed to estimate the prevalence of zoonotic pathogens Giardia duodenalis, Cryptosporidium spp., Toxoplasma gondii and Erysipelothrix in muskoxen (Ovibos moschatus) and sheep (Ovis aries) from Greenland. In 2017 and 2018, faecal samples were collected from wild muskoxen from three distinct populations (Zackenberg, Kangerlussuaq, and Ivittuut) and from domestic sheep from southwest Greenland. Blood samples were collected from muskoxen from Kangerlussuaq and Ivittuut and from sheep. Faecal samples were tested for specific DNA of G. duodenalis and Cryptosporidium spp., and blood samples were tested for antibodies against T. gondii and Erysipelothrix. The estimated prevalence of G. duodenalis was 0% (0/58), 17% (7/41) and 0% (0/55) in muskoxen from Zackenberg, Kangerlussuaq and Ivittuut, respectively, and 37% (16/43) in sheep. The estimated prevalence of Cryptosporidium was 0% (0/58), 2% (1/41), 7% (4/55) in muskoxen from Zackenberg, Kangerlussuaq, Ivittuut, respectively, and 2% (1/43) in sheep. Neither Giardia nor Cryptosporidium were detected in winter samples (0/78). Of the positive samples, Giardia from one muskox sample only was successfully typed as G. duodenalis assemblage A, and Cryptosporidium from two muskoxen was successfully typed as C. parvum, subtype IIdA20G1e. The estimated T. gondii seroprevalence was 2% (1/44) and 0% (0/8) in muskoxen from Kangerlussuaq and Ivittuut, respectively, and 1% (1/155) in sheep. The estimated Erysipelothrix seroprevalence was 2% (1/45) and 13% (1/8) in muskoxen from Kangerlussuaq and Ivittuut, respectively, and 7% (10/150) in sheep. The results of this study add to the scarce knowledge on zoonotic pathogens in the Arctic.

Corynebacterium freneyi Bacterial Septicemia Secondary to Contagious Ecthyma in a Wild Muskox (Ovibos moschatus)

A wild muskox (Ovibos moschatus) with dermatitis typical of contagious ecthyma had secondary bacterial septicemia with Corynebacterium freneyi that included laminitis, hepatitis, and suppurative encephalitis. This case supports the association between orf virus infection and fatal secondary infections, which may have contributed to population declines on Victoria Island, Canada.

Qiviut cortisol reflects hypothalamic-pituitary-adrenal axis activity in muskoxen (Ovibos moschatus)

Muskoxen (Ovibos moschatus) are increasingly exposed to a broad diversity of stressors in their rapidly changing Arctic environment. There is an urgent need to develop validated tools to monitor the impact of these stressors on the hypothalamic-pituitary-adrenal (HPA) axis activity of muskoxen to help inform conservation actions. Here, we evaluated whether muskox qiviut (dense wooly undercoat) cortisol accurately reflects changes in HPA axis activity. Two repeated pharmacological challenges, involving weekly administrations of saline (control group) or adrenocorticotropic hormone (ACTH) during five consecutive weeks, were done on captive muskoxen, in winter (no hair growth) and summer (maximum hair growth). Pre-challenge qiviut cortisol levels were significantly higher in the shoulder than in the neck, but neither differed from rump concentrations. Qiviut cortisol levels significantly increased (p < 0.001) in response to the administration of ACTH during the hair growth phase, but not in the absence of growth (p = 0.84). Cortisol levels in the qiviut segment grown during the summer challenge increased significantly over a six-month period in the ACTH-injected muskoxen with a similar trend occurring in the control animals. Finally, cortisol levels in shed qiviut were significantly higher and not correlated to those of fully grown qiviut shaved three months earlier. Our results show that cortisol is deposited in qiviut during its growth and that qiviut cortisol can thus be used as an integrated measure of HPA axis activity over the period of the hair's growth. Differences in qiviut cortisol across body regions, significant differences in qiviut segments over time, and differences between shed qiviut versus unshed qiviut, highlight the importance of consistent design and methodology for sample collection and analyses in order to account for sources of variation when using qiviut cortisol as a biomarker of HPA axis activity in muskoxen.

Fecal glucocorticoid metabolites reflect hypothalamic-pituitary-adrenal axis activity in muskoxen (Ovibos moschatus)

Muskoxen (Ovibos moschatus), a taxonomically unique Arctic species, are increasingly exposed to climate and other anthropogenic changes. It is critical to develop and validate reliable tools to monitor their physiological stress response in order to assess the impacts of these changes. Here, we measured fecal glucocorticoid metabolite (FGM) levels in response to the administration of adrenocorticotropic hormone (ACTH) in the winter (1 IU/kg) and summer (2 IU/kg) using two enzyme immunoassays, one targeting primarily cortisol and the other targeting primarily corticosterone. Fecal cortisol levels varied substantially within and among individuals, and none of the animals in either challenge showed an increase in fecal cortisol following the injection of ACTH. By contrast, two of six (winter) and two of five (summer) muskoxen showed a clear response in fecal corticosterone levels (i.e., maximal percentage increase as compared to time 0 levels > 100%). Increases in fecal corticosterone post-ACTH injection occurred earlier and were of shorter duration in the summer than in the winter and fecal corticosterone levels were, in general, lower during the summer. These seasonal differences in FGM responses may be related to the use of different individuals (i.e., influence of sex, age, social status, etc.) and to seasonal variations in the metabolism and excretion of glucocorticoids, intestinal transit time, voluntary food intake, and fecal output and moisture content. Results from this study support using FGMs as a biomarker of hypothalamic–pituitary–adrenal axis activity in muskoxen, advance our understanding of the physiological adaptations of mammals living in highly seasonal and extreme environments such as the Arctic, and emphasize the importance of considering seasonality in other species when interpreting FGM levels.

Range expansion of muskox lungworms track rapid arctic warming: implications for geographic colonization under climate forcing

Rapid climate warming in the Arctic results in multifaceted disruption of biodiversity, faunal structure, and ecosystem health. Hypotheses have linked range expansion and emergence of parasites and diseases to accelerating warming globally but empirical studies demonstrating causality are rare. Using historical data and recent surveys as baselines, we explored climatological drivers for Arctic warming as determinants of range expansion for two temperature-dependent lungworms, Umingmakstrongylus pallikuukensis and Varestrongylus eleguneniensis, of muskoxen (Ovibos moschatus) and caribou (Rangifer tarandus), in the Canadian Arctic Archipelago from 1980 through 2017. Our field data shows a substantial northward shift of the northern edge of the range for both parasites and increased abundance across the expanded ranges during the last decade. Mechanistic models parameterized with parasites’ thermal requirements demonstrated that geographical colonization tracked spatial expansion of permissive environments, with a temporal lag. Subtle differences in life histories, thermal requirements of closely related parasites, climate oscillations and shifting thermal balances across environments influence faunal assembly and biodiversity. Our findings support that persistence of host-parasite assemblages reflects capacities of parasites to utilize host and environmental resources in an ecological arena of fluctuating opportunity (alternating trends in exploration and exploitation) driving shifting boundaries for distribution across spatial and temporal scales.

Environmental conditions alter behavioural organization and rhythmicity of a large Arctic ruminant across the annual cycle

The existence and persistence of rhythmicity in animal activity during phases of environmental change is of interest in ecology, evolution and chronobiology. A wide diversity of biological rhythms in response to exogenous conditions and internal stimuli have been uncovered, especially for polar vertebrates. However, empirical data supporting circadian organization in behaviour of large ruminating herbivores remains inconclusive. Using year-round tracking data of the largest Arctic ruminant, the muskox (Ovibos moschatus), we modelled rhythmicity as a function of behaviour and environmental conditions. Behavioural states were classified based on patterns in hourly movements, and incorporated within a periodicity analyses framework. Although circadian rhythmicity in muskox behaviour was detected throughout the year, ultradian rhythmicity was most prevalent, especially when muskoxen were foraging and resting in mid-winter (continuous darkness). However, when combining circadian and ultradian rhythmicity together, the probability of behavioural rhythmicity declined with increasing photoperiod until largely disrupted in mid-summer (continuous light). Individuals that remained behaviourally rhythmic during mid-summer foraged in areas with lower plant productivity (NDVI) than individuals with arrhythmic behaviour. Based on our study, we conclude that muskoxen may use an interval timer to schedule their behavioural cycles when forage resources are low, but that the importance and duration of this timer are reduced once environmental conditions allow energetic reserves to be replenished ad libitum. We argue that alimentary function and metabolic requirements are critical determinants of biological rhythmicity in muskoxen, which probably applies to ruminating herbivores in general.

Novel insights into serodiagnosis and epidemiology of Erysipelothrix rhusiopathiae, a newly recognized pathogen in muskoxen (Ovibos moschatus)

BACKGROUND

Muskoxen are a key species of Arctic ecosystems and are important for food security and socio-economic well-being of many Indigenous communities in the Arctic and Subarctic. Between 2009 and 2014, the bacterium Erysipelothrix rhusiopathiae was isolated for the first time in this species in association with multiple mortality events in Canada and Alaska, raising questions regarding the spatiotemporal occurrence of the pathogen and its potential impact on muskox populations.

MATERIALS AND METHODS

We adapted a commercial porcine E. rhusiopathiae enzyme-linked immunosorbent assay to test 958 blood samples that were collected from muskoxen from seven regions in Alaska and the Canadian Arctic between 1976 and 2017. The cut-off between negative and positive results was established using mixture-distribution analysis, a data-driven approach. Based on 818 samples for which a serological status could be determined and with complete information, we calculated trends in sample seroprevalences in population time-series and compared them with population trends in the investigated regions.

RESULTS

Overall, 219/818 (27.8%, 95% Confidence Interval: 24.7-31.0) samples were classified as positive for exposure to E. rhusiopathiae. There were large variations between years and regions. Seropositive animals were found among the earliest serum samples tested; 1976 in Alaska and 1991 in Canada. In Alaskan muskoxen, sample seroprevalence increased after 2000 and, in two regions, peak seroprevalences occurred simultaneously with population declines. In one of these regions, concurrent unusual mortalities were observed and E. rhusiopathiae was isolated from muskox carcasses. In Canada, there was an increase in sample seroprevalence in two muskox populations following known mortality events that had been attributed to E. rhusiopathiae.

CONCLUSION

Our results indicate widespread exposure of muskoxen to E. rhusiopathiae in western Canada and Alaska. Although not new to the Arctic, we documented an increased exposure to the pathogen in several regions concurrent with population declines. Understanding causes for the apparent increased occurrence of this pathogen and its association with large scale mortality events for muskoxen is critical to evaluate the implications for wildlife and wildlife-dependent human populations in the Arctic.

Arctic terrestrial biodiversity status and trends: A synopsis of science supporting the CBMP State of Arctic Terrestrial Biodiversity Report

This review provides a synopsis of the main findings of individual papers in the special issue Terrestrial Biodiversity in a Rapidly Changing Arctic. The special issue was developed to inform the State of the Arctic Terrestrial Biodiversity Report developed by the Circumpolar Biodiversity Monitoring Program (CBMP) of the Conservation of Arctic Flora and Fauna (CAFF), Arctic Council working group. Salient points about the status and trends of Arctic biodiversity and biodiversity monitoring are organized by taxonomic groups: (1) vegetation, (2) invertebrates, (3) mammals, and (4) birds. This is followed by a discussion about commonalities across the collection of papers, for example, that heterogeneity was a predominant pattern of change particularly when assessing global trends for Arctic terrestrial biodiversity. Finally, the need for a comprehensive, integrated, ecosystem-based monitoring program, coupled with targeted research projects deciphering causal patterns, is discussed.

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.

Already at the bottom? Demographic declines are unlikely further to undermine genetic diversity of a large Arctic ungulate: muskox, Ovibos moschatus (Artiodactyla: Bovidae)

Low genetic diversity is associated with low fitness and evolutionary potential, yet the demographic and life-history traits of some species contribute to low genetic diversity, without empirical evidence of negative impacts on fitness. Modelling past and future trajectories of genetic diversity under different demographic scenarios can provide insight into how genetic variation might impact population fitness. The muskox is an Arctic species that has undergone multiple population bottlenecks and, although populations have rebounded repeatedly, two large populations have recently declined by &gt; 50%. It is unclear how these demographic patterns influence muskox genetic diversity and fitness. We compared the genetic diversity of Canadian muskox populations undergoing opposing population trends. Genotyping 84 mainland and 244 Victoria Island individuals at ten microsatellite loci revealed low genetic variation (Victoria Island, mean allelic richness 1.66, expected heterozygosity 0.16; mainland, mean allelic richness 2.58, expected heterozygosity 0.41), with no evidence of further reductions in diversity subsequent to recent demographic declines. Bayesian modelling showed that a 1900s bottleneck contributed to the lack of diversity in contemporary populations, and forward-in-time simulations suggested little effect on genetic diversity over the next 100 years. Muskoxen might have reached a genetic diversity minimum, and additional research will be needed to determine their capacity to adapt to rapid changes in selective pressures in a rapidly changing Arctic.

A Transdisciplinary Approach to Brucella in Muskoxen of the Western Canadian Arctic 1989-2016

Brucella serostatus was evaluated in 3189 muskoxen sampled between 1989 and 2016 from various locations of the Canadian Arctic archipelago and mainland, near the communities of Sachs Harbour and Ulukhaktok, Northwest Territories, and Cambridge Bay and Kugluktuk, Nunavut. Brucella antibodies were found only in muskoxen sampled around Cambridge Bay, both on southern Victoria Island and on the adjacent mainland (Kent Peninsula). Consistent with participatory epidemiology data documented from local harvesters describing increased Brucella-like syndromes (swollen joints and lameness) and a decreased proportion of juveniles, the apparent Brucella seroprevalence in the sampled muskoxen of the Cambridge Bay area increased from 0.9% (95% CI 0.3-2.1) in the period of 1989-2001 to 5.6% (95% CI 3.3-8.9) in 2010-2016. The zoonotic bacteria Brucella suis biovar 4 was also cultured from tissues of muskoxen sampled on Victoria Island near Ulukhaktok in 1996 (n = 1) and Cambridge Bay in 1998, 2014, and 2016 (n = 3). Overall, our data demonstrate that B. suis biovar 4 is found in muskoxen that are harvested for food and by guided hunts on Victoria Island and Kent Peninsula, adding an important public health dimension to this study. Robust participatory epidemiology data on muskox health and diseases greatly enhanced the interpretation of our Cambridge Bay data and, combined with the serological and microbiological data, provide compelling evidence that the prevalence of B. suis biovar 4 has increased in this area since the late 1990s. This study enhances the available knowledge on Brucella exposure and infection in muskoxen and provides an example of how scientific knowledge and local knowledge can work together to better understand disease status in wildlife.