Watershed characteristics shape the landscape genetics of brook stickleback (Culaea inconstans) in shallow prairie lakes

Investigating the consequences of landscape features on population genetic patterns is increasingly important to elucidate the ecological factors governing connectivity between populations and predicting the evolutionary consequences of landscapes. Small prairie lakes in Alberta, Canada, and the brook stickleback (Culaea inconstans) that inhabit them, provide a unique aquatic system whereby populations are highly isolated from one another. These heterogeneous and extreme environments are prone to winterkills, an event whereby most of the fish die and frequent bottlenecks occur. In this study, we characterized the genetic population structure of brook stickleback among several lakes, finding that the species is hierarchically influenced by within‐lake characteristics in small‐scale watersheds. Landscape genetic analyses of the role of spatial features found support for basin characteristics associated with genetic diversity and bottlenecks in 20% of the sampled lakes. These results suggest that brook stickleback population genetic patterns may be driven, at least in part, by ecological processes that accelerate genetic drift and landscape patterns associated with reduced dispersal. Collectively, these results reinforce the potential importance of connectivity in the maintenance of genetic diversity, especially in fragmented landscapes.     

Tracking climate change in a dispersal-limited species: reduced spatial and genetic connectivity in a montane salamander

Tropical montane taxa are often locally adapted to very specific climatic conditions, contributing to their lower dispersal potential across complex landscapes. Climate and landscape features in montane regions affect population genetic structure in predictable ways, yet few empirical studies quantify the effects of both factors in shaping genetic structure of montane-adapted taxa. Here, we considered temporal and spatial variability in climate to explain contemporary genetic differentiation between populations of the montane salamander, Pseudoeurycea leprosa. Specifically, we used ecological niche modelling (ENM) and measured spatial connectivity and gene flow (using both mtDNA and microsatellite markers) across extant populations of P. leprosa in the Trans-Mexican Volcanic Belt (TVB). Our results indicate significant spatial and genetic isolation among populations, but we cannot distinguish between isolation by distance over time or current landscape barriers as mechanisms shaping population genetic divergences. Combining ecological niche modelling, spatial connectivity analyses, and historical and contemporary genetic signatures from different classes of genetic markers allows for inference of historical evolutionary processes and predictions of the impacts future climate change will have on the genetic diversity of montane taxa with low dispersal rates. Pseudoeurycea leprosa is one montane species among many endemic to this region and thus is a case study for the continued persistence of spatially and genetically isolated populations in the highly biodiverse TVB of central Mexico.     

Landscape genetics in a changing world: disentangling historical and contemporary influences and inferring change

Landscape genetics seeks to determine the effect of landscape features on gene flow and genetic structure. Often, such analyses are intended to inform conservation and management. However, depending on the many factors that influence the time to reach equilibrium, genetic structure may more strongly represent past rather than contemporary landscapes. This well‐known lag between current demographic processes and population genetic structure often makes it challenging to interpret how contemporary landscapes and anthropogenic activity shape gene flow. Here, we review the theoretical framework for factors that influence time lags, summarize approaches to address this temporal disconnect in landscape genetic studies, and evaluate ways to make inferences about landscape change and its effects on species using genetic data alone or in combination with other data. Those approaches include comparing correlation of genetic structure with historical versus contemporary landscapes, using molecular markers with different rates of evolution, contrasting metrics of genetic structure and gene flow that reflect population genetic processes operating at different temporal scales, comparing historical and contemporary samples, combining genetic data with contemporary estimates of species distribution or movement, and controlling for phylogeographic history. We recommend using simulated data sets to explore time lags in genetic structure, and argue that time lags should be explicitly considered both when designing and interpreting landscape genetic studies. We conclude that the time lag problem can be exploited to strengthen inferences about recent landscape changes and to establish conservation baselines, particularly when genetic data are combined with other data.     

Microgeographic population structure of brook charr: a comparison of microsatellite and mark-recapture data

Polymorphism at five microsatellite genetic markers (genotyped n  = 496) and mark-recapture tagging data (tagged n  = 9813) were used to define the population structure of brook charr, Salvelinus fontinalis from the Indian Bay watershed, Newfoundland, Canada. Despite the absence of physical barriers to migration among lakes, both genetic and tagging data suggest that brook charr in each lake represent reproductively isolated populations. Exact tests comparing allele frequencies, θ (global value = 0·063), R st (global value = 0·052), individual assignment tests, and Nei's genetic distance provided congruent estimates of population subdivision in agreement with the tagging data (only 2·2% of recaptures were lake-to-lake). The genetic structure of the brook charr populations corresponded with the geographic structure of the drainage basin on a qualitative level, although linear distance over water was not significantly correlated with the tagging data or the genetic distance measures. The agreement between the tagging and the genetic data suggest that microsatellite markers can be useful tools for defining real biological units. The results also suggest that brook charr exhibit microgeographic population structure at the watershed scale, and that this is the scale at which conservation and management of this salmonid might best be implemented.     

Effects of structural connectivity on fine scale population genetic structure of muskrat,Ondatra zibethicus

In heterogeneous landscapes, physical barriers and loss of structural connectivity have been shown to reduce gene flow and therefore lead to population structuring. In this study, we assessed the influence of landscape features on population genetic structure and gene flow of a semiaquatic species, the muskrat. A total of 97 muskrats were sampled from three watersheds near Sudbury, Ontario, Canada. We estimated population genetic structure using 11 microsatellite loci and identified a single genetic cluster and no genetic differences were found among the watersheds as a result of high levels of gene flow. At finer scales, we assessed the correlation between individual pairwise genetic distances and Euclidean distance as well as different models of least cost path (LCP). We used a range of cost values for the landscape types in order to build our LCP models. We found a positive relationship between genetic distance and least cost distance when we considered roads as corridors for movements. Open landscapes and urban areas seemed to restrict but not prevent gene flow within the study area. Our study underlines the high‐dispersal ability of generalist species in their use of landscape and highlights how landscape features often considered barriers to animal movements are corridors for other species.     

Disentangling the effects of historic vs. contemporary landscape structure on population genetic divergence

Increasing habitat fragmentation poses an immediate threat to population viability, as gene flow patterns are changed in these altered landscapes. Patterns of genetic divergence can potentially reveal the impact of these shifts in landscape connectivity. However, divergence patterns not only carry the signature of altered contemporary landscapes, but also historical ones. When considered separately, both recent and historical landscape structure appear to significantly affect connectivity among 51 wood frog ( Rana sylvatica ) populations. However, by controlling for correlations among landscape structure from multiple time periods, we show that patterns of genetic divergence reflect recent landscape structure as opposed to landscape structure prior to European settlement of the region (before 1850s). At the same time, within-population genetic diversities remain high and a genetic signature of population bottlenecks is lacking. Together, these results suggest that metapopulation processes – not drift-induced divergence associated with strong demographic bottlenecks following habitat loss – underlie the strikingly rapid consequences of temporally shifting landscape structure on these amphibians. We discuss the implications of these results in the context of understanding the role of population demography in the adaptive variation observed in wood frog populations.     

A portrait of a sucker using landscape genetics: how colonization and life history undermine the idealized dendritic metapopulation

Dendritic metapopulations have been attributed unique properties by in silico studies, including an elevated genetic diversity relative to a panmictic population of equal total size. These predictions have not been rigorously tested in nature, nor has there been full consideration of the interacting effects among contemporary landscape features, colonization history and life history traits of the target species. We tested for the effects of dendritic structure as well as the relative importance of life history, environmental barriers and historical colonization on the neutral genetic structure of a longnose sucker (Catostomus catostomus) metapopulation in the Kogaluk watershed of northern Labrador, Canada. Samples were collected from eight lakes, genotyped with 17 microsatellites, and aged using opercula. Lakes varied in differentiation, historical and contemporary connectivity, and life history traits. Isolation by distance was detected only by removing two highly genetically differentiated lakes, suggesting a lack of migration–drift equilibrium and the lingering influence of historical factors on genetic structure. Bayesian analyses supported colonization via the Kogaluk's headwaters. The historical concentration of genetic diversity in headwaters inferred by this result was supported by high historical and contemporary effective sizes of the headwater lake, T‐Bone. Alternatively, reduced allelic richness in headwaters confirmed the dendritic structure's influence on gene flow, but this did not translate to an elevated metapopulation effective size. A lack of equilibrium and upstream migration may have dampened the effects of dendritic structure. We suggest that interacting historical and contemporary factors prevent the achievement of the idealized traits of a dendritic metapopulation in nature.     

INTEGRATIVE TESTING OF HOW ENVIRONMENTS FROM THE PAST TO THE PRESENT SHAPE GENETIC STRUCTURE ACROSS LANDSCAPES

Tests of the genetic structure of empirical populations typically focus on the correlative relationships between population connectivity and geographic and/or environmental factors in landscape genetics. However, such tests may overlook or misidentify the impact of candidate factors on genetic structure, especially when connectivity patterns differ between past and present populations because of shifting environmental conditions over time. Here we account for the underlying demographic component of population connectivity associated with a temporarily dynamic landscape in tests of the factors structuring population genetic variation in an Australian lizard, Lerista lineopunctulata, from 24 nuclear loci. Correlative tests did not support significant effect from factors associated with a static contemporary landscape. However, spatially explicit demographic modeling of genetic differentiation shows that changes in environmental conditions (as estimated from paleoclimatic data) and corresponding distributional shifts from the past to present landscape significantly structures genetic variation. Results from model‐based inference (i.e., from an integrative modeling approach that generates spatially explicit expectations that are tested with approximate Bayesian computation) contrasts with those from correlative analyses, highlighting the importance of expanding the landscape genetic perspective to tests the links between pattern and process, revealing how factors shape patterns of genetic variation within species.     

Past forest cover explains current genetic differentiation in the Carpathian newt (Lissotriton montandoni), but not in the smooth newt (L. vulgaris)

Aim

Current genetic variation and differentiation are expected to reflect the effects of past rather than present landscapes due to time lags, that is, the time necessary for genetic diversity to reach equilibrium and reflect demography. Time lags can affect our ability to infer landscape use and model connectivity and also obscure the genetic consequences of recent landscape changes. In this work, we test whether past forest cover better explains contemporary patterns of genetic differentiation in two closely related but ecologically distinct newt species—Lissotriton montandoni and L. vulgaris.

Location

Carpathian Mountains and foothills.

Methods

Genetic differentiation between populations was related to landscape resistance optimized with tools from landscape genetics, for multiple timeframes, using forest-cover data from 1963 to 2015. Analyses were conducted for pairs of populations at distances from 1 to 50 km.

Results

We find evidence for a time lag in L. montandoni, with forest cover from 40 years ago (ca. 10 newt generations) better explaining current genetic differentiation. In L. vulgaris, current genetic differentiation was better predicted by present land-cover models with lower resistance given to open forests. This result may reflect the generalist ecology of L. vulgaris, its lower effective population sizes and exposure to habitat destruction and fragmentation.

Main Conclusions

Our study provides evidence for time lags in L. montandoni, showing that the genetic consequences of landscape change for some species are not yet evident. Our findings highlight the interspecific variation in time lag prevalence and demonstrate that current patterns of genetic differentiation should be interpreted in the context of historical landscape changes.     

Effects of landscape and history on diversification of a montane, stream-breeding amphibian

Aim The aim of this study was to understand the roles of landscape features in shaping patterns of contemporary and historical genetic diversification among populations of the Andean tree frog (Hypsiboas andinus) across spatial scales. Location Andes mountains, north‐western Argentina, South America. Methods Mitochondrial DNA control region sequences were utilized to assess genetic differentiation among populations and calculate population pair‐wise genetic distances. Three models of movement, namely traditional straight‐line distance and two effective distances based on habitat classification, were examined to determine which of these explained the most variation in pair‐wise population genetic differentiation. The two habitat classifications were based on digital vegetation and hydrology layers that were generated from a 90‐m resolution digital elevation model (DEM) and known relationships between elevation and habitat. Mantel tests were conducted to test for correlations between geographic and genetic distance matrices and to estimate the percentage variation explained by each type of geographic distance. To investigate the location of possible barriers to gene flow, we used Monmonier’s maximum difference algorithm as implemented in barrier 2.2. Results At both geographic scales, effective distances explained more variation in genetic differentiation than did straight‐line distance. The least‐cost distances based on the simple classification performed better than the more detailed habitat classification. We controlled for the effects of historical range fragmentation determined from previous nested clade analyses, and therefore evaluated the effect of different distances on the genetic variation attributable to more recent factors. Effective distances identified populations that were highly divergent as a result of isolation in unsuitable habitats. The proposed locations of barriers to gene flow identified using Monmonier’s maximum difference algorithm corresponded well with earlier analyses and supported findings from our partial Mantel tests. Main conclusions Our results indicate that landscape features have been important in both historical and contemporary genetic structuring of populations of H. andinus at both large and small spatial scales. A landscape genetic perspective offers novel insights not provided by traditional phylogeographic studies: (1) effective distances can better explain patterns of differentiation in populations, especially in heterogeneous landscapes where barriers to dispersal may be common; and (2) least‐cost path analysis can help to identify corridors of movement between populations that are biologically more realistic.     

Landscape genetics, historical isolation and cross-Andean gene flow in the wax palm, Ceroxylon echinulatum (Arecaceae)

Knowledge of the role of landscapes in shaping genetic connectivity and divergence is essential for understanding patterns of biogeography and diversity. This is particularly relevant for the Andes region, a major biodiversity hotspot of relatively recent origin. We examined the phylogeography and landscape genetics of the Andean wax palm Ceroxylon echinulatum (Arecaceae) that occurs in two narrow bands of montane forests on each side of the Andes in Ecuador and northeastern Peru. First, we tested the hypothesis of C. echinulatum being a geographic cline species crossing the Andes in the Amotape–Huancabamba zone (AHZ) of southern Ecuador/northern Peru, as indicated by observations on fruit morphology. Second, we assessed the timeframe of cross-Andean divergence, and third, we investigated the impact of contemporary and historical landscape features on observed spatio-genetic patterns. Individual-based Bayesian clustering (BC) identified a northeastern, southeastern, southwestern, and northwestern cluster, with areas of genetic discontinuity coinciding with the Andes and the Giron–Paute deflection. F -statistics derived from BC suggested an east-to-west dispersal history. Population-based analyses revealed strong genetic structuring at both small and large geographic scales. Interpopulation relationships and Mantel tests strongly supported the cline model with cross-Andean dispersal in the AHZ. Along the cline, gene flow measured as F _ST was mainly limited by distance, with less but significant impact of climatic friction. Coalescent analysis revealed that cross-Andean divergence took place during the Quaternary. Significant historical isolation ( R _ST >  F _ST) was found in the southwestern population. The current study illustrates a joint effect of founder dynamics, divergence by distance and historical isolation on patterns of Andean diversity and distribution.     

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.     

The influence of landscape on gene flow in the eastern massasauga rattlesnake (Sistrurus c. catenatus): insight from computer simulations

Understanding how gene flow shapes contemporary population structure requires the explicit consideration of landscape composition and configuration. New landscape genetic approaches allow us to link such heterogeneity to gene flow within and among populations. However, the attribution of cause is difficult when landscape features are spatially correlated, or when genetic patterns reflect past events. We use spatial Bayesian clustering and landscape resistance analysis to identify the landscape features that influence gene flow across two regional populations of the eastern massasauga rattlesnake, Sistrurus c. catenatus. Based on spatially explicit simulations, we inferred how habitat distribution modulates gene flow and attempted to disentangle the effects of spatially confounded landscape features. We found genetic clustering across one regional landscape but not the other, and also local differences in the effect of landscape on gene flow. Beyond the effects of isolation‐by‐distance, water bodies appear to underlie genetic differentiation among individuals in one regional population. Significant effects of roads were additionally detected locally, but these effects are possibly confounded with the signal of water bodies. In contrast, we found no signal of isolation‐by‐distance or landscape effects on genetic structure in the other regional population. Our simulations imply that these local differences have arisen as a result of differences in population density or tendencies for juvenile rather than adult dispersal. Importantly, our simulations also demonstrate that the ability to detect the consequences of contemporary anthropogenic landscape features (e.g. roads) on gene flow may be compromised when long‐standing natural features (e.g. water bodies) co‐exist on the landscape.     

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.     

The scale of genetic differentiation in the Dunes Sagebrush-Lizard (Sceloporus arenicolus), an endemic habitat specialist

The Dunes Sagebrush-Lizard (Sceloporus arenicolus) is a North American species endemic to sand-shinnery oak habitats of the Mescalero and Monahans sand dunes in eastern New Mexico and western Texas. This lizard is listed as Endangered in New Mexico and exhibits habitat specificity at several geographic scales. Dunes Sagebrush-Lizards are only found in topographically complex shinnery oak (Quercus havardii) dominated landscapes within their small geographic distribution and are not found in surrounding human-altered landscapes. Within suitable sand-shinnery oak habitat, individuals predominantly occupy non-vegetated sand dune blowouts and utilize blowouts with particular physical characteristics due to thermoregulatory, reproduction, and foraging requirements. Here, we examined historical and contemporary patterns of genetic differentiation with respect to the current distribution of suitable habitat at multiple spatial scales using mitochondrial DNA sequences and microsatellite data from individuals throughout the entire range. We found three genetic clusters of individuals generally concordant with geographic regions and low sequence divergence at mitochondrial loci suggesting a recent origin of these populations. We also found high levels of genetic structure at microsatellite loci among populations within each of these groups indicating restricted gene flow at intermediate scales. Despite high habitat specificity, we did not detect genetic structure among sand blowouts at finer spatial scales. Within each population, matrices comprised of both sand blowouts and vegetated shinnery oak patches are necessary for genetic connectivity, but the fine scale spatial arrangement of blowouts may not be as critical. We discuss our results with respect to the scale of landscape heterogeneity and habitat connectivity and consider the conservation implications for this threatened taxon.     

Fragmentation of Atlantic Forest has not affected gene flow of a widespread seed‐dispersing bat

Habitat loss and resultant fragmentation are major threats to biodiversity, particularly in tropical and subtropical ecosystems. It is increasingly urgent to understand fragmentation effects, which are often complex and vary across taxa, time and space. We determined whether recent fragmentation of Atlantic forest is causing population subdivision in a widespread and important Neotropical seed disperser: Artibeus lituratus (Chiroptera: Phyllostomidae). Genetic structure within highly fragmented forest in Paraguay was compared to that in mostly contiguous forest in neighbouring Misiones, Argentina. Further, observed genetic structure across the fragmented landscape was compared with expected levels of structure for similar time spans in realistic simulated landscapes under different degrees of reduction in gene flow. If fragmentation significantly reduced successful dispersal, greater population differentiation and stronger isolation by distance would be expected in the fragmented than in the continuous landscape, and genetic structure in the fragmented landscape should be similar to structure for simulated landscapes where dispersal had been substantially reduced. Instead, little genetic differentiation was observed, and no significant correlation was found between genetic and geographic distance in fragmented or continuous landscapes. Furthermore, comparison of empirical and simulated landscapes indicated empirical results were consistent with regular long‐distance dispersal and high migration rates. Our results suggest maintenance of high gene flow for this relatively mobile and generalist species, which could be preventing or significantly delaying reduction in population connectivity in fragmented habitat. Our conclusions apply to A. lituratus in Interior Atlantic Forest, and do not contradict broad evidence that habitat fragmentation is contributing to extinction of populations and species, and poses a threat to biodiversity worldwide.     

Contemporary and historic factors influence differently genetic differentiation and diversity in a tropical palm

Population genetics theory predicts loss in genetic variability because of drift and inbreeding in isolated plant populations; however, it has been argued that long-distance pollination and seed dispersal may be able to maintain gene flow, even in highly fragmented landscapes. We tested how historical effective population size, historical migration and contemporary landscape structure, such as forest cover, patch isolation and matrix resistance, affect genetic variability and differentiation of seedlings in a tropical palm (Euterpe edulis) in a human-modified rainforest. We sampled 16 sites within five landscapes in the Brazilian Atlantic forest and assessed genetic variability and differentiation using eight microsatellite loci. Using a model selection approach, none of the covariates explained the variation observed in inbreeding coefficients among populations. The variation in genetic diversity among sites was best explained by historical effective population size. Allelic richness was best explained by historical effective population size and matrix resistance, whereas genetic differentiation was explained by matrix resistance. Coalescence analysis revealed high historical migration between sites within landscapes and constant historical population sizes, showing that the genetic differentiation is most likely due to recent changes caused by habitat loss and fragmentation. Overall, recent landscape changes have a greater influence on among-population genetic variation than historical gene flow process. As immediate restoration actions in landscapes with low forest amount, the development of more permeable matrices to allow the movement of pollinators and seed dispersers may be an effective strategy to maintain microevolutionary processes.     

Genetic signature of population fragmentation varies with mobility in seven bird species of a fragmented Kenyan cloud forest

Habitat fragmentation can restrict geneflow, reduce neighbourhood effective population size, and increase genetic drift and inbreeding in small, isolated habitat remnants. The extent to which habitat fragmentation leads to population fragmentation, however, differs among landscapes and taxa. Commonly, researchers use information on the current status of a species to predict population effects of habitat fragmentation. Such methods, however, do not convey information on species-specific responses to fragmentation. Here, we compare levels of past population differentiation, estimated from microsatellite genotypes, with contemporary dispersal rates, estimated from multi-strata capture-recapture models, to infer changes in mobility over time in seven sympatric, forest-dependent bird species of a Kenyan cloud forest archipelago. Overall, populations of sedentary species were more strongly differentiated and clustered compared to those of vagile ones, while geographical patterning suggested an important role of landscape structure in shaping genetic variation. However, five of seven species with broadly similar levels of genetic differentiation nevertheless differed substantially in their current dispersal rates. We conclude that post-fragmentation levels of vagility, without reference to past population connectivity, may not be the best predictor of how forest fragmentation affects the life history of forest-dependent species. As effective conservation strategies often hinge on accurate prediction of shifts in ecological and genetic relationships among populations, conservation practices based solely upon current population abundances or movements may, in the long term, prove to be inadequate.     

Microsatellite variation and genetic structure of brook trout (Salvelinus fontinalis) populations in Labrador and neighboring Atlantic Canada: evidence for ongoing gene flow and dual routes of post-Wisconsinan colonization

In conservation genetics and management, it is important to understand the contribution of historical and contemporary processes to geographic patterns of genetic structure in order to characterize and preserve diversity. As part of a 10-year monitoring program by the Government of Newfoundland and Labrador, Canada, we measured the population genetic structure of the world's most northern native populations of brook trout (Salvelinus fontinalis) in Labrador to gather baseline data to facilitate monitoring of future impacts of the recently opened Trans-Labrador Highway. Six-locus microsatellite profiles were obtained from 1130 fish representing 32 populations from six local regions. Genetic diversity in brook trout populations in Labrador (average H(E)= 0.620) is within the spectrum of variability found in other brook trout across their northeastern range, with limited ongoing gene flow occurring between populations (average pairwise F(ST)= 0.139). Evidence for some contribution of historical processes shaping genetic structure was inferred from an isolation-by-distance analysis, while dual routes of post-Wisconsinan recolonization were indicated by STRUCTURE analysis: K= 2 was the most likely number of genetic groups, revealing a separation between northern and west-central Labrador from all remaining populations. Our results represent the first data from the nuclear genome of brook trout in Labrador and emphasize the usefulness of microsatellite data for revealing the extent to which genetic structure is shaped by both historical and contemporary processes.