skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 5:00 PM ET until 11:00 PM ET on Friday, June 21 due to maintenance. We apologize for the inconvenience.

Title: Area Not Geographic Isolation Mediates Biodiversity Responses of Alpine Refugia to Climate Change
Climate refugia, where local populations of species can persist through periods of unfavorable regional climate, play a key role in the maintenance of regional biodiversity during times of environmental change. However, the ability of refugia to buffer biodiversity change may be mediated by the landscape context of refugial habitats. Here, we examined how plant communities restricted to refugial sky islands of alpine tundra in the Colorado Rockies are changing in response to rapid climate change in the region (increased temperature, declining snowpack, and earlier snow melt-out) and if these biodiversity changes are mediated by the area or geographic isolation of the sky island. We resampled plant communities in 153 plots at seven sky islands distributed across the Colorado Rockies at two time points separated by 12 years (2007/2008–2019/2020) and found changes in taxonomic, phylogenetic, and functional diversity over time. Specifically, we found an increase in species richness, a trend toward increased phylogenetic diversity, a shift toward leaf traits associated with the stress-tolerant end of leaf economics spectrum (e.g., lower specific leaf area, higher leaf dry matter content), and a decrease in the functional dispersion of specific leaf area. Importantly, these changes were partially mediated by refugial area but not by geographic isolation, suggesting that dispersal from nearby areas of tundra does not play a strong role in mediating these changes, while site characteristics associated with a larger area (e.g., environmental heterogeneity, larger community size) may be relatively more important. Taken together, these results suggest that considering the landscape context (area and geographic isolation) of refugia may be critical for prioritizing the conservation of specific refugial sites that provide the most conservation value.  more » « less
Award ID(s):
Author(s) / Creator(s):
Date Published:
Journal Name:
Frontiers in Ecology and Evolution
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 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.




    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.


    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
  3. Abstract

    The paradigm of past climate-driven range shifts structuring the distribution of marine intraspecific biodiversity lacks replication in biological models exposed to comparable limiting conditions in independent regions. This may lead to confounding effects unlinked to climate drivers. We aim to fill in this gap by asking whether the global distribution of intraspecific biodiversity of giant kelp (Macrocystis pyrifera) is explained by past climate changes occurring across the two hemispheres. We compared the species’ population genetic diversity and structure inferred with microsatellite markers, with range shifts and long-term refugial regions predicted with species distribution modelling (SDM) from the last glacial maximum (LGM) to the present. The broad antitropical distribution ofMacrocystis pyriferais composed by six significantly differentiated genetic groups, for which current genetic diversity levels match the expectations of past climate changes. Range shifts from the LGM to the present structured low latitude refugial regions where genetic relics with higher and unique diversity were found (particularly in the Channel Islands of California and in Peru), while post-glacial expansions following ~ 40% range contraction explained extensive regions with homogenous reduced diversity. The estimated effect of past climate-driven range shifts was comparable between hemispheres, largely demonstrating that the distribution of intraspecific marine biodiversity can be structured by comparable evolutionary forces across the global ocean. Additionally, the differentiation and endemicity of regional genetic groups, confers high conservation value to these localized intraspecific biodiversity hotspots of giant kelp forests.

    more » « less
  4. Abstract Background

    Antarctica and its unique biodiversity are increasingly at risk from the effects of global climate change and other human influences. A significant recent element underpinning strategies for Antarctic conservation has been the development of a system of Antarctic Conservation Biogeographic Regions (ACBRs). The datasets supporting this classification are, however, dominated by eukaryotic taxa, with contributions from the bacterial domain restricted to Actinomycetota and Cyanobacteriota. Nevertheless, the ice-free areas of the Antarctic continent and the sub-Antarctic islands are dominated in terms of diversity by bacteria. Our study aims to generate a comprehensive phylogenetic dataset of Antarctic bacteria with wide geographical coverage on the continent and sub-Antarctic islands, to investigate whether bacterial diversity and distribution is reflected in the current ACBRs.


    Soil bacterial diversity and community composition did not fully conform with the ACBR classification. Although 19% of the variability was explained by this classification, the largest differences in bacterial community composition were between the broader continental and maritime Antarctic regions, where a degree of structural overlapping within continental and maritime bacterial communities was apparent, not fully reflecting the division into separate ACBRs. Strong divergence in soil bacterial community composition was also apparent between the Antarctic/sub-Antarctic islands and the Antarctic mainland. Bacterial communities were partially shaped by bioclimatic conditions, with 28% of dominant genera showing habitat preferences connected to at least one of the bioclimatic variables included in our analyses. These genera were also reported as indicator taxa for the ACBRs.


    Overall, our data indicate that the current ACBR subdivision of the Antarctic continent does not fully reflect bacterial distribution and diversity in Antarctica. We observed considerable overlap in the structure of soil bacterial communities within the maritime Antarctic region and within the continental Antarctic region. Our results also suggest that bacterial communities might be impacted by regional climatic and other environmental changes. The dataset developed in this study provides a comprehensive baseline that will provide a valuable tool for biodiversity conservation efforts on the continent. Further studies are clearly required, and we emphasize the need for more extensive campaigns to systematically sample and characterize Antarctic and sub-Antarctic soil microbial communities.

    more » « less
  5. Abstract

    A central goal at the interface of ecology and conservation is understanding how the relationship between biodiversity and ecosystem function (B–EF) will shift with changing climate. Despite recent theoretical advances, studies which examine temporal variation in the functional traits and mechanisms (mass ratio effects and niche complementarity effects) that underpin the B–EF relationship are lacking.

    Here, we use 13 years of data on plant species composition, plant traits, local‐scale abiotic variables, above‐ground net primary productivity (ANPP), and climate from the alpine tundra of Colorado (USA) to investigate temporal dynamics in the B–EF relationship. To assess how changing climatic conditions may alter the B–EF relationship, we built structural equation models (SEMs) for 11 traits across 13 years and evaluated the power of different trait SEMs to predict ANPP, as well as the relative contributions of mass ratio effects (community‐weighted mean trait values; CWM), niche complementarity effects (functional dispersion; FDis) and local abiotic variables. Additionally, we coupled linear mixed effects models with Multimodel inference methods to assess how inclusion of trait–climate interactions might improve our ability to predict ANPP through time.

    In every year, at least one SEM exhibited good fit, explaining between 19.6% and 57.2% of the variation in ANPP. However, the identity of the trait which best explained ANPP changed depending on winter precipitation, with leaf area, plant height and foliar nitrogen isotope content (δ15N) SEMs performing best in high, middle and low precipitation years, respectively. Regardless of trait identity, CWMs exerted a stronger influence on ANPP than FDis and total biotic effects were always greater than total abiotic effects. Multimodel inference reinforced the results of SEM analysis, with the inclusion of climate–trait interactions marginally improving our ability to predict ANPP through time.

    Synthesis. Our results suggest that temporal variation in climatic conditions influences which traits, mechanisms and abiotic variables were most responsible for driving the B–EF relationship. Importantly, our findings suggest that future research should consider temporal variability in the B–EF relationship, particularly how the predictive power of individual functional traits and abiotic variables may fluctuate as conditions shift due to climate change.

    more » « less