Phylogeography and the origins of range disjunctions in a north temperate fish,the pygmy whitefish (Prosopium coulterii), inferred from mitochondrial and nuclear DNA sequence analysis

Aim To investigate the degree of phylogeographical divergence within pygmy whitefish (Prosopium coulterii) and to test hypotheses concerning the origin of disjunct populations within North America. Location North America from western Alaska to Lake Superior. Methods Mitochondrial (ATPase subunit VI) and nuclear (ITS‐1, ITS‐2) DNA sequence variation was assessed across the species’ North American range to test for the existence of distinct phylogeographical groupings of pygmy whitefish associated with known glacial refugia. Coalescent simulations of the mitochondrial DNA (mtDNA) data were used to test hypotheses of population structure. Results This species is composed of two monophyletic mitochondrial clades across its North American range. The two mtDNA clades differed by an average 3.3% nucleotide sequence divergence. These clades were also distinguished by ITS‐2, but the relationships among lineages were not resolved by the ITS‐1 analysis. Coalescent analyses rejected the null hypothesis of the current disjunct distributions being a result of fragmentation of a single widespread ancestral lineage across a variety of effective population sizes and divergence times. Main conclusions The current range disjunctions of pygmy whitefish in North America probably resulted from isolation, genetic divergence, and selective dispersal from at least two major Pleistocene glacial refugia: Beringia and Cascadia. More recent isolation and dispersal from an upper Mississippi refugium is suggested by relationships within one of the clades and by distributional evidence from co‐distributed species. The Beringian and Cascadian refugia have played major roles in the zoogeography of Nearctic temperate aquatics, but the roles of smaller refugia appear more variable among other species.     

Pleistocene persistence and expansion in tarantulas on the Colorado Plateau and the effects of missing data on phylogeographical inferences from RADseq

Montane species endemic to the “sky islands” of the North American southwest were significantly impacted by changing climates during the Pleistocene. We combined mitochondrial and genomic data with species distribution modelling to determine whether Aphonopelma marxi, a large tarantula from the nearby Colorado Plateau, was similarly impacted by glacial climates. Genetic analyses revealed that the species comprises three main clades that diverged in the Pleistocene. A clade distributed along the Mogollon Rim appears to have persisted in place during glacial conditions, whereas the other two clades probably colonized central and northeastern portions of the species' range from refugia in canyons. Climate models support this hypothesis for the Mogollon Rim, but late glacial climate data appear too coarse to detect suitable areas in canyons. Locations of canyon refugia could not be inferred from genomic analyses due to missing data, encouraging us to explore the effect of missing loci in phylogeographical inferences using RADseq. Results from analyses with varying amounts of missing data suggest that samples with large amounts of missing data can still improve inferences, and the specific loci that are missing matters more than the number of missing loci. This study highlights the profound impact of Pleistocene climates on tarantulas endemic to the Colorado Plateau, as well as the mixed nature of the region's fauna. Some animals recently colonized from nearby deserts as glacial climates receded, whereas others, like tarantulas, appear to have persisted on the Mogollon Rim and in refugia associated with the region's famous river‐cut canyons.     

Are glacial refugia hotspots of speciation and cytonuclear discordances? Answers from the genomic phylogeography of Spanish common frogs

Subdivided Pleistocene glacial refugia, best known as “refugia within refugia”, provided opportunities for diverging populations to evolve into incipient species and/or to hybridize and merge following range shifts tracking the climatic fluctuations, potentially promoting extensive cytonuclear discordances and “ghost” mtDNA lineages. Here, we tested which of these opposing evolutionary outcomes prevails in northern Iberian areas hosting multiple historical refugia of common frogs (Rana cf. temporaria), based on a genomic phylogeography approach (mtDNA barcoding and RAD‐sequencing). We found evidence for both incipient speciation events and massive cytonuclear discordances. On the one hand, populations from northwestern Spain (Galicia and Asturias, assigned to the regional endemic R. parvipalmata), are deeply‐diverged at mitochondrial and nuclear genomes (~4 My of independent evolution), and barely admix with northeastern populations (assigned to R. temporaria sensu stricto) across a narrow hybrid zone (~25 km) located in the Cantabrian Mountains, suggesting that they represent distinct species. On the other hand, the most divergent mtDNA clade, widespread in Cantabria and the Basque country, shares its nuclear genome with other R. temporaria s. s. lineages. Patterns of population expansions and isolation‐by‐distance among these populations are consistent with past mitochondrial capture and/or drift in generating and maintaining this ghost mitochondrial lineage. This remarkable case study emphasizes the complex evolutionary history that shaped the present genetic diversity of refugial populations, and stresses the need to revisit their phylogeography by genomic approaches, in order to make informed taxonomic inferences.     

Phylogeographical lineages of Arctic grayling (Thymallus arcticus) in North America: divergence, origins and affinities with Eurasian Thymallus

The number and location of Arctic glacial refugia utilized by taxa during the Pleistocene are continuing uncertainties in Holarctic phylogeography. Arctic grayling (Thymallus arcticus) are widely distributed in freshwaters from the eastern side of Hudson Bay (Canada) west to central Asia. We studied mitochondrial DNA (mtDNA) and microsatellite DNA variation in North American T. arcticus to test for genetic signatures of survival in, and postglacial dispersal from, multiple glacial refugia, and to assess their evolutionary affinities with Eurasian Thymallus. In samples from 32 localities, we resolved 12 mtDNA haplotypes belonging to three assemblages that differed from each other in sequence by between 0.75 and 2.13%: a 'South Beringia' lineage found from western Alaska to northern British Columbia, Canada; a 'North Beringia' lineage found on the north slope of Alaska, the lower Mackenzie River, and to eastern Saskatchewan; and a 'Nahanni' lineage confined to the Nahanni River area of the upper Mackenzie River drainage. Sequence analysis of a portion of the control region indicated monophyly of all North American T. arcticus and their probable origin from eastern Siberian T. arcticus at least 3 Mya. Arctic grayling sampled from 25 localities displayed low allelic diversity and expected heterozygosity (H(E)) across five microsatellite loci (means of 2.1 alleles and 0.27 H(E), respectively) and there were declines in these measures of genetic diversity with distance eastward from the lower Yukon River Valley. Assemblages defined by mtDNA divergences were less apparent at microsatellite loci, but again the Nahanni lineage was the most distinctive. Analysis of molecular variance indicated that between 24% (microsatellite DNA) and 81% (mtDNA) of the variance was attributable to differences among South Beringia, North Beringia and Nahanni lineages. Our data suggest that extant North American Arctic grayling are more diverse phylogeographically than previously suspected and that they consist of at least three major lineages that originated in distinct Pleistocene glacial refugia. T. arcticus probably originated and dispersed from Eurasia to North America in the late to mid-Pliocene, but our data also suggest more recent (mid-late Pleistocene) interactions between lineages across Beringia.     

A parapatric propensity for breeding precludes the completion of speciation in common teal (Anas crecca, sensu lato)

Speciation is a process in which genetic drift and selection cause divergence over time. However, there is no rule dictating the time required for speciation, and even low levels of gene flow hinder divergence, so that taxa may be poised at the threshold of speciation for long periods of evolutionary time. We sequenced mitochondrial DNA (mtDNA) and eight nuclear introns (nuDNA) to estimate genomic levels of differentiation and gene flow between the Eurasian common teal (Anas crecca crecca) and the North American green‐winged teal (Anas crecca carolinensis). These ducks come into contact in Beringia (north‐eastern Asia and north‐western North America) and have probably done so, perhaps cyclically, since the Pliocene–Pleistocene transition, ～2.6 Ma, when they apparently began diverging. They have diagnosable differences in male plumage and are 6.9% divergent in the mtDNA control region, with only 1 of 58 crecca and 2 of 86 carolinensis having haplotypes grouping with the other. Two nuclear loci were likewise strongly structured between these teal (Φ_st ≥ 0.35), but six loci were undifferentiated or only weakly structured (Φ_st = 0.0–0.06). Gene flow between crecca and carolinensis was ～1 individual per generation in both directions in mtDNA, but was asymmetrical in nuDNA, with ～1 and ～20 individuals per generation immigrating into crecca and carolinensis, respectively. This study illustrates that species delimitation using a single marker oversimplifies the complexity of the speciation process, and it suggests that even with divergent selection, moderate levels of gene flow may stall the speciation process short of completion.     

Multilocus phylogeography and population structure of common eiders breeding in North America and Scandinavia

Aim Glacial refugia during the Pleistocene had major impacts on the levels and spatial apportionment of genetic diversity of species in northern latitude ecosystems. We characterized patterns of population subdivision, and tested hypotheses associated with locations of potential Pleistocene refugia and the relative contribution of these refugia to the post‐glacial colonization of North America and Scandinavia by common eiders (Somateria mollissima). Specifically, we evaluated localities hypothesized as ice‐free areas or glacial refugia for other Arctic vertebrates, including Beringia, the High Arctic Canadian Archipelago, Newfoundland Bank, Spitsbergen Bank and north‐west Norway. Location Alaska, Canada, Norway and Sweden. Methods Molecular data from 12 microsatellite loci, the mitochondrial DNA (mtDNA) control region, and two nuclear introns were collected and analysed for 15 populations of common eiders (n = 716) breeding throughout North America and Scandinavia. Population genetic structure, historical population fluctuations and gene flow were inferred using F‐statistics, analyses of molecular variance, and multilocus coalescent analyses. Results Significant inter‐population variation in allelic and haplotypic frequencies were observed (nuclear DNA F_ST = 0.004–0.290; mtDNA Φ_ST = 0.051–0.927). Whereas spatial differentiation in nuclear genes was concordant with subspecific designations, geographic proximity was more predictive of inter‐population variance in mitochondrial DNA haplotype frequency. Inferences of historical population demography were consistent with restriction of common eiders to four geographic areas during the Last Glacial Maximum: Belcher Islands, Newfoundland Bank, northern Alaska and Svalbard. Three of these areas coincide with previously identified glacial refugia: Newfoundland Bank, Beringia and Spitsbergen Bank. Gene‐flow and clustering analyses indicated that the Beringian refugium contributed little to common eider post‐glacial colonization of North America, whereas Canadian, Scandinavian and southern Alaskan post‐glacial colonization is likely to have occurred in a stepwise fashion from the same glacial refugium. Main conclusions Concordance of proposed glacial refugia used by common eiders and other Arctic species indicates that Arctic and subarctic refugia were important reservoirs of genetic diversity during the Pleistocene. Furthermore, suture zones identified at MacKenzie River, western Alaska/Aleutians and Scandinavia coincide with those identified for other Arctic vertebrates, suggesting that these regions were strong geographic barriers limiting dispersal from Pleistocene refugia.     

Genomic divergence between Spanish Littorina saxatilis ecotypes unravels limited admixture and extensive parallelism associated with population history

The rough periwinkle, Littorina saxatilis, is a model system for studying parallel ecological speciation in microparapatry. Phenotypically parallel wave‐adapted and crab‐adapted ecotypes that hybridize within the middle shore are replicated along the northwestern coast of Spain and have likely arisen from two separate glacial refugia. We tested whether greater geographic separation corresponding to reduced opportunity for contemporary or historical gene flow between parallel ecotypes resulted in less parallel genomic divergence. We sequenced double‐digested restriction‐associated DNA (ddRAD) libraries from individual snails from upper, mid, and low intertidal levels of three separate sites colonized from two separate refugia. Outlier analysis of 4256 SNP markers identified 34.4% sharing of divergent loci between two geographically close sites; however, these sites each shared only 9.9%–15.1% of their divergent loci with a third more‐distant site. STRUCTURE analysis revealed that genotypes from only three of 166 phenotypically intermediate mid‐shore individuals appeared to result from recent hybridization, suggesting that hybrids cannot be reliably identified using shell traits. Hierarchical AMOVA indicated that the primary source of genomic differentiation was geographic separation, but also revealed greater similarity of the same ecotype across the two geographically close sites than previously estimated with dominant markers. These results from a model system for ecological speciation suggest that genomic parallelism is affected by the opportunity for historical or contemporary gene flow between populations.     

A mirage of cryptic species: Genomics uncover striking mitonuclear discordance in the butterfly Thymelicus sylvestris

Mitochondrial DNA (mtDNA) sequencing has led to an unprecedented rise in the identification of cryptic species. However, it is widely acknowledged that nuclear DNA (nuDNA) sequence data are also necessary to properly define species boundaries. Next generation sequencing techniques provide a wealth of nuclear genomic data, which can be used to ascertain both the evolutionary history and taxonomic status of putative cryptic species. Here, we focus on the intriguing case of the butterfly Thymelicus sylvestris (Lepidoptera: Hesperiidae). We identified six deeply diverged mitochondrial lineages; three distributed all across Europe and found in sympatry, suggesting a potential case of cryptic species. We then sequenced these six lineages using double‐digest restriction‐site associated DNA sequencing (ddRADseq). Nuclear genomic loci contradicted mtDNA patterns and genotypes generally clustered according to geography, i.e., a pattern expected under the assumption of postglacial recolonization from different refugia. Further analyses indicated that this strong mtDNA/nuDNA discrepancy cannot be explained by incomplete lineage sorting, sex‐biased asymmetries, NUMTs, natural selection, introgression or Wolbachia‐mediated genetic sweeps. We suggest that this mitonuclear discordance was caused by long periods of geographic isolation followed by range expansions, homogenizing the nuclear but not the mitochondrial genome. These results highlight T. sylvestris as a potential case of multiple despeciation and/or lineage fusion events. We finally argue, since mtDNA and nuDNA do not necessarily follow the same mechanisms of evolution, their respective evolutionary history reflects complementary aspects of past demographic and biogeographic events.     

HOT SPOTS OF GENETIC DIVERSITY DESCENDED FROM MULTIPLE PLEISTOCENE REFUGIA IN AN ALPINE UNGULATE

Species that inhabit naturally fragmented environments are expected to be spatially structured and exhibit reduced genetic diversity at the periphery of their range. Patterns of differentiation may also reflect historical processes such as recolonization from glacial refugia. We examined the relative importance of these factors in shaping the spatial patterns of genetic differentiation across the range of an alpine specialist, the North American mountain goat (Oreamnos americanus). Contrary to fossil evidence that suggests a single southern refugium, we detected evidence for additional refugia in northern British Columbia and the Alaskan coast using both mitochondrial and microsatellite DNA. A core area of elevated genetic diversity characterized both regions, and molecular dating suggested a recent Pleistocene split was followed by demographic expansion. Across their range, mountain goats were highly genetically structured and displayed the expected pattern of declining diversity toward the periphery. Gene flow was high within contiguous mountain ranges, but cross‐assignments paradoxically suggest that long‐distance contemporary dispersal movements are not uncommon. These results improve our understanding of how historical vicariance and contemporary fragmentation influence population differentiation, and have implications for conserving the adaptive potential of alpine populations and habitat.     

Dispersal differences predict population genetic structure in Mormon crickets

Research investigating the geographical context of speciation has primarily focused on abiotic factors such as the role of Pleistocene glacial cycles, or geotectonic events. Few study systems allow a direct comparison of how biological differences, such as dispersal behaviour, affect population genetic structure of organisms that were subdivided during the Pleistocene. Mormon crickets exist in solitary and gregarious 'phases', which broadly correspond with an east-west mtDNA division across the Rocky Mountains. Gregarious individuals form bands that can move up to 2 km daily. This study assessed whether population genetic structure results mainly from deep Pleistocene vicariance or if we can also detect more recent genetic patterns due to phase and dispersal differences superimposed on the older, deeper divisions. We found that separation in refugia was a more important influence on genetic divergence than phase, with the Rockies acting as a barrier that separated Mormon cricket populations into eastern and western refugia during Pleistocene glacial cycles. However, patterns of isolation by distance differ between eastern and western clades for both mitochondrial and nuclear DNA, with greater divergence within the eastern, solitary clade. An mtDNA haplotype mismatch distribution is compatible with historical population expansion in the western clade but not in the eastern clade. A persistent (and possibly sex-biased) difference in dispersal ability has most likely influenced the greater population genetic structure seen in the eastern clade, emphasizing the importance of the interaction of Quaternary climate fluctuations and geography with biotic factors in producing the patterns of genetic subdivision observed today.     

Population genomic evidence for plant glacial survival in Scandinavia

Quaternary glaciations have played a major role in shaping the genetic diversity and distribution of plant species. Strong palaeoecological and genetic evidence supports a postglacial recolonization of most plant species to northern Europe from southern, eastern and even western glacial refugia. Although highly controversial, the existence of small in situ glacial refugia in northern Europe has recently gained molecular support. We used genomic analyses to examine the phylogeography of a species that is critical in this debate. Carex scirpoidea Michx subsp. scirpoidea is a dioecious, amphi‐Atlantic arctic–alpine sedge that is widely distributed in North America, but absent from most of Eurasia, apart from three extremely disjunct populations in Norway, all well within the limits of the Weichselian ice sheet. Range‐wide population sampling and variation at 5,307 single nucleotide polymorphisms show that the three Norwegian populations comprise unique evolutionary lineages divergent from Greenland with high between‐population divergence. The Norwegian populations have low within‐population genetic diversity consistent with having experienced genetic bottlenecks in glacial refugia, and host private alleles that probably accumulated in long‐term isolated populations. Demographic analyses support a single, pre‐Weichselian colonization into Norway from East Greenland, and subsequent divergence of the three populations in separate refugia. Other refugial areas are identified in North‐east Greenland, Minnesota/Michigan, Colorado and Alaska. Admixed populations in British Columbia and West Greenland indicate postglacial contact. Taken together, evidence from this study strongly indicates in situ glacial survival in Scandinavia.     

Phylogeography of <Emphasis Type="Italic">Leiopelma hochstetteri</Emphasis> reveals strong genetic structure and suggests new conservation priorities

Leiopelma hochstetteri belongs to a singular anuran lineage endemic to New Zealand that diverged from other frogs about 200 million years ago. The species now is reduced to a series of isolated populations in the northern half of the North Island and Great Barrier Island. We have used mitochondrial and nuclear sequence data to examine the genetic affinities of extant populations of L. hochstetteri. Phylogenetic reconstructions reveal that populations are highly structured. Each of the geographically isolated populations harbours independent mtDNA lineages as well as different degrees of nuDNA differentiation. Moreover, molecular dating reveals that this structure originated in the early Pleistocene. This pattern of genetic structure likely results from unfavourable climatic conditions during the Pleistocene combined with the low dispersal ability of the species. Isolated populations in forested refugia existed even in the southern part of the distribution of the species during glacial cycles. Previously published variation in chromosome numbers and isozyme data are consistent with the new evidence. We identified 13 evolutionary significant units (ESUs) that should serve as the focus for future management and conservation of this species.     

GLOBAL SURVEY OF MITOCHONDRIAL DNA SEQUENCES IN THE THREESPINE STICKLEBACK: EVIDENCE FOR RECENT MIGRATIONS

Phylogenetic analyses of mitochondrial DNA (mtDNA) sequences were used to assess the matriarchal genetic structure of the threespine stickleback, Gasterosteus aculeatus. A 747 base-pair (bp) fragment of the cytochrome b was sequenced from 36 individuals collected from 25 localities in Europe, North America, and Japan. Two major divergent clades were revealed: one widespread in Japan but with representatives in some Alaskan and British Columbian lakes and the other common in Europe and North America. A simple diagnostic test using the polymerase chain reaction (PCR) and a restriction enzyme was used to assay additional individuals, confirming the absence of the Japanese clade in the Atlantic basin. Geographic distribution of mtDNA variation suggests (1) a recent origin of the Atlantic populations, and (2) support for previous hypotheses about the existence of Pleistocene refugia for freshwater fishes in Alaska and British Columbia. Silent substitution rates were used to date the colonization of the Atlantic at 90,000 to 260,000 yr before present, which conflicts with earlier dates implied by the fossil record. The recent replacement of Atlantic mitochondrial lineages suggested by our data may be explained by severe reduction or extinction of northern Atlantic populations during the Pleistocene, followed by a recent reinvasion from the Pacific. With a global perspective of the distribution of genetic variation as a framework, meaningful comparisons at a smaller geographical scale will now be possible.     

Surviving glacial ages within the Biotic Gap: phylogeography of the New Zealand cicada Maoricicada campbelli

Aim New Zealand is an ideal location in which to investigate the roles of landscape and climate change on speciation and biogeography. An earlier study of the widespread endemic cicada Maoricicada campbelli (Myers) found two phylogeographically distinguishable major clades – northern South Island plus North Island (northern‐SI + NI) and Otago. These two clades appeared to have diverged on either side of an area of the South Island known as the Biotic Gap. We sampled more intensively to test competing theories for this divergence. We aimed to discover if M. campbelli had survived within the Biotic Gap during recent glacial maxima, and if predicted areas of secondary contact between the two major clades existed. Location New Zealand. Methods We analysed mitochondrial DNA sequences (1520 bp; 212 individuals; 91 populations) using phylogenetic (maximum likelihood, Bayesian), population genetic (analysis of molecular variance) and molecular dating methods (Bayesian relaxed clock with improved priors). Results We found strong geographical structuring of genetic variation. Our dating analyses suggest that M. campbelli originated 1.83–2.58 Ma, and split into the two major clades 1.45–2.09 Ma. The main subclades in the northern‐SI + NI clade arose almost simultaneously at 0.69–1.03 Ma. Most subclades are supported by long internal branches and began to diversify 0.40–0.78 Ma. We found four narrow areas of secondary contact between the two major clades. We also found a difference between calling songs of the Otago vs. northern‐SI + NI clades. Main conclusions Phylogeographical patterns within M. campbelli indicate an early Pleistocene split into two major clades, followed by late Pleistocene range expansion and in situ population differentiation of subclades. The northern‐SI + NI clade diversified so rapidly that the main subclade relationships cannot be resolved, and we now have little evidence for a disjunction across the Biotic Gap. Structure within the main subclades indicates rapid divergence after a common bottlenecking event, perhaps attributable to an extremely cold glacial maximum at c. 0.43 Ma. Clade structure and dating analyses indicate that M. campbelli survived in many refugia during recent glacial maxima, including within the Biotic Gap. The narrow overlap between the two major clades is attributed to recent contact during the current interglacial and slow gene diffusion. The two major clades appear to be in the early stages of speciation based on genetic and behavioural differences.     

The distribution of divergent mitochondrial DNA lineages of threespine stickleback (Gasterosteus aculeatus) in the northeastern Pacific Basin: post-glacial dispersal and lake accessibility

Aim This study furthers the documentation of the geographical distribution of two divergent (c. 3%) mitochondrial DNA clades in the threespine stickleback (Gasterosteus aculeatus) and tests the hypotheses that the northeastern Pacific distribution has been influenced by post‐glacial colonization and lake elevation and that clade identity is associated with certain morphological attributes such as reduction in body armour. Location Lakes and nearshore marine environments of the eastern Pacific Basin from southcentral Alaska to southeastern British Columbia, (BC) Canada. Methods Restriction enzyme analysis of polymerase chain reaction‐amplified mitochondrial DNA fragments (cytochrome b) from a total of 45 new populations combined with existing data for a further 45 populations. Lake elevation data were collected for 78 localities and tested for an association with mtDNA clade by contingency table analyses. Morphological data were collected on sticklebacks from eight samples representing four lake‐stream systems and tested for differentiation among populations with different mtDNA clade identities using analyses of variance. Results We extend the known distribution of the haplotypes diagnostic of the Trans‐North Pacific Clade (TNPC) southward to mid‐Vancouver Island and, for the first time, on mainland BC, in other island populations far from putative refugia, and in nearby anadromous populations. A morphological analysis indicated that the mainland population with the TNPC was not characterized by reduced spine or lateral plate (‘armour’) traits that characterize some putative relict populations on the Queen Charlotte Islands. We found a significant association between lake elevation and the presence of the TNPC; the TNPC was present more often in lakes located at or lower than 42 m than in higher elevation lakes. Main conclusions Our data support the hypothesis that post‐glacial colonization by TNPC‐bearing marine sticklebacks and aspects of lake ‘accessibility’ were important in determining the distribution of mtDNA clades in the eastern Pacific Ocean basin. More generally, our study demonstrates how processes acting both across immense geographic scales (e.g. pan‐Pacific dispersal) and local scales (lake accessibility contingent on timing and extent of isostatic rebound) may interact to explain biogeographical patterns.     

Population genomic evidence for multiple Pliocene refugia in a montane‐restricted harvestman (Arachnida,Opiliones, Sclerobunus robustus) from the southwestern United States

The integration of ecological niche modelling into phylogeographic analyses has allowed for the identification and testing of potential refugia under a hypothesis‐based framework, where the expected patterns of higher genetic diversity in refugial populations and evidence of range expansion of nonrefugial populations are corroborated with empirical data. In this study, we focus on a montane‐restricted cryophilic harvestman, Sclerobunus robustus, distributed throughout the heterogeneous Southern Rocky Mountains and Intermontane Plateau of southwestern North America. We identified hypothetical refugia using ecological niche models (ENMs) across three time periods, corroborated these refugia with population genetic methods using double‐digest RAD‐seq data and conducted population‐level phylogenetic and divergence dating analyses. ENMs identify two large temporally persistent regions in the mid‐latitude highlands. Genetic patterns support these two hypothesized refugia with higher genetic diversity within refugial populations and evidence for range expansion in populations found outside hypothesized refugia. Phylogenetic analyses identify five to six genetically divergent, geographically cohesive clades of S. robustus. Divergence dating analyses suggest that these separate refugia date to the Pliocene and that divergence between clades pre‐dates the late Pleistocene glacial cycles, while diversification within clades was likely driven by these cycles. Population genetic analyses reveal effects of both isolation by distance (IBD) and isolation by environment (IBE), with IBD more important in the continuous mountainous portion of the distribution, while IBE was stronger in the populations inhabiting the isolated sky islands of the south. Using model‐based coalescent approaches, we find support for postdivergence migration between clades from separate refugia.     

Across the southern Andes on fin: glacial refugia,drainage reversals and a secondary contact zone revealed by the phylogeographical signal of Galaxias platei in Patagonia

We employed DNA sequence variation at two mitochondrial (control region, COI) regions from 212 individuals of Galaxias platei (Pisces, Galaxiidae) collected throughout Patagonia (25 lakes/rivers) to examine how Andean orogeny and the climatic cycles throughout the Quaternary affected the genetic diversity and phylogeography of this species. Phylogenetic analyses revealed four deep genealogical lineages which likely represent the initial division of G. platei into eastern and western lineages by Andean uplift, followed by further subdivision of each lineage into separate glacial refugia by repeated Pleistocene glacial cycles. West of the Andes, refugia were likely restricted to the northern region of Patagonia with small relicts in the south, whereas eastern refugia appear to have been much larger and widespread, consisting of separate northern and southern regions that collectively spanned most of Argentinean Patagonia. The retreat of glacial ice following the last glacial maximum allowed re‐colonization of central Chile from nonlocal refugia from the north and east, representing a region of secondary contact between all four glacial lineages. Northwestern glacial relicts likely followed pro‐glacial lakes into central Chilean Patagonia, whereas catastrophic changes in drainage direction (Atlantic → Pacific) for several eastern palaeolakes were the likely avenues for invasions from the east. These mechanisms, combined with evidence for recent, rapid and widespread population growth could explain the extensive contemporary distribution of G. platei throughout Patagonia.     

Chromosome polymorphisms track trans‐Atlantic divergence and secondary contact in Atlantic salmon

Pleistocene glaciations drove repeated range contractions and expansions shaping contemporary intraspecific diversity. Atlantic salmon (Salmo salar) in the western and eastern Atlantic diverged >600,000 years before present, with the two lineages isolated in different southern refugia during glacial maxima, driving trans‐Atlantic genomic and karyotypic divergence. Here, we investigate the genomic consequences of glacial isolation and trans‐Atlantic secondary contact using 108,870 single nucleotide polymorphisms genotyped in 80 North American and European populations. Throughout North America, we identified extensive interindividual variation and discrete linkage blocks within and between chromosomes with known trans‐Atlantic differences in rearrangements: Ssa01/Ssa23 translocation and Ssa08/Ssa29 fusion. Spatial genetic analyses suggest independence of rearrangements, with Ssa01/Ssa23 showing high European introgression (>50%) in northern populations indicative of post‐glacial trans‐Atlantic secondary contact, contrasting with low European ancestry genome‐wide (3%). Ssa08/Ssa29 showed greater intrapopulation diversity, suggesting a derived chromosome fusion polymorphism that evolved within North America. Evidence of potential selection on both genomic regions suggests that the adaptive role of rearrangements warrants further investigation in Atlantic salmon. Our study highlights how Pleistocene glaciations can influence large‐scale intraspecific variation in genomic architecture of northern species.     

Molecular markers provide insights into contemporary and historic gene flow for a non‐migratory species

Hairy woodpeckers Picoides villosus are a common, year round resident with distinct plumage and morphological variation across North America. We genotyped 335 individuals at six variable microsatellite loci and analyzed 322 mtDNA control region sequences in order to examine the role of contemporary and historical barriers to gene flow. In addition we combined genetic analyses with ecological niche modelling to test if hairy woodpeckers were isolated in northern refugia (Alaska, Newfoundland and the Queen Charlotte Islands) during the last glacial maximum. Genetic analyses revealed that gene flow among North American hairy woodpecker populations is restricted, but not to the extent predicted for a sedentary species. Populations clustered into two main genetic groups, east and west of the Great Plains in the south and the Rocky Mountains in the north. Contact zones between the two main genetic groups exist in central British Columbia and Washington, but are narrow. Within each group we found additional population structure with genetic breaks between subgroups in the geographic west corresponding to breaks in forested habitat and physical barriers like open expanses of water. Population genetic patterns for hairy woodpeckers have resulted from isolation in multiple southern refugia with the current distribution of genetic groups resulting from post‐glacial expansion and subsequent reduction in gene flow. While populations in Alaska, Newfoundland and the Queen Charlotte Islands are genetically distinct from other populations, we found no evidence of these areas acting as refugia throughout the Pleistocene. Atlantic Canada populations contained unique haplotypes raising the possibility of a separate colonization from the rest of eastern Canada. The endemic subspecies on the island of Newfoundland is not genetically distinct from their closest mainland population unlike the Queen Charlotte Island subspecies.     

Biogeographic barriers,Pleistocene refugia,and climatic gradients in the southeastern Nearctic drive diversification in cornsnakes (Pantherophis guttatus complex)

The southeastern Nearctic is a biodiversity hotspot that is also rich in cryptic species. Numerous hypotheses (e.g., vicariance, local adaptation, and Pleistocene speciation in glacial refugia) have been tested in an attempt to explain diversification and the observed pattern of extant biodiversity. However, previous phylogeographic studies have both supported and refuted these hypotheses. Therefore, while data support one or more of these diversification hypotheses, it is likely that taxa are forming within this region in species‐specific ways. Here, we generate a genomic data set for the cornsnakes (Pantherophis guttatus complex), which are widespread across this region, spanning both biogeographic barriers and climatic gradients. We use phylogeographic model selection combined with hindcast ecological niche models to determine regions of habitat stability through time. This combined approach suggests that numerous drivers of population differentiation explain the current diversity of this group of snakes. The Mississippi River caused initial speciation in this species complex, with more recent divergence events linked to adaptations to ecological heterogeneity and allopatric Pleistocene refugia. Lastly, we discuss the taxonomy of this group and suggest there may be additional cryptic species in need of formal recognition.