Measuring probabilistic reaction norms for age and size at maturation

We present a new probabilistic concept of reaction norms for age and size at maturation that is applicable when observations are carried out at discrete time intervals. This approach can also be used to estimate reaction norms for age and size at metamorphosis or at other ontogenetic transitions. Such estimations are critical for understanding phenotypic plasticity and life-history changes in variable environments, assessing genetic changes in the presence of phenotypic plasticity, and calibrating size- and age-structured population models. We show that previous approaches to this problem, based on regressing size against age at maturation, give results that are systematically biased when compared to the probabilistic reaction norms. The bias can be substantial and is likely to lead to qualitatively incorrect conclusions; it is caused by failing to account for the probabilistic nature of the maturation process. We explain why, instead, robust estimations of maturation reaction norms should be based on logistic regression or on other statistical models that treat the probability of maturing as a dependent variable. We demonstrate the utility of our approach with two examples. First, the analysis of data generated for a known reaction norm highlights some crucial limitations of previous approaches. Second, application to the northeast arctic cod (Gadus morhua) illustrates how our approach can be used to shed new light on existing real-world data.     

The evolution of reaction norms: simple models for age and size at maturity

This paper presents a simple model for the evolution of reaction norms for age and size at maturity that predicts reaction norms with a variety of shapes. Using realistic parameter values the model predicts reaction norms close to those observed in Drosophila. The major assumptions of the model are: 1) that net reproductive rate is maximized, 2) that growth is determinate, and 3) that mortality rates are independent of age and size at maturity. If, additionally, juvenile mortality is uncorrelated with a growth coefficient, k, the model predicts that selection favors maturation later at a smaller size when k is reduced by environmental factors and that decreased juvenile mortality leads to delayed maturity. These two predictions conform with those found by previous models using other measures of fitness. Correlations between k and juvenile mortality can change the shape of the predicted reaction norm. Depending on the precise form of the correlation, the model can predict done- or bowl-shaped reaction norms and can predict delayed or earlier maturity as k decreases. These shapes are qualitatively different from those predicted by previous models that used different fitness measures. Systematic estimates of the parameter values for this and for related models are required to determine the appropriate fitness measure for models of reaction norms.     

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.     

Exploring patterns of variation in clutch size–density reaction norms in a wild passerine bird

Negative density dependence of clutch size is a ubiquitous characteristic of avian populations and is partly due to within‐individual phenotypic plasticity. Yet, very little is known about the extent to which individuals differ in their degree of phenotypic plasticity, whether such variation has a genetic basis and whether level of plasticity can thus evolve in response to selection. Using 18 years of data of a Dutch great tit population (Parus major), we show that females reduced clutch size with increasing population density (slopes of the reaction norms), differed strongly in their average clutch size (elevations of the reaction norms) at the population‐mean density and that the latter variation was partly heritable. In contrast, we could not detect individual variation in phenotypic plasticity (‘I × E’). Level of plasticity is thus not likely to evolve in response to selection in this population. Observed clutch sizes deviated more from the estimated individual reaction norms in certain years and densities, implying that the within‐individual between‐year variance (so‐called residual variance) of clutch size was heterogeneous with respect to these factors. Given the observational nature of this study, experimental manipulation of density is now warranted to confirm the causality of the observed density effects. Our analyses demonstrate that failure to acknowledge this heterogeneity would have inflated the estimate of ‘I × E’ and led to misinterpretation of the data. This paper thereby emphasizes the fact that heterogeneity in residuals can provide biologically insightful information about the ecological processes underlying the data.     

Local adaptation of a Drosophila parasitoid: habitat‐specific differences in thermal reaction norms

Local climate is an important source of selection on thermal reaction norms that has been well investigated in cline studies, where populations sampled along altitudinal or latitudinal gradients are compared. Several biotic factors vary with climate, but are rarely integrated as alternative agents of selection to climatic factors. We tested the hypothesis that habitat may select for thermal reaction norms and magnitude of phenotypic plasticity in a drosophila parasitoid, independently of the climate of origin. We sampled populations of Leptopilina boulardi, a Drosophila parasitoid in two different habitats, orchards and forests. Orchards offer laying opportunities over small distances for parasitoids, with a low variability in the number of hosts per patch, while forests offer more dispersed and more variable patches. The sampling was realized in a temperate and a Mediterranean climate. We measured egg load, volume of eggs, longevity and lipid content for parasitoids reared at two temperatures. Reaction norms were opposite for populations from forests and orchards for investment in reproduction, independently of the climate of origin. The maximal investment of resources in reproduction occurred at the lower temperature in orchards and the higher temperature in forests. Host distribution differences between habitats may explain these opposite reaction norms. We also observed a flatter reaction norm for egg load in forests than in orchards. This relative canalization may have been selected in response to the higher variability in laying opportunities observed in forests. Our results demonstrate the potential role of resource distribution in evolution of thermal plasticity.     

Squeezing out the last egg—annual fish increase reproductive efforts in response to a predation threat

Both constitutive and inducible antipredator strategies are ubiquitous in nature and serve to maximize fitness under a predation threat. Inducible strategies may be favored over constitutive defenses depending on their relative cost and benefit and temporal variability in predator presence. In African temporary ponds, annual killifish of the genus Nothobranchius are variably exposed to predators, depending on whether larger fish invade their habitat from nearby rivers during floods. Nonetheless, potential plastic responses to predation risk are poorly known. Here, we studied whether Nothobranchius furzeri individuals adjust their life history in response to a predation threat. For this, we monitored key life history traits in response to cues that signal the presence of predatory pumpkinseed sunfish (Lepomis gibbosus). While growth rate, adult body size, age at maturation, and initial fecundity were not affected, peak and total fecundity were higher in the predation risk treatment. This contrasts with known life history strategies of killifish from permanent waters, which tend to reduce reproduction in the presence of predators. Although our results show that N. furzeri individuals are able to detect predators and respond to their presence by modulating their reproductive output, these responses only become evident after a few clutches have been deposited. Overall our findings suggest that, in the presence of a predation risk, it can be beneficial to increase the production of life stages that can persist until the predation risk has faded.     

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.     

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.     

Adaptive divergence in body size overrides the effects of plasticity across natural habitats in the brown trout

The evolution of life-history traits is characterized by trade-offs between different selection pressures, as well as plasticity across environmental conditions. Yet, studies on local adaptation are often performed under artificial conditions, leaving two issues unexplored: (i) how consistent are laboratory inferred local adaptations under natural conditions and (ii) how much phenotypic variation is attributed to phenotypic plasticity and to adaptive evolution, respectively, across environmental conditions? We reared fish from six locally adapted (domesticated and wild) populations of anadromous brown trout (Salmo trutta) in one semi-natural and three natural streams and recorded a key life-history trait (body size at the end of first growth season). We found that population-specific reaction norms were close to parallel across different streams and Q_ST was similar – and larger than F_ST – within all streams, indicating a consistency of local adaptation in body size across natural environments. The amount of variation explained by population origin exceeded the variation across stream environments, indicating that genetic effects derived from adaptive processes have a stronger effect on phenotypic variation than plasticity induced by environmental conditions. These results suggest that plasticity does not “swamp” the phenotypic variation, and that selection may thus be efficient in generating genetic change.     

Phenotypic convergence of artificially reared and wild trout is mediated by shape plasticity

Phenotypic plasticity can be viewed as the first level of defense of organism homeostasis against environmental stress and therefore represents the potential to deal with rapid environmental changes. Transitions between low complexity, artificial environments and complex, natural habitats can promote phenotypic plasticity. Here, we conducted an experimental introduction with juvenile brown trout to evaluate the plasticity of shape in response to a transition between contrasting environments. We released 202 juvenile trout reared under hatchery conditions in a natural stream and analyzed changes in shape and morphological variability after 5 months. A geometric morphometrics approach based on 14 landmarks was used to compare changes in body shape for 37 fish recaptured at the end of the experiment. A similar number of hatchery and wild fish caught at the receptor stream were used as controls for shape in the two environments. After 5‐months, fish showed significant change in shape, shifting from elongated to robust shapes, and affecting to the relative position of the caudal peduncle. These new shapes were closer to wild than to the hatchery shapes, suggesting a process of rapid phenotype change. Moreover, these changes were concomitant with a marked increase in morphological variability. Our results support the hypothesis that phenotypic plasticity is a major potential for adjustment to environmental change but not the idea that shape can be constrained by initial shapes. We confirmed the “increased” variance hypothesis and phenotype convergence with wild morphs. This has important implications because stresses the role of phenotypic plasticity as a buffer that allows organisms to cope with important environmental discontinuities at time scales that preclude the onset of adaptive adjustments. We suggest that environmental conditioning and shape plasticity can overcome both reduced morphological diversity and phenotype uncoupling with habitat characteristics resulting from initial rearing in low complexity artificial environments.     

The evolution of life histories in spatially heterogeneous environments: Optimal reaction norms revisited

Summary Natural populations live in heterogeneous environments, where habitat variation drives the evolution of phenotypic plasticity. The key feature of population structure addressed in this paper is the net flow of individuals from source (good) to sink (poor) habitats. These movements make it necessary to calculate fitness across the full range of habitats encountered by the population, rather than independently for each habitat. As a consequence, the optimal phenotype in a given habitat not only depends on conditions there but is linked to the performance of individuals in other habitats. We generalize the Euler-Lotka equation to define fitness in a spatially heterogeneous environment in which individuals disperse among habitats as newborn and then stay in a given habitat for life. In this case, maximizing fitness (the rate of increase over all habitats) is equivalent to maximizing the reproductive value of newborn in each habitat but not to maximizing the rate of increase that would result if individuals in each habitat were an isolated population. The new equation can be used to find optimal reaction norms for life history traits, and examples are calculated for age at maturity and clutch size. In contrast to previous results, the optimal reaction norm differs from the line connecting local adaptations of isolated populations each living in only one habitat. Selection pressure is higher in good and frequent habitats than in poor and rare ones. A formula for the relative importance of these two factors allows predictions of the habitat in which the genetic variance about the optimal reaction norm should be smallest.     

Genomic consequences of a recent three‐way admixture in supplemented wild brown trout populations revealed by local ancestry tracts

Understanding the evolutionary consequences of human‐mediated introductions of domesticated strains into the wild and their subsequent admixture with natural populations is of major concern in conservation biology. However, the genomic impacts of stocking from distinct sources (locally derived vs. divergent) on the genetic integrity of wild populations remain poorly understood. We designed an approach based on estimating local ancestry along individual chromosomes to provide a detailed picture of genomic admixture in supplemented populations. We used this approach to document admixture consequences in the brown trout Salmo trutta, for which decades of stocking practices have profoundly impacted the genetic make‐up of wild populations. In southern France, small local Mediterranean populations have been subject to successive introductions of domestic strains derived from the Atlantic and Mediterranean lineages. To address the impact of stocking, we evaluate the extent of admixture from both domestic strains within populations, using 75,684 mapped SNPs obtained from double‐digested restriction site‐associated DNA sequencing. Then, the chromosomal ancestry profiles of admixed individuals reveal a wider diversity of hybrid and introgressed genotypes than estimated using classical methods for inferring ancestry and hybrid pedigrees. In addition, the length distribution of introgressed tracts retained different timings of introgression between the two domestic strains. We finally reveal opposite consequences of admixture on the level of polymorphism of the recipient populations between domestic strains. Our study illustrates the potential of using the information contained in the genomic mosaic of ancestry tracts in combination with classical methods based on allele frequencies for analysing multiple‐way admixture with population genomic data.     

Early exposure to nonlethal predation risk by size-selective predators increases somatic growth and decreases size at adulthood in three-spined sticklebacks

Predation has an important influence on life history traits in many organisms, especially when they are young. When cues of trout were present, juvenile sticklebacks grew faster. The increase in body size as a result of exposure to cues of predators was adaptive because larger individuals were more likely to survive predation. However, sticklebacks that had been exposed to cues of predators were smaller at adulthood. This result is consistent with some life history theory. However, these results prompt an alternative hypothesis, which is that the decreased size at adulthood reflects a deferred cost of early rapid growth. Compared to males, females were more likely to survive predation, but female size at adulthood was more affected by cues of predators than male size at adulthood, suggesting that size at adulthood might be more important to male fitness than to female fitness.     

Flowering time QTL in natural populations of Arabidopsis thaliana and implications for their adaptive value

The genetic basis of phenotypic traits is of great interest to evolutionary biologists, but their contribution to adaptation in nature is often unknown. To determine the genetic architecture of flowering time in ecologically relevant conditions, we used a recombinant inbred line population created from two locally adapted populations of Arabidopsis thaliana from Sweden and Italy. Using these RILs, we identified flowering time QTL in growth chambers that mimicked the natural temperature and photoperiod variation across the growing season in each native environment. We also compared the genomic locations of flowering time QTL to those of fitness (total fruit number) QTL from a previous three‐year field study. Ten total flowering time QTL were found, and in all cases, the Italy genotype caused early flowering regardless of the conditions. Two QTL were consistent across chamber environments, and these had the largest effects on flowering time. Five of the fitness QTL colocalized with flowering time QTL found in the Italy conditions, and in each case, the local genotype was favoured. In contrast, just two flowering time QTL found in the Sweden conditions colocalized with fitness QTL and in only one case was the local genotype favoured. This implies that flowering time may be more important for adaptation in Italy than Sweden. Two candidate genes (FLC and VIN3) underlying the major flowering time QTL found in the current study are implicated in local adaptation.     

Timing matters: species‐specific interactions between spawning time,substrate quality,and recruitment success in three salmonid species

Substratum quality and oxygen supply to the interstitial zone are crucial for the reproductive success of salmonid fishes. At present, degradation of spawning grounds due to fine sediment deposition and colmation are recognized as main factors for reproductive failure. In addition, changes in water temperatures due to climate change, damming, and cooling water inlets are predicted to reduce hatching success. We tested the hypothesis that the biological effects of habitat degradation depend strongly on the species‐specific spawning seasons and life‐history strategies (e.g., fall‐ vs. spring‐spawners, migratory vs. resident species) and assessed temperature as an important species‐specific factor for hatching success within river substratum. We studied the species‐specific differences in their responses to such disturbances using egg‐to‐fry survival of Danube Salmon (Hucho hucho), resident brown trout (Salmo trutta fario), and migratory brown trout (Salmo trutta lacustris) as biological endpoint. The egg incubation and hatching success of the salmonids and their dependence on temperature and stream substratum quality were compared. Hatching rates of Danube salmon were lower than of brown trout, probably due to higher oxygen demands and increased interstitial respiration in spring. Increases in maximum water temperature reduced hatching rates of resident and migratory brown trout (both fall‐spawners) but were positively correlated with hatching rates of Danube salmon (a spring‐spawner). Significantly longer incubation periods of resident and migratory brown trout coincided with relatively low stream substratum quality at the end of the egg incubation. Danube salmon seem to avoid low oxygen concentrations in the hyporheic zone by faster egg development favored by higher water temperatures. Consequently, the prediction of effects of temperature changes and altered stream substratum properties on gravel‐spawning fishes and biological communities should consider the observed species‐specific variances in life‐history strategies to increase conservation success.     

Environmental determinants of population divergence in life‐history traits for an invasive species: climate,seasonality and natural enemies

Invasive species cope with novel environments through both phenotypic plasticity and evolutionary change. However, the environmental factors that cause evolutionary divergence in invasive species are poorly understood. We developed predictions for how different life‐history traits, and plasticity in those traits, may respond to environmental gradients in seasonal temperatures, season length and natural enemies. We then tested these predictions in four geographic populations of the invasive cabbage white butterfly (Pieris rapae) from North America. We examined the influence of two rearing temperatures (20 and 26.7 °C) on pupal mass, pupal development time, immune function and fecundity. As predicted, development time was shorter and immune function was greater in populations adapted to longer season length. Also, phenotypic plasticity in development time was greater in regions with shorter growing seasons. Populations differed significantly in mean and plasticity of body mass and fecundity, but these differences were not associated with seasonal temperatures or season length. Our study shows that some life‐history traits, such as development time and immune function, can evolve rapidly in response to latitudinal variation in season length and natural enemies, whereas others traits did not. Our results also indicate that phenotypic plasticity in development time can also diverge rapidly in response to environmental conditions for some traits.     

Phenotypic integration plasticity in Daphnia magna: an integral facet of G × E interactions

Phenotypic integration can be defined as the network of multivariate relationships among behavioural, physiological and morphological traits that describe the organism. Phenotypic integration plasticity refers to the change in patterns of phenotypic integration across environments or ontogeny. Because studies of phenotypic plasticity have predominantly focussed on single traits, a G × E interaction is typically perceived as differences in the magnitude of trait expression across two or more environments. However, many plastic responses involve coordinated responses in multiple traits, raising the possibility that relative differences in trait expression in different environments are an important, but often overlooked, source of G × E interaction. Here, we use phenotypic change vectors to statistically compare the multivariate life‐history plasticity of six Daphnia magna clones collected from four disparate European populations. Differences in the magnitude of plastic responses were statistically distinguishable for two of the six clones studied. However, differences in phenotypic integration plasticity were statistically distinguishable for all six of the clones studied, suggesting that phenotypic integration plasticity is an important component of G × E interactions that may be missed unless appropriate multivariate analyses are used.     

Insulin‐like growth factor 1 relieves the constraints on the growth of young wild passerines

To gain a selective advantage for survival in stochastic environments, the growth of different body parameters of juvenile animals should be constantly adjusted according to prevailing conditions. Hormones, especially insulin‐like growth factor 1 (IGF‐1), are an important part of physiological mechanisms mediating life‐history variability in free‐living animals when connecting available resources (e.g. food) with pathways of somatic growth. We used an IGF‐1 injection treatment in free‐living European Pied Flycatcher Ficedula hypoleuca nestlings to mimic experimentally the differentiation of growth conditions for chicks with a similar genetic background. We showed that there is probably a physiological trade‐off for young animals between the growth rates of structural size and body mass, where IGF‐1 could be part of the physiological modulatory system of this trade‐off. By weakening internal constraints that limit growth, IGF‐1 could help to relieve the trade‐off between these competing body size parameters.     

Age and size at maturity: A quantitative review of diet‐induced reaction norms in insects

Optimality models predict that diet‐induced bivariate reaction norms for age and size at maturity can have diverse shapes, with the slope varying from negative to positive. To evaluate these predictions, we perform a quantitative review of relevant data, using a literature‐derived database of body sizes and development times for over 200 insect species. We show that bivariate reaction norms with a negative slope prevail in nearly all taxonomic and ecological categories of insects as well as in some other ectotherm taxa with comparable life histories (arachnids and amphibians). In insects, positive slopes are largely limited to species, which feed on discrete resource items, parasitoids in particular. By contrast, with virtually no meaningful exceptions, herbivorous and predatory insects display reaction norms with a negative slope. This is consistent with the idea that predictable resource depletion, a scenario selecting for positively sloped reaction norms, is not frequent for these insects. Another source of such selection—a positive correlation between resource levels and juvenile mortality rates—should similarly be rare among insects. Positive slopes can also be predicted by models which integrate life‐history evolution and population dynamics. As bottom‐up regulation is not common in most insect groups, such models may not be most appropriate for insects.     

On the fate of seasonally plastic traits in a rainforest butterfly under relaxed selection

Many organisms display phenotypic plasticity as adaptation to seasonal environmental fluctuations. Often, such seasonal responses entails plasticity of a whole suite of morphological and life‐history traits that together contribute to the adaptive phenotypes in the alternative environments. While phenotypic plasticity in general is a well‐studied phenomenon, little is known about the evolutionary fate of plastic responses if natural selection on plasticity is relaxed. Here, we study whether the presumed ancestral seasonal plasticity of the rainforest butterfly Bicyclus sanaos (Fabricius, 1793) is still retained despite the fact that this species inhabits an environmentally stable habitat. Being exposed to an atypical range of temperatures in the laboratory revealed hidden reaction norms for several traits, including wing pattern. In contrast, reproductive body allocation has lost the plastic response. In the savannah butterfly, B. anynana (Butler, 1879), these traits show strong developmental plasticity as an adaptation to the contrasting environments of its seasonal habitat and they are coordinated via a common developmental hormonal system. Our results for B. sanaos indicate that such integration of plastic traits – as a result of past selection on expressing a coordinated environmental response – can be broken when the optimal reaction norms for those traits diverge in a new environment.