skip to main content

Title: Genotyping‐by‐sequencing and ecological niche modeling illuminate phylogeography, admixture, and Pleistocene range dynamics in quaking aspen ( Populus tremuloides )

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
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
Ecology and Evolution
Page Range / eLocation ID:
p. 4609-4629
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    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 ofS. 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.

    more » « less
  2. Abstract

    Understanding the processes that shape genetic diversity by either promoting or preventing population divergence can help identify geographic areas that either facilitate or limit gene flow. Furthermore, broadly distributed species allow us to understand how biogeographic and ecogeographic transitions affect gene flow. We investigated these processes using genomic data in the Northern Alligator Lizard (Elgaria coerulea), which is widely distributed in Western North America across diverse ecoregions (California Floristic Province and Pacific Northwest) and mountain ranges (Sierra Nevada, Coastal Ranges, and Cascades). We collected single-nucleotide polymorphism data from 120 samples of E. coerulea. Biogeographic analyses of squamate reptiles with similar distributions have identified several shared diversification patterns that provide testable predictions for E. coerulea, including deep genetic divisions in the Sierra Nevada, demographic stability of southern populations, and recent post-Pleistocene expansion into the Pacific Northwest. We use genomic data to test these predictions by estimating the structure, connectivity, and phylogenetic history of populations. At least 10 distinct populations are supported, with mixed-ancestry individuals situated at most population boundaries. A species tree analysis provides strong support for the early divergence of populations in the Sierra Nevada Mountains and recent diversification into the Pacific Northwest. Admixture and migration analyses detect gene flow among populations in the Lower Cascades and Northern California, and a spatial analysis of gene flow identified significant barriers to gene flow across both the Sierra Nevada and Coast Ranges. The distribution of genetic diversity in E. coerulea is uneven, patchy, and interconnected at population boundaries. The biogeographic patterns seen in E. coerulea are consistent with predictions from co-distributed species.

    more » « less
  3. Abstract

    Phylogeographic studies can uncover robust details about the population structure, demographics, and diversity of species. The smooth greensnake, Opheodrys vernalis, is a small, cryptic snake occupying mesic grassland and sparsely wooded habitats. Although O. vernalis has a wide geographical range, many metapopulations are patchy and some are declining. We used mitochondrial DNA and double digest restriction-site associated DNA sequencing to construct the first phylogeographic assessment of O. vernalis. Genomic analysis of 119 individuals (mitochondrial DNA) and a subset of another 45 smooth greensnakes (nuclear DNA; N = 3031 single nucleotide polymorphisms) strongly supports two longitudinally separated lineages, with admixture in the Great Lakes region. Post-Pleistocene secondary contact best explains admixture from populations advancing northwards. Overall, populations expressed low heterozygosity, variable inbreeding rates, and moderate to high differentiation. Disjunct populations in the Rocky Mountains and central Great Plains regions might be contracting relicts, whereas northerly populations in more continuous mesic habitats (e.g., Prairie Pothole region, southern Canada) had signals of population expansion. Broadly, conservation management efforts should be focused on local populations, because habitat connectivity may facilitate gene flow and genetic diversity.

    more » « less
  4. Abstract Aim

    We used genome‐scale sampling to assess the phylogeography of a group of topminnows in theFundulus notatusspecies complex. Two of the species have undergone extensive range expansions resulting in broadly overlapping distributions, and sympatry within drainages has provided opportunities for hybridization and introgression. We assessed the timing and pattern of range expansion in the context of late Pleistocene–Holocene drainage events and evaluated the evidence for introgressive hybridization between species.


    Central and southern United States including drainages of the Gulf of Mexico Coastal Plain and portions of the Mississippi River drainage in and around the Central Highlands.


    Topminnows, GenusFundulus, subgenusZygonectesFundulus notatus, Fundulus olivaceus, Fundulus euryzonus.


    We sampled members of theF. notatusspecies complex throughout their respective ranges, including numerous drainage systems where species co‐occur. We collected genome‐wide single nucleotide polymorphisms (SNPs) using the genotype‐by‐sequencing (GBS) method and subjected data to population genetic analyses to infer the population histories of both species, including explicit tests for admixture and introgression. The methods employed includedSTRUCTURE, principal coordinates analysis, TreeMix and approximate Bayesian computation.


    Genetic data are presented for 749 individuals sampled from 14F. notatus, 20F. olivaceusand 2F. euryzonuspopulations. Members of the species complex differed in phylogeographic structure, withF. notatusexhibiting geographic clusters corresponding to Pleistocene coastal drainages andF. olivaceuscomparatively lacking in phylogeographic structure. Evidence for interspecific introgression varied by drainage.

    Main conclusions

    Populations ofF. notatusandF. olivaceusexhibited contrasting patterns of lineage diversity among coastal drainages, indicating interspecific differences in their Pleistocene southern refugia. Phylogeographic patterns in both species indicated that range expansions into the northern limits of contemporary distributions coincided with and continued subsequent to the Last Glacial Maximum. There was evidence of introgression between species in some, but not all drainages where the species co‐occur, in a pattern that is correlated with previous estimates of hybridization rates.

    more » « less
  5. Abstract

    Across eastern North America, glacial cycles of the Pleistocene drove episodic latitudinal range shifts by temperate species. Isolation of populations within low-latitude refugia during glacial maxima was enhanced by physiographic barriers, leading to patterns of phylogeographic differentiation that are shared across diverse taxa. Postglacial population expansion created opportunities for differentiated lineages to come into contact, with various potential population-genetic outcomes. Northern short-tailed shrews (Blarina brevicauda) exhibit three mitochondrial phylogroups that probably originated via glacial-age range restriction and isolation. We investigate the history of postglacial expansion and interlineage contact between historically isolated regional populations of B. brevicauda. Morphological differences between skulls of shrews representing a Western lineage and those representing Central and Eastern lineages are consistent with past subspecies delineations. However, we demonstrate broad range overlap between neighboring phylogroups across the Upper Peninsula and Lower Peninsula in Michigan. Further, incongruence between phylogroup association and morphology among individuals in Upper Peninsula populations suggests that genetic admixture between shrews representing the Western and Central groups has occurred in the past and may be ongoing. We show that across most cranial measurements, shrews within the contact zone are morphologically most similar to the Central group regardless of mitochondrial identity, but one measurement in these contact zone shrews (depth of skull) is more similar to that seen in the Western group. These results suggest that hybridization between historically isolated populations has resulted in the origin of a novel skull phenotype that is proportionally deeper, narrower, and shorter than those seen in core Western and Central populations.

    more » « less