Mixed population genomics support for the central marginal hypothesis across the invasive range of the cane toad (Rhinella marina) in Australia

Understanding factors that cause species' geographic range limits is a major focus in ecology and evolution. The central marginal hypothesis (CMH) predicts that species cannot adapt to conditions beyond current geographic range edges because genetic diversity decreases from core to edge due to smaller, more isolated edge populations. We employed a population genomics framework using 24 235–33 112 SNP loci to test major predictions of the CMH in the ongoing invasion of the cane toad (Rhinella marina) in Australia. Cane toad tissue samples were collected along broad‐scale, core‐to‐edge transects across their invasive range. Geographic and ecological core areas were identified using GIS and habitat suitability indices from ecological niche modelling. Bayesian clustering analyses revealed three genetic clusters, in the northwest invasion‐front region, northeast precipitation‐limited region and southeast cold temperature‐limited region. Core‐to‐edge patterns of genetic diversity and differentiation were consistent with the CMH in the southeast, but were not supported in the northeast and showed mixed support in the northwest. Results suggest cold temperatures are a likely contributor to southeastern range limits, consistent with CMH predictions. In the northeast and northwest, ecological processes consisting of a steep physiological barrier and ongoing invasion dynamics, respectively, are more likely explanations for population genomic patterns than the CMH.     

Disentangling the effects of geographic peripherality and habitat suitability on neutral and adaptive genetic variation in Swiss stone pine

It is generally accepted that the spatial distribution of neutral genetic diversity within a species’ native range mostly depends on effective population size, demographic history, and geographic position. However, it is unclear how genetic diversity at adaptive loci correlates with geographic peripherality or with habitat suitability within the ecological niche. Using exome‐wide genomic data and distribution maps of the Alpine range, we first tested whether geographic peripherality correlates with four measures of population genetic diversity at > 17,000 SNP loci in 24 Alpine populations (480 individuals) of Swiss stone pine (Pinus cembra) from Switzerland. To distinguish between neutral and adaptive SNP sets, we used four approaches (two gene diversity estimates, F_ST outlier test, and environmental association analysis) that search for signatures of selection. Second, we established ecological niche models for P. cembra in the study range and investigated how habitat suitability correlates with genetic diversity at neutral and adaptive loci. All estimates of neutral genetic diversity decreased with geographic peripherality, but were uncorrelated with habitat suitability. However, heterozygosity (H_e) at adaptive loci based on Tajima's D declined significantly with increasingly suitable conditions. No other diversity estimates at adaptive loci were correlated with habitat suitability. Our findings suggest that populations at the edge of a species' geographic distribution harbour limited neutral genetic diversity due to demographic properties. Moreover, we argue that populations from suitable habitats went through strong selection processes, are thus well adapted to local conditions, and therefore exhibit reduced genetic diversity at adaptive loci compared to populations at niche margins.     

Spatial genetic structure in Beta vulgaris subsp. maritima and Beta macrocarpa reveals the effect of contrasting mating system,influence of marine currents,and footprints of postglacial recolonization routes

Understanding the factors that contribute to population genetic divergence across a species' range is a long‐standing goal in evolutionary biology and ecological genetics. We examined the relative importance of historical and ecological features in shaping the present‐day spatial patterns of genetic structure in two related plant species, Beta vulgaris subsp. maritima and Beta macrocarpa. Using nuclear and mitochondrial markers, we surveyed 93 populations from Brittany (France) to Morocco – the southern limit of their species' range distribution. Whereas B. macrocarpa showed a genotypic structure and a high level of genetic differentiation indicative of selfing, the population genetic structure of B. vulgaris subsp. maritima was consistent with an outcrossing mating system. We further showed (1) a strong geographic clustering in coastal B. vulgaris subsp. maritima populations that highlighted the influence of marine currents in shaping different lineages and (2) a peculiar genetic structure of inland B. vulgaris subsp. maritima populations that could indicate the admixture of distinct evolutionary lineages and recent expansions associated with anthropogenic disturbances. Spatial patterns of nuclear diversity and differentiation also supported a stepwise recolonization of Europe from Atlantic‐Mediterranean refugia after the last glacial period, with leading‐edge expansions. However, cytoplasmic diversity was not impacted by postglacial recolonization: stochastic long‐distance seed dispersal mediated by major oceanic currents may mitigate the common patterns of reduced cytoplasmic diversity observed for edge populations. Overall, the patterns we documented here challenge the general view of reduced genetic diversity at the edge of a species' range distribution and provide clues for understanding how life‐history and major geographic features interact to shape the distribution of genetic diversity.     

Losing your edge: climate change and the conservation value of range‐edge populations

Populations occurring at species' range edges can be locally adapted to unique environmental conditions. From a species' perspective, range‐edge environments generally have higher severity and frequency of extreme climatic events relative to the range core. Under future climates, extreme climatic events are predicted to become increasingly important in defining species' distributions. Therefore, range‐edge genotypes that are better adapted to extreme climates relative to core populations may be essential to species' persistence during periods of rapid climate change. We use relatively simple conceptual models to highlight the importance of locally adapted range‐edge populations (leading and trailing edges) for determining the ability of species to persist under future climates. Using trees as an example, we show how locally adapted populations at species' range edges may expand under future climate change and become more common relative to range‐core populations. We also highlight how large‐scale habitat destruction occurring in some geographic areas where many species range edge converge, such as biome boundaries and ecotones (e.g., the arc of deforestation along the rainforest‐cerrado ecotone in the southern Amazonia), can have major implications for global biodiversity. As climate changes, range‐edge populations will play key roles in helping species to maintain or expand their geographic distributions. The loss of these locally adapted range‐edge populations through anthropogenic disturbance is therefore hypothesized to reduce the ability of species to persist in the face of rapid future climate change.     

Ecological niche models display nonlinear relationships with abundance and demographic performance across the latitudinal distribution of Astragalus utahensis (Fabaceae)

The potential for ecological niche models (ENMs) to accurately predict species' abundance and demographic performance throughout their geographic distributions remains a topic of substantial debate in ecology and biogeography. Few studies simultaneously examine the relationship between ENM predictions of environmental suitability and both a species' abundance and its demographic performance, particularly across its entire geographic distribution. Yet, studies of this type are essential for understanding the extent to which ENMs are a viable tool for identifying areas that may promote high abundance or performance of a species or how species might respond to future climate conditions. In this study, we used an ensemble ecological niche model to predict climatic suitability for the perennial forb Astragalus utahensis across its geographic distribution. We then examined relationships between projected climatic suitability and field‐based measures of abundance, demographic performance, and forecasted stochastic population growth (λ_s). Predicted climatic suitability showed a J‐shaped relationship with A. utahensis abundance, where low‐abundance populations were associated with low‐to‐intermediate suitability scores and abundance increased sharply in areas of high predicted climatic suitability. A similar relationship existed between climatic suitability and λ_s from the center to the northern edge of the latitudinal distribution. Patterns such as these, where density or demographic performance only increases appreciably beyond some threshold of climatic suitability, support the contention that ENM‐predicted climatic suitability does not necessarily represent a reliable predictor of abundance or performance across large geographic regions.     

Waterscape genetics of the yellow perch ( Perca flavescens): patterns across large connected ecosystems and isolated relict populations

Comparisons of a species' genetic diversity and divergence patterns across large connected populations vs. isolated relict areas provide important data for understanding potential response to global warming, habitat alterations and other perturbations. Aquatic taxa offer ideal case studies for interpreting these patterns, because their dispersal and gene flow often are constrained through narrow connectivity channels that have changed over geological time and/or from contemporary anthropogenic perturbations. Our research objective is to better understand the interplay between historic influences and modern‐day factors (fishery exploitation, stocking supplementation and habitat loss) in shaping population genetic patterns of the yellow perch Perca flavescens (Percidae: Teleostei) across its native North American range. We employ a modified landscape genetics approach, analysing sequences from the entire mitochondrial DNA control region and 15 nuclear DNA microsatellite loci of 664 spawning adults from 24 populations. Results support that perch from primary glacial refugium areas (Missourian, Mississippian and Atlantic) founded contemporary northern populations. Genetic diversity today is highest in southern (never glaciated) populations and also is appreciable in northern areas that were founded from multiple refugia. Divergence is greater among isolated populations, both north and south; the southern Gulf Coast relict populations are the most divergent, reflecting their long history of isolation. Understanding the influence of past and current waterway connections on the genetic structure of yellow perch populations may help us to assess the roles of ongoing climate change and habitat disruptions towards conserving aquatic biodiversity.     

History matters more when explaining genetic diversity within the context of the core–periphery hypothesis

The core–periphery hypothesis (CPH) predicts that populations located at the periphery of a species' range should have lower levels of genetic variation than those at the centre of the range. However, most of the research on the CPH focuses on geographic distance and not on ecological distance, or uses categorical definitions of core and periphery to explain the distribution of genetic diversity. We use current climate data and historical climate data from the last glacial maxima to develop quantitative estimates of contemporary and historical ecological suitability using ecological niche models. We analysed genetic diversity using 12 polymorphic microsatellites to estimate changes in heterozygosity, allelic richness and population differentiation in 31 populations of the wood frog (Lithobates sylvaticus) spanning the species’ entire eastern clade (33^o to 45^o latitude) from Alabama, USA, to Nova Scotia, Canada. Our data support predictions based on the CPH. Populations showed significant differences in genetic diversity across the range, with lower levels of genetic variation at the geographic range edge and in areas with lower levels of historical and contemporary ecological suitability. However, history and geography (not current ecological suitability) best explain the patterns. This study highlights the importance of examining more than just geography when assessing the CPH, and the importance of historical ecological suitability in the maintenance of genetic diversity and population differentiation.     

An approach for identifying cryptic barriers to gene flow that limit species' geographic ranges

Species' geographic range limits are most often not demarcated by obvious dispersal barriers. Poor‐quality habitat at the edge of a species' range can prevent range expansion by preventing outward migration or through reducing adaptive potential resulting from decreased genetic diversity. We identified habitat variables that constrain gene flow across the entire geographic range of an endemic salamander (Ambystoma barbouri) in the eastern United States, and we tested whether increased resistance resulting from these variables provides cryptic dispersal barriers at the range edges. Using polymorphic microsatellite loci, we first identified three genetic clusters that are separated by the Ohio and Kentucky rivers. Through a combination of landscape genetic analyses and generalized dissimilarity modelling, we then classified variables that (i) restrict gene flow in each of the genetic clusters across the geographic distribution of A. barbouri and (ii) become more common towards the peripheries of the distribution. A decrease in limestone availability and an increase in growing season precipitation were correlated with high resistance to gene flow across the range, and both became more common at the edges of the species' distribution. However, other landscape variables were more important for explaining variation in geneflow rates in different portions of the range, such as increased mean annual temperature and frost‐free period in the south vs. growing season precipitation in the north. Taken together, these results suggest that there are both range‐wide and regionally specific cryptic habitat barriers preventing geographic range expansion. Species ‘geographic range limits are probably governed by a set of ecological and evolutionary factors, and our landscape genetic approach could be applied to gain additional insight into many systems.     

Rangewide landscape genetics of an endemic Pacific northwestern salamander

A species' genetic structure often varies in response to ecological and landscape processes that differ throughout the species' geographic range, yet landscape genetics studies are rarely spatially replicated. The Cope's giant salamander (Dicamptodon copei) is a neotenic, dispersal‐limited amphibian with a restricted geographic range in the Pacific northwestern USA. We investigated which landscape factors affect D. copei gene flow in three regions spanning the species' range, which vary in climate, landcover and degree of anthropogenic disturbance. Least cost paths and Circuitscape resistance analyses revealed that gene flow patterns vary across the species' range, with unique combinations of landscape variables affecting gene flow in different regions. Populations in the northern coastal portions of the range had relatively high gene flow, largely facilitated by stream and river networks. Near the southeastern edge of the species' range, gene flow was more restricted overall, with relatively less facilitation by streams and more limitation by heat load index and fragmented forest cover. These results suggested that the landscape is more difficult for individuals to disperse through at the southeastern edge of the species' range, with terrestrial habitat desiccation factors becoming more limiting to gene flow. We suggest that caution be used when attempting to extrapolate landscape genetic models and conservation measures from one portion of a species' range to another.     

Habitat fragmentation in coastal southern California disrupts genetic connectivity in the cactus wren (Campylorhynchus brunneicapillus)

Achieving long‐term persistence of species in urbanized landscapes requires characterizing population genetic structure to understand and manage the effects of anthropogenic disturbance on connectivity. Urbanization over the past century in coastal southern California has caused both precipitous loss of coastal sage scrub habitat and declines in populations of the cactus wren (Campylorhynchus brunneicapillus). Using 22 microsatellite loci, we found that remnant cactus wren aggregations in coastal southern California comprised 20 populations based on strict exact tests for population differentiation, and 12 genetic clusters with hierarchical Bayesian clustering analyses. Genetic structure patterns largely mirrored underlying habitat availability, with cluster and population boundaries coinciding with fragmentation caused primarily by urbanization. Using a habitat model we developed, we detected stronger associations between habitat‐based distances and genetic distances than Euclidean geographic distance. Within populations, we detected a positive association between available local habitat and allelic richness and a negative association with relatedness. Isolation‐by‐distance patterns varied over the study area, which we attribute to temporal differences in anthropogenic landscape development. We also found that genetic bottleneck signals were associated with wildfire frequency. These results indicate that habitat fragmentation and alterations have reduced genetic connectivity and diversity of cactus wren populations in coastal southern California. Management efforts focused on improving connectivity among remaining populations may help to ensure population persistence.     

Seascape drivers of Macrocystis pyrifera population genetic structure in the northeast Pacific

At small spatial and temporal scales, genetic differentiation is largely controlled by constraints on gene flow, while genetic diversity across a species' distribution is shaped on longer temporal and spatial scales. We assess the hypothesis that oceanographic transport and other seascape features explain different scales of genetic structure of giant kelp, Macrocystis pyrifera. We followed a hierarchical approach to perform a microsatellite‐based analysis of genetic differentiation in Macrocystis across its distribution in the northeast Pacific. We used seascape genetic approaches to identify large‐scale biogeographic population clusters and investigate whether they could be explained by oceanographic transport and other environmental drivers. We then modelled population genetic differentiation within clusters as a function of oceanographic transport and other environmental factors. Five geographic clusters were identified: Alaska/Canada, central California, continental Santa Barbara, California Channel Islands and mainland southern California/Baja California peninsula. The strongest break occurred between central and southern California, with mainland Santa Barbara sites forming a transition zone between the two. Breaks between clusters corresponded approximately to previously identified biogeographic breaks, but were not solely explained by oceanographic transport. An isolation‐by‐environment (IBE) pattern was observed where the northern and southern Channel Islands clustered together, but not with closer mainland sites, despite the greater distance between them. The strongest environmental association with this IBE pattern was observed with light extinction coefficient, which extends suitable habitat to deeper areas. Within clusters, we found support for previous results showing that oceanographic connectivity plays an important role in the population genetic structure of Macrocystis in the Northern hemisphere.     

Predicting the effects of climate change on population connectivity and genetic diversity of an imperiled freshwater mussel,Cumberlandia monodonta (Bivalvia: Margaritiferidae), in riverine systems

In the face of global climate change, organisms may respond to temperature increases by shifting their ranges poleward or to higher altitudes. However, the direction of range shifts in riverine systems is less clear. Because rivers are dendritic networks, there is only one dispersal route from any given location to another. Thus, range shifts are only possible if branches are connected by suitable habitat, and stream‐dwelling organisms can disperse through these branches. We used Cumberlandia monodonta (Bivalvia: Unionoida: Margaritiferidae) as a model species to investigate the effects of climate change on population connectivity because a majority of contemporary populations are panmictic. We combined ecological niche models (ENMs) with population genetic simulations to investigate the effects of climate change on population connectivity and genetic diversity of C. monodonta. The ENMs were constructed using bioclimatic and landscape data to project shifts in suitable habitat under future climate scenarios. We then used forward‐time simulations to project potential changes in genetic diversity and population connectivity based on these range shifts. ENM results under current conditions indicated long stretches of highly suitable habitat in rivers where C. monodonta persists; populations in the upper Mississippi River remain connected by suitable habitat that does not impede gene flow. Future climate scenarios projected northward and headwater‐ward range contraction and drastic declines in habitat suitability for most extant populations throughout the Mississippi River Basin. Simulations indicated that climate change would greatly reduce genetic diversity and connectivity across populations. Results suggest that a single, large population of C. monodonta will become further fragmented into smaller populations, each of which will be isolated and begin to differentiate genetically. Because C. monodonta is a widely distributed species and purely aquatic, our results suggest that persistence and connectivity of stream‐dwelling organisms will be significantly altered in response to future climate change.     

Spatial differences in genetic diversity and northward migration suggest genetic erosion along the boreal caribou southern range limit and continued range retraction

With increasing human activities and associated landscape changes, distributions of terrestrial mammals become fragmented. These changes in distribution are often associated with reduced population sizes and loss of genetic connectivity and diversity (i.e., genetic erosion) which may further diminish a species' ability to respond to changing environmental conditions and lead to local population extinctions. We studied threatened boreal caribou (Rangifer tarandus caribou) populations across their distribution in Ontario/Manitoba (Canada) to assess changes in genetic diversity and connectivity in areas of high and low anthropogenic activity. Using data from >1,000 caribou and nine microsatellite loci, we assessed population genetic structure, genetic diversity, and recent migration rates using a combination of network and population genetic analyses. We used Bayesian clustering analyses to identify population genetic structure and explored spatial and temporal variation in those patterns by assembling networks based on R_ST and F_ST as historical and contemporary genetic edge distances, respectively. The Bayesian clustering analyses identified broad‐scale patterns of genetic structure and closely aligned with the R_ST network. The F_ST network revealed substantial contemporary genetic differentiation, particularly in areas presenting contemporary anthropogenic disturbances and habitat fragmentation. In general, relatively lower genetic diversity and greater genetic differentiation were detected along the southern range limit, differing from areas in the northern parts of the distribution. Moreover, estimation of migration rates suggested a northward movement of animals away from the southern range limit. The patterns of genetic erosion revealed in our study suggest ongoing range retraction of boreal caribou in central Canada.     

High genetic diversity declines towards the geographic range periphery of Adonis vernalis,a Eurasian dry grassland plant

Genetic diversity is important for species' fitness and evolutionary processes but our knowledge on how it varies across a species' distribution range is limited. The abundant centre hypothesis (ACH) predicts that populations become smaller and more isolated towards the geographic range periphery – a pattern that in turn should be associated with decreasing genetic diversity and increasing genetic differentiation. We tested this hypothesis in Adonis vernalis, a dry grassland plant with an extensive Eurasian distribution. Its life‐history traits and distribution characteristics suggest a low genetic diversity that decreases and a high genetic differentiation that increases towards the range edge. We analysed AFLP fingerprints in 28 populations along a 4698‐km transect from the geographic range core in Russia to the western range periphery in Central and Western Europe. Contrary to our expectation, our analysis revealed high genetic diversity (range of proportion of polymorphic bands = 56–81%, He = 0.168–0.238) and low genetic differentiation across populations (Φ_ST = 0.18). However, in congruence with the genetic predictions of the ACH, genetic diversity decreased and genetic differentiation increased towards the range periphery. Spanish populations were genetically distinct, suggesting a divergent post‐glacial history in this region. The high genetic diversity and low genetic differentiation in the remaining A. vernalis populations is surprising given the species' life‐history traits and points to the possibility that the species has been widely distributed in the studied region or that it has migrated from a diverse source in an East–West direction, in the past.     

Population dynamics can be more important than physiological limits for determining range shifts under climate change

Evidence is accumulating that species' responses to climate changes are best predicted by modelling the interaction of physiological limits, biotic processes and the effects of dispersal‐limitation. Using commercially harvested blacklip (Haliotis rubra) and greenlip abalone (Haliotis laevigata) as case studies, we determine the relative importance of accounting for interactions among physiology, metapopulation dynamics and exploitation in predictions of range (geographical occupancy) and abundance (spatially explicit density) under various climate change scenarios. Traditional correlative ecological niche models (ENM) predict that climate change will benefit the commercial exploitation of abalone by promoting increased abundances without any reduction in range size. However, models that account simultaneously for demographic processes and physiological responses to climate‐related factors result in future (and present) estimates of area of occupancy (AOO) and abundance that differ from those generated by ENMs alone. Range expansion and population growth are unlikely for blacklip abalone because of important interactions between climate‐dependent mortality and metapopulation processes; in contrast, greenlip abalone should increase in abundance despite a contraction in AOO. The strongly non‐linear relationship between abalone population size and AOO has important ramifications for the use of ENM predictions that rely on metrics describing change in habitat area as proxies for extinction risk. These results show that predicting species' responses to climate change often require physiological information to understand climatic range determinants, and a metapopulation model that can make full use of this data to more realistically account for processes such as local extirpation, demographic rescue, source‐sink dynamics and dispersal‐limitation.     

Demographic and population‐genetic tests provide mixed support for the abundant centre hypothesis in the endemic plant Leavenworthia stylosa

The abundant centre hypothesis (ACH) assumes that population abundance, population size, density and per‐capita reproductive output should peak at the centre of a species' geographic range and decline towards the periphery. Increased isolation among and decreased reproductive output within edge populations should reduce within‐population genetic diversity and increase genetic differentiation among edge relative to central populations. The ACH also predicts asymmetrical gene flow, with net movement of migrants from the centre to edges. We evaluated these ecological assumptions and population‐genetic predictions in the endemic flowering plant Leavenworthia stylosa. Although populations were more spatially isolated near range edges, the geographic centre was surrounded by and not coincident with areas of peak population abundance, and plant density increased towards range edges. Per‐capita seed number was not associated with distance to the range centre, but seed number/m² increased near range edges. In support of ACH predictions, allelic diversity at 12 microsatellite loci declined with distance from the range centre, and pairwise F_ST values were higher between edge populations than between central populations. Coalescent analyses confirmed that gene flow was most infrequent between edge populations, but there was not an asymmetric pattern of gene flow predicted by the ACH. This study shows that among‐population demographic variability largely did not support the ACH, while patterns of genetic diversity, differentiation and gene flow were generally consistent with its predictions. Such mixed support has frequently been observed in tests of the ACH and raises concerns regarding the generality of this hypothesis for species range limits.     

Postglacial recolonization in a cold climate specialist in western Europe: patterns of genetic diversity in the adder (Vipera berus) support the central–marginal hypothesis

Understanding the impact of postglacial recolonization on genetic diversity is essential in explaining current patterns of genetic variation. The central–marginal hypothesis (CMH) predicts a reduction in genetic diversity from the core of the distribution to peripheral populations, as well as reduced connectivity between peripheral populations. While the CMH has received considerable empirical support, its broad applicability is still debated and alternative hypotheses predict different spatial patterns of genetic diversity. Using microsatellite markers, we analysed the genetic diversity of the adder (Vipera berus) in western Europe to reconstruct postglacial recolonization. Approximate Bayesian Computation (ABC) analyses suggested a postglacial recolonization from two routes: a western route from the Atlantic Coast up to Belgium and a central route from the Massif Central to the Alps. This cold‐adapted species likely used two isolated glacial refugia in southern France, in permafrost‐free areas during the last glacial maximum. Adder populations further from putative glacial refugia had lower genetic diversity and reduced connectivity; therefore, our results support the predictions of the CMH. Our study also illustrates the utility of highly variable nuclear markers, such as microsatellites, and ABC to test competing recolonization hypotheses.     

Distance,elevation and environment as drivers of diversity and divergence in bumble bees across latitude and altitude

Identifying drivers of dispersal limitation and genetic differentiation is a key goal in biogeography. We examine patterns of population connectivity and genetic diversity using restriction site‐associated DNA sequencing (RADseq) in two bumble bee species, Bombus vosnesenskii and Bombus bifarius, across latitude and altitude in mountain ranges from California, Oregon and Washington, U.S.A. Bombus vosnesenskii, which occurs across a broader elevational range at most latitudes, exhibits little population structure while B. bifarius, which occupies a relatively narrow higher elevation niche across most latitudes, exhibits much stronger population differentiation, although gene flow in both species is best explained by isolation with environmental niche resistance. A relationship between elevational habitat breadth and genetic diversity is also apparent, with B. vosnesenskii exhibiting relatively consistent levels of genetic diversity across its range, while B. bifarius has reduced genetic diversity at low latitudes, where it is restricted to high‐elevation habitat. The results of this study highlight the importance of the intersect between elevational range and habitat suitability in influencing population connectivity and suggest that future climate warming will have a fragmenting effect even on populations that are presently well connected, as they track their thermal niches upward in montane systems.     

Climate change reduces genetic diversity of Canada lynx at the trailing range edge

Shifts in species distributions due to environmental change may affect the spatial pattern of genetic structure within a species' range, including possible changes to the adaptive potential of populations. We investigated spatial patterns of neutral genetic diversity and differentiation at the southern edge of the Canada lynx Lynx canadensis distribution in Ontario, Canada. We analyzed provincial fur harvest records (1972–2010) and collected and genotyped lynx pelt samples (2007–2009) from 702 lynx at 14 microsatellite loci. We show that the southern range boundary of lynx in central Canada has contracted northward by > 175 km since the 1970s, and that high winter temperature, low snow depth, and low proportion of suitable habitat are strongly correlated with low neutral genetic diversity and high genetic differentiation at the trailing range edge. Our work tests fundamental ideas about species range limits and demonstrates that environmental conditions can have a marked influence on neutral genetic structure. Our results suggest that changes in environmental conditions will result in further loss of genetic diversity and possibly reduce adaptive potential in southern peripheral lynx populations.     

Colonization from divergent ancestors: glaciation signatures on contemporary patterns of genomic variation in Collared Pikas (Ochotona collaris)

Identifying the genetic structure of a species and the factors that drive it is an important first step in modern population management, in part because populations evolving from separate ancestral sources may possess potentially different characteristics. This is especially true for climate‐sensitive species such as pikas, where the delimitation of distinct genetic units and the characterization of population responses to contemporary and historical environmental pressures are of particular interest. We combined a restriction site‐associated DNA sequencing (RADSeq) data set containing 4156 single nucleotide polymorphisms with ecological niche models (ENMs) of present and past habitat suitability to characterize population composition and evaluate the effects of historical range shifts, contemporary climates and landscape factors on gene flow in Collared Pikas, which are found in Alaska and adjacent regions of northwestern Canada and are the lesser‐studied of North America's two pika species. The results suggest that contemporary environmental factors contribute little to current population connectivity. Instead, genetic diversity is strongly shaped by the presence of three ancestral lineages isolated during the Pleistocene (~148 and 52 kya). Based on ENMs and genetic data, populations originating from a northern refugium experienced longer‐term stability, whereas both southern lineages underwent population expansion – contradicting the southern stability and northern expansion patterns seen in many other taxa. Current populations are comparable with respect to generally low diversity within populations and little‐to‐no recent admixture. The predominance of divergent histories structuring populations implies that if we are to understand and manage pika populations, we must specifically assess and accurately account for the forces underlying genetic similarity.