Aaron B. A. Shafer

Of glaciers and refugia: a decade of study sheds new light on the phylogeography of northwestern North America

Glacial cycles have played a dominant role in shaping the genetic structure and distribution of biota in northwestern North America. The two major ice age refugia of Beringia and the Pacific Northwest were connected by major mountain chains and bordered by the Pacific Ocean. As a result, numerous refugial options were available for the regions taxa during glacial advances. We reviewed the importance of glaciations and refugia in shaping northwestern North America’s phylogeographic history. We also tested whether ecological variables were associated with refugial history. The recurrent phylogeographic patterns that emerged were the following: (i) additional complexity, i.e. refugia within refugia, in both Beringia and the Pacific Northwest; and (ii) strong evidence for cryptic refugia in the Alexander Archipelago and Haida Gwaii, the Canadian Arctic and within the ice‐sheets. Species with contemporary ranges that covered multiple refugia, or those with high dispersal ability, were significantly more likely to have resided in multiple refugia. Most of the shared phylogeographic patterns can be attributed to multiple refugial locales during the last glacial maximum or major physiographic barriers like rivers and glaciers. However, some of the observed patterns are much older and appear connected to the orogeny of the Cascade‐Sierra chain or allopatric differentiation during historic glacial advances. The emergent patterns from this review suggest we should refine the classic Beringian‐southern refugial paradigm for northwestern North American biota and highlight the ecological and evolutionary consequences of colonization from multiple refugia.

HOT SPOTS OF GENETIC DIVERSITY DESCENDED FROM MULTIPLE PLEISTOCENE REFUGIA IN AN ALPINE UNGULATE

Species that inhabit naturally fragmented environments are expected to be spatially structured and exhibit reduced genetic diversity at the periphery of their range. Patterns of differentiation may also reflect historical processes such as recolonization from glacial refugia. We examined the relative importance of these factors in shaping the spatial patterns of genetic differentiation across the range of an alpine specialist, the North American mountain goat (Oreamnos americanus). Contrary to fossil evidence that suggests a single southern refugium, we detected evidence for additional refugia in northern British Columbia and the Alaskan coast using both mitochondrial and microsatellite DNA. A core area of elevated genetic diversity characterized both regions, and molecular dating suggested a recent Pleistocene split was followed by demographic expansion. Across their range, mountain goats were highly genetically structured and displayed the expected pattern of declining diversity toward the periphery. Gene flow was high within contiguous mountain ranges, but cross‐assignments paradoxically suggest that long‐distance contemporary dispersal movements are not uncommon. These results improve our understanding of how historical vicariance and contemporary fragmentation influence population differentiation, and have implications for conserving the adaptive potential of alpine populations and habitat.

Habitat selection predicts genetic relatedness in an alpine ungulate

Landscape heterogeneity plays an integral role in shaping ecological and evolutionary processes. Despite links between the two disciplines, ecologists and population geneticists have taken different approaches to evaluating habitat selection, animal movement, and gene flow across the landscape. Ecologists commonly use statistical models such as resource selection functions (RSFs) to identify habitat features disproportionately selected by animals, whereas population genetic approaches model genetic differentiation according to the distribution of habitat variables. We combined ecological and genetic approaches by using RSFs to predict genetic relatedness across a heterogeneous landscape. We constructed sex- and season-specific resistance surfaces based on RSFs estimated using data from 102 GPS (global positioning system) radio-collared mountain goats (Oreamnos americanus) in southeast Alaska, USA. Based on mountain goat ecology, we hypothesized that summer and male surfaces would be the best predictors of relatedness. All individuals were genotyped at 22 microsatellite loci, which we used to estimate genetic relatedness. Summer resistance surfaces derived from RSFs were the best predictors of genetic relatedness, and winter models the poorest. Mountain goats generally selected for areas close to escape terrain and with a high heat load (a metric related to vegetative productivity and snow depth), while avoiding valleys. Male- and female-specific surfaces were similar, except for winter, for which male habitat selection better predicted genetic relatedness. The null models of isolation-by-distance and barrier only outperformed the winter models. This study merges high-resolution individual locations through GPS telemetry and genetic data, that can be used to validate and parameterize landscape genetics models, and further elucidates the relationship between landscape heterogeneity and genetic differentiation.

Cophylogeny of Nosema (Microsporidia: Nosematidae) and Bees (Hymenoptera: Apidae) Suggests Both Cospeciation and a Host-switch

Abstract Some microsporidian parasites belonging to the genus Nosema infect bees. Previous phylogenies of these parasites have produced alternative, conflicting relationships. We analyzed separately, and in combination, large and small subunit ribosomal DNA sequences of Nosema species infecting bees under neighbor-joining, maximum parsimony, maximum likelihood, and Bayesian frameworks. We observed a sister relationship between Nosema ceranae and Nosema bombi, with Nosema apis as a basal member to this group. When compared to their respective hosts (Apis cerana, Bombus spp., and A. mellifera), 2 plausible evolutionary scenarios emerged. The first hypothesis involves a common ancestor of N. bombi host-switching from a historical Bombus lineage to A. cerana. The second suggests an ancestral N. ceranae host-switching to a species of Bombus. The reported events offer insight into the evolutionary history of these organisms and may explain host specificity and virulence of Nosema in these economically important insects.

Temporal dynamics of genetic variability in a mountain goat (Oreamnos americanus) population

The association between population dynamics and genetic variability is of fundamental importance for both evolutionary and conservation biology. We combined long‐term population monitoring and molecular genetic data from 123 offspring and their parents at 28 microsatellite loci to investigate changes in genetic diversity over 14 cohorts in a small and relatively isolated population of mountain goats (Oreamnos americanus) during a period of demographic increase. Offspring heterozygosity decreased while parental genetic similarity and inbreeding coefficients (FIS) increased over the study period (1995–2008). Immigrants introduced three novel alleles into the population and matings between residents and immigrants produced more heterozygous offspring than local crosses, suggesting that immigration can increase population genetic variability. The population experienced genetic drift over the study period, reflected by a reduced allelic richness over time and an ‘isolation‐by‐time’ pattern of genetic structure. The temporal decline of individual genetic diversity despite increasing population size probably resulted from a combination of genetic drift due to small effective population size, inbreeding and insufficient counterbalancing by immigration. This study highlights the importance of long‐term genetic monitoring to understand how demographic processes influence temporal changes of genetic diversity in long‐lived organisms.

Demographic inferences using short-read genomic data in an approximate Bayesian computation framework: in silico evaluation of power, biases and proof of concept in Atlantic walrus

Approximate Bayesian computation (ABC) is a powerful tool for model‐based inference of demographic histories from large genetic data sets. For most organisms, its implementation has been hampered by the lack of sufficient genetic data. Genotyping‐by‐sequencing (GBS) provides cheap genome‐scale data to fill this gap, but its potential has not fully been exploited. Here, we explored power, precision and biases of a coalescent‐based ABC approach where GBS data were modelled with either a population mutation parameter (θ) or a fixed site (FS) approach, allowing single or several segregating sites per locus. With simulated data ranging from 500 to 50 000 loci, a variety of demographic models could be reliably inferred across a range of timescales and migration scenarios. Posterior estimates were informative with 1000 loci for migration and split time in simple population divergence models. In more complex models, posterior distributions were wide and almost reverted to the uninformative prior even with 50 000 loci. ABC parameter estimates, however, were generally more accurate than an alternative composite‐likelihood method. Bottleneck scenarios proved particularly difficult, and only recent bottlenecks without recovery could be reliably detected and dated. Notably, minor‐allele‐frequency filters – usual practice for GBS data – negatively affected nearly all estimates. With this in mind, we used a combination of FS and θ approaches on empirical GBS data generated from the Atlantic walrus (Odobenus rosmarus rosmarus), collectively providing support for a population split before the last glacial maximum followed by asymmetrical migration and a high Arctic bottleneck. Overall, this study evaluates the potential and limitations of GBS data in an ABC‐coalescence framework and proposes a best‐practice approach.

Does Learning or Instinct Shape Habitat Selection

Habitat selection is an important behavioural process widely studied for its population-level effects. Models of habitat selection are, however, often fit without a mechanistic consideration. Here, we investigated whether patterns in habitat selection result from instinct or learning for a population of grizzly bears (Ursus arctos) in Alberta, Canada. We found that habitat selection and relatedness were positively correlated in female bears during the fall season, with a trend in the spring, but not during any season for males. This suggests that habitat selection is a learned behaviour because males do not participate in parental care: a genetically predetermined behaviour (instinct) would have resulted in habitat selection and relatedness correlations for both sexes. Geographic distance and home range overlap among animals did not alter correlations indicating that dispersal and spatial autocorrelation had little effect on the observed trends. These results suggest that habitat selection in grizzly bears are partly learned from their mothers, which could have implications for the translocation of wildlife to novel environments.

Genome-wide cross-amplification of domestic sheep microsatellites in bighorn sheep and mountain goats.

We tested for cross‐species amplification of microsatellite loci located throughout the domestic sheep (Ovis aries) genome in two north American mountain ungulates (bighorn sheep, Ovis canadensis, and mountain goats, Oreamnos americanus). We identified 247 new polymorphic markers in bighorn sheep (≥ 3 alleles in one of two study populations) and 149 in mountain goats (≥ 2 alleles in a single study population) using 648 and 576 primer pairs, respectively. Our efforts increased the number of available polymorphic microsatellite markers to 327 for bighorn sheep and 180 for mountain goats. The average distance between successive polymorphic bighorn sheep and mountain goat markers inferred from the Australian domestic sheep genome linkage map (mean ± 1 SD) was 11.9 ± 9.2 and 15.8 ± 13.8 centimorgans, respectively. The development of genomic resources in these wildlife species enables future studies of the genetic architecture of trait variation.

Lack of) Genetic Diversity in Immune Genes Predates Glacial Isolation in the North American Mountain Goat (Oreamnos americanus)

The major histocompatibility complex (MHC) plays an important role in an organisms ability to respond to pathogens. Immunogenetic diversity is advantageous as it permits the recognition of more external antigens. For this reason, MHC and immune gene variation are considered a barometer for the genetic health of wild populations. Mountain goats (Oreamnos americanus) were previously shown to have little variation at the MHC Class II Oram-DRB locus, which was attributed to population bottlenecks during the last glacial maximum (LGM). In this paper, we extended the analysis of immunogenetic variability in mountain goats to 5 genes representing the 3 classes of MHC gene (Class I OLA, Class II DRA and DRB, and Class III TNF-α) and the natural resistance-associated macrophage protein. We sequenced approximately 3000 bp from 31 individuals sampled across the range of mountain goats and found very low levels of diversity (1-3 polymorphic sites per gene) with the exception of the Class I Oram-OLA gene. Oram-OLA was nearly 30 times more diverse than the other immune genes and appears to represent a source of increased immunogenetic diversity. This diversity may be attributed to multiple loci, mediated by pathogen exposure, or potentially influenced by social factors. The distribution of SNPs was not associated with refugial history, suggesting that the current distribution of immunogenetic diversity was present prior to the LGM. These data suggest that although they have low levels of diversity at the 4 of 5 immune loci, mountain goats may be better equipped for future climate oscillations and pathogen exposure than previously thought.

LetterReply to Garner et al.

Genomics in Conservation : Case Studies and Bridging the Gap between Data and Application Reply