Patterns of covariation in the hominoid craniofacial skeleton: implications for paleoanthropological models

Living species are often used as analogues for fossil ones. When this is done, the implicit assumption is made that hominids and living hominoids vary in the same way. This paper addresses the validity of this assumption by comparing patterns of facial variation among humans and African apes. In particular, it addresses three major questions that underlie approaches to reconstructing hominid relationships. First, is phenotypic variation similar between closely related species? Second, if it is dissimilar, why? Third, is it feasible to use analogue species for modeling purposes? Measurements are obtained from 542 crania of adult apes and humans. Care is taken to choose homologous data, and account for differences in population size and structure. Variance/covariance and correlation matrices among the species are compared using common principal component (CPC) analysis, random skewers methods and matrix correlations. Morphological distances (D(2)) are calculated between population means, and between randomized pairs of individuals within each population, to evaluate intraspecific variation. Morphological distances are also calculated between randomized pairs of individuals using the variation patterns of analogue populations, in order to evaluate the efficacy of such substitutions. Results show that while the hominoids share a similar pattern of facial variation overall, the patterns do diverge. This difference generally corresponds to the phylogenetic relationships among these species, suggesting that patterns of variation may have diverged through time in the large bodied hominoids. Because interpretation of relationships in the fossil record is confounded by a lack of understanding of how variation changes through time, exploration of such patterns of divergence can provide important clues to understanding human evolution. Additionally, neglecting to account for this divergence when using living analogues as variation "yardsticks" can give rise to interpretations of the fossil record that are more speciose than is warranted.     

Patterns of phenotypic covariation and correlation in modern humans as viewed from morphological integration

Proportionality of phenotypic and genetic distance is of crucial importance to adequately focus on population history and structure, and it depends on the proportionality of genetic and phenotypic covariance. Constancy of phenotypic covariances is unlikely without constancy of genetic covariation if the latter is a substantial component of the former. If phenotypic patterns are found to be relatively stable, the most probable explanation is that genetic covariance matrices are also stable. Factors like morphological integration account for such stability. Morphological integration can be studied by analyzing the relationships among morphological traits. We present here a comparison of phenotypic correlation and covariance structure among worldwide human populations. Correlation and covariance matrices between 47 cranial traits were obtained for 28 populations, and compared with design matrices representing functional and developmental constraints. Among-population differences in patterns of correlation and covariation were tested for association with matrices of genetic distances (obtained after an examination of 10 Alu-insertions) and with Mahalanobis distances (computed after craniometrical traits). All matrix correlations were estimated by means of Mantel tests. Results indicate that correlation and covariance structure in our species is stable, and that among-group correlation/covariance similarity is not related to genetic or phenotypic distance. Conversely, genetic and morphological distance matrices were highly correlated. Correlation and covariation patterns were largely associated with functional and developmental factors, which probably account for the stability of covariance patterns.     

Comparison of cranial ontogenetic trajectories among great apes and humans

Molecular data suggest that humans are more closely related to chimpanzees than either is to the gorillas, yet one finds the closest similarity in craniofacial morphology to be among the great apes to the exclusion of humans. To clarify how and when these differences arise in ontogeny, we studied ontogenetic trajectories for Homo sapiens, Pan paniscus, Pan troglodytes, Gorilla gorilla and Pongo pygmaeus. A total of 96 traditional three-dimensional landmarks and semilandmarks on the face and cranial base were collected on 268 adult and sub-adult crania for a geometric morphometric analysis. The ontogenetic trajectories are compared by various techniques, including a new method, relative warps in size-shape space. We find that adult Homo sapiens specimens are clearly separated from the great apes in shape space and size-shape space. Around birth, Homo sapiens infants are already markedly different from the great apes, which overlap at this age but diverge among themselves postnatally. The results suggest that the small genetic differences between Homo and Pan affect early human ontogeny to induce the distinct adult human craniofacial morphology. Pure heterochrony does not sufficiently explain the human craniofacial morphology nor the differences among the African apes.     

Ontogeny and phylogeny of the pelvis in Gorilla, Pongo, Pan, Australopithecus and Homo

To examine the evolutionary differences between hominoid locomotor systems, a number of observations concerning the growth of the pelvis among the great apes as compared to modern and fossil hominids are reported. We are interested in the size and shape of the coxal bones at different developmental stages across species that may elucidate the relationship between ontogeny and phylogeny (i.e., heterochrony) in the hominoid pelvis. Our hypotheses are: (1) do rates of absolute growth differ?, (2) do rates of relative growth differ?, and (3) does heterochrony explain these differences? Bivariate and multivariate analyses of pelvic dimensions demonstrate both the diversity of species-specific ontogenetic patterns among hominoids, and an unequivocal separation of hominids and the great apes. Heterochrony alone fails to account for the ontogenetic differences between hominids and the great apes. Compared to recent Homo,Australopithecus can be described as 'hyper-human' from the relative size of the ischium, and short but broad ilium. Australopithecus afarensis differs from Australopithecus africanus by its relatively long pubis. In multivariate analyses of ilium shape, the most complete coxal bone attributed to Homo erectus, KNM-ER 3228, falls within the range of juvenile and adult Australopithecus, whereas Broken Hill falls within the range of modern Homo, suggesting that the modern human ilium shape arose rather recently. Among the great apes, patterns of pelvic ontogeny do not exclusively separate the African apes from Pongo.     

Morphological Integration and the Interpretation of Fossil Hominin Diversity

The fossil record of primate and human evolution cannot provide accurate estimates of within species variation and integration. This means that we cannot directly observe how patterns of integration have evolved over time in this lineage. And yet, our interpretations of fossil diversity are awash with assumptions about variation patterning in precisely these fossil taxa. Most commonly, researchers rely on extant models of variation for interpreting past diversity, by assuming equality of variation (and occasionally covariation) among extant and fossil populations. Yet one of the things we know from studies of integration in primates is that patterns of morphological covariation can differ among even closely related taxa, indicating that they have diverged over evolutionary time, either in response to selection or as the result of neutral evolution. At the same time, overall patterns of integration remain remarkably similar, meaning that in many respects they are highly conserved evolutionarily. Taken together, these seemingly contradictory observations offer an important conceptual framework for interpreting patterns that we observe in the fossil past. This framework dictates that while we can use patterns of covariation in extant taxa as proxies for extinct diversity, and indeed their conserved nature makes them superior to approaches that rely on variation alone, we also need to account for the fact that such patterns change over time, and incorporate that into our models. Here I provide examples using covariation patterns estimated from modern humans and African great apes to demonstrate the extent to which divergence in covariance structure might affect our interpretations of hominin diversity.

HOW SIMILAR ARE GENETIC CORRELATION STRUCTURES? DATA FROM MICE AND RATS

An integral assumption of many models of morphometric evolution is the equality of the genetic variance-covariance structure across evolutionary time. To examine this assumption, the quantitative-genetic aspects of morphometric form are examined for eight pelvic traits in laboratory rats (Rattus norvegicus) and random-bred ICR mice (Mus musculus). In both species, all traits are significantly heritable, and there are significant phenotypic and genetic correlations among traits, although environmental correlations among the eight traits are low. The size relations among the pelvic variables are isometric. Three matrix-permutation tests are used to examine similarity of phenotypic, genetic, and environmental covariance and correlation matrices within and between species. Independent patterns of morphometric covariation and correlation arise from genetic and environmental effects within each species and from environmental effects between species. The patterns of phenotypic and genetic covariation and correlation are similar within each species, and the phenotypic and genetic correlations are also similar between these species. However, genetic covariance matrices show no significant statistical association between species. It is suggested that the assumption of equality of genetic variance-covariance structures across divergent taxa should be approached with caution.     

Covariance structure in the skull of Catarrhini: a case of pattern stasis and magnitude evolution

The study of the genetic variance/covariance matrix (G-matrix) is a recent and fruitful approach in evolutionary biology, providing a window of investigating for the evolution of complex characters. Although G-matrix studies were originally conducted for microevolutionary timescales, they could be extrapolated to macroevolution as long as the G-matrix remains relatively constant, or proportional, along the period of interest. A promising approach to investigating the constancy of G-matrices is to compare their phenotypic counterparts (P-matrices) in a large group of related species; if significant similarity is found among several taxa, it is very likely that the underlying G-matrices are also equivalent. Here we study the similarity of covariance and correlation structure in a broad sample of Old World monkeys and apes (Catarrhini). We made phylogenetically structured comparisons of correlation and covariance matrices derived from 39 skull traits, ranging from between species to the superfamily level. We also compared the overall magnitude of integration between skull traits (r²) for all Catarrhini genera. Our results show that P-matrices were not strictly constant among catarrhines, but the amount of divergence observed among taxa was generally low. There was significant and positive correlation between the amount of divergence in correlation and covariance patterns among the 30 genera and their phylogenetic distances derived from a recently proposed phylogenetic hypothesis. Our data demonstrate that the P-matrices remained relatively similar along the evolutionary history of catarrhines, and comparisons with the G-matrix available for a New World monkey genus (Saguinus) suggests that the same holds for all anthropoids. The magnitude of integration, in contrast, varied considerably among genera, indicating that evolution of the magnitude, rather than the pattern of inter-trait correlations, might have played an important role in the diversification of the catarrhine skull.     

Evolution of a complex phenotype with biphasic ontogeny: Contribution of development versus function and climatic variation to skull modularity in toads

The theory of morphological integration and modularity predicts that if functional correlations among traits are relevant to mean population fitness, the genetic basis of development will be molded by stabilizing selection to match functional patterns. Yet, how much functional interactions actually shape the fitness landscape is still an open question. We used the anuran skull as a model of a complex phenotype for which we can separate developmental and functional modularity. We hypothesized that functional modularity associated to functional demands of the adult skull would overcome developmental modularity associated to bone origin at the larval phase because metamorphosis would erase the developmental signal. We tested this hypothesis in toad species of the Rhinella granulosa complex using species phenotypic correlation pattern (P‐matrices). Given that the toad species are distributed in very distinct habitats and the skull has important functions related to climatic conditions, we also hypothesized that differences in skull trait covariance pattern are associated to differences in climatic variables among species. Functional and hormonal‐regulated modules are more conspicuous than developmental modules only when size variation is retained on species P‐matrices. Without size variation, there is a clear modularity signal of developmental units, but most species have the functional model as the best supported by empirical data without allometric size variation. Closely related toad species have more similar climatic niches and P‐matrices than distantly related species, suggesting phylogenetic niche conservatism. We infer that the modularity signal due to embryonic origin of bones, which happens early in ontogeny, is blurred by the process of growth that occurs later in ontogeny. We suggest that the species differing in the preferred modularity model have different demands on the orbital functional unit and that species contrasting in climate are subjected to divergent patterns of natural selection associated to neurocranial allometry and T3 hormone regulation.     

Ontogenetic variation in the mandibular ramus of great apes and humans

Considerable variation exists in mandibular ramus form among primates, particularly great apes and humans. Recent analyses of adult ramal morphology have suggested that features on the ramus, especially the coronoid process and sigmoid notch, can be treated as phylogenetic characters that can be used to reconstruct relationships among great ape and fossil hominin taxa. Others have contended that ramal morphology is more influenced by function than phylogeny. In addition, it remains unclear how ontogeny of the ramus contributes to adult variation in great apes and humans. Specifically, it is unclear whether differences among adults appear early and are maintained throughout ontogeny, or if these differences appear, or are enhanced, during later development. To address these questions, the present study examined a broad ontogenetic sample of great apes and humans using two‐dimensional geometric morphometric analysis. Variation within and among species was summarized using principal component and thin plate spline analyses, and Procrustes distances and discriminant function analyses were used to statistically compare species and age classes. Results suggest that morphological differences among species in ramal morphology appear early in ontogeny and persist into adulthood. Morphological differences among adults are particularly pronounced in the height and angulation of the coronoid process, the depth and anteroposterior length of the sigmoid notch, and the inclination of the ramus. In all taxa, the ascending ramus of the youngest specimens is more posteriorly inclined in relation to the occlusal plane, shifting to become more upright in adults. These results suggest that, although there are likely functional influences over the form of the coronoid process and ramus, the morphology of this region can be profitably used to differentiate among great apes, modern humans, and fossil hominid taxa. J. Morphol. 275:661–677, 2014. © 2013 Wiley Periodicals, Inc.     

Developmental regulation of skull morphology II: ontogenetic dynamics of covariance

Canalization may play a critical role in molding patterns of integration when variability is regulated by the balance between processes that generate and remove variation. Under these conditions, the interaction among those processes may produce a dynamic structure of integration even when the level of variability is constant. To determine whether the constancy of variance in skull shape throughout most of postnatal growth results from a balance between processes generating and removing variation, we compare covariance structures from age to age in two rodent species, cotton rats (Sigmodon fulviventer) and house mice (Mus musculus domesticus). We assess the overall similarity of covariance matrices by the matrix correlation, and compare the structures of covariance matrices using common subspace analysis, a method related to common principal components (PCs) analysis but suited to cases in which variation is so nearly spherical that PCs are ambiguous. We find significant differences from age to age in covariance structure and the more effectively canalized ones tend to be least stable in covariance structure. We find no evidence that canalization gradually and preferentially removes deviations arising early in development as we might expect if canalization results from compensatory differential growth. Our results suggest that (co)variation patterns are continually restructured by processes that equilibrate variance, and thus that canalization plays a critical role in molding patterns of integration.     

ONTOGENETIC VARIATION IN PATTERNS OF DEVELOPMENTAL AND FUNCTIONAL INTEGRATION IN SKULLS OF SIGMODON FULVIVENTER

Patterns of variation and covariation within populations can influence how characters respond to natural selection and random genetic drift and so constrain the ability of natural selection to modify the phenotype. We examined several potential developmental and functional explanations of character covariation throughout ontogeny using known-age samples of the cotton rat (Sigmodon fulviventer) to identify the causes of covariation and to assess the variability of patterns of covariation throughout postnatal growth. Competing developmental and functional models were fit to samples of orofacial and neurocranial measures by confirmatory factor analysis and evaluated for their ability to reconstruct observed variance-covariance matrices. Samples of successive ages were simultaneously fit to a common model to test the hypothesis that the patterns of developmental and functional integration were invariant between ages. Orofacial characters derived from the same branchial-arch primordium covary early in ontogeny. Subsequently, there is a repatterning of integration that may reflect a transition from developmental to functional sources of integration. Neurocranial characters exhibit even more variation in patterns of covariation: initially, characters appear to comprise a single integrated unit; before puberty, they appear to respond to localized bone growth; after puberty, they form separate calvarial and basicranial components. This ontogenetic variation in patterns of covariation suggests that developmental constraints are transient and flexible and that the consequences of selection may depend upon the age at which it acts.     

Causes of covariation of phenotypic traits among populations

Morphological and life-history traits often vary among populations of a species. Traits generally do not vary independently, but show patterns of covariation that can arise from genetic and environmental influences on phenotype. Covariance of traits may arise at an among-population level when genetically influenced traits diverge among populations in a correlated manner. Genetic correlations caused by pleiotropy and/or gene linkage can cause traits to evolve together, but among-population covariance can also arise among traits that are not genetically correlated. For example, “selective covariance” can arise when natural selection directly causes correlated change in a suite of traits. Similarly, mutation, migration, and drift may also sometimes cause correlated genetic changes among populations. Because covariation of traits among populations can arise by several different processes, the evolution of suites of traits must be interpreted with great caution. We discuss the sources of among-population covariance and illustrate one approach to identifying the sources' using data on floral traits of Dalechampia scandens (Euphorbiaceae).     

ONTOGENETIC VARIATION IN PATTERNS OF PHENOTYPIC INTEGRATION IN THE LABORATORY RAT

I used confirmatory factor analysis to evaluate the ability of causal developmental models to predict observed phenotypic integration in limb and skull measures at five stages of postnatal ontogeny in the laboratory rat. To analyze the dynamics of phenotypic integration, I fit successive age-classes simultaneously to a common model. Growth was the principal developmental explanation of observed phenotypic covariation in the limb and skull. No complex morphogenetic model more adequately reconstructed observed covariance structure. Models that could not be interpreted in embryological terms, coupled with a growth component, provide the best models for observed phenotypic integration. During postnatal growth, some aspects of integration vary in both the skull and limb. The covariance between factors and the proportion of variance unique to each character differ between some sequential age-classes. The factor-pattern is invariant in the limb; however, repatterning in the skull occurs in the interval between eye-opening and weaning. The temporal variation in the structure of covariation suggests that functional interactions among characters may create observed patterns of phenotypic integration. The developmental constraints responsible for evolutionary modification of phenotypes might be equally dynamic and responsive to embryonic functional interactions.     

Genetic covariance structure of growth in the salamander Ambystoma macrodactylum

The size of an organism at any point during ontogeny often has fitness consequences through either direct selection on size or through selection on size-related morphological, performance, or life history traits. However, the evolutionary response to selection on size across ontogeny (a growth trajectory) may be limited by genetic correlations across ages. Here we characterize the phenotypic and genetic covariance structure of length and mass growth trajectories in a natural population of larval Ambystoma macrodactylum using function-valued quantitative genetic analyses and principal component decomposition. Most of the phenotypic and genetic variation in both growth trajectories appears to be confined to a single principal component describing a pattern of positive covariation among sizes across all ages. Higher order principal components with no significant associated genetic variation were identified for both trajectories, suggesting that evolution towards certain patterns of negative covariation between sizes across ages is constrained. The well-characterized positive relationship between size at metamorphosis and fitness in pond-breeding amphibians predicts that the across-age covariance structure will strongly limit evolution only if there is negative selection on size prior to metamorphosis. The pattern of genetic covariation observed in this study is similar to that observed in other vertebrate taxa, indicating that size may often be highly genetically and phenotypically integrated across ontogeny. Additionally, we find that phenotypic and genetic analyses of growth trajectories can yield qualitatively similar patterns of covariance structure.     

Hominins do not share a common postnatal facial ontogenetic shape trajectory

This paper examines the hypothesis raised by recent studies that postnatal trajectories of shape change in the facial skeleton are parallel between, at least, chimpanzees, modern humans and also fossil hominins, specifically australopithecines and possibly Neanderthals. In contrast, other studies point to divergences in postnatal shape trajectories within diverse groups of primates. As such there is some debate regarding the relative contributions of pre and postnatal ontogeny to adult morphological differences. This paper presents a series of geometric morphometric studies of the ontogeny of facial shape in hominins with the specific aim of resolving these issues. The results indicate that many differences in facial shape between hominins are established prenatally, however highly significant divergences of postnatal facial ontogeny are found among living hominins. Our studies point to possible differences between the shape ontogeny of the Australopithecus africanus face and that of African apes on the one hand and humans on the other. However, sampling experiments indicate that the small sample size of available specimens of A. africanus does not permit any conclusions to be drawn regarding comparative shape ontogeny of the face.     

Enthesis bilateral asymmetry in humans and African apes

Entheses (skeletal muscle and tendon attachment sites) have often been used to infer handedness and activity variability among human populations. However, the specific roles that intensity vs. frequency of muscle contractions play in modifying entheses are not well understood and the assumption that entheses reflect muscle activity levels has been challenged. This study explores the effect of habitual muscular activity on enthesis morphology in humans and African apes by investigating bilateral asymmetry in the forelimbs and hindlimbs of these taxa. Humans have generally more developed entheses in the lower limb while African apes have generally more developed entheses in the forelimbs. All species studied have more asymmetric forelimbs than hindlimbs except humans that show more asymmetrical expression of bony spurs in the lower limbs than in the upper limbs. When comparing species, humans are always more asymmetric in ethesis development than apes for both the forelimbs and hindlimbs, which reflects the relatively greater asymmetry in limb use in humans and the more symmetric use in apes. Enthesis development may reflect cross-symmetry patterns in humans and, more subtly, a moderate handedness in apes during manipulative activities. This study suggests that enthesis morphology provides information on muscle activity levels, with greater development of entheses associated with more habitual or powerful muscle use. The general similarity of ape and human responses to muscle activity suggests that muscle activity influenced enthesis development in Plio-Pleistocene hominins and that interpretation of muscle markings in these fossils can provide data for functional inferences in these extinct species.     

A comparison of phenotypic variation and covariation patterns and the role of phylogeny, ecology, and ontogeny during cranial evolution of new world monkeys

Similarity of genetic and phenotypic variation patterns among populations is important for making quantitative inferences about past evolutionary forces acting to differentiate populations and for evaluating the evolution of relationships among traits in response to new functional and developmental relationships. Here, phenotypic co variance and correlation structure is compared among Platyrrhine Neotropical primates. Comparisons range from among species within a genus to the superfamily level. Matrix correlation followed by Mantel's test and vector correlation among responses to random natural selection vectors (random skewers) were used to compare correlation and variance/covariance matrices of 39 skull traits. Sampling errors involved in matrix estimates were taken into account in comparisons using matrix repeatability to set upper limits for each pairwise comparison. Results indicate that covariance structure is not strictly constant but that the amount of variance pattern divergence observed among taxa is generally low and not associated with taxonomic distance. Specific instances of divergence are identified. There is no correlation between the amount of divergence in covariance patterns among the 16 genera and their phylogenetic distance derived from a conjoint analysis of four already published nuclear gene datasets. In contrast, there is a significant correlation between phylogenetic distance and morphological distance (Mahalanobis distance among genus centroids). This result indicates that while the phenotypic means were evolving during the last 30 millions years of New World monkey evolution, phenotypic covariance structures of Neotropical primate skulls have remained relatively consistent. Neotropical primates can be divided into four major groups based on their feeding habits (fruit-leaves, seed-fruits, insect-fruits, and gum-insect-fruits). Differences in phenotypic covariance structure are correlated with differences in feeding habits, indicating that to some extent changes in interrelationships among skull traits are associated with changes in feeding habits. Finally, common patterns and levels of morphological integration are found among Platyrrhine primates, suggesting that functional/developmental integration could be one major factor keeping covariance structure relatively stable during evolutionary diversification of South American monkeys.     

Postnatal Cranial Development in Papionin Primates: An Alternative Model for Hominin Evolutionary Development

The evolution of hominin growth and life history has long been a subject of intensive research, but it is only recently that paleoanthropologists have considered the ontogenetic basis of human morphological evolution. To date, most human EvoDevo studies have focused on developmental patterns in extant African apes and humans. However, the Old World monkey tribe Papionini, a diverse clade whose members resemble hominins in their ecology and population structure, has been proposed as an alternative model for human craniofacial evolution. This paper reviews prior studies of papionin development and socioecology and presents new analyses of juvenile shape variation and ontogeny to address fundamental questions concerning primate cranial development, including: (1) When are cranial shape differences between species established? (2) How do epigenetic influences modulate early-arising pattern differences? (3) How much do postnatal developmental trajectories vary? (4) What is the impact of developmental variation on adult cranial shape? and, (5) What role do environmental factors play in establishing adult cranial form? Results of this inquiry suggest that species differences in cranial morphology arise during prenatal or earliest postnatal development. This is true even for late-arising features that develop under the influence of epigenetic factors such as mechanical loading. Papionins largely retain a shared, ancestral pattern of ontogenetic shape change, but large size and sexual dimorphism are associated with divergent developmental trajectories, suggesting differences in cranial integration. Developmental simulation studies indicate that postnatal ontogenetic variation has a limited influence on adult cranial morphology, leaving early morphogenesis as the primary determinant of cranial shape. The ability of social factors to influence craniofacial development in Mandrillus suggests a possible role for phentotypic plasticity in the diversification of primate cranial form. The implications of these findings for taxonomic attribution of juvenile fossils, the developmental basis of early hominin characters, and hominin cranial diversity are discussed.     

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.