skip to main content


Title: The role of multiple Pleistocene refugia in promoting diversification in the Pacific Northwest
Abstract

Pleistocene glacial cycles drastically changed the distributions of taxa endemic to temperate rainforests in the Pacific Northwest, with many experiencing reduced habitat suitability during glacial periods. In this study, we investigate whether glacial cycles promoted intraspecific divergence and whether subsequent range changes led to secondary contact and gene flow. For seven invertebrate species endemic to the PNW, we estimated species distribution models (SDMs) and projected them onto current and historical climate conditions to assess how habitat suitability changed during glacial cycles. Using single nucleotide polymorphism (SNP) data from these species, we assessed population genetic structure and used a machine‐learning approach to compare models with and without gene flow between populations upon secondary contact after the last glacial maximum (LGM). Finally, we estimated divergence times and rates of gene flow between populations. SDMs suggest that there was less suitable habitat in the North Cascades and Northern Rocky Mountains during glacial compared to interglacial periods, resulting in reduced habitat suitability and increased habitat fragmentation during the LGM. Our genomic data identify population structure in all taxa, and support gene flow upon secondary contact in five of the seven taxa. Parameter estimates suggest that population divergences date to the later Pleistocene for most populations. Our results support a role of refugial dynamics in driving intraspecific divergence in the Cascades Range. In these invertebrates, population structure often does not correspond to current biogeographic or environmental barriers. Rather, population structure may reflect refugial lineages that have since expanded their ranges, often leading to secondary contact between once isolated lineages.

 
more » « less
Award ID(s):
2009989
NSF-PAR ID:
10369530
Author(s) / Creator(s):
 ;  ;  ;  ;  
Publisher / Repository:
Wiley-Blackwell
Date Published:
Journal Name:
Molecular Ecology
Volume:
31
Issue:
16
ISSN:
0962-1083
Page Range / eLocation ID:
p. 4402-4416
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Aim

    Refugial isolation during glaciation is an established driver of speciation; however, the opposing role of interglacial population expansion, secondary contact, and gene flow on the diversification process remains less understood. The consequences of glacial cycling on diversity are complex and especially so for archipelago species, which experience dramatic fluctuations in connectivity in response to both lower sea levels during glacial events and increased fragmentation during glacial recession. We test whether extended refugial isolation has led to the divergence of genetically and morphologically distinct species within Holarctic ermine (Mustela erminea), a small cosmopolitan carnivore species that harbours 34 extant subspecies, 14 of which are insular endemics.

    Location

    Holarctic.

    Methods

    We use genetic sequences (complete mitochondrial genomes, four nuclear genes) from >100 ermine (stoats) and geometric morphometric data for >200 individuals (27 of the 34 extant subspecies) from across their Holarctic range to provide an integrative perspective on diversification and endemism across this complex landscape. Multiple species delimitation methods (iBPP,bPTP) assessed congruence between morphometric and genetic data.

    Results

    Our results support the recognition of at least three species within theM. ermineacomplex, coincident with three of four genetic clades, tied to diversification in separate glacial refugia. We found substantial geographic variation within each species, with geometric morphometric results largely consistent with historical infraspecific taxonomy.

    Main conclusions

    Phylogeographic structure mirrors patterns of diversification in other Holarctic species, with a major Nearctic‐Palearctic split, but with greater intraspecific morphological diversity. Recognition of insular endemic speciesM. haidarumis consistent with a deep history of refugial persistence and highlights the urgency of mindful management of island populations along North America's North Pacific Coast. Significant environmental modification (e.g. industrial‐scale logging, mining) has been proposed for a number of these islands, which may elevate the risk of extinction of insular palaeoendemics.

     
    more » « less
  2. Abstract

    Secondary sympatry amongst sister lineages is strongly associated with genetic and ecological divergence. This pattern suggests that for closely related species to coexist in secondary sympatry, they must accumulate differences in traits that mediate ecological and/or reproductive isolation. Here, we characterized inter‐ and intraspecific divergence in three giant tree frog species whose distributions stretch across West and Central Africa. Using genome‐wide single‐nucleotide polymorphism data, we demonstrated that species‐level divergence coincides temporally and geographically with a period of large‐scale forest fragmentation during the late Pliocene. Our environmental niche models further supported a dynamic history of climatic suitability and stability, and indicated that all three species occupy distinct environmental niches. We found modest morphological differentiation amongst the species with significant divergence in tympanum diameter and male advertisement call. In addition, we confirmed that two species occur in secondary sympatry in Central Africa but found no evidence of hybridization. These patterns support the hypothesis that cycles of genetic exchange and isolation across West and Central Africa have contributed to globally significant biodiversity. Furthermore, divergence in both ecology and reproductive traits appear to have played important roles in maintaining distinct lineages. At the intraspecific level, we found that climatic refugia, precipitation gradients, marine incursions, and potentially riverine barriers generated phylogeographic structure throughout the Pleistocene and into the Holocene. Further studies examining phenotypic divergence and secondary contact amongst these geographically structured populations may demonstrate how smaller scale and more recent biogeographic barriers contribute to regional diversification.

     
    more » « less
  3. Abstract

    Determining how genetic diversity is structured between populations that span the divergence continuum from populations to biological species is key to understanding the generation and maintenance of biodiversity. We investigated genetic divergence and gene flow in eight lineages of birds with a trans‐Beringian distribution, where Asian and North American populations have likely been split and reunited through multiple Pleistocene glacial cycles. Our study transects the speciation process, including eight pairwise comparisons in three orders (ducks, shorebirds and passerines) at population, subspecies and species levels. Using ultraconserved elements (UCEs), we found that these lineages represent conditions from slightly differentiated populations to full biological species. Although allopatric speciation is considered the predominant mode of divergence in birds, all of our best divergence models included gene flow, supporting speciation with gene flow as the predominant mode in Beringia. In our eight lineages, three were best described by a split‐migration model (divergence with gene flow), three best fit a secondary contact scenario (isolation followed by gene flow), and two showed support for both models. The lineages were not evenly distributed across a divergence space defined by gene flow (M) and differentiation (FST), instead forming two discontinuous groups: one with relatively shallow divergence, no fixed single nucleotide polymorphisms (SNPs), and high rates of gene flow between populations; and the second with relatively deeply divergent lineages, multiple fixed SNPs, and low gene flow. Our results highlight the important role that gene flow plays in avian divergence in Beringia.

     
    more » « less
  4. Abstract

    Populus tremuloidesis the widest‐ranging tree species in North America and an ecologically important component of mesic forest ecosystems displaced by the Pleistocene glaciations. Using phylogeographic analyses of genome‐wide SNPs (34,796 SNPs, 183 individuals) and ecological niche modeling, we inferred population structure, ploidy levels, admixture, and Pleistocene range dynamics ofP. tremuloides, and tested several historical biogeographical hypotheses. We found three genetic lineages located mainly in coastal–Cascades (cluster 1), east‐slope Cascades–Sierra Nevadas–Northern Rockies (cluster 2), and U.S. Rocky Mountains through southern Canadian (cluster 3) regions of theP. tremuloidesrange, with tree graph relationships of the form ((cluster 1, cluster 2), cluster 3). Populations consisted mainly of diploids (86%) but also small numbers of triploids (12%) and tetraploids (1%), and ploidy did not adversely affect our genetic inferences. The main vector of admixture was from cluster 3 into cluster 2, with the admixture zone trending northwest through the Rocky Mountains along a recognized phenotypic cline (Utah to Idaho). Clusters 1 and 2 provided strong support for the “stable‐edge hypothesis” that unglaciated southwestern populations persisted in situ since the last glaciation. By contrast, despite a lack of clinal genetic variation, cluster 3 exhibited “trailing‐edge” dynamics from niche suitability predictions signifying complete northward postglacial expansion. Results were also consistent with the “inland dispersal hypothesis” predicting postglacial assembly of Pacific Northwestern forest ecosystems, but rejected the hypothesis that Pacific‐coastal populations were colonized during outburst flooding from glacial Lake Missoula. Overall, congruent patterns between our phylogeographic and ecological niche modeling results and fossil pollen data demonstrate complex mixtures of stable‐edge, refugial locations, and postglacial expansion withinP. tremuloides. These findings confirm and refine previous genetic studies, while strongly supporting a distinct Pacific‐coastal genetic lineage of quaking aspen.

     
    more » « less
  5. Abstract

    Divergence is often ephemeral, and populations that diverge in response to regional topographic and climatic factors may not remain reproductively isolated when they come into secondary contact. We investigated the geographical structure and evolutionary history of population divergence withinSceloporus occidentalis(western fence lizard), a habitat generalist with a broad distribution that spans the major biogeographical regions of Western North America. We used double digest RAD sequencing to infer population structure, phylogeny and demography. Population genetic structure is hierarchical and geographically structured with evidence for gene flow between biogeographical regions. Consistent with the isolation–expansion model of divergence during Quaternary glacial–interglacial cycles, gene flow and secondary contact are supported as important processes explaining the demographic histories of populations. Although populations may have diverged as they spread northward in a ring‐like manner around the Sierra Nevada and southern Cascade Ranges, there is strong evidence for gene flow among populations at the northern terminus of the ring. We propose the concept of an “ephemeral ring species” and contrastS. occidentaliswith the classic North American ring species,Ensatina eschscholtzii. Contrary to expectations of lower genetic diversity at northern latitudes following post‐Quaternary‐glaciation expansion, the ephemeral nature of divergence inS. occidentalishas produced centres of high genetic diversity for different reasons in the south (long‐term stability) vs. the north (secondary contact).

     
    more » « less