Evolutionary selection and morphological integration in the foot of modern humans

Objectives

To advance our understanding of the evolution of the hominin foot by quantifying integration and responses to selection in the foot of modern humans.

Materials and Methods

The sample includes 247 female and male adult individuals from Euro-American, Afro-American, European, and Amerindian populations. We collected 190 linear measurements from the 26 skeletal elements that constitute the modern human foot. With these data, we calculated the magnitudes of integration and the ability of the foot to respond to selection demands.

Results

The results revealed that distal phalanges are less integrated, more evolvable, and more flexible than proximal elements (i.e., proximal phalanges and metatarsals). Also, bones from the medial ray (e.g., hallux) show stronger integration and weaker evolvability than their counterparts from the lateral column (e.g., fifth ray), following this trend from medial to lateral positions. Among the tarsals, the talus and calcaneus are the most integrated, least evolvable, and flexible elements from that module.

Discussion

These results suggest that selection for bipedalism would have reorganized the variance/covariance matrix of the foot. The hallux might have been under strong functional selection pressures for bipedal requirements, resulting in a strong integration and low evolvability. Also, differences in the developmental process of each bone seem to have played an essential role in the degree of evolvability, showing those elements that develop earlier have less ability to respond to selection demands.     

Evolution of adaptive phenotypic variation patterns by direct selection for evolvability

A basic assumption of the Darwinian theory of evolution is that heritable variation arises randomly. In this context, randomness means that mutations arise irrespective of the current adaptive needs imposed by the environment. It is broadly accepted, however, that phenotypic variation is not uniformly distributed among phenotypic traits, some traits tend to covary, while others vary independently, and again others barely vary at all. Furthermore, it is well established that patterns of trait variation differ among species. Specifically, traits that serve different functions tend to be less correlated, as for instance forelimbs and hind limbs in bats and humans, compared with the limbs of quadrupedal mammals. Recently, a novel class of genetic elements has been identified in mouse gene-mapping studies that modify correlations among quantitative traits. These loci are called relationship loci, or relationship Quantitative Trait Loci (rQTL), and affect trait correlations by changing the expression of the existing genetic variation through gene interaction. Here, we present a population genetic model of how natural selection acts on rQTL. Contrary to the usual neo-Darwinian theory, in this model, new heritable phenotypic variation is produced along the selected dimension in response to directional selection. The results predict that selection on rQTL leads to higher correlations among traits that are simultaneously under directional selection. On the other hand, traits that are not simultaneously under directional selection are predicted to evolve lower correlations. These results and the previously demonstrated existence of rQTL variation, show a mechanism by which natural selection can directly enhance the evolvability of complex organisms along lines of adaptive change.     

Genotype-Phenotype Maps Maximizing Evolvability: Modularity Revisited

The mechanisms translating genetic to phenotypic variation determine the distribution of heritable phenotypic variance available to selection. Pleiotropy is an aspect of this structure that limits independent variation of characters. Modularization of pleiotropy has been suggested to promote evolvability by restricting genetic covariance among unrelated characters and reducing constraints due to correlated response. However, modularity may also reduce total genetic variation of characters. We study the properties of genotype-phenotype maps that maximize average conditional evolvability, measured as the amount of unconstrained genetic variation in random directions of phenotypic space. In general, maximal evolvability occurs by maximizing genetic variance and minimizing genetic covariance. This does not necessarily require modularity, only patterns of pleiotropy that cancel on average. The detailed structure of the most evolvable genotype-phenotype maps depends on the distribution of molecular variance. When molecular variance is determined by mutation-selection equilibrium either highly pleiotropic or highly modular genotype-phenotype maps can be optimal, depending on the mutation rate and the relative strengths of stabilizing selection on the characters.     

Evolution of morphological integration in the skull of Carnivora (Mammalia): Changes in Canidae lead to increased evolutionary potential of facial traits

Morphological integration refers to the fact that different phenotypic traits of organisms are not fully independent from each other, and tend to covary to different degrees. The covariation among traits is thought to reflect properties of the species' genetic architecture and thus can have an impact on evolutionary responses. Furthermore, if morphological integration changes along the history of a group, inferences of past selection regimes might be problematic. Here, we evaluated the stability and evolution of the morphological integration of skull traits in Carnivora by using evolutionary simulations and phylogenetic comparative methods. Our results show that carnivoran species are able to respond to natural selection in a very similar way. Our comparative analyses show that the phylogenetic signal for pattern of integration is lower than that observed for morphology (trait averages), and that integration was stable throughout the evolution of the group. That notwithstanding, Canidae differed from other families by having higher integration, evolvability, flexibility, and allometric coefficients on the facial region. These changes might have allowed canids to rapidly adapt to different food sources, helping to explain not only the phenotypic diversification of the family, but also why humans were able to generate such a great diversity of dog breeds through artificial selection.     

Evolvability and genetic constraint in Dalechampia blossoms: components of variance and measures of evolvability

Many evolutionary arguments are based on the assumption that quantitative characters are highly evolvable entities that can be rapidly moulded by changing selection pressures. The empirical evaluation of this assumption depends on having an operational measure of evolvability that reflects the ability of a trait to respond to a given external selection pressure. We suggest short-term evolvability be measured as expected proportional response in a trait to a unit strength of directional selection, where strength of selection is defined independently of character variation and in units of the strength of selection on fitness itself. We show that the additive genetic variance scaled by the square of the trait mean, IA, is such a measure. The heritability, h2, does not measure evolvability in this sense. Based on a diallel analysis, we use IA to assess the evolvability of floral characters in a population of the neotropical vine Dalechampia scandens (Euphorbiaceae). Although we are able to demonstrate that there is additive genetic variation in a number of floral traits, we also find that most of the traits are not expected to change by more than a fraction of a percent per generation. We provide evidence that the degree of among-population divergence of traits is related to their predicted evolvabilities, but not to their heritabilities.     

Is modularity necessary for evolvability? Remarks on the relationship between pleiotropy and evolvability

Evolvability is the ability to respond to a selective challenge. This requires the capacity to produce the right kind of variation for selection to act upon. To understand evolvability we therefore need to understand the variational properties of biological organisms. Modularity is a variational property, which has been linked to evolvability. If different characters are able to vary independently, selection will be able to optimize each character separately without interference. But although modularity seems like a good design principle for an evolvable organism, it does not therefore follow that it is the only design that can achieve evolvability. In this essay I analyze the effects of modularity and, more generally, pleiotropy on evolvability. Although, pleiotropy causes interference between the adaptation of different characters, it also increases the variational potential of those characters. The most evolvable genetic architectures may often be those with an intermediate level of integration among characters, and in particular those where pleiotropic effects are variable and able to compensate for each other's constraints.     

THE COEVOLUTION OF HUMAN HANDS AND FEET

Human hands and feet have longer, more robust first digits, and shorter lateral digits compared to African apes. These similarities are often assumed to be independently evolved adaptations for manipulative activities and bipedalism, respectively. However, hands and feet are serially homologous structures that share virtually identical developmental blueprints, raising the possibility that digital proportions coevolved in human hands and feet because of underlying developmental linkages that increase phenotypic covariation between them. Here we show that phenotypic covariation between serially homologous fingers and toes in Homo and Pan is not only higher than expected, it also causes these digits to evolve along highly parallel trajectories under episodes of simulated directional selection, even when selection pressures push their means in divergent directions. Further, our estimates of the selection pressures required to produce human‐like fingers and toes from an African ape‐like ancestor indicate that selection on the toes was substantially stronger, and likely led to parallel phenotypic changes in the hands. Our data support the hypothesis that human hands and feet coevolved, and suggest that the evolution of long robust big toes and short lateral toes for bipedalism led to changes in hominin fingers that may have facilitated the emergence of stone tool technology.     

Limited plasticity in the phenotypic variance‐covariance matrix for male advertisement calls in the black field cricket,Teleogryllus commodus

Phenotypic integration and plasticity are central to our understanding of how complex phenotypic traits evolve. Evolutionary change in complex quantitative traits can be predicted using the multivariate breeders’ equation, but such predictions are only accurate if the matrices involved are stable over evolutionary time. Recent study, however, suggests that these matrices are temporally plastic, spatially variable and themselves evolvable. The data available on phenotypic variance‐covariance matrix ( P ) stability are sparse, and largely focused on morphological traits. Here, we compared P for the structure of the complex sexual advertisement call of six divergent allopatric populations of the Australian black field cricket, Teleogryllus commodus. We measured a subset of calls from wild‐caught crickets from each of the populations and then a second subset after rearing crickets under common‐garden conditions for three generations. In a second experiment, crickets from each population were reared in the laboratory on high‐ and low‐nutrient diets and their calls recorded. In both experiments, we estimated P for call traits and used multiple methods to compare them statistically (Flury hierarchy, geometric subspace comparisons and random skewers). Despite considerable variation in means and variances of individual call traits, the structure of P was largely conserved among populations, across generations and between our rearing diets. Our finding that P remains largely stable, among populations and between environmental conditions, suggests that selection has preserved the structure of call traits in order that they can function as an integrated unit.     

Deciphering the Palimpsest: Studying the Relationship Between Morphological Integration and Phenotypic Covariation

Organisms represent a complex arrangement of anatomical structures and individuated parts that must maintain functional associations through development. This integration of variation between functionally related body parts and the modular organization of development are fundamental determinants of their evolvability. This is because integration results in the expression of coordinated variation that can create preferred directions for evolutionary change, while modularity enables variation in a group of traits or regions to accumulate without deleterious effects on other aspects of the organism. Using our own work on both model systems (e.g., lab mice, avians) and natural populations of rodents and primates, we explore in this paper the relationship between patterns of phenotypic covariation and the developmental determinants of integration that those patterns are assumed to reflect. We show that integration cannot be reliably studied through phenotypic covariance patterns alone and argue that the relationship between phenotypic covariation and integration is obscured in two ways. One is the superimposition of multiple determinants of covariance in complex systems and the other is the dependence of covariation structure on variances in covariance-generating processes. As a consequence, we argue that the direct study of the developmental determinants of integration in model systems is necessary to fully interpret patterns of covariation in natural populations, to link covariation patterns to the processes that generate them, and to understand their significance for evolutionary explanation.     

Sexual selection and population divergence II. Divergence in different sexual traits and signal modalities in field crickets (Teleogryllus oceanicus)

Sexual selection can target many different types of traits. However, the relative influence of different sexually selected traits during evolutionary divergence is poorly understood. We used the field cricket Teleogryllus oceanicus to quantify and compare how five traits from each of three sexual signal modalities and components diverge among allopatric populations: male advertisement song, cuticular hydrocarbon (CHC) profiles and forewing morphology. Population divergence was unexpectedly consistent: we estimated the among‐population (genetic) variance‐covariance matrix, D , for all 15 traits, and D_max explained nearly two‐thirds of its variation. CHC and wing traits were most tightly integrated, whereas song varied more independently. We modeled the dependence of among‐population trait divergence on genetic distance estimated from neutral markers to test for signatures of selection versus neutral divergence. For all three sexual trait types, phenotypic variation among populations was largely explained by a neutral model of divergence. Our findings illustrate how phenotypic integration across different types of sexual traits might impose constraints on the evolution of mating isolation and divergence via sexual selection.     

Phenotypic integration and the evolution of signal repertoires: A case study of treefrog acoustic communication

Animal signals are inherently complex phenotypes with many interacting parts combining to elicit responses from receivers. The pattern of interrelationships between signal components reflects the extent to which each component is expressed, and responds to selection, either in concert with or independently of others. Furthermore, many species have complex repertoires consisting of multiple signal types used in different contexts, and common morphological and physiological constraints may result in interrelationships extending across the multiple signals in species’ repertoires. The evolutionary significance of interrelationships between signal traits can be explored within the framework of phenotypic integration, which offers a suite of quantitative techniques to characterize complex phenotypes. In particular, these techniques allow for the assessment of modularity and integration, which describe, respectively, the extent to which sets of traits covary either independently or jointly. Although signal and repertoire complexity are thought to be major drivers of diversification and social evolution, few studies have explicitly measured the phenotypic integration of signals to investigate the evolution of diverse communication systems. We applied methods from phenotypic integration studies to quantify integration in the two primary vocalization types (advertisement and aggressive calls) in the treefrogs Hyla versicolor, Hyla cinerea, and Dendropsophus ebraccatus. We recorded male calls and calculated standardized phenotypic variance–covariance ( P ) matrices for characteristics within and across call types. We found significant integration across call types, but the strength of integration varied by species and corresponded with the acoustic similarity of the call types within each species. H. versicolor had the most modular advertisement and aggressive calls and the least acoustically similar call types. Additionally, P was robust to changing social competition levels in H. versicolor. Our findings suggest new directions in animal communication research in which the complex relationships among the traits of multiple signals are a key consideration for understanding signal evolution.     

Morphological integration and natural selection in the postcranium of wild verreaux's sifaka (Propithecus verreauxi verreauxi)

Morphological integration manifests as strong phenotypic covariation among interacting traits. In this study, a graph-theory approach is used to analyze patterns of morphological integration in a wild population of Verreaux's sifaka (Propithecus verreauxi verreauxi). The motivation for this study is to determine the relative roles of development versus function in shaping patterns of morphological integration in the sifaka postcranium. A developmental and a functional hypothesis of integration are compared with the observed pattern of integration and the fit of these hypotheses is assessed using information theoretic statistics. Correlational selection is also estimated on limb elements. Information theoretic statistics indicate that the developmental hypothesis fits the observed pattern of integration slightly better than the functional hypothesis. Only two pairs of traits experience correlational selection but neither of the traits within each pair are morphologically integrated. The observed pattern of integration contains several trait-trait associations that are specified by both the functional and developmental hypotheses. These results likely reflect the nested covariation structure in which a novel locomotor mode, vertical clinging and leaping, is derived from a primitive quadrupedal morphotype.     

Morphological evolution through integration: a quantitative study of cranial integration in Homo, Pan, Gorilla and Pongo

Morphological integration refers to coordinated variation among traits that are closely related in development and/or function. Patterns of integration can offer important insight into the structural relationship between phenotypic units, providing a framework to address questions about phenotypic evolvability and constraints. Integrative features of the primate cranium have recently become a popular subject of study. However, an important question that still remains under-investigated is: what is the pattern of cranial shape integration among closely related hominoids? To address this question, we conducted a Procrustes-based geometric morphometrics study to quantify and analyze shape covariation patterns between different cranial regions in Homo, Pan, Gorilla and Pongo. A total of fifty-six 3D landmarks were collected on 407 adult individuals. We then sub-divided the landmarks corresponding to cranial units as outlined in the ‘functional matrix hypothesis.’ Sub-dividing the cranium in this manner allowed us to explore patterns of covariation between the face, basicranium and cranial vault, using the two-block partial least squares approach. Our results suggest that integrated shape changes in the hominoid cranium are complex, but that the overall pattern of integration is similar among human and non-human apes. Thus, despite having very distinct morphologies the way in which the face, basicranium and cranial vault covary is shared among these taxa. These results imply that the pattern of cranial integration among hominoids is conserved.     

Local adaptation within a hybrid species

Ecological divergence among populations may be strongly influenced by their genetic background. For instance, genetic admixture through introgressive hybridization or hybrid speciation is likely to affect the genetic variation and evolvability of phenotypic traits. We studied geographic variation in two beak dimensions and three other phenotypic traits of the Italian sparrow (Passer italiae), a young hybrid species formed through interbreeding between house sparrows (P. domesticus) and Spanish sparrows (P. hispaniolensis). We found that beak morphology was strongly influenced by precipitation regimes and that it appeared to be the target of divergent selection within Italian sparrows. Interestingly, however, the degree of parental genetic contribution in the hybrid species had no effect on phenotypic beak variation. Moreover, beak height divergence may mediate genetic differentiation between populations, consistent with isolation-by-adaptation within this hybrid species. The study illustrates how hybrid species may be relatively unconstrained by their admixed genetic background, allowing them to adapt rapidly to environmental variation.     

Constraints on evolution and postcopulatory sexual selection: trade-offs among ejaculate characteristics

Ejaculates function as an integrated unit to ensure male fertility and paternity, can have a complex structure, and can experience multiple episodes of selection. Current studies on the evolution of ejaculates typically focus on phenotypic variation in sperm number, size, or related traits such as testes size as adaptations to postcopulatory male-male competition. However, the evolution of the integrated nature of ejaculate structure and function depends on genetic variation in and covariation between the component parts. Here we report a quantitative genetic study of the components of the ejaculate of the cockroach Nauphoeta cinerea, including those we know to experience postcopulatory sexual selection, in the context of functional integration of ejaculate characters. We use the patterns of genetic variation and covariation to infer how the integration of the functions of the ejaculate constrain and shape its evolution. Ejaculate components were highly variable, showed significant additive genetic variance, and moderate to high evolvability. The level of genetic variation in these characters, despite strong directional or truncating selection, may reflect the integration of multiple episodes of selection that occur in N. cinerea. There were few significant phenotypic correlations, but all the genetic correlations among ejaculate characters were significantly different from zero. The patterns of genetic variation and covariation suggest that there are important trade-offs among individual traits of the ejaculate and that evolution of ejaculate characteristics will not proceed unconstrained. Fully describing the genetic relationships among traits that perform as an integrated unit helps us understand how functional relationships constrain or facilitate the evolution of the complex structure that is the ejaculate.     

Sexually dimorphic levels of color trait integration and the resolution of sexual conflict in Lake Malawi cichlids

East African cichlids are renowned for their propensity to radiate, and variation in color patterns accounts for much of endemic cichlid diversity. Sexual dimorphism in color among cichlid species likely represents the outcome of different selective regimes acting on each sex, and is a classic example of sexual conflict. It is generally assumed that this conflict has been mitigated through the evolution of sex-linked color polymorphisms. Here, we propose that the evolution of sex-specific differences in levels of color trait integration may represent an additional mechanism through which sexual conflict has been resolved in this group. Specifically, we predict: (1) that general patterns of integration are influenced by early developmental events and thus conserved across sexes and (2) that male color is less integrated than females, and thus more evolvable in terms of producing an elaborate palette (i.e., in response to sexual selection), whereas female color is more integrated, facilitating wholesale shifts in color for background matching (i.e., in response to natural selection for crypsis). We tested these hypotheses using an F(2) design to compare the segregation of male and female color patterns. Both exploratory methods and hypothesis-driven analyses of integration demonstrate that the covariance structure of color traits in males and females is distinct, and that males are significantly less integrated than females. We suggest that the ability of species to promote different levels, and to a lesser extent patterns, of phenotypic integration between males and females may have contributed to the evolutionary success of this group.     

Genetic variation in pleiotropy: differential epistasis as a source of variation in the allometric relationship between long bone lengths and body weight

Pleiotropy is an aspect of genetic architecture underlying the phenotypic covariance structure. The presence of genetic variation in pleiotropy is necessary for natural selection to shape patterns of covariation between traits. We examined the contribution of differential epistasis to variation in the intertrait relationship and the nature of this variation. Genetic variation in pleiotropy was revealed by mapping quantitative trait loci (QTLs) affecting the allometry of mouse limb and tail length relative to body weight in the mouse-inbred strain LG/J by SM/J intercross. These relationship QTLs (rQTLs) modify relationships between the traits affected by a common pleiotropic locus. We detected 11 rQTLs, mostly affecting allometry of multiple bones. We further identified epistatic interactions responsible for the observed allometric variation. Forty loci that interact epistatically with the detected rQTLs were identified. We demonstrate how these epistatic interactions differentially affect the body size variance and the covariance of traits with body size. We conclude that epistasis, by differentially affecting both the canalization and mean values of the traits of a pleiotropic domain, causes variation in the covariance structure. Variation in pleiotropy maintains evolvability of the genetic architecture, in particular the evolvability of its modular organization.     

Pervasive genetic integration directs the evolution of human skull shape

It has long been unclear whether the different derived cranial traits of modern humans evolved independently in response to separate selection pressures or whether they resulted from the inherent morphological integration throughout the skull. In a novel approach to this issue, we combine evolutionary quantitative genetics and geometric morphometrics to analyze genetic and phenotypic integration in human skull shape. We measured human skulls in the ossuary of Hallstatt (Austria), which offer a unique opportunity because they are associated with genealogical data. Our results indicate pronounced covariation of traits throughout the skull. Separate simulations of selection for localized shape changes corresponding to some of the principal derived characters of modern human skulls produced outcomes that were similar to each other and involved a joint response in all of these traits. The data for both genetic and phenotypic shape variation were not consistent with the hypothesis that the face, cranial base, and cranial vault are completely independent modules but relatively strongly integrated structures. These results indicate pervasive integration in the human skull and suggest a reinterpretation of the selective scenario for human evolution where the origin of any one of the derived characters may have facilitated the evolution of the others.     

Functional and Genetic Integration in the Skulls of Lake Malawi Cichlids

The level of integration present among organismal traits is thought to influence evolutionary potential, and this potential should be affected by the type or types of integration displayed (e.g., functional, developmental, or genetic). Morphological integration is generally high among functionally related traits, but whether this is predominantly determined by genetic architecture, or is instead a result of biomechanical remodeling during development remains poorly understood. We examine this question in Lake Malawi cichlid fishes by combining a finite-element analysis (FEA) of bite force transmission with quantitative genetic analyses of skull morphology in order to test the hypothesis that functionally coupled traits share a common genetic basis. FEA modeling indicates that the profile of the neurocranium affects its ability to resist forces transmitted from the jaws during biting, and suggests a novel role for skull shape in fish feeding mechanics. Quantitative trait loci mapping demonstrates that the functional integration between jaw and neurocranial shape has a genetic basis, and that this association is being driven by alleles inherited from the specialized biting species. Notably, the co-inheritance of these two functionally related traits in our F₂ matches patterns of covariation within and between Lake Malawi cichlid species. Across species, jaw and neurocranial shapes covary, but the trend appears strongest among biting species. Similarly, within populations of biting species, the dimensions of the jaw and neurocranium are tightly linked, whereas this correlation disappears within populations of omnivorous and suction feeding fish. These data suggest (1) that either pleiotropy, or physical linkage maintained by selection, underlies the phenotypic integration of these two functionally related traits, and (2) that this pattern of integration may have influenced the radiation of craniofacial morphology in Lake Malawi cichlids.     

Sifaka positional behavior: ontogenetic and quantitative genetic approaches

In many primate species, hands and feet are large relative to neonatal body weight, and they subsequently exhibit negative allometric growth during ontogeny. Here, data are presented showing that this pattern holds for a wild population of lemur, Verreaux's sifaka (Propithecus verreauxi verreauxi). Using morphometric data collected on this population, it is shown that younger animals possess relatively large hands and feet. This ontogenetic pattern suggests a simple behavioral test: do juvenile animals with their larger, almost adult‐sized hands and feet locomote on similarly sized substrates as adult animals? Using locomotor bout sampling, this question was tested by collecting positional behavior data on this population. Results from this test find no differences in locomotor behaviors or substrate use between yearlings and adult animals. To place these results in a broader evolutionary context, heritabilities and selection gradients of hands, feet, and other limb elements for animals in this population were estimated. Among limb elements, heritabilities range from 0.16–0.44, with the foot having the lowest value. Positive directional selection acts most strongly on the foot (directional selection gradient = 0.119). The low heritability and positive selection coefficient indicate that selection has acted, and continues to act, on foot size in young animals. These results are interpreted within a functional context with respect to the development of locomotor coordination: larger feet enable young animals to use “adult‐sized” substrates when they move through their habitat. It is suggested that the widespread pattern of negative allometry of the extremities in sifaka and other primates is maintained by selection, and does not simply reflect a primitive developmental pathway that has no adaptive basis. Am J Phys Anthropol 131:261–271, 2006. © 2006 Wiley‐Liss, Inc.