Genetic structure of the world's polar bear populations

We studied genetic structure in polar bear (Ursus maritimus) populations by typing a sample of 473 individuals spanning the species distribution at 16 highly variable microsatellite loci. No genetic discontinuities were found that would be consistent with evolutionarily significant periods of isolation between groups. Direct comparison of movement data and genetic data from the Canadian Arctic revealed a highly significant correlation. Genetic data generally supported existing population (management unit) designations, although there were two cases where genetic data failed to differentiate between pairs of populations previously resolved by movement data. A sharp contrast was found between the minimal genetic structure observed among populations surrounding the polar basin and the presence of several marked genetic discontinuities in the Canadian Arctic. The discontinuities in the Canadian Arctic caused the appearance of four genetic clusters of polar bear populations. These clusters vary in total estimated population size from 100 to over 10 000, and the smallest may merit a relatively conservative management strategy in consideration of its apparent isolation. We suggest that the observed pattern of genetic discontinuities has developed in response to differences in the seasonal distribution and pattern of sea ice habitat and the effects of these differences on the distribution and abundance of seals.     

Drawing inferences about the coancestry coefficient

The coancestry coefficient, also known as the population structure parameter, is of great interest in population genetics. It can be thought of as the intraclass correlation of pairs of alleles within populations and it can serve as a measure of genetic distance between populations. For a general class of evolutionary models it determines the distribution of allele frequencies among populations. Under more restrictive models it can be regarded as the probability of identity by descent of any pair of alleles at a locus within a random mating population. In this paper we review estimation procedures that use the method of moments or are maximum likelihood under the assumption of normally distributed allele frequencies. We then consider the problem of testing hypotheses about this parameter. In addition to parametric and non-parametric bootstrap tests we present an asymptotically-distributed chi-square test. This test reduces to the contingency-table test for equal sample sizes across populations. Our new test appears to be more powerful than previous tests, especially for loci with multiple alleles. We apply our methods to HapMap SNP data to confirm that the coancestry coefficient for humans is strictly positive.     

Deciphering ecological barriers to North American river otter (Lontra canadensis) gene flow in the Louisiana landscape

For North American river otters (Lontra canadensis) in Louisiana, statewide distribution, availability of aquatic habitats, and the absence of physical barriers to dispersal might suggest that they exist as a large, panmictic population. However, the wide variety of habitat types in this region, and the dynamic nature of these habitats over time, could potentially structure river otter populations in accordance with cryptic landscape features. Recently developed landscape genetic models offer a spatially explicit approach that could be useful in identifying potential barriers to the movement of river otters through the dynamic aquatic landscape of Louisiana. We used georeferenced multilocus microsatellite genotypes in spatially implicit (STRUCTURE) and spatially explicit (GENELAND) models to characterize patterns of landscape genetic structure. All models identified 3 subpopulations of river otters in Louisiana, corresponding to Inland, Atchafalaya River, and Mississippi River regions. Variation in breeding seasonality, brought about by variation in prey abundance between inland and coastal populations, may have contributed to genetic differentiation among populations. It is also possible that the genetic discontinuities we observed indicate a correlation between otter distribution and access to freshwater. Regardless of the mechanism, it is likely that any genetic differentiation among subpopulations is exacerbated by relatively poor dispersal.     

A new method to uncover signatures of divergent and stabilizing selection in quantitative traits

While it is well understood that the pace of evolution depends on the interplay between natural selection, random genetic drift, mutation, and gene flow, it is not always easy to disentangle the relative roles of these factors with data from natural populations. One popular approach to infer whether the observed degree of population differentiation has been influenced by local adaptation is the comparison of neutral marker gene differentiation (as reflected in FST) and quantitative trait divergence (as reflected in QST). However, this method may lead to compromised statistical power, because FST and QST are summary statistics which neglect information on specific pairs of populations, and because current multivariate tests of neutrality involve an averaging procedure over the traits. Further, most FST-QST comparisons actually replace QST by its expectation over the evolutionary process and are thus theoretically flawed. To overcome these caveats, we derived the statistical distribution of population means generated by random genetic drift and used the probability density of this distribution to test whether the observed pattern could be generated by drift alone. We show that our method can differentiate between genetic drift and selection as a cause of population differentiation even in cases with FST=QST and demonstrate with simulated data that it disentangles drift from selection more accurately than conventional FST-QST tests especially when data sets are small.     

Congruent population structure inferred from dispersal behaviour and intensive genetic surveys of the threatened Florida scrub-jay (Aphelocoma coerulescens)

The delimitation of populations, defined as groups of individuals linked by gene flow, is possible by the analysis of genetic markers and also by spatial models based on dispersal probabilities across a landscape. We combined these two complimentary methods to define the spatial pattern of genetic structure among remaining populations of the threatened Florida scrub-jay, a species for which dispersal ability is unusually well-characterized. The range-wide population was intensively censused in the 1990s, and a metapopulation model defined population boundaries based on predicted dispersal-mediated demographic connectivity. We subjected genotypes from more than 1000 individual jays screened at 20 microsatellite loci to two Bayesian clustering methods. We describe a consensus method for identifying common features across many replicated clustering runs. Ten genetically differentiated groups exist across the present-day range of the Florida scrub-jay. These groups are largely consistent with the dispersal-defined metapopulations, which assume very limited dispersal ability. Some genetic groups comprise more than one metapopulation, likely because these genetically similar metapopulations were sundered only recently by habitat alteration. The combined reconstructions of population structure based on genetics and dispersal-mediated demographic connectivity provide a robust depiction of the current genetic and demographic organization of this species, reflecting past and present levels of dispersal among occupied habitat patches. The differentiation of populations into 10 genetic groups adds urgency to management efforts aimed at preserving what remains of genetic variation in this dwindling species, by maintaining viable populations of all genetically differentiated and geographically isolated populations.     

Phylogeography of the montane caddisfly Drusus discolor: evidence for multiple refugia and periglacial survival

We studied the genetic population structure and phylogeography of the montane caddisfly Drusus discolor across its entire range in central and southern Europe. The species is restricted to mountain regions and exhibits an insular distribution across the major mountain ranges. Mitochondrial sequence data (COI) of 254 individuals from the entire species range is analysed to reveal population genetic structure. The data show little molecular variation within populations and regions, but distinct genetic differentiation between mountain ranges. Most populations are significantly differentiated based on F(ST) and exact tests of population differentiation and most haplotypes are unique to a single mountain range. Phylogenetic analyses reveal deep divergence between geographically isolated lineages. Combined, these results suggest that past fragmentation is the prominent process structuring the populations across Europe. We use tests of selective neutrality and mismatch distributions, to study the demographic population history of regions with haplotype overlap. The high level of genetic differentiation between mountain ranges and estimates of demographic history provide evidence for the existence of multiple glacial refugia, including several in central Europe. The study shows that these aquatic organisms reacted differently to Pleistocene cooling than many terrestrial species. They persisted in numerous refugia over multiple glacial cycles, allowing many local endemic clades to form.     

Stronger spatial genetic structure in recolonized areas than in refugia in the European beech

Extant rear‐edge populations located in former glacial refugia remain understudied despite their high conservation value. These populations should have experienced strong genetic drift due to their small size and long isolation. Moreover, the prolonged action of isolation by distance in refugial areas should result in stronger regional spatial genetic structure (SGS) than in recolonized areas, but empirical tests of this prediction are scarce. To fill this gap, we first used a set of 16 microsatellite markers to investigate the genetic structure of European beech in France in 65 populations from three refugial areas and one control recolonized (nonrefugial) area. Then, using the same approach, we reanalysed published isozyme data from 375 populations distributed across the entire species range. We found stronger genetic differentiation among populations in refugia than in recolonized areas. However, contrary to expectations, regional SGS was lower within refugia than within recolonized areas. Published studies presenting similar analyses suggest that our results could have generality across different biogeographical settings and types of organisms. Strong and prolonged genetic drift in refugial areas could have erased the signature of range expansions that is still visible in recolonized areas. Our results therefore suggest that Pleistocene population isolation has played a key role in increasing the genetic complexity of extant rear‐edge populations.     

Integrative species delimitation reveals cryptic diversity in the southern Appalachian Antrodiaetus unicolor (Araneae: Antrodiaetidae) species complex

Although species delimitation can be highly contentious, the development of reliable methods to accurately ascertain species boundaries is an imperative step in cataloguing and describing Earth's quickly disappearing biodiversity. Spider species delimitation remains largely based on morphological characters; however, many mygalomorph spider populations are morphologically indistinguishable from each other yet have considerable molecular divergence. The focus of our study, the Antrodiaetus unicolor species complex containing two sympatric species, exhibits this pattern of relative morphological stasis with considerable genetic divergence across its distribution. A past study using two molecular markers, COI and 28S, revealed that A. unicolor is paraphyletic with respect to A. microunicolor. To better investigate species boundaries in the complex, we implement the cohesion species concept and use multiple lines of evidence for testing genetic exchangeability and ecological interchangeability. Our integrative approach includes extensively sampling homologous loci across the genome using a RADseq approach (3RAD), assessing population structure across their geographic range using multiple genetic clustering analyses that include structure , principal components analysis and a recently developed unsupervised machine learning approach (Variational Autoencoder). We evaluate ecological similarity by using large‐scale ecological data for niche‐based distribution modelling. Based on our analyses, we conclude that this complex has at least one additional species as well as confirm species delimitations based on previous less comprehensive approaches. Our study demonstrates the efficacy of genomic‐scale data for recognizing cryptic species, suggesting that species delimitation with one data type, whether one mitochondrial gene or morphology, may underestimate true species diversity in morphologically homogenous taxa with low vagility.     

Isolation by resistance across a complex coral reef seascape

A detailed understanding of the genetic structure of populations and an accurate interpretation of processes driving contemporary patterns of gene flow are fundamental to successful spatial conservation management. The field of seascape genetics seeks to incorporate environmental variables and processes into analyses of population genetic data to improve our understanding of forces driving genetic divergence in the marine environment. Information about barriers to gene flow (such as ocean currents) is used to define a resistance surface to predict the spatial genetic structure of populations and explain deviations from the widely applied isolation-by-distance model. The majority of seascape approaches to date have been applied to linear coastal systems or at large spatial scales (more than 250 km), with very few applied to complex systems at regional spatial scales (less than 100 km). Here, we apply a seascape genetics approach to a peripheral population of the broadcast-spawning coral Acropora spicifera across the Houtman Abrolhos Islands, a high-latitude complex coral reef system off the central coast of Western Australia. We coupled population genetic data from a panel of microsatellite DNA markers with a biophysical dispersal model to test whether oceanographic processes could explain patterns of genetic divergence. We identified significant variation in allele frequencies over distances of less than 10 km, with significant differentiation occurring between adjacent sites but not between the most geographically distant ones. Recruitment probabilities between sites based on simulated larval dispersal were projected into a measure of resistance to connectivity that was significantly correlated with patterns of genetic divergence, demonstrating that patterns of spatial genetic structure are a function of restrictions to gene flow imposed by oceanographic currents. This study advances our understanding of the role of larval dispersal on the fine-scale genetic structure of coral populations across a complex island system and applies a methodological framework that can be tailored to suit a variety of marine organisms with a range of life-history characteristics.     

Number of SNPS loci needed to detect population structure

The study of the association of polymorphic genetic markers with common diseases is one of the most powerful tools in modern genetics. Interest in single nucleotide polymorphisms (SNPs) has steadily grown over the last decade. SNPs are currently the most developed markers in the human genome because they have a number of advantages over other marker types. One of the critical problems responsible for 'spurious' association findings in case-control studies is population stratification. There are many statistical approaches developed for detecting population heterogeneity. However the power to detect population structure by known methods is highly dependent on the number of loci utilised. We performed an analysis of SNPs data available in the public domain from The Single Nucleotide Consortia Ltd. (TSCL). Three populations, Afro-American, Asian and Caucasian, were compared. Estimation of the minimum number of SNPs loci necessary for detection of the population structure was performed. Two clustering approaches, distance-based and model-based, were compared. The model-based approach was superior when compared with the distance-based method. We found more than 65 random SNPs loci are required for identifying distinct geographically separated populations. Increasing the number of markers to over 100 raises the probability of correct assignment of a particular individual to an origin group to over 90%, even with conventional clustering methods.     

Migratory Flyways May Affect Population Structure in Double-Crested Cormorants

Double-crested cormorants (Phalacrocorax auritus) recovered from a demographic bottleneck so well that they are now considered a nuisance species at breeding and wintering grounds across the United States and Canada. Management of this species could be improved by refining genetic population boundaries and assigning individuals to their natal population. Further, recent radio-telemetry data suggest the existence of Interior and Atlantic migratory flyways, which could reduce gene flow and result in substantial genetic isolation. In this study, we used 1,784 individuals collected across the eastern United States, a large panel of microsatellite markers developed for this species, and individuals banded as chicks and recaptured as adults to explore the effects of migratory flyways on population structure, quantify the genetic effects of demographic bottlenecks, and determine whether individuals could be assigned to their natal population based on genotype. We found evidence for genetic population division only along migratory flyways, no evidence of genetic bottlenecks, and mixed effectiveness of assignment tests. Our population structure findings suggest that gene flow is high across large scales; for example, individuals from New York, Minnesota, and Alabama are all in panmixia. We also found that traditional subspecies ranges may not be valid because >1 subspecies was present in single genetic populations. The lack of evidence for genetic bottlenecks also likely underscores the vagility of this species, suggesting that even during demographic bottlenecks, populations were not isolated from allelic exchange. Finally, the failure of assignment tests to consistently perform is likely due in part to imperfect a priori sampling of Atlantic and Interior chicks and the high vagility of adults. We conclude that the demographic bottleneck is not likely to have reduced genetic diversity, and that assignment tests remain unreliable for this species. We recommend double-crested cormorants be managed by flyway. Further development of genomic resources in this species could improve population subdivision resolution, improve assignment tests, and reveal further information on demographic histories. © 2020 The Wildlife Society.     

Integration of genetic and demographic data to assess population risk in a continuously distributed species

The identification and demographic assessment of biologically meaningful populations is fundamental to species’ ecology and management. Although genetic tools are used frequently to identify populations, studies often do not incorporate demographic data to understand their respective population trends. We used genetic data to define subpopulations in a continuously distributed species. We assessed demographic independence and variation in population trends across the distribution. Additionally, we identified potential barriers to gene flow among subpopulations. We sampled greater sage-grouse (Centrocercus urophasianus) leks from across their range (≈175,000 Km²) in Wyoming and amplified DNA at 14 microsatellite loci for 1761 samples. Subsequently, we assessed population structure in unrelated individuals (n = 872) by integrating results from multiple Bayesian clustering approaches and used the boundaries to inform our assessment of long-term population trends and lek activity over the period of 1995–2013. We identified four genetic clusters of which two northern ones showed demographic independence from the others. Trends in population size for the northwest subpopulation were statistically different from the other three genetic clusters and the northeast and southwest subpopulations demonstrated a general trend of increasing proportion of inactive leks over time. Population change from 1996 to 2012 suggested population growth in the southern subpopulations and decline, or neutral, change in the northern subpopulations. We suggest that sage-grouse subpopulations in northern Wyoming are at greater risk of extirpation than the southern subpopulations due to smaller census and effective population sizes and higher variability within subpopulations. Our research is an example of incorporating genetic and demographic data and provides guidance on the identification of subpopulations of conservation concern.     

A consensus tree approach for reconstructing human evolutionary history and detecting population substructure

The random accumulation of variations in the human genome over time implicitly encodes a history of how human populations have arisen, dispersed, and intermixed since we emerged as a species. Reconstructing that history is a challenging computational and statistical problem but has important applications both to basic research and to the discovery of genotype-phenotype correlations. We present a novel approach to inferring human evolutionary history from genetic variation data. We use the idea of consensus trees, a technique generally used to reconcile species trees from divergent gene trees, adapting it to the problem of finding robust relationships within a set of intraspecies phylogenies derived from local regions of the genome. Validation on both simulated and real data shows the method to be effective in recapitulating known true structure of the data closely matching our best current understanding of human evolutionary history. Additional comparison with results of leading methods for the problem of population substructure assignment verifies that our method provides comparable accuracy in identifying meaningful population subgroups in addition to inferring relationships among them. The consensus tree approach thus provides a promising new model for the robust inference of substructure and ancestry from large-scale genetic variation data.     

Defining population structure for the Mojave desert tortoise

We used highly variable microsatellite markers to identify population structure, movement, and biological boundaries for populations of the desert tortoise, Gopherus agassizii, in the Mojave and Colorado Deserts of the southwestern United States. The Mojave desert tortoise (listed as “threatened” by the U.S. Fish and Wildlife Service) has a large geographic range, long generation time, low population densities, and little above-ground activity. Additionally, the dispersal patterns of individual tortoises are virtually unknown, making indirect methods to assess movement among populations valuable. Using Bayesian assignment tests, we detected hierarchical structuring within the Mojave desert tortoise. Three basal groups were identified, and these corresponded to the mitochondrial DNA haplotypes reported in 1989. Additional population structure was evident within each basal unit, and this structure corresponds with major geographic barriers. Our analyses suggest that gene flow among populations was historically high because levels of population differentiation were low across the range. Geographic distance explained a large proportion of variation in genetic distance (68%), which pinpoints that dispersal is limited only on a regional scale. In light of these new analyses of the genetic population structure of the Mojave desert tortoise, we make new recommendations for the number and locations of recovery units for conservation of this species.     

Genetic Structuring across Marine Biogeographic Boundaries in Rocky Shore Invertebrates

Biogeography investigates spatial patterns of species distribution. Discontinuities in species distribution are identified as boundaries between biogeographic areas. Do these boundaries affect genetic connectivity? To address this question, a multifactorial hierarchical sampling design, across three of the major marine biogeographic boundaries in the central Mediterranean Sea (Ligurian-Tyrrhenian, Tyrrhenian-Ionian and Ionian-Adriatic) was carried out. Mitochondrial COI sequence polymorphism of seven species of Mediterranean benthic invertebrates was analysed. Two species showed significant genetic structure across the Tyrrhenian-Ionian boundary, as well as two other species across the Ionian Sea, a previously unknown phylogeographic barrier. The hypothesized barrier in the Ligurian-Tyrrhenian cannot be detected in the genetic structure of the investigated species. Connectivity patterns across species at distances up to 800 km apart confirmed that estimates of pelagic larval dispersal were poor predictors of the genetic structure. The detected genetic discontinuities seem more related to the effect of past historical events, though maintained by present day oceanographic processes. Multivariate statistical tools were used to test the consistency of the patterns across species, providing a conceptual framework for across-species barrier locations and strengths. Additional sequences retrieved from public databases supported our findings. Heterogeneity of phylogeographic patterns shown by the 7 investigated species is relevant to the understanding of the genetic diversity, and carry implications for conservation biology.     

Examining genetic structure in a bogus yucca moth: a sequential approach to phylogeography

Understanding the phylogeography of a species requires not only elucidating patterns of genetic structure among populations, but also identifying the possible evolutionary events creating that structure. The use of a single phylogeographic test or analysis, however, usually provides a picture of genetic structure without revealing the possible underlying evolutionary causes. We used current analytical techniques in a sequential approach to examine genetic structure and its underlying causes in the bogus yucca moth Prodoxus decipiens (Lepidoptera: Prodoxidae). Both historical biogeography and recent human transplantations of the moth's host plants provided a priori expectations of the pattern of genetic structure and its underlying causes. We evaluated these expectations by using a progression of phylogenetic, demographic, and population genetic analyses of mtDNA sequence data from 476 individuals distributed across 25 populations that encompassed the range of P. decipiens. The combination of these analyses revealed that much of the genetic structure has evolved more recently than suggested by historical biogeography, has been influenced by changes in demography, and can be best explained by long distance dispersal and isolation by distance. We suggest that performing a suite of analyses that focus on different temporal scales may be an effective approach to investigating the patterns and causes of genetic structure within species.     

Matching genetics with oceanography: directional gene flow in a Mediterranean fish species

Genetic connectivity and geographic fragmentation are two opposing mechanisms determining the population structure of species. While the first homogenizes the genetic background across populations the second one allows their differentiation. Therefore, knowledge of processes affecting dispersal of marine organisms is crucial to understand their genetic distribution patterns and for the effective management of their populations. In this study, we use genetic analyses of eleven microsatellites in combination with oceanographic satellite and dispersal simulation data to determine distribution patterns for Serranus cabrilla, a ubiquitous demersal broadcast spawner, in the Mediterranean Sea. Pairwise population F_ST values ranged between ?0.003 and 0.135. Two genetically distinct clusters were identified, with a clear division located between the oceanographic discontinuities at the Ibiza Channel (IC) and the Almeria‐Oran Front (AOF), revealing an admixed population in between. The Balearic Front (BF) also appeared to dictate population structure. Directional gene flow on the Spanish coast was observed as S. cabrilla dispersed from west to east over the AOF, from north to south on the IC and from south of the IC towards the Balearic Islands. Correlations between genetic and oceanographic data were highly significant. Seasonal changes in current patterns and the relationship between ocean circulation patterns and spawning season may also play an important role in population structure around oceanographic fronts.     

Bayesian clustering algorithms ascertaining spatial population structure: a new computer program and a comparison study

On the basis of simulated data, this study compares the relative performances of the Bayesian clustering computer programs structure , geneland , geneclust and a new program named tess . While these four programs can detect population genetic structure from multilocus genotypes, only the last three ones include simultaneous analysis from geographical data. The programs are compared with respect to their abilities to infer the number of populations, to estimate membership probabilities, and to detect genetic discontinuities and clinal variation. The results suggest that combining analyses using tess and structure offers a convenient way to address inference of spatial population structure.     

An approximate likelihood for genetic data under a model with recombination and population splitting

We describe a new approximate likelihood for population genetic data under a model in which a single ancestral population has split into two daughter populations. The approximate likelihood is based on the ‘Product of Approximate Conditionals’ likelihood and ‘copying model’ of Li and Stephens [Li, N., Stephens, M., 2003. Modeling linkage disequilibrium and identifying recombination hotspots using single-nucleotide polymorphism data. Genetics 165 (4), 2213–2233]. The approach developed here may be used for efficient approximate likelihood-based analyses of unlinked data. However our copying model also considers the effects of recombination. Hence, a more important application is to loosely-linked haplotype data, for which efficient statistical models explicitly featuring non-equilibrium population structure have so far been unavailable. Thus, in addition to the information in allele frequency differences about the timing of the population split, the method can also extract information from the lengths of haplotypes shared between the populations. There are a number of challenges posed by extracting such information, which makes parameter estimation difficult. We discuss how the approach could be extended to identify haplotypes introduced by migrants.     

Assessing changes in functional connectivity in a desert bighorn sheep metapopulation after two generations

Determining how species move across complex and fragmented landscapes and interact with human‐made barriers is a major research focus in conservation. Studies estimating functional connectivity from movement, dispersal or gene flow usually rely on a single study period and rarely consider variation over time. We contrasted genetic structure and gene flow across barriers for a metapopulation of desert bighorn sheep (Ovis canadensis nelsoni) using genotypes collected 2000–2003 and 2013–2015. Based on the recently observed but unexpected spread of a respiratory pathogen across an interstate highway previously identified as a barrier to gene flow, we hypothesized that bighorn sheep changed how they interacted with that barrier, and that shifts in metapopulation structure influenced gene flow, genetic diversity and connectivity. Population assignment tests, genetic structure and genetic recapture demonstrated that bighorn sheep crossed the interstate highway in at least one location in 2013–2015, sharply reducing genetic structure between two populations, but supported conclusions of an earlier study that such crossings were very infrequent or unknown in 2000–2003. A recently expanded population established new links and caused decreases in genetic structure among multiple populations. Genetic diversity showed only slight increases in populations linked by new connections. Genetic structure and assignments revealed other previously undetected changes in movements and distribution, but much was consistent. Thus, we observed changes in both structural and functional connectivity over just two generations, but only in specific locations. Movement patterns of species should be revisited periodically to enable informed management, particularly in dynamic and fragmented systems.