When three traits make a line: evolution of phenotypic plasticity and genetic assimilation through linear reaction norms in stochastic environments

Genetic assimilation emerges from selection on phenotypic plasticity. Yet, commonly used quantitative genetics models of linear reaction norms considering intercept and slope as traits do not mimic the full process of genetic assimilation. We argue that intercept–slope reaction norm models are insufficient representations of genetic effects on linear reaction norms and that considering reaction norm intercept as a trait is unfortunate because the definition of this trait relates to a specific environmental value (zero) and confounds genetic effects on reaction norm elevation with genetic effects on environmental perception. Instead, we suggest a model with three traits representing genetic effects that, respectively, (i) are independent of the environment, (ii) alter the sensitivity of the phenotype to the environment and (iii) determine how the organism perceives the environment. The model predicts that, given sufficient additive genetic variation in environmental perception, the environmental value at which reaction norms tend to cross will respond rapidly to selection after an abrupt environmental change, and eventually becomes equal to the new mean environment. This readjustment of the zone of canalization becomes completed without changes in genetic correlations, genetic drift or imposing any fitness costs of maintaining plasticity. The asymptotic evolutionary outcome of this three‐trait linear reaction norm generally entails a lower degree of phenotypic plasticity than the two‐trait model, and maximum expected fitness does not occur at the mean trait values in the population.     

Selection on plasticity of seasonal life-history traits using random regression mixed model analysis

Theory considers the covariation of seasonal life-history traits as an optimal reaction norm, implying that deviating from this reaction norm reduces fitness. However, the estimation of reaction-norm properties (i.e., elevation, linear slope, and higher order slope terms) and the selection on these is statistically challenging. We here advocate the use of random regression mixed models to estimate reaction-norm properties and the use of bivariate random regression to estimate selection on these properties within a single model. We illustrate the approach by random regression mixed models on 1115 observations of clutch sizes and laying dates of 361 female Ural owl Strix uralensis collected over 31 years to show that (1) there is variation across individuals in the slope of their clutch size-laying date relationship, and that (2) there is selection on the slope of the reaction norm between these two traits. Hence, natural selection potentially drives the negative covariance in clutch size and laying date in this species. The random-regression approach is hampered by inability to estimate nonlinear selection, but avoids a number of disadvantages (stats-on-stats, connecting reaction-norm properties to fitness). The approach is of value in describing and studying selection on behavioral reaction norms (behavioral syndromes) or life-history reaction norms. The approach can also be extended to consider the genetic underpinning of reaction-norm properties.     

Phenological mismatch drives selection on elevation,but not on slope,of breeding time plasticity in a wild songbird

Phenotypic plasticity is an important mechanism for populations to respond to fluctuating environments, yet may be insufficient to adapt to a directionally changing environment. To study whether plasticity can evolve under current climate change, we quantified selection and genetic variation in both the elevation (RN_E) and slope (RN_S) of the breeding time reaction norm in a long‐term (1973–2016) study population of great tits (Parus major). The optimal RN_E (the caterpillar biomass peak date regressed against the temperature used as cue by great tits) changed over time, whereas the optimal RN_S did not. Concordantly, we found strong directional selection on RN_E, but not RN_S, of egg‐laying date in the second third of the study period; this selection subsequently waned, potentially due to increased between‐year variability in optimal laying dates. We found individual and additive genetic variation in RN_E but, contrary to previous studies on our population, not in RN_S. The predicted and observed evolutionary change in RN_E was, however, marginal, due to low heritability and the sex limitation of laying date. We conclude that adaptation to climate change can only occur via micro‐evolution of RN_E, but this will necessarily be slow and potentially hampered by increased variability in phenotypic optima.     

Adaptation to an extraordinary environment by evolution of phenotypic plasticity and genetic assimilation

Adaptation to a sudden extreme change in environment, beyond the usual range of background environmental fluctuations, is analysed using a quantitative genetic model of phenotypic plasticity. Generations are discrete, with time lag τ between a critical period for environmental influence on individual development and natural selection on adult phenotypes. The optimum phenotype, and genotypic norms of reaction, are linear functions of the environment. Reaction norm elevation and slope (plasticity) vary among genotypes. Initially, in the average background environment, the character is canalized with minimum genetic and phenotypic variance, and no correlation between reaction norm elevation and slope. The optimal plasticity is proportional to the predictability of environmental fluctuations over time lag τ. During the first generation in the new environment the mean fitness suddenly drops and the mean phenotype jumps towards the new optimum phenotype by plasticity. Subsequent adaptation occurs in two phases. Rapid evolution of increased plasticity allows the mean phenotype to closely approach the new optimum. The new phenotype then undergoes slow genetic assimilation, with reduction in plasticity compensated by genetic evolution of reaction norm elevation in the original environment.     

Natural selection and inheritance of breeding time and clutch size in the collared flycatcher

Many characteristics of organisms in free-living populations appear to be under directional selection, possess additive genetic variance, and yet show no evolutionary response to selection. Avian breeding time and clutch size are often-cited examples of such characters. We report analyses of inheritance of, and selection on, these traits in a long-term study of a wild population of the collared flycatcher Ficedula albicollis. We used mixed model analysis with REML estimation ("animal models") to make full use of the information in complex multigenerational pedigrees. Heritability of laying date, but not clutch size, was lower than that estimated previously using parent-offspring regressions, although for both traits there was evidence of substantial additive genetic variance (h2 = 0.19 and 0.29, respectively). Laying date and clutch size were negatively genetically correlated (rA = -0.41 +/- 0.09), implying that selection on one of the traits would cause a correlated response in the other, but there was little evidence to suggest that evolution of either trait would be constrained by correlations with other phenotypic characters. Analysis of selection on these traits in females revealed consistent strong directional fecundity selection for earlier breeding at the level of the phenotype (beta = -0.28 +/- 0.03), but little evidence for stabilising selection on breeding time. We found no evidence that clutch size was independently under selection. Analysis of fecundity selection on breeding values for laying date, estimated from an animal model, indicated that selection acts directly on additive genetic variance underlying breeding time (beta = -0.20 +/- 0.04), but not on clutch size (beta = 0.03 +/- 0.05). In contrast, selection on laying date via adult female survival fluctuated in sign between years, and was opposite in sign for selection on phenotypes (negative) and breeding values (positive). Our data thus suggest that any evolutionary response to selection on laying date is partially constrained by underlying life-history trade-offs, and illustrate the difficulties in using purely phenotypic measures and incomplete fitness estimates to assess evolution of life-history trade-offs. We discuss some of the difficulties associated with understanding the evolution of laying date and clutch size in natural populations.     

Plasticity results in delayed breeding in a long‐distant migrant seabird

A major question for conservationists and evolutionary biologists is whether natural populations can adapt to rapid environmental change through micro‐evolution or phenotypic plasticity. Making use of 17 years of data from a colony of a long‐distant migratory seabird, the common tern (Sterna hirundo), we examined phenotypic plasticity and the evolutionary potential of breeding phenology, a key reproductive trait. We found that laying date was strongly heritable (0.27 ± 0.09) and under significant fecundity selection for earlier laying. Paradoxically, and in contrast to patterns observed in most songbird populations, laying date became delayed over the study period, by about 5 days. The discrepancy between the observed changes and those predicted from selection on laying date was explained by substantial phenotypic plasticity. The plastic response in laying date did not vary significantly among individuals. Exploration of climatic factors showed individual responses to the mean sea surface temperature in Senegal in December prior to breeding: Common terns laid later following warmer winters in Senegal. For each 1°C of warming of the sea surface in Senegal, common terns delayed their laying date in northern Germany by 6.7 days. This suggests that warmer waters provide poorer wintering resources. We therefore found that substantial plastic response to wintering conditions can oppose natural selection, perhaps constraining adaptation.     

The genetics of phenotypic plasticity. VII. Evolution in a spatially-structured environment

Previous models of the evolution of phenotypic plasticity have, for the most part, not considered the effects of genetic architecture and spatial structure. I examine those factors with an individual-based simulation model. With regard to genetic architecture, I considered how the presence of different types of loci would affect medium-term evolutionary outcomes. The types of loci differed in how the environment determined phenotypic expression and included loci that were insensitive to the environment (non-plastic loci), sensitive in a linear fashion, and sensitive in a quadratic fashion (both plastic loci). With regard to spatial structure, I investigated the affects of migration patterns. These simulations demonstrated that two general conditions are necessary for phenotypic plasticity to be selected. (1) The environment must have a strong influence on genotypic expression. (2) The between-generation changes in the environment must be large and predictable, in the current instance because of migration in a spatially-structured (clinal) environment. Responses to selection were not simple, however. Rarely were pure strategies — genetic specialization or phenotypic plasticity — selected for. Instead, the existence of multiple types of loci led to mixed genetic outcomes. The result of this mixed outcome were individuals with reaction norms that were less steep than the optimal reaction norm (when non-plastic and linear-plastic loci were present) or individuals with curved reaction norms when the optimal reaction norms was linear (when all three types of loci were present). A pure plasticity strategy had the highest global fitness because plastic individuals would match the optimal phenotype everywhere. The reason that the metapopulation did not achieve this global fitness optimum is that local selection is stronger than global selection. Each deme is driven to a local fitness peak based on the combined, locally additive effects of the non-plastic and plastic loci. Plasticity is only selected globally, so plasticity becomes more highly favored with high migration rates. This effect was greatest in parts of the cline where the plasticity loci were not being expressed and, thus, not locally selected upon. That is, in these demes local selection was weak or absent allowing global fitness effects to predominate.     

Ontogenetic changes in genetic variances of age-dependent plasticity along a latitudinal gradient

The expression of phenotypic plasticity may differ among life stages of the same organism. Age-dependent plasticity can be important for adaptation to heterogeneous environments, but this has only recently been recognized. Whether age-dependent plasticity is a common outcome of local adaptation and whether populations harbor genetic variation in this respect remains largely unknown. To answer these questions, we estimated levels of additive genetic variation in age-dependent plasticity in six species of damselflies sampled from 18 populations along a latitudinal gradient spanning 3600 km. We reared full sib larvae at three temperatures and estimated genetic variances in the height and slope of thermal reaction norms of body size at three points in time during ontogeny using random regression. Our data show that most populations harbor genetic variation in growth rate (reaction norm height) in all ontogenetic stages, but only some populations and ontogenetic stages were found to harbor genetic variation in thermal plasticity (reaction norm slope). Genetic variances in reaction norm height differed among species, while genetic variances in reaction norm slope differed among populations. The slope of the ontogenetic trend in genetic variances of both reaction norm height and slope increased with latitude. We propose that differences in genetic variances reflect temporal and spatial variation in the strength and direction of natural selection on growth trajectories and age-dependent plasticity. Selection on age-dependent plasticity may depend on the interaction between temperature seasonality and time constraints associated with variation in life history traits such as generation length.     

Phenotypic plasticity in morphological traits in two populations of Drosophila melanogaster

Isofemale lines of two populations of Drosophila melanogaster, originating from France and Tanzania, were examined over a range of temperatures. Morphological traits showed distinct patterns in phenotypic plasticity; flies of the two populations differed in shape. Genotype-by-Environment (G*E) interactions were frequently found in the Tanzania population, but were hardly present in the France population. If G*E interaction was present over temperature, estimates of additive genetic variance and additive genetic covariance were made to compare theoretical models with our data. The conclusion is that in France Drosophila melanogaster has been selected over a wider range of temperatures, resulting in parallel reaction norms of more optimal slope. In contrast, selection must have taken place over a narrower temperature range in Tanzanian flies, and will have exerted no direct influence on the slope of the reaction norm.     

A modular concept of phenotypic plasticity in plants

Based on empirical evidence from the literature we propose that, in nature, phenotypic plasticity in plants is usually expressed at a subindividual level. While reaction norms (i.e. the type and the degree of plant responses to environmental variation) are a property of genotypes, they are expressed at the level of modular subunits in most plants. We thus contend that phenotypic plasticity is not a whole-plant response, but a property of individual meristems, leaves, branches and roots, triggered by local environmental conditions. Communication and behavioural integration of interconnected modules can change the local responses in different ways: it may enhance or diminish local plastic effects, thereby increasing or decreasing the differences between integrated modules exposed to different conditions. Modular integration can also induce qualitatively different responses, which are not expressed if all modules experience the same conditions. We propose that the response of a plant to its environment is the sum of all modular responses to their local conditions plus all interaction effects that are due to integration. The local response rules to environmental variation, and the modular interaction rules may be seen as evolving traits targeted by natural selection. Following this notion, whole-plant reaction norms are an integrative by-product of modular plasticity, which has far-reaching methodological, ecological and evolutionary implications.     

CONTRASTING PATTERNS OF PHENOTYPIC PLASTICITY IN REPRODUCTIVE TRAITS IN TWO GREAT TIT (PARUS MAJOR) POPULATIONS

Phenotypic plasticity is an important mechanism via which populations can respond to changing environmental conditions, but we know very little about how natural populations vary with respect to plasticity. Here we use random‐regression animal models to understand the multivariate phenotypic and genetic patterns of plasticity variation in two key life‐history traits, laying date and clutch size, using data from long‐term studies of great tits in The Netherlands (Hoge Veluwe [HV]) and UK (Wytham Woods [WW]). We show that, while population‐level responses of laying date and clutch size to temperature were similar in the two populations, between‐individual variation in plasticity differed markedly. Both populations showed significant variation in phenotypic plasticity (IxE) for laying date, but IxE was significantly higher in HV than in WW. There were no significant genotype‐by‐environment interactions (GxE) for laying date, yet differences in GxE were marginally nonsignificant between HV and WW. For clutch size, we only found significant IxE and GxE in WW but no significant difference between populations. From a multivariate perspective, plasticity in laying date was not correlated with plasticity in clutch size in either population. Our results suggest that generalizations about the form and cause of any response to changing environmental conditions across populations may be difficult.     

Genotype by environment interaction for litter size in pigs as quantified by reaction norms analysis

A Bayesian procedure was used to estimate linear reaction norms (i.e. individual G × E plots) on 297 518 litter size records of 121 104 sows, daughters of 2040 sires, recorded on 144 farms in North and Latin America, Europe, Asia and Australia. The method allowed for simultaneous estimation of all parameters involved. The analysis was carried out on three subsets, comprising (i) parity 1 records of 33 641 sows of line B, (ii) all parity records of 52 120 sows of line B and (iii) all parity records of 121 104 sows of lines A, B and A × B. Estimated heritabilities ranged from 0.09 to 0.10 (smallest to largest subset) for the intercept of the reaction norms, and were 0.15, 0.08 and 0.02 (ditto) for the slope. Estimated genetic correlations between intercept and slope were -0.09, +0.26 and +0.69 (ditto). The three subsets therefore showed a progressively lower genetic component to environmental sensitivity, and progressively less re-ranking of genotypes across the environmental (herd-year-season) range. In a genetic evaluation that does not include reaction norms in the statistical model, part of the G × E effect remains confounded with the additive genetic effect, which may lead to errors in the estimates of the additive genetic effect; the reaction norms model removes this confounding. The intercept estimates from the largest data subset show correlations with litter size estimated breeding values (EBV) from routine genetic evaluation (without reaction norms included) of 0.78 to 0.85 for sows with one to seven litter records, and 0.75 for sires. Hence, including reaction norms in genetic evaluation would increase the reliability of the EBV of young selection candidates without own performance or progeny data by considerably more than 100 × (1/0.75-1) = 33%. Reaction norm slope estimates turn out to be very demanding statistics; environmental sensitivity must therefore be classified as a 'hard-to-measure' trait.     

Reaction norms of Arabidopsis. I. Plasticity of characters and correlations across water,nutrient and light gradients

The univariate and multivariate study of variation for phenotypic plasticity is central to providing a clear understanding of hypotheses about the genetic control and evolution of reaction norms in natural populations. Arabidopsis thaliana is an ideal organism for the study of Genotype × Environment interactions (i.e., genetic variation for plasticity), because of the ease with which it can be grown in large numbers and due to the amount of information already available on its genetics, physiology and developmental biology. In this paper, we report on the plasticity, genetic variation and G × E interactions of four populations of A. thaliana in response to three environmental gradients (water, light and nutrients), each characterized by four levels of the controlled parameter. We measured nine traits and obtained their reaction norms. Path analysis was used to study the plasticity of character correlations. We found a tendency for A. thaliana reaction norms to be linear (either flat, i.e. no plasticity, or with a significant slope), in accordance with previous studies. We detected substantial amounts of genetic variation for plasticity in the light and nutrient gradients, but not in the water gradient. Dramatic restructuring of character correlations was induced by changes in environmental conditions, although some paths tended to be stable irrespective of the environment, thereby suggesting some degree of canalization.     

Testing Local Host Adaptation and Phenotypic Plasticity in a Herbivore When Alternative Related Host Plants Occur Sympatrically

Host race formation in phytophagous insects can be an early stage of adaptive speciation. However, the evolution of phenotypic plasticity in host use is another possible outcome. Using a reciprocal transplant experiment we tested the hypothesis of local adaptation in the aphid Brevicoryne brassicae. Aphid genotypes derived from two sympatric host plants, Brassica oleracea and B. campestris, were assessed in order to measure the extent of phenotypic plasticity in morphological and life history traits in relation to the host plants. We obtained an index of phenotypic plasticity for each genotype. Morphological variation of aphids was summarized by principal components analysis. Significant effects of recipient host on morphological variation and life history traits (establishment, age at first reproduction, number of nymphs, and intrinsic growth rate) were detected. We did not detected genotype × host plant interaction; in general the genotypes developed better on B. campestris, independent of the host plant species from which they were collected. Therefore, there was no evidence to suggest local adaptation. Regarding plasticity, significant differences among genotypes in the index of plasticity were detected. Furthermore, significant selection on PC1 (general aphid body size) on B. campestris, and on PC1 and PC2 (body length relative to body size) on B. oleracea was detected. The elevation of the reaction norm of PC1 and the slope of the reaction norm for PC2 (i.e., plasticity) were under directional selection. Thus, host plant species constitute distinct selective environments for B. brassicae. Aphid genotypes expressed different phenotypes in response to the host plant with low or nil fitness costs. Phenotypic plasticity and gene flow limits natural selection for host specialization promoting the maintenance of genetic variation in host exploitation.     

Variation in plant thermal reaction norms along a latitudinal gradient – more than adaptation to season length

Little is known about the extent to which observed phenological responses to changes in climate are the result of phenotypic plasticity or genetic changes. We also know little about how plasticity, in terms of thermal reaction norms, vary spatially. We investigated if the thermal reaction norms for flower development of five crucifer species (Brassicaceae) differed among three regions along a south–north latitudinal gradient in replicated experiments. The mean response (elevation) of thermal reaction norms of flowering differed among regions in all study species, while sensitivity of flower development to temperature (slope) differed in only one of the species. Differences in mean responses corresponded to cogradient patterns in some species, but countergradient patterns in other. This suggests that differences among regions were not solely the result of adaptation to differences in the length of the growing season, but that other factors, such as herbivory, play an important role. Differences in development rate within species were mainly explained by variation in early phases of bud formation in some species but by variation in later phases of bud formation in other species. The differences in latitudinal patterns of thermal reaction norms among species observed in this study are important, both to identify agents of selection and to predict short‐ and long‐term responses to a warming climate.     

Developmental constraints on an adaptive plasticity: reaction norms of pigmentation in adult segments of Drosophila melanogaster

SUMMARY Variation of dark pigmentation according to developmental temperature was investigated in two geographic populations (France and India) with the isofemale line technique (20 lines for each population). The response curves called the reaction norms, were established in females for seven different segments: the mesothorax and abdomen segments 2–7 (Abd 2–7). In all cases the response curves were non-linear and had to be described either by a quadratic convex polynomial for thorax and Abd 2–5, or by a cubic polynomial for Abd 6 and 7. Among abdomen segments, increasing antero-posterior gradients were observed for several traits, including average pigmentation, overall phenotypic plasticity, the temperature of minimum pigmentation, and the curvature parameter of quadratic norms. Genetic correlations between abdomen segments were high when adjacent segments were considered, but became nil when more distant segments were correlated, suggesting that different pigmentation genes are expressed in the anterior and the posterior part of the abdomen. Characteristic values of reaction norms provided information either on trait value (i.e., the extension of pigmentation) or on plasticity. Correlations between plasticity and pigmentation were generally low and non-significant, suggesting their genetic independence. The overall darker pigmentation which is observed at low temperatures is assumed to be an adaptive plasticity. However, the differences which are evidenced among segments reveal strong interactions with developmental genes. These interactions are less likely to be a consequence of natural selection and are better interpreted as developmental constraints. The reaction norms analysis reveals the complexity of these interactions and should help, in the future, in the identification of the responsible thermosensitive genes.     

Plasticity in laying dates of Canada Geese in response to spring phenology

Observed phenological changes can be explained either by individual phenotypic plasticity or by evolutionary changes, but there is more evidence pointing towards phenotypic plasticity to explain the mechanism behind changes in bird phenology. However, most studies on phenology have been conducted on insectivorous bird species for which breeding is closely tied to temperature and insect emergence. In this study, we examined the consequences of climatic conditions on the nesting phenology of temperate breeding Canada Geese Branta canadensis maxima, which rely on a continuous food supply, during a 14‐year period (2003–16). We determined whether laying dates were plastically adjusted to spring environmental conditions, and whether this adjustment resulted in a laying date advancement. We further estimated the strength and shape of selection acting on breeding timing, by looking at the effect of laying date on the relative number of young successfully hatched in a nest. We found that Geese plastically adjusted their laying date to spring maximum temperature (and not to precipitation or ice break‐up), resulting in a 9‐day advancement of laying date in the population for that period. Laying date was also moderately repeatable (r = 0.23) and subject to directional selection, but stabilizing selection was negligible. We thus demonstrate how Canada Geese plastically adjust laying dates to temperature, which may further be beneficial to nesting success. Evolutionary change of laying date to selection related to climate change, however, is still possible.     

Old wine in new bottles: reaction norms in salmonid fishes

Genetic variability in reaction norms reflects differences in the ability of individuals, populations and ultimately species to respond to environmental change. By increasing our understanding of how genotype × environment interactions influence evolution, studies of genetic variation in phenotypic plasticity serve to refine our capacity to predict how populations will respond to natural and anthropogenic environmental variability, including climate change. Given the extraordinary variability in morphology, behaviour and life history in salmonids, one might anticipate the research milieu on reaction norms in these fishes to be empirically rich and intellectually engaging. Here, I undertake a review of genetic variability in continuous and discontinuous (threshold) norms of reaction in salmonid fishes, as determined primarily (but not exclusively) by common-garden experiments. Although in its infancy from a numerical publication perspective, there is taxonomically broad evidence of genetic differentiation in continuous, threshold and bivariate reaction norms among individuals, families and populations (including inter-population hybrids and backcrosses) for traits as divergent as embryonic development, age and size at maturity, and gene expression. There is compelling inferential evidence that plasticity is heritable and that population differences in reaction norms can reflect adaptive responses, by natural selection, to local environments. As a stimulus for future work, a series of 20 research questions are identified that focus on reaction-norm variability, selection, costs and constraints, demographic and conservation consequences, and genetic markers and correlates of phenotypic plasticity.     

Variation in phenotypic plasticity and selection patterns in blue tit breeding time: between- and within-population comparisons

1.?Phenotypic plasticity, the response of individual phenotypes to their environment, can allow organisms to cope with spatio-temporal variation in environmental conditions. Recent studies have shown that variation exists among individuals in their capacity to adjust their traits to environmental changes and that this individual plasticity can be under strong selection. Yet, little is known on the extent and ultimate causes of variation between populations and individuals in plasticity patterns. 2.?In passerines, timing of breeding is a key life-history trait strongly related to fitness and is known to vary with the environment, but few studies have investigated the within-species variation in individual plasticity. 3.?Here, we studied between- and within-population variation in breeding time, phenotypic plasticity and selection patterns for this trait in four Mediterranean populations of blue tits (Cyanistes caeruleus) breeding in habitats varying in structure and quality. 4.?Although there was no significant warming over the course of the study, we found evidence for earlier onset of breeding in warmer years in all populations, with reduced plasticity in the less predictable environment. In two of four populations, there was significant inter-individual variation in plasticity for laying date. Interestingly, selection for earlier laying date was significant only in populations where there was no inter-individual differences in plasticity. 5.?Our results show that generalization of plasticity patterns among populations of the same species might be challenging even at a small spatial scale and that the amount of within-individual variation in phenotypic plasticity may be linked to selective pressures acting on these phenotypic traits.     

Habitat-specific differences in thermal plasticity in natural populations of a soil arthropod

When populations experience substantial variation in environmental conditions, they may evolve phenotypic plasticity in response to these varying selection pressures. Evolutionary theory predicts differentiation in the level of phenotypic plasticity among different habitats. We evaluated temperature-induced phenotypic responses in juvenile growth rate in natural populations of the springtail Orchesella cincta , inhabiting forest and heathland. These habitats typically co-occur but differ strongly with respect to, for example, thermal regime, relative humidity, and structure. Offspring of females from the two habitats were reared at different temperatures in climate rooms and the temperature response of juvenile growth rate and egg size was measured. We found a habitat-specific difference in plasticity of juvenile growth rate. The reaction norms of the forest populations were steeper than the reaction norms for heath populations at two replicated sampling sites. Egg weight itself was demonstrated to be a plastic trait with a higher egg weight at low temperatures, but the thermal response did not differ between habitats. We conclude that these populations have diverged due to strong local natural selection. Our results support the argument that the level of phenotypic plasticity itself can be under selection and that differentiation in reaction norms can occur even in neighbouring habitats with no barrier to gene flow.  © 2008 The Linnean Society of London, Biological Journal of the Linnean Society , 2008, 94 , 265–271.