Probabilistic maturation reaction norms assessed from mark–recaptures of wild fish in their natural habitat

Reaction norms are a valuable tool in evolutionary biology. Lately, the probabilistic maturation reaction norm approach, describing probabilities of maturing at combinations of age and body size, has been much applied for testing whether phenotypic changes in exploited populations of fish are mainly plastic or involving an evolutionary component. However, due to typical field data limitations, with imperfect knowledge about individual life histories, this demographic method still needs to be assessed. Using 13 years of direct mark–recapture observations on individual growth and maturation in an intensively sampled population of brown trout (Salmo trutta), we show that the probabilistic maturation reaction norm approach may perform well even if the assumption of equal survival of juvenile and maturing fish does not hold. Earlier studies have pointed out that growth effects may confound the interpretation of shifts in maturation reaction norms, because this method in its basic form deals with body size rather than growth. In our case, however, we found that juvenile body size, rather than annual growth, was more strongly associated with maturation. Viewed against earlier studies, our results also underscore the challenges of generalizing life‐history patterns among species and populations.     

DOES SIZE MATTER MOST? THE EFFECT OF GROWTH HISTORY ON PROBABILISTIC REACTION NORM FOR SALMON MATURATION

Body size is widely believed to affect the occurrence of sexual maturation. Recent studies have used changes in the age-specific body size at which the probability of maturing is 50%, a feature of probabilistic reaction norms, to quantify purported evolution of life histories. However, body size results from a combination of growth rates during successive developmental stages. Therefore, to understand the evolution of the maturation schedule, it is necessary to comprehend the relationships among body size, growth history, and maturation schedule. We examined the relationships among body size, previous growth history, and maturation probability in chum salmon (Oncorhynchus keta). In this study, previous growth history was estimated from yearly specific growth increments that provide information describing body size. Previous growth history was found to be more closely linked to maturation probability than body size. The most recent growth condition was the most important factor affecting whether a fish matured during the subsequent breeding season. Because individuals of similar body size and same age can have different growth histories, the relationship between body size and maturation probability could be plastically modified by growth history. This may violate an assumption required to infer evolution, namely that size-related maturation trends in probabilistic reaction norms are immune to growth history.     

Phenotypic convergence along a gradient of predation risk

A long-standing question in ecology is whether phenotypic plasticity, rather than selection per se, is responsible for phenotypic variation among populations. Plasticity can increase or decrease variation, but most previous studies have been limited to single populations, single traits and a small number of environments assessed using univariate reaction norms. Here, examining two genetically distinct populations of Daphnia pulex with different predation histories, we quantified predator-induced plasticity among 11 traits along a fine-scale gradient of predation risk by a predator (Chaoborus) common to both populations. We test the hypothesis that plasticity can be responsible for convergence in phenotypes among different populations by experimentally characterizing multivariate reaction norms with phenotypic trajectory analysis (PTA). Univariate analyses showed that all genotypes increased age and size at maturity, and invested in defensive spikes (neckteeth), but failed to quantitatively describe whole-organism response. In contrast, PTA quantified and qualified the phenotypic strategy the organism mobilized against the selection pressure. We demonstrate, at the whole-organism level, that the two populations occupy different areas of phenotypic space in the absence of predation but converge in phenotypic space as predation threat increases.     

Reduction of Delia radicum attack in field brassicas using a vertical barrier

According to life‐history theory, longer development time may result in bigger adults. However, reaction norms describing age and size at maturity often follow an L‐shaped form. This relationship is attributable to the simple notion that slowly growing individuals may not lengthen their development excessively after the maturation decision has been made, for example, when development is time limited in seasonal environments. In arthropods, growth occurs within instars, and thus the optimal growth strategy might be mediated by the phenotypic adjustment of instar numbers. We studied the relationship between age and size at maturity of a lichen‐feeding moth, Eilema depressum (Esper) (Lepidoptera: Arctiidae: Lithosiinae), and the variability of instar numbers in relation to achieved adult body mass and time used for maturation. A positive relationship between age and size at maturity was found across developmental pathways and a negative one within the developmental pathways. Directly developing larvae had higher growth rates, attained smaller pupal mass, and passed fewer instars than larvae maturing after overwintering. Host quality did not affect whether larvae matured during the remaining or the next season. High variation in the number of instars together with variable growth rates indicates high plasticity in adaptation to varying environmental conditions. Our results also confirm previous results that instar number variability may be a key characteristic mediating age and size at maturity in insects.     

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.     

Impact of environmental covariation in growth and mortality on evolving maturation reaction norms

Maturation age and size have important fitness consequences through their effects on survival probabilities and body sizes. The evolution of maturation reaction norms in response to environmental covariation in growth and mortality is therefore a key subject of life-history theory. The eco-evolutionary model we present and analyze here incorporates critical features that earlier studies of evolving maturation reaction norms have often neglected: the trade-off between growth and reproduction, source-sink population structure, and population regulation through density-dependent growth and fecundity. We report the following findings. First, the evolutionarily optimal age at maturation can be decomposed into the sum of a density-dependent and a density-independent component. These components measure, respectively, the hypothetical negative age at which an individual's length would be 0 and the delay in maturation relative to this offset. Second, along any growth trajectory, individuals mature earlier when mortality is higher. This allows us to deduce, third, how the shapes of evolutionarily optimal maturation reaction norms depend on the covariation between growth and mortality (positive or negative, linear or curvilinear, and deterministic or probabilistic). Providing eco-evolutionary explanations for many alternative reaction-norm shapes, our results appear to be in good agreement with current empirical knowledge on maturation dynamics.     

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.     

THE EVOLUTION OF PHENOTYPIC PLASTICITY IN LIFE-HISTORY TRAITS: PREDICTIONS OF REACTION NORMS FOR AGE AND SIZE AT MATURITY

We used life-history theory to predict reaction norms for age and size at maturation. We assumed that fecundity increases with size and that juvenile mortality rates of offspring decrease as ages-at-maturity of parents increase, then calculated the reaction norm by varying growth rate and calculating an optimal age at maturity for each growth rate. The reaction norm for maturation should take one of at least four shapes that depend on specific relations between changes in growth rates and changes in adult mortality rates, juvenile mortality rates, or both. Most organisms should mature neither at a fixed size nor at a fixed age, but along an age-size trajectory. The model makes possible a clear distinction between the genetic and phenotypic components of variation. The evolved response to selection is reflected in the shape and position of the reaction norm. The phenotypic response of a single organism to rapid or slow growth is defined by the location of its maturation event as a point on the reaction norm. A quantitative test with data from 19 populations and species of fish showed that predictions were in good agreement with observations (r = 0.93, P < 0.0001). The predictions of the model also agreed qualitatively with observed phenotypic variation in age and size at maturity in humans, platyfish, fruit flies, and red deer. This preliminary success suggests that experiments designed to test the predictions directly will be worthwhile.     

Gene regulation,quantitative genetics and the evolution of reaction norms

Summary The ideas of phenotypic plasticity and of reaction norm are gaining prominence as important components of theories of phenotypic evolution. Our understanding of the role of phenotypic plasticity as an adaptation of organisms to variable environments will depend on (1) the form(s) of genetic and developmental control exerted on the shape of the reaction norm and (2) the nature of the constraints on the possible evolutionary trajectories in multiple environments. In this paper we identify two categories of genetic control of plasticity: allelic sensitivity and gene regulation. These correspond generally to two classes of response by the developmental system to environmental change: phenotypic modulation, in which plastic responses are a continuous and proportional function of environmental stimuli and developmental conversion, where responses tend to be not simply proportional to the stimuli. We propose that control of plasticity by regulatory actions has distinct advantages over simple allelic sensitivity: stability of phenotypic expression, capacity for anticipatory response and relaxation of constraints due to genetic correlations. We cite examples of the extensive molecular evidence for the existence of environmentally-cued gene regulation leading to developmental conversion. The results of quantitative genetic investigations on the genetics and evolution of plasticity, as well as the limits of current approaches are discussed. We suggest that evolution of reaction norms would be affected by the ecological context (i.e. spatial versus temporal variation, hard versus soft selection, and fine versus coarse environmental grain). We conclude by discussing some empirical approaches to address fundamental questions about plasticity evolution.     

Experimental assessment of the probabilistic maturation reaction norm: condition matters

The probabilistic maturation reaction norm (PMRN) describes an individual''s probability of maturing at a given age as a function of size and other relevant phenotypic traits. Population-level shifts in the PMRN are often interpreted to indicate genetic as opposed to phenotypic changes in maturation in fish. Inferences derived from trends in the PMRN have been challenged, warranting an experimental assessment of the method. This was accomplished in a laboratory experiment using zebrafish (Danio rerio). Fish were reared under different food levels to induce variation in growth and maturation. Plasticity in maturation was not entirely captured by the demographic age- and length-based PMRN. Adding condition to the PMRN captured a greater amount of environmental variation in maturation probability. Nevertheless, significant differences in the PMRNs among the food levels remained after accounting for the influences of age, size and condition on maturation probability indicating plasticity of the PMRN. This was particularly pronounced between fish held on low food levels as compared with fish experiencing abundant resources, with the latter experiencing higher size-specific maturation probabilities. Our analysis emphasizes the need for incorporating salient physiological traits influencing maturation, such as condition, to make accurate inferences about documented shifts observed in the position of PMRNs on maturation trends in wild fish stocks.     

Genetic variation in threshold reaction norms for alternative reproductive tactics in male Atlantic salmon, Salmo salar

Alternative reproductive tactics may be a product of adaptive phenotypic plasticity, such that discontinuous variation in life history depends on both the genotype and the environment. Phenotypes that fall below a genetically determined threshold adopt one tactic, while those exceeding the threshold adopt the alternative tactic. We report evidence of genetic variability in maturation thresholds for male Atlantic salmon (Salmo salar) that mature either as large (more than 1 kg) anadromous males or as small (10-150 g) parr. Using a common-garden experimental protocol, we find that the growth rate at which the sneaker parr phenotype is expressed differs among pure- and mixed-population crosses. Maturation thresholds of hybrids were intermediate to those of pure crosses, consistent with the hypothesis that the life-history switch points are heritable. Our work provides evidence, for a vertebrate, that thresholds for alternative reproductive tactics differ genetically among populations and can be modelled as discontinuous reaction norms for age and size at maturity.     

An empirical test for a zone of canalization in thermal reaction norms

Theoretical models on the evolution of phenotypic plasticity predict a zone of canalization where reaction norms cross, and genetic variation is minimized in the environment a population most frequently encounter. Empirical tests of this prediction are largely missing, in particular for life‐history traits. We addressed this prediction by quantifying thermal reaction norms of three life‐history traits (somatic growth rate, age and size at maturation) of a Norwegian population of Daphnia magna and testing for the occurrence of an intermediate temperature (T_m) at which genetic variance in the traits is minimized. Size at maturation changed relatively little with temperature compared to the other traits, and there was no genetic variance in the shape of the reaction norm. Consequently, age at maturation and somatic growth rate were strongly negatively correlated. Both traits showed a strong genotype–environment interaction, and the estimated T_m was 14 °C for both age at maturation and growth rate. This value of T_m corresponds well with mean summer temperatures experienced by the population and suggests that the population has evolved under stabilizing selection in temperatures that fluctuate around this mean temperature. These results suggest local adaptation to temperature in the studied population and allow predicting evolutionary trajectories of thermal reaction norms under changing thermal regimes.     

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.     

The evolution of human phenotypic plasticity: age and nutritional status at maturity

Several evolutionary optimal models of human plasticity in age and nutritional status at reproductive maturation are proposed and their dynamics examined. These models differ from previously published models because fertility is not assumed to be a function of body size or nutritional status. Further, the models are based on explicitly human demographic patterns, that is, model human life-tables, model human fertility tables, and, a nutrient flow-based model of maternal nutritional status. Infant survival (instead of fertility as in previous models) is assumed to be a function of maternal nutritional status. Two basic models are examined. In the first the cost of reproduction is assumed to be a constant proportion of total nutrient flow. In the second the cost of reproduction is constant for each birth. The constant proportion model predicts a negative slope of age and nutritional status at maturation. The constant cost per birth model predicts a positive slope of age and nutritional status at maturation. Either model can account for the secular decline in menarche observed over the last several centuries in Europe. A search of the growth literature failed to find definitive empirical documentation of human phenotypic plasticity in age and nutritional status at maturation. Most research strategies confound genetics with phenotypic plasticity. The one study that reports secular trends suggests a marginally insignificant, but positive slope. This view tends to support the constant cost per birth model.     

Phenotypic plasticity and constancy of life-history traits in laboratory clones of Daphnia magna Straus: effects of neonatal length

1. The phenotypic constancy of four laboratory Daphnia magna clones in fitness-related life-history traits, such as age and clutch size at maturity, was studied among consecutive experimental runs in differing food environments.
2. A significant part of the observed clonal and genetic-by-environmental variation in age and clutch size at maturity was explained by experimentally uncontrollable variations in neonatal body length.
3. Despite food availability, neonatal length determined the number of instars invested to maturity and thus maturation age. Clonal differences in neonatal length and thus in maturation instar occurrence across environments explained most of the clonal variability observed in maturation age.
4. Although interclonal differences in clutch size existed, most of the phenotypic plasticity observed for clutch size was mediated by clonal differences in neonatal length.
5. Clonal differences in neonatal length and in the occurrence of maturation instars across environments dramatically affected the body length of instar IM-2 where provisioning of eggs take place. Since clutch size is determined from clutch mass and clutch mass was strongly related to the body length of instar IM-2, clonal differences across environments in body length of instar IM-2 mirrored clonal differences across environments in clutch size.
6. The results reported in the present study show that maternally mediated traits such as neonatal length affect how genotypes respond to different environmental selection regimes (genetic-by-environmental interaction). Future research needs to focus on the effects of neonatal length on the heritability or genetic variation of the reaction norms, since prediction of the response to selection is a key research objective in quantitative genetic studies.     

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.     

HOW TO MEASURE MATURATION: A COMPARISON OF PROBABILISTIC METHODS USED TO TEST FOR GENOTYPIC VARIATION AND PLASTICITY IN THE DECISION TO MATURE

Maturation is a developmental trait that plays a key role in shaping organisms’ life‐history. However, progress in understanding how maturation phenotypes evolve has been held back by confusion over how best to model maturation decisions and a lack of studies comparing genotypic variation in maturation. Here, we fitted probabilistic maturation reaction norms (PMRNs) to data collected from five clones of Daphnia magna and five of Daphnia pulex collected from within and between different populations. We directly compared the utility of modeling approaches that assume maturation to be a process with an instantaneous rate with those that do not by fitting maturation rate and logistic regression models, and emphasize similarities and differences between them. Our results demonstrate that in Daphnia, PMRNs using a logistic regression approach were simpler to use and provided a better fit to the data. The decision to mature was plastic across a range of growth trajectories and dependent upon both body size and age. However, the age effect was stronger in D. magna than D. pulex and varied considerably between clones. Our results support the idea that maturation thresholds can evolve but also suggest that the notion of a threshold based on a single fixed state is an oversimplification that underestimates the adaptability of these important traits.     

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.