Genetics

Currently there are large areas of the North American landscape that are occupied by Canis spp. hybrids of several varieties, leading to the logical question as to the genetic structure and ecological function of Canis populations across the continent, and to what extent hybrids reflect contemporary landscapes. This study illustrated patterns of niche differentiation between parental canid species and their hybrids using individual high quality genetic profile and species distribution models to support the intermediate phenotype hypothesis. In general, hybrids demonstrated an intermediate habitat suitability compared to its parental species, across most environmental variables used. A similar trend was observed in the niche metric analysis, where we found that hybrids exhibit intermediate niche breadth, with eastern coyotes and eastern wolves exhibiting the broader and narrower niche, respectively. Our results demonstrate that the intermediate phenotype hypothesis is supported even at a large scale and when involving highly mobile large mammal species.

Author Keywords: canid, ecological niche modelling, hybridization, intermediate phenotype, microsatellite genotype, niche differentiation

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

Giardia intestinalis possesses four isotypes of cytochrome b5 (gCYTB-I-IV) that differ from their mammalian counterparts, suggesting different functions in this protozoan parasite. Although the recently discovered Giardia flavoenzyme, GiOR-1, reduces these cytochromes, its properties have not been thoroughly studied, owing to the difficulty in its expression. Here I describe successful conditions for expression of GiOR-1 using autoinduction. GiOR-1 is obtained with flavins bound as indicated by its UV-visible spectrum. Its ability to catalyze electron transfer from donors (NADH, NADPH) to acceptors (oxygen, ferricyanide, cytochrome c, gCYTB5-III) were studied in spectrophotometric rate assays. NADPH is the preferred electron donor, while cytochromes are the preferred electron acceptors. Interestingly, the His-tag used to purify gCYTB5-III decreases its reaction rate with GiOR-1, as an untagged version has slightly faster rates. These findings establish the appropriate conditions for further studies on GiOR-1, including the identification of endogenous electron acceptors.

Author Keywords: Autoinduction, Cytochrome b5, Cytochrome P450 oxidoreductase, Giardia intestinalis, GiOR-1, Polyhistidine tag

Evaluating the impact of climate change is arguably one of the main goals of conservation biology, which can be addressed in part by studying the demographic history of species in the region of interest. In North America, landscape and species composition during the most recent Pleistocene epoch was primarily influenced by glaciation cycles. Glacial advance and retreat caused species ranges to shift as well, leaving signatures of past population bottlenecks in the genetic code of most species. Genomic tools have shown to be important tools for understanding these demographic events to enhance conservation biology measures in several species. In my thesis I first reviewed the state of ungulate genomics, with a focus on how such data sets can be used in understand demography, adaptation, and inform conservation and management. Importantly, the review introduces key analyses like the pairwise sequentially Markovian coalescent and features like variation in antlers and horns and selection pressures that are used throughout subsequent chapters. Using the North American mountain goat as a model species, I then explored the genomic and phenotypic variation in this alpine specialist mammal. Starting with the generation of the first genome assembly for the mountain goat, I identified genes unique to the mountain goat and modeled demographic history going back millions of years using a pairwise sequentially Markovian coalescent approach. Species' effective population size generally paralleled climatic trends over the past one hundred thousand years and severely declined to under a thousand individuals during the last glacial maximum. Given the biological importance of horns in mountain goats and the recent scientific interest in genetic basis of headgear, I analyzed over 23,000 horn records from goats harvested in British Columbia, Alaska and Northwest Territories from 1980 to 2017. Overall, variation in horn size over space and time was low; goats harvested further North had shorter horn lengths and smaller horn circumferences in one year old and 4 years and older age classes and 4 years and older age class, respectively. Proximity of roads, which was used as an indicator of artificial selection, had a small effect on horn size, with larger horns being harvested closer to major roads. Finally, I used two range-wide genomic data sets sequenced with a whole genome re-sequencing and reduced representation approaches to provide estimates of genetic diversity, contemporary effective population sizes and population structure. These insights can help inform management and will potentially make an impact in preserving the mountain goat.

Author Keywords: genome assembly, horn size, Oreamnos americanus, population demography, reduced representation sequencing, whole genome resequencing

Glycine max (soybean) is an economically important plant species that registers a relatively low yield/seed weight compared to other food and oil seed crops due to higher rates of flower and pod abortion. Alleviation of this abortion rate can be achieved by altering the sink strength of the reproductive organs of soybeans. Cytokinin (CK) plays a fundamental role in promoting growth of sink organ (flowers and seeds) by increasing the assimilate demand. Cytokinin oxidase/dehydrogenase (CKX) is an enzyme that catalyses the irreversible breakdown of active CKs and hence reduce the cytokinin content. The current thesis uncovers the members of CKX gene family in soybeans and the natural variations among CKX genes within soybean varieties with different yield characteristics. The identification of null variants of OsCKX2 that resulted in large yield increases by Ashikari et al. (2005) provided a rationale for current thesis. The soybean CKX genes along with the ones from Arabidopsis, Rice and Maize were used to construct a phylogenetic tree. Using comparative phylogeny, protein properties and bioinformatic programs, the potential effect of the identified natural variations on soybean yield was predicted. Five genes among the seventeen soybean CKXs identified, showed polymorphisms. One of the natural variations, A159G, in the gene GmCKX16 occurred close to the active site of the protein and was predicted to affect the activity of enzyme leading to higher accumulation of CKs and hence increased seed weight. Use of such natural variations in marker assisted breeding could lead to the development of higher yielding soybean varieties.

Author Keywords: CKX, Cytokinins, Seed weight, Seed Yield, SNPs, Soybeans

The structure and composition of the landscape can facilitate or impede gene flow, which can have important consequences because genetically isolated groups of individuals may be prone to inbreeding depression and possible extinction. My dissertation examines how landscape structure influences spatial patterns of genetic differentiation and diversity of American marten (Martes americana) and Canada lynx (Lynx canadensis) in Ontario, Canada, and provides methodological advances useful for landscape geneticists. First, I identified the effects of map boundaries on estimates of landscape resistance, and proposed a solution to the bias: a buffer around the map boundary. Second, I assessed the sensitivity of a network-based estimate of genetic distance, conditional genetic distance, to incomplete sampling. I then used these landscape genetic tools in a pairwise, distance-based analysis of 653 martens genotyped at 12 microsatellite loci. I evaluated whether forest management in Ontario has influenced the genetic structure of martens. Although forest management practices had some impact, isolation by distance best described marten gene flow. Our results suggest that managed forests in Ontario are well connected for marten and do not impede marten gene flow. Finally, I used a site-based analysis of 702 lynx genotyped at 14 microsatellite loci to investigate spatial patterns of genetic diversity and differentiation at the trailing (contracting) edge of the lynx distribution in Ontario. I analyzed harvest records and found that the southern edge of lynx range has contracted by >175 km since the 1970s. I also found that neutral genetic diversity decreased towards the trailing edge, whereas genetic differentiation increased. Furthermore, I found strong correlations between gradients of lynx genetic structure and gradients of climate and land cover in Ontario. My findings suggest that increases in winter air temperature, decreases in snow depth, and loss of suitable habitat will result in further loss of genetic diversity in peripheral populations of lynx. Consequently, the adaptive potential of lynx populations on the southern range periphery could decline. In conclusion, my dissertation demonstrates the varying influences that contemporary landscape structure and climate gradients can have on genetic diversity and differentiation of different species.

Author Keywords: Circuitscape, genetic network, landscape genetics, Lynx canadensis, Martes americana, range shift

Interspecific hybridization rates can vary depending on genomic compatibilities between progenitors, while subsequent hybrid spread can vary depending on hybrid performance and habitat availability for hybrid establishment and persistence. As a result, hybridization rates and hybrid frequencies can vary across regions of parental sympatry. In areas around the Laurentian Great Lakes, Typha × glauca is an invasive plant hybrid of native Typha latifolia and introduced Typha angustifolia. In areas of parental sympatry in Atlantic Canada and outside of North America, T. × glauca has been reported as either rare or non-existent. I investigated whether the low frequencies of hybrids documented in Nova Scotia, Atlantic Canada, are influenced by reproductive barriers that prevent hybrid formation or environmental factors (salinity) that reduce hybrid performance. I identified an abundance of hybrids in the Annapolis Valley (inland) and a scarcity of hybrids in coastal wetlands through preliminary site surveys throughout Nova Scotia. In Annapolis Valley populations, flowering times of progenitor species overlapped, indicating that asynchronous flowering times do not limit hybrid formation in this region. Viable progeny were created from interspecific crosses of T. latifolia and T. angustifolia from Nova Scotia, indicating that there are no genomic barriers to fertilization and germination of hybrid seeds. Typha × glauca germination in high salinity was significantly lower than that of T. latifolia, but there was no difference at lower salinities. Therefore, while germination of hybrid seeds may be impeded in the coastal wetlands where salinity is high, inland sites have lower salinity and thus an environment conducive to hybrid germination. However, I found that once established as seedlings, hybrids appear to have greater performance over T. latifolia across all salinities through higher ramet production. Moreover, I found that T. latifolia sourced from Ontario had reduced germination and lower survivorship in high salinities compared to T. latifolia sourced from Nova Scotia, which could indicate local adaptation by T. latifolia to increased salinity. These findings underline that interactions between environment and local progenitor lineages can influence the viability and the consequent distribution and abundance of hybrids. This, in turn, can help explain why hybrids demonstrate invasiveness in some areas of parental sympatry but remain largely absent from other areas.

Author Keywords: flowering phenology, Hybridization, invasive species, physiology, pollen compatibility, salinity tolerance

Understanding the complex genomic architecture underlying quantitative traits can provide valuable insight for the conservation and management of wildlife. Despite improvements in sequencing technologies, few empirical studies have identified quantitative trait loci (QTL) via whole genome sequencing in free-ranging mammal populations outside a few well-studied systems. This thesis uses high-depth whole genome pooled re-sequencing to characterize the molecular basis of the natural variation observed in two sexually selected, heritable traits in white-tailed deer (Odocoileus virginianus, WTD). Specifically, sampled individuals representing the phenotypic extremes from an island population of WTD for antler and body size traits. Our results showed a largely homogenous genome between extreme phenotypes for each trait, with many highly differentiated regions throughout the genome, indicative of a quantitative model for polygenic traits. We identified and validated several potential QTL of putatively small-to-moderate effect for each trait, and discuss the potential for real-world application to conservation and management.

Author Keywords: evolution, extreme phenotypes, genetics, genomics, quantitative traits, sexual selection