skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Pleistocene persistence and expansion in tarantulas on the Colorado Plateau and the effects of missing data on phylogeographical inferences from RADseq
Abstract 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 whetherAphonopelma 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.  more » « less
Award ID(s):
1754030
PAR ID:
10452369
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
Wiley-Blackwell
Date Published:
Journal Name:
Molecular Ecology
Volume:
29
Issue:
19
ISSN:
0962-1083
Page Range / eLocation ID:
p. 3684-3701
Format(s):
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 Many coastal marine species experienced Pleistocene gene flow between the North Pacific and Atlantic. Understanding historical connectivity between ocean basins should aid in predicting how regional faunas will respond to recent warming that has intensified trans‐Arctic dispersal. Wetland fauna of the Northwestern Atlantic may have survived in estuarine refugia throughout glacial cycles, or recolonised from the southern coast, North Pacific or Northeastern Atlantic. Here, we used multilocus genetic markers and historical climate data to investigate lineage distribution and connectivity among populations of the nominally cosmopolitan sea slugAlderia modesta, sampled from mudflats on both coasts of the North Pacific and North Atlantic. Mitochondrial DNA clades from European and North American populations were deeply divergent and reciprocally monophyletic; differences at seven polymorphic nuclear loci indicated prolonged absence of trans‐Atlantic gene flow. A Pacific ancestor likely first colonised the Atlantic during the marine biotic interchange of the middle Pliocene ~3.5 Ma. Both mtDNA phylogenetics and nuclear genotype assignments support repeated trans‐Arctic colonisation of the Northwestern Atlantic from the Pacific during inter‐glacial cycles; no gene flow was evident since the last glacial maximum, however. Time‐calibrated coalescent phylogenies, Bayesian skyline plots and haplotype networks all indicated recent population expansions in the Pacific and Europe, but not Northwestern Atlantic. In both the Pacific and Northwestern Atlantic, older lineages persisted in patchy refugia north of glacial margins, while a derived clade of Pacific haplotypes indicates northward post‐LGM expansion. The biogeographical history ofAlderiacontrasts with rocky‐shore taxa that were largely extirpated by glacial advance and recolonised from refugia following the last glacial maximum. Based on molecular differences and distinctions in radular and penial stylet morphology, we resurrect the nameAlderia harvardiensisGould 1870 forAlderiafrom the Northwestern Atlantic and North Pacific;A. modestarefers exclusively to European slugs. 
    more » « less
  3. Abstract Quaternary climate change has been strongly linked to distributional shifts and recent species diversification. Montane species, in particular, have experienced enhanced isolation and rapid genetic divergence during glacial fluctuations, and these processes have resulted in a disproportionate number of neo‐endemic species forming in high‐elevation habitats. In temperate montane environments, a general model of alpine population history is well supported, where cold‐specialized species track favourable climate conditions downslope during glacial episodes and upslope during warmer interglacial periods, which leads to a climate‐driven population or species diversification pump. However, it remains unclear how geography mediates distributional changes and whether certain episodes of glacial history have differentially impacted rates of diversification. We address these questions by examining phylogenomic data in a North American clade of flightless, cold‐specialized insects, the ice crawlers (Insecta: Grylloblattodea: Grylloblattidae:Grylloblatta). These low‐vagility organisms have the potential to reveal highly localized refugia and patterns of spatial recolonization, as well as a longer history of in situ diversification. Using continuous phylogeographic analysis of species groups, we show that all species tend to retreat to nearby low‐elevation habitats across western North America during episodes of glaciation, but species at high latitude exhibit larger distributional shifts. Lineage diversification was examined over the course of the Neogene and Quaternary periods, with statistical analysis supporting a direct association between climate variation and diversification rate. Major increases in lineage diversification appear to be correlated with warm and dry periods, rather than with extreme glacial events. Finally, we identify substantial cryptic diversity among ice crawlers, leading to high endemism across their range. This diversity provides new insights into highly localized glacial refugia for cold‐specialized species across western North America. 
    more » « less
  4. Abstract AimArctic plants survived the Pleistocene glaciations in unglaciated refugia. The number, ages, and locations of these refugia are often unclear. We use high‐resolution genomic data from present‐day and Little‐Ice‐Age populations of Arctic Bell‐Heather to re‐evaluate the biogeography of this species and determine whether it had multiple independent refugia or a single refugium in Beringia. LocationCircumpolar Arctic and Coastal British Columbia (BC) alpine. TaxonCassiope tetragonaL., subspeciessaximontanaandtetragona, outgroupC. mertensiana(Ericaceae). MethodsWe built genotyping‐by‐sequencing (GBS) libraries usingCassiope tetragonatissue from 36 Arctic locations, including two ~250‐ to 500‐year‐old populations collected under glacial ice on Ellesmere Island, Canada. We assembled a de novo GBS reference to call variants. Population structure, genetic diversity and demography were inferred from PCA, ADMIXTURE, fastsimcoal2, SplitsTree, and several population genomics statistics. ResultsPopulation structure analyses identified 4–5 clusters that align with geographic locations. Nucleotide diversity was highest in Beringia and decreased eastwards across Canada. Demographic coalescent analyses dated the following splits with Alaska: BC subspeciessaximontana(5 mya), Russia (~1.4 mya), Europe (>200–600 kya), and Greenland (~60 kya). Northern Canada populations appear to have formed during the current interglacial (7–9 kya). Admixture analyses show genetic variants from Alaska appear more frequently in present‐day than historic plants on Ellesmere Island. ConclusionsPopulation and demographic analyses support BC, Alaska, Russia, Europe and Greenland as all having had independent Pleistocene refugia. Northern Canadian populations appear to be founded during the current interglacial with genetic contributions from Alaska, Europe and Greenland. We found evidence, on Ellesmere Island, for continued recent gene flow in the last 250–500 years. These results suggest that a re‐analysis of other Arctic species with shallow population structure using higher resolution genomic markers and demographic analyses may help reveal deeper structure and other circumpolar glacial refugia. 
    more » « less
  5. Although many studies have examined how taxa responded to Pleistocene climate fluctuations in the Appalachian Mountains, impacts on high-elevation endemics of Central Appalachia are not yet understood. We use mitochondrial (ND4 & Cytb) and nuclear (GAPD) DNA sequences to investigate the phylogeography of the Cow Knob Salamander (Plethodon punctatus), a woodland species from Central Appalachian highlands thought to have origins in the Pleistocene. Data from 72 tail tips representing 25 sites revealed that the species comprises two geographically cohesive mitochondrial clades with a narrow, putative contact zone on Shenandoah Mountain. Molecular clock estimates indicate the clades diverged in the Middle Pleistocene. The population size of the Southern clade appears to have remained stable for at least 50,000 years. Despite spanning several isolated mountain systems, the Northern clade has exceptionally low genetic diversity, probably due to recent demographic expansion. Palaeodemographic hypothesis testing supported a scenario in which a founder effect characterized the Northern clade as it diverged from the Southern clade. Species distribution models predicted no suitable habitat for the species during the Last Glacial Maximum. Ultimately, Pleistocene glacial climates may have driven the species from the northern half of its current range, with recolonization events by members of the Northern clade as climates warmed. Density dependent processes may now maintain a narrow contact zone between the two clades. 
    more » « less