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.     

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.     

The phenotypic variance gradient – a novel concept

Evolutionary ecologists commonly use reaction norms, which show the range of phenotypes produced by a set of genotypes exposed to different environments, to quantify the degree of phenotypic variance and the magnitude of plasticity of morphometric and life‐history traits. Significant differences among the values of the slopes of the reaction norms are interpreted as significant differences in phenotypic plasticity, whereas significant differences among phenotypic variances (variance or coefficient of variation) are interpreted as differences in the degree of developmental instability or canalization. We highlight some potential problems with this approach to quantifying phenotypic variance and suggest a novel and more informative way to plot reaction norms: namely “a plot of log (variance) on the y‐axis versus log (mean) on the x‐axis, with a reference line added”. This approach gives an immediate impression of how the degree of phenotypic variance varies across an environmental gradient, taking into account the consequences of the scaling effect of the variance with the mean. The evolutionary implications of the variation in the degree of phenotypic variance, which we call a “phenotypic variance gradient”, are discussed together with its potential interactions with variation in the degree of phenotypic plasticity and canalization.     

Neuroendocrine control of life histories: what do we need to know to understand the evolution of phenotypic plasticity?

Almost all life histories are phenotypically plastic: that is, life-history traits such as timing of breeding, family size or the investment in individual offspring vary with some aspect of the environment, such as temperature or food availability. One approach to understanding this phenotypic plasticity from an evolutionary point of view is to extend the optimality approach to the range of environments experienced by the organism. This approach attempts to understand the value of particular traits in terms of the selection pressures that act on them either directly or owing to trade-offs due to resource allocation and other factors such as predation risk. Because these selection pressures will between environments, the predicted optimal phenotype will too. The relationship expressing the optimal phenotype for different environments is the optimal reaction norm and describes the optimal phenotypic plasticity. However, this view of phenotypic plasticity ignores the fact that the reaction norm must be underlain by some sort of control system: cues about the environment must be collected by sense organs, integrated into a decision about the appropriate life history, and a message sent to the relevant organs to implement that decision. In multicellular animals, this control mechanism is the neuroendocrine system. The central question that this paper addresses is whether the control system affects the reaction norm that evolves. This might happen in two different ways: first, the control system will create constraints on the evolution of reaction norms if it cannot be configured to produce the optimal reaction norm and second, the control system will create additional selection pressures on reaction norms if the neuroendocrine system is costly. If either of these happens, a full understanding of the way in which selection shapes reaction norms must include details of the neuroendocrine control system. This paper presents the conceptual framework needed to explain what is meant by a constraint or cost being created by the neuroendocrine system and discusses the extent to which this occurs and some possible examples. The purpose of doing this is to encourage endocrinologists to take a fresh look at neuroendocrine mechanisms and help identify the properties of the system and situations in which these generate constraints and costs that impinge on the evolution of phenotypic plasticity.     

Stressful environments induce novel phenotypic variation: hierarchical reaction norms for sperm performance of a pervasive invader

Genetic variation for phenotypic plasticity is ubiquitous and important. However, the scale of such variation including the relative variability present in reaction norms among different hierarchies of biological organization (e.g., individuals, populations, and closely related species) is unknown. Complicating interpretation is a trade‐off in environmental scale. As plasticity can only be inferred over the range of environments tested, experiments focusing on fine tuned responses to normal or benign conditions may miss cryptic phenotypic variation expressed under novel or stressful environments. Here, we sought to discern the presence and shape of plasticity in the performance of brown trout sperm as a function of optimal to extremely stressful river pH, and demarcate if the reaction norm varies among genotypes. Our overarching goal was to determine if deteriorating environmental quality increases expressed variation among individuals. A more applied aim was to ascertain whether maintaining sperm performance over a wide pH range could help explain how brown trout are able to invade diverse river systems when transplanted outside of their native range. Individuals differed in their reaction norms of phenotypic expression of an important trait in response to environmental change. Cryptic variation was revealed under stressful conditions, evidenced through increasing among‐individual variability. Importantly, data on population averages masked this variability in plasticity. In addition, canalized reaction norms in sperm swimming velocities of many individuals over a very large range in water chemistry may help explain why brown trout are able to colonize a wide variety of habitats.     

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.     

Quantitative genetics of behavioural reaction norms: genetic correlations between personality and behavioural plasticity vary across stickleback populations

Behavioural ecologists have proposed various evolutionary mechanisms as to why different personality types coexist. Our ability to understand the evolutionary trajectories of personality traits requires insights from the quantitative genetics of behavioural reaction norms. We assayed > 1000 pedigreed stickleback for initial exploration behaviour of a novel environment, and subsequent changes in exploration over a few hours, representing their capacity to adjust their behaviour to changes in perceived novelty and risk. We found heritable variation in both the average level of exploration and behavioural plasticity, and population differences in the sign of the genetic correlation between these two reaction norm components. The phenotypic correlation was not a good indicator of the genetic correlation, implying that quantitative genetics are necessary to appropriately evaluate evolutionary hypotheses in cases such as these. Our findings therefore have important implications for future studies concerning the evolution of personality and plasticity.     

Adaptation to spatially heterogeneous modifying and adaptive environments

The development of an individual's phenotype is influenced by environmental factors (the modifying environment) which may differ from those factors (the adaptive environment) that decide on the adaptational value of the developed phenotype. The shapes of adaptationally optimal norms of reaction are therefore essentially determined by associations between these two environmental components together with the degree of adaptational sensitivity of the developed phenotypes. Two complementary aspects of optimality are accounted for: (a) environments can be optimal for a given norm of reaction and (b) norms of reaction can be optimal for a given environment. The results are obtained for random distribution of genotypes over environmental conditions and under the physiologically reasonable premise that fitness is a function of the costs of modification and adaptation. It turned out that the associations of adaptive and modifying environments are the primary sources of adaptational optimization. More specifically, it is shown that (i) independence between the two environmental components constitutes an adaptationally optimal environment only for norms of reaction in which all phenotypes are adaptively insensitive; (ii) if costs of modification do not depend on the environment, and if the two environmental components are not associated, adaptationally optimal norms of reaction can always be realized through phenogenetic invariance; (iii) as a rule, adaptively sensitive phenotypes developed under strong environmental associations necessitate phenogenetic plasticity for the optimal norm of reaction; (iv) a norm of reaction which is adaptationally optimal in its adaptationally optimal environment can always be realized through phenogenetic invariance, if costs of modification do not vary with the environment. These results reveal an important role of patterns of adaptive sensitivity of phenotypes, which may even surpass that of shapes of norms of reaction in adaptational processes.     

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.     

Variation in reaction norms: Statistical considerations and biological interpretation

Analysis of reaction norms, the functions by which the phenotype produced by a given genotype depends on the environment, is critical to studying many aspects of phenotypic evolution. Different techniques are available for quantifying different aspects of reaction norm variation. We examine what biological inferences can be drawn from some of the more readily applicable analyses for studying reaction norms. We adopt a strongly biologically motivated view, but draw on statistical theory to highlight strengths and drawbacks of different techniques. In particular, consideration of some formal statistical theory leads to revision of some recently, and forcefully, advocated opinions on reaction norm analysis. We clarify what simple analysis of the slope between mean phenotype in two environments can tell us about reaction norms, explore the conditions under which polynomial regression can provide robust inferences about reaction norm shape, and explore how different existing approaches may be used to draw inferences about variation in reaction norm shape. We show how mixed model‐based approaches can provide more robust inferences than more commonly used multistep statistical approaches, and derive new metrics of the relative importance of variation in reaction norm intercepts, slopes, and curvatures.     

Plasticity as a developing trait: exploring the implications

Individual differences in plasticity have been classically framed as genotype-by-environment interactions, with different genotypes showing different reaction norms in response to environmental conditions. However, research has shown that early experience can be a critical factor in shaping an individual's plasticity to later environmental factors. In other words, plasticity itself can be investigated as a developing trait that reflects the combined action of an individual's genes and previous interactions with the environment. In this paper I explore some implications of the idea that the early environment modulates long-term plasticity, with an emphasis on plasticity in behavioral traits. I begin by focusing on the mechanisms that mediate plasticity at the proximate level, and discussing the possibility that some traits may work as generalized mediators of plasticity by affecting the sensitivity of multiple phenol types across developmental contexts. I then tackle the complex problem of the evolution of reaction norms for plasticity. Next, I consider a number of potential implications for research on parental effects and phenotypic matching, and conclude by discussing how plasticity may become a target of evolutionary conflict between parents and offspring. In total, I aim to show how the idea of plasticity as a developing trait offers a rich source of questions and insights that may inform future research in this area.     

Diapause and bet-hedging strategies in insects: a meta-analysis of reaction norm shapes

Many organisms escape from lethal climatological conditions by entering a resistant resting stage called diapause, which needs to be optimally timed with seasonal change. As climate change exerts selection pressure on phenology, the evolution of mean diapause timing, but also of phenotypic plasticity and bet-hedging strategies is expected. The potential of the latter strategy as a means of coping with environmental unpredictability has received little attention in the climate change literature. Populations should be adapted to spatial variation in local conditions; contemporary patterns of phenological strategies across a geographic range may hence provide information about their evolvability. We thus extracted 458 diapause reaction norms from 60 studies. First, we correlated mean diapause timing with mean winter onset. Then we partitioned the reaction norm variance into a temporal component (phenotypic plasticity) and among-offspring variance (diversified bet-hedging) and correlated this variance composition with variability of winter onset. Mean diapause timing correlated reasonably well with mean winter onset, except for populations at high latitudes, which apparently failed to track early onsets. Variance among offspring was, however, limited and correlated only weakly with environmental variability, indicating little scope for bet-hedging. The apparent lack of phenological bet-hedging strategies may pose a risk in a less predictable climate, but we also highlight the need for more data on alternative strategies.     

Natural selection and genetic variation for reproductive reaction norms in a wild bird population

Many morphological and life-history traits show phenotypic plasticity that can be described by reaction norms, but few studies have attempted individual-level analyses of reaction norms in the wild. We analyzed variation in individual reaction norms between laying date and three climatic variables (local temperature, local rainfall, and North Atlantic Oscillation) of 1126 female collared flycatchers (Ficedula albicollis) with a restricted maximum likehood linear mixed model approach using random-effect best linear unbiased predictor estimates for the elevation (i.e., expected laying date in the average environment) and slope (i.e., adjustment in laying date as a function of environment) of females' reaction norms. Variation in laying date was best explained by local temperature, and individual females differed in both the elevation and the slope of their laying date-temperature reaction norms. As revealed by animal model analyses, there was weak evidence for additive genetic variance of elevation (h2 +/- SE = 0.09 +/- 0.09), whereas there was no evidence for heritability of slope (h2 +/- SE = 0.00 +/- 0.01). Selection analysis, using a female's lifetime production of fledglings or recruits as an estimate of her fitness, revealed significant selection for a lower phenotypic value and breeding value for elevation (i.e., earlier laying date at the average temperature). There was selection for steeper phenotypic values of slope (i.e., greater plasticity in the adjustment of laying date to temperature), but no significant selection on the breeding values of slope. Although these results suggest that phenotypic laying date is influenced by additive genetic factors, as well as by an interaction with the environment, selection on plasticity would not produce an evolutionary response.     

Evolution and maintenance of the environmental component of the phenotypic variance: benefit of plastic traits under changing environments

How environmental variances in quantitative traits are influenced by variable environments is an important problem in evolutionary biology. In this study, the evolution and maintenance of phenotypic variance in a plastic trait under stabilizing selection are investigated. The mapping from genotypic value to phenotypic value of the quantitative trait is approximated by a linear reaction norm, with genotypic effects on its phenotypic mean and sensitivity to environment. The environmental deviation is assumed to be decomposed into environmental quality, which interacts with genotypic value, and residual developmental noise, which is independent of genotype. Environmental quality and the optimal phenotype of stabilizing selection are allowed to randomly fluctuate in both space and time, and individuals migrate equally before development and reproduction among different niches. Analyses show that phenotypic plasticity is adaptive within variable environments if correlations have become established between the optimal phenotype and environmental quality in space and/or time. The evolved plasticity increases with variances in optimal phenotypes and correlations between optimal phenotype and environmental quality; this further induces increases in mean fitness and the environmental variance in the trait. Under certain circumstances, however, the environmental variance may decrease with increase in variation in environmental quality.     

Effects of predation environment and food availability on somatic growth in the Livebearing Fish Brachyrhaphis rhabdophora (Pisces: Poeciliidae)

Variation in somatic growth rates is of great interest to biologists because of the relationship between growth and other fitness‐determining traits, and it results from both genetic and environmentally induced variation (i.e. plasticity). Theoretical predictions suggest that mean somatic growth rates and the shape of the reaction norm for growth can be influenced by variation in predator‐induced mortality rates. Few studies have focused on variation in reaction norms for growth in response to resource availability between high‐predation and low‐predation environments. We used juvenile Brachyrhaphis rhabdophora from high‐predation and low‐predation environments to test for variation in mean growth rates and for variation in reaction norms for growth at two levels of food availability in a common‐environment experiment. To test for variation in growth rates in the field, we compared somatic growth rates in juveniles in high‐predation and low‐predation environments. In the common‐environment experiment, mean growth rates did not differ between fish from differing predation environments, but the interaction between predation environment and food level took the form of a crossing reaction norm for both growth in length and mass. Fish from low‐predation environments exhibited no significant difference in growth rate between high and low food treatments. In contrast, fish from high‐predation environments exhibited variation in growth rates between high and low food treatments, with higher food availability resulting in higher growth rates. In the field, individuals in the high‐predation environment grow at a faster rate than those in low‐predation environments at the smallest sizes (comparable to sizes in the common‐environment experiment). These data provide no evidence for evolved differences in mean growth rates between predation environments. However, fish from high‐predation environments exhibited greater plasticity in growth rates in response to resource availability suggesting that predation environments may exhibit increased variation in food availability for prey fish and consequent selection for plasticity.     

Evolutionary changes of nonlinear reaction norms according to thermal adaptation: a comparison of two Drosophila species

While the adaptive significance of discontinuous reaction norms is generally accepted, the evolutionary interpretation of continuous response curves remains speculative, and the occurrence of internal constraints is often suggested as an explanation of experimental observations. In Drosophila melanogaster, various morphometrical traits exhibit convex reaction norms to growth temperature, with a maximum value within the developmental thermal range. We compared a cold-adapted species (D. subobscura) with a mid thermal range at 16 °C, to the warm-adapted D. melanogaster (mid thermal range at 22 °C) for three different morphometrical traits: wing and thorax length in both sexes and ovariole number in females. Maximum value temperatures were ordered in the same way for the three traits in both species: ovariole number > thorax length > wing length. Significant differences were also observed between the two species for the curvature parameter of the quadratic adjustment. The major observation was a significant lateral shift in the reaction norms: maximum values were observed at much lower temperatures in the cold-adapted species than in the warm-adapted one. The parallelism between mid thermal range variation and the position of the maximum value strongly suggests an adaptive displacement of the response curves. Natural selection may thus act not only on trait mean values but also on phenotypic plasticity and on the shape of reaction norms.     

PLASTICITY VERSUS ENVIRONMENTAL CANALIZATION: POPULATION DIFFERENCES IN THERMAL RESPONSES ALONG A LATITUDINAL GRADIENT IN DROSOPHILA SERRATA

The phenotypic plasticity of traits, defined as the ability of a genotype to express different phenotypic values of the trait across a range of environments, can vary between habitats depending on levels of temporal and spatial heterogeneity. Other traits can be insensitive to environmental perturbations and show environmental canalization. We tested levels of phenotypic plasticity in diverse Drosophila serrata populations along a latitudinal cline ranging from a temperate, variable climate to a tropical, stable climate by measuring developmental rate and size-related traits at three temperatures (16°C, 22°C, and 28°C). We then compared the slopes of the thermal reaction norms among populations. The 16–22°C part of the reaction norms for developmental rate was flatter (more canalized) for the temperate populations than for the tropical populations. However, slopes for the reaction norms of the two morphological traits (wing size, wing:thorax ratio), were steeper (more plastic) in the temperate versus the tropical populations over the entire thermal range. The different latitudinal patterns in plasticity for developmental rate and the morphological traits may reflect contrasting selection pressures along the tropical–temperate thermal gradient.