Global phylogeography of sailfish: deep evolutionary lineages with implications for fisheries management

Hydrobiologia - Since the Miocene profound climatic changes have influenced the biology and ecology of species worldwide, such as their connectivity, genetic population structure, and biogeography....     

DNA from bird-dispersed seed and wind-disseminated pollen provides insights into postglacial colonization and population genetic structure of whitebark pine (Pinus albicaulis)

Uniparentally inherited mitochondrial (mt)DNA and chloroplast (cp)DNA microsatellites (cpSSRs) were used to examine population genetic structure and biogeographic patterns of bird-dispersed seed and wind-disseminated pollen of whitebark pine (Pinus albicaulis Engelm.). Sampling was conducted from 41 populations throughout the range of the species. Analyses provide evidence for an ancestral haplotype and two derived mtDNA haplotypes with distinct regional distributions. An abrupt contact zone between mtDNA haplotypes in the Cascade Range suggests postglacial biogeographic movements. Among three cpSSR loci, 42 haplotypes were detected within 28 cpSSR sample populations that were aggregated into six regions. Analysis of molecular variance (amova) was used to determine the hierarchical genetic structure of cpSSRs. amova and population pairwise comparisons (FST ) of cpSSR, and geographical distribution of mtDNA haplotypes provide insights into historical changes in biogeography. The genetic data suggest that whitebark pine has been intimately tied to climatic change and associated glaciation, which has led to range movements facilitated by seed dispersal by Clark's nutcracker (Nucifraga columbiana Wilson). The two hypotheses proposed to explain the genetic structure are: (i) a northward expansion into Canada and the northern Cascades in the early Holocene; and (ii) historical gene flow between Idaho and the Oregon Cascades when more continuous habitat existed in Central Oregon during the late Pleistocene. Genetic structure and insights gained from historical seed movements provide a basis on which to develop recovery plans for a species that is at risk from multiple threats.     

Plio‐Pleistocene diversification and genetic population structure of an endangered lizard (the Blue Mountains water skink,Eulamprus leuraensis) in south‐eastern Australia

Aim Although climatic fluctuations occurred world‐wide during the Pleistocene, the severity of glacial and drought events – and hence their influence on animal and plant biogeography – differed among regions. Many Holarctic species were forced to warmer‐climate refugia during glacial periods, leaving the genetic signature of recent expansion and gene flow among modern‐day populations. Montane south‐eastern Australia experienced less extreme glaciation, but the effects of drier and colder climatic conditions over this period on biotic distributions, and hence on the present‐day genetic structure of animal and plant populations, are poorly known. Location South‐eastern Australia. Methods The endangered Blue Mountains water skink (Eulamprus leuraensis) is a viviparous lizard known from fewer than 40 isolated small swamps at 560–1060 m elevation in south‐eastern Australia. We conducted molecular phylogenetic, dating and population genetics analyses using the mitochondrial NADH dehydrogenase 4 (ND4) of 224 individuals of E. leuraensis sampled across the species’ distribution. Results Ancient divergences in haplotype groups between lizards from the Blue Mountains and the Newnes Plateau, and strong genetic differences, even between swamps separated by only a few kilometres, suggest that the species has persisted as a series of relatively isolated populations within its current distribution for about a million years. Presumably, habitat patches similar to current‐day swamps persisted throughout glacial–interglacial cycles in this region, allowing the development of high levels of genetic structuring within and among present‐day populations. Main conclusions Our results suggest that less extreme glacial conditions occurred in the Southern Hemisphere compared with the Northern Hemisphere, allowing cold‐adapted species (such as E. leuraensis) to persist in montane areas. However, additional studies are needed before we can assemble a comprehensive view of the impact of Pleistocene climatic variation on the phylogeography of Southern Hemisphere taxa.     

Phylogeography and species biogeography of montane Great Basin stoneflies

Sky islands are ideal systems for determining the effects of climatic oscillations on species distributions and genetic structure. Our study focused on montane stonefly populations in the Great Basin of western North America. We used niche-based distribution modelling, phylogeography and traditional species-based biogeography to test several hypotheses as follows: (i) genetic differentiation among Doroneuria baumanni populations will be independent of hydrologic connectivity (headwater model); (ii) Sky islands were colonized when habitat was more continuous and populations likely experienced multiple expansions and contractions; (iii) Colonization events were coincident with the late Pleistocene and Holocene; and (iv) Shared topography and climate history will result in concordant patterns of genetic differentiation in D. baumanni and occurrences of 32 stonefly species across the region. Overall, Φ(ST) 's and coalescent-based estimates of migration were consistent with the headwater model. Maximum likelihood and Bayesian gene trees identified three major nonoverlapping east-west clades. Distribution modelling indicated more suitable habitat in the Great Basin during the Last Glacial Maximum than at present, but none during the last interglacial period. Demographic analyses showed evidence of population expansion in one of the three major east-west clades. Intra-clade divergence times (60,000-183,000ybp) were well within the late Pleistocene while among-clade divergence times (499.000-719,000ybp) were deeper. Genetic differentiation in D. baumanni and distributions of stonefly species were significantly concordant. These results imply that climatic oscillations have played major roles in shaping the genetic structure and distributions of Great Basin stoneflies, but that divergence among clades occurred much earlier than our late Pleistocence/early Holocene predictions.     

Landscape discontinuities influence gene flow and genetic structure in a large, vagile Australian mammal, Macropus fuliginosus

Large vagile mammals typically exhibit little genetic structuring across their range, particularly when their habitat is essentially continuous. We investigated the population genetic structure of a large vagile Australian macropodid, Macropus fuliginosus, which is continuously distributed across most of southern Australia, using nine highly polymorphic nuclear microsatellite loci. Five distinct genetic units were identified across the range, four on the mainland and one on Kangaroo Island. In addition to the predicted historic Nullarbor Plain Barrier, two unexpected mainland barriers to gene flow were identified. Both were associated with landscape discontinuities (Swan River, Flinders Ranges), which appear within the dispersal capabilities of M. fuliginosus. Typical of large vagile mammals, M. fuliginosus displays high genetic diversity (with the exception of an insular population) and weak genetic structuring (within genetic units). However, the expansion of M. fuliginosus from southwestern Australia during the Pleistocene has resulted in significantly reduced genetic diversity in eastern populations. No significant sex-biased dispersal was detected, although differences in habitat, densities and climatic conditions between the eastern and western regions of the range appear to influence dispersal with the effects of isolation by distance only evident in the west. These results suggest that the biogeography of southern Australia is more complex than previously thought and reveal that seemingly minor landscape features can significantly impact genetic structuring in large vagile mammals.     

Historical perspectives on population genetics and conservation of three marine turtle species

Considerable phylogeographic and population genetic research has been conducted on marine turtles. Less attention, however, has been paid to the historical patterns and processes that have led to present patterns of genetic structure, and particularly, how these populations have responded to major climatic changes in the past. To address these questions, we analyzed previously published mitochondrial haplotype data independently for three marine turtle species, the loggerhead (Caretta caretta), hawksbill (Eretmochelys imbricata), and green turtle (Chelonia mydas). Considering all three species, we conducted analyses on a total of 657 individuals from 36 nesting beaches in the Atlantic and Mediterranean. Our results suggest that much of the contemporary genetic structure has been significantly affected by complex patterns of historical population subdivision, long-distance dispersal, and restricted geneflow. These inferences also imply that the climatic and sea level fluctuations during the Pleistocene may have had contrasting effects on genetic structure (e.g., fragmenting versus homogenizing) and on population sizes. Estimates of historical and current effective population sizes further highlight differential demographic responses across species to historical climatic cycles. Collectively, our results provide evidence for the occurrence of historical refugia through climatic cycles and complex historical metapopulation dynamics, and identify common and unique patterns of metapopulation structure across species. These historical patterns provide a basis for predictive estimates of metapopulation responses to habitat loss, population extirpation, and global climatic change.     

Range expansion and retraction along a moving contact zone has no effect on the genetic diversity of two passerine birds

Disentangling the factors shaping species distributions remains a central goal in biogeography, ecology and evolutionary biology. The extrinsic pressures that may facilitate range shifts, such as climatic factors or biotic interactions are well known. However, in contrast, the possible intrinsic factors are manifold and hard to generalize across taxa. Recently, several theoretical studies have investigated the consequences of moving range borders on genetic diversity. However, empirical studies that support or refute these theoretical predictions are scarce. Moving contact zones between parapatric sister species are suitable models to test these hypotheses. Changes in genetic diversity can be tested simultaneously along the expanding and receding edges of two species of the contact zone while accounting for intra‐specific effects (e.g. introgression). The two Old World warblers Hippolais polyglotta and H. icterina form a narrow moving contact zone, where interspecific interactions are suspected to be the main factor shaping this zone. We investigated the population genetic structure of both species along a transect ranging from the core range of the expanding H. polyglotta across the contact zone and far into the range of the receding H. icterina. The theoretical predictions of changes in genetic diversity at the range edges were tested. No gradual change in genetic diversity was detected for both the expanding and the receding range margin. Furthermore, no genetic structure was found in either species supporting the hypothesis that long distance dispersal (LDD) occurs frequently due to the high mobility of these long‐distance migrants. The results suggest that when dispersal propensity is high and accompanied by frequent LDD events, then neither an enrichment nor a depletion of alleles along moving range edges would be detected. This these species as the probability to retain genetic diversity during exogenous induced range shifts is high in such mobile species.     

A trait‐based approach to predict population genetic structure in bees

Understanding population genetic structure is key to developing predictions about species susceptibility to environmental change, such as habitat fragmentation and climate change. It has been theorized that life‐history traits may constrain some species in their dispersal and lead to greater signatures of population genetic structure. In this study, we use a quantitative comparative approach to assess if patterns of population genetic structure in bees are driven by three key species‐level life‐history traits: body size, sociality, and diet breadth. Specifically, we reviewed the current literature on bee population genetic structure, as measured by the differentiation indices Nei's G_ST, Hedrick's G′_ST, and Jost's D. We then used phylogenetic generalised linear models to estimate the correlation between the evolution of these traits and patterns of genetic differentiation. Our analyses revealed a negative and significant effect of body size on genetic structure, regardless of differentiation index utilized. For Hedrick's G′_ST and Jost's D, we also found a significant impact of sociality, where social species exhibited lower levels of differentiation than solitary species. We did not find an effect of diet specialization on population genetic structure. Overall, our results suggest that physical dispersal or other functions related to body size are among the most critical for mediating population structure for bees. We further highlight the importance of standardizing population genetic measures to more easily compare studies and to identify the most susceptible species to landscape and climatic changes.     

A watershed study on genetic diversity: Phylogenetic analysis of the Platypleura plumosa (Hemiptera: Cicadidae) complex reveals catchment-specific lineages

Historical biogeography studies have at their disposal a small suite of vicariance models to explain genetic differentiation within and between species. One of these processes involves the role of river catchments and their associated watersheds, in driving diversification and is applicable to both aquatic and terrestrial organisms. Although the idea of catchments structuring the genetic history of aquatic organisms is reasonably well understood, their effect on terrestrial organisms has largely been overlooked, with relevant studies being limited in scope. South Africa presents a perfect test-bed for elucidating this mechanism of diversification due to its rich biodiversity, range of climatic environments and many large river catchments. Here we use the cicadas of the Platypleura plumosa complex to highlight the importance of catchments and their associated watersheds in driving diversification of terrestrial invertebrates that lack an aquatic life-stage. Population structure was found to correspond to primary and in some cases secondary catchments; highlighting the need to include information on catchment structure when formulating hypotheses of population diversification. Recognizing that climate change in the near future is likely to alter the environment, and particularly precipitation patterns, insight into recent patterns of population change related to catchments may be useful in a conservation context.     

Genetic and phylogenetic consequences of island biogeography

Abstract.— Island biogeography theory predicts that the number of species on an island should increase with island size and decrease with island distance to the mainland. These predictions are generally well supported in comparative and experimental studies. These ecological, equilibrium predictions arise as a result of colonization and extinction processes. Because colonization and extinction are also important processes in evolution, we develop methods to test evolutionary predictions of island biogeography. We derive a population genetic model of island biogeography that incorporates island colonization, migration of individuals from the mainland, and extinction of island populations. The model provides a means of estimating the rates of migration and extinction from population genetic data. This model predicts that within an island population the distribution of genetic divergences with respect to the mainland source population should be bimodal, with much of the divergence dating to the colonization event. Across islands, this model predicts that populations on large islands should be on average more genetically divergent from mainland source populations than those on small islands. Likewise, populations on distant islands should be more divergent than those on close islands. Published observations of a larger proportion of endemic species on large and distant islands support these predictions.     

Montane refugia predict population genetic structure in the Large‐blotched Ensatina salamander

Understanding the biotic consequences of Pleistocene range shifts and fragmentation remains a fundamental goal in historical biogeography and evolutionary biology. Here, we combine species distribution models (SDM) from the present and two late Quaternary time periods with multilocus genetic data (mitochondrial DNA and microsatellites) to evaluate the effect of climate‐induced habitat shifts on population genetic structure in the Large‐blotched Ensatina (Ensatina eschscholtzii klauberi), a plethodontid salamander endemic to middle and high‐elevation conifer forest in the Transverse and Peninsular Ranges of southern California and northern Baja California. A composite SDM representing the range through time predicts two disjunct refugia, one in southern California encompassing the core of the species range and the other in the Sierra San Pedro Mártir of northern Baja California at the southern limit of the species range. Based on our spatial model, we would expect a pattern of high connectivity among populations within the northern refugium and, conversely, a pattern of isolation due to long‐term persistence of the Sierra San Pedro Mártir population. Our genetic results are consistent with these predictions based on the hypothetical refugia in that (i) historical measures of population connectivity among stable areas are correlated with gene flow estimates; and (ii) there is strong geographical structure between separate refugia. These results provide evidence for the role of recent climatic change in shaping patterns of population persistence and connectivity within the Transverse and Peninsular Ranges, an evolutionary hotspot.     

Population fragmentation causes randomly fixed genotypes in populations of <Emphasis Type="Italic">Arabidopsis kamchatica</Emphasis> in the Japanese Archipelago

Populations of arctic alpine plants likely disappeared and re-colonised several times at the southern edge of their distributions during glacial and interglacial cycles throughout the Quaternary. Range shift and population fragmentation after a glacial period would affect the genetic structure of such plants in southernmost populations. We aimed to elucidate how climatic oscillations influenced the population subsistence of alpine plants in the Japanese Archipelago as one of the southernmost populations, by inferring the genetic structure of Arabidopsis kamchatica subsp. kamchatica and the intraspecific littoral taxon, subsp. kawasakiana. We identified genotypes based on the haplotypes of five nuclear genes and two chloroplast DNA spacers for 164 individuals from 24 populations. Most populations harboured only one private genotype, whereas few polymorphisms were found in each population. Two genetic genealogies were found, suggesting that northern Japanese populations of alpine subsp. kamchatica, subsp. kawasakiana and the northerly subsp. kamchatica in eastern Russia and Alaska clustered and differentiated from populations in central Honshu, western Japan and Taiwan. During climatic oscillations, the genetic structure of extant southernmost populations would have been shaped by strong genetic drift under population fragmentation and randomly fixed to a single genotype among their ancestral polymorphisms.     

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.     

Impact of Quaternary climatic changes and interspecific competition on the demographic history of a highly mobile generalist carnivore, the coyote

Recurrent cycles of climatic change during the Quaternary period have dramatically affected the population genetic structure of many species. We reconstruct the recent demographic history of the coyote (Canis latrans) through the use of Bayesian techniques to examine the effects of Late Quaternary climatic perturbations on the genetic structure of a highly mobile generalist species. Our analysis reveals a lack of phylogeographic structure throughout the range but past population size changes correlated with climatic changes. We conclude that even generalist carnivorous species are very susceptible to environmental changes associated with climatic perturbations. This effect may be enhanced in coyotes by interspecific competition with larger carnivores.     

Revisiting Amazonian phylogeography: insights into diversification hypotheses and novel perspectives

The Amazon Basin harbors one of the richest biotas on Earth, such that a number of diversification hypotheses have been formulated to explain patterns of Amazonian biodiversity and biogeography. For nearly two decades, phylogeographic approaches have been applied to better understand the underlying causes of genetic differentiation and geographic structure among Amazonian organisms. Although this research program has made progress in elucidating several aspects of species diversification in the region, recent methodological and theoretical developments in the discipline of phylogeography will provide new perspectives through more robust hypothesis testing. Herein, we outline central aspects of Amazonian geology and landscape evolution as well as climate and vegetation dynamics through the Neogene and Quaternary to contextualize the historical settings considered by major hypotheses of diversification. We address each of these hypotheses by reviewing key phylogeographic papers and by expanding their respective predictions. We also propose future directions for devising and testing hypotheses. Specifically, combining the exploratory power of phylogeography with the statistical rigor of coalescent methods will greatly expand analytical inferences on the evolutionary history of Amazonian biota. Incorporation of non-genetic data from Earth science disciplines into the phylogeographic approach is key to a better understanding of the influence of climatic and geophysical events on patterns of Amazonian biodiversity and biogeography. In addition, achieving such an integrative enterprise must involve overcoming issues such as limited geographic and taxonomic sampling. These future challenges likely will be accomplished by a combination of extensive collaborative research and incentives for conducting basic inventories.     

Extensive population expansion of Pedicularis longiflora (Orobanchaceae) on the Qinghai‐Tibetan Plateau and its correlation with the Quaternary climate change

The Qinghai‐Tibetan Plateau (QTP) is thought to be more strongly affected by the Quaternary glaciations than most other regions of the same latitude. It would be of great interest to investigate the population genetic structure of organisms distributed on the platform and its correlation with the Quaternary climatic oscillations. Here we used the chloroplast (cp)DNA trnT‐trnF sequence to study genetic variation and phylogeography of Pedicularis longiflora, an alpine herb with extensive distribution on the QTP. Based on a range‐wide sampling comprising 41 populations and 910 individuals, we detected 30 cpDNA haplotypes that were divided into five clades by phylogenetic and network analyses and a strong phylogeographical structure. All haplotypes but one in the three basal clades occur exclusively in the southeast QTP, whereas haplotypes in the young clade V occupy almost the whole species range. In particular, the young haplotype H18 occurs in 420 individuals, even at a frequency of 100% in some QTP platform populations and the Altai population. The haplotype distribution pattern, together with molecular clock estimation and mismatch distribution analysis, suggests that the southeast QTP was either a refuge for P. longiflora during the Quaternary climatic change or is the place of origin of the species. The present wide distribution of the species on the QTP platform has resulted from recent population expansions which could be dated back to 120 000–17 000 years ago, a period mostly before the last glacial maximum. The possible relationships among geographic genetic structure, climatic change and species diversification in Pedicularis are also discussed.     

Differential effects of historical migration,glaciations and human impact on the genetic structure and diversity of the mountain pasture weed Veratrum album L.

Aim Today’s genetic population structure and diversity of species can be understood as the result of range expansion from the area of origin, past climatic oscillations and contemporary processes. We examined the relative importance of these factors in Veratrum album L., a toxic weed of mountain grasslands. Location Continental Europe. Methods Forty populations from the Asian border (Urals and Caucasus) to Portugal were studied using amplified fragment length polymorphisms (AFLPs) combined with selected plant and population measures. The data were analysed with phylogenetic, population genetic and regression methods inferring both genetic structure and diversity from geographic and ecological factors. Results Fragment frequency clines together with genetic distance clustering and principal coordinates analysis indicated an east–west direction in the genetic structure of V. album, suggesting ancient migration into Europe from a proposed Asian origin. However, the strong geographic pattern in the genetic structure, pronounced isolation by distance (R² = 0.74) and moderate overall population differentiation (F_ST = 0.13) suggests high historical gene flow, possibly during glacials, and vicariance into mountainous regions during interglacials. Occurrence of V. album during the last glaciation in several areas along the periphery of the Alps and recolonization of this mountain range from both eastern and central–western areas was indicated. Genetic diversity was highest in central Europe, a pattern that did not agree with the expectations from east–west migration into Europe. Furthermore, managed habitats showed higher levels of genetic diversity compared to unmanaged habitats. Stepwise linear regression determined shoot density and soil phosphorus as the main predictors of within‐population genetic diversity (R² = 0.40). Main conclusions Our results showed that V. album retained genetic imprints of historical range expansion into Europe, although this was alleviated by the influence of climatic oscillations and contemporary processes. For example, genetic population structure was strongly affected by post‐glacial vicariance while patterns of genetic diversity seemed mainly to be influenced by human land use. Our findings highlight the importance of applying a synthetic approach, testing the influence of both historical and contemporary processes on genetic structure and diversity in order to understand complex phylogeographic patterns. This may especially apply to widespread species, such as weeds. Implications of our findings for biological control are briefly discussed.     

Genetic and paleomodelling evidence of the population expansion of the cattle egret Bubulcus ibis in Africa during the climatic oscillations of the Late Pleistocene

Increasing aridity during glacial periods produced the retraction of forests and the expansion of arid and semi‐arid environments in Africa, with consequences for birds. Cattle egret Bubulcus ibis is a dispersive species that prefers semi‐arid environments and requires proximity to bodies of water. We expected that climatic oscillations led to the expansion of the range of the cattle egret during arid periods, such as the Last Maximum Glacial (LGM) and contraction of distribution during the Last Interglacial (LIG) period, resulting in contact of populations previously isolated. We investigated this hypothesis by evaluating the genetic structure and population history of 15 cattle egret breeding colonies located in west and South Africa using the mitochondrial DNA (mtDNA) control region, mtDNA ATPase 8 and 6, and an intron of nuclear gene transforming growth factor‐beta 2. Occurrence data and bioclimatic information were used to generate ecological niche models of three periods (present, LGM and LIG). We used the genetic and paleomodelling data to assess the responses of the cattle egret from Africa to the climatic oscillations during the late Pleistocene. Genetic data revealed low levels of genetic differentiation, signs of isolation‐by‐distance, as well as recent increases in effective population size that started during the LGM. The observed low genetic structure may be explained by recent colonization events due to the demographic expansion following the last glacial period and by dispersal capacity of this species. The paleomodels corroborated the expansion during the LGM, and a more restricted potential distribution during the LIG. Our findinds supports the hypothesis that the species range of the cattle egret expanded during arid periods and contracted during wet periods.     

Genetic structure of Partamona helleri (Apidae, Meliponini) from Neotropical Atlantic rainforest

Stingless bees of the genus Partamona are distributed from southern Mexico to southern Brazil. This genus has been subject to different approaches to solve questions concerning general biology, taxonomy, systematics and biogeography, but population studies applying molecular techniques are inexistent. We analyzed the genetic structure of P. helleri across its geographic distribution along the coastal Atlantic tropical rainforest in Brazil. Ten mtDNA haplotypes were observed in 47 colonies of P. helleri of which some were exclusive and others shared among geographic sub-groups. Statistical analysis showed high genetic differentiation between geographic areas sampled. Fragmentation of the Atlantic forest during Pleistocene glaciations is discussed as a possible cause of the present haplotype distribution and frequency.     

Genetic structure in the European endemic seabird,Phalacrocorax aristotelis,shaped by a complex interaction of historical and contemporary,physical and nonphysical drivers

Geographically separated populations tend to be less connected by gene flow, as a result of physical or nonphysical barriers preventing dispersal, and this can lead to genetic structure. In this context, highly mobile organisms such as seabirds are interesting because the small effect of physical barriers means nonphysical ones may be relatively more important. Here, we use microsatellite and mitochondrial data to explore the genetic structure and phylogeography of Atlantic and Mediterranean populations of a European endemic seabird, the European shag, Phalacrocorax aristotelis, and identify the primary drivers of their diversification. Analyses of mitochondrial markers revealed three phylogenetic lineages grouping the North Atlantic, Spanish/Corsican and eastern Mediterranean populations, apparently arising from fragmentation during the Pleistocene followed by range expansion. These traces of historical fragmentation were also evident in the genetic structure estimated by microsatellite markers, despite significant contemporary gene flow among adjacent populations. Stronger genetic structure, probably promoted by landscape, philopatry and local adaptation, was found among distant populations and those separated by physical and ecological barriers. This study highlights the enduring effect of Pleistocene climatic changes on shag populations, especially within the Mediterranean Basin, and suggests a role for cryptic northern refugia, as well as known southern refugia, on the genetic structure of European seabirds. Finally, it outlines how contemporary ecological barriers and behavioural traits may maintain population divergence, despite long‐distance dispersal triggered by extreme environmental conditions (e.g. population crashes).