Plasticity and local adaptation explain lizard cold tolerance

How does climate variation limit the range of species and what does it take for species to colonize new regions? In this issue of Molecular Ecology, Campbell‐Staton et al. ( 2018 ) address these broad questions by investigating cold tolerance adaptation in the green anole lizard (Anolis carolinensis) across a latitudinal transect. By integrating physiological data, gene expression data and acclimation experiments, the authors disentangle the mechanisms underlying cold adaptation. They first establish that cold tolerance adaptation in Anolis lizards follows the predictions of the oxygen‐ and capacity‐limited thermal tolerance hypothesis, which states that organisms are limited by temperature thresholds at which oxygen supply cannot meet demand. They then explore the drivers of cold tolerance at a finer scale, finding evidence that northern populations are adapted to cooler thermal regimes and that both phenotypic plasticity and heritable genetic variation contribute to cold tolerance. The integration of physiological and gene expression data further highlights the varied mechanisms that drive cold tolerance adaptation in Anolis lizards, including both supply‐side and demand‐side adaptations that improve oxygen economy. Altogether, their work provides new insight into the physiological and genetic mechanisms underlying adaptation to new climatic niches and demonstrates that cold tolerance in northern lizard populations is achieved through the synergy of physiological plasticity and local genetic adaptation for thermal performance.     

Genomic and phenotypic differentiation of Arabidopsis thaliana along altitudinal gradients in the North Italian Alps

Altitudinal gradients in mountain regions are short‐range clines of different environmental parameters such as temperature or radiation. We investigated genomic and phenotypic signatures of adaptation to such gradients in five Arabidopsis thaliana populations from the North Italian Alps that originated from 580 to 2350 m altitude by resequencing pools of 19–29 individuals from each population. The sample includes two pairs of low‐ and high‐altitude populations from two different valleys. High‐altitude populations showed a lower nucleotide diversity and negative Tajima's D values and were more closely related to each other than to low‐altitude populations from the same valley. Despite their close geographic proximity, demographic analysis revealed that low‐ and high‐altitude populations split between 260 000 and 15 000 years before present. Single nucleotide polymorphisms whose allele frequencies were highly differentiated between low‐ and high‐altitude populations identified genomic regions of up to 50 kb length where patterns of genetic diversity are consistent with signatures of local selective sweeps. These regions harbour multiple genes involved in stress response. Variation among populations in two putative adaptive phenotypic traits, frost tolerance and response to light/UV stress was not correlated with altitude. Taken together, the spatial distribution of genetic diversity reflects a potentially adaptive differentiation between low‐ and high‐altitude populations, whereas the phenotypic differentiation in the two traits investigated does not. It may resemble an interaction between adaptation to the local microhabitat and demographic history influenced by historical glaciation cycles, recent seed dispersal and genetic drift in local populations.     

Intraspecific variation in thermal acclimation and tolerance between populations of the winter ant,Prenolepis imparis

Thermal phenotypic plasticity, otherwise known as acclimation, plays an essential role in how organisms respond to short‐term temperature changes. Plasticity buffers the impact of harmful temperature changes; therefore, understanding variation in plasticity in natural populations is crucial for understanding how species will respond to the changing climate. However, very few studies have examined patterns of phenotypic plasticity among populations, especially among ant populations. Considering that this intraspecies variation can provide insight into adaptive variation in populations, the goal of this study was to quantify the short‐term acclimation ability and thermal tolerance of several populations of the winter ant, Prenolepis imparis. We tested for correlations between thermal plasticity and thermal tolerance, elevation, and body size. We characterized the thermal environment both above and below ground for several populations distributed across different elevations within California, USA. In addition, we measured the short‐term acclimation ability and thermal tolerance of those populations. To measure thermal tolerance, we used chill‐coma recovery time (CCRT) and knockdown time as indicators of cold and heat tolerance, respectively. Short‐term phenotypic plasticity was assessed by calculating acclimation capacity using CCRT and knockdown time after exposure to both high and low temperatures. We found that several populations displayed different chill‐coma recovery times and a few displayed different heat knockdown times, and that the acclimation capacities of cold and heat tolerance differed among most populations. The high‐elevation populations displayed increased tolerance to the cold (faster CCRT) and greater plasticity. For high‐temperature tolerance, we found heat tolerance was not associated with altitude; instead, greater tolerance to the heat was correlated with increased plasticity at higher temperatures. These current findings provide insight into thermal adaptation and factors that contribute to phenotypic diversity by revealing physiological variance among populations.     

The genetic architecture of freezing tolerance varies across the range of Arabidopsis thaliana

The capacity to tolerate freezing temperatures limits the geographical distribution of many plants, including several species of agricultural importance. However, the genes involved in freezing tolerance remain largely unknown. Here, we describe the variation in constitutive freezing tolerance that occurs among worldwide accessions of Arabidopsis thaliana. We found that although plants from high latitudes tend to be more freezing tolerant than plants from low latitudes, the environmental factors that shape cold adaptation differ across the species range. Consistent with this, we found that the genetic architecture of freezing tolerance also differs across its range. Conventional genome‐wide association studies helped identify a priori and other promising candidate genes. However, simultaneously modelling climate variables and freezing tolerance together pinpointed other excellent a priori candidate genes. This suggests that if the selective factor underlying phenotypic variation is known, multi‐trait mixed models may aid in identifying the genes that underlie adaptation.     

Local adaptations to frost in marginal and central populations of the dominant forest tree Fagus sylvatica L. as affected by temperature and extreme drought in common garden experiments

Local adaptations to environmental conditions are of high ecological importance as they determine distribution ranges and likely affect species responses to climate change. Increased environmental stress (warming, extreme drought) due to climate change in combination with decreased genetic mixing due to isolation may lead to stronger local adaptations of geographically marginal than central populations. We experimentally observed local adaptations of three marginal and four central populations of Fagus sylvatica L., the dominant native forest tree, to frost over winter and in spring (late frost). We determined frost hardiness of buds and roots by the relative electrolyte leakage in two common garden experiments. The experiment at the cold site included a continuous warming treatment; the experiment at the warm site included a preceding summer drought manipulation. In both experiments, we found evidence for local adaptation to frost, with stronger signs of local adaptation in marginal populations. Winter frost killed many of the potted individuals at the cold site, with higher survival in the warming treatment and in those populations originating from colder environments. However, we found no difference in winter frost tolerance of buds among populations, implying that bud survival was not the main cue for mortality. Bud late frost tolerance in April differed between populations at the warm site, mainly because of phenological differences in bud break. Increased spring frost tolerance of plants which had experienced drought stress in the preceding summer could also be explained by shifts in phenology. Stronger local adaptations to climate in geographically marginal than central populations imply the potential for adaptation to climate at range edges. In times of climate change, however, it needs to be tested whether locally adapted populations at range margins can successfully adapt further to changing conditions.     

The trade-off between cold resistance and growth determines the <Emphasis Type="Italic">Nothofagus pumilio</Emphasis> treeline

The upper and poleward limit of tree distribution are usually determined by abiotic factors such as low temperature and strong winds. Thus, cold resistance is a key element for survival in high altitudes and latitudes where conditions can reduce plant growth. A trade-off between resource allocation to cold resistance and growth could emerge in populations frequently exposed to low temperatures like those in the treeline zone. We studied annual height growth and ice nucleation temperature in Nothofagus pumilio (Nothofagaceae) populations growing in its extremes of altitudinal distribution and in 3 sites situated on a latitudinal gradient in the Chilean Andes. Additionally, gas exchange, water and nitrogen use efficiency and total soluble sugar (TSS) were also measured as possible mechanisms for survival in high altitudes. Individuals from the treeline populations showed lower annual height growth and lower ice nucleation temperatures compared with those from lower populations. In the same way, individuals from more poleward populations showed lower annual height growth and lower ice nucleation temperatures. Gas exchange, water and nitrogen use efficiency and TSS were also higher in the high altitude populations. The results obtained support the hypothesis of trade-off, because the upper and poleward populations showed more cold resistance but a lower height growth. Additionally, we show that cold resistance mechanisms do not impact the physiological performance, suggesting possible adaptation of the high altitude populations. Low temperatures may be affecting cellular growth instead of photosynthesis, creating a pool of carbohydrates that could participate in cold tolerance. Other abiotic and biotic factors should be also assessed to fully understand the distributional range of Nothofagus species.     

Eco‐genomic analysis of the poleward range expansion of the wasp spider Argiope bruennichi shows rapid adaptation and genomic admixture

Poleward range expansions are commonly attributed to global change, but could alternatively be driven by rapid evolutionary adaptation. A well‐documented example of a range expansion during the past decades is provided by the European wasp spider Argiope bruennichi. Using ecological niche modeling, thermal tolerance experiments and a genome‐wide analysis of gene expression divergence, we show that invasive populations have adapted to novel climatic conditions in the course of their expansion. Their climatic niche shift is mirrored in an increased cold tolerance and a population‐specific and functionally differentiated gene expression response. We generated an Argiope reference genome sequence and used population genome resequencing to assess genomic changes associated with the new climatic adaptations. We find clear genetic differentiation and a significant admixture with alleles from East Asian populations in the invasive Northern European populations. Population genetic modeling suggests that at least some of these introgressing alleles have contributed to the new adaptations during the expansion. Our results thus confirm the notion that range expansions are not a simple consequence of climate change, but are accompanied by fast genetic changes and adaptations that may be fuelled through admixture between long separated lineages.     

Indirect selection of thermal tolerance during experimental evolution of Drosophila melanogaster

Natural selection alters the distribution of a trait in a population and indirectly alters the distribution of genetically correlated traits. Long‐standing models of thermal adaptation assume that trade‐offs exist between fitness at different temperatures; however, experimental evolution often fails to reveal such trade‐offs. Here, we show that adaptation to benign temperatures in experimental populations of Drosophila melanogaster resulted in correlated responses at the boundaries of the thermal niche. Specifically, adaptation to fluctuating temperatures (16–25°C) decreased tolerance of extreme heat. Surprisingly, flies adapted to a constant temperature of 25°C had greater cold tolerance than did flies adapted to other thermal conditions, including a constant temperature of 16°C. As our populations were never exposed to extreme temperatures during selection, divergence of thermal tolerance likely reflects indirect selection of standing genetic variation via linkage or pleiotropy. We found no relationship between heat and cold tolerances in these populations. Our results show that the thermal niche evolves by direct and indirect selection, in ways that are more complicated than assumed by theoretical models.     

INDUCER OF CBF EXPRESSION 1 integrates cold signals into FLOWERING LOCUS C‐mediated flowering pathways in Arabidopsis

Plants constantly monitor changes in photoperiod and temperature throughout the year to synchronize flowering with optimal environmental conditions. In the temperate zones, both photoperiod and temperature fluctuate in a somewhat predictable manner through the seasons, although a transient shift to low temperature is also encountered during changing seasons, such as early spring. Although low temperatures are known to delay flowering by inducing the floral repressor FLOWERING LOCUS C (FLC), it is not fully understood how temperature signals are coordinated with photoperiodic signals in the timing of seasonal flowering. Here, we show that the cold signaling activator INDUCER OF CBF EXPRESSION 1 (ICE1), FLC and the floral promoter SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1 (SOC1) constitute an elaborate signaling network that integrates cold signals into flowering pathways. The cold‐activated ICE1 directly induces the gene encoding FLC, which represses SOC1 expression, resulting in delayed flowering. In contrast, under floral promotive conditions, SOC1 inhibits the binding of ICE1 to the promoters of the FLC gene, inducing flowering with a reduction of freezing tolerance. These observations indicate that the ICE1‐FLC‐SOC1 signaling network contributes to the fine‐tuning of flowering during changing seasons.     

Nucleotide diversity patterns of local adaptation at drought‐related candidate genes in wild tomatoes

We surveyed nucleotide diversity at two candidate genes LeNCED1 and pLC30‐15, involved in an ABA (abscisic acid) signalling pathway, in two closely related tomato species Solanum peruvianum and Solanum chilense. Our six population samples (three for each species) cover a range of mesic to very dry habitats. The ABA pathway plays an important role in the plants’ response to drought stress. LeNCED1 is an upstream gene involved in ABA biosynthesis, and pLC30‐15 is a dehydrin gene positioned downstream in the pathway. The two genes show very different patterns of nucleotide variation. LeNCED1 exhibits very low nucleotide diversity relative to the eight neutral reference loci that were previously surveyed in these populations. This suggests that strong purifying selection has been acting on this gene. In contrast, pLC30‐15 exhibits higher levels of nucleotide diversity and, in particular in S. chilense, higher genetic differentiation between populations than the reference loci, which is indicative of local adaptation. In the more drought‐tolerant species S. chilense, one population (from Quicacha) shows a significant haplotype structure, which appears to be the result of positive (diversifying) selection.     

Local adaptation with high gene flow: temperature parameters drive adaptation to altitude in the common frog (Rana temporaria)

Both environmental and genetic influences can result in phenotypic variation. Quantifying the relative contributions of local adaptation and phenotypic plasticity to phenotypes is key to understanding the effect of environmental variation on populations. Identifying the selective pressures that drive divergence is an important, but often lacking, next step. High gene flow between high‐ and low‐altitude common frog (Rana temporaria) breeding sites has previously been demonstrated in Scotland. The aim of this study was to assess whether local adaptation occurs in the face of high gene flow and to identify potential environmental selection pressures that drive adaptation. Phenotypic variation in larval traits was quantified in R. temporaria from paired high‐ and low‐altitude sites using three common temperature treatments. Local adaptation was assessed using Q_ST–F_ST analyses, and quantitative phenotypic divergence was related to environmental parameters using Mantel tests. Although evidence of local adaptation was found for all traits measured, only variation in larval period and growth rate was consistent with adaptation to altitude. Moreover, this was only evident in the three mountains with the highest high‐altitude sites. This variation was correlated with mean summer and winter temperatures, suggesting that temperature parameters are potentially strong selective pressures maintaining local adaptation, despite high gene flow.     

Use of intraspecific variation in thermal responses for estimating an elevational cline in the timing of cold hardening in a sub‐boreal conifer

To avoid winter frost damage, evergreen coniferous species develop cold hardiness with suitable phenology for the local climate regime. Along the elevational gradient, a genetic cline in autumn phenology is often recognised among coniferous populations, but further quantification of evolutionary adaptation related to the local environment and its responsible signals generating the phenological variation are poorly understood. We evaluated the timing of cold hardening among populations of Abies sachalinensis, based on time series freezing tests using trees derived from four seed source populations × three planting sites. Furthermore, we constructed a model to estimate the development of hardening from field temperatures and the intraspecific variations occurring during this process. An elevational cline was detected such that high‐elevation populations developed cold hardiness earlier than low‐elevation populations, representing significant genetic control. Because development occurred earlier at high‐elevation planting sites, the genetic trend across elevation overlapped with the environmental trend. Based on the trade‐off between later hardening to lengthen the active growth period and earlier hardening to avoid frost damage, this genetic cline would be adaptive to the local climate. Our modelling approach estimated intraspecific variation in two model components: the threshold temperature, which was the criterion for determining whether the trees accumulated the thermal value, and the chilling requirement for trees to achieve adequate cold hardiness. A higher threshold temperature and a lower chilling requirement could be responsible for the earlier phenology of the high‐elevation population. These thermal responses may be one of the important factors driving the elevation‐dependent adaptation of A. sachalinensis.     

Genomewide scan for adaptive differentiation along altitudinal gradient in the Andrew's toad Bufo andrewsi

Recent studies of humans, dogs and rodents have started to discover the genetic underpinnings of high altitude adaptations, yet amphibians have received little attention in this respect. To identify possible signatures of adaptation to altitude, we performed a genome scan of 15 557 single nucleotide polymorphisms (SNPs) obtained with restriction site‐associated DNA sequencing of pooled samples from 11 populations of Andrew's toad (Bufo andrewsi) from the edge of the Tibetan Plateau, spanning an altitudinal gradient from 1690 to 2768 m.a.s.l. We discovered significant geographic differentiation among all sites, with an average F_ST = 0.023 across all SNPs. Apart from clear patterns of isolation by distance, we discovered numerous outlier SNPs showing strong associations with variation in altitude (1394 SNPs), average annual temperature (1859 SNPs) or both (1051 SNPs). Levels and patterns of genetic differentiation in these SNPs were consistent with the hypothesis that they have been subject to directional selection and reflect adaptation to altitudinal variation among the study sites. Genes with footprints of selection were significantly enriched in binding and metabolic processes. Several genes potentially related to high altitude adaptation were identified, although the identity and functional significance of most genomic targets of selection remain unknown. In general, the results provide genomic support for results of earlier common garden and low coverage genetic studies that have uncovered substantial adaptive differentiation along altitudinal and latitudinal gradients in amphibians.     

Melanic through nature or nurture: genetic polymorphism and phenotypic plasticity in Harmonia axyridis

Abstract Individuals can adapt to heterogeneity in their environment through either local adaptation or phenotypic plasticity. Colour forms of the ladybird Harmonia axyridis are a classic example of local adaptation, in which the frequency of melanic forms varies greatly between populations. In some populations, there are also large seasonal changes in allele frequency, with melanism being costly in summer and beneficial in winter. We report that the non‐melanic morph of H. axyridis dramatically increases its degree of melanization at cold temperatures. Furthermore, there is genetic variation in reaction norms, with different families responding to temperature in different ways. Variation at different spatial and temporal scales appears to have selected for either genetic or phenotypically plastic adaptations, which may be important in thermoregulation. As melanism is known to have a large effect on fitness in H. axyridis, this plasticity of melanization may have hastened its spread as an invasive species.