Ecology

The objective of this thesis is to better understand the demography and habitat selection of Newfoundland caribou. Chapter 1 provides a general introduction of elements of population ecology and behavioural ecology discussed in the thesis. In Chapter 2, I examine the causes of long-term fluctuations among caribou herds. My findings indicate that winter severity and density-dependent degradation of summer range quality offer partial explanations for the observed patterns of population change. In Chapter 3, I investigate the influence of climate, predation and density-dependence on cause-specific neonate survival. I found that when caribou populations are in a period of increase, predation from coyotes and bears is most strongly influenced by the abiotic conditions that precede calving. However, when populations begin to decline, weather conditions during calving also influenced survival. I build on this analysis in Chapter 4 by determining the influence of climate change on the interplay between predation risk and neonate survival. I found that the relative equilibrium between bears and coyotes may not persist in the future as risk from coyotes could increase due to climate change. In Chapter 5, I investigate the relationships in niche overlap between caribou and their predators and how this may influence differential predation risk by affecting encounter rates. For coyotes, seasonal changes in niche overlap mirrored variation in caribou calf risk, but had less association with the rate of encounter with calves. In contrast, changes in niche overlap during the calving season for black bears had little association with these parameters. In Chapter 6, I examine broad-level habitat selection of caribou to study trade-offs between predator avoidance and foraging during the calving season. The results suggest that caribou movements are oriented towards increased access to foraging and the reduction of encounter risk with bears, and to a lesser extent, coyotes. Finally, I synthesize the major findings from this thesis and their relevance to caribou conservation in Chapter 7, to infer that Newfoundland caribou decline is ultimately driven by extrinsic and intrinsic elements related to density-dependence. Reduction in neonate survival emerged from nutritionally-stressed caribou females producing calves with lower survival.

Author Keywords: Behavioural ecology, Black bear (Ursus americanus), Coyote (Canis latrans), Population ecology, Predator-prey interactions, Woodland caribou (Rangifer tarandus)

Effective population size (Ne) is a key concept in population genetics, evolutionary biology and conservation biology that describes an important facet of genetic diversity and the capacity of populations to respond to future evolutionary pressures. The importance of Ne in management and conservation of wild populations encouraged the development of numerous genetic estimators which rely on a variety of methods. Despite the number and diversity of available Ne methods, however, tests of estimator performance have largely relied on simulations, with relatively few tests based on empirical data. I used well-studied wild populations of brook trout (Salvelinus fontinalis) in Algonquin Park, Ontario as a model system to assess the comparative performance of multiple Ne estimation methods and programs, comparing the resultant Ne estimates against demographic population size estimates.

As a first step, the genetic diversity and ancestry of wild brook trout populations was determined using 14 microsatellite loci. Genetic structure of brook trout populations showed variable contributions from historical supplemental stocking and also identified localized gene pools within and between watersheds, reflecting variable levels of connectivity and gene flow. Once the genetic ancestry and connectivity of populations had been resolved, single sample (point) and two samples (temporal) genetic estimators were used to estimate Ne of populations with pure native ancestry. Values obtained from genetic estimators utilizing both methods were variable within as well as among populations. Single sample (point) estimators were variable within individual populations, but substantially less than was observed among the temporal methods. The ratios of Ne to the estimated demographic population size (N) in small populations were substantially higher than in larger populations. Variation among estimates obtained from the different methods reflects varying assumptions that underlay the estimation algorithms. This research further investigated the effect of sampling effort and number of microsatellite loci used on Ne values obtained using the linkage disequilibrium (LD) estimation method. Ne estimates varied substantially among values generated from subsets of loci and genotyped individuals, highlighting the necessity for proper sampling design for efforts aiming to measure Ne. Despite the variation observed among and within estimation methods, the Ne concept is a valuable for the conservation and management of both exploited and endangered species.

Author Keywords: Brook Trout, Effective population size, Genetic Diversity, Genetic Structure

ABSTRACT

Calving site selection and fidelity in a restored elk (Cervus elaphus) herd in Bancroft

Ontario, Canada.

Michael R. Allan

Parturition site selection by ungulates is believed to be influenced by forage abundance and concealment from predators. In 2011 and 2012, I used vaginal implant transmitters and movements to identify calving sites for 23 GPS collared elk (Cervus elaphus) from a restored herd. I tested the hypothesis that maternal elk used sites with higher forage and denser concealment compared to pre-calving sites at micro and macrohabitat levels. I detected no significant microhabitat differences from direct measurements of vegetation. At the macrohabitat scale, based on proximity of landcover classes, mean distances to hardwood forests was significantly less for calving (153 m) than pre-calving sites (198 m).

Site fidelity is hypothesized to offer security in terms of familiarity to an area. I tested the hypothesis that females demonstrated fidelity to their previous year's location during pre-partum, parturition, post partum, breeding and winter periods. Elk were more philopatric during parturition and post partum than during breeding. Compared to winter elk were more philopatric during pre-partum, parturition and post-partum periods. Expressed as distance between consecutive-year calving locations, site fidelity varied with 27% of females exhibiting high (<1 km), 18% moderate and 55% (>2.9 km) low fidelity.

I measured nearest-neighbour distances at calving time, exploring the hypothesis that females distance themselves from conspecifics. Elk increased the average distances

to collared conspecifics during parturition; however, sample sizes were small. This strategy might influence calving site selection.

Rapid movement prior to parturition, low site fidelity and spacing-out of females during parturition appear to be strategies to minimize predator risk and detection. Little evidence of selection for vegetation structure suggests this may not be limiting to these elk.

Author Keywords: calving, elk, fidelity, movement, parturition, selection

Benthic macroinvertebrtes (BMI) are functionally important in aquatic ecosystems; as such, knowledge of their community structure and function is critical for understanding these systems. BMI were sampled from ten lakes in each of two regions of south-central Ontario to investigate which chemical and physical variables could be shaping their community structure and function. Ten Precambrian Shield lakes in the Muskoka-Haliburton region, and ten St Lawrence lowland lakes in the Kawartha lakes region were sampled. These lakes are geologically and chemically distinct, creating natural chemical and physical gradients within and between both regions. Community function was assessed using stable isotope analysis to elucidate carbon transfer dynamics (δ13C) and food web interactions (δ15N). It was predicted that the BMI from Shield lakes would have a δ13C signature indicative of allochthonous carbon subsidies, whereas the lowland lake BMI signatures would reflect autochthonous production. Additionally, it was predicted that the food web length (measured in δ15N units) would be different in Shield and lowland lakes. Both of these predictions were supported; however, the data indicate that δ13C signatures are more likely influenced by catchment geology (represented by bicarbonate concentration) than the extent of allochthony. The best predictor of food web length was found to be region. To assess BMI community structure, taxonomic richness, %EPT (% Ephemeroptera, Plecoptera, Trichoptera; a water quality index), and distribution of functional feeding groups were examined. Based on chemistry it was expected that the Shield lakes would be more speciose, and of greater water quality (relatively lower nutrient levels). These predictions were rejected; since there were no significant regional differences in taxonomic richness or biologically inferred water quality (%EPT). However, sediment size was found to best explain the variability in both metrics, with greater richness and %EPT found at sites with medium and small substrates than those with large substrates. Significant regional differences were found in the distribution of functional feeding groups. Most notably, there were significantly greater proportions of scrapers and shredders in the lowland and Shield lakes, respectively. Based on the feeding mechanisms of these invertebrates it can be inferred that allochthonous subsidies are likely of greater importance to Shield lake BMI communities than those of the lowland lakes; supporting the carbon transfer prediction. These findings provide insight about the structure and function of BMI communities from two dominant lake types in Ontario, and could be useful when determining how future chemical and physical changes will impact these communities.

Author Keywords: benthic macroinvertebrates, community function, community structure, Precambrian Shield, stable isotopes, St. Lawrence lowlands

Eastern wolves (Canis lycaon) have hybridized extensively with coyotes (C. latrans) and gray wolves (C. lupus) and are listed as a `species of special concern' in Canada. Previous studies have not linked genetic analysis with field data to investigate the mechanisms underlying Canis hybridization. Accordingly, I studied genetics, morphology, mortality, and behavior of wolves, coyotes, and hybrids in and adjacent to Algonquin Provincial Park (APP), Ontario. I documented 3 genetically distinct Canis types within the APP region that also differed morphologically, corresponding to putative gray wolves, eastern wolves, and coyotes. I also documented a substantial number of hybrids (36%) that exhibited intermediate morphology relative to parental types. I found that individuals with greater wolf ancestry occupied areas of higher moose density and fewer roads. Next, I studied intrinsic and extrinsic factors influencing survival and cause-specific mortality of canids in the hybrid zone. I found that survival was poor and harvest mortality was high for eastern wolves in areas adjacent to APP compared with other sympatric Canis types outside of APP and eastern wolves within APP. Contrary to previous studies of wolves and coyotes elsewhere, I hypothesized that all Canis types exhibit a high degree of spatial segregation in the Ontario hybrid zone. My hypothesis was supported as home range overlap and shared space use between neighboring Canis packs of all ancestry classes were low. Territoriality among Canis may increase the likelihood of eastern wolves joining coyote and hybrid packs and exacerbate hybridization. Canids outside APP modified their use of roads between night and day strongly at high road densities (selecting roads more at night), whereas they responded weakly at lower road densities (generally no selection). Individuals that survived exhibited a highly significant relationship between the difference in their night and day selection of roads and availability of roads, whereas those that died showed a weaker, non-significant response. My results suggest that canids in the unprotected landscape outside APP must balance trade-offs between exploiting benefits associated with secondary roads while mitigating risk of human-caused mortality. Overall, my results suggest that the distinct eastern wolf population of APP is unlikely to expand numerically and/or geographically under current environmental conditions and management regulations. If expansion of the APP eastern wolf population (numerically and in terms of its geographic distribution) is a conservation priority for Canada and Ontario, additional harvest protection in areas outside of APP may be required. If additional harvest protection is enacted, a detailed study within the new areas of protection would be important to document specific effects on eastern wolf population growth.

Author Keywords: Canis, coyotes, eastern wolves, hybridization, resource selection, survival

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

Coldwater streams are becoming increasingly impacted due to urbanization. Using environmental surveys, mark-recapture and telemetry, I assessed factors influencing seasonal brook trout (Salvelinus fontinalis) habitat use and movement in urban headwater streams in central Ontario between 2017-18. Generalized additive models were used to assess which habitat variables best explained seasonal yearling and older brook trout abundance, while generalized least squares models were used to assess overall trends in radio-tagged brook trout movement. My research demonstrated dynamic patterns in habitat use and movement by urban stream-dwelling brook trout. Yearlings were primarily influenced by water quality (stream temperature, conductivity), while older brook trout were most strongly influenced by stream morphology (depth, undercut bank). Movement occurred disproportionately around the spawning season and was more limited in the smaller, more altered stream. These findings may be used to inform fisheries managers on crucial timing and location of brook trout habitat refugia within urbanized environments.

Author Keywords: Brook trout, coldwater stream, groundwater, habitat use, radiotelemetry, urbanization

Invasive Phragmites australis ssp. australis (herein "Phragmites") has established and rapidly spread throughout many coastal areas of the Great Lakes. Known to displace native vegetation communities as it forms large, monotypic stands, Phragmites has a bad reputation when it comes to losses of biodiversity and habitat provision for wildlife. However, the extent to which Phragmites provides habitat for fish and invertebrates in coastal freshwater wetlands remains relatively unquantified. Thus, this study assessed whether fish assemblages and invertebrate communities in stands of Phragmites differ from those in stands of two native emergent vegetation communities, Typha spp. and Schoenoplectus spp. The findings showed significant differences in habitat variables among the vegetation communities in terms of water depth, macrophyte species richness, stem density and water quality. While abundance of the functional feeding group filterer-collectors was found to be significantly less in stands of Phragmites when compared to Schoenoplectus, no difference was observed in invertebrate taxa richness among vegetation communities. Lastly, no difference in fish assemblage or invertebrate community was detected when using multivariate analyses, implying that invasive Phragmites provides habitat that appears to be as valuable for fish and invertebrates as other emergent vegetation types in the St. Clair River Delta. The findings of this study will ultimately benefit the literature on invasive Phragmites and its role as fish habitat in freshwater wetlands, and aid management agencies in decisions regarding control of the invasive species.

Author Keywords: aquatic invasive species, aquatic macroinvertebrates, freshwater fish, freshwater wetlands, nMDS, Phragmites

Predation risk can elicit a range of responses in prey, but to date little is known about breadth of potential responses that may arise under realistic field conditions and how such responses are linked, leaving a fragmented picture of risk-related consequences on individuals. We increased predation risk in free-ranging snowshoe hares (Lepus americanus) during two consecutive summers by simulating natural chases using a model predator (i.e., domestic dog), and monitored hare stress physiology, energy expenditure, behaviour, condition, and habitat use. We show that higher levels of risk elicited marked changes in physiological stress metrics including sustained high levels of free plasma cortisol which had cascading effects on glucose, and immunology, but not condition. Risk-augmented hares also had lowered daily energy expenditure, spent more time foraging, and decreased rest, vigilance, and travel. It is possible that these alterations allowed risk-exposed hares to increase their condition at the same rate as controls. Additionally, risk-augmented hares selected, had high fidelity to, and were more mobile in structurally dense habitat (i.e., shrubs) which provided them additional cover from predators. They also used more open habitat (i.e., conifer) differently based on locale within the home range, using familiar conifer areas within cores for rest while moving through unfamiliar conifer areas in the periphery. Overall, these findings show that prey can have a multi-faceted, highly plastic response in the face of risk and can mitigate the effects of their stress physiology given the right environmental conditions.

Author Keywords: behaviour, condition, daily energy expenditure, predator-prey interactions, snowshoe hare, stress physiology

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