Genetic sex assignment in wild populations using genotyping‐by‐sequencing data: A statistical threshold approach

Establishing the sex of individuals in wild systems can be challenging and often requires genetic testing. Genotyping‐by‐sequencing (GBS) and other reduced‐representation DNA sequencing (RRS) protocols (e.g., RADseq, ddRAD) have enabled the analysis of genetic data on an unprecedented scale. Here, we present a novel approach for the discovery and statistical validation of sex‐specific loci in GBS data sets. We used GBS to genotype 166 New Zealand fur seals (NZFS, Arctocephalus forsteri) of known sex. We retained monomorphic loci as potential sex‐specific markers in the locus discovery phase. We then used (i) a sex‐specific locus threshold (SSLT) to identify significantly male‐specific loci within our data set; and (ii) a significant sex‐assignment threshold (SSAT) to confidently assign sex in silico the presence or absence of significantly male‐specific loci to individuals in our data set treated as unknowns (98.9% accuracy for females; 95.8% for males, estimated via cross‐validation). Furthermore, we assigned sex to 86 individuals of true unknown sex using our SSAT and assessed the effect of SSLT adjustments on these assignments. From 90 verified sex‐specific loci, we developed a panel of three sex‐specific PCR primers that we used to ascertain sex independently of our GBS data, which we show amplify reliably in at least two other pinniped species. Using monomorphic loci normally discarded from large SNP data sets is an effective way to identify robust sex‐linked markers for nonmodel species. Our novel pipeline can be used to identify and statistically validate monomorphic and polymorphic sex‐specific markers across a range of species and RRS data sets.     

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

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

Comparison of historical bottleneck effects and genetic consequences of re‐introduction in a critically endangered island passerine

Re‐introduction is an important tool for recovering endangered species; however, the magnitude of genetic consequences for re‐introduced populations remains largely unknown, in particular the relative impacts of historical population bottlenecks compared to those induced by conservation management. We characterize 14 microsatellite loci developed for the Seychelles paradise flycatcher and use them to quantify temporal and spatial measures of genetic variation across a 134‐year time frame encompassing a historical bottleneck that reduced the species to ~28 individuals in the 1960s, through the initial stages of recovery and across a second contemporary conservation‐introduction‐induced bottleneck. We then evaluate the relative impacts of the two bottlenecks, and finally apply our findings to inform broader re‐introduction strategy. We find a temporal trend of significant decrease in standard measures of genetic diversity across the historical bottleneck, but only a nonsignificant downward trend in number of alleles across the contemporary bottleneck. However, accounting for the different timescales of the two bottlenecks (~40 historical generations versus <1 contemporary generation), the loss of genetic diversity per generation is greater across the contemporary bottleneck. Historically, the flycatcher population was genetically structured; however, extinction on four of five islands has resulted in a homogeneous contemporary population. We conclude that severe historical bottlenecks can leave a large footprint in terms of sheer quantity of genetic diversity lost. However, severely depleted genetic diversity does not render a species immune to further genetic erosion upon re‐introduction. In some cases, the loss of genetic diversity per generation can, initially at least, be greater across re‐introduction‐induced bottlenecks.     

High gene flow and metapopulation dynamics detected for three species in a dryland river system

1. The highly variable hydrology of dryland rivers has important implications for population dynamics in these systems. In western Queensland, fluctuations in sub‐population size are likely to lead to local bottlenecks and extinctions, increasing the need for connectivity and gene flow to maintain population viability. 2. Using microsatellite markers, we explored evidence for this metapopulation structure in two species of freshwater fish (Maquaria ambigua and Tandanus tandanus) and one crustacean (Macrobrachium australiense) in a sub‐catchment of the upper Murray–Darling Basin, Australia. 3. Overall, we found very weak genetic structure for all three species. Two species (M. ambigua and M. australiense) showed some significant genetic structure that did not correlate with geographic distance. However, decomposed pairwise regression analysis revealed evidence for intense genetic drift at the waterhole scale, suggesting that local bottlenecks are driving what little genetic structure does exist for these species. 4. The results identify the local impact of bottlenecks on genetic diversity, but highlight the importance of gene flow in maintaining population viability in these highly variable systems. As the impacts of bottlenecks are likely to be tempered by gene flow, it is suggested that the maintenance of connectivity is of paramount importance in this dryland system.     

Comparative phylogeography of mainland and insular species of Neotropical molossid bats (Molossus)

Historical events, habitat preferences, and geographic barriers might result in distinct genetic patterns in insular versus mainland populations. Comparison between these two biogeographic systems provides an opportunity to investigate the relative role of isolation in phylogeographic patterns and to elucidate the importance of evolution and demographic history in population structure. Herein, we use a genotype‐by‐sequencing approach (GBS) to explore population structure within three species of mastiff bats (Molossus molossus, M. coibensis, and M. milleri), which represent different ecological histories and geographical distributions in the genus. We tested the hypotheses that oceanic straits serve as barriers to dispersal in Caribbean bats and that isolated island populations are more likely to experience genetic drift and bottlenecks in comparison with highly connected ones, thus leading to different phylogeographic patterns. We show that population structures vary according to general habitat preferences, levels of population isolation, and historical fluctuations in climate. In our dataset, mainland geographic barriers played only a small role in isolation of lineages. However, oceanic straits posed a partial barrier to the dispersal for some populations within some species (M. milleri), but do not seem to disrupt gene flow in others (M. molossus). Lineages on distant islands undergo genetic bottlenecks more frequently than island lineages closer to the mainland, which have a greater exchange of haplotypes.     

Defining evolutionary boundaries across parapatric ecomorphs of Black Salamanders (Aneides flavipunctatus) with conservation implications

The accurate delimitation of evolutionary population units represents an important component in phylogeographic and conservation genetic studies. Here, we used a combined population assignment and historical demographic approach to study a complex of ecomorphologically distinctive populations of Black Salamanders (Aneides flavipunctatus) that are parapatrically distributed and meet at a three‐way contact zone in north‐western California. We used mitochondrial tree‐based and multilocus clustering methods to evaluate a priori two‐ (Northern and Southern) and three (Northern, Coast and Inland) population hypotheses derived from previous studies. Mitochondrial results were consistent with the two‐ and three‐population hypotheses, while the nDNA clustering results supported only the two‐population hypothesis. Historical demographic analyses and mtDNA gene divergence estimates revealed that the Northern and Southern populations split during the Pliocene (2–5 Ma). Subdivision of the Southern population into Coast and Inland populations was estimated to be late Pleistocene (0.24 Ma), although our mtDNA results suggested a Pliocene divergence. Effective gene flow estimates (2N_em) suggest that either the two‐ or three‐population hypotheses remain valid. However, our results unexpectedly revealed that the Northern population might instead represent two parapatric populations that separated nearly 4 Ma. These results are surprising because the Pliocene divergence between these ecomorphologically conservative forms is similar or older than for the ecomorphologically divergent Coast and Inland sister populations. We conclude that Black Salamanders in north‐western California belong to at least three or four populations or species, and these all meet criteria for being Evolutionary Significant Units or ‘ESUs’ and therefore warrant conservation consideration.     

Drought Impacts Survivorship and Reproductive Strategies of an Epiphyllous Leafy Liverwort in Central Amazonia

Understanding how changing precipitation patterns impact the population dynamics of Amazonian plants is necessary for predicting their long‐term survival. Most terrestrial plants are characterized by life cycles intractably slow for evaluating such demographic consequences. Due to the demands of surviving on an ephemeral substrate such as a leaf surface, epiphyllous (leaf‐inhabiting) bryophytes have some of the fastest generation times for terrestrial plants. Considering these advantages, we investigated the demographic consequences of seasonality in precipitation during one‐year of the epiphyllous bryophyte Radula flaccida Gott. (Radulaceae) in a central Amazonian rain forest. In a 1‐ha sized study area, 154 epiphyll colonies from 18 phorophytes of the understory shrub Naucleopsis ulei (Warb.) Ducke (Moraceae) were marked and colony growth, extinction, recolonization, as well as rates of sexual and asexual expression were measured in nine demographic censuses. The probability of survivorship decreased in the dry season due to the synergistic effects of both increased leaf fall (e.g., loss of occupied patches) and colony mortality on persisting leaves. Asexual expression also declined significantly during the driest months; however, sexual expression, as measured by sporophyte density, was not related to seasonality. Logistic regression probabilities also indicate a minimum threshold colony size for optimal sporophyte expression. Lower survivorship and colony growth in the dry season along with reduced sporophyte production in smaller colonies implies that both demographic stability as well as population structure of Amazonian epiphyll species may be compromised whether climate models projecting prolonged droughts for the region are accurate.     

Outlier analyses to test for local adaptation to breeding grounds in a migratory arctic seabird

Investigating the extent (or the existence) of local adaptation is crucial to understanding how populations adapt. When experiments or fitness measurements are difficult or impossible to perform in natural populations, genomic techniques allow us to investigate local adaptation through the comparison of allele frequencies and outlier loci along environmental clines. The thick‐billed murre (Uria lomvia) is a highly philopatric colonial arctic seabird that occupies a significant environmental gradient, shows marked phenotypic differences among colonies, and has large effective population sizes. To test whether thick‐billed murres from five colonies along the eastern Canadian Arctic coast show genomic signatures of local adaptation to their breeding grounds, we analyzed geographic variation in genome‐wide markers mapped to a newly assembled thick‐billed murre reference genome. We used outlier analyses to detect loci putatively under selection, and clustering analyses to investigate patterns of differentiation based on 2220 genomewide single nucleotide polymorphisms (SNPs) and 137 outlier SNPs. We found no evidence of population structure among colonies using all loci but found population structure based on outliers only, where birds from the two northernmost colonies (Minarets and Prince Leopold) grouped with birds from the southernmost colony (Gannet), and birds from Coats and Akpatok were distinct from all other colonies. Although results from our analyses did not support local adaptation along the latitudinal cline of breeding colonies, outlier loci grouped birds from different colonies according to their non‐breeding distributions, suggesting that outliers may be informative about adaptation and/or demographic connectivity associated with their migration patterns or nonbreeding grounds.     

Supportive breeding boosts natural population abundance with minimal negative impacts on fitness of a wild population of Chinook salmon

While supportive breeding programmes strive to minimize negative genetic impacts to populations, case studies have found evidence for reduced fitness of artificially produced individuals when they reproduce in the wild. Pedigrees of two complete generations were tracked with molecular markers to investigate differences in reproductive success (RS) of wild and hatchery‐reared Chinook salmon spawning in the natural environment to address questions regarding the demographic and genetic impacts of supplementation to a natural population. Results show a demographic boost to the population from supplementation. On average, fish taken into the hatchery produced 4.7 times more adult offspring, and 1.3 times more adult grand‐offspring than naturally reproducing fish. Of the wild and hatchery fish that successfully reproduced, we found no significant differences in RS between any comparisons, but hatchery‐reared males typically had lower RS values than wild males. Mean relative reproductive success (RRS) for hatchery F₁ females and males was 1.11 (P = 0.84) and 0.89 (P = 0.56), respectively. RRS of hatchery‐reared fish (H) that mated in the wild with either hatchery or wild‐origin (W) fish was generally equivalent to W × W matings. Mean RRS of H × W and H × H matings was 1.07 (P = 0.92) and 0.94 (P = 0.95), respectively. We conclude that fish chosen for hatchery rearing did not have a detectable negative impact on the fitness of wild fish by mating with them for a single generation. Results suggest that supplementation following similar management practices (e.g. 100% local, wild‐origin brood stock) can successfully boost population size with minimal impacts on the fitness of salmon in the wild.     

Contrasting Demographic Structure of Short‐ and Long‐lived Pioneer Tree Species on Amazonian Forest Edges

Although tropical forests have been rapidly converted into human‐modified landscapes, tree species response to forest edges remains poorly examined. In this study, we addressed four pioneer tree species to document demographic shifts experienced by this key ecological group and make inferences about pioneer response to forest edges. All individuals with dbh ≥ 1 cm of two short‐lived (Bellucia grossularioides and Cecropia sciadophylla) and two long‐lived species (Goupia glabra and Laetia procera) were sampled in 20 1‐ha forest edge plots and 20 1‐ha forest interior plots in Oiapoque and Manaus, Northeast and Central Amazon, respectively. As expected, pioneer stem density with dbh ≥ 1 cm increased by around 10–17‐fold along forest edges regardless of species, lifespan, and study site. Edge populations of long‐lived pioneers presented 84–94 percent of their individuals in sapling/subadult size classes, whereas edge populations of short‐lived pioneers showed 56–97 percent of their individuals in adult size classes. These demographic biases were associated with negative and positive net adult recruitment of long‐ and short‐lived pioneers, respectively. Our population‐level analyses support three general statements: (1) native pioneer tree species proliferate along forest edges (i.e., increased density), at least in terms of non‐reproductive individuals; (2) pioneer response to edge establishment is not homogeneous as species differ in terms of demographic structure and net adult recruitment; and (3) some pioneer species, particularly long‐lived ones, may experience population decline due to adult sensitivity to edge‐affected habitats.     

Genetic depletion at adaptive but not neutral loci in an endangered bird species

Many endangered species suffer from the loss of genetic diversity, but some populations may be able to thrive even if genetically depleted. To investigate the underlying genetic processes of population bottlenecks, we apply an innovative approach for assessing genetic diversity in the last known population of the endangered Pale‐headed Brushfinch (Atlapetes pallidiceps) in Ecuador. First, we measure genetic diversity at eleven neutral microsatellite loci and adaptive SNP variation in five Toll‐like receptor (TLR) immune system genes. Bottleneck tests confirm genetic drift as the main force shaping genetic diversity in this species and indicate a 99 % reduction in population size dating back several hundred years. Second, we compare contemporary microsatellite diversity with historic museum samples of A. pallidiceps, finding no change in genetic diversity. Third, we compare genetic diversity in the Pale‐headed Brushfinch with two co‐occurring‐related brushfinch species (Atlapetes latinuchus, Buarremon torquatus), finding a reduction of up to 91% diversity in the immune system genes but not in microsatellites. High TLR diversity is linked to decreased survival probabilities in A. pallidiceps. Low TLR diversity is thus probably an adaptation to the specific selection regime within its currently very restricted distribution (approximately 200 ha), but could severely restrict the adaptive potential of the species in the long run. Our study illustrates the importance of investigating both neutral and adaptive markers to assess the effect of population bottlenecks and for recommending specific management plans in endangered species.     

Why to account for finite sites in population genetic studies and how to do this with Jaatha 2.0

With the advent of next‐generation sequencing technologies, large data sets of several thousand loci from multiple conspecific individuals are available. Such data sets should make it possible to obtain accurate estimates of population genetic parameters, even for complex models of population history. In the analyses of large data sets, it is difficult to consider finite‐sites mutation models (FSMs). Here, we use extensive simulations to demonstrate that the inclusion of FSMs is necessary to avoid severe biases in the estimation of the population mutation rate θ, population divergence times, and migration rates. We present a new version of Jaatha, an efficient composite‐likelihood method for estimating demographic parameters from population genetic data and evaluate the usefulness of Jaatha in two biological examples. For the first application, we infer the speciation process of two wild tomato species, Solanum chilense and Solanum peruvianum. In our second application example, we demonstrate that Jaatha is readily applicable to NGS data by analyzing genome‐wide data from two southern European populations of Arabidopsis thaliana. Jaatha is now freely available as an R package from the Comprehensive R Archive Network (CRAN).     

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.     

Postglacial history of a widespread conifer produces inverse clines in selective neutrality tests

Deviations of the site frequency spectrum of mutations (SFS) from neutral expectations may be caused by natural selection or by demographic processes such as population subdivision or temporal changes in population size. As most widespread temperate and boreal tree species have expanded from glacial refugia in the past 13 000 years, colonization bottlenecks associated with this migration may have left variable demographic signatures among geographic populations corresponding to distance from the refugia. To determine whether the signature of postglacial re‐colonization has skewed the SFS in the widely distributed conifer Sitka spruce (Picea sitchensis (Bong.) Carr.), we re‐sequenced 153 nuclear genes in six populations from across the species range. We found that while the SFS for the pooled sample produced negative values for Tajima’s D and Fay and Wu’s H, these statistics exhibited strong clinal variation when populations were analysed separately (R² = 0.84, P = 0.007 for Tajima’s D and R² = 0.65, P = 0.033 for Fay and Wu’s H). When historical bottlenecks of varying age were simulated using approximate Bayesian computation, distance of populations from the southern range limit explained most of the variation in bottleneck timing among populations (R² = 0.89, P = 0.003). These data suggest that sequential population bottlenecks during postglacial re‐colonization have resulted in diverse among‐population signatures within the contemporary SFS in Sitka spruce, with rare variants more common in the south, and medium‐frequency variants more common in the north. Our results also emphasize the need to consider sampling strategy and to explore population‐specific null demographic models in surveys of nucleotide variation in widely distributed species.     

Two genetically distinct units of Sinomanglietia glauca (Magnoliaceae) detected by chloroplast PCR‐SSCP

Sinomanglietia glauca is a critically endangered species described from Jiangxi Province in the 1990s. Recently two populations were discovered from Yongshun County of west Hunan Province, about 450 km away from those in Jiangxi. Because of the new findings and the poor reproducibility inherent to RAPD and ISSR markers of previous studies, the population structure of this rare species was reanalyzed with chloroplast PCR‐SSCP (single‐stranded conformation polymorphism), including all of four recorded populations. The results showed that two distinct haplotypes characterized Jiangxi and Hunan populations separately, with no genetic variation occurring within regions. We postulated that this surprising pattern might result from habitat fragmentation and demographic bottlenecks during and/or after the Quaternary glaciation. On the basis of the pronounced genetic structure, two evolutionarily significant units (ESUs) were recommended for effective conservation of S. glauca.     

Patterns of genetic diversity in Great Lakes bloaters (Coregonus hoyi) with a view to future reintroduction in Lake Ontario

The originally diverse ciscoe fish fauna of the Laurentian Great Lakes has suffered many extinctions and local extirpations. Bloaters (Coregonus hoyi) are presumed extirpated from Lake Ontario and the reintroduction of this deepwater fish is under consideration. Given the demographic fluctuations of this species in the other Great Lakes and its recent intralacustrine origin, we sought to identify a genetically diverse and similar source of C. hoyi via an analysis of genetic diversity and population structure using 10 microsatellite loci. Despite well-documented demographic declines, we found no genetic evidence of bottlenecks in 12 C. hoyi samples from the four potential donor lakes (Huron, Michigan, Superior and Nipigon). By contrast, evidence of bottlenecks in historical samples of C. artedi from Lake Ontario suggested that standard genetic methods frequently used to identify population bottlenecks can only detect very severe and long-lasting demographic declines in naturally large populations. Patterns of genetic differentiation and assignment tests indicated that C. hoyi from Lake Huron and Lake Michigan, which are not differentiated, are genetically most similar to Lake Ontario ciscoes. The small available sample of deepwater ciscoes recently caught in Lake Ontario did not allow determining if these represent a small undetected C. hoyi population or a recent invasion of the deep section by C. artedi. On the basis of genetic criteria, we conclude that C. hoyi from any location within Lake Huron or Lake Michigan would be judicious sources of breeders for reintroducing C. hoyi in Lake Ontario.     

Duplication and population dynamics shape historic patterns of selection and genetic variation at the major histocompatibility complex in rodents

Genetic variation at the major histocompatibility complex (MHC) is vitally important for wildlife populations to respond to pathogen threats. As natural populations can fluctuate greatly in size, a key issue concerns how population cycles and bottlenecks that could reduce genetic diversity will influence MHC genes. Using 454 sequencing, we characterized genetic diversity at the DRB Class II locus in montane voles (Microtus montanus), a North American rodent that regularly undergoes high‐amplitude fluctuations in population size. We tested for evidence of historic balancing selection, recombination, and gene duplication to identify mechanisms maintaining allelic diversity. Counter to our expectations, we found strong evidence of purifying selection acting on the DRB locus in montane voles. We speculate that the interplay between population fluctuations and gene duplication might be responsible for the weak evidence of historic balancing selection and strong evidence of purifying selection detected. To further explore this idea, we conducted a phylogenetically controlled comparative analysis across 16 rodent species with varying demographic histories and MHC duplication events (based on the maximum number of alleles detected per individual). On the basis of phylogenetic generalized linear model‐averaging, we found evidence that the estimated number of duplicated loci was positively related to allelic diversity and, surprisingly, to the strength of purifying selection at the DRB locus. Our analyses also revealed that species that had undergone population bottlenecks had lower allelic richness than stable species. This study highlights the need to consider demographic history and genetic structure alongside patterns of natural selection to understand resulting patterns of genetic variation at the MHC.     

Introduction history overrides social factors in explaining genetic structure of females in Mediterranean mouflon

Fine‐scale spatial genetic structure of populations results from social and spatial behaviors of individuals such as sex‐biased dispersal and philopatry. However, the demographic history of a given population can override such socio‐spatial factors in shaping genetic variability when bottlenecks or founder events occurred in the population. Here, we investigated whether socio‐spatial organization determines the fine‐scale genetic structure for both sexes in a Mediterranean mouflon (Ovis gmelini musimon × Ovis sp.) population in southern France 60 years after its introduction. Based on multilocus genotypes at 16 loci of microsatellite DNA (n = 230 individuals), we identified three genetic groups for females and two for males, and concurrently defined the same number of socio‐spatial units using both GPS‐collared individuals (n = 121) and visual resightings of marked individuals (n = 378). The socio‐spatial and genetic structures did not match, indicating that the former was not the main driver of the latter for both sexes. Beyond this structural mismatch, we found significant, yet low, genetic differentiation among female socio‐spatial groups, and no genetic differentiation in males, with this suggesting female philopatry and male‐biased gene flow, respectively. Despite spatial disconnection, females from the north of the study area were genetically closer to females from the south, as indicated by the spatial analysis of the genetic variability, and this pattern was in accordance with the common genetic origin of their founders. To conclude, more than 14 generations later, genetic signatures of first introduction are not only still detectable among females, but they also represent the main factor shaping their present‐time genetic structure.     

A range-wide genetic bottleneck overwhelms contemporary landscape factors and local abundance in shaping genetic patterns of an alpine butterfly (Lepidoptera: Pieridae: Colias behrii)

Spatial and environmental heterogeneity are major factors in structuring species distributions in alpine landscapes. These landscapes have also been affected by glacial advances and retreats, causing alpine taxa to undergo range shifts and demographic changes. These nonequilibrium population dynamics have the potential to obscure the effects of environmental factors on the distribution of genetic variation. Here, we investigate how demographic change and environmental factors influence genetic variation in the alpine butterfly Colias behrii. Data from 14 microsatellite loci provide evidence of bottlenecks in all population samples. We test several alternative models of demography using approximate Bayesian computation (ABC), with the results favouring a model in which a recent bottleneck precedes rapid population growth. Applying independent calibrations to microsatellite loci and a nuclear gene, we estimate that this bottleneck affected both northern and southern populations 531–281 years ago, coinciding with a period of global cooling. Using regression approaches, we attempt to separate the effects of population structure, geographical distance and landscape on patterns of population genetic differentiation. Only 40% of the variation in F_ST is explained by these models, with geographical distance and least‐cost distance among meadow patches selected as the best predictors. Various measures of genetic diversity within populations are also decoupled from estimates of local abundance and habitat patch characteristics. Our results demonstrate that demographic change can have a disproportionate influence on genetic diversity in alpine species, contrasting with other studies that suggest landscape features control contemporary demographic processes in high‐elevation environments.