Microsatellite analysis of population structure in Canadian polar bears

Attempts to study the genetic population structure of large mammals are often hampered by the low levels of genetic variation observed in these species. Polar bears have particularly low levels of genetic variation with the result that their genetic population structure has been intractable. We describe the use of eight hypervariable microsatellite loci to study the genetic relationships between four Canadian polar bear populations: the northern Beaufort Sea, southern Beaufort Sea, western Hudson Bay, and Davis Strait - Labrador Sea. These markers detected considerable genetic variation, with average heterozygosity near 60% within each population. Interpopulation differences in allele frequency distribution were significant between all pairs of populations, including two adjacent populations in the Beaufort Sea. Measures of genetic distance reflect the geographic distribution of populations, but also suggest patterns of gene flow which are not obvious from geography and may reflect movement patterns of these animals. Distribution of variation is sufficiently different between the Beaufort Sea populations and the two more eastern ones that the region of origin for a given sample can be predicted based on its expected genotype frequency using an assignment test. These data indicate that gene flow between local populations is restricted despite the long-distance seasonal movements undertaken by polar bears.     

Life-history and habitat features influence the within-river genetic structure of Atlantic salmon

Defining populations and identifying ecological and life-history characteristics affecting genetic structure is important for understanding species biology and hence, for managing threatened or endangered species or populations. In this study, populations of the world's largest indigenous Atlantic salmon (Salmo salar) stock were first inferred using model-based clustering methods, following which life-history and habitat variables best predicting the genetic diversity of populations were identified. This study revealed that natal homing of Atlantic salmon within the Teno River system is accurate at least to the tributary level. Generally, defining populations by main tributaries was observed to be a reasonable approach in this large river system, whereas in the mainstem of the river, the number of inferred populations was fewer than the number of distinct sampling sites. Mainstem and headwater populations were genetically more diverse and less diverged, while each tributary fostered a distinct population with high genetic differentiation and lower genetic diversity. Population structure and variation in genetic diversity among populations were poorly explained by geographical distance. In contrast, age-structure, as estimated by the proportion of multisea-winter spawners, was the most predictive variable in explaining the variation in the genetic diversity of the populations. This observation, being in agreement with theoretical predictions, emphasizes the essence of large multisea-winter females in maintaining the genetic diversity of populations. In addition, the unique genetic diversity of populations, as estimated by private allele richness, was affected by the ease of accessibility of a site, with more difficult to access sites having lower unique genetic diversity. Our results show that despite this species' high capacity for migration, tributaries foster relatively closed populations with little gene flow which will be important to consider when developing management strategies for the system.     

Development of microsatellite markers in Acer capillipes

Acer capillipes is an insect‐pollinated tree species that grows in temperate regions of Japan. We isolated seven polymorphic microsatellite loci from this species using a dual‐suppression polymerase chain reaction (PCR) technique. The number of alleles per locus ranged from two to 10, and the expected heterozygosity ranged from 0.042 to 0.828. Cross‐species amplification from 14 other Acer species was successful for the majority of the isolated loci, suggesting that these loci may be useful for the characterization of other maple species.     

Panmictic population structure in the hooded seal (Cystophora cristata)

Two putative populations of hooded seals (Cystophora cristata) occur in the North Atlantic. The Greenland Sea population pup and breed on the pack ice near Jan Mayen ('West Ice') while the Northwest Atlantic population is thought to pup in the Davis Strait, in the Gulf of St. Lawrence (the 'Gulf'), and off southern Labrador or northeast Newfoundland (the 'Front'). We used microsatellite profiling of 300 individuals at 13 loci and mitochondrial DNA sequencing of the control region of 123 individuals to test for genetic differentiation between these four breeding herds. We found no significant genetic differences between breeding areas, nor evidence for cryptic nor higher level genetic structure in this species. The Greenland Sea breeding herd was genetically most distant from the Northwest Atlantic breeding areas; however, the differences were statistically nonsignificant. Our data therefore suggest that the world's hooded seals comprise a single panmictic genetic population.     

In situ population structure and ex situ representation of the endangered Amur tiger

The Amur tiger ( Panthera tigris altaica ) is a critically endangered felid that suffered a severe demographic contraction in the 1940s. In this study, we sampled 95 individuals collected throughout their native range to investigate questions relative to population genetic structure and demographic history. Additionally, we sampled targeted individuals from the North American ex situ population to assess the genetic representation found in captivity. Population genetic and Bayesian structure analyses clearly identified two populations separated by a development corridor in Russia. Despite their well-documented 20th century decline, we failed to find evidence of a recent population bottleneck, although genetic signatures of a historical contraction were detected. This disparity in signal may be due to several reasons, including historical paucity in population genetic variation associated with postglacial colonization and potential gene flow from a now extirpated Chinese population. Despite conflicting signatures of a bottleneck, our estimates of effective population size ( N _e = 27–35) and N _e /N ratio (0.07–0.054) were substantially lower than the only other values reported for a wild tiger population. Lastly, the extent and distribution of genetic variation in captive and wild populations were similar, yet gene variants persisted ex situ that were lost in situ . Overall, our results indicate the need to secure ecological connectivity between the two Russian populations to minimize loss of genetic diversity and overall susceptibility to stochastic events, and support a previous study suggesting that the captive population may be a reservoir of gene variants lost in situ .     

Sampling related individuals within ponds biases estimates of population structure in a pond‐breeding amphibian

Effective conservation and management of pond‐breeding amphibians depends on the accurate estimation of population structure, demographic parameters, and the influence of landscape features on breeding‐site connectivity. Population‐level studies of pond‐breeding amphibians typically sample larval life stages because they are easily captured and can be sampled nondestructively. These studies often identify high levels of relatedness between individuals from the same pond, which can be exacerbated by sampling the larval stage. Yet, the effect of these related individuals on population genetic studies using genomic data is not yet fully understood. Here, we assess the effect of within‐pond relatedness on population and landscape genetic analyses by focusing on the barred tiger salamanders (Ambystoma mavortium) from the Nebraska Sandhills. Utilizing genome‐wide SNPs generated using a double‐digest RADseq approach, we conducted standard population and landscape genetic analyses using datasets with and without siblings. We found that reduced sample sizes influenced parameter estimates more than the inclusion of siblings, but that within‐pond relatedness led to the inference of spurious population structure when analyses depended on allele frequencies. Our landscape genetic analyses also supported different models across datasets depending on the spatial resolution analyzed. We recommend that future studies not only test for relatedness among larval samples but also remove siblings before conducting population or landscape genetic analyses. We also recommend alternative sampling strategies to reduce sampling siblings before sequencing takes place. Biases introduced by unknowingly including siblings can have significant implications for population and landscape genetic analyses, and in turn, for species conservation strategies and outcomes.     

Population genetics of the people of Iran II. Genetic differentiation and population structure

The genetic structure of six populations of Iran (Turks, Kurds, Lurs, Zabolis, Baluchis and Zoroastrians) was examined using data on blood groups, serum proteins and cell enzymes. Our results show conclusively that there are genetic differences among the six populations and the analysis of superimposed R and S matrices defined Harpending & Jenkins (1973) show by that the dispersion of some of the alleles correspond to the dispersion of the populations. The F_ST estimates are not large enough to favour selection on any of the loci studied. The F_IT and F_IS estimates are positive and moderately high suggesting that the genetic differentiation to some extent is influenced by inbreeding.     

The phylogeographical and management implications of genetic population structure in the imperiled snuffbox mussel, Epioblasma triquetra (Bivalvia: Unionidae)

Seven populations of the imperiled snuffbox mussel, Epioblasma triquetra , were sampled from across the central basin of North America. Samples were genotyped using 15 microsatellite DNA loci, and maternal history was inferred using mitochondrial DNA (mtDNA) cytochrome  c oxidase subunit-I (COI) sequences. Populations in the Clinch and St Francis rivers were quite distinct in their mtDNA. The population in the St Francis River had a unique, fixed haplotype. Among a suite of haplotypes, the population in the Clinch River had two unique haplotypes of common ancestry. The other populations were dominated by a common haplotype, which also occurred in the Clinch River population. Analysis of DNA microsatellites revealed much greater divergences and showed significant genetic structure between populations in the formerly glaciated regions. Divergence has occurred between the populations, as evidenced by moderate to high fixation indices ( F _ST and R _ST values) and nearly perfect assignment tests. These results indicate the occurrence of three glacial refugia for E. triquetra : the Tennessee River, rivers south of the Ozark Crest, and the lower Ohio River drainage near the confluence with the Mississippi. Populations in the lower Ohio River were likely to be responsible for the postglacial reinvasion into formerly glaciated regions, and into the upper Tennessee River drainage. The population of the St Francis River may constitute a distinct taxonomic entity. Conservation efforts, if necessary for this imperiled species, should not mix populations.  © 2008 The Linnean Society of London, Biological Journal of the Linnean Society , 2008, 93 , 371–384.     

Population structure of North American beluga whales (Delphinapterus leucas) based on nuclear DNA microsatellite variation and contrasted with the population structure revealed by mitochondrial DNA variation

Beluga whales ( Delphinapterus leucas ) in North American waters migrate seasonally between wintering areas in broken pack ice and summering locations in estuaries and other open water areas in the Arctic and sub-Arctic. Results from our previous investigation of beluga whale mitochondrial DNA (mtDNA) revealed genetic heterogeneity among beluga from different summering locations that was interpreted as representing a high degree of summering site philopatry. However, mtDNA is maternally inherited and does not reflect mating that may occur among beluga from different summering locations in wintering areas or during annual migrations. To test the possibility that breeding occurs among beluga from different summering locations, genetic variability at five nuclear DNA (nDNA) microsatellite loci was examined in the same animals tested in the mtDNA study. Beluga samples ( n = 640) were collected between 1984 and 1994 from 24 sites across North America, mostly during the summer. Whales from the various sites were categorized into eight summering locations as identified by mtDNA analysis, as well as four hypothesized wintering areas: Bering Sea, Hudson Strait (Hudson Strait, Labrador Sea, southwest Davis Strait), Baffin Bay (North Water, east Davis Strait), and St Lawrence River. Microsatellite allele frequencies indicated genetic homogeneity among animals from summering sites believed to winter together but differentiation among whales from some of the wintering areas. In particular, beluga from western North America (Bering Sea) were clearly distinguished from beluga from eastern North America (Hudson Strait, Baffin Bay, and St Lawrence River). Based upon the combined data set, the population of North American beluga whales was divided into two evolutionarily significant units. However, the population may be further subdivided into management units to reflect distinct groups of beluga at summering locations.     

Population structure of North American beluga whales (Delphinapterus leucas) based on nuclear DNA microsatellite variation and contrasted with the population structure revealed by mitochondrial DNA variation

Beluga whales ( Delphinapterus leucas ) in North American waters migrate seasonally between wintering areas in broken pack ice and summering locations in estuaries and other open water areas in the Arctic and sub-Arctic. Results from our previous investigation of beluga whale mitochondrial DNA (mtDNA) revealed genetic heterogeneity among beluga from different summering locations that was interpreted as representing a high degree of summering site philopatry. However, mtDNA is maternally inherited and does not reflect mating that may occur among beluga from different summering locations in wintering areas or during annual migrations. To test the possibility that breeding occurs among beluga from different summering locations, genetic variability at five nuclear DNA (nDNA) microsatellite loci was examined in the same animals tested in the mtDNA study. Beluga samples ( n = 640) were collected between 1984 and 1994 from 24 sites across North America, mostly during the summer. Whales from the various sites were categorized into eight summering locations as identified by mtDNA analysis, as well as four hypothesized wintering areas: Bering Sea, Hudson Strait (Hudson Strait, Labrador Sea, southwest Davis Strait), Baffin Bay (North Water, east Davis Strait), and St Lawrence River. Microsatellite allele frequencies indicated genetic homogeneity among animals from summering sites believed to winter together but differentiation among whales from some of the wintering areas. In particular, beluga from western North America (Bering Sea) were clearly distinguished from beluga from eastern North America (Hudson Strait, Baffin Bay, and St Lawrence River). Based upon the combined data set, the population of North American beluga whales was divided into two evolutionarily significant units. However, the population may be further subdivided into management units to reflect distinct groups of beluga at summering locations.     

The landscape genetics of yellow perch (Perca flavescens) in a large fluvial ecosystem

Landscape genetics is being increasingly applied to elucidate the role of environmental features on the population structure of terrestrial organisms. However, the potential of this framework has been little explored in aquatic ecosystems such as large rivers. Here, we used a landscape genetics approach in order to (i) document the population structure of the yellow perch ( Perca flavescens ) by means of genetic variation at microsatellite markers, (ii) assess to what extent the structure was explained by landscape heterogeneity, and (iii) interpret the relevance of interactions between genetics and landscape for management and conservation. Analysis of the genetic variation among 1715 individuals from 16 localities and distributed over 310 km in the freshwater section of the Saint Lawrence River (Québec, Canada) revealed a relatively modest level of genetic structuring ( F _ST = 0.039). Application of the Monmonier's algorithm combining geographical and genetic information identified three zones of restricted gene flow defining four distinct populations. Physical barriers played a more important role on gene flow and genetic structure than waterway geographical distance. We found correlations between genetic differentiation and presence of distinct water masses in the sector of Lake Saint-Louis ( r =  0.7177, P  = 0.0340) and with fragmentation of spawning habitats in the sector of Lake Saint-Pierre ( r =  0.8578, P  = 0.0095). Our results support the treatment of four distinct biological units, which is in contrast with the current basis for yellow perch management. Finally, this study showed that landscape genetics is a powerful means to identify environmental barriers to gene flow causing genetic discontinuities in apparently highly connected aquatic landscapes.     

Landscape genetics and least-cost path analysis reveal unexpected dispersal routes in the California tiger salamander (Ambystoma californiense)

A major goal of landscape genetics is to understand how landscapes structure genetic variation in natural populations. However, landscape genetics still lacks a framework for quantifying the effects of landscape features, such as habitat type, on realized gene flow. Here, we present a methodology for identifying the costs of dispersal through different habitats for the California tiger salamander ( Ambystoma californiense ), an endangered species restricted to grassland/vernal pool habitat mosaics. We sampled larvae from all 16 breeding ponds in a geographically restricted area of vernal pool habitat at the Fort Ord Natural Reserve, Monterey County, California. We estimated between-pond gene flow using 13 polymorphic microsatellite loci and constructed GIS data layers of habitat types in our study area. We then used least-cost path analysis to determine the relative costs of movement through each habitat that best match rates of gene flow measured by our genetic data. We identified four measurable rates of gene flow between pairs of ponds, with between 10.5% and 19.9% of larvae having immigrant ancestry. Although A. californiense is typically associated with breeding ponds in grassland habitat, we found that dispersal through grassland is nearly twice as costly as dispersal through chaparral and that oak woodland is by far the most costly habitat to traverse. With the increasing availability of molecular resources and GIS data, we anticipate that these methods could be applied to a broad range of study systems, particularly those with cryptic life histories that make direct observation of movement challenging.     

Subtle population structuring within a highly vagile marine invertebrate, the veined squid Loligo forbesi, demonstrated with microsatellite DNA markers

Microsatellite DNA markers were applied for the first time in a population genetic study of a cephalopod and compared with previous estimates of genetic differentiation obtained using allozyme and mitochondrial DNA (mtDNA) markers. Levels of genetic variation detected with microsatellites were much higher than found with previous markers (mean number of alleles per locus=10.6, mean expected heterozygosity ( H _E)=0.79; allozyme H _E=0.08; mtDNA restriction fragment length polymorphism (RFLP) H _E=0.16). In agreement with previous studies, microsatellites demonstrated genetic uniformity across the population occupying the European shelf seas of the North East Atlantic, and extreme genetic differentiation of the Azores population ( R _ST/ F _ST=0.252/0.245; allozyme F _ST=0.536; mtDNA F _ST=0.789). In contrast to other markers, microsatellites detected more subtle, and significant, levels of differentiation between the populations of the North East Atlantic offshore banks (Rockall and Faroes) and the shelf population ( R _ST=0.048 and 0.057). Breakdown of extensive gene flow among these populations is indicated, with hydrographic (water depth) and hydrodynamic (isolating current regimes) factors suggested as possible barriers to migration. The demonstration of genetic subdivision in an abundant, highly mobile marine invertebrate has implications for the interpretation of dispersal and population dynamics, and consequent management, of such a commercially exploited species. Relative levels of differentiation indicated by the three different marker systems, and the use of measures of differentiation (assuming different mutation models), are discussed.     

Influence of landscape on the population genetic structure of the alpine butterfly parnassius smintheus (Papilionidae)

Four microsatellite DNA markers were developed which were used to examine the relationship between landscape and population genetic structure among a set of populations of the butterfly Parnassius smintheus located in the foothills of the Canadian Rockies. Detailed information on the dispersal of adult butterflies among this same set of populations was available. Simple and partial Mantel tests were used to examine the relationships between genetic distances, predicted rates of dispersal, and a number of landscape variables, all measured pairwise for 17 sample sites. Nei's standard genetic distance was negatively correlated with predicted dispersal. We observed a significant pattern of isolation by distance at a very small spatial scale. The distance between sites that was through forest was a stronger predictor of genetic distance than the distance through open meadow, indicating a significant effect of landscape on population genetic structure beyond that of simple isolation by distance. Our results suggest that rises in the tree-line in alpine areas, caused by global warming, will lead to reduced gene flow among populations of P. smintheus.     

Landscape structure and the genetic effects of a population collapse

Both landscape structure and population size fluctuations influence population genetics. While independent effects of these factors on genetic patterns and processes are well studied, a key challenge is to understand their interaction, as populations are simultaneously exposed to habitat fragmentation and climatic changes that increase variability in population size. In a population network of an alpine butterfly, abundance declined 60–100% in 2003 because of low over-winter survival. Across the network, mean microsatellite genetic diversity did not change. However, patch connectivity and local severity of the collapse interacted to determine allelic richness change within populations, indicating that patch connectivity can mediate genetic response to a demographic collapse. The collapse strongly affected spatial genetic structure, leading to a breakdown of isolation-by-distance and loss of landscape genetic pattern. Our study reveals important interactions between landscape structure and temporal demographic variability on the genetic diversity and genetic differentiation of populations. Projected future changes to both landscape and climate may lead to loss of genetic variability from the studied populations, and selection acting on adaptive variation will likely occur within the context of an increasing influence of genetic drift.     

Fine‐scale population structure of Collichtys lucidus populations inferred from microsatellite markers

The spinyhead croaker Collichthys lucidus (Richardson) is a small sciaenid species distributed along the inshore waters of northwestern Pacific Ocean, and now has been listed as Key Protected Commercial Sources of Aquatic Animals and Plants in China. To delineate stock boundaries and inform conservation policy for its management, samples were collected from eight locations across the Chinese coastal waters and analyzed at nine microsatellite loci. C. lucidus populations showed low genetic diversity (expected heterozygosity = 0.445–0.542; observed heterozygosity = 0.392–0.539; Polymorphism Information Content = 0.268–0.684). Strong genetic fdifferentiation (F_st = 0.065–0.510, all significant after Bonferroni correction) among all populations and high levels of self‐recruitment (89.2%–91.5%) were observed, which suggested limited genetic exchange for this species. Clustering results of discriminant analysis of principal components and STRUCTURE found strong support for obvious genetic clusters (populations FZ, XM and SZ vs. populations SH, YRE, ZS, WZ and ND). The results of the present study not only supported the phylogeographic pattern of north‐south differentiation, but also suggested that C. lucidus populations may be predominantly sustained by self‐replenishment rather than by recruitment from distant populations.     

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.     

Demographics and landscape features determine intrariver population structure in Atlantic salmon (Salmo salar L.): the case of the River Moy in Ireland

Contemporary genetic structure of Atlantic salmon (Salmo salar L.) in the River Moy in Ireland is shown here to be strongly related to landscape features and population demographics, with populations being defined largely by their degree of physical isolation and their size. Samples of juvenile salmon were collected from the 17 major spawning areas on the river Moy and from one spawning area in each of five smaller nearby rivers. No temporal allele frequency differences were observed within locations for 12 microsatellite loci, whereas nearly all spatial samples differed significantly, suggesting that each was a separate population. Bayesian clustering and landscape genetic analyses suggest that these populations can be combined hierarchically into five genetically informative larger groupings. Lakes were found to be the single most important determinant of the observed population structure. Spawning area size was also an important factor. The salmon population of the closest nearby river resembled genetically the largest Moy population grouping. In addition, we showed that anthropogenic influences on spawning habitats, in this case arterial drainage, can affect relationships between populations. Our results show that Atlantic salmon biodiversity can be largely defined by geography, and thus, knowledge of landscape features (for example, as characterized within Geographical Information Systems) has the potential to predict population structure in other rivers without an intensive genetic survey, or at least to help direct sampling. This approach of combining genetics and geography, for sampling and in subsequent statistical analyses, has wider application to the investigation of population structure in other freshwater/anadromous fish species and possibly in marine fish and other organisms.     

Landscape genetics and hierarchical genetic structure in Atlantic salmon: the interaction of gene flow and local adaptation

Disentangling evolutionary forces that may interact to determine the patterns of genetic differentiation within and among wild populations is a major challenge in evolutionary biology. The objective of this study was to assess the genetic structure and the potential influence of several ecological variables on the extent of genetic differentiation at multiple spatial scales in a widely distributed species, the Atlantic salmon, Salmo salar . A total of 2775 anadromous fish were sampled from 51 rivers along the North American Atlantic coast and were genotyped using 13 microsatellites. A Bayesian analysis clustered these populations into seven genetically and geographically distinct groups, characterized by different environmental and ecological factors, mainly temperature. These groups were also characterized by different extent of genetic differentiation among populations. Dispersal was relatively high and of the same magnitude within compared to among regional groups, which contrasted with the maintenance of a regional genetic structure. However, genetic differentiation was lower among populations exchanging similar rates of local as opposed to inter-regional migrants, over the same geographical scale. This raised the hypothesis that gene flow could be constrained by local adaptation at the regional scale. Both coastal distance and temperature regime were found to influence the observed genetic structure according to landscape genetic analyses. The influence of other factors such as latitude, river length and altitude, migration tactic, and stocking was not significant at any spatial scale. Overall, these results suggested that the interaction between gene flow and thermal regime adaptation mainly explained the hierarchical genetic structure observed among Atlantic salmon populations.     

Isolation and characterization of microsatellite loci in Geum urbanum (Rosaceae) and their transferability within the genus Geum

Thirteen novel polymorphic microsatellite loci are presented for Geum urbanum (Rosaceae). The microsatellites will be useful tools to analyse the influence of landscape structure and land‐use intensity in agricultural landscapes on genetic diversity within and among populations of Geum urbanum. Transferability was tested in 19 other Geum species and two Waldsteinia species. In most species polymerase chain reaction (PCR) products of the expected range were obtained, therefore the markers reported here appear to be applicable across the whole genus.