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.

Genome-environment associations along elevation gradients in two snowbed species of the North-Eastern Calcareous Alps.

BACKGROUND: Anthropogenic climate change leads to increasing temperatures and altered precipitation and snowmelt patterns, especially in alpine ecosystems. To understand species' responses to climate change, assessment of genetic structure and diversity is crucial as the basis for the evaluation of migration patterns, genetic adaptation potential as well as the identification of adaptive alleles. RESULTS: We studied genetic structure, diversity and genome-environment associations of two snowbed species endemic to the Eastern Alps with a large elevational range, Achillea clusiana Tausch and Campanula pulla L. Genotyping-by-sequencing was employed to assemble loci de novo, call variants and perform population genetic analyses. Populations of either species were distinguishable by mountain, and to some extent by elevation. We found evidence for gene flow between elevations. Results of genome-environment associations suggested similar selective pressures acting on both species, emanating mainly from precipitation and exposition rather than temperature. CONCLUSIONS: Given their genetic structure and amount of gene flow among populations the two study species are suitable to serve as a model for genetic monitoring of climate change adaptation along an elevation gradient. Consequences of climate change will predominantly manifest via changes in precipitation and, thus, duration of snow cover in the snowbeds and indirectly via shrub encroachment accompanied by increasing shading of snowbeds at lower range margins. Assembling genomes of the study species and studying larger sample sizes and time series will be necessary to functionally characterize and validate the herein identified genomic loci putatively involved in adaptive processes.

Diversified chemical profiles of cuticular wax on alpine meadow plants of the Qinghai-Tibet Plateau.

MAIN CONCLUSION: The alpine meadow plants showed great intra- and inter-genera variations of chemical profiles of cuticular waxes. Developing an understanding of wax structure-function relationships that will help us tackle global climate change requires a detailed understanding of plant wax chemistry. The goal in this study was to provide a catalog of wax structures, abundances, and compositions on alpine meadow plants. Here, leaf waxes from 33 plant species belonging to 11 families were sampled from alpine meadows of the east side of the Qinghai-Tibet Plateau. Across these species, total wax coverage varied from 2.30 µg cm-2 to 40.70 µg cm-2, showing variation both within as well as between genera and suggesting that wax variation is subject to both environmental and genetic effects. Across all wax samples, more than 140 wax compounds belonging to 13 wax compound classes were identified, including both ubiquitous wax compounds and lineage-specific compounds. Among the ubiquitous compounds (primary alcohols, alkyl esters, aldehydes, alkanes, and fatty acids), chain length profiles across a wide range of species point to key differences in the chain length specificity of alcohol and alkane formation machinery. The lineage-specific wax compound classes (diols, secondary alcohols, lactones, iso-alkanes, alkyl resorcinols, phenylethyl esters, cinnamate esters, alkyl benzoates, and triterpenoids) nearly all consisted of isomers with varying chain lengths or functional group positions, making the diversity of specialized wax compounds immense. The comparison of species relationships between chemical data and genetic data highlighted the importance of inferring phylogenetic relationships from data sets that contain a large number of variables that do not respond to environmental stimuli.

Extent of intraspecific trait variability in ecologically central and marginal populations of a dominant alpine plant across European mountains.

BACKGROUND AND AIMS: Studying trait variability and restricted gene flow between populations of species can reveal species dynamics. Peripheral populations commonly exhibit lower genetic diversity and trait variability due to isolation and ecological marginality, unlike central populations experiencing gene flow and optimal conditions. This study focused on Carex curvula, the dominant species in alpine acidic meadows of European mountain regions. The species is sparser in dry areas such as the Pyrenees and Balkans, compared to the Central-Eastern Alps and Carpathians. We hypothesized that distinct population groups could be identified based on their mean functional trait values and their correlation with the environment; we predicted that ecologically marginal populations would have stronger trait correlations, lower within-population trait variability (intraspecific trait variability, ITV) and lower genetic diversity than populations of optimal habitats. METHODS: Sampling was conducted in 34 populations that spanned the entire distribution range of C. curvula. We used hierarchical clustering to identify emergent functional groups of populations, defined by combinations of multiple traits associated with nutrient economy and drought tolerance (e.g. specific leaf area, anatomy). We contrasted the geographical distribution of these groups in relation to environment and genetic structure. We compared pairwise trait relationships, within-population trait variation (ITV) and neutral genetic diversity between groups. KEY RESULTS: Our study identified emergent functional groups of populations. Those in the southernmost ranges, specifically the Pyrenees and Balkan region, showed drought-tolerant trait syndromes and correlated with indicators of limited water availability. While we noted a decline in population genetic diversity, we did not observe any significant changes in ITV in ecologically marginal (peripheral) populations. CONCLUSIONS: Our research exemplifies the relationship between ecological marginality and geographical peripherality, which in this case study is linked to genetic depauperation but not to reduced ITV. Understanding these relationships is crucial for understanding the biogeographical factors shaping trait variation.

Development of SSR markers for Carex curvula (Cyperaceae) and their importance in investigating the species genetic structure.

BACKGROUND: Microsatellite primers were developed and tested to genotype several populations of Carex curvula s. l. (Cyperaceae), in order to infer the phylogeographic relationships of the populations within species and the boundaries between the two described subspecies: C. curvula subsp. curvula and C. curvula subsp. rosae. METHODS AND RESULTS: Candidate microsatellite loci were isolated based on next-generation sequencing. We tested 18 markers for polymorphism and replicability in seven C. curvula s. l. populations and identified 13 polymorphic loci with dinucleotide repeats. Genotyping results showed the total number of alleles per locus varied from four to 23 (including both infrataxa), and the observed and expected heterozygosity ranged between 0.1 to 0.82 and 0.219 to 0.711, respectively. Furthermore, the NJ tree showed a clear separation between C. curvula subsp. curvula and C. curvula subsp. rosae. CONCLUSION: The development of these highly polymorphic markers proved to be very efficient not only in delineating between the two subspecies, but also in genetic discriminating at population level within each infrataxon. They are promising tools for evolutionary studies in Cariceae section, as well as in providing knowledge on patterns of the species phylogeography.

Genomic convergence underlying high-altitude adaptation in alpine plants.

Evolutionary convergence is one of the most striking examples of adaptation driven by natural selection. However, genomic evidence for convergent adaptation to extreme environments remains scarce. Here, we assembled reference genomes of two alpine plants, Saussurea obvallata (Asteraceae) and Rheum alexandrae (Polygonaceae), with 37,938 and 61,463 annotated protein-coding genes. By integrating an additional five alpine genomes, we elucidated genomic convergence underlying high-altitude adaptation in alpine plants. Our results detected convergent contractions of disease-resistance genes in alpine genomes, which might be an energy-saving strategy for surviving in hostile environments with only a few pathogens present. We identified signatures of positive selection on a set of genes involved in reproduction and respiration (e.g., MMD1, NBS1, and HPR), and revealed signatures of molecular convergence on genes involved in self-incompatibility, cell wall modification, DNA repair and stress resistance, which may underlie adaptation to extreme cold, high ultraviolet radiation and hypoxia environments. Incorporating transcriptomic data, we further demonstrated that genes associated with cuticular wax and flavonoid biosynthetic pathways exhibit higher expression levels in leafy bracts, shedding light on the genetic mechanisms of the adaptive "greenhouse" morphology. Our integrative data provide novel insights into convergent evolution at a high-taxonomic level, aiding in a deep understanding of genetic adaptation to complex environments.

Ecological lags govern the pace and outcome of plant community responses to 21st-century climate change.

Forecasting the trajectories of species assemblages in response to ongoing climate change requires quantifying the time lags in the demographic and ecological processes through which climate impacts species' abundances. Since experimental climate manipulations are typically abrupt, the observed species responses may not match their responses to gradual climate change. We addressed this problem by transplanting alpine grassland turfs to lower elevations, recording species' demographic responses to climate and competition, and using these data to parameterise community dynamics models forced by scenarios of gradual climate change. We found that shifts in community structure following an abrupt climate manipulation were not simply accelerated versions of shifts expected under gradual warming, as the former missed the transient rise of species benefiting from moderate warming. Time lags in demography and species interactions controlled the pace and trajectory of changing species' abundances under simulated 21st-century climate change, and thereby prevented immediate diversity loss.

Acclimation to water stress improves tolerance to heat and freezing in a common alpine grass.

Alpine plants in Australia are increasingly exposed to more frequent drought and heatwaves, with significant consequences for physiological stress responses. Acclimation is a critical feature that allows plants to improve tolerance to environmental extremes by directly altering their physiology or morphology. Yet it is unclear how plant performance, tolerance, and recovery are affected when heat and water stress co-occur, and whether prior exposure affects responses to subsequent climate extremes. We grew a common alpine grass species under high or low watering treatments for three weeks before exposure to either none, one, or two heat stress events. We determined photosynthetic heat and freezing tolerance (LT50, mean temperature causing 50% irreversible damage to photosystem II) and growth. Physiological adjustments to low watering, including more negative water potentials and reduced growth, were also characterised by improved tolerance to high and low-temperature extremes. Shifts to higher heat tolerance were also evident with increasing exposure to heat stress events, though freezing tolerance was not affected. Acclimation effects were mostly short-term, however; prior exposure to heat and/or water stress had little to no effect on growth and thermal tolerance following the six-week recovery period. We conclude that rapid acclimation to water and heat stress that co-occur during summer enhances the capacity of alpine plants to tolerate increasingly frequent temperature extremes.

Decades-long phylogeographic issues: complex historical processes and ecological factors on genetic structure of alpine plants in the Japanese Archipelago.

Mountain regions are important places for biodiversity, where organisms could persist throughout prolonged periods and accumulate genetic divergence as well as promote speciation. Roles of mountains for biodiversity have been exclusively discussed in regions that have specifically diverse species or covered with ice-sheets during the Pleistocene glacial periods, whereas the importance of mountainous regions in East Asia has been less disputed. High mountains in the Japanese Archipelago, located at the eastern edge of the Eurasia continent, have one of southernmost populations of alpine and arctic-alpine plants that are also distributed in the northern Pacific and/or the circumarctic regions. Phylogeographic studies on the Japanese alpine plants have excluded their possible ephemeral occurrence during the current warm period, and rather, suggest persistence of alpine plants throughout several cycles of climate changes in the Pleistocene on high mountains in central Honshu, the main island of the Japanese Archipelago. In this review, I look through decade long phylogeographic studies and show complicated patterns of range dynamics of Japanese alpine plants. In addition, I note recent findings of genetic relationships of Japanese populations of alpine and/or arctic-alpine plants with those in northern regions and their possible ecological divergence in the Japanese Archipelago. Taken together, I provide several issues for understanding historical processes that established distribution of alpine plants following climate changes as well as their diversification and propose importance of Japanese populations of alpine plants on biodiversity in alpine communities across broader range, especially in the northern Pacific region.

Divergence in red light responses associated with thermal reversion of phytochrome B between high- and low-latitude species.

Phytochromes play a central role in mediating adaptive responses to light and temperature throughout plant life cycles. Despite evidence for adaptive importance of natural variation in phytochromes, little information is known about molecular mechanisms that modulate physiological responses of phytochromes in nature. We show evolutionary divergence in physiological responses relevant to thermal stability of a physiologically active form of phytochrome (Pfr) between two sister species of Brassicaceae growing at different latitudes. The higher latitude species (Cardamine bellidifolia; Cb) responded more strongly to light-limited conditions compared with its lower latitude sister (C. nipponica; Cn). Moreover, CbPHYB conferred stronger responses to both light-limited and warm conditions in the phyB-deficient mutant of Arabidopsis thaliana than CnPHYB: that is Pfr CbphyB was more stable in nuclei than CnphyB. Our findings suggest that fine tuning Pfr stability is a fundamental mechanism for plants to optimise phytochrome-related traits in their evolution and adapt to spatially varying environments, and open a new avenue to understand molecular mechanisms that fine tune phytochrome responses in nature.

Climate-driven range shifts reduce persistence of competitors in a perennial plant community.

Forecasting the impacts of climate change on species persistence in diverse natural communities requires a way to account for indirect effects mediated through species interactions. In particular, we expect species to experience major changes in competition as they track favorable climates. Here, we combine experimental data with a recently developed theoretical framework based on coexistence theory to measure the impact of climate-driven range shifts on alpine plant persistence under climate change. We transplanted three co-dominant alpine perennial species to five elevations, creating a maximum of 5°C increase in average growing-season temperature. We statistically modeled species' demographic rates in response to the environment and interpolated species' intrinsic ranges-the environmental mapping of reproduction in the absence of competition. We used low-density population growth rates-species' initial rate of invasion into an established community-as a metric of persistence. Further analysis of low-density growth rates (LGRs) allowed us to parse the direct impacts of climate change from indirect impacts mediated by shifting competition. Our models predict qualitatively different range shifts for each species based on the climate conditions under which growth rates are maximized and where they are zero. Overall, climate change is predicted to increase the intrinsic (competition free) growth rates of all species, as warmer and wetter conditions increase the favorability of alpine habitat. However, these benefits are entirely negated by increased competition arising from greater overlap between competitors in their intrinsic ranges. Species were highly dispersal limited, which can prevent species from tracking shifting intrinsic ranges by reducing population spread rates. Yet dispersal limitation also promoted species' persistence because it promotes persistence mechanisms. Our study demonstrates the complex pathways by which climate change impacts species' persistence by altering their competitive environment, and highlights how a persistence framework based on LGRs can help disentangle impacts.

Genetic structure of a widespread alpine shrub Rhododendron aureum (Ericaceae) across East Asia.

The vast territory of East Asia, including southwestern Beringia, is considered to have been almost ice free during the Pleistocene. Cold-resistant flora may have persisted in this region expanding or contracting its range during the climate cooling. Only a few plant genera have been studied with a sampling area across their entire geographic range in East Asia; therefore, the understanding of the biogeographic history of alpine flora in this region remains limited. In the present study, genetic variation and population structure in 21 populations of the alpine shrub Rhododendron aureum across its range in East Asia were assessed using 18 microsatellite loci. Phylogenetic analyses revealed three main genetic groups: Siberia, Northeast, and North Pacific. According to the geographical pattern of genetic diversity, the North Pacific group includes populations from Kamchatka, south of Russian Far East, and territories close to central Japan. This group is the most diverse and likely diverged earlier than the Siberia and Northeast groups. Ecological niche modeling predicts range expansion of this species during the period of cooling and, together with demographic history, suggests that the divergence between the three main genetic groups predated the Last Glacial Maximum. Similar to other cold-resistant species such as Larix sibirica and Juniperus communis, the pattern of genetic diversity of R. aureum supports the survival of the species at high latitudes during the Pleistocene with limited contribution of the southern populations to expansion of the species range to the Northeast region and Siberia.

Environmental and geographic variables are effective surrogates for genetic variation in conservation planning.

Protected areas buffer species from anthropogenic threats and provide places for the processes that generate and maintain biodiversity to continue. However, genetic variation, the raw material for evolution, is difficult to capture in conservation planning, not least because genetic data require considerable resources to obtain and analyze. Here we show that freely available environmental and geographic distance variables can be highly effective surrogates in conservation planning for representing adaptive and neutral intraspecific genetic variation. We obtained occurrence and genetic data from the IntraBioDiv project for 27 plant species collected over the European Alps using a gridded sampling scheme. For each species, we identified loci that were potentially under selection using outlier loci methods, and mapped their main gradients of adaptive and neutral genetic variation across the grid cells. We then used the cells as planning units to prioritize protected area acquisitions. First, we verified that the spatial patterns of environmental and geographic variation were correlated, respectively, with adaptive and neutral genetic variation. Second, we showed that these surrogates can predict the proportion of genetic variation secured in randomly generated solutions. Finally, we discovered that solutions based only on surrogate information secured substantial amounts of adaptive and neutral genetic variation. Our work paves the way for widespread integration of surrogates for genetic variation into conservation planning.

Winter Nights during Summer Time: Stress Physiological Response to Ice and the Facilitation of Freezing Cytorrhysis by Elastic Cell Wall Components in the Leaves of a Nival Species.

Ranunculus glacialis grows and reproduces successfully, although the snow-free time period is short (2-3 months) and night frosts are frequent. At a nival site (3185 m a.s.l.), we disentangled the interplay between the atmospheric temperature, leaf temperatures, and leaf freezing frequency to assess the actual strain. For a comprehensive understanding, the freezing behavior from the whole plant to the leaf and cellular level and its physiological after-effects as well as cell wall chemistry were studied. The atmospheric temperatures did not mirror the leaf temperatures, which could be 9.3 °C lower. Leaf freezing occurred even when the air temperature was above 0 °C. Ice nucleation at on average -2.6 °C started usually independently in each leaf, as the shoot is deep-seated in unfrozen soil. All the mesophyll cells were subjected to freezing cytorrhysis. Huge ice masses formed in the intercellular spaces of the spongy parenchyma. After thawing, photosynthesis was unaffected regardless of whether ice had formed. The cell walls were pectin-rich and triglycerides occurred, particularly in the spongy parenchyma. At high elevations, atmospheric temperatures fail to predict plant freezing. Shoot burial prevents ice spreading, specific tissue architecture enables ice management, and the flexibility of cell walls allows recurrent freezing cytorrhysis. The peculiar patterning of triglycerides close to ice rewards further investigation.

Identification of the Volatile Components of Galium verum L. and Cruciata leavipes Opiz from the Western Italian Alps.

The chemical composition of the volatile fraction from Galium verum L. (leaves and flowers) and Cruciata laevipes Opiz (whole plant), Rubiaceae, was investigated. Samples from these two plant species were collected at full bloom in Val di Susa (Western Alps, Turin, Italy), distilled in a Clevenger-type apparatus, and analyzed by GC/FID and GC/MS. A total of more than 70 compounds were identified, making up 92%-98% of the total oil. Chemical investigation of their essential oils indicated a quite different composition between G. verum and C. laevipes, both in terms of the major constituents and the dominant chemical classes of the specialized metabolites. The most abundant compounds identified in the essential oils from G. verum were 2-methylbenzaldheyde (26.27%, corresponding to 11.59 µg/g of fresh plant material) in the leaves and germacrene D (27.70%; 61.63 µg/g) in the flowers. C. laevipes essential oils were instead characterized by two sesquiterpenes, namely ß-caryophyllene (19.90%; 15.68 µg/g) and trans-muurola-4(15),5-diene (7.60%; 5.99 µg/g); two phenylpropanoids, benzyl alcohol (8.30%; 6.71 µg/g), and phenylacetaldehyde (7.74%; 6.26 µg/g); and the green-leaf alcohol cis-3-hexen-1-ol (9.69%; 7.84 µg/g). The ecological significance of the presence of such compounds is discussed.

Microclimate and demography interact to shape stable population dynamics across the range of an alpine plant.

Heterogeneous terrain in montane systems results in a decoupling of climatic gradients. Population dynamics across species' ranges in these heterogeneous landscapes are shaped by relationships between demographic rates and these interwoven climate gradients. Linking demography and climate variables across species' ranges refines our understanding of the underlying mechanisms of species' current and future ranges. We explored the importance of multiple microclimatic gradients in shaping individual demographic rates and population growth rates in 16 populations across the elevational distribution of an alpine plant (Ivesia lycopodioides var. scandularis). Using integral projection modeling, we ask how each rate varies across three microclimate gradients: accumulated degree-days, growing-season soil moisture, and days of snow cover. Range-wide variation in demographic rates was best explained by the combined influence of multiple microclimatic variables. Different pairs of demographic rates exhibited both similar and inverse responses to the same microclimatic gradient, and the microclimatic effects often varied with plant size. These responses resulted in range-wide projected population persistence, with no declining populations at either elevational range edge or at the extremes of the microclimate gradients. The complex relationships between topography, microclimate and demography suggest that populations across a species' range may have unique demographic pathways to stable population dynamics.

Environmental effects on fine-scale spatial genetic structure in four Alpine keystone forest tree species.

Genetic responses to environmental changes take place at different spatial scales. While the effect of environment on the distribution of species' genetic diversity at large geographical scales has been the focus of several recent studies, its potential effects on genetic structure at local scales are understudied. Environmental effects on fine-scale spatial genetic structure (FSGS) were investigated in four Alpine conifer species (five to eight populations per species) from the eastern Italian Alps. Significant FSGS was found for 11 of 25 populations. Interestingly, we found no significant differences in FSGS across species but great variation among populations within species, highlighting the importance of local environmental factors. Interannual variability in spring temperature had a small but significant effect on FSGS of Larix decidua, probably related to species-specific life history traits. For Abies alba, Picea abies and Pinus cembra, linear models identified spring precipitation as a potentially relevant climate factor associated with differences in FSGS across populations; however, models had low explanatory power and were strongly influenced by a P. cembra outlier population from a very dry site. Overall, the direction of the identified effects is according to expectations, with drier and more variable environments increasing FSGS. Underlying mechanisms may include climate-related changes in the variance of reproductive success and/or environmental selection of specific families. This study provides new insights on potential changes in local genetic structure of four Alpine conifers in the face of environmental changes, suggesting that new climates, through altering FSGS, may also have relevant impacts on plant microevolution.

Unchanged risk of frost exposure for subalpine and alpine plants after snowmelt in Switzerland despite climate warming.

The length of the snow-free season is a key factor regulating plant phenology and shaping plant community composition in cold regions. While global warming has significantly advanced the time of snowmelt and the growth period at all elevations in the Swiss Alps, it remains unclear if it has altered the likelihood of frost risk for alpine plants. Here, we analyzed the influence of the snowmelt timing on the risk of frost exposure for subalpine and alpine plants shortly after snowmelt, i.e., during their most vulnerable period to frost at the beginning of their growth period. Furthermore, we tested whether recent climate warming has changed the risk of exposure of plants to frost after snowmelt. We analyzed snow and air temperature data in the Swiss Alps using six weather stations covering the period 1970-2016 and 77 weather stations covering the period 1998-2016, spanning elevations from 1418 to 2950 m asl. When analyzed across all years within each station, our results showed strong negative relationships between the time of snowmelt and the frequency and intensity of frost during the most vulnerable period to frost for subalpine and alpine plants, indicating a higher frost risk damage for plants during years with earlier snowmelt. However, over the last 46 years, the time of snowmelt and the last spring frost date have advanced at similar rates, so that the frequency and intensity of frost during the vulnerable period for plants remained unchanged.

Arctic and boreal plant species decline at their southern range limits in the Rocky Mountains.

Climate change is predicted to cause a decline in warm-margin plant populations, but this hypothesis has rarely been tested. Understanding which species and habitats are most likely to be affected is critical for adaptive management and conservation. We monitored the density of 46 populations representing 28 species of arctic-alpine or boreal plants at the southern margin of their ranges in the Rocky Mountains of Montana, USA, between 1988 and 2014 and analysed population trends and relationships to phylogeny and habitat. Marginal populations declined overall during the past two decades; however, the mean trend for 18 dicot populations was -5.8% per year, but only -0.4% per year for the 28 populations of monocots and pteridophytes. Declines in the size of peripheral populations did not differ significantly among tundra, fen and forest habitats. Results of our study support predicted effects of climate change and suggest that vulnerability may depend on phylogeny or associated anatomical/physiological attributes.