It has been claimed that hundreds of researchers use nested clade phylogeographic analysis (NCPA) based on what the method promises rather than requiring objective validation of the method. The supposed failure of NCPA is based upon the argument that validating it by using positive controls ignored type I error, and that computer simulations have shown a high type I error. The first argument is factually incorrect: the previously published validation analysis fully accounted for both type I and type II errors. The simulations that indicate a 75% type I error rate have serious flaws and only evaluate outdated versions of NCPA. These outdated type I error rates fall precipitously when the 2003 version of single-locus NCPA is used or when the 2002 multilocus version of NCPA is used. It is shown that the tree-wise type I errors in single-locus NCPA can be corrected to the desired nominal level by a simple statistical procedure, and that multilocus NCPA reconstructs a simulated scenario used to discredit NCPA with 100% accuracy. Hence, NCPA is a not a failed method at all, but rather has been validated both by actual data and by simulated data in a manner that satisfies the published criteria given by its critics. The critics have come to different conclusions because they have focused on the pre-2002 versions of NCPA and have failed to take into account the extensive developments in NCPA since 2002. Hence, researchers can choose to use NCPA based upon objective critical validation that shows that NCPA delivers what it promises.

Evolutionary quantitative genetics has recently advanced in two distinct streams. Many biologists address evolutionary questions by estimating phenotypic selection and genetic (co)variances (G matrices). Simultaneously, an increasing number of studies have applied quantitative trait locus (QTL) mapping methods to dissect variation. Both conceptual and practical difficulties have isolated these two foci of quantitative genetics. A conceptual integration follows from the recognition that QTL allele frequencies are the essential variables relating the G-Matrix to marker-based mapping experiments. Breeding designs initiated from randomly selected parental genotypes can be used to estimate QTL-specific genetic (co)variances. These statistics appropriately distill allelic variation and provide an explicit population context for QTL mapping estimates. Within this framework, one can parse the G-Matrix into a set of mutually exclusive genomic components and ask whether these parts are similar or dissimilar in their respective features, for example the magnitude of phenotypic effects and the extent and nature of pleiotropy. As these features are critical determinants of sustained response to selection, the integration of QTL mapping methods into G-Matrix estimation can provide a concrete, genetically based experimental program to investigate the evolutionary potential of natural populations.

The mean fitness of a population, often equal to its growth rate, measures its level of adaptation to particular environmental conditions. A better understanding of the evolution of mean fitness could thus provide a natural link between evolution and demography. Yet, after the seminal work of Fisher and its renowned “fundamental theorem of natural selection,” the dynamics of mean fitness has attracted little attention, and mostly from theoretical population geneticists. Here we analyze the dynamics of mean fitness in the context of host-parasite interactions. We illustrate the potential relevance of this analysis under different scenarios ranging from a simple situation in which a parasite evolves in a homogeneous host population to a more complex one with host-parasite coevolution. In each case, we contrast the effects of natural selection, recurrent mutations, and the change of the biotic environment, on the dynamics of adaptation. Decoupling these three components helps elucidate the interplay between evolutionary and ecological dynamics. In particular, it offers new perspectives on situations leading to evolutionary suicide. As mean fitness is an easily measurable quantity in microbial systems, this analysis provides new ways to track the dynamics of adaptation in experimental evolution and coevolution studies.

Mutualisms often involve reciprocal adaptations of both partners. Acacia ant-plants defended by symbiotic Pseudomyrmex ant mutualists secrete sucrose-free extrafloral nectar, which is unattractive to generalists. We aimed to investigate whether this extrafloral nectar can also exclude exploiters, that is nondefending ant species. Mutualist workers discriminated against sucrose whereas exploiters and generalists with no affinity toward Acacia myrmecophytes preferred sucrose, because mutualist workers lacked the sucrose-cleaving enzyme invertase, which is present in workers of the other two groups. Sucrose uptake induced invertase activity in workers of parasites and generalists, but not mutualists, and in larvae of all species: the mutualists loose invertase during their ontogeny. This reduced metabolic capacity ties the mutualists to their plant hosts, but it does not completely prevent the mutualism from exploitation. We therefore investigated whether the exploiters studied here are cheaters (i.e., have evolved from former mutualists) or parasites (exploiters with no mutualistic ancestor). A molecular phytogeny demonstrates that the exploiter species did not evolve from former mutualists, and no evidence for cheaters was found. We conclude that being specialized to their partner can prevent mutualists from becoming cheaters, whereas other mechanisms are required to stabilize a mutualism against the exploitation by parasites.

Diversification of insect herbivores is often associated with coevolution between plant toxins and insect countermeasures, resulting in a specificity that restricts host plant shifts. Gall inducers, however, bypass plant toxins and the factors influencing host plant associations in these specialized herbivores remain unclear. We reconstructed the evolution of host plant associations in Western Palaearctic oak gallwasps (Cynipidae: Cynipini), a species-rich lineage of specialist herbivores on oak (Quercus). (1) Bayesian analyses of sequence data for three genes revealed extreme host plant conservatism, with inferred shifts between major oak lineages (sections Cerris and Quercus) closely matching the minimum required to explain observed diversity. It thus appears that the coevolutionary demands of gall induction constrain host plant shifts, both in cases of mutualism (e.g., fig wasps, yucca moths) and parasitism (oak gallwasps). (2) Shifts between oak sections occurred independently in sexual and asexual generations of the gallwasp lifecycle, implying that these can evolve independently. (3) Western Palaearctic gallwasps associated with sections Cerris and Quercus diverged at least 20 million years ago (mya), prior to the arrival of oaks in the Western Palaearctic from Asia 5–7 mya. This implies an Asian origin for Western Palaearctic gallwasps, with independent westwards range expansion by multiple lineages.

The wild progenitors of cultivated rice, Oryza nivara and Oryza rufipogon, provide an experimental system for characterizing the genetic basis of adaptation. The evolution of annual O. nivara from a perennial ancestor resembling its sister species, O. rufipogon, was associated with an ecological shift from persistently wet to seasonally dry habitats. Here we report a quantitative trait locus (QTL) analysis of phenotypic differentiation in life history, mating system, and flowering time between O. nivara and O. rufipogon. The exponential distribution of effect sizes of QTL fits the prediction of a recently proposed population genetic model of adaptation. More than 80% of QTL alleles of O. nivara acted in the same direction of phenotypic evolution, suggesting that they were fixed under directional selection. The loss of photoperiod sensitivity, which might be essential to the survival of the ancestral populations of O. nivara in the new environment, was controlled by QTL of relatively large effect. Mating system evolution from cross- to self-fertilization through the modification of panicle and floral morphology was controlled by QTL of small-to-moderate effect. The lack of segregation of the recessive annual habit in the F₂ mapping populations suggested that the evolution of annual from perennial life form had a complex genetic basis. The study captured the genetic architecture for the adaptive origin of O. nivara and provides a foundation for rigorous experimental tests of population genetic theories of adaptation.

Interspecific hybridization can generate transgressive hybrid phenotypes with extreme trait values exceeding the combined range of the parental species. Such variation can enlarge the working surface for natural selection, and may facilitate the evolution of novel adaptations where ecological opportunity exists. The number of quantitative trait loci fixed for different alleles in different species should increase with time since speciation. If transgression is caused by complementary gene action or epistasis, hybrids between more distant species should be more likely to display transgressive phenotypes. To test this prediction we collected data on transgression frequency from the literature, estimated genetic distances between the hybridizing species from gene sequences, and calculated the relationship between the two using phylogenetically controlled methods. We also tested if parental phenotypic divergence affected the occurrence of transgression. We found a highly significant positive correlation between transgression frequency and genetic distance in eudicot plants explaining 43% of the variance in transgression frequency. In total, 36% of the measured traits were transgressive. The predicted effect of time since speciation on transgressive segregation was unconfounded by the potentially conflicting effects of phenotypic differentiation between species. Our analysis demonstrates that the potential impact hybridization may have on phenotypic evolution is predictable from the genetic distance between species.

Despite long-standing interest in character displacement, little is known of its underlying proximate causes. Here, we explore the role of maternal effects in character displacement. We specifically investigated whether differences in maternal body condition mediate divergence in the expression of resource-use traits between populations of spadefoot toads (Spea multiplicata) that occur in sympatry with a heterospecific competitor and those that occur in allopatry. In sympatry, S. multiplicata is forced by its competitor onto a less profitable resource. As a result, sympatric females mature in poorer condition and invest less into offspring. Consequently, their offspring produce a resource-use phenotype that minimizes competition with the other species and that also differs from the phenotype produced in allopatry. These differences in trait expression between allopatry and sympatry disappear once mothers are equilibrated in body condition in the laboratory. Thus, a condition-dependent maternal effect mediates population divergence and character displacement. Such effects potentially buffer populations from extinction (via competitive exclusion) while genetic changes accumulate, which produce divergent traits in the absence of the maternal effect. Maternal effects may therefore often be important in determining the initial direction and rate of evolution during the early stages of character displacement.

The origin of sexual size dimorphisms (SSD) has long been a central topic in evolutionary biology. However, there is little agreement as to which factors are most important in driving the evolution of SSD, and several hypotheses concerning SSD evolution have never been tested empirically. Emydid turtles include species with both male and female-biased SSD, and some emydids exhibit among the most extreme SSD in tetrapods. Here, we use a comparative phylogenetic approach in emydids to analyze the origins of SSD and test several hypotheses for the evolution of SSD, some for the first time. We test the Fairbairn-Preziosi hypothesis for the origin of Rensch's rule, and support it in lineages with male-biased SSD but not those with female-biased SSD. We also find support for the secondary ecological dimorphism hypothesis, which proposes that selection for ecological divergence between sexes exaggerates preexisting SSD. Finally, we find only equivocal support for the Bolnick-Doebeli hypothesis, which relates intersexual ecological divergence to interspecific ecological divergence. Our results also illustrate how global analyses of SSD may mislead in groups in which the factors that drive the evolution of SSD vary among clades.

Sexual selection theory makes clear predictions regarding male spermatogenic investment. To test these predictions we used experimental sexual selection in Drosophila pseudoobscura, a sperm heteromorphic species in which males produce both fertile and sterile sperm, the latter of which may function in postmating competition. Specifically, we determined whether the number and size of both sperm types, as well as relative testis mass and accessory gland size, increased with increased sperm competition risk and whether any fitness benefits could accrue from such changes. We found no effect of sexual selection history on either the number or size of either sperm morph, or on relative testis mass. However, males experiencing a greater opportunity for sexual selection evolved the largest accessory glands, had the greatest mating capacity, and sired the most progeny. These findings suggest that sterile sperm are not direct targets of sexual selection and that accessory gland size, rather than testis mass, appears to be an important determinant of male reproductive success. We briefly review the data from experimental sexual selection studies and find that testis mass may not be a frequent target of postcopulatory sexual selection and, even when it is, the resulting changes do not always improve fitness.

Numerous theoretical studies have investigated how limited dispersal may provide an explanation for the evolution of cooperation, by leading to interactions between relatives. However, despite considerable theoretical attention, there has been a lack of empirical tests. In this article, we test how patterns of dispersal influence the evolution of cooperation, using iron-scavenging in the bacterium Pseudomonas aeruginosa as our cooperative trait. We found that relatively limited dispersal does not favor cooperation. The reason for this is that although limited dispersal increases the relatedness between interacting individuals, it also leads to increased local competition for resources between relatives. This result supports Taylor's prediction that in the simplest possible scenario, the effects of increased relatedness and local competition exactly cancel out. In contrast, we show that one way for cooperation to be favored is if individuals disperse in groups (budding dispersal), because this maintains high relatedness while reducing local competition between relatives (relatively global competition).

Clonal interference refers to the competition that arises in asexual populations when multiple beneficial mutations segregate simultaneously. A large body of theoretical and experimental work now addresses this issue. Although much of the experimental work is performed in populations that grow exponentially between periodic population bottlenecks, the theoretical work to date has addressed only populations of a constant size. We derive an analytical approximation for the rate of adaptation in the presence of both clonal interference and bottlenecks, and compare this prediction to the results of an individual-based simulation, showing excellent agreement in the parameter regime in which clonal interference prevails. We also derive an appropriate definition for the effective population size for adaptive evolution experiments in the presence of population bottlenecks. This ‘adaptation effective population size’ allows for a good approximation of the expected rate of adaptation, either in the strong-selection weak-mutation regime, or when clonal interference comes into play. In the multiple mutation regime, when the product of the population size and mutation rate is extremely large, these results no longer hold.

Substantial insular speciation has resulted in exceptionally high levels of endemism in Madagascar, creating locally restricted species' ranges that remain poorly understood. The contributions of alternative processes that could influence patterns of local endemism—including speciation by geographic isolation or adaptation to environmental gradients—are widely debated, both for Madagascar and elsewhere. A recently proposed hypothesis (the “watershed hypothesis”) suggests that allopatric speciation driven by isolation in watersheds during Quaternary climate shifts provides a general explanation for patterns of local endemism across taxa in Madagascar. Here we tested coincidence between species' distributions and areas of endemism predicted by two contrasting biogeographic hypotheses: (1) the watershed hypothesis, and (2) an alternative hypothesis driven by climatic gradients (the “current climate hypothesis”). Statistical significance of coincidence was assessed by comparing against a null model. Surprisingly, we found that extant distributions of lemurs, geckos, and chameleons reveal species patterns that are significantly coincident with the watershed and current climate hypotheses. These results strongly support local endemism developing from multiple processes, even among closely related species. Our findings thus indicate that pluralistic approaches will offer the best option both for understanding processes that generate local endemism, and for incorporating endemism within conservation priority setting.

Discrete behavioral strategies comprise a suite of traits closely integrated in their expression with consistent natural selection for such coexpression leading to developmental and genetic integration of their components. However, behavioral traits are often also selected to respond rapidly to changing environments, which should both favor their context-dependent expression and inhibit evolution of genetic integration with other, less flexible traits. Here we use a multigeneration pedigree and longterm data on lifetime fitness to test whether behaviors comprising distinct dispersal strategies of western bluebirds—a species in which the propensity to disperse is functionally integrated with aggressive behavior—are genetically correlated. We further investigated whether selection favors flexibility in the expression of aggression in relation to current social context. We found a significant genetic correlation between aggression and dispersal that is concordant with consistent selection for coexpression of these behaviors. To a limited extent, individuals modified their aggression to match their mate; however, we found no fitness consequences on such adjustments. These results introduce a novel way of viewing behavioral strategies, where flexibility of behavior, while often aiding an organism's fit in its current environment, may be limited and thereby enable integration with less flexible traits.

The causes of craniofacial variation among human populations have been the subject of controversy. In this work, we studied aboriginal populations from southern South America, the last continental region peopled by humans and with a wide range of ecological conditions. Because of these characteristics, southern South America provides a unique opportunity to study the relative importance of random and nonrandom factors in human diversification. Previous craniometric studies recognized remarkable differences among populations from this region, usually resorting to random factors as the main explanation. In contrast, here we suggest, using tests based on quantitative genetic models, that: (1) the rate of craniofacial divergence among these populations is too high and (2) the patterns of variation within and between populations are too different to be explained by genetic drift alone. In addition, the among-sample craniofacial variation is correlated with climate and diet but not with mtDNA variation. We suggest that the influence of nonrandom factors (e.g., plasticity, selection) on human craniofacial diversification in regions with large ecological variation is more important than generally acknowledged and capable to generate large craniofacial divergence in a short period of time. These results bring nonrandom factors into focus for the interpretation of human craniofacial variation.

Ontogenetic and behavioral studies using birds currently do not document the early evolution of flight because birds (including juveniles) used in such studies employ forelimb oscillation frequencies over 10 Hz, forelimb stroke-angles in excess of 130°, and possess uniquely avian flight musculatures. Living birds are an advanced morphological stage in the development of flapping flight. To gain insight into the early stages of flight evolution (i.e., prebird), in the absence of a living analogue, a new approach using Strouhal number (St = f f · A/U) was used. Strouhal number is a nondimensional number that describes the relationship between wing-stroke amplitude (A), wing-beat frequency (f), and flight speed (U). Calculations indicated that even moderate wing movements are enough to generate rudimentary thrust and that a propulsive flapping flight-stroke could have evolved via gradual incremental changes in wing movement and wing morphology. More fundamental to the origin of the avian flapping flight-stroke is the question of how a symmetrical forelimb posture—required for gliding and flapping flight—evolved from an alternating forelimb motion, evident in all extant bipeds when running except birds.

Quantifying morphological shape is a fundamental issue in evolutionary biology. Recent technological advances (e.g., confocal microscopy, laser scanning, computer tomography) have made the capture of detailed three-dimensional (3D) morphological structure easy and cost-effective. In this article, we develop a 3D analytic framework (SPHARM—spherical harmonics) for modeling the shapes of complex morphological structures from continuous surface maps that can be produced by these technologies. Because the traditional SPHARM methodology has limitations in several of its processing steps, we present new algorithms for two SPHARM processing steps: spherical parameterization and SPHARM registration. These new algorithms allow for the numerical characterization of a much larger class of 3D models. We demonstrate the effectiveness of the method by applying it to modeling the cerci of Enallagma damselflies.

G.C. Williams's 1957 hypothesis famously argues that higher age-independent, or ‘extrinsic,’ mortality should select for faster rates of senescence. Long-lived species should therefore show relatively few deaths from extrinsic causes such as predation and starvation. Theoretical explorations and empirical tests of Williams's hypothesis have flourished in the past decade but it has not yet been tested empirically among humans. We test Williams's hypothesis using mortality data from subsistence populations and from historical cohorts from Sweden and England/Wales, and examine whether rates of actuarial aging declined over the past two centuries. We employ three aging measures: mortality rate doubling time (MRDT), Ricklefs's ω;, and the slope of mortality hazard from ages 60–70, m′_50–70, and model mortality using both Weibull and Gompertz-Makeham hazard models. We find that (1) actuarial aging in subsistence societies is similar to that of early Europe, (2) actuarial senescence has slowed in later European cohorts, (3) reductions in extrinsic mortality associate with slower actuarial aging in longitudinal samples, and (4) men senesce more rapidly than women, especially in later cohorts. To interpret these results, we attempt to bridge population-based evolutionary analysis with individual-level proximate mechanisms.

If immune functions confer obvious benefits to hosts, life-history theory assumes that they also induce costs, leading to trade-offs between immunity and other fitness components. However, whether substantial fitness costs are associated with immune systems in the wild is debatable, as numerous factors may influence the costs and benefits associated with immune activation. Here, we explore the survival cost of immune deployment in postfledging birds. We injected Eurasian collared dove nestlings (Streptopelia decaocto) with antigens from Escherichia coll, and examined whether this immune challenge affected survival after fledging. To assess survival, birds were fitted with radlotags and the fate of each individual was monitored regularly. Our results show that mimicking a bacterial infection in nestlings lowered their survival prospects after fledging, in comparison to controls. The main identified cause of mortality (by examination of dead birds) was presumed to be predation. This study provides experimental evidence that immune activation may entail dramatic survival costs in a free-ranging vertebrate, and emphasizes the potential role that environmental factors such as predation may play in this interaction.

Phenotypic selection is widely recognized as the primary cause of adaptive evolution in natural populations, a fact that has been documented frequently over the last few decades, mainly in morphological and life-history traits. The energetic definition of fitness predicts that natural selection will maximize the residual energy available for growth and reproduction, suggesting that energy metabolism could be a target of selection. To address this problem, we chose the garden snail. Helix asperse (Cornu aspersum). We performed a seminatural experiment for measuring phenotypic selection on standard metabolic rate (SMR), the minimum cost of maintenance in ectotherm organisms. To discount selection on correlated traits, we included two additional whole-organism performance traits (mean speed and maximum force of dislodgement). We found a combination of linear (negative directional selection, β = -0.106 ± 0.06; P = 0.001) and quadratic (stabilizing selection, γ = -0.012 ± 0.033; P = 0.061) selection on SMR. Correlational selection was not significant for any possible pair of traits. This suggests that individuals with average-to-reduced SMRs were promoted by selection. To the best of our knowledge, this is the first study showing significant directional selection on the obligatory cost of maintenance in an animal, providing support for the energetic definition of fitness.

Estimating quantitative genetic parameters ideally takes place in natural populations, but relatively few studies have overcome the inherent logistical difficulties. For this reason, no estimates currently exist for the genetic basis of life-history traits in natural populations of large marine vertebrates. And yet such estimates are likely to be important given the exposure of this taxon to changing selection pressures, and the relevance of life-history traits to population productivity. We report such estimates from a long-term (1995–2007) study of lemon sharks (Negaprion brevirostris) conducted at Bimini, Bahamas. We obtained these estimates by genetically reconstructing a population pedigree (117 dams, 487 sires, and 1351 offspring) and then using an “animal model” approach to estimate quantitative genetic parameters. We find significant additive genetic (co)variance, and hence moderate heritability, for juvenile length and mass. We also find substantial maternal effects for these traits at age-0, but not age-1, confirming that genotype-phenotype interactions between mother and offspring are strongest at birth; although these effects could not be parsed into their genetic and nongenetic components. Our results suggest that human-imposed selection pressures (e.g., size-selective harvesting) might impose noteworthy evolutionary change even in large marine vertebrates. We therefore use our findings to explain how maternal effects may sometimes promote maladaptive juvenile traits, and how lemon sharks at different nursery sites may show “constrained local adaptation.” We also show how single-generation pedigrees, and even simple marker-based regression methods, can provide accurate estimates of quantitative genetic parameters in at least some natural systems.

During the process of coevolution, social parasites have evolved sophisticated strategies to exploit the brood care behavior of their social hosts. Slave-making ant queens invade host colonies and kill or eject all adult host ants. Host workers, which eclose from the remaining brood, are tricked into caring for the parasite brood. Due to their high prevalence and frequent raids, following which stolen host broods are similarly enslaved, slave-making ants exert substantial selection upon their hosts, leading to the evolution of antiparasite adaptations. However, all host defenses shown to date are active before host workers are parasitized, whereas selection was thought to be unable to act on traits of already enslaved hosts. Yet, here we demonstrate the rebellion of enslaved Temnothorax workers, which kill two-thirds of the female pupae of the slave-making ant Protomognathus americanus. Thereby, slaves decrease the long-term parasite impact on surrounding related host colonies. This novel antiparasite strategy of enslaved workers constitutes a new level in the coevolutionary battle after host colony defense has failed. Our discovery is analogous to recent findings in hosts of avian brood parasites where perfect mimicry of parasite eggs leads to the evolution of chick recognition as a second line of defense.

Most models of Fisherian sexual selection assume haploidy. However, analytical models that focus on dynamics near fixation boundaries and simulations show that the resulting behavior depends on ploidy. Here we model sexual selection in a diploid to characterize behaviour away from fixation boundaries. The model assumes two di-allelic loci, a male-limited trait locus subject to viability selection, and a preference locus that determines a female's tendency to mate with males based on their genotype at the trait locus. Using a quasi-linkage equilibrium (QLE) approach, we find a general equation for the curves of quasi-neutral equilibria, and the conditions under which they are attracting or repelling. Unlike in the haploid model, the system can move away from the internal curve of equilibria in the diploid model. We show that this is the case when the combined forces of natural and sexual selection induce underdominance at the trait locus.

Many essential organelles and endosymbionts exhibit a strict matrilineal pattern of inheritance. The absence of paternal transmission of such extranuclear components is thought to preclude a response to selection on their effects on male viability and fertility. We overturn this dogma by showing that two mechanisms, inbreeding and kin selection, allow mitochondria to respond to selection on both male viability and fertility. Even modest levels of inbreeding allow such a response to selection when there are direct fitness effects of mitochondria on male fertility because inbreeding associates male fertility traits with mitochondrial matrilines. Male viability effects of mitochondria are also selectable whenever there are indirect fitness effects of males on the fitness of their sisters. When either of these effects is sufficiently strong, we show that there are conditions that allow the spread of mitochondria with direct effects that are harmful to females, contrary to standard expectation. We discuss the implications of our findings for the evolution of organelles and endosymbionts and genomic conflict.

Many evolutionary processes can lead to a change in the correlation between continuous characters over time or on different branches of a phylogenetic tree. Shifts in genetic or functional constraint, in the selective regime, or in some combination there of can influence both the evolution of continuous traits and their relation to each other. These changes can often be mapped on a phylogenetic tree to examine their influence on multivariate phenotypic diversification. We propose a new likelihood method to fit multiple evolutionary rate matrices (also called evolutionary variance-covariance matrices) to species data for two or more continuous characters and a phytogeny. The evolutionary rate matrix is a matrix containing the evolutionary rates for individual characters on its diagonal, and the covariances between characters (of which the evolutionary correlations are a function) elsewhere. To illustrate our approach, we apply the method to an empirical dataset consisting of two features of feeding morphology sampled from 28 centrarchid fish species, as well as to data generated via phylogenetic numerical simulations. We find that the method has appropriate type I error, power, and parameter estimation. The approach presented herein is the first to allow for the explicit testing of how and when the evolutionary covariances between characters have changed in the history of a group.