Adaptive divergence in plasticity in natural populations of Impatiens capensis and its consequences for performance in novel habitats

We tested for adaptive differentiation between two natural populations of Impatiens capensis from sites known to differ in selection on plasticity to density. We also determined the degree to which plasticity to density within a site was correlated with plastic responses of experimental immigrants to foreign sites. Inbred lines, derived from natural populations in an open-canopy site and a woodland site, were planted reciprocally in both original sites at naturally occurring high densities and at low density. The density manipulation represents environmental variation typically experienced within the site of a given population, and the transplant manipulation represents environmental differences between sites of different populations. Internode elongation, meristem allocation, leaf length, flowering date, and total lifetime fitness were measured. Genotypes originating in the open site, where selection favored plasticity of first internode length and flowering time (Donohue et al. 2000a), were more plastic in those characters than genotypes originating from the woodland site, where plasticity was maladaptive. Therefore, these two populations appear to have responded to divergent selection on plasticity. Plasticity to density strongly resembled plasticity to site differences for many characters, suggesting that similar environmental factors elicit plasticity both to density and to overhead canopy. Thus, plasticity that evolved in response to density variation within a site influenced phenotypic expression in the foreign site. Plastic responses to site caused immigrants from foreign populations to resemble native genotypes more closely. In particular, immigrants from the open site converged toward the selectively favored early-flowering phenotype of native genotypes in the woodland site, thereby reducing potential fitness differences between foreign and native genotypes. However, because genotypes from the woods population were less plastic than genotypes from the sun population, phenotypic differences between populations were greatest in the open site at low density. Therefore, population differences in plasticity can cause genotypes from foreign populations to be more strongly selected against in some environments than in others. However, genetic constraints and limits to plasticity prevented complete convergence of immigrants to the native phenotype in any environment.     

Evidence of adaptive divergence in plasticity: density- and site-dependent selection on shade-avoidance responses in Impatiens capensis

We investigated the conditions under which plastic responses to density are adaptive in natural populations of Impatiens capensis and determined whether plasticity has evolved differently in different selective environments. Previous studies showed that a population that evolved in a sunny site exhibited greater plasticity in response to density than did a population that evolved in a woodland site. Using replicate inbred lines in a reciprocal transplant that included a density manipulation, we asked whether such population differentiation was consistent with the hypothesis of adaptive divergence. We hypothesized that plasticity would be more strongly favored in the sunny site than in the woodland site; consequently, we predicted that selection would be more strongly density dependent in the sunny site, favoring the phenotype that was expressed at each density. Selection on internode length and flowering date was consistent with the hypothesis of adaptive divergence in plasticity. Few costs or benefits of plasticity were detected independently from the expressed phenotype, so plasticity was selected primarily through selection on the phenotype. Correlations between phenotypes and their plasticity varied with the environment and would cause indirect selection on plasticity to be environment dependent. We showed that an appropriate plastic response even to a rare environment can greatly increase genotypic fitness when that environment is favorable. Selection on the measured characters contributed to local adaptation and fully accounted for fitness differences between populations in all treatments except the woodland site at natural density.     

HOMOGENEITY OF THE GENETIC VARIANCE-COVARIANCE MATRIX FOR ANTIPREDATOR TRAITS IN TWO NATURAL POPULATIONS OF THE GARTER SNAKE THAMNOPHIS ORDINOIDES

Quantitative genetic models of evolution rely on the genetic variance-covariance matrix to predict the phenotypic response to selection. Both prospective and retrospective studies of phenotypic evolution across generations rely on assumptions about the constancy of patterns of genetic covariance through time. In the absence of robust theoretical predictions about the stability of genetic covariances, this assumption must be tested with empirical comparisons of genetic parameters among populations and species. Genetic variance-covariance matrices were estimated for a suite of antipredator traits in two populations of the northwestern garter snake, Thamnophis ordinoides. The characters studied include color pattern and antipredator behaviors that interact to facilitate escape from predators. Significant heritabilities for all traits were detected in both populations. Genetic correlations and covariances were found among behaviors in both populations and between color pattern and behavior in one of the populations. Phenotypic means differed among populations, but pairwise comparisons revealed no heterogeneity of genetic parameters between the populations. The structure of the genetic variance-covariance matrix has apparently not changed significantly during the divergence of these two populations.     

The evolution of sexual dimorphism in the house finch. I. Population divergence in morphological covariance structure

Abstract Patterns of genetic variation and covariation strongly affect the rate and direction of evolutionary change by limiting the amount and form of genetic variation available to natural selection. We studied evolution of morphological variance-covariance structure among seven populations of house finches (Carpodacus mexicanus) with a known phylogenetic history. We examined the relationship between within- and among-population covariance structure and, in particular, tested the concordance between hierarchical changes in morphological variance-covariance structure and phylogenetic history of this species. We found that among-population morphological divergence in either males or females did not follow the within-population covariance patterns. Hierarchical patterns of similarity in morphological covariance matrices were not congruent with a priori defined historical pattern of population divergence. Both of these results point to the lack of proportionality in morphological covariance structure of finch populations, suggesting that random drift alone is unlikely to account for observed divergence. Furthermore, drift alone cannot explain the sex differences in within- and among-population covariance patterns or sex-specific patterns of evolution of covariance structure. Our results suggest that extensive among-population variation in sexual dimorphism in morphological covariance structure was produced by population differences in local selection pressures acting on each sex.     

The quantitative genetics of floral trait variation in Lobelia: potential constraints on adaptive evolution

Abstract Although pollinator-mediated natural selection has been measured on many floral traits and in many species, the extent to which selection is constrained from producing optimal floral phenotypes is less frequently studied. In particular, negative correlations between flower size and flower number are hypothesized to be a major constraint on the evolution of floral displays, yet few empirical studies have documented such a trade-off. To determine the potential for genetic constraints on the adaptive evolution of floral displays, I estimated the quantitative genetic basis of floral trait variation in two populations of Lobelia siphilitica . Restricted maximum likelihood (REML) analyses of greenhouse-grown half-sib families were used to estimate genetic variances and covariances for flower number and six measures of flower size. There was significant genetic variation for all seven floral traits in both populations. Flower number was negatively genetically correlated with four measures of flower size in one population and three measures in the other. When the genetic variance-covariance matrices were combined with field estimates of phenotypic selection gradients, the predicted multivariate evolutionary response was less than or opposite in sign to the selection gradient for flower number and five of six measures of flower size, suggesting genetic constraints on the evolution of these traits. More generally, my results indicate that the adaptive evolution of floral displays can be constrained by tradeoffs between flower size and number, as has been assumed by many theoretical models of floral evolution.     

Constraints on the evolution of adaptive phenotypic plasticity in plants

The high potential fitness benefit of phenotypic plasticity tempts us to expect phenotypic plasticity as a frequent adaptation to environmental heterogeneity. Examples of proven adaptive plasticity in plants, however, are scarce and most plastic responses actually may be 'passive' rather than adaptive. This suggests that frequently requirements for the evolution of adaptive plasticity are not met or that such evolution is impeded by constraints. Here we outline requirements and potential constraints for the evolution of adaptive phenotypic plasticity, identify open questions, and propose new research approaches. Important open questions concern the genetic background of plasticity, genetic variation in plasticity, selection for plasticity in natural habitats, and the nature and occurrence of costs and limits of plasticity. Especially promising tools to address these questions are selection gradient analysis, meta-analysis of studies on genotype-by-environment interactions, QTL analysis, cDNA-microarray scanning and quantitative PCR to quantify gene expression, and two-dimensional gel electrophoresis to quantify protein expression. Studying plasticity along the pathway from gene expression to the phenotype and its relationship with fitness will help us to better understand why adaptive plasticity is not more universal, and to more realistically predict the evolution of plastic responses to environmental change.     

THE GENETIC ARCHITECTURE OF PLASTICITY TO DENSITY IN IMPATIENS CAPENSIS

Plant responses to crowding may be mediated by resource availability and/or by a specific environmental cue, the ratio of red:far red wavelengths (R:FR) perceived by phytochrome. This study examined the contribution of phytochrome-mediated photomorphogenesis to genetic variation in plastic responses to density in the annual plant Impatiens capensis. Inbred lines derived from open and woodland populations were grown under low density high density, and high density with selective removal of FR wavelengths to block phytochrome-mediated perception of neighbor proximity. Genetic variation in plasticity to density and to the R:FR cue was detected for several traits Plants grown at high density displayed increased internode elongation; decreased branch, flower, and node production; increased menstem dormancy; and decreased leaf area and specific leaf weight compared to plants grown at low density. Stem elongation responses to density were suppressed when phytochrome perception was blocked at high density. For these phytochrome-mediated traits, a genotype's plasticity to density was strongly correlated with its response to R:FR. Phytochrome-mediated traits were tightly correlated with one another, regardless of the density environment. However, the responses to density of meristem allocation to branching and leaf traits were less strongly phytochrome-mediated. These traits differed in patterns of plasticity, and their genetic correlations often differed across environments. In particular, genetic trade-offs involving meristem allocation to branching were expressed only at low density. The observed density dependence of phenotypic and genetic correlations implies that indirect selection and the potential for correlated response to selection will depend upon the competitive environment. Thus the differential sensitivity of characters to the R:FR cue can influence the evolution of integrated plastic responses to density.     

Genetic constraints on the evolution of mate recognition under natural selection

Field populations of Drosophila serrata display reproductive character displacement in cuticular hydrocarbons (CHCs) when sympatric with Drosophila birchii. We have previously shown that the naturally occurring pattern of reproductive character displacement can be experimentally replicated by exposing field allopatric populations of D. serrata to experimental sympatry with D. birchii. Here, we tested whether the repeated evolution of reproductive character displacement in natural and experimental populations was a consequence of genetic constraints on the evolution of CHCs. The genetic variance-covariance (G) matrices for CHCs were determined for populations of D. serrata that had evolved in either the presence or absence of D. birchii under field and experimental conditions. Natural selection on mate recognition under both field and experimental sympatric conditions increased the genetic variance in CHCs consistent with a response to selection based on rare alleles. A close association between G eigenstructure and the eigenstructure of the phenotypic divergence (D) matrix in natural and experimental populations suggested that G matrix eigenstructure may have determined the direction in which reproductive character displacement evolved during the reinforcement of mate recognition.     

MUTATIONAL CONTRIBUTIONS TO GENETIC VARIANCE-COVARIANCE MATRICES: AN EXPERIMENTAL APPROACH USING INDUCED MUTATIONS IN ARABIDOPSIS THALIANA

Genetic potential for evolutionary change and covariational constraints are typically summarized as the genetic variance-covariance matrix G , and there is currently debate over the extent to which G remains effectively constant during the course of adaptive evolution. However, G provides only a temporally restricted view of constraints that ignores possible biases in how new mutations affect multivariate phenotypes. We used chemical mutagenesis to study the effect of mutations as summarized by the mutational covariance matrix, M , in Arabidopsis thaliana. By introducing mutations into three isogenic strains of A. thaliana, we were able to quantify M directly as the genetic variance-covariance matrix of mutagenized lines. Induced mutations generally did not alter the means of the six morphology and life-history traits we measured, but they did affect the levels of available genetic variation and the covariances among traits. However, these effects were not consistent among the three isogenic lines; that is, there were significant differences among the lines in both the number of mutations produced by ethyl-methane-sulfonate treatment and the M matrices they induced. The evolutionary implications of the dependence of M on the number of mutations, the particular genetic background, and the mutagenic sampling of loci in the genome are discussed in light of commonly applied models of multivariate evolution and the potential for the genetic architecture itself to change in ways that facilitate the coordinated evolution of complex phenotypes.     

Population divergence along lines of genetic variance and covariance in the invasive plant Lythrum salicaria in eastern North America

Evolution during biological invasion may occur over contemporary timescales, but the rate of evolutionary change may be inhibited by a lack of standing genetic variation for ecologically relevant traits and by fitness trade-offs among them. The extent to which these genetic constraints limit the evolution of local adaptation during biological invasion has rarely been examined. To investigate genetic constraints on life-history traits, we measured standing genetic variance and covariance in 20 populations of the invasive plant purple loosestrife (Lythrum salicaria) sampled along a latitudinal climatic gradient in eastern North America and grown under uniform conditions in a glasshouse. Genetic variances within and among populations were significant for all traits; however, strong intercorrelations among measurements of seedling growth rate, time to reproductive maturity and adult size suggested that fitness trade-offs have constrained population divergence. Evidence to support this hypothesis was obtained from the genetic variance-covariance matrix (G) and the matrix of (co)variance among population means (D), which were 79.8% (95% C.I. 77.7-82.9%) similar. These results suggest that population divergence during invasive spread of L. salicaria in eastern North America has been constrained by strong genetic correlations among life-history traits, despite large amounts of standing genetic variation for individual traits.     

Sexually dimorphic inflorescence traits in a wind-pollinated species: heritabilities and genetic correlations in Schiedea adamantis (Caryophyllaceae)

Sexual dimorphism may be especially pronounced in wind-pollinated species because they lack the constraints of biotically pollinated species that must present their pollen and stigmas in similar positions to ensure pollen transfer. Lacking these constraints, the sexes of wind-pollinated species may diverge in response to the different demands of pollen dispersal and receipt, depending on the magnitude of genetic correlations preventing divergence between sexes. Patterns of sexual dimorphism and genetic variation were investigated for inflorescence traits in Schiedea adamantis (Caryophyllaceae), a species well adapted to wind-pollination, and compared to S. salicaria, a species with fewer adaptations to wind pollination. For S. adamantis, dimorphism was pronounced for inflorescence condensation and its components, including lateral flower number and pedicel length. Within sexes, genetic correlations between traits may constrain the relative shape of the inflorescence. Correlations detected across sexes may retard the evolution of sexual dimorphism in inflorescence structure, including features favoring enhanced dispersal and receipt of pollen. Despite genetic correlations across sexes, common principal components analysis showed that genetic variance-covariance matrices (G matrices) differed significantly between the sexes, in part because of greater genetic variation for flower number in hermaphrodites than in females. G matrices also differed between closely related S. adamantis and S. salicaria, indicating the potential for divergent evolution of inflorescence structure despite general similarities in morphology and pollination biology.     

Plasticity versus canalization: population differences in the timing of shade-avoidance responses

The reliability of environmental cues and costs of a fixed phenotype are two factors determining whether selection favors phenotypic plasticity or environmental specialization. This study examines the relationship between these two factors and the evolution of plant competitive strategies (plastic vs. fixed morphologies). In natural plant populations, shifts in light quality associated with foliar shade reliably indicate the presence of neighbors. These cues mediate plastic stem-elongation responses that often increase competitive ability and access to light. Using experimental light treatments (full sun, neutral shade, and foliar shade), genetic differences among populations of Abutilon theophrasti (velvetleaf) in average elongation and plasticity to foliar-shade cues were examined. Six populations, two from each of three site types (fields in continuous corn cultivation, fields undergoing corn-soy rotation, and weedy sites), were exposed to the light treatments at two stages in their life history. At the seedling stage, populations derived from cornfield sites exhibited higher, average elongation than populations from either rotating corn-soy fields or weedy areas. Because seedling elongation may delay shading of velvetleaf by corn, population differences may reflect adaptive responses to directional selection imposed by competitive conditions. However, the effects of simulated foliar shade on elongation were three times as great as the average population differences, and these comparatively higher levels of elongation were associated with an allocation cost. These results are consistent with the hypothesis that phenotypic plasticity may limit the evolution of specialists; reliable environmental cues enable individuals to facultatively adopt highly elongated, costly phenotypes in crowded patches while avoiding the costs of that phenotype in less crowded microsites. At later life-history stages, populations experiencing competition with corn exhibited lower plasticity to light quality than populations derived from weedy areas. Elongation at later nodes is maladaptive in cornfields because velvetleaf is ultimately incapable of overtopping corn; individuals that elongate therefore experience the cost of allocating to stems but fail to improve leaf exposure. The decreased responsiveness of cornfield populations to light quality is consistent with theoretical predictions in which reduced plasticity is favored when environmental cues fail to mediate an adaptive response.     

Variability in plant fitness influences range expansion of Leptospermum scoparium

A reciprocal transplant experiment was conducted to uncouple genetic and environmental effects on plant fitness of a range-expanding indigenous shrub, Leptospermum scoparium . Patterns of encroachment differed significantly between three regions of Eucalyptus camaldulensis woodland in the Grampians National Park, south-eastern Australia, with greater shrub densities found in the Victoria Valley region. Higher shrub densities in this region are likely due to greater fruit production/shrub, increased field germination and increased early seedling survival. The results of the reciprocal transplant experiment showed increased fitness appears to be due to phenotypic plasticity in response to favourable environmental conditions rather than genotypic variation. Field germination and seedling survival were influenced by region of planting and not seed origin. However, differences in shoot length and biomass of seedlings between the regions were detected after five months of growth in uniform conditions which may indicate genetic differentiation. The results, however, suggest range expansion is partly due to increased plant fitness which is likely a plastic response to favourable environmental conditions, especially soil moisture.     

THE EVOLUTIONARY STABILITY OF CROSS‐SEX,CROSS‐TRAIT GENETIC COVARIANCES

Although knowledge of the selective agents behind the evolution of sexual dimorphism has advanced considerably in recent years, we still lack a clear understanding of the evolutionary durability of cross‐sex genetic covariances that often constrain its evolution. We tested the relative stability of cross‐sex genetic covariances for a suite of homologous contact pheromones of the fruit fly Drosophila serrata, along a latitudinal gradient where these traits have diverged in mean. Using a Bayesian framework, which allowed us to account for uncertainty in all parameter estimates, we compared divergence in the total amount and orientation of genetic variance across populations, finding divergence in orientation but not total variance. We then statistically compared orientation divergence of within‐sex ( G ) to cross‐sex ( B ) covariance matrices. In line with a previous theoretical prediction, we find that the cross‐sex covariance matrix, B , is more variable than either within‐sex G matrix. Decomposition of B matrices into their symmetrical and nonsymmetrical components revealed that instability is linked to the degree of asymmetry. We also find that the degree of asymmetry correlates with latitude suggesting a role for spatially varying natural selection in shaping genetic constraints on the evolution of sexual dimorphism.     

Constraints on the evolution of males and sexual dimorphism: field estimates of genetic architecture of reproductive traits in three populations of gynodioecious Fragaria virginiana

Abstract To understand how genetic constraints may limit the evolution of males and sexual dimorphism in a gynodioecious species, I conducted a quantitative genetic experiment in a gynodioecious wild strawberry, Fragaria virginiana . I estimated and compared genetic parameters (narrow-sense heritabilities, between-trait and between-sex genetic correlations, as well as phenotypic and genetic variance-covariance matrices) in the two sex morphs from three populations grown in a common field garden. I measured pollen and ovule production per flower, petal size, fruit set, and flower number. My major findings are as follows. (1) The presence of a phenotypic trade-off between pollen production and fruit set in hermaphrodites reflects a negative genetic correlation in the narrow sense that is statistically significant when pooled across populations. (2) The main constraints on the evolution of males are low genetic variation for pollen per flower and strong positive correlations associated with ovule number (e.g., between pollen and ovules in hermaphrodites, and between ovules in hermaphrodites and females). (3) Traits with the lowest levels of sexual dimorphism (ovule number and flower number) have the highest between-sex genetic correlations suggesting that overlap in the expression of genes in the sex morphs constrains their independent evolution. (4) There are significant differences in G matrices between sex morphs but not among populations. However, evidence that male-female trait correlations in hermaphrodites were lower in populations with higher frequencies of females may indicate subtle changes in genetic architecture.     

The constancy of the G matrix through species divergence and the effects of quantitative genetic constraints on phenotypic evolution: a case study in crickets

Long-term phenotypic evolution can be modeled using the response-to-selection equation of quantitative genetics, which incorporates information about genetic constraints (the G matrix). However, little is known about the evolution of G and about its long-term importance in constraining phenotypic evolution. We first investigated the degree of conservation of the G matrix across three species of crickets and qualitatively compared the pattern of variation of G to the phylogeny of the group. Second, we investigated the effect of G on phenotypic evolution by comparing the direction of greatest quantitative genetic variation within species (g(max)) to the direction of phenotypic divergence between species (Delta(z)). Each species, Gryllus veletis, G. firmus, and G. pennsylvanicus, was reared in the laboratory using a full-sib breeding design to extract quantitative genetic information. Five morphological traits related to size were measured. G matrices were compared using three statistical approaches: the T method, the Flury hierarchy, and the MANOVA method. Results revealed that the differences between matrices were small and mostly caused by differences in the magnitude of the genetic variation, not by differences in principal component structure. This suggested that the G matrix structure of this group of species was preserved, despite significant phenotypic divergence across species. The small observed differences in G matrices across species were qualitatively consistent with genetic distances, whereas ecological information did not provide a good prediction of G matrix variation. The comparison of g(max) and Delta(z) revealed that the angle between these two vectors was small in two of three species comparisons, whereas the larger angle corresponding to the third species comparison was caused in large part by one of the five traits. This suggests that multivariate phenotypic divergence occurred mostly in a direction predicted by the direction of greatest genetic variation, although it was not possible to demonstrate the causal relationship from G to Delta(z). Overall, this study provided some support for the validity of the predictive power of quantitative genetics over evolutionary time scales.     

Phenotypic Plasticity in the Annual Weed Polygonum aviculare

Reaction norms of fourteen life history and morphological traits were investigated in four tetra- and two hexaploid genotypes of the annual weed species complex, Polygonum aviculare. The plants were cultivated in six treatments consisting of factorial combinations of three pot sizes and two fertility levels. All characters, except life span, were plastic but the relative importance of genotype (G), treatment (T) and interaction (G × T) to total variance was strongly trait-specific. Consistent genetic differentiation, not correlated with ploidy level, was found in metamer size and life history: genotypes originating from trampled sites had smaller metamers and shorter shoots while those originating from sites with a short growing season, due to weeding activities, had a shorter life span, an earlier flowering date and a higher biomass allocation to reproduction compared to genotypes from less disturbed sites. Significant variation was found in reaction norms for all characters, including a lower amount of plasticity in metamer size in genotypes with numerous metamers and a lower amount of plasticity in total weight in shortlived genotypes. This suggested that variation in phenotypic plasticity reflected developmental constraints imposed by contrasting life span and metamer size in different genotypes. There was no evidence for niche differentiation along the soil resource gradient, suggesting that the species is comprised of “general purpose” genotypes with respect to soil fertility. It is concluded that the Polygonum aviculare complex has evolved a “dual” adaptive strategy i.e. a combination of genetic polymorphism and high phenotypic plasticity.     

The role of phenotypic plasticity in driving genetic evolution

Models of population divergence and speciation are often based on the assumption that differences between populations are due to genetic factors, and that phenotypic change is due to natural selection. It is equally plausible that some of the differences among populations are due to phenotypic plasticity. We use the metaphor of the adaptive landscape to review the role of phenotypic plasticity in driving genetic evolution. Moderate levels of phenotypic plasticity are optimal in permitting population survival in a new environment and in bringing populations into the realm of attraction of an adaptive peak. High levels of plasticity may increase the probability of population persistence but reduce the likelihood of genetic change, because the plastic response itself places the population close to a peak. Moderate levels of plasticity arise whenever multiple traits, some of which are plastic and others not, form a composite trait involved in the adaptive response. For example, altered behaviours may drive selection on morphology and physiology. Because there is likely to be a considerable element of chance in which behaviours become established, behavioural change followed by morphological and physiological evolution may be a potent force in driving evolution in novel directions. We assess the role of phenotypic plasticity in stimulating evolution by considering two examples from birds: (i) the evolution of red and yellow plumage coloration due to carotenoid consumption; and (ii) the evolution of foraging behaviours on islands. Phenotypic plasticity is widespread in nature and may speed up, slow down, or have little effect on evolutionary change. Moderate levels of plasticity may often facilitate genetic evolution but careful analyses of individual cases are needed to ascertain whether plasticity has been essential or merely incidental to population differentiation.     

Plasticity of functional traits varies clinally along a rainfall gradient in Eucalyptus tricarpa

Widespread species often occur across a range of climatic conditions, through a combination of local genetic adaptations and phenotypic plasticity. Species with greater phenotypic plasticity are likely to be better positioned to cope with rapid anthropogenic climate changes, while those displaying strong local adaptations might benefit from translocations to assist the movement of adaptive genes as the climate changes. Eucalyptus tricarpa occurs across a climatic gradient in south‐eastern Australia, a region of increasing aridity, and we hypothesized that this species would display local adaptation to climate. We measured morphological and physiological traits reflecting climate responses in nine provenances from sites of 460 to 1040 mm annual rainfall, in their natural habitat and in common gardens near each end of the gradient. Local adaptation was evident in functional traits and differential growth rates in the common gardens. Some traits displayed complex combinations of plasticity and genetic divergence among provenances, including clinal variation in plasticity itself. Provenances from drier locations were more plastic in leaf thickness, whereas leaf size was more plastic in provenances from higher rainfall locations. Leaf density and stomatal physiology (as indicated by δ¹³C and δ¹⁸O) were highly and uniformly plastic. In addition to variation in mean trait values, genetic variation in trait plasticity may play a role in climate adaptation.     

Mixed support for an alignment between phenotypic plasticity and genetic differentiation in damselfly wing shape

The relationship between genetic differentiation and phenotypic plasticity can provide information on whether plasticity generally facilitates or hinders adaptation to environmental change. Here, we studied wing shape variation in a damselfly (Lestes sponsa) across a latitudinal gradient in Europe that differed in time constraints mediated by photoperiod and temperature. We reared damselflies from northern and southern populations in the laboratory using a reciprocal transplant experiment that simulated time-constrained (i.e. northern) and unconstrained (southern) photoperiods and temperatures. After emergence, adult wing shape was analysed using geometric morphometrics. Wings from individuals in the northern and southern populations differed significantly in shape when animals were reared in their respective native environment. Comparing wing shape across environments, we found evidence for phenotypic plasticity in wing shape, and this response differed across populations (i.e. G × E interactions). This interaction was driven by a stronger plastic response by individuals from the northern population and differences in the direction of plastic wing shape changes among populations. The alignment between genetic and plastic responses depended on the specific combination of population and rearing environment. For example, there was an alignment between plasticity and genetic differentiation under time-constrained, but not under non-time-constrained conditions for forewings. We thus find mixed support for the hypothesis that environmental plasticity and genetic population differentiation are aligned. Furthermore, although our laboratory treatments mimicked the natural climatic conditions at northern and southern latitudes, the effects of population differences on wing shape were two to four times stronger than plastic effects. We discuss our results in terms of time constraints and the possibility that natural and sexual selection is acting differently on fore- and hindwings.