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.     

Disintegrating the fly: A mutational perspective on phenotypic integration and covariation

The structure of environmentally induced phenotypic covariation can influence the effective strength and magnitude of natural selection. Yet our understanding of the factors that contribute to and influence the evolutionary lability of such covariation is poor. Most studies have either examined environmental variation without accounting for covariation, or examined phenotypic and genetic covariation without distinguishing the environmental component. In this study, we examined the effect of mutational perturbations on different properties of environmental covariation, as well as mean shape. We use strains of Drosophila melanogaster bearing well‐characterized mutations known to influence wing shape, as well as naturally derived strains, all reared under carefully controlled conditions and with the same genetic background. We find that mean shape changes more freely than the covariance structure, and that different properties of the covariance matrix change independently from each other. The perturbations affect matrix orientation more than they affect matrix eccentricity or total variance. Yet, mutational effects on matrix orientation do not cluster according to the developmental pathway that they target. These results suggest that it might be useful to consider a more general concept of “decanalization,” involving all aspects of variation and covariation.     

A COMPARISON OF COVARIANCE STRUCTURE IN WILD AND LABORATORY MUROID CRANIA

Mutations have the ability to produce dramatic changes to covariance structure by altering the variance of covariance-generating developmental processes. Several evolutionary mechanisms exist that may be acting interdependently to stabilize covariance structure, despite this developmental potential for variation within species. We explore covariance structure in the crania of laboratory mouse mutants exhibiting mild-to-significant developmental perturbations of the cranium, and contrast it with covariance structure in related wild muroid taxa. Phenotypic covariance structure is conserved among wild muroidea, but highly variable and mutation-dependent within the laboratory group. We show that covariance structures in natural populations of related species occupy a more restricted portion of covariance structure space than do the covariance structures resulting from single mutations of significant effect or the almost nonexistent genetic differences that separate inbred mouse strains. Our results suggest that developmental constraint is not the primary mechanism acting to stabilize covariance structure, and imply a more important role for other mechanisms.     

Impact of Hybridization on Shape,Variation and Covariation of the Mouse Molar

Hybridization can generate phenotypes departing from the parental range through many processes: increased variance due to increased heterozygosity, decreased phenotypic integration and altered covariation structure due to changes in the expression of (co)variation generating developmental processes, and transgression that may arise from epistasis and compensatory genes. Morphometric assessment of shape differentiation, variance and covariance may shed light on these processes, especially for complex characters for which the genetic basis has not yet been assessed. The shape of the first upper molar was thus quantified in a cross between inbred strains of the two European subspecies of the house mouse (Mus musculus domesticus and M. m. musculus). Hybrids appeared as moderately transgressive. Morphological variance was increased in F2s, back to levels observed in wild populations. The pattern of variance–covariance was different between the two parental strains, but restored to a wild-type structure in F2s. Finally, F2s displayed a degree of morphological integration comparable to wild populations but lower than observed in the parental strains. This is interpreted as the result of the reshuffling of the standing genetic variation in hybrids that should have restored the expression of (co)variance generating processes made ineffective in parents due to random allele fixation of some loci. Inter-parental differentiation was more important in regions of the tooth developing early during embryogenesis, whereas transgression was more pronounced in late developing regions. Mutations may more easily generate important geometric differences early during tooth development, occurring as a cascade of signalization from the first cusp to initiate onward. Epistasis and constraints of neighboring parts may be more important later, promoting transgression and impeding differentiation.     

Morphological integration and developmental progress during fish ontogeny in two contrasting habitats

SUMMARY Morphological integration can respond to environmental conditions, a response that may be dynamic through ontogeny. Among fishes, brook charrs ( Salvelinus fontinalis ) display a trophic polymorphism that makes it a good species for analyzing the ontogeny of morphological integration. To better understand the processes regulating variation and integration, we assess the ontogenetic dynamics of covariances and developmental progress for populations of S. fontinalis from two habitats that differ in water velocity; lake and stream. Geometric morphometrics and developmental progress were evaluated on 751 and 198 specimens, respectively. In both habitats, most ossification events occur before the transition from alevin to juvenile. This threshold defines two distinct periods. During the first period representing free-embryos and alevins, there are important shape changes and rapid ossification, integration tends to be relatively low and decreasing and the variance of shape drastically decreases. During the juvenile period, the rate of shape change decreases and the onset of ossification is nearly complete, plus integration increases and shape variance stabilizes. While we find two distinct developmental periods, we nonetheless find a notable stability underlying the ontogenetic dynamics of variability as well as gradual change in the structure of covariation within each habitat. Our results imply that the variability of juvenile body shape does not seem to retain signals of variability determined early in ontogeny and warrants caution in using juvenile as guides to the earlier causes of variability. Overall, this study highlights the difficulty of inferring causes of integration from studies of static covariance.     

How to Explore Morphological Integration in Human Evolution and Development?

Most studies in evolutionary developmental biology focus on large-scale evolutionary processes using experimental or molecular approaches, whereas evolutionary quantitative genetics provides mathematical models of the influence of heritable phenotypic variation on the short-term response to natural selection. Studies of morphological integration typically are situated in-between these two styles of explanation. They are based on the consilience of observed phenotypic covariances with qualitative developmental, functional, or evolutionary models. Here we review different forms of integration along with multiple other sources of phenotypic covariances, such as geometric and spatial dependencies among measurements. We discuss one multivariate method [partial least squares analysis (PLS)] to model phenotypic covariances and demonstrate how it can be applied to study developmental integration using two empirical examples. In the first example we use PLS to study integration between the cranial base and the face in human postnatal development. Because the data are longitudinal, we can model both cross-sectional integration and integration of growth itself, i.e., how cross-sectional variance and covariance is actually generated in the course of ontogeny. We find one factor of developmental integration (connecting facial size and the length of the anterior cranial base) that is highly canalized during postnatal development, leading to decreasing cross-sectional variance and covariance. A second factor (overall cranial length to height ratio) is less canalized and leads to increasing (co)variance. In a second example, we examine the evolutionary significance of these patterns by comparing cranial integration in humans to that in chimpanzees.     

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.     

No relationship between canalization and developmental stability of the skull in a natural population of Mastomys natalensis (Rodentia: Muridae)

The aim of the present work was to investigate the relationship between canalization and developmental stability under varying environmental conditions. Three different cohorts of Mastomys natalensis (Rodentia, Muridae), displaying different growth trajectories, were analysed by means of geometric morphometrics. A set of 23 landmarks was digitalized on the dorsal skull of 292 specimens from Morogoro (Tanzania). Patterns of among‐ and within‐individual (measured as fluctuating asymmetry, FA) variation were assessed and compared among and within the three groups to test for the presence of a common mechanism between canalization and developmental stability. Results showed that there was no congruence between canalization and developmental stability: (1) levels of FA and among‐individual variation varied in a discordant fashion, (2) no correspondence between the variance–covariance matrix of among‐ and within individual variation was found, and (3) environmental effects were able to alter the covariance structure of among‐individual variation leaving patterns associated with fluctuating asymmetry unaffected. These findings support the view of multiple mechanisms underlying developmental buffering of shape variation. © 2011 The Linnean Society of London, Biological Journal of the Linnean Society, 2011, 104 , 207–216.     

Ontogeny of additive and maternal genetic effects: lessons from domestic mammals

Evolution of size and growth depends on heritable variation arising from additive and maternal genetic effects. Levels of heritable (and nonheritable) variation might change over ontogeny, increasing through "variance compounding" or decreasing through "compensatory growth." We test for these processes using a meta-analysis of age-specific weight traits in domestic ungulates. Generally, mean standardized variance components decrease with age, consistent with compensatory growth. Phenotypic convergence among adult sheep occurs through decreasing environmental and maternal genetic variation. Maternal variation similarly declines in cattle. Maternal genetic effects are thus reduced with age (both in absolute and relative terms). Significant trends in heritability (decreasing in cattle, increasing in sheep) result from declining maternal and environmental components rather than from changing additive variation. There was no evidence for increasing standardized variance components. Any compounding must therefore be masked by more important compensatory processes. While extrapolation of these patterns to processes in natural population is difficult, our results highlight the inadequacy of assuming constancy in genetic parameters over ontogeny. Negative covariance between direct and maternal genetic effects was common. Negative correlations with additive and maternal genetic variances indicate that antagonistic pleiotropy (between additive and maternal genetic effects) may maintain genetic variance and limit responses to selection.     

Oceanographic variation influences spatial genomic structure in the sea scallop,Placopecten magellanicus

Environmental factors can influence diversity and population structure in marine species and accurate understanding of this influence can both improve fisheries management and help predict responses to environmental change. We used 7163 SNPs derived from restriction site‐associated DNA sequencing genotyped in 245 individuals of the economically important sea scallop, Placopecten magellanicus, to evaluate the correlations between oceanographic variation and a previously identified latitudinal genomic cline. Sea scallops span a broad latitudinal area (>10 degrees), and we hypothesized that climatic variation significantly drives clinal trends in allele frequency. Using a large environmental dataset, including temperature, salinity, chlorophyll a, and nutrient concentrations, we identified a suite of SNPs (285–621, depending on analysis and environmental dataset) potentially under selection through correlations with environmental variation. Principal components analysis of different outlier SNPs and environmental datasets revealed similar northern and southern clusters, with significant associations between the first axes of each (R²_adj = .66–.79). Multivariate redundancy analysis of outlier SNPs and the environmental principal components indicated that environmental factors explained more than 32% of the variance. Similarly, multiple linear regressions and random‐forest analysis identified winter average and minimum ocean temperatures as significant parameters in the link between genetic and environmental variation. This work indicates that oceanographic variation is associated with the observed genomic cline in this species and that seasonal periods of extreme cold may restrict gene flow along a latitudinal gradient in this marine benthic bivalve. Incorporating this finding into management may improve accuracy of management strategies and future predictions.     

The opportunity for canalization and the evolution of genetic networks

There has been a recent revival of interest in how genetic interactions evolve, spurred on by an increase in our knowledge of genetic interactions at the molecular level. Empirical work on genetic networks has revealed a surprising amount of robustness to perturbations, suggesting that robustness is an evolved feature of genetic networks. Here, we derive a general model for the evolution of canalization that can incorporate any form of perturbation. We establish an upper bound to the strength of selection on canalization that is approximately equal to the fitness load in the system. This method makes it possible to compare different forms of perturbation, including genetic, developmental, and environmental effects. In general, load that arises from mutational processes is low because the mutation rate is itself low. Mutation load can create selection for canalization in a small network that can be achieved through dominance evolution or gene duplication, and in each case selection for canalization is weak at best. In larger genetic networks, selection on genetic canalization can be reasonably strong because larger networks have higher mutational load. Because load induced through migration, segregation, developmental noise, and environmental variance is not mutation limited, each can cause strong selection for canalization.     

The brachymorph mouse and the developmental-genetic basis for canalization and morphological integration

Although it is well known that many mutations influence phenotypic variability as well as the mean, the underlying mechanisms for variability effects are very poorly understood. The brachymorph (bm) phenotype results from an autosomal recessive mutation in the phosphoadenosine-phosphosulfate synthetase 2 gene (Papps2). A major cranial manifestation is a dramatic reduction in the growth of the chondrocranium which results from undersulfation of glycosaminoglycans (GAGs) in the cartilage matrix. We found that this reduction in the growth of the chondrocranium is associated with an altered pattern of craniofacial shape variation, a significant increase in phenotypic variance and a dramatic increase in morphological integration for craniofacial shape. Both effects are largest in the basicranium. The altered variation pattern indicates that the mutation produces developmental influences on shape that are not present in the wildtype. As the mutation dramatically reduces sulfation of GAGs, we infer that this influence is variation among individuals in the degree of sulfation, or variable expressivity of the mutation. This variation may be because of genetic variation at other loci that influence sulfation, environmental effects, or intrinsic effects. We infer that chondrocranial development exhibits greater sensitivity to variation in the sulfation of chondroitin sulfate when the degree of sulfation is low. At normal levels, sulfation probably contributes minimally to phenotypic variation. This case illustrates canalization in a particular developmental-genetic context.     

Morphological intergration between development compartments in the Drosophila wing

Abstract. Developmental integration is the covariation among morphological structures due to connections between the developmental processes that built them. Here we use the methods of geometric morphometrics to study integration in the wing of Drosophila melanogaster . In particular, we focus on the hypothesis that the anterior and posterior wing compartments are separate developmental units that vary independently. We measured both variation among genetically diverse individuals and random differences between body sides of single individuals (fluctuating asymmetry, FA). For both of these sources of variation, the patterns of variation identified by principal component analyses all involved landmarks in both the anterior and posterior compartments simultaneously. Analyses focusing exclusively on the covariation between the anterior and posterior compartments, by the partial least-squares method, revealed pervasive integration of the two compartments, for both individual variation and FA. These analyses clearly indicate that the anterior and posterior compartments are not separate units of variation, but that the covariation between compartments is sufficient to account for nearly all the variation throughout the entire wing. We conclude that variation among individuals as well as the developmental perturbations responsible for FA generate shape variation primarily through developmental processes that are integrated across both compartments. In contrast, much less of the shape variation in our sample can be attributed to the localized processes that establish the identity of particular wing veins.     

Developmental associations between traits: covariance and beyond

Statistical associations between phenotypic traits often result from shared developmental processes and include both covariation between the trait values and more complex associations between higher moments of the joint distribution of traits. In this article, an analytical technique for calculating the covariance between traits is presented on the basis of (1). the distribution of underlying genetic and environmental variation that jointly influences the traits and (2). the mechanics of how these underlying factors influence the development of each trait. It is shown that epistasis can produce patterns of covariation between traits that are not seen in additive models. Applying this approach to a trait in parents and the same trait in their offspring allows us to study the consequences of epistasis for the evolution of additive genetic variance and heritability. This analysis is then extended to the study of more complicated associations between traits. It is shown that even traits that are not correlated may exhibit developmental associations that influence their joint evolution.     

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.     

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.     

Studying morphological integration and modularity at multiple levels: concepts and analysis

Although most studies on integration and modularity have focused on variation among individuals within populations or species, this is not the only level of variation for which integration and modularity exist. Multiple levels of biological variation originate from distinct sources: genetic variation, phenotypic plasticity resulting from environmental heterogeneity, fluctuating asymmetry from random developmental variation and, at the interpopulation or interspecific levels, evolutionary change. The processes that produce variation at all these levels can impart integration or modularity on the covariance structure among morphological traits. In turn, studies of the patterns of integration and modularity can inform about the underlying processes. In particular, the methods of geometric morphometrics offer many advantages for such studies because they can characterize the patterns of morphological variation in great detail and maintain the anatomical context of the structures under study. This paper reviews biological concepts and analytical methods for characterizing patterns of variation and for comparing across levels. Because research comparing patterns across level has only just begun, there are relatively few results, generalizations are difficult and many biological and statistical questions remain unanswered. Nevertheless, it is clear that research using this approach can take advantage of an abundance of new possibilities that are so far largely unexplored.     

Epigenetic Effects on Integration of Limb Lengths in a Mouse Model: Selective Breeding for High Voluntary Locomotor Activity

Mammals exhibit a similar pattern of integration among homologous limb elements, the strength of which is believed to vary in response to selection for functional coordination or similarity. Although integration is hypothesized to primarily reflect the effect of genes intrinsic to limbs, extrinsic genetic or epigenetic factors may also affect the strength of integration through their impact on the magnitude and direction of skeletal variance or covariance. Such factors as neuromuscular coordination or bone-muscle interactions may therefore play a role in both canalization and the structure or magnitude of limb integration. If this were the case, then increased levels of locomotor activity would be predicted to increase canalization and the magnitude of covariation between limbs. To investigate whether postnatal activity levels can have a significant effect on variance within or covariance among homologous limb elements, we compared four groups of male mice from a long-term selective breeding experiment: (1) mice from lines bred for increased voluntary activity on running wheels and allowed free access to a wheel for 8 weeks beginning at weaning (“active”), (2) selected mice that did not have wheel access (“sedentary”), (3) active mice from non-selected control lines, and (4) sedentary control mice. Mice from selected lines that had wheel access ran significantly more than control-line mice. However, when controlled for activity, linetype, and body mass, results indicate few significant differences in means, variance, or covariation structure, and no significant differences in integration between limbs, suggesting that postnatal activity levels do not significantly affect canalization or integration of limb lengths. A possible explanation for this result is that whereas baseline levels of postnatal activity may help to maintain patterns of variance and integration, increased levels of activity do not further increase these measures. Investigations into disrupted epigenetic processes (e.g., via models in which neuromuscular coordination is impaired) are required to further test hypotheses about how canalization or integration of limb variation is affected by epigenetic factors.     

QUANTITATIVE GENETICS OF SHAPE IN CRICKET WINGS: DEVELOPMENTAL INTEGRATION IN A FUNCTIONAL STRUCTURE

The role of developmental and genetic integration for evolution is contentious. One hypothesis states that integration acts as a constraint on evolution, whereas an alternative is that developmental and genetic systems evolve to match the functional modularity of organisms. This study examined a morphological structure, the cricket wing, where developmental and functional modules are discordant, making it possible to distinguish the two alternatives. Wing shape was characterized with geometric morphometrics, quantitative genetic information was extracted using a full‐sibling breeding design, and patterns of developmental integration were inferred from fluctuating asymmetry of wing shape. The patterns of genetic, phenotypic, and developmental integration were clearly similar, but not identical. Heritabilities for different shape variables varied widely, but no shape variables were devoid of genetic variation. Simulated selection for specific shape changes produced predicted responses with marked deflections due to the genetic covariance structure. Three hypotheses of modularity according to the wing structures involved in sound production were inconsistent with the genetic, phenotypic, or developmental covariance structure. Instead, there appears to be strong integration throughout the wing. The hypothesis that genetic and developmental integration evolve to match functional modularity can therefore be rejected for this example.     

Multilevel assessment of the Lacertid lizard cranial modularity

Different factors and processes that produce phenotypic variation at the individual, population, or interspecific level can influence or alter the covariance structure among morphological traits. Therefore, studies of the patterns of integration and modularity at multiple levels—static, ontogenetic, and evolutionary, can provide invaluable data on underlying factors and processes that structured morphological variation, directed, or constrained evolutionary changes. Our dataset, consisting of cranium shape data for 14 lizard species from the family Lacertidae, with substantial samples of hatchlings and adults along with their inferred evolutionary relationships, enabled us to assess modularity and morphological integration at all three levels. Five, not mutually exclusive modularity hypotheses of lizard cranium, were tested, and the effects of allometry on intensity and the pattern of integration and modularity were estimated. We used geometric morphometrics to extract symmetric and asymmetric, as well as allometric and nonallometric, components of shape variation. At the static level, firm confirmation of cranial modularity was found for hypotheses which separate anterior and posterior functional compartments of the skull. At the ontogenetic level, two alternative hypotheses (the “anteroposterior” and “neurodermatocranial” hypotheses) of ventral cranial modularity were confirmed. At the evolutionary level, the “neurodermatocranial” hypothesis was confirmed for the ventral cranium, which is in accordance with the pattern observed at the ontogenetic level. The observed pattern of static modularity could be driven by functional demands and can be regarded as adaptive. Ontogenetic modularity and evolutionary modularity show the same developmental origin, indicating conservatism of modularity patterns driven by developmental constraints.