Effect of snowmelt timing on the genetic structure of an Erythronium grandiflorum population in an alpine environment

Genetic diversity and differentiation of populations of Anthyllis vulneraria along elevational and latitudinal gradients

The abundant centre model (ACM) predicts that the suitability of environmental conditions for a species decreases from the centre of its distribution toward its range periphery and, consequently, its populations will become scarcer, smaller and more isolated, resulting in lower genetic diversity and increased differentiation. However, little is known about whether genetic diversity shows similar patterns along elevational and latitudinal gradients with similar changes in important environmental conditions. Using microsatellite markers, we studied the genetic diversity and structure of 20 populations each of Anthyllis vulneraria along elevational gradients in the Alps from the valleys to the elevational limit (2500 m) and along a latitudinal gradient (2500 km) from Central Europe to the range margin in northern Scandinavia. Both types of gradients corresponded to an 11.5°C difference in mean annual temperature. Genetic diversity strongly declined and differentiation increased with latitude in line with the predictions of the ACM. However, as population size did not decline with latitude and genetic diversity was not related to population size in A. vulneraria, this pattern is not likely to be due to less favorable conditions in the North, but due to serial founder effects during the post-glacial recolonization process. Genetic diversity was not related to elevation, but we found significant isolation by distance along both gradients, although the elevational gradient was shorter by orders of magnitude. Subarctic populations differed genetically from alpine populations indicating that the northern populations did not originate from high elevational Alpine ones. Our results support the notion that postglacial latitudinal colonization over large distances resulted in a larger loss of genetic diversity than elevational range shifts. The lack of genetic diversity in subarctic populations may threaten their long-term persistence in the face of climate change, whereas alpine populations could benefit from gene flow from low-elevation populations.

Does Phenological Plasticity Help or Hinder Range Shifts Under Climate Change?

Climate warming is predicted to shift species’ ranges as previously uninhabitable environments just beyond the leading range edges become suitable habitat and trailing range edges become increasingly unsuitable. Understanding which aspects of the environment and species traits mediate these range shifts is critical for understanding species’ possible redistributions under global change, yet we have a limited understanding of the ecological and evolutionary responses underlying population spread or extinction at species’ range edges. Within plant populations, shifts in flowering phenology have been one of the strongest and most consistent responses to climate change, and are likely to play an important role in mediating population dynamics within and beyond species’ ranges. However, the role of phenological shifts, and particularly phenological plasticity, in species’ range shifts remains relatively unstudied. Here, we synthesize literature on phenology, plasticity, and adaptation to suggest ways in which phenological responses to climate may vary across species’ ranges and review the empirical evidence for and against these hypotheses. We then outline how phenological plasticity could facilitate or hinder persistence and potential consequences of phenological plasticity in range expansions, including phenological cues, shifts in correlated traits, altered species interactions, and effects on gene flow. Finally, we suggest future avenues for research, such as characterizing reaction norms for phenology across a species’ range and in beyond-the-range transplant experiments. Given the prevalence and magnitude of phenological shifts, future work should carefully dissect its costs and benefits for population persistence, and incorporate phenological plasticity into models predicting species’ persistence and geographic range shifts under climate change.

Predicting Thermal Adaptation by Looking Into Populations' Genomic Past

Molecular evolution offers an insightful theory to interpret the genomic consequences of thermal adaptation to previous events of climate change beyond range shifts. However, disentangling often mixed footprints of selective and demographic processes from those due to lineage sorting, recombination rate variation, and genomic constrains is not trivial. Therefore, here we condense current and historical population genomic tools to study thermal adaptation and outline key developments (genomic prediction, machine learning) that might assist their utilization for improving forecasts of populations' responses to thermal variation. We start by summarizing how recent thermal-driven selective and demographic responses can be inferred by coalescent methods and in turn how quantitative genetic theory offers suitable multi-trait predictions over a few generations via the breeder's equation. We later assume that enough generations have passed as to display genomic signatures of divergent selection to thermal variation and describe how these footprints can be reconstructed using genome-wide association and selection scans or, alternatively, may be used for forward prediction over multiple generations under an infinitesimal genomic prediction model. Finally, we move deeper in time to comprehend the genomic consequences of thermal shifts at an evolutionary time scale by relying on phylogeographic approaches that allow for reticulate evolution and ecological parapatric speciation, and end by envisioning the potential of modern machine learning techniques to better inform long-term predictions. We conclude that foreseeing future thermal adaptive responses requires bridging the multiple spatial scales of historical and predictive environmental change research under modern cohesive approaches such as genomic prediction and machine learning frameworks.

Progressive deterioration of pollination service detected in a 17-year study vanishes in a 26-year study

Widespread reports of declining populations of pollinators have raised concerns that plant populations may be incurring increasing shortfalls in pollination, but few studies have measured pollination deficits over enough seasons to detect such changes. I have conducted pollen-supplementation experiments in a wild population of the glacier lily (Erythronium grandiflorum, Liliaceae) from 1993 to 2018. Pollination deficits were estimated by comparing the fruit set of hand-pollinated, single-flowered plants to that of open-pollinated controls. For a subset of years, seed set data were also available. A previous publication reported a significant deterioration of pollination in this population from 1993 to 2009, and suggested phenological dislocation as a possible cause. That deterioration is no longer evident in the longer-term data set. Very long time series may be necessary to detect temporal trends in pollination service. This population consistently experiences stronger pollination deficits before its flowering peak than after. This heterogeneity suggests caution in characterizing a population as pollination-limited or not, even within a single season.

Interannual variation in season length is linked to strong co-gradient plasticity of phenology in a montane annual plant

Species are commonly distributed along latitudinal and elevational gradients of growing season length to which they might respond via phenotypic plasticity and/or adaptive genetic differentiation. However, the relative contribution of these processes and whether plasticity, if it occurs, facilitates expansion along season-length gradients remain unclear, but are important for predicting species fates during anthropogenic change. We quantified phenological trait variation in the montane annual Rhinanthus minor for three generations at 12 sites across 900 m of elevation in the Canadian Rocky Mountains and conducted a reciprocal transplant experiment for two generations among nine sites. We compared clines and interannual variation of phenological traits between natural and transplanted individuals. Season length declined by c. 37% along our elevational gradient and, as expected, plants emerged, reached first flower and made their first seed in c. 41% fewer growing degree days under shorter growing seasons. Although reciprocal transplants revealed modest genetic differentiation across elevation, trait clines primarily were due to striking co-gradient plasticity that paralleled genetic differentiation. Co-gradient plasticity likely evolved in response to considerable interannual variation in season length across our elevational transect, and should prepare R. minor to make adaptive changes to phenology in response to ongoing climate change predicted for montane environments.

Ecotypic divergences of the alpine herb Potentilla matsumurae adapted to fellfield-snowbed habitats across a series of mountain sky islands

PREMISE

Divergent selection due to environmental heterogeneity can lead to local adaptation. However, the ecological and evolutionary processes of local adaptation that occurs across multiple regions are often unknown. Our previous studies reported on the ecotypic divergence within a local area of variation of Potentilla matsumurae, an alpine herb adapted to the fellfield-snowbed environment. Here we investigated large-scale geographic patterns of ecotypic differentiation in this species to infer local adaptation and selective forces across multiple regions.

METHODS

We compiled information on the overall distributions of fellfield and snowbed habitats on the mountains in Japan across the distribution of the species. Next, we conducted common garden experiments to test the adaptive divergence of the fellfield-snowbed plants derived from multiple regions. Finally, we evaluated phylogeographic structures based on cpDNA and allozyme variations and inferred the evolutionary history of ecotype differentiation.

RESULTS

The mosaic distribution of the fellfield-snowbed ecotypes across isolated mountaintops constitutes indirect evidence for habitat-specific natural selection. The significant difference in survivorship between the ecotypes observed in a controlled snow environment provides more substantial evidence of local selection. Phylogeographic structures support the hypothesis that ecotypic divergence events from fellfield to snowbed populations occurred independently in at least two distinct regions.

CONCLUSIONS

Ecotypic divergence of P. matsumurae has occurred across a series of mountain sky islands. Local selection in snowy environments is a driving force that maintains the divergent ecotypes across multiple mountain regions and can contribute to the diversification of plants in heavy-snow regions.

Evolutionary lessons from California plant phylogeography

Phylogeography documents the spatial distribution of genetic lineages that result from demographic processes, such as population expansion, population contraction, and gene movement, shaped by climate fluctuations and the physical landscape. Because most phylogeographic studies have used neutral markers, the role of selection may have been undervalued. In this paper, we contend that plants provide a useful evolutionary lesson about the impact of selection on spatial patterns of neutral genetic variation, when the environment affects which individuals can colonize new sites, and on adaptive genetic variation, when environmental heterogeneity creates divergence at specific loci underlying local adaptation. Specifically, we discuss five characteristics found in plants that intensify the impact of selection: sessile growth form, high reproductive output, leptokurtic dispersal, isolation by environment, and the potential to evolve longevity. Collectively, these traits exacerbate the impact of environment on movement between populations and local selection pressures-both of which influence phylogeographic structure. We illustrate how these unique traits shape these processes with case studies of the California endemic oak, Quercus lobata, and the western North American lichen, Ramalina menziesii Obviously, the lessons we learn from plant traits are not unique to plants, but they highlight the need for future animal, plant, and microbe studies to incorporate its impact. Modern tools that generate genome-wide sequence data are now allowing us to decipher how evolutionary processes affect the spatial distribution of different kinds of genes and also to better model future spatial distribution of species in response to climate change.

Population structure of Erythronium dens-canis L. (Liliaceae) in the northern Apennines (Italy)

Relationships between age, time of emergence, and leaf traits of individuals were investigated in a population of Erythroniumdens-canis L. in a hilly woodland area named Farneto-C, near Bologna, Italy. In 2015, 591 individuals were counted, 19 of which were flowering (FLO), 442 were mature non-flowering (MNF) and 130 were juveniles (JUV). FLO emerged at the end of February, whereas most MNF and JUV appeared at the middle and end of March, respectively. The mean aboveground survivorship of MNF was 24 days. Most MNF had large, oval to shield-shaped leaves with red-brown mottling, whereas most JUV leaves were smaller, usually oblong or lanceolate with a rough maculation or none. These results suggest that both timing of emergence and leaf shape are related to the age of the bulb. Based on leaf background, plants were classified into three major types with a likely genetic basis in the 2015 and 2016 surveys (the latter limited to FLO): a dominant silvery type (SLV, 62–74%), silvery-and-green type (S&G, 23–32%), and a less frequent vivid-green type (GRN, 3–5%). Several subtypes were also identified, but only one was dominant within each type. The three basic patterns appear to be phenotypically stable and no differences between MNF and FLO were found; once the juvenile stage has passed, each plant produces the same leaf type year after year. In addition, our results on the discoloration time-course of red-brown spots suggest that the functional role of leaf mottling is not related to pollinator attraction. Instead, leaf mottling could play a role in camouflage against herbivores.The observed massive grazing on flowers, more than leaves, could explain why the frequency of mature individuals was biased towards the non-flowering ones.

Limited Pollen Dispersal Contributes to Population Genetic Structure but Not Local Adaptation in Quercus oleoides Forests of Costa Rica

Background

Quercus oleoides Cham. and Schlect., tropical live oak, is a species of conservation importance in its southern range limit of northwestern Costa Rica. It occurs in high-density stands across a fragmented landscape spanning a contrasting elevation and precipitation gradient. We examined genetic diversity and spatial genetic structure in this geographically isolated and genetically distinct population. We characterized population genetic diversity at 11 nuclear microsatellite loci in 260 individuals from 13 sites. We monitored flowering time at 10 sites, and characterized the local environment in order to compare observed spatial genetic structure to hypotheses of isolation-by-distance and isolation-by-environment. Finally, we quantified pollen dispersal distances and tested for local adaptation through a reciprocal transplant experiment in order to experimentally address these hypotheses.

Results

High genetic diversity is maintained in the population and the genetic variation is significantly structured among sampled sites. We identified 5 distinct genetic clusters and average pollen dispersal predominately occurred over short distances. Differences among sites in flowering phenology and environmental factors, however, were not strictly associated with genetic differentiation. Growth and survival of upland and lowland progeny in their native and foreign environments was expected to exhibit evidence of local adaptation due to the more extreme dry season in the lowlands. Seedlings planted in the lowland garden experienced much higher mortality than seedlings in the upland garden, but we did not identify evidence for local adaptation.

Conclusion

Overall, this study indicates that the Costa Rican Q. oleoides population has a rich population genetic history. Despite environmental heterogeneity and habitat fragmentation, isolation-by-distance and isolation-by-environment alone do not explain spatial genetic structure. These results add to studies of genetic structure by examining a common, tropical tree over multiple habitats and provide information for managers of a successional forest in a protected area.

Impact of intra- versus inter-annual snow depth variation on water relations and photosynthesis for two Great Basin Desert shrubs

Snowfall provides the majority of soil water in certain ecosystems of North America. We tested the hypothesis that snow depth variation affects soil water content, which in turn drives water potential (Ψ) and photosynthesis, over 10 years for two widespread shrubs of the western USA. Stem Ψ (Ψ stem) and photosynthetic gas exchange [stomatal conductance to water vapor (g s), and CO2 assimilation (A)] were measured in mid-June each year from 2004 to 2013 for Artemisia tridentata var. vaseyana (Asteraceae) and Purshia tridentata (Rosaceae). Snow fences were used to create increased or decreased snow depth plots. Snow depth on +snow plots was about twice that of ambient plots in most years, and 20 % lower on -snow plots, consistent with several down-scaled climate model projections. Maximal soil water content at 40- and 100-cm depths was correlated with February snow depth. For both species, multivariate ANOVA (MANOVA) showed that Ψ stem, g s, and A were significantly affected by intra-annual variation in snow depth. Within years, MANOVA showed that only A was significantly affected by spatial snow depth treatments for A. tridentata, and Ψ stem was significantly affected by snow depth for P. tridentata. Results show that stem water relations and photosynthetic gas exchange for these two cold desert shrub species in mid-June were more affected by inter-annual variation in snow depth by comparison to within-year spatial variation in snow depth. The results highlight the potential importance of changes in inter-annual variation in snowfall for future shrub photosynthesis in the western Great Basin Desert.

Small-scale patterns in snowmelt timing affect gene flow and the distribution of genetic diversity in the alpine dwarf shrub Salix herbacea

Current threats to biodiversity, such as climate change, are thought to alter the within-species genetic diversity among microhabitats in highly heterogeneous alpine environments. Assessing the spatial organization and dynamics of genetic diversity within species can help to predict the responses of organisms to environmental change. In this study, we evaluated whether small-scale heterogeneity in snowmelt timing restricts gene flow between microhabitats in the common long-lived dwarf shrub Salix herbacea L. We surveyed 273 genets across 12 early- and late-snowmelt sites (that is, ridges and snowbeds) in the Swiss Alps for phenological variation over 2 years and for genetic variation using seven SSR markers. Phenological differentiation triggered by differences in snowmelt timing did not correlate with genetic differentiation between microhabitats. On the contrary, extensive gene flow appeared to occur between microhabitats and slightly less extensively among adjacent mountains. However, ridges exhibited significantly lower levels of genetic diversity than snowbeds, and patterns of effective population size (Ne) and migration (Nem) between microhabitats were strongly asymmetric, with ridges acting as sources and snowbeds as sinks. As no recent genetic bottlenecks were detected in the studied sites, this asymmetry is likely to reflect current meta-population dynamics of the species dominated by gene flow via seeds rather than ancient re-colonization after the last glacial period. Overall, our results suggest that seed dispersal prevents snowmelt-driven genetic isolation, and snowbeds act as sinks of genetic diversity. We discuss the consequences of such small-scale variation in gene flow and diversity levels for population responses to climate change.

Seasonal change in a pollinator community and the maintenance of style length variation in Mertensia fusiformis (Boraginaceae)

BACKGROUND AND AIMS

In sub-alpine habitats, patchiness in snowpack produces marked, small-scale variation in flowering phenology. Plants in early- and late-melting patches are therefore likely to experience very different conditions during their flowering periods. Mertensia fusiformis is an early-flowering perennial that varies conspicuously in style length within and among populations. The hypothesis that style length represents an adaptation to local flowering time was tested. Specifically, it was hypothesized that lower air temperatures and higher frost risk would favour short-styled plants (with stigmas more shielded by corollas) in early-flowering patches, but that the pollen-collecting behaviour of flower visitors in late-flowering patches would favour long-styled plants.

METHODS

Floral morphology was measured, temperatures were monitored and pollinators were observed in several matched pairs of early and late populations. To evaluate effects of cold temperatures on plants of different style lengths, experimental pollinations were conducted during mornings (warm) and evenings (cool), and on flowers that either had or had not experienced a prior frost. The effectiveness of different pollinators was quantified as seed set following single visits to plants with relatively short or long styles.

KEY RESULTS

Late-flowering populations experienced warmer temperatures than early-flowering populations and a different suite of pollinators. Nectar-foraging bumble-bee queens and male solitary bees predominated in early populations, whereas pollen-collecting female solitary bees were more numerous in later sites. Pollinators differed significantly in their abilities to transfer pollen to stigmas at different heights, in accordance with our prediction. However, temperature and frost sensitivity did not differ between long- and short-styled plants. Although plants in late-flowering patches tended to have longer styles than those in early patches, this difference was not consistent.

CONCLUSIONS

Seasonal change in pollinator-mediated selection on style length may help maintain variation in this trait in M. fusiformis, but adaptation to local flowering time is not apparent. The prevalence of short styles in these populations requires further explanation.

Flowering phenology, fruiting success and progressive deterioration of pollination in an early-flowering geophyte

Spatio-temporal patterns of snowmelt and flowering times affect fruiting success in Erythronium grandiflorum Pursh (Liliaceae) in subalpine western Colorado, USA. From 1990 to 1995, I measured the consistency across years of snowmelt patterns and flowering times along a permanent transect. In most years since 1993, I have monitored fruit set in temporal cohorts (early- to late-flowering groups of plants) at one site. To assess 'pollination limitation', I have also conducted supplemental hand-pollination experiments at various times through the blooming season. The onset of blooming is determined by snowmelt, with the earliest years starting a month before the latest years owing to variation in winter snowpack accumulation. Fruit set is diminished or prevented entirely by killing frosts in some years, most frequently but not exclusively for the earlier cohorts. When frosts do not limit fruit set, pollination limitation is frequent, especially in the earlier cohorts. Pollination limitation is strongest for middle cohorts: it tends to be negated by frost in early cohorts and ameliorated by continuing emergence of bumble-bee queens in later cohorts. This lily appears to be poorly synchronized with its pollinators. Across the years of the study, pollination limitation appears to be increasing, perhaps because the synchronization is getting worse.

Changes in snowmelt date and summer precipitation affect the flowering phenology of Erythronium grandiflorum (glacier lily; Liliaceae)

PREMISE OF THE STUDY

Climate change has affected species worldwide, including alterations in phenology, migration patterns, distribution, and survival. Because Erythronium grandiflorum is an early-season bloomer, alterations in its phenology may have serious implications for many North American Rocky Mountain communities, including changes in resource availability for pollinators and herbivores. •

METHODS

We investigated whether changes in the snowmelt date, summer temperature, and summer precipitation have altered the timing and abundance of flowering in E. grandiflorum by collecting long-term data on floral abundance from 1975-2008 in a series of 2 × 2 m plots at the Rocky Mountain Biological Laboratory (RMBL) in Gothic, Colorado in the United States. •

KEY RESULTS

Snowmelt date and mean summer temperature were negatively correlated. Over the 30-yr study, the snowmelt date advanced by 4.14 d/decade, and mean summer temperature increased by 0.38°C/decade. Summer precipitation was variable, showing no change. The first, peak, and last flowering dates of E. grandiflorum advanced an average of 3.2 d/decade. Furthermore, earlier snowmelt and greater summer precipitation in the previous year led to earlier flowering in E. grandiflorum. There was no change in flowering abundance in this species, indicating it may be controlled by a complex set of abiotic and biotic variables. •

CONCLUSIONS

Our study indicates that snowmelt is arriving earlier at the RMBL, which has caused earlier flowering in E. grandiflorum. Because alterations in phenology can disrupt important ecological interactions, information on potential phenological shifts in species that interact with E. grandiflorum is essential in determining the net effect of climate-driven alterations in phenology.

Morphological and genetic variations of Potentilla matsumurae (Rosaceae) between fellfield and snowbed populations

Identifying ecological factors associated with local differentiation of populations is important for understanding microevolutionary processes. Alpine environments offer a unique opportunity to investigate the effects of habitat-specific selective forces and gene flow limitations among populations at a microscale on local adaptation because the heterogeneous snowmelt patterns in alpine ecosystems provide steep environmental changes. We investigated the variation in morphological traits and enzyme loci between fellfield and snowbed populations of Potentilla matsumurae, a common alpine herb with a wide distribution along snowmelt gradients in northern Japan. We found significant differences in morphological traits between fellfield and snowbed habitats in a northern distribution region. These differences were maintained when plants were grown under uniform conditions in a greenhouse. Allozyme variations among 15 populations from geographically separated regions with different historical backgrounds showed that the populations are more genetically differentiated between the fellfield and snowbed habitats within a region than between populations occupying the same habitat type in different regions. These results suggest that variation in snowmelt regimes could be a driving force creating local adaptation and genetic differentiation of alpine plant populations.

The effect of segregation of flowering time on fine-scale spatial genetic structure in an alpine-snowbed herb Primula cuneifolia

The flowering phenology of alpine-snowbed plants varies widely depending on the time of snowmelt. This variation may cause spatial and temporal heterogeneity in pollen dispersal, which in turn may influence genetic structure. We used spatial autocorrelation analyses to evaluate relative effect of segregation in flowering time and physical distance on fine-scale spatial genetic structure (SGS) of a snowbed herb Primula cuneifolia sampled in 10-m grids within a continuous snow patch (110 x 250 m) using nine allozyme loci. Although the individual flower lasts for </=10 days, flowering season varied over 50 days from late June to the middle of August within the plot. The effect of flowering phenology on SGS was assessed using spatial autocorrelation analyses based on the pairwise kinship coefficients for all sampled plants (control pairs), plants with flowering overlap (co-flowering pairs) and plants with separate flowering season (non-co-flowering pairs). The degree of SGS increased as the extent of flowering segregation increased: co-flowering pairs < control pairs < non-co-flowering pairs, indicating substantial effect of restriction in gene flow due to phenological heterogeneity. Flowering segregation caused by snowmelt timing is a critical factor for reinforcing the fine-scale SGS in this species.