Chromosomal rearrangements,phenotypic variation and modularity: a case study from a contact zone between house mouse Robertsonian races in Central Italy

The Western European house mouse, Mus musculus domesticus, is well‐known for the high frequency of Robertsonian fusions that have rapidly produced more than 50 karyotipic races, making it an ideal model for studying the mechanisms of chromosomal speciation. The mouse mandible is one of the traits studied most intensively to investigate the effect of Robertsonian fusions on phenotypic variation within and between populations. This complex bone structure has also been widely used to study the level of integration between different morphogenetic units. Here, with the aim of testing the effect of different karyotypic assets on the morphology of the mouse mandible and on its level of modularity, we performed morphometric analyses of mice from a contact area between two highly metacentric races in Central Italy. We found no difference in size, while the mandible shape was found to be different between the two Robertsonian races, even after accounting for the genetic relationships among individuals and geographic proximity. Our results support the existence of two modules that indicate a certain degree of evolutionary independence, but no difference in the strength of modularity between chromosomal races. Moreover, the ascending ramus showed more pronounced interpopulation/race phenotypic differences than the alveolar region, an effect that could be associated to their different polygenic architecture. This study suggests that chromosomal rearrangements play a role in the house mouse phenotypic divergence, and that the two modules of the mouse mandible are differentially affected by environmental factors and genetic makeup.     

Integration and modularity of quantitative trait locus effects on geometric shape in the mouse mandible

The mouse mandible has long served as a model system for complex morphological structures. Here we use new methodology based on geometric morphometrics to test the hypothesis that the mandible consists of two main modules, the alveolar region and the ascending ramus, and that this modularity is reflected in the effects of quantitative trait loci (QTL). The shape of each mandible was analyzed by the positions of 16 morphological landmarks and these data were analyzed using Procrustes analysis. Interval mapping in the F(2) generation from intercrosses of the LG/J and SM/J strains revealed 33 QTL affecting mandible shape. The QTL effects corresponded to a variety of shape changes, but ordination or a parametric bootstrap test of clustering did not reveal any distinct groups of QTL that would affect primarily one module or the other. The correlations of landmark positions between the two modules tended to be lower than the correlations between arbitrary subsets of landmarks, indicating that the modules were relatively independent of each other and confirming the hypothesized location of the boundary between them. While these results are in agreement with the hypothesis of modularity, they also underscore that modularity is a question of the relative degrees to which QTL contribute to different traits, rather than a question of discrete sets of QTL contributing to discrete sets of traits.     

Evolutionary modularity of the mouse mandible: dissecting the effect of chromosomal reorganizations and isolation by distance in a Robertsonian system of Mus musculus domesticus

The mouse mandible consists of several morphogenetic units that are usually grouped into two main modules: the alveolar region and the ascending ramus. The genetic/ontogenetic modularity of the two regions implies that they might evolve independently to some extent. In particular, evolutionary modularity in quantitative traits could arise during chromosomal speciation due to lower gene flow in rearranged chromosomes. With the aim of uncovering the autonomous evolution of the mandible modules, the form variation of each of them was assessed in the house mouse Robertsonian system from Barcelona, in which chromosomal variation and geographical distance may act as isolation factors. The association between these factors and morphological changes was analysed to determine their contribution to the differentiation of each module. Although size changes in the two modules were highly correlated, shape changes were not, and their association with karyotype differences, but not geographical distance, was dependent on the module. The results support the existence of two evolutionary modules and highlight the importance of size in morphological integration of the mandible. They also suggest that geographical distance and chromosomal reorganizations reduce gene flow between karyotypically divergent populations, but although geographical distance represents a global barrier to gene flow, the isolation produced by a set of chromosomal reorganizations only affects particular modules, probably depending on the number and location of loci with effects on a particular morphological region.     

Morphological variation in house mice from the Robertsonian polymorphism area of Barcelona

Morphometric variation in the Robertsonian polymorphism zone of Barcelona of Mus musculus domesticus was studied by geometric morphometrics. This system is characterized by populations of reduced diploid number (2 n  = 27–39) surrounded by standard populations (2 n  = 40). We investigated the morphological variation in mice from this area, as well as the effect of geographical distance and karyotype on this variation. We also investigated the degree of co-variation between the two functional units of the mandible to explore the origin of this system (primary intergradation or secondary contact). The size and shape of the cranium, mandible and scapula were analysed for 226 specimens grouped by population, chromosome number and structural heterozygosity. Size was estimated as the centroid size, and shape was estimated after Procrustes superimposition. No significant differences in size between populations or chromosomal groups were detected. Diploid number, structural heterozygosity and local geographical isolation contributed to the differentiation in shape. Morphological differentiation between standard mice and Robertsonian specimens was observed, suggesting genetic isolation between these groups. Co-variation between the ascending ramus and alveolar region of the mandible was quantified by the trace correlation between landmark subsets of these modules. The trace values showed an ascending trend, correlated with the distance from the centre of the polymorphism area, a pattern consistent with a primary intergradation scenario.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 97 , 555–570.     

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.     

Morphological integration and modularity in the hyperkinetic feeding system of aquatic‐foraging snakes

The kinetic skull is a key innovation that allowed snakes to capture, manipulate, and swallow prey exclusively using their heads using the coordinated movement of eight bones. Despite these unique feeding behaviors, patterns of evolutionary integration and modularity within the feeding bones of snakes in a phylogenetic framework have yet to be addressed. Here, we use a dataset of 60 μCT‐scanned skulls and high‐density geometric morphometric methods to address the origin and patterns of variation and integration in the feeding bones of aquatic‐foraging snakes. By comparing alternate superimposition protocols allowing us to analyze the entire kinetic feeding system simultaneously, we find that the feeding bones are highly integrated, driven predominantly by functional selective pressures. The most supported pattern of modularity contains four modules, each associated with distinct functional roles: the mandible, the palatopterygoid arch, the maxilla, and the suspensorium. Further, the morphological disparity of each bone is not linked to its magnitude of integration, indicating that integration within the feeding system does not strongly constrain morphological evolution, and that adequate biomechanical solutions to a wide range of feeding ecologies and behaviors are readily evolvable within the constraint due to integration in the snake feeding system.     

Modularity of a Cambrian ptychoparioid trilobite cranidium

SUMMARY Modularity of the cranidium of Crassifimbra? metalaspis, a Cambrian ptychoparioid trilobite, is investigated using landmark‐based geometric morphometric methods to gain insight into the integration among morphogenetic processes responsible for shaping the head of an ancient arthropod. Of particular interest is the extent to which the structure of phenotypic integration was governed by direct interactions among developmental pathways, because these interactions may generate long‐term constraints on evolutionary innovation. A modified two‐way ANOVA decomposes cranidial shape variation into components representing symmetric variation among individuals and fluctuating asymmetry (FA). The structure of integration of each of these components is inferred from correlated deviations in shape among nine partitions of the cranidium. Significant correlation among partitions in FA indicates direct interactions among their respective developmental pathways. An a priori hypothesis that modularity was determined by functional association among partitions is not well supported by the among‐partition correlation structure for either component of variation. Instead, exploratory analyses reveal that phenotypic integration was strongly influenced by spatially localized morphogenetic controls. Comparison of the structures of the Individuals and FA components of variation reveals that the two share relatively few commonalities: the structure of phenotypic integration was only weakly influenced by direct interactions. The large contribution of parallel variation to phenotypic integration suggests that modularity was unlikely to have imposed a long‐term constraint on evolutionary innovation in these early trilobites.     

Invasive house mice facing a changing environment on the Sub‐Antarctic Guillou Island (Kerguelen Archipelago)

Adaptation to new environments is a key feature in evolution promoting divergence in morphological structures under selection. The house mouse (Mus musculus domesticus) introduced on the Sub‐Antarctic Guillou Island (Kerguelen Archipelago) had and still has to face environmental conditions that likely shaped the pattern and pace of its insular evolution. Since mouse arrival on the island, probably not more than two centuries ago, ecological conditions dramatically differed from those available to their Western European commensal source populations. In addition, over the last two decades, the plant and animal communities of Guillou Island were considerably modified by the eradication of rabbits, the effects of climate change and the spread of invasive species detrimental to native communities. Under such a changing habitat, the mouse response was investigated using a morphometric quantification of mandible and molar tooth, two morphological structures related to food processing. A marked differentiation of the insular mice compared with their relatives from Western Europe was documented for both mandibles and molar shapes. Moreover, these shapes changed through the 16 years of the record, in agreement with expectations of drift for the molar, but more than expected by chance for the mandible. These results suggest that mice responded to the recent changes in food resources, possibly with a part of plastic variation for the mandible prone to bone remodelling. This pattern exemplifies the intricate interplay of evolution, ecology and plasticity that is a probable key of the success of such an invasive rodent facing pronounced shifts in food resources exploitation under a changing environment.     

Genome‐wide association analyses reveal polygenic genomic architecture underlying divergent shell morphology in Spanish Littorina saxatilis ecotypes

Gene flow between diverging populations experiencing dissimilar ecological conditions can theoretically constrain adaptive evolution. To minimize the effect of gene flow, alleles underlying traits essential for local adaptation are predicted to be located in linked genome regions with reduced recombination. Local reduction in gene flow caused by selection is expected to produce elevated divergence in these regions. The highly divergent crab‐adapted and wave‐adapted ecotypes of the marine snail Littorina saxatilis present a model system to test these predictions. We used genome‐wide association (GWA) analysis of geometric morphometric shell traits associated with microgeographic divergence between the two L. saxatilis ecotypes within three separate sampling sites. A total of 477 snails that had individual geometric morphometric data and individual genotypes at 4,066 single nucleotide polymorphisms (SNPs) were analyzed using GWA methods that corrected for population structure among the three sites. This approach allowed dissection of the genomic architecture of shell shape divergence between ecotypes across a wide geographic range, spanning two glacial lineages. GWA revealed 216 quantitative trait loci (QTL) with shell size or shape differences between ecotypes, with most loci explaining a small proportion of phenotypic variation. We found that QTL were evenly distributed across 17 linkage groups, and exhibited elevated interchromosomal linkage, suggesting a genome‐wide response to divergent selection on shell shape between the two ecotypes. Shell shape trait‐associated loci showed partial overlap with previously identified outlier loci under divergent selection between the two ecotypes, supporting the hypothesis of diversifying selection on these genomic regions. These results suggest that divergence in shell shape between the crab‐adapted and wave‐adapted ecotypes is produced predominantly by a polygenic genomic architecture with positive linkage disequilibrium among loci of small effect.     

Shape variation in the mole dentary (Talpidae: Mammalia)

The clade Talpidae consists of specialized fossorial forms, shrew‐like moles and semi‐aquatic desmans. As with all higher jawed vertebrates, different functional, phylogenetic and developmental constraints act on different parts of dentary influencing its shape. In order to determine whether morphological variation in the dentary was unified or dispersed into an integrated complex of structural–functional components, a morphometric analysis of the mole dentary was undertaken. The dentary was subdivided into component parts – horizonal ramus; coronoid, condylar, angular processes of the ascending ramus – and outline‐based geometric morphometric methods used to quantify, compare and contrast modes of shape variation within the clade. These were successful in revealing subtle differences and aspects of shape important in distinguishing between mole genera. Closer examination of shape variation within the two fully fossorial mole clades (Talpini and Scalopini) revealed several similarities in ascending ramus shapes between genera from each clade. For example, the broad, truncated appearance of the coronoid process in the talpine genera Talpa and Parascalops was shared with the scalopine genus Scapanus. Also, the more slender, hook‐shaped coronoid process of Euroscaptor and Parascaptor (Talpini) closely resembles that of Scalopus (Scalopini). Interestingly, subspecies (one from each clade) more closely resembled genera other than their own in coronoid process shape. Important distinctions in horizontal ramus shape were found to exist between the two clades, such as the extent of curvature of the ventral margin and relative depth of the horizontal ramus. Results show shape variation in this region is correlated with dental formulae and the relative sizes of the teeth. The taxonomically important dentition differences characteristic of mammals are also reflected in the horizontal ramus results. Moreover, these results suggest size may be affecting shape and the extent of variation in, for example, the coronoid and condylar processes between the semi‐aquatic moles Desmana and Galemys. It is likely that the effects of morphological integration seen at this level of analysis – covariation between shapes of dentary components – may exist because interacting traits are evolving together. Horizontal ramus and coronoid process shape, for example, are similar across Scapanus and Parascalops, but both these shapes have diverged in Scalopus. © 2008 Trustees of the Natural History Museum (London). Journal compilation © 2008 The Linnean Society of London, Zoological Journal of the Linnean Society, 2008, 153 , 187–211.     

PLEIOTROPIC EFFECTS OF INDIVIDUAL GENE LOCI ON MANDIBULAR MORPHOLOGY

The genotypic basis of morphological variation is largely unknown. In this study we examine patterns of pleiotropic effects on mandibular morphology at individual gene loci to determine whether the pleiotropic effects of individual genes are restricted to functionally and developmentally related traits. Mandibular measurements were obtained from 480 mice from the F₂ generation of an intercross between the LG/J and SM/J mouse strains. DNA was also extracted from these animals, and 76 microsatellite loci covering the autosomes were scored. Interval mapping was used to detect chromosomal locations with significant effects on various mandibular measurements. Sets of traits mapping to a common chromosomal region were considered as being affected by a single quantitative trait locus (QTL) for mandibular morphology. Thirty-seven such chromosomal regions were identified spread throughout the autosomes. Gene effects were small to moderate with the allele derived from the LG/J strain typically leading to larger size. When dominance was present, the LG/J allele was typically dominant to the SM/J allele. Most loci affected restricted functional and developmental regions of the mandible. Of the 26 chromosomal regions affecting more than two traits, 50% affect the muscular processes of the ascending ramus, 27% affect the alveolar processes carrying the teeth, and 23% affect the whole mandible. Four additional locations affecting two traits had effects significantly associated with alveolar regions. Pleiotropic effects are typically restricted to morphologically integrated complexes.     

A MODULAR FRAMEWORK CHARACTERIZES MICRO‐ AND MACROEVOLUTION OF OLD WORLD MONKEY DENTITIONS

The study of modularity can provide a foundation for integrating development into studies of phenotypic evolution. The dentition is an ideal phenotype for this as it is developmentally relatively simple, adaptively highly significant, and evolutionarily tractable through the fossil record. Here, we use phenotypic variation in the dentition to test a hypothesis about genetic modularity. Quantitative genetic analysis of size variation in the baboon dentition indicates a genetic modular framework corresponding to tooth type categories. We analyzed covariation within the dentitions of six species of Old World monkeys (OWMs) to assess the macroevolutionary extent of this framework: first by estimating variance–covariance matrices of linear tooth size, and second by performing a geometric morphometric (GM) analysis of tooth row shape. For both size and shape, we observe across OWMs a framework of anterior and postcanine modules, as well as submodularity between the molars and premolars. Our results of modularity by tooth type suggest that adult variation in the OWM dentition is influenced by early developmental processes such as odontogenesis and jaw patterning. This study presents a comparison of genotypic modules to phenotypic modules, which can be used to better understand their action across evolutionary time scales.     

Is the genotype-phenotype map modular? A statistical approach using mouse quantitative trait loci data.

Various theories about the evolution of complex characters make predictions about the statistical distribution of genetic effects on phenotypic characters, also called the genotype-phenotype map. With the advent of QTL technology, data about these distributions are becoming available. In this article, we propose simple tests for the prediction that functionally integrated characters have a modular genotype-phenotype map. The test is applied to QTL data on the mouse mandible. The results provide statistical support for the notion that the ascending ramus region of the mandible is modularized. A data set comprising the effects of QTL on a more extensive portion of the phenotype is required to determine if the alveolar region of the mandible is also modularized.     

Semi‐supervised determination of pseudocryptic morphotypes using observer‐free characterizations of anatomical alignment and shape

Accurate, quantitative characterization of complex shapes is recognized as a key methodological challenge in biology. Recent development of automated three‐dimensional geometric morphometric protocols (auto3dgm) provides a promising set of tools to help address this challenge. While auto3dgm has been shown to be useful in characterizing variation across clades of morphologically very distinct mammals, it has not been adequately tested in more problematic cases where pseudolandmark placement error potentially confounds interpretation of true shape variation. Here, we tested the sensitivity of auto3dgm to the degree of variation and various parameterization settings using a simulation and three microCT datasets that characterize mammal tooth crown morphology as biological examples. The microCT datasets vary in degree of apparent morphological differentiation, with two that include grossly similar morphospecies and one that includes two laboratory strains of a single species. Resulting alignments are highly sensitive to the number of pseudolandmarks used to quantify shapes. The degree to which the surfaces were downsampled and the apparent degree of morphological differentiation across the dataset also influenced alignment repeatability. We show that previous critiques of auto3dgm were based on poorly parameterized alignments and suggest that sample‐specific sensitivity analyses should be added to any research protocol including auto3dgm. Auto3dgm is a useful tool for studying samples when pseudolandmark placement error is small relative to the true differences between specimens. This method therefore represents a promising avenue forward in morphometric studies at a wide range of scales, from samples that differ by a single genetic locus to samples that represent multiple phylogenetically diverse clades.     

Does nasal echolocation influence the modularity of the mammal skull?

In vertebrates, changes in cranial modularity can evolve rapidly in response to selection. However, mammals have apparently maintained their pattern of cranial integration throughout their evolutionary history and across tremendous morphological and ecological diversity. Here, we use phylogenetic, geometric morphometric and comparative analyses to test the hypothesis that the modularity of the mammalian skull has been remodelled in rhinolophid bats due to the novel and critical function of the nasal cavity in echolocation. We predicted that nasal echolocation has resulted in the evolution of a third cranial module, the ‘nasal dome’, in addition to the braincase and rostrum modules, which are conserved across mammals. We also test for similarities in the evolution of skull shape in relation to habitat across rhinolophids. We find that, despite broad variation in the shape of the nasal dome, the integration of the rhinolophid skull is highly consistent with conserved patterns of modularity found in other mammals. Across their broad geographical distribution, cranial shape in rhinolophids follows two major divisions that could reflect adaptations to dietary and environmental differences in African versus South Asian distributions. Our results highlight the potential of a relatively simple modular template to generate broad morphological and functional variation in mammals.     

Chondrocranial development in larval Rana sylvatica (Anura: Ranidae): morphometric analysis of cranial allometry and ontogenetic shape change

This study provides baseline quantitative data on the morphological development of the chondrocranium in a larval anuran. Both linear and geometric morphometric methods are used to quantitatively analyze size-related shape change in a complete developmental series of larvae of the wood frog, Rana sylvatica. The null hypothesis of isometry was rejected in all geometric morphometric and most linear morphometric analyses. Reduced major axis regressions of 11 linear chondrocranial measurements on size indicate a mixture of allometric and isometric scaling. Measurements in the otic and oral regions tend to scale with negative allometry and those associated with the palatoquadrate and muscular process scale with isometry or positive allometry. Geometric morphometric analyses, based on a set of 11 chondrocranial landmarks, include linear regression of relative warp scores and multivariate regression of partial warp scores and uniform components on log centroid size. Body size explains about one-quarter to one-third of the total shape variation found in the sample. Areas of regional shape transformation (e.g., palatoquadrate, otic region, trabecular horns) are identified by thin-plate spline deformation grids and are concordant with linear morphometric results. Thus, the anuran chondrocranium is not a static structure during premetamorphic stages and allometric patterns generally follow scaling predictions for tetrapod cranial development. Potential implications regarding larval functional morphology, cranial development, and chondrocranial evolution in anurans are discussed.     

Does divergence from normal patterns of integration increase as chromosomal fusions increase in number? A test on a house mouse hybrid zone

Chromosomal evolution is widely considered an important driver of speciation because it can promote the establishment of reproductive barriers. Karyotypic reorganization is also expected to affect the mean phenotype, as well as its development and patterns of phenotypic integration, through processes such as variation in genetic linkage between quantitative trait loci or between regulatory regions and their targets. Here we explore the relationship between chromosomal evolution and phenotypic integration by analyzing a well-known house mouse parapatric contact zone between a highly derived Robertsonian (Rb) race (2n = 22) and populations with standard karyotype (2n = 40). Populations with hybrid karyotypes are scattered throughout the hybrid zone connecting the two parental races. Using mandible shape data and geometric morphometrics, we test the hypothesis that patterns of integration progressively diverge from the “normal” integration pattern observed in the standard race as they accumulate Rb fusions. We find that the main pattern of integration observed between the posterior and anterior part of the mandible can be largely attributed to allometry. We find no support for a gradual increase in divergence from normal patterns of integration as fusions accumulate. Surprisingly, however, we find that the derived Rb race (2n = 22) has a distinct allometric trajectory compared with the standard race. Our results suggest that either individual fusions disproportionately affect patterns of integration or that there are mechanisms which “purge” extreme variants in hybrids (e.g. reduced fitness of hybrid shape).     

Heritability of male mandible length in the stag beetle Cyclommatus metallifer

Numerous coleopteran species express male‐specific “weapon traits” that often show size variations among males, even within a single population. Many empirical studies have demonstrated that environmental conditions during development affect absolute weapon size. However, relatively few studies in horned beetles support the hypothesis that the relationship between weapon size and body size, also referred to as a “scaling relationship” or “static allometry”, is largely determined by genetic factors. In this study, the heritability of absolute mandible length and static allometry between mandible length and body size were estimated in the stag beetle Cyclommatus metallifer. While no significant heritable variation was observed in absolute mandible length, high heritability (h² = 0.57 ± 0.25) was detected in the static allometry between mandible length and body size. This is the first report on the genetic effect on male mandible size in Lucanidae, suggesting that absolute mandible size is largely determined by environmental conditions while the static allometry between weapon size and body size is primarily determined by genetic factors.     

Rapid morpho‐functional changes among insular populations of the greater white‐toothed shrew

Islands are often considered to be natural laboratories where repeated ‘evolutionary experiments’ have taken place. Consequently, islands have been key model systems in our understanding of evolutionary theory. The greater white‐toothed shrew (Crocidura russula) is of interest as it has invaded French Atlantic islands within the last few thousand years and is considered to be morphologically and genetically stable in this area. In this article, we study the shape of the mandible of the greater white‐toothed shrew on four islands and compare it with that of individuals from populations on the mainland to quantify the effects of insularity. The degree of insularity (i.e. island size and distance to the continent) is thought to be linked to differences in ecological characteristics of islands compared with the mainland. We used geometric morphometric analyses to quantify differences in size and shape between populations and employed a simple biomechanical model to evaluate the potential effects of shape differences on bite force. Specimens from island populations are different from continental populations in shape and mechanical potential of the mandible. Among islands, the mandible shows various shapes that are correlated with both the distance from the coast and island area. The shape differences are located on different parts of the mandible, suggesting different ecological constraints on each island. Moreover, these shapes are linked to the ‘mechanical potential’, which is markedly different between islands. Mechanical potential has been suggested to evolve in response to prey size and or mechanical properties. In conclusion, our results show that, in spite of the relatively recent colonization of the Atlantic Islands, the mandible of C. russula possesses a distinct shape. Moreover, the shape differs among islands and is probably linked to the consumption of different prey. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, ?? , ??–??.     

The Influence of Modularity on Cranial Morphological Disparity in Carnivora and Primates (Mammalia)

Background

Although variation provides the raw material for natural selection and evolution, few empirical data exist about the factors controlling morphological variation. Because developmental constraints on variation are expected to act by influencing trait correlations, studies of modularity offer promising approaches that quantify and summarize patterns of trait relationships. Modules, highly-correlated and semi-autonomous sets of traits, are observed at many levels of biological organization, from genes to colonies. The evolutionary significance of modularity is considerable, with potential effects including constraining the variation of individual traits, circumventing pleiotropy and canalization, and facilitating the transformation of functional structures. Despite these important consequences, there has been little empirical study of how modularity influences morphological evolution on a macroevolutionary scale. Here, we conduct the first morphometric analysis of modularity and disparity in two clades of placental mammals, Primates and Carnivora, and test if trait integration within modules constrains or facilitates morphological evolution.

Principal Findings

We used both randomization methods and direct comparisons of landmark variance to compare disparity in the six cranial modules identified in previous studies. The cranial base, a highly-integrated module, showed significantly low disparity in Primates and low landmark variance in both Primates and Carnivora. The vault, zygomatic-pterygoid and orbit modules, characterized by low trait integration, displayed significantly high disparity within Carnivora. 14 of 24 results from analyses of disparity show no significant relationship between module integration and morphological disparity. Of the ten significant or marginally significant results, eight support the hypothesis that integration within modules constrains morphological evolution in the placental skull. Only the molar module, a highly-integrated and functionally important module, showed significantly high disparity in Carnivora, in support of the facilitation hypothesis.

Conclusions

This analysis of within-module disparity suggested that strong integration of traits had little influence on morphological evolution over large time scales. However, where significant results were found, the primary effect of strong integration of traits was to constrain morphological variation. Thus, within Primates and Carnivora, there was some support for the hypothesis that integration of traits within cranial modules limits morphological evolution, presumably by limiting the variation of individual traits.