Sexual variation in heritability and genetic correlations of morphological traits in house sparrow (Passer domesticus)

Estimates of genetic components are important for our understanding of how individual characteristics are transferred between generations. We show that the level of heritability varies between 0.12 and 0.68 in six morphological traits in house sparrows (Passer domesticus L.) in northern Norway. Positive and negative genetic correlations were present among traits, suggesting evolutionary constraints on the evolution of some of these characters. A sexual difference in the amount of heritable genetic variation was found in tarsus length, wing length, bill depth and body condition index, with generally higher heritability in females. In addition, the structure of the genetic variance-covariance matrix for the traits differed between the sexes. Genetic correlations between males and females for the morphological traits were however large and not significantly different from one, indicating that sex-specific responses to selection will be influenced by intersexual differences in selection differentials. Despite this, some traits had heritability above 0.1 in females, even after conditioning on the additive genetic covariance between sexes and the additive genetic variances in males. Moreover, a meta-analysis indicated that higher heritability in females than in males may be common in birds. Thus, this indicates sexual differences in the genetic architecture of birds. Consequently, as in house sparrows, the evolutionary responses to selection will often be larger in females than males. Hence, our results suggest that sex-specific additive genetic variances and covariances, although ignored in most studies, should be included when making predictions of evolutionary changes from standard quantitative genetic models.     

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.     

Cranial evolution in sakis (Pithecia, Platyrrhini). II: Evolutionary processes and morphological integration

Patterns of interspecific differentiation in saki monkeys (Pithecia) are quantitatively described and possible evolutionary processes producing them are examined. The comparison of species correlation matrices to expected patterns of morphological integration reveal significant and similar patterns of development-based cranial integration among species. Aspects of the facial region are more heavily influenced by general size variation than features of the neural region. The comparison of pooled within- and between-groups V/CV matrices suggests that genetic drift might be a sufficient explanation for saki cranial evolution. Differential natural selection gradients are also reconstructed because selection may also have caused population differentiation through evolutionary time. These gradients illustrate the inherent multivariate nature of selection, being a consequence of the interaction between existing morphological integration (correlation) among traits and the action of natural selection. Yet, our attempt to interpret selection gradients in terms of their functional significance did not result in any clear association between selection and function. Perhaps this is also an indication that morphological evolution in sakis was mostly neutral.     

The G matrix under fluctuating correlational mutation and selection

Theoretical quantitative genetics provides a framework for reconstructing past selection and predicting future patterns of phenotypic differentiation. However, the usefulness of the equations of quantitative genetics for evolutionary inference relies on the evolutionary stability of the additive genetic variance-covariance matrix (G matrix). A fruitful new approach for exploring the evolutionary dynamics of G involves the use of individual-based computer simulations. Previous studies have focused on the evolution of the eigenstructure of G. An alternative approach employed in this paper uses the multivariate response-to-selection equation to evaluate the stability of G. In this approach, I measure similarity by the correlation between response-to-selection vectors due to random selection gradients. I analyze the dynamics of G under several conditions of correlational mutation and selection. As found in a previous study, the eigenstructure of G is stabilized by correlational mutation and selection. However, over broad conditions, instability of G did not result in a decreased consistency of the response to selection. I also analyze the stability of G when the correlation coefficients of correlational mutation and selection and the effective population size change through time. To my knowledge, no prior study has used computer simulations to investigate the stability of G when correlational mutation and selection fluctuate. Under these conditions, the eigenstructure of G is unstable under some simulation conditions. Different results are obtained if G matrix stability is assessed by eigenanalysis or by the response to random selection gradients. In this case, the response to selection is most consistent when certain aspects of the eigenstructure of G are least stable and vice versa.     

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.     

The alignment between phenotypic plasticity,the major axis of genetic variation and the response to selection

Phenotypic plasticity is the ability of a genotype to produce more than one phenotype in order to match the environment. Recent theory proposes that the major axis of genetic variation in a phenotypically plastic population can align with the direction of selection. Therefore, theory predicts that plasticity directly aids adaptation by increasing genetic variation in the direction favoured by selection and reflected in plasticity. We evaluated this theory in the freshwater crustacean Daphnia pulex, facing predation risk from two contrasting size-selective predators. We estimated plasticity in several life-history traits, the G matrix of these traits, the selection gradients on reproduction and survival, and the predicted responses to selection. Using these data, we tested whether the genetic lines of least resistance and the predicted response to selection aligned with plasticity. We found predator environment-specific G matrices, but shared genetic architecture across environments resulted in more constraint in the G matrix than in the plasticity of the traits, sometimes preventing alignment of the two. However, as the importance of survival selection increased, the difference between environments in their predicted response to selection increased and resulted in closer alignment between the plasticity and the predicted selection response. Therefore, plasticity may indeed aid adaptation to new environments.     

Early diversification of sperm size in the evolutionary history of the old world leaf warblers (Phylloscopidae)

Sperm morphological traits are highly variable among species and are commonly thought to evolve by post‐copulatory sexual selection. However, little is known about the evolutionary dynamics of sperm morphology, and whether rates of evolutionary change are variable over time and among taxonomic groups. Here, we examine sperm morphology from 21 species of Old World leaf warblers (Phylloscopidae), a group of generally dull, sexually monochromatic birds, which are known to have high levels of extra‐pair paternity. We found that sperm length differs markedly across species, spanning about 40% of the range observed across a larger selection of passerine birds. Furthermore, we found strong support for an ‘early‐burst’ model of trait evolution, implying that the majority of divergence in sperm length has occurred early in the evolutionary history of this clade with subsequent evolutionary stasis. This large early divergence matches the early divergence reported in ecological traits (i.e. body size and feeding behaviour). Our findings demonstrate that rates of evolution in sperm morphology can change over time in passerine taxa, and that evolutionary stasis in sperm traits can occur even in species exhibiting characteristics consistent with moderate‐to‐high levels of sperm competition. It remains a major challenge to identify the selection mechanisms and possible constraints responsible for these variable rates of sperm evolution.     

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.     

MULTIYEAR STUDY OF MULTIVARIATE LINEAR AND NONLINEAR PHENOTYPIC SELECTION ON FLORAL TRAITS OF HUMMINGBIRD-POLLINATED SILENE VIRGINICA

Pollination syndromes suggest that convergent evolution of floral traits and trait combinations reflects similar selection pressures. Accordingly, a pattern of selection on floral traits is expected to be consistent with increasing the attraction and pollen transfer of the important pollinator. We measured individual variation in six floral traits and yearly and lifetime total plant seed and fruit production of 758 plants across nine years of study in natural populations of Ruby-Throated Hummingbird-pollinated Silene virginica. The type, strength, and direction of selection gradients were observed by year, and for two cohorts selection was estimated through lifetime maternal fitness. Positive directional selection was detected on floral display height in all years of study and stigma exsertion in all years but one. Significant quadratic and correlational selection gradients were rare. However, a canonical analysis of the gamma matrix indicated nonlinear selection was common; if significant curvature was detected it was convex with one exception. Our analyses demonstrated selection favored trait combinations and the integration of floral features of attraction and pollen transfer efficiency that were consistent with the hummingbird pollination syndrome.     

Divergence with gene flow and fine-scale phylogeographical structure in the wedge-billed woodcreeper, Glyphorynchus spirurus, a Neotropical rainforest bird

Determining the relative roles of vicariance and selection in restricting gene flow between populations is of central importance to the evolutionary process of population divergence and speciation. Here we use molecular and morphological data to contrast the effect of isolation (by mountains and geographical distance) with that of ecological factors (altitudinal gradients) in promoting differentiation in the wedge-billed woodcreeper, Glyphorynchus spirurus , a tropical forest bird, in Ecuador. Tarsus length and beak size increased relative to body size with altitude on both sides of the Andes, and were correlated with the amount of moss on tree trunks, suggesting the role of selection in driving adaptive divergence. In contrast, molecular data revealed a considerable degree of admixture along these altitudinal gradients, suggesting that adaptive divergence in morphological traits has occurred in the presence of gene flow. As suggested by mitochondrial DNA sequence data, the Andes act as a barrier to gene flow between ancient subspecific lineages. Genome-wide amplified fragment length polymorphism markers reflected more recent patterns of gene flow and revealed fine-scale patterns of population differentiation that were not detectable with mitochondrial DNA, including the differentiation of isolated coastal populations west of the Andes. Our results support the predominant role of geographical isolation in driving genetic differentiation in G. spirurus , yet suggest the role of selection in driving parallel morphological divergence along ecological gradients.     

A test of the hypothesis that correlational selection generates genetic correlations

Theory predicts that correlational selection on two traits will cause the major axis of the bivariate G matrix to orient itself in the same direction as the correlational selection gradient. Two testable predictions follow from this: for a given pair of traits, (1) the sign of correlational selection gradient should be the same as that of the genetic correlation, and (2) the correlational selection gradient should be positively correlated with the value of the genetic correlation. We test this hypothesis with a meta-analysis utilizing empirical estimates of correlational selection gradients and measures of the correlation between the two focal traits. Our results are consistent with both predictions and hence support the underlying hypothesis that correlational selection generates a genetic correlation between the two traits and hence orients the bivariate G matrix.     

Evidence for strong intralocus sexual conflict in the Indian meal moth, Plodia interpunctella

Males and females share a genome and express many shared phenotypic traits, which are often selected in opposite directions. This generates intralocus sexual conflict that may constrain trait evolution by preventing the sexes from reaching their optimal phenotype. Furthermore, if present across multiple loci, intralocus sexual conflict can result in a gender load that may diminish the benefits of sexual selection and help maintain genetic variation for fitness. Despite the importance of intralocus sexual conflict, surprisingly few empirical studies conclusively demonstrate its operation. We show that the pattern of multivariate selection acting on three sexually dimorphic life-history traits (development time, body size, and longevity) in the Indian meal moth, Plodia interpunctella, is opposing for the sexes. Moreover, we combined our estimates of selection with the additive genetic variance-covariance matrix (G) to predict the evolutionary response of the life-history traits in the sexes and showed that the angle between the vector of responses and the vector of sexually antagonistic selection was almost orthogonal at 84.70°. Thus, G biases the predicted response of life-history traits in the sexes away from the direction of sexually antagonistic selection, confirming the presence of strong intralocus sexual conflict in this species. Despite this, sexual dimorphism has evolved in all of the life-history traits examined suggesting that mechanism(s) have evolved to resolve this conflict and allow the sexes to reach their life-history optima. We argue that intralocus sexual conflict is likely to play an important role in the evolution of divergent life-history strategies between the sexes in this species.     

Evolution of body size in Galapagos marine iguanas

Body size is one of the most important traits of organisms and allows predictions of an individual's morphology, physiology, behaviour and life history. However, explaining the evolution of complex traits such as body size is difficult because a plethora of other traits influence body size. Here I review what we know about the evolution of body size in a group of island reptiles and try to generalize about the mechanisms that shape body size. Galapagos marine iguanas occupy all 13 larger islands in this Pacific archipelago and have maximum island body weights between 900 and 12 000g. The distribution of body sizes does not match mitochondrial clades, indicating that body size evolves independently of genetic relatedness. Marine iguanas lack intra- and inter-specific food competition and predators are not size-specific, discounting these factors as selective agents influencing body size. Instead I hypothesize that body size reflects the trade-offs between sexual and natural selection. We found that sexual selection continuously favours larger body sizes. Large males establish display territories and some gain over-proportional reproductive success in the iguanas' mating aggregations. Females select males based on size and activity and are thus responsible for the observed mating skew. However, large individuals are strongly selected against during El Ni?o-related famines when dietary algae disappear from the intertidal foraging areas. We showed that differences in algae sward ('pasture') heights and thermal constraints on large size are causally responsible for differences in maximum body size among populations. I hypothesize that body size in many animal species reflects a trade-off between foraging constraints and sexual selection and suggest that future research could focus on physiological and genetic mechanisms determining body size in wild animals. Furthermore, evolutionary stable body size distributions within populations should be analysed to better understand selection pressures on individual body size.     

Environmental effects on the structure of the G‐matrix

Genetic correlations between traits determine the multivariate response to selection in the short term, and thereby play a causal role in evolutionary change. Although individual studies have documented environmentally induced changes in genetic correlations, the nature and extent of environmental effects on multivariate genetic architecture across species and environments remain largely uncharacterized. We reviewed the literature for estimates of the genetic variance–covariance ( G ) matrix in multiple environments, and compared differences in G between environments to the divergence in G between conspecific populations (measured in a common garden). We found that the predicted evolutionary trajectory differed as strongly between environments as it did between populations. Between‐environment differences in the underlying structure of G (total genetic variance and the relative magnitude and orientation of genetic correlations) were equal to or greater than between‐population differences. Neither environmental novelty, nor the difference in mean phenotype predicted these differences in G . Our results suggest that environmental effects on multivariate genetic architecture may be comparable to the divergence that accumulates over dozens or hundreds of generations between populations. We outline avenues of future research to address the limitations of existing data and characterize the extent to which lability in genetic correlations shapes evolution in changing environments.     

Evolutionary potential of morphological traits across different life‐history stages

Despite accumulating examples of selection acting on heritable traits in the wild, predicted evolutionary responses are often different from observed phenotypic trends. Various explanations have been suggested for these mismatches. These include within‐individual changes across lifespan that can create important variation in genetic architecture of traits and selection acting on them, but also potential problems with the methodological approach used to predict evolutionary responses of traits. Here, we used an 8‐year data set on tree swallow (Tachycineta bicolor) to first assess the effects of differences among three nestling life‐history stages on the genetic (co)variances of two morphological traits (body mass and primary feather length) and the selection acting on them over three generations. We then estimated the evolutionary potential of these traits by predicting their evolutionary responses using the breeder's equation and the secondary theorem of selection approaches. Our results showed variation in strength and direction of selection and slight changes in trait variance across ages. Predicted evolutionary responses differed importantly between both approaches for half of the trait–age combinations we studied, suggesting the presence of environmentally induced correlations between focal traits and fitness possibly biasing breeder's equation predictions. Our results emphasize that predictions of evolutionary potential for morphological traits are likely to be highly variable, both in strength and direction, depending on the life stage and method used, thus mitigating our capacity to predict adaptation and persistence of wild populations.     

Rainfall can explain adaptive phenotypic variation with high gene flow in the New Holland Honeyeater (Phylidonyris novaehollandiae)

Identifying environmentally driven changes in traits that serve an ecological function is essential for predicting evolutionary outcomes of climate change. We examined population genetic structure, sex‐specific dispersal patterns, and morphology in relation to rainfall patterns across an island and three peninsulas in South Australia. The study system was the New Holland Honeyeater (Phylidonyris novaehollandiae), a nectarivorous passerine that is a key pollinator species. We predicted that rainfall‐related mechanisms would be driving local adaptation of morphological traits, such that in areas of lower rainfall, where nectar is less available, more insectivorous traits – shorter, deeper bills, longer tarsi, and longer wings – would be favored. The study populations differed in phenotype across the Eyre, Yorke, and Fleurieu Peninsulas and Kangaroo Island despite high gene flow (single continuous population) and sex‐biased dispersal (males were philopatric and females dispersed). We tested the role of rainfall in shaping the observed phenotypic differences, and found strong support for our predicted relationships: birds in areas of higher rainfall had higher condition indices, as well as longer bill‐head length, deeper bills, and shorter tarsi. Bill depth in males in high‐rainfall sites showed signals of stabilizing selection, suggesting local adaptation. In addition to these local indications of selection, a global pattern of directional selection toward larger size for bill‐head length, bill‐nostril length, and wing length was also observed. We suggest this pattern may reflect an adaptive response to the relatively dry conditions that South Australia has experienced over the last decade. We conclude that rainfall has shaped aspects of phenology in P. novaehollandiae, both locally, with different patterns of stabilizing and directional selection, and globally, with evidence of adaptive divergence at a landscape scale.     

Effects of migration on the genetic covariance matrix

In 1996, Schluter showed that the direction of morphological divergence of closely related species is biased toward the line of least genetic resistance, represented by g_max , the leading eigenvector of the matrix of genetic variance–covariance (the G -matrix). G is used to predict the direction of evolutionary change in natural populations. However, this usage requires that G is sufficiently constant over time to have enough predictive significance. Here, we explore the alternative explanation that G can evolve due to gene flow to conform to the direction of divergence between incipient species. We use computer simulations in a mainland–island migration model with stabilizing selection on two quantitative traits. We show that a high level of gene flow from a mainland population is required to significantly affect the orientation of the G -matrix in an island population. The changes caused by the introgression of the mainland alleles into the island population affect all aspects of the shape of G (size, eccentricity, and orientation) and lead to the alignment of g_max with the line of divergence between the two populations' phenotypic optima. Those changes decrease with increased correlation in mutational effects and with a correlated selection. Our results suggest that high migration rates, such as those often seen at the intraspecific level, will substantially affect the shape and orientation of G , whereas low migration (e.g., at the interspecific level) is unlikely to substantially affect the evolution of G .     

Dwarfism in insular sloths: biogeography,selection, and evolutionary rate

The islands of Bocas del Toro, Panama, were sequentially separated from the adjacent mainland by rising sea levels during the past 10,000 years. Three-toed sloths (Bradypus) from five islands are smaller than their mainland counterparts, and the insular populations themselves vary in mean body size. We first examine relationships between body size and physical characteristics of the islands, testing hypotheses regarding optimal body size, evolutionary equilibria, and the presence of dispersal in this system. To do so, we conduct linear regressions of body size onto island area, distance from the mainland, and island age. Second, we retroactively calculate two measures of the evolutionary rate of change in body size (haldanes and darwins) and the standardized linear selection differential, or selection intensity (i). We also test the observed morphological changes against models of evolution by genetic drift. The results indicate that mean body size decreases linearly with island age, explaining up to 97% of the variation among population means. Neither island area nor distance from the mainland is significant in multiple regressions that include island age. Thus, we find no evidence for differential optimal body size among islands, or for dispersal in the system. In contrast, the dependence of body size on island age suggests uniform directional selection for small body size in the insular populations. Although genetic drift cannot be discounted as the cause for this evolution in body size, the probability is small given the consistent direction of evolution (repeated dwarfism). The insular sloths show a sustained rate of evolution similar to those measured in haldanes over tens of generations, appearing to unite micro- and macroevolutionary time scales. Furthermore, the magnitude and rate of this example of rapid differentiation fall within predictions of theoretical models from population genetics. However, the linearity of the relationship between body size and island age is not predicted, suggesting that either more factors are involved than those considered here, or that theoretical advances are necessary to explain constant evolutionary rates over long time spans in new selective environments.     

Evolutionary rates for multivariate traits: the role of selection and genetic variation

A fundamental question in evolutionary biology is the relative importance of selection and genetic architecture in determining evolutionary rates. Adaptive evolution can be described by the multivariate breeders'' equation (), which predicts evolutionary change for a suite of phenotypic traits () as a product of directional selection acting on them (β) and the genetic variance–covariance matrix for those traits (G). Despite being empirically challenging to estimate, there are enough published estimates of G and β to allow for synthesis of general patterns across species. We use published estimates to test the hypotheses that there are systematic differences in the rate of evolution among trait types, and that these differences are, in part, due to genetic architecture. We find some evidence that sexually selected traits exhibit faster rates of evolution compared with life-history or morphological traits. This difference does not appear to be related to stronger selection on sexually selected traits. Using numerous proposed approaches to quantifying the shape, size and structure of G, we examine how these parameters relate to one another, and how they vary among taxonomic and trait groupings. Despite considerable variation, they do not explain the observed differences in evolutionary rates.     

Stability of the G-matrix in a population experiencing pleiotropic mutation, stabilizing selection, and genetic drift

Abstract. Quantitative genetics theory provides a framework that predicts the effects of selection on a phenotype consisting of a suite of complex traits. However, the ability of existing theory to reconstruct the history of selection or to predict the future trajectory of evolution depends upon the evolutionary dynamics of the genetic variance-covariance matrix (G-matrix). Thus, the central focus of the emerging field of comparative quantitative genetics is the evolution of the G-matrix. Existing analytical theory reveals little about the dynamics of G, because the problem is too complex to be mathematically tractable. As a first step toward a predictive theory of G-matrix evolution, our goal was to use stochastic computer models to investigate factors that might contribute to the stability of G over evolutionary time. We were concerned with the relatively simple case of two quantitative traits in a population experiencing stabilizing selection, pleiotropic mutation, and random genetic drift. Our results show that G-matrix stability is enhanced by strong correlational selection and large effective population size. In addition, the nature of mutations at pleiotropic loci can dramatically influence stability of G. In particular, when a mutation at a single locus simultaneously changes the value of the two traits (due to pleiotropy) and these effects are correlated, mutation can generate extreme stability of G. Thus, the central message of our study is that the empirical question regarding G-matrix stability is not necessarily a general question of whether G is stable across various taxonomic levels. Rather, we should expect the G-matrix to be extremely stable for some suites of characters and unstable for others over similar spans of evolutionary time.