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.     

QUANTITATIVE GENETICS OF MORPHOLOGICAL DIFFERENTIATION IN PEROMYSCUS. I. TESTS OF THE HOMOGENEITY OF GENETIC COVARIANCE STRUCTURE AMONG SPECIES AND SUBSPECIES

Phenotypic and additive genetic covariance matrices were estimated for 15 morphometric characters in three species and subspecies of Peromyscus. Univariate and multivariate ANOVAs indicate these groups are highly diverged in all characters, P. leucopus having the largest body size, P. maniculatus bairdii the smallest, and P. maniculatus nebrascensis being intermediate. Comparing the structure of P and G within each taxon revealed significant similarities in all three cases. This proportionality was strong enough to justify using P in the place of G to analyze evolutionary processes using quantitative genetic models when G can not be estimated, as in fossil material. However, the similarity between genetic and phenotypic covariance structures is sufficiently low that estimates of the genetic parameters should be used when possible. The additive genetic covariance matrices were compared to examine the assumption that they remain constant during evolution, an assumption which underlies many applications of quantitative-genetic models. While matrix permutation tests indicated statistically significant proportionality between the genetic covariance structures of the two P. maniculatus subspecies, there is no evidence of significant genetic structural similarity between species. This result suggests that the assumption of constant genetic covariance structure may be valid only within species. (It does not, however, necessarily imply a causal relationship between speciation and heterogeneity of genetic covariance structures.) The low matrix correlation for the two P. maniculatus subspecies' genetic covariance matrices indicates G may not be functionally constant, even within species. The lack of similarity observed here may be due partly to sampling variation.     

Effects of selection and drift on G matrix evolution in a heterogeneous environment: a multivariate Qst-Fst Test with the freshwater snail Galba truncatula

Unraveling the effect of selection vs. drift on the evolution of quantitative traits is commonly achieved by one of two methods. Either one contrasts population differentiation estimates for genetic markers and quantitative traits (the Q(st)-F(st) contrast) or multivariate methods are used to study the covariance between sets of traits. In particular, many studies have focused on the genetic variance-covariance matrix (the G matrix). However, both drift and selection can cause changes in G. To understand their joint effects, we recently combined the two methods into a single test (accompanying article by Martin et al.), which we apply here to a network of 16 natural populations of the freshwater snail Galba truncatula. Using this new neutrality test, extended to hierarchical population structures, we studied the multivariate equivalent of the Q(st)-F(st) contrast for several life-history traits of G. truncatula. We found strong evidence of selection acting on multivariate phenotypes. Selection was homogeneous among populations within each habitat and heterogeneous between habitats. We found that the G matrices were relatively stable within each habitat, with proportionality between the among-populations (D) and the within-populations (G) covariance matrices. The effect of habitat heterogeneity is to break this proportionality because of selection for habitat-dependent optima. Individual-based simulations mimicking our empirical system confirmed that these patterns are expected under the selective regime inferred. We show that homogenizing selection can mimic some effect of drift on the G matrix (G and D almost proportional), but that incorporating information from molecular markers (multivariate Q(st)-F(st)) allows disentangling the two effects.     

TYPE I ERROR RATES FOR TESTING GENETIC DRIFT WITH PHENOTYPIC COVARIANCE MATRICES: A SIMULATION STUDY

Studies of evolutionary divergence using quantitative genetic methods are centered on the additive genetic variance–covariance matrix ( G ) of correlated traits. However, estimating G properly requires large samples and complicated experimental designs. Multivariate tests for neutral evolution commonly replace average G by the pooled phenotypic within‐group variance–covariance matrix ( W ) for evolutionary inferences, but this approach has been criticized due to the lack of exact proportionality between genetic and phenotypic matrices. In this study, we examined the consequence, in terms of type I error rates, of replacing average G by W in a test of neutral evolution that measures the regression slope between among‐population variances and within‐population eigenvalues (the Ackermann and Cheverud [AC] test) using a simulation approach to generate random observations under genetic drift. Our results indicate that the type I error rates for the genetic drift test are acceptable when using W instead of average G when the matrix correlation between the ancestral G and P is higher than 0.6, the average character heritability is above 0.7, and the matrices share principal components. For less‐similar G and P matrices, the type I error rates would still be acceptable if the ratio between the number of generations since divergence and the effective population size (t/N_e) is smaller than 0.01 (large populations that diverged recently). When G is not known in real data, a simulation approach to estimate expected slopes for the AC test under genetic drift is discussed.     

QUANTITATIVE GENETICS OF BRYOZOAN PHENOTYPIC EVOLUTION. I. RATE TESTS FOR RANDOM CHANGE VERSUS SELECTION IN DIFFERENTIATION OF LIVING SPECIES

The possible roles of random genetic change and natural selection in bryozoan speciation were analyzed using quantitative genetic methods on breeding data for traits of skeletal morphology in two closely related species of the cheilostome Stylopoma. The hypothesis that morphologic differences between the species are caused entirely by mutation and genetic drift could not be rejected for reasonable rates of mutation maintained for as few as 10³ to 10⁴ generations. Divergence times this short or shorter are consistent with the abrupt appearances of many invertebrate species in the fossil record, commonly followed by millions of years of morphologic stasis. To produce these differences over 10³ generations or fewer, directional selection acting alone would require unrealistically high levels of minimum selective mortality throughout divergence. Thus, selection is unnecessary to explain the divergence of these species, except as a means of accelerating the effects of random genetic change on shorter time scales (directional selection), or decelerating them over longer ones (stabilizing selection). These results are consistent with a variety of models of phenotypic evolution involving random shifts between multiple adaptive peaks. Similar results were obtained by substituting trait heritabilities and genetic covariances reconstructed by partitioning within- and among-colony phenotypic variance in place of the values based on breeding data. Quantitative genetic analysis of speciation in fossil bryozoan lineages is thus justified.     

BRINGING HABITAT INFORMATION INTO STATISTICAL TESTS OF LOCAL ADAPTATION IN QUANTITATIVE TRAITS: A CASE STUDY OF NINE‐SPINED STICKLEBACKS

Detection of footprints of historical natural selection on quantitative traits in cross‐sectional data sets is challenging, especially when the number of populations to be compared is small and the populations are subject to strong random genetic drift. We extend a recent Bayesian multivariate approach to differentiate between selective and neutral causes of population differentiation by the inclusion of habitat information. The extended framework allows one to test for signals of selection in two ways: by comparing the patterns of population differentiation in quantitative traits and in neutral loci, and by comparing the similarity of habitats and phenotypes. We illustrate the framework using data on variation of eight morphological and behavioral traits among four populations of nine‐spined sticklebacks (Pungitius pungitius). In spite of the strong signal of genetic drift in the study system (average F_ST = 0.35 in neutral markers), strong footprints of adaptive population differentiation were uncovered both in morphological and behavioral traits. The results give quantitative support for earlier qualitative assessments, which have attributed the observed differentiation to adaptive divergence in response to differing ecological conditions in pond and marine habitats.     

Population divergence and morphometric integration in the greenfinch (Carduelis chloris) – evolution against the trajectory of least resistance?

Theory predicts that genetic and phenotypic correlations among traits may direct the process of short-term evolution by limiting the directions of variation available to natural selection to act on. We studied correlations between 14 skeletal traits in 10 geographically distinct and relatively young greenfinch (Carduelis chloris) populations to unravel whether the divergence among populations has occurred into directions predicted by the within-population correlations (cf. drift/correlated responses models), or whether it is better explained by ‘adaptive’ models, which predict no necessary association between within- and among-population correlations (allometries). We found that the within-population character correlations (or covariances) did not predict character divergence between populations. This was because the first eigenvector of the among-population correlation/covariance matrix, summarizing the major dimension of divergence, was a bipolar body:beak dimension, and distinct from the (≈ isometric) first eigenvector of within-population matrix. Hence, as the divergence among greenfinch populations cannot be satisfactorily accommodated by drift/correlated response models, an adaptive basis for divergence is suggested. The second major axis of within-population variation was a classical ‘group size’ factor revealing that beak size was more or less free to vary independently of body size. Consequently, even if the divergence among populations cannot be simply accommodated to expectations of drift and correlated response models, it is striking that the most pronounced size-independent (nonallometric) changes had occurred along the second largest dimension of variance. This could mean that selection pressures which shape integration within populations are the same as those that cause divergence among populations. A relaxed beak:body integration could also occur as a result of species level selection favouring taxa in which independent evolution of beak and body is made possible.     

Microevolution in island forms: the roles of drift and directional selection in morphological divergence of a passerine bird

Abstract.— Theory predicts that in small isolated populations random genetic drift can lead to phenotypic divergence; however this prediction has rarely been tested quantitatively in natural populations. Here we utilize natural repeated island colonization events by members of the avian species complex, Zosterops lateralis , to assess whether or not genetic drift alone is an adequate explanation for the observed patterns of microevolutionary divergence in morphology. Morphological and molecular genetic characteristics of island and mainland populations are compared to test three predictions of drift theory: (1) that the pattern of morphological change is idiosyncratic to each island; (2) that there is concordance between morphological and neutral genetic shifts across island populations; and (3) for populations whose time of colonization is known, that the rate of morphological change is sufficiently slow to be accounted for solely by genetic drift. Our results are not consistent with these predictions. First, the direction of size shifts was consistently towards larger size, suggesting the action of a nonrandom process. Second, patterns of morphological divergence among recently colonized populations showed little concordance with divergence in neutral genetic characters. Third, rate tests of morphological change showed that effective population sizes were not small enough for random processes alone to account for the magnitude of microevolutionary change. Altogether, these three lines of evidence suggest that drift alone is not an adequate explanation of morphological differentiation in recently colonized island Zosterops and therefore we suggest that the observed microevolutionary changes are largely a result of directional natural selection.     

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.     

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.     

A genetic model of interpopulation variation and covariation of quantitative characters

Evolutionary consequences of natural selection, migration, genotype-environment interaction, and random genetic drift on interpopulation variation and covariation of quantitative characters are analysed in terms of a selection model that partitions natural selection into directional and stabilizing components. Without migration, interpopulation variation and covariation depend mainly on the pattern and intensities of selection among populations and the harmonic mean of effective population sizes. Both transient and equilibrium covariance structures are formulated with suitable approximations. Migration reduces the differentiation among populations, but its effect is less with genotype-environment interaction. In some special cases of genotype-environment interaction, the equilibrium interpopulation variation and covariation is independent of migration.     

QUANTITATIVE GENETICS OF BRYOZOAN PHENOTYPIC EVOLUTION. II. ANALYSIS OF SELECTION AND RANDOM CHANGE IN FOSSIL SPECIES USING RECONSTRUCTED GENETIC PARAMETERS

The roles of natural selection and random genetic change in the punctuated phenotypic evolution of eight Miocene-Pliocene tropical American species of the cheilostome bryozoan Metrarabdotos are analyzed by quantitative genetic methods. Trait heritabilities and genetic covariances reconstructed by partitioning within- and among-colony phenotypic variance are similar to those previously obtained for living species of the cheilostome Stylopoma using breeding data. The hypothesis that differences in skeletal morphology between species of Metrarabdotos are entirely due to mutation and genetic drift cannot be rejected for reasonable rates of mutation maintained for periods brief enough to account for the geologically abrupt appearances of these species in the fossil record. Except for one pair of species, separated by the largest morphologic distance, directional selection acting alone would require unrealistically high rates of selective mortality to be maintained for these periods. Thus, directional selection is not strongly implicated in the divergence of Metrarabdotos species. Within species, rates of net phenotypic change are slow enough to require stabilizing selection, but mask large, relatively rapid fluctuations, all of which, however, can be attributed to chance departures from the mean phenotype by mutation and genetic drift, rather than to tracking environmental fluctuation by directional selection. The results are consistent with genetic models involving shifts between multiple adaptive peaks on which phenotypes remain more or less static through long-term stabilizing selection. Regardless of the degree to which directional selection may be involved in peak shifts, phenotypic differentiation is thus related to processes different than the pervasive stabilizing selection acting within species.     

Diversifying selection drives parallel evolution of gill raker number and body size along the speciation continuum of European whitefish

Adaptive radiation is the evolution of ecological and phenotypical diversity. It arises via ecological opportunity that promotes the exploration of underutilized or novel niches mediating specialization and reproductive isolation. The assumed precondition for rapid local adaptation is diversifying natural selection, but random genetic drift could also be a major driver of this process. We used 27 populations of European whitefish (Coregonus lavaretus) from nine lakes distributed in three neighboring subarctic watercourses in northern Fennoscandia as a model to test the importance of random drift versus diversifying natural selection for parallel evolution of adaptive phenotypic traits. We contrasted variation for two key adaptive phenotypic traits correlated with resource utilization of polymorphic fish; the number of gill rakers and the total length of fish, with the posterior distribution of neutral genetic differentiation from 13 microsatellite loci, to test whether the observed phenotypic divergence could be achieved by random genetic drift alone. Our results show that both traits have been under diversifying selection and that the evolution of these morphs has been driven by isolation through habitat adaptations. We conclude that diversifying selection acting on gill raker number and body size has played a significant role in the ongoing adaptive radiation of European whitefish morphs in this region.     

Contrasting patterns of body shape and neutral genetic divergence in marine and lake populations of threespine sticklebacks

Comparisons of neutral marker and quantitative trait divergence can provide important insights into the relative roles of natural selection and neutral genetic drift in population differentiation. We investigated phenotypic and genetic differentiation among Fennoscandian threespine stickleback (Gasterosteus aculeatus) populations, and found that the highest degree of differentiation occurred between sea and freshwater habitats. Within habitats, morphological divergence was highest among the different freshwater populations. Pairwise phenotypic and neutral genetic distances among populations were positively correlated, suggesting that genetic drift may have contributed to the morphological differentiation among habitats. On the other hand, the degree of phenotypic differentiation (PST) clearly surpassed the neutral expectation set by FST, suggesting a predominant role for natural selection over genetic drift as an explanation for the observed differentiation. However, separate PST/FST comparisons by habitats revealed that body shape divergence between lake and marine populations, and even among marine populations, can be strongly influenced by natural selection. On the other hand, genetic drift can play an important role in the differentiation among lake populations.     

High evolutionary constraints limited adaptive responses to past climate changes in toad skulls

Interactions among traits that build a complex structure may be represented as genetic covariation and correlation. Genetic correlations may act as constraints, deflecting the evolutionary response from the direction of natural selection. We investigated the relative importance of drift, selection, and constraints in driving skull divergence in a group of related toad species. The distributional range of these species encompasses very distinct habitats with important climatic differences and the species are primarily distinguished by differences in their skulls. Some parts of the toad skull, such as the snout, may have functional relevance in reproductive ecology, detecting water cues. Thus, we hypothesized that the species skull divergence was driven by natural selection associated with climatic variation. However, given that all species present high correlations among skull traits, our second prediction was of high constraints deflecting the response to selection. We first extracted the main morphological direction that is expected to be subjected to selection by using within- and between-species covariance matrices. We then used evolutionary regressions to investigate whether divergence along this direction is explained by climatic variation between species. We also used quantitative genetics models to test for a role of random drift versus natural selection in skull divergence and to reconstruct selection gradients along species phylogeny. Climatic variables explained high proportions of between-species variation in the most selected axis. However, most evolutionary responses were not in the direction of selection, but aligned with the direction of allometric size, the dimension of highest phenotypic variance in the ancestral population. We conclude that toad species have responded to selection related to climate in their skulls, yet high evolutionary constraints dominated species divergence and may limit species responses to future climate change.     

Rapid divergence of genetic variance-covariance matrix within a natural population

The matrix of genetic variances and covariances (G matrix) represents the genetic architecture of multiple traits sharing developmental and genetic processes and is central for predicting phenotypic evolution. These predictions require that the G matrix be stable. Yet the timescale and conditions promoting G matrix stability in natural populations remain unclear. We studied stability of the G matrix in a 20-year evolution field experiment, where a population of the cosmopolitan parthenogenetic soil nematode Acrobeloides nanus was subjected to drift and divergent selection (benign and stress environments). Selection regime did not influence the level of absolute genetic constraints: under both regimes, two genetic dimensions for three life-history traits were identified. A substantial response to selection in principal components structure and in general matrix pattern was indicated by three statistical methods. G structure was also influenced by drift, with higher divergence under benign conditions. These results show that the G matrix might evolve rapidly in natural populations. The observed high dynamics of G structure probably represents the general feature of asexual species and limits the predictive power of G in phenotypic evolution analyses.     

Inbreeding changes the shape of the genetic covariance matrix in Drosophila melanogaster.

The pattern of genetic covariation among traits (the G matrix) plays a central role in determining the pattern of evolutionary change from both natural selection and random genetic drift. Here we measure the effect of genetic drift on the shape of the G matrix using a large data set on the inheritance of wing characteristics in Drosophila melanogaster. Fifty-two inbred lines with a total of 4680 parent-offspring families were generated by one generation of brother-sister mating and compared to an outbred control population of 1945 families. In keeping with the theoretical expectation for a correlated set of additively determined traits, the average G matrix of the inbred lines remained proportional to the outbred control G matrix with a proportionality constant approximately equal to (1 - F), where F is the inbreeding coefficient. Further, the pattern of covariance among the means of the inbred lines induced by inbreeding was also proportional to the within-line G matrix of the control population with a constant very close to the expectation of 2F. Although the average G of the inbred lines did not show change in overall structure relative to the outbred controls, separate analysis revealed a great deal of variation among inbred lines around this expectation, including changes in the sign of genetic correlations. Since any given line can be quite different from the outbred control, it is likely that in nature unreplicated drift will lead to changes in the G matrix. Thus, the shape of G is malleable under genetic drift, and the evolutionary response of any particular population is likely to depend on the specifics of its evolutionary history.     

Statistical methods for detecting natural selection from genomic data

In the study of molecular and phenotypic evolution, understanding the relative importance of random genetic drift and positive selection as the mechanisms for driving divergences between populations and maintaining polymorphisms within populations has been a central issue. A variety of statistical methods has been developed for detecting natural selection operating at the amino acid and nucleotide sequence levels. These methods may be largely classified into those aimed at detecting recurrent and/or recent/ongoing natural selection by utilizing the divergence and/or polymorphism data. Using these methods, pervasive positive selection has been identified for protein-coding and non-coding sequences in the genomic analysis of some organisms. However, many of these methods have been criticized by using computer simulation and real data analysis to produce excessive false-positives and to be sensitive to various disturbing factors. Importantly, some of these methods have been invalidated experimentally. These facts indicate that many of the statistical methods for detecting natural selection are unreliable. In addition, the signals that have been believed as the evidence for fixations of advantageous mutations due to positive selection may also be interpreted as the evidence for fixations of deleterious mutations due to random genetic drift. The genomic diversity data are rapidly accumulating in various organisms, and detection of natural selection may play a critical role for clarifying the relative role of random genetic drift and positive selection in molecular and phenotypic evolution. It is therefore important to develop reliable statistical methods that are unbiased as well as robust against various disturbing factors, for inferring natural selection.     

Does natural selection alter genetic architecture? An evaluation of quantitative genetic variation among populations of Allonemobiussocius and A. fasciatus

To make long-term predictions using present quantitative genetic theory it is necessary to assume that the genetic variance–covariance matrix ( G ) remains constant or at least changes by a constant fraction. In this paper we examine the stability of the genetic architecture of two traits known to be subject to natural selection; femur length and ovipositor length in two species of the cricket Allonemobius. Previous studies have shown that in A. fasciatus and A. socius natural selection favours an increased body size southwards but a decreased ovipositor length. Such countergradient selection should tend to favour a change in G . In the total sample of eight populations of A. socius and one of A. fasciatus we show that there is significant variation in all genetic covariance components, i.e. V_A for body size, V_A for ovipositor length, and Cov_A. This variation results entirely from an increase in the covariances of A. fasciatus. However, although larger, these components are approximately proportionally increased, thereby leading to no statistically significant change in the genetic correlation. A proportional increase in the covariance components is consistent with changes resulting from genetic drift. On the other hand, the genetic covariance components are significantly correlated with the length of the growing season suggesting that the change in the genetic architecture is the result of selection and drift.     

Cranial phenotypic evolution in Proechimys iheringi Thomas (Rodentia: Echimyidae)

Models of evolutionary quantitative genetics were employed to analyze cranial phenotypic evolution in the echimyid rodent Proechimys iheringi Thomas. The hypothesis that cranial differentiation among populations and subspecies of P. iheringi resulted from random drift was rejected. The completely selective model was then used to reconstruct the forces of selection that might have acted to produce the observed differences in cranial morphometric traits in P. iheringi . The net gradient of selection indicates that some cranial traits evolved in the direction contrary to the forces of selection, acting to increase mean values for cranial traits in the evolutionary transitions between populations of P. iheringi . Minimum intensities of selection were calculated and showed that selective mortality per generation was low, of the order of l0^-3, suggesting that weak selection can explain the observed cranial morphometric differences among the populations of P. iheringi .