Shafer, Aaron B.A.

Natural populations are often difficult and costly to study, due to the plethora of confounding processes and variables present. This is of particular importance when dealing with managed species. Ungulates, for example, act as both consumers and prey sources; they also provide economic benefit through harvest, and as such, are of high ecological and economic value. I addressed conservation and management concerns by quantifying subdivision in wild populations and combined movement with non-invasive sampling to provide novel insight on the physiological drivers of space use in multiple species. This thesis explored biological patterns in ungulates using two distinct approaches: the first used molecular genetics to quantify gene flow, while the second examined the relationship between movement and the gut microbiome using high-throughput sequencing and GPS tracking. The goal of the first chapter was to quantify gene flow and assess the population structure of mountain goats (Oreamnos americanus) in northern British Columbia (BC) to inform management. I used microsatellites to generate genotype data and used a landscape genetics framework to evaluate the possible drivers behind genetic differentiation. The same analyses were performed at both a broad and fine scale, assessing genetic differentiation between populations in all of northern BC and in a case management study area northeast of Smithers BC. The results indicated panmixia among mountain goats regardless of scale, suggesting distance and landscape resistance were minimally inhibiting gene flow. Therefore, management at local scales can continue with little need for genetically informed boundaries, but regulations should be tailored to specific regions incorporating data on local access and harvest pressure. My second chapter aimed to determine the extent to which the gut microbiome drives space-use patterns in a specialist (mountain goat) and generalist (white-tailed deer, Odocoileus virginianus) ungulate. Using fecal samples, we generated genomic data using 16S rRNA high-throughput sequencing to evaluate gut diversity and gut microbiome characteristics. Additionally, individuals were fitted with GPS collars so that we could gain insight into movement patterns. Gut microbiome metrics were stronger predictors of space use and movement patterns with respect to home range size, whereas they were weaker predictors of habitat use. Notably, factors of both the gut microbiome and age of a given species were correlated with changes in space use and habitat use. Ultimately, this research linked high-throughput sequencing and GPS data to better understand ecological processes in wild ungulates.

Author Keywords: gene flow, genomics, gut microbiome, home range, population genetic structure, ungulates

Genetic variation is the raw material and basis for evolutionary changes in nature. The loss of genetic diversity is a challenge many species are facing, with genomics being a potential tool to inform and prioritize decision making. Whole genome analysis can be an asset to conservation biology and the management of species through the generation of more precise and novel metrics. This thesis uses whole genome re-sequencing to characterize the demographic history and quantify genomic metrics relevant to conservation of caribou (Rangifer tarandus) in Québec, Canada. We calculated the ancestral and contemporary patterns of genomic diversity of five representative caribou populations and applied a comparative population genomics framework to assess the interplay between demographic events and genomic diversity. When compared to the census size, NC, the endangered Gaspésie Mountain caribou population had the highest ancestral Ne:NC ratio which is consistent with recent work suggesting high ancestral Ne:NC is of conservation concern. These ratios were highly correlated with genomic signatures (i.e. Tajima's D) of recent population declines and explicit demographic model parameters. Values of contemporary Ne, estimated from linkage-disequilibrium showed Gaspêsie having among the highest contemporary Ne:NC ratio. Importantly, classic conservation genetics theory would predict this population to be of less concern based off this metric alone. Inbreeding measures suggested nuanced patterns of inbreeding and correlated to the demographic models. This study suggests that while the Québec populations are all under decline, they harbour enough ancestral genetic variation to replenish any lost diversity, if conservation decisions are made in favour of these populations, specifically supporting NC.

The integration of genetic and environmental information can help wildlife managers better understand the factors affecting a species' population structure and their response to disturbance. This thesis uses genetic techniques to assess the broad and fine scale population structure of mountain goats in Alaska. The first chapter aims to determine the number of genetically distinct subpopulations and model the demographic history of mountain goats in Alaska. The second chapter investigates the population structure and demographic history of mountain goats in Glacier Bay National Park and examines the impact that climate change will have on these mountain goats. My results indicate that Alaska has eight subpopulations which diverged during the Wisconsin glaciation. In Glacier Bay, population structure is reflective of the landscape during colonization, and mountain goat population density and movement corridors are likely to decline due to future climate change.

Author Keywords: Alaska, biogeography, gene flow, landscape genetics, mountain goat, population genetic structure

Assessments of the adaptive potential of natural populations are essential for understanding and predicting responses to environmental stressors like climate change and infectious disease. The range of stressors species face in a human-dominated landscape, often have contrasting effects. White-tailed deer (Odocoileus virginianus, deer) are expanding in the northern part of their range following decreasing winter severity and increasing forage availability, caused by climate change. Chronic wasting disease (CWD), a prion disease affecting cervids, is likewise expanding and represents a major threat to deer and other cervids We obtained tissue samples from free-ranging deer across their native range in Ontario, Canada which has yet to detect CWD in wild populations of cervids. High throughput sequencing was used to assess neutral genomic variation and variation in the gene responsible for the protein that misfolds into prions when deer contract CWD, known as the PRNP gene. Neutral variation revealed a high number of rare alleles and no population structure, consistent with an expanding population of deer. Functional genetic variation revealed that the frequencies of variants associated to CWD susceptibility and disease progression were evenly distributed across the landscape and the frequencies were consistent with deer populations not infected with CWD. These findings suggest that an observable shift in PRNP allele frequencies likely coincides with the start of a novel CWD epidemic. Sustained surveillance of genomic and genetic variation can be a useful tool for CWD-free regions where deer are managed for ecological and economic benefits.

Author Keywords: Canadian wildlife, population genetics, prion, PRNP, RADseq, ungulate

Wild populations are notoriously difficult to study due to confounding stochastic variables. This thesis tackles two components of investigating wild populations. The first examines the use of niche modeling to quantify macro-scale predator-prey relationships in canid populations across eastern North America, while the second examines range-wide molecular structure in Canada lynx. The goal of the first chapter is to quantify niche characteristics in a Canis hybrid zone of C. lupus, C. lycaon, and C. latrans to better understand the ecological differentiation of these species, and to assess the impacts of incorporating biotic interactions into species distribution models. The goal of the second chapter is to determine if DNA methylation, an epigenetic marker that modifies the structure of DNA, can be used to differentiate populations, and might be a signature of local adaptation. Our results indicated that canids across the hybrid zone in eastern North America exhibit low levels of genetic and ecological differentiation, and that the importance of biotic interactions are largely lost at large spatial scales. We also identified cryptic structure in methylation patterns in Canada lynx populations, which suggest signatures of local adaptation, and indicate the utility of DNA methylation as a marker for investigating adaptive divergence.

Author Keywords: Ecological Epigenetics, Ecological Genetics, SDM