The genetic variance but not the genetic covariance of life‐history traits changes towards the north in a time‐constrained insect

Seasonal time constraints are usually stronger at higher than lower latitudes and can exert strong selection on life‐history traits and the correlations among these traits. To predict the response of life‐history traits to environmental change along a latitudinal gradient, information must be obtained about genetic variance in traits and also genetic correlation between traits, that is the genetic variance‐covariance matrix, G . Here, we estimated G for key life‐history traits in an obligate univoltine damselfly that faces seasonal time constraints. We exposed populations to simulated native temperatures and photoperiods and common garden environmental conditions in a laboratory set‐up. Despite differences in genetic variance in these traits between populations (lower variance at northern latitudes), there was no evidence for latitude‐specific covariance of the life‐history traits. At simulated native conditions, all populations showed strong genetic and phenotypic correlations between traits that shaped growth and development. The variance–covariance matrix changed considerably when populations were exposed to common garden conditions compared with the simulated natural conditions, showing the importance of environmentally induced changes in multivariate genetic structure. Our results highlight the importance of estimating variance–covariance matrixes in environments that mimic selection pressures and not only trait variances or mean trait values in common garden conditions for understanding the trait evolution across populations and environments.     

QUANTITATIVE GENETIC DIVERGENCE AND STANDING GENETIC (CO)VARIANCE IN THERMAL REACTION NORMS ALONG LATITUDE

Although the potential to adapt to warmer climate is constrained by genetic trade‐offs, our understanding of how selection and mutation shape genetic (co)variances in thermal reaction norms is poor. Using 71 isofemale lines of the fly Sepsis punctum, originating from northern, central, and southern European climates, we tested for divergence in juvenile development rate across latitude at five experimental temperatures. To investigate effects of evolutionary history in different climates on standing genetic variation in reaction norms, we further compared genetic (co)variances between regions. Flies were reared on either high or low food resources to explore the role of energy acquisition in determining genetic trade‐offs between different temperatures. Although the latter had only weak effects on the strength and sign of genetic correlations, genetic architecture differed significantly between climatic regions, implying that evolution of reaction norms proceeds via different trajectories at high latitude versus low latitude in this system. Accordingly, regional genetic architecture was correlated to region‐specific differentiation. Moreover, hot development temperatures were associated with low genetic variance and stronger genetic correlations compared to cooler temperatures. We discuss the evolutionary potential of thermal reaction norms in light of their underlying genetic architectures, evolutionary histories, and the materialization of trade‐offs in natural environments.     

Similarity in G matrix structure among natural populations of Arabidopsis lyrata

Understanding the stability of the G matrix in natural populations is fundamental for predicting evolutionary trajectories; yet, the extent of its spatial variation and how this impacts responses to selection remain open questions. With a nested paternal half‐sib crossing design and plants grown in a field experiment, we examined differences in the genetic architecture of flowering time, floral display, and plant size among four Scandinavian populations of Arabidopsis lyrata. Using a multivariate Bayesian framework, we compared the size, shape, and orientation of G matrices and assessed their potential to facilitate or constrain trait evolution. Flowering time, floral display and rosette size varied among populations and significant additive genetic variation within populations indicated potential to evolve in response to selection. Yet, some characters, including flowering start and number of flowers, may not evolve independently because of genetic correlations. Using a multivariate framework, we found few differences in the genetic architecture of traits among populations. G matrices varied mostly in size rather than shape or orientation. Differences in multivariate responses to selection predicted from differences in G were small, suggesting overall matrix similarity and shared constraints to trait evolution among populations.     

Characterizing the evolution of genetic variance using genetic covariance tensors

Determining how genetic variance changes under selection in natural populations has proved to be a very resilient problem in evolutionary genetics. In the same way that understanding the availability of genetic variance within populations requires the simultaneous consideration of genetic variance in sets of functionally related traits, determining how genetic variance changes under selection in natural populations will require ascertaining how genetic variance–covariance (G) matrices evolve. Here, we develop a geometric framework using higher-order tensors, which enables the empirical characterization of how G matrices have diverged among populations. We then show how divergence among populations in genetic covariance structure can then be associated with divergence in selection acting on those traits using key equations from evolutionary theory. Using estimates of G matrices of eight male sexually selected traits from nine geographical populations of Drosophila serrata, we show that much of the divergence in genetic variance occurred in a single trait combination, a conclusion that could not have been reached by examining variation among the individual elements of the nine G matrices. Divergence in G was primarily in the direction of the major axes of genetic variance within populations, suggesting that genetic drift may be a major cause of divergence in genetic variance among these populations.     

Population differentiation in G matrix structure due to natural selection in Rana temporaria

The additive genetic variance-covariance matrix (G) is a concept central to discussions about evolutionary change over time in a suite of traits. However, at the moment we do not know how fast G itself changes as a consequence of selection or how sensitive it is to environmental influences. We investigated possible evolutionary divergence and environmental influences on G using data from a factorial common-garden experiment where common frog (Rana temporaria) tadpoles from two divergent populations were exposed to three different environmental treatments. G-matrices were estimated using an animal model approach applied to data from a NCII breeding design. Matrix comparisons using both Flury and multivariate analysis of variance methods revealed significant differences in G matrices both between populations and between treatments within populations, the former being generally larger than the latter. Comparison of levels of population differentiation in trait means using Q(ST) indices with that observed in microsatellite markers (F(ST)) revealed that the former values generally exceeded the neutral expectation set by F(ST). Hence, the results suggest that intraspecific divergence in G matrix structure has occurred mainly due to natural selection.     

Parallel and nonparallel ecological,morphological and genetic divergence in lake–stream stickleback from a single catchment

Parallel phenotypic evolution in similar environments has been well studied in evolutionary biology; however, comparatively little is known about the influence of determinism and historical contingency on the nature, extent and generality of this divergence. Taking advantage of a novel system containing multiple lake–stream stickleback populations, we examined the extent of ecological, morphological and genetic divergence between three‐spined stickleback present in parapatric environments. Consistent with other lake–stream studies, we found a shift towards a deeper body and shorter gill rakers in stream fish. Morphological shifts were concurrent with changes in diet, indicated by both stable isotope and stomach contents analysis. Performing a multivariate test for shared and unique components of evolutionary response to the distance gradient from the lake, we found a strong signature of parallel adaptation. Nonparallel divergence was also present, attributable mainly to differences between river locations. We additionally found evidence of genetic substructuring across five lake–stream transitions, indicating that some level of reproductive isolation occurs between populations in these habitats. Strong correlations between pairwise measures of morphological, ecological and genetic distance between lake and stream populations supports the hypothesis that divergent natural selection between habitats drives adaptive divergence and reproductive isolation. Lake–stream stickleback divergence in Lough Neagh provides evidence for the deterministic role of selection and supports the hypothesis that parallel selection in similar environments may initiate parallel speciation.     

Quantitative genetic and translocation experiments reveal genotype‐by‐environment effects on juvenile life‐history traits in two populations of Chinook salmon (Oncorhynchus tshawytscha)

Understanding the genetic architecture of phenotypic plasticity is required to assess how populations might respond to heterogeneous or changing environments. Although several studies have examined population‐level patterns in environmental heterogeneity and plasticity, few studies have examined individual‐level variation in plasticity. Here, we use the North Carolina II breeding design and translocation experiments between two populations of Chinook salmon to detail the genetic architecture and plasticity of offspring survival and growth. We followed the survival of 50 800 offspring through the larval stage and used parentage analysis to examine survival and growth through freshwater rearing. In one population, we found that additive genetic, nonadditive genetic and maternal effects explained 25%, 34% and 55% of the variance in larvae survival, respectively. In the second population, these effects explained 0%, 24% and 61% of the variance in larvae survival. In contrast, fry survival was regulated primarily by additive genetic effects, which indicates a shift from maternal to genetic effects as development proceeds. Fry growth also showed strong additive genetic effects. Translocations between populations revealed that offspring survival and growth varied between environments, the degree of which differed among families. These results indicate genetic differences among individuals in their degree of plasticity and consequently their ability to respond to environmental variation.     

Environmentally induced changes in correlated responses to selection reveal variable pleiotropy across a complex genetic network

Selection in novel environments can lead to a coordinated evolutionary response across a suite of characters. Environmental conditions can also potentially induce changes in the genetic architecture of complex traits, which in turn could alter the pattern of the multivariate response to selection. We describe a factorial selection experiment using the nematode Caenorhabditis remanei in which two different stress‐related phenotypes (heat and oxidative stress resistance) were selected under three different environmental conditions. The pattern of covariation in the evolutionary response between phenotypes or across environments differed depending on the environment in which selection occurred, including asymmetrical responses to selection in some cases. These results indicate that variation in pleiotropy across the stress response network is highly sensitive to the external environment. Our findings highlight the complexity of the interaction between genes and environment that influences the ability of organisms to acclimate to novel environments. They also make clear the need to identify the underlying genetic basis of genetic correlations in order understand how patterns of pleiotropy are distributed across complex genetic networks.     

Genetic correlations and their dependence on environmental similarity—Insights from livestock data

Genetic correlations for a trait across environments are predicted to decrease as environments diverge. However, estimates of genetic correlations from natural populations are typically defined across a limited environmental range and prone to very large standard errors, making it difficult to test this prediction. We address the importance of environmental distance on genetic correlations by employing data from domestic cattle in which abundant and accurate estimates are available from a wide range of environments. Three production traits related to milk yield show a clear decrease in genetic correlations with increasing environmental divergence. This pattern was also evident for growth traits and other yield traits but not for traits related to reproduction, morphology, physiology, or disease. We suspect that this reflects weaker selection on these latter trait classes compared to production traits, or alternatively the effects of selection are constrained by unfavorable genetic correlations between traits. The results support the notion that traits that historically have been under strong directional selection in a small range of frequently encountered environments will evolve high genetic correlations across these environments, while exposure to uncommon (and dissimilar) environments lead to a reranking of gene effects and a decrease in genetic correlations across environments.     

Genetic architecture of survival and fitness-related traits in two populations of Atlantic salmon

The additive genetic effects of traits can be used to predict evolutionary trajectories, such as responses to selection. Non-additive genetic and maternal environmental effects can also change evolutionary trajectories and influence phenotypes, but these effects have received less attention by researchers. We partitioned the phenotypic variance of survival and fitness-related traits into additive genetic, non-additive genetic and maternal environmental effects using a full-factorial breeding design within two allopatric populations of Atlantic salmon (Salmo salar). Maternal environmental effects were large at early life stages, but decreased during development, with non-additive genetic effects being most significant at later juvenile stages (alevin and fry). Non-additive genetic effects were also, on average, larger than additive genetic effects. The populations, generally, did not differ in the trait values or inferred genetic architecture of the traits. Any differences between the populations for trait values could be explained by maternal environmental effects. We discuss whether the similarities in architectures of these populations is the result of natural selection across a common juvenile environment.     

Protein deprivation facilitates the independent evolution of behavior and morphology

Ecological conditions such as nutrition can change genetic covariances between traits and accelerate or slow down trait evolution. As adaptive trait correlations can become maladaptive following rapid environmental change, poor or stressful environments are expected to weaken genetic covariances, thereby increasing the opportunity for independent evolution of traits. Here, we demonstrate the differences in genetic covariance among multiple behavioral and morphological traits (exploration, aggression, and body weight) between southern field crickets (Gryllus bimaculatus) raised in favorable (free‐choice) versus stressful (protein‐deprived) nutritional environments. We also quantify the extent to which differences in genetic covariance structures contribute to the potential for the independent evolution of these traits. We demonstrate that protein‐deprived environments tend to increase the potential for traits to evolve independently, which is caused by genetic covariances that are significantly weaker for crickets raised on protein‐deprived versus free‐choice diets. The weakening effects of stressful environments on genetic covariances tended to be stronger in males than in females. The weakening of the genetic covariance between traits under stressful nutritional environments was expected to facilitate the opportunity for adaptive evolution across generations. Therefore, the multivariate gene‐by‐environment interactions revealed here may facilitate behavioral and morphological adaptations to rapid environmental change.     

The extent and genetic basis of phenotypic divergence in life history traits in Mimulus guttatus

Differential natural selection acting on populations in contrasting environments often results in adaptive divergence in multivariate phenotypes. Multivariate trait divergence across populations could be caused by selection on pleiotropic alleles or through many independent loci with trait‐specific effects. Here, we assess patterns of association between a suite of traits contributing to life history divergence in the common monkey flower, Mimulus guttatus, and examine the genetic architecture underlying these correlations. A common garden survey of 74 populations representing annual and perennial strategies from across the native range revealed strong correlations between vegetative and reproductive traits. To determine whether these multitrait patterns arise from pleiotropic or independent loci, we mapped QTLs using an approach combining high‐throughput sequencing with bulk segregant analysis on a cross between populations with divergent life histories. We find extensive pleiotropy for QTLs related to flowering time and stolon production, a key feature of the perennial strategy. Candidate genes related to axillary meristem development colocalize with the QTLs in a manner consistent with either pleiotropic or independent QTL effects. Further, these results are analogous to previous work showing pleiotropy‐mediated genetic correlations within a single population of M. guttatus experiencing heterogeneous selection. Our findings of strong multivariate trait associations and pleiotropic QTLs suggest that patterns of genetic variation may determine the trajectory of adaptive divergence.     

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.     

Structure and stability of genetic variance–covariance matrices: A Bayesian sparse factor analysis of transcriptional variation in the three‐spined stickleback

The genetic variance–covariance matrix ( G ) is a quantity of central importance in evolutionary biology due to its influence on the rate and direction of multivariate evolution. However, the predictive power of empirically estimated G ‐matrices is limited for two reasons. First, phenotypes are high‐dimensional, whereas traditional statistical methods are tuned to estimate and analyse low‐dimensional matrices. Second, the stability of G to environmental effects and over time remains poorly understood. Using Bayesian sparse factor analysis (BSFG) designed to estimate high‐dimensional G ‐matrices, we analysed levels variation and covariation in 10,527 expressed genes in a large (n = 563) half‐sib breeding design of three‐spined sticklebacks subject to two temperature treatments. We found significant differences in the structure of G between the treatments: heritabilities and evolvabilities were higher in the warm than in the low‐temperature treatment, suggesting more and faster opportunity to evolve in warm (stressful) conditions. Furthermore, comparison of G and its phenotypic equivalent P revealed the latter is a poor substitute of the former. Most strikingly, the results suggest that the expected impact of G on evolvability—as well as the similarity among G ‐matrices—may depend strongly on the number of traits included into analyses. In our results, the inclusion of only few traits in the analyses leads to underestimation in the differences between the G ‐matrices and their predicted impacts on evolution. While the results highlight the challenges involved in estimating G , they also illustrate that by enabling the estimation of large G ‐matrices, the BSFG method can improve predicted evolutionary responses to selection.     

Revealing hidden evolutionary capacity to cope with global change

The extent to which global change will impact the long‐term persistence of species depends on their evolutionary potential to adapt to future conditions. While the number of studies that estimate the standing levels of adaptive genetic variation in populations under predicted global change scenarios is growing all the time, few studies have considered multiple environments simultaneously and even fewer have considered evolutionary potential in multivariate context. Because conditions will not be constant, adaptation to climate change is fundamentally a multivariate process so viewing genetic variances and covariances over multivariate space will always be more informative than relying on bivariate genetic correlations between traits. A multivariate approach to understanding the evolutionary capacity to cope with global change is necessary to avoid misestimating adaptive genetic variation in the dimensions in which selection will act. We assessed the evolutionary capacity of the larval stage of the marine polychaete Galeolaria caespitosa to adapt to warmer water temperatures. Galeolaria is an important habitat‐forming species in Australia, and its earlier life‐history stages tend to be more susceptible to stress. We used a powerful quantitative genetics design that assessed the impacts of three temperatures on subsequent survival across over 30 000 embryos across 204 unique families. We found adaptive genetic variation in the two cooler temperatures in our study, but none in the warmest temperature. Based on these results, we would have concluded that this species has very little capacity to evolve to the warmest temperature. However, when we explored genetic variation in multivariate space, we found evidence that larval survival has the potential to evolve even in the warmest temperatures via correlated responses to selection across thermal environments. Future studies should take a multivariate approach to estimating evolutionary capacity to cope with global change lest they misestimate a species’ true adaptive potential.     

Expression of additive genetic variances and covariances for wild radish floral traits: comparison between field and greenhouse environments

Measurements of the genetic variation and covariation underlying quantitative traits are crucial to our understanding of current evolutionary change and the mechanisms causing this evolution. This fact has spurred a large number of studies estimating heritabilities and genetic correlations in a variety of organisms. Most of these studies have been done in laboratory or greenhouse settings, but it is not well known how accurately these measurements estimate genetic variance and covariance expressed in the field. We conducted a quantitative genetic half-sibling analysis on six floral traits in wild radish. Plants were grown from seed in the field and were exposed to natural environmental variation throughout their lives, including herbivory and intra- and interspecific competition. The estimates of heritabilities and the additive genetic variance-covariance matrix (G) obtained from this analysis were then compared to previous greenhouse estimates of the same floral traits from the same natural population. Heritabilities were much lower in the field for all traits, and this was due to both large increases in environmental variance and decreases in additive genetic variance. Additive genetic covariance expressed was also much lower in the field. These differences resulted in highly significant differences in the G matrix between the greenhouse and field environments using two complementary testing methods. Although the G matrices shared some principal components in common, they were not simply proportional to each other. Therefore, the greenhouse results did not accurately depict how the floral traits would respond to natural selection in the field.     

Density dependence and population differentiation of genetic architecture in Impatiens capensis in natural environments

We identified environment-dependent constraints on the evolution of plasticity to density under natural conditions in two natural populations of Impatiens capensis. We also examined the expression of population divergence in genetic variance-covariance matrices in these natural environments. Inbred lines, originally collected from a sunny site with high seedling densities and a woodland site with low seedling densities, were planted in both original sites at natural high densities and at low density. Morphological and life-history characters were measured. More genetic variation for plastic responses to density was expressed in the sun site than in the woodland site, so the evolutionary potential of plasticity was greater in the sun site. Strong genetic correlations between the same character expressed at different densities and correlations among different characters could constrain the evolution of plasticity in both sites. Genetically based trade-offs in meristem allocation to vegetative growth and reproduction were apparent only in the high-resource environment with no overhead canopy and no intraspecific competition. Therefore, genetic constraints on the evolution of plasticity depended on the site and density in which plants were grown, and correlated responses to selection on plastic characters are also expected to differ between sites and densities. Population differentiation in genetic variance-covariance matrices was detected, but matrix structural differences, as opposed to proportional differences, were detected between populations only in the sun site at natural high density. Thus, population divergence in genetic architecture can occur rapidly and on a fine spatial scale, but the expression of such divergence may depend on the environment.     

Phenotypic variation and covariation indicate high evolvability of acoustic communication in crickets

Studying the genetic architecture of sexual traits provides insight into the rate and direction at which traits can respond to selection. Traits associated with few loci and limited genetic and phenotypic constraints tend to evolve at high rates typically observed for secondary sexual characters. Here, we examined the genetic architecture of song traits and female song preferences in the field crickets Gryllus rubens and Gryllus texensis. Song and preference data were collected from both species and interspecific F1 and F2 hybrids. We first analysed phenotypic variation to examine interspecific differentiation and trait distributions in parental and hybrid generations. Then, the relative contribution of additive and additive‐dominance variation was estimated. Finally, phenotypic variance–covariance ( P ) matrices were estimated to evaluate the multivariate phenotype available for selection. Song traits and preferences had unimodal trait distributions, and hybrid offspring were intermediate with respect to the parents. We uncovered additive and dominance variation in song traits and preferences. For two song traits, we found evidence for X‐linked inheritance. On the one hand, the observed genetic architecture does not suggest rapid divergence, although sex linkage may have allowed for somewhat higher evolutionary rates. On the other hand, P matrices revealed that multivariate variation in song traits aligned with major dimensions in song preferences, suggesting a strong selection response. We also found strong covariance between the main traits that are sexually selected and traits that are not directly selected by females, providing an explanation for the striking multivariate divergence in male calling songs despite limited divergence in female preferences.     

Genetic covariances promote climatic adaptation in Australian Drosophila*

Evolutionary potential for adaptation hinges upon the orientation of genetic variation for traits under selection, captured by the additive genetic variance-covariance matrix (G), as well as the evolutionary stability of G. Yet studies that assess both the stability of G and its alignment with selection are extraordinarily rare. We evaluated the stability of G in three Drosophila melanogaster populations that have adapted to local climatic conditions along a latitudinal cline. We estimated population- and sex-specific G matrices for wing size and three climatic stress-resistance traits that diverge adaptively along the cline. To determine how G affects evolutionary potential within these populations, we used simulations to quantify how well G aligns with the direction of trait divergence along the cline (as a proxy for the direction of local selection) and how genetic covariances between traits and sexes influence this alignment. We found that G was stable across the cline, showing no significant divergence overall, or in sex-specific subcomponents, among populations. G also aligned well with the direction of clinal divergence, with genetic covariances strongly elevating evolutionary potential for adaptation to climatic extremes. These results suggest that genetic covariances between both traits and sexes should significantly boost evolutionary responses to environmental change.