How Do Cold-Adapted Plants Respond to Climatic Cycles? Interglacial Expansion Explains Current Distribution and Genomic Diversity in Primula farinosa L

Unveiling the Genome-Wide Consequences of Range Expansion and Mating System Transitions in Primula vulgaris

Genetic diversity is heterogeneously distributed among populations of the same species, due to the joint effects of multiple demographic processes, including range contractions and expansions, and mating systems shifts. Here, we ask how both processes shape genomic diversity in space and time in the classical Primula vulgaris model. This perennial herb originated in the Caucasus region and was hypothesized to have expanded westward following glacial retreat in the Quaternary. Moreover, this species is a long-standing model for mating system transitions, exemplified by shifts from heterostyly to homostyly. Leveraging a high-quality reference genome of the closely related Primula veris and whole-genome resequencing data from both heterostylous and homostylous individuals from populations encompassing a wide distribution of P. vulgaris, we reconstructed the demographic history of P. vulgaris. Results are compatible with the previously proposed hypothesis of range expansion from the Caucasus region approximately 79,000 years ago and suggest later shifts to homostyly following rather than preceding postglacial colonization of England. Furthermore, in accordance with population genetic theoretical predictions, both processes are associated with reduced genetic diversity, increased linkage disequilibrium, and reduced efficacy of purifying selection. A novel result concerns the contrasting effects of range expansion versus shift to homostyly on transposable elements, for the former, process is associated with changes in transposable element genomic content, while the latter is not. Jointly, our results elucidate how the interactions among range expansion, transitions to selfing, and Quaternary climatic oscillations shape plant evolution.

Alpine Extremophytes in Evolutionary Turmoil: Complex Diversification Patterns and Demographic Responses of a Halophilic Grass in a Central Asian Biodiversity Hotspot

Diversification and demographic responses are key processes shaping species evolutionary history. Yet we still lack a full understanding of ecological mechanisms that shape genetic diversity at different spatial scales upon rapid environmental changes. In this study, we examined genetic differentiation in an extremophilic grass Puccinellia pamirica and factors affecting its population dynamics among the occupied hypersaline alpine wetlands on the arid Pamir Plateau in Central Asia. Using genomic data, we found evidence of fine-scale population structure and gene flow among the localities established across the high-elevation plateau as well as fingerprints of historical demographic expansion. We showed that an increase in the effective population size could coincide with the Last Glacial Period, which was followed by the species demographic decline during the Holocene. Geographic distance plays a vital role in shaping the spatial genetic structure of P. pamirica alongside with isolation-by-environment and habitat fragmentation. Our results highlight a complex history of divergence and gene flow in this species-poor alpine region during the Late Quaternary. We demonstrate that regional climate specificity and a shortage of nonclimate data largely impede predictions of future range changes of the alpine extremophile using ecological niche modeling. This study emphasizes the importance of fine-scale environmental heterogeneity for population dynamics and species distribution shifts.

Survival in nunatak and peripheral glacial refugia of three alpine plant species is partly predicted by altitudinal segregation

Mountain biota survived the Quaternary cold stages most probably in peripheral refugia and/or ice-free peaks within ice-sheets (nunataks). While survival in peripheral refugia has been broadly demonstrated, evidence for nunatak refugia is still scarce. We generated RADseq data from three mountain plant species occurring at different elevations in the southeastern European Alps to investigate the role of different glacial refugia during the Last Glacial Maximum (LGM). We tested the following hypotheses. (i) The deep Piave Valley forms the deepest genetic split in the species distributed across it, delimiting two peripheral refugia. (ii) The montane to alpine species Campanula morettiana and Primula tyrolensis survived the LGM in peripheral refugia, while high-alpine to subnival Saxifraga facchinii likely survived in several nunatak refugia. (iii) The lower elevation species suffered a strong population decline during the LGM. By contrast, the higher elevation species shows long-term stability of population sizes due to survival on permanently ice-free peaks and small population sizes at present. We found peripheral refugia on both sides of the Piave Valley, which acted as a major genetic barrier. Demographic modelling confirmed nunatak survival not only for S. facchinii but also for montane to alpine C. morettiana. Altitudinal segregation influenced the species' demographic fluctuations, with the lower elevation species showing a significant population increase at the end of the LGM, and the higher elevation species either showing decrease towards the present or stable population sizes with a short bottleneck. Our results highlight the role of nunatak survival and species ecology in the demographic history of mountain species.

The Role of the Hercynian Mountains of Central Europe in Shaping Plant Migration Patterns in the Pleistocene-A Review

The climatic changes that took place in Europe during the Quaternary period influenced plant habitats as well as their species and vegetation composition. In this article, biogeographical studies on Hercynian mountain plants that include data for the Alps, Carpathians, and European lowlands are reviewed in order to discuss the phylogeographical structure and divergence of the Hercynian populations from those in other European mountain ranges, Scandinavia, and lowlands. The analyzed studies show specific phylogeographical relations between the Hercynian mountains, Alps, Scandinavia, Carpathians, and European lowlands. The results also indicate that the genetic patterns of plant populations in the Hercynian Mountains may differ significantly in terms of origin. The main migration routes of species to the Hercynian ranges began in the Alps or Carpathians. Some species, such as Rubus chamaemorus L., Salix lapponum L., and Salix herbacea L., are glacial relics that may have arrived and settled in the Hercynian Mountains during the Ice Age and that survived in isolated habitats. The Hercynian Mountains are composed of various smaller mountain ranges and are a crossroads of migration routes from different parts of Europe; thus, intensive hybridization has occurred between the plant populations therein, which is indicated by the presence of several divergent genetic lines.

Climate has contributed to population diversification of Daphnia galeata across Eurasia

Climate is a fundamental abiotic factor that plays a key role in driving the evolution, distribution and population diversification of species. However, there have been few investigations of genomic signatures of adaptation to local climatic conditions in cladocerans. Here, we have provided the first high-quality chromosome-level genome assembly (~143 Mb, scaffold N50 12.6 Mb) of the waterflea, Daphnia galeata, and investigated genomic variation in 22 populations from Central Europe and Eastern China. Our ecological-niche models suggested that the historic distribution of D. galeata in Eurasia was significantly affected by Quaternary climate fluctuations. We detected pronounced genomic and morphometric divergences between European and Chinese D. galeata populations. Such divergences could be partly explained by genomic signatures of thermal adaptation to distinct climate regimes: a set of candidate single-nucleotide polymorphisms (SNPs) potentially associated with climate were detected. These SNPs were in genes significantly enriched in the Gene ontology terms "determination of adult lifespan" and "translation repressor activity", and especially, mthl5 and SOD1 involved in the IIS pathway, and EIF4EBP2 involved in the target of the rapamycin signalling pathway. Our study indicates that certain alleles might be associated with particular temperature regimes, playing a functional role in shaping the population structure of D. galeata at a large geographical scale. These results highlight the potential role of molecular variation in the response to climate variation, in the context of global climate change.

Latitudinal Diversity Gradient in the Changing World: Retrospectives and Perspectives

The latitudinal diversity gradient (LDG) is one of the most extensive and important biodiversity patterns on the Earth. Various studies have established that species diversity increases with higher taxa numbers from the polar to the tropics. Studies of multicellular biotas have supported the LDG patterns from land (e.g., plants, animals, forests, wetlands, grasslands, fungi, and so forth) to oceans (e.g., marine organisms from freshwater invertebrates, continental shelve, open ocean, even to the deep sea invertebrates). So far, there are several hypotheses proposed to explore the diversity patterns and mechanisms of LDG, however, there has been no consensus on the underlying causes of LDG over the past few decades. Thus, we reviewed the progress of LDG studies in recent years. Although several explanations for the LDG have been proposed, these hypotheses are only based on species richness, evolution and the ecosystems. In this review, we summarize the effects of evolution and ecology on the LDG patterns to synthesize the formation mechanisms of the general biodiversity distribution patterns. These intertwined factors from ecology and evolution in the LDG are generally due to the wider distribution of tropical areas, which hinders efforts to distinguish their relative contributions. However, the mechanisms of LDG always engaged controversies, especially in such a context that the human activity and climate change has affected the biodiversity. With the development of molecular biology, more genetic/genomic data are available to facilitate the estimation of global biodiversity patterns with regard to climate, latitude, and other factors. Given that human activity and climate change have inevitably impacted on biodiversity loss, biodiversity conservation should focus on the change in LDG pattern. Using large-scale genetic/genomic data to disentangle the diversity mechanisms and patterns of LDG, will provide insights into biodiversity conservation and management measures. Future perspectives of LDG with integrative genetic/genomic, species, evolution, and ecosystem diversity patterns, as well as the mechanisms that apply to biodiversity conservation, are discussed. It is imperative to explore integrated approaches for recognizing the causes of LDG in the context of rapid loss of diversity in a changing world.

Congruent evolutionary responses of European steppe biota to late Quaternary climate change

Quaternary climatic oscillations had a large impact on European biogeography. Alternation of cold and warm stages caused recurrent glaciations, massive vegetation shifts, and large-scale range alterations in many species. The Eurasian steppe biome and its grasslands are a noteworthy example; they underwent climate-driven, large-scale contractions during warm stages and expansions during cold stages. Here, we evaluate the impact of these range alterations on the late Quaternary demography of several phylogenetically distant plant and insect species, typical of the Eurasian steppes. We compare three explicit demographic hypotheses by applying an approach combining convolutional neural networks with approximate Bayesian computation. We identified congruent demographic responses of cold stage expansion and warm stage contraction across all species, but also species-specific effects. The demographic history of the Eurasian steppe biota reflects major paleoecological turning points in the late Quaternary and emphasizes the role of climate as a driving force underlying patterns of genetic variance on the biome level.

Glacial Expansion or Interglacial Expansion? Contrasting Demographic Models of Four Cold-Adapted Fir Species in North America and East Asia

Understanding and forecasting species’ response to climate change is a critical need for future conservation and management. Two expansion hypotheses, the glacial expansion versus the interglacial expansion, have been proposed to interpret how cold-adapted organisms in the northern hemisphere respond to Quaternary climatic fluctuations. To test these two hypotheses, we originally used two pairs of high-low elevation firs from North America (Abies lasiocarpa and Abies balsamea) and East Asia (Abies chensiensis and Abies nephrolepis). Abies lasiocarpa and Abies chensiensis are widely distributed in high-elevation regions of western North America and central China. Abies balsamea and Abies nephrolepis occur in central North America and northeast China, with much lower elevations. These fir species are typical cold-adapted species and sensitive to climate fluctuations. Here, we integrated the mtDNA and cpDNA polymorphisms involving 44 populations and 585 individuals. Based on phylogeographic analyses, recent historical range expansions were indicated in two high-elevation firs (Abies lasiocarpa and Abies chensiensis) during the last glaciation (43.8–28.4 or 21.9–14.2 kya, 53.1–34.5 or 26.6–17.2 kya). Such glacial expansions in high-elevation firs were further confirmed by the evidence of species distribution modelling, geographic-driven genetic patterns, palynological records, and current distribution patterns. Unlike the north American firs, the SDM models indicated unremarkable expansion or contraction in East Asia firs for its much more stable conditions during different historical periods. Taken together, our findings highly supported that high-elevation firs experienced glacial expansion during the Quaternary climate change in East Asia and North America, as interglacial expansion within low-elevation firs. Under this situation, the critically endangered fir species distributed in high elevation would have no enough higher elevational space to migrate. Facing the increasing global warming, thus we proposed ex-situ conservation of defining conservation units as the most meaningful strategy.

Climate-Change Impacts on the Southernmost Mediterranean Arctic-Alpine Plant Populations

Human-induced climate- and land-use change have been affecting biogeographical and biodiversity patterns for the past two centuries all over the globe, resulting in increased extinction and biotic homogenization rates. High mountain ecosystems are more sensitive to these changes, which have led to physiological and phenological shifts, as well as to ecosystem processes’ deformation. Glacial relicts, such as arctic-alpine taxa, are sensitive indicators of the effects of global warming and their rear-edge populations could include warm-adapted genotypes that might prove—conservation-wise—useful in an era of unprecedented climate regimes. Despite the ongoing thermophilization in European and Mediterranean summits, it still remains unknown how past and future climate-change might affect the distributional patterns of the glacial relict, arctic-alpine taxa occurring in Greece, their European southernmost distributional limit. Using species distribution models, we investigated the impacts of past and future climate changes on the arctic-alpine taxa occurring in Greece and identified the areas comprising arctic-alpine biodiversity hotspots in Greece. Most of these species will be faced with severe range reductions in the near future, despite their innate resilience to a multitude of threats, while the species richness hotspots will experience both altitudinal and latitudinal shifts. Being long-lived perennials means that there might be an extinction-debt present in these taxa, and a prolonged stability phase could be masking the deleterious effects of climate change on them. Several ex situ conservation measures (e.g., seed collection, population augmentation) should be taken to preserve the southernmost populations of these rare arctic-alpine taxa and a better understanding of their population genetics is urgently needed.

Addressing alpine plant phylogeography using integrative distributional, demographic and coalescent modeling

Phylogeographic studies of alpine plants have evolved considerably in the last two decades from ad hoc interpretations of genetic data to statistical model-based approaches. In this review we outline the developments in alpine plant phylogeography focusing on the recent approach of integrative distributional, demographic and coalescent (iDDC) modeling. By integrating distributional data with spatially explicit demographic modeling and subsequent coalescent simulations, the history of alpine species can be inferred and long-standing hypotheses, such as species-specific responses to climate change or survival on nunataks during the last glacial maximum, can be efficiently tested as exemplified by available case studies. We also discuss future prospects and improvements of iDDC.

Population Genetic Structure and Demographic History of Primula fasciculata in Southwest China

Understanding the factors that drive the genetic structure of a species and its responses to past climatic changes is an important first step in modern population management. The response to the last glacial maximum (LGM) has been well studied, however, the effect of previous glaciation periods on plant demographic history is still not well studied. Here we investigated the population structure and demographic history of Primula fasciculata that widely occurs in the Hengduan Mountains and Qinghai-Tibetan Plateau. We obtained genomic data for 234 samples of the species using restriction site-associated DNA (RAD) sequencing and combined approximate Bayesian computation (ABC) and species distribution modeling (SDM) to evaluate the effects of multiple glaciation periods by testing several population divergence models and demographic scenarios. The analyses of population structure showed that P. fasciculata displays a striking population structure with six groups that could be identified genetically. Our ABC modeling suggested that the current groups diverged from ancestral populations located in the eastern Hengduan Mountains after the largest glaciation occurred in the region (~ 0.8-0.5 million years ago), which is consistent with the result of SDMs. Each current group has survived in different glacial refugia during the LGM and experienced expansions and/or bottlenecks since their divergence during or across the following Quaternary glacial cycles. Our study demonstrates the usefulness of population genomics for evaluating the effects of past climatic changes in alpine plant species with shallow population structure.

Evolutionary history and past climate change shape the distribution of genetic diversity in terrestrial mammals

Knowledge of global patterns of biodiversity, ranging from intraspecific genetic diversity (GD) to taxonomic and phylogenetic diversity, is essential for identifying and conserving the processes that shape the distribution of life. Yet, global patterns of GD and its drivers remain elusive. Here we assess existing biodiversity theories to explain and predict the global distribution of GD in terrestrial mammal assemblages. We find a strong positive covariation between GD and interspecific diversity, with evolutionary time, reflected in phylogenetic diversity, being the best predictor of GD. Moreover, we reveal the negative effect of past rapid climate change and the positive effect of inter-annual precipitation variability in shaping GD. Our models, explaining almost half of the variation in GD globally, uncover the importance of deep evolutionary history and past climate stability in accumulating and maintaining intraspecific diversity, and constitute a crucial step towards reducing the Wallacean shortfall for an important dimension of biodiversity.

The drivers of genetic diversity (GD) are poorly understood at the global scale. Here the authors show, for terrestrial mammals, that within-species GD covaries with phylogenetic diversity and is higher in locations with more stable past climates. They also interpolate GD for data-poor locations such as the tropics.

Ecotypic differentiation, hybridization and clonality facilitate the persistence of a cold-adapted sedge in European bogs

Recent research has shown that many cold-adapted species survived the last glacial maximum (LGM) in northern refugia. Whether this evolutionary history has had consequences for their genetic diversity and adaptive potential remains unknown. We sampled 14 populations of Carex limosa, a sedge specialized to bog ecosystems, along a latitudinal gradient from its Scandinavian core to the southern lowland range-margin in Germany. Using microsatellite and experimental common-garden data, we evaluated the impacts of global climate change along this gradient and assessed the conservation status of the southern marginal populations. Microsatellite data revealed two highly distinct genetic groups and hybrid individuals. In our common-garden experiment, the two groups showed divergent responses to increased nitrogen/phosphorus (N/P) availability, suggesting ecotypic differentiation. Each group formed genetically uniform populations at both northern and southern sampling areas. Mixed populations occurred throughout our sampling area, an area that was entirely glaciated during the LGM. The fragmented distribution implies allopatric divergence at geographically separated refugia that putatively differed in N/P availability. Molecular data and an observed low hybrid fecundity indicate the importance of clonal reproduction for hybrid populations. At the southern range-margin, however, all populations showed effects of clonality, lowered fecundity and low competitiveness, suggesting abiotic and biotic constraints to population persistence.

Is the incidence of survival in interior Pleistocene refugia (nunataks) underestimated? Phylogeography of the high mountain plant Androsace alpina (Primulaceae) in the European Alps revisited

Temperate mountain ranges such as the European Alps have been strongly affected by the Pleistocene glaciations. Glacial advances forced biota into refugia, which were situated either at the periphery of mountain ranges or in their interior. Whereas in the Alps peripheral refugia have been repeatedly and congruently identified, support for the latter scenario, termed "nunatak hypothesis," is still limited and no general pattern is recognizable yet. Here, we test the hypothesis of nunatak survival for species growing in the high alpine to subnival zones on siliceous substrate using the cushion plant Androsace alpina (Primulaceae), endemic to the European Alps, as our model species. To this end, we analyzed AFLP and plastid DNA sequence data obtained from a dense and range-wide sampling. Both AFLPs and plastid sequence data identified the southwestern-most population as the most divergent one. AFLP data did not allow for discrimination of interior and peripheral populations, but rather identified two to three longitudinally separated major gene pools. In contrast, in the eastern half of the Alps several plastid haplotypes of regional or local distribution in interior ranges-the Alpine periphery mostly harbored a widespread haplotype-were indicative for the presence of interior refugia. Together with evidence from other Alpine plant species, this study shows that in the eastern Alps silicicolous species of open habitats in the alpine and subnival zone survived, also or exclusively so, in interior refugia. As the corresponding genetic structure may be lost in mostly nuclear-derived, rapidly homogenizing marker systems such as AFLPs or RAD sequencing tags, markers not prone to homogenization, as is the case for plastid sequences (Sanger-sequenced or extracted from an NGS data set) will continue to be important for detecting older, biogeographically relevant patterns.