Developmental plasticity, morphological variation and evolvability: a multilevel analysis of morphometric integration in the shape of compound leaves

The structure of compound leaves provides flexibility for morphological change by variation in the shapes, sizes and arrangement of leaflets. Here, we conduct a multilevel analysis of shape variation in compound leaves to explore the developmental plasticity and evolutionary potential that are the basis of diversification in leaf shape. We use the methods of geometric morphometrics to study the shapes of individual leaflets and whole leaves in 20 taxa of Potentilla (sensu lato). A newly developed test based on the bootstrap approach suggests that uncertainty in the molecular phylogeny precludes firm conclusions whether there is a phylogenetic signal in the data on leaf shape. For variation among taxa, variation within taxa, as well as fluctuating asymmetry, there is evidence of strong morphological integration. The patterns of variation are similar across all three levels, suggesting that integration within taxa may act as a constraint on evolutionary change.     

Developmental constraints revealed by co-variation within and among molar rows in two murine rodents

SUMMARY Morphological integration corresponds to interdependency between characters that can arise from several causes. Proximal causes of integration include that different phenotypic features may share common genetic sets and/or interact during their development. Ultimate causes may be the prolonged effect of selection favoring integration of functionally interacting characters, achieved by the molding of these proximal causes. Strong and direct interactions among successive teeth of a molar row are predicted by genetic and developmental evidences. Functional constraints related to occlusion, however, should have selected more strongly for a morphological integration of occluding teeth and a corresponding evolution of the underlying developmental and genetic pathways. To investigate how these predictions match the patterns of phenotypic integration, we studied the co‐variation among the six molars of the murine molar row, focusing on two populations of house mice (Mus musculus domesticus) and wood mice (Apodemus sylvaticus). The size and shape of the three upper and lower molars were quantified and compared. Our results evidenced similar patterns in both species, size being more integrated than shape among all the teeth, and both size and shape co‐varying strongly between adjacent teeth, but also between occluding teeth. Strong co‐variation within each molar row is in agreement with developmental models showing a cascade influence of the first molar on the subsequent molars. In contrast, the strong co‐variation between molars of the occluding tooth rows confirms that functional constraints molded patterns of integration and probably the underlying developmental pathways despite the low level of direct developmental interactions occurring among molar rows. These patterns of co‐variation are furthermore conserved between the house mouse and the wood mouse that diverged >10 Ma, suggesting that they may constitute long‐running constraints to the diversification of the murine rodent dentition.     

Phylogenetic signal,function and integration in the subunits of the carnivoran mandible

Complex phenotypes could be interpreted as the result of functional integration between identifiable subunits. Common developmental or ecological factors may favour macroevolutionary morphological integration so that functional subunits also covary above the species level. We investigate shape variation and functional integration in two subunits of the mammalian mandible: the corpus and the ramus in a subset of extant terrestrial Carnivora using geometric morphometric and comparative methods. More specifically, we test if corpus and ramus shape exhibit similar degree of homoplasy and if these traits covary above species level. Additionally, broad functional categorisations (predaceous and non predaceous) are investigated to test if differences in morphological variation and integration at macroevolutionary scale occur. Principal components of shape data show a significant phylogenetic signal in both mandibular subunits, with the ramus exhibiting a higher degree of homoplasy than the corpus. Functional groups (predators and non-predators) are significantly distinct in corpus shape, while in the ramus significance emerges only after removing the phylogenetic signal. Partial Least Square shows that mandibular corpus and ramus region covaries above species level even if this trend is not supported when employing comparative methods. Only in a subset of predaceous species covariation still hold. We conclude that mandibular subunits of Carnivora differ considerably in shape among predaceous and non-predaceous species because of the adaptive selection pressure imposed by catching and hold of live prey. This selective process also favoured macroevolutionary integration in predaceous carnivorans.     

Phylogenetic and environmental components of morphological variation: skull, mandible, and molar shape in marmots (Marmota, Rodentia)

The phenotype is a product of its phylogenetic history and its recent adaptation to local environments, but the relative importance of the two factors is controversial. We assessed the effects of diet, habitat, elevation, temperature, precipitation, body size, and mtDNA genetic divergence on shape variation in skulls, mandibles, and molars, structures that differ in their genetic and functional control. We asked whether these structures have adapted to environment to the same extent and whether they retain the same amount of phylogenetic signal. We studied these traits in intra- and interspecific populations of Eurasian marmots whose last common ancestor lived 2-5 million years ago. Path Analysis revealed that body size explained 10% of variation in skulls, 7% in mandibles, and 15% in molars. Local vegetation explained 7% of variation in skulls, 11% in mandibles, and 12% in molars. Dietary category explained 25% of variation in skulls, 11% in mandibles, and 9% in molars. Cyt b mtDNA divergence (phylogeny) explained 15% of variation in skulls, 7% in mandibles, and 5% in molars. Despite the percentages of phylogenetic variance, maximum-likelihood trees based on molar and skull shape recovered most phylogenetic groupings correctly, but mandible shape did not. The good performance of molars and skulls was probably due to different factors. Skulls are genetically and functionally more complicated than teeth, and they had more mathematically independent components of variation (5-6-in skulls compared to 3-in molars). The high proportion of diet-related variance was not enough to mask the phylogenetic signal. Molars had fewer independent components, but they also have less ecophenotypic variation and evolve more slowly, giving each component a proportionally stronger phylogenetic signal. Molars require larger samples for each operational taxonomic unit than the other structures because the proportion of within-taxon to between-taxon variation was higher. Good phylogenetic signal in quantitative skeletal morphology is likely to be found only when the taxa have a common ancestry no older than hundreds of thousands or millions of years (1% to 10% mtDNA divergence)--under these conditions skulls and molars provide stronger signal than mandibles.     

A geometric morphometric analysis of the shape of the first upper molar in mice of the genus Mus (Muridae, Rodentia)

Phenotypic variation in the shape of the first upper molar among 595 mice, representing nine extant and three extinct taxa of the genus Mus , was studied with thin-plate spline analysis. The reliability of classification of individual specimens into known groups based on their molars varied from 75 to 100%, depending on group and method used. Including 13 sliding semilandmarks to the analysis improved the detection of different kinds of size and shape variation as well as visualization of shape differences between studied groups. Correlation between phylogenetic and morphometric distances suggested about 80% contribution of phylogenetic inertia to the molar shape variation; moreover, the importance of localized versus global shape changes was similar in the detection of phylogenetic signals. Finally, shape changes along individual evolutionary lineages were revealed, suggesting a few cases of reversals, convergence and/or retention of ancestral shape. The evolution of mouse molars has thus been driven by random effects of drift together with stabilizing selection and convergence.     

Morphological integration in the hominin dentition: evolutionary, developmental, and functional factors

As the most common and best preserved remains in the fossil record, teeth are central to our understanding of evolution. However, many evolutionary analyses based on dental traits overlook the constraints that limit dental evolution. These constraints are diverse, ranging from developmental interactions between the individual elements of a homologous series (the whole dentition) to functional constraints related to occlusion. This study evaluates morphological integration in the hominin dentition and its effect on dental evolution in an extensive sample of Plio- and Pleistocene hominin teeth using geometric morphometrics and phylogenetic comparative methods. Results reveal that premolars and molars display significant levels of covariation; that integration is stronger in the mandibular dentition than in the maxillary dentition; and that antagonist teeth, especially first molars, are strongly integrated. Results also show an association of morphological integration and evolution. Stasis is observed in elements with strong functional and/or developmental interactions, namely in first molars. Alternatively, directional evolution (and weaker integration) occurs in the elements with marginal roles in occlusion and mastication, probably in response to other direct or indirect selective pressures. This study points to the need to reevaluate hypotheses about hominin evolution based on dental characters, given the complex scenario in which teeth evolve.     

Assessing Trait Covariation and Morphological Integration on Phylogenies Using Evolutionary Covariance Matrices

Morphological integration describes the degree to which sets of organismal traits covary with one another. Morphological covariation may be evaluated at various levels of biological organization, but when characterizing such patterns across species at the macroevolutionary level, phylogeny must be taken into account. We outline an analytical procedure based on the evolutionary covariance matrix that allows species-level patterns of morphological integration among structures defined by sets of traits to be evaluated while accounting for the phylogenetic relationships among taxa, providing a flexible and robust complement to related phylogenetic independent contrasts based approaches. Using computer simulations under a Brownian motion model we show that statistical tests based on the approach display appropriate Type I error rates and high statistical power for detecting known levels of integration, and these trends remain consistent for simulations using different numbers of species, and for simulations that differ in the number of trait dimensions. Thus, our procedure provides a useful means of testing hypotheses of morphological integration in a phylogenetic context. We illustrate the utility of this approach by evaluating evolutionary patterns of morphological integration in head shape for a lineage of Plethodon salamanders, and find significant integration between cranial shape and mandible shape. Finally, computer code written in R for implementing the procedure is provided.     

Functional and phylogenetic implications of molar flare variation in Miocene hominoids

Comparative analyses of molar shape figure prominently in Miocene hominoid evolutionary studies, and incomplete understanding of functional and phylogenetic influences on molar shape variation can have direct consequences for the interpretation of fossil taxa. Molar flare is a shape trait whose polarity, phylogenetic distribution, and functional significance have been sources of contention. To clarify the determinants of molar flare variation in the hominoid radiation, a combination of statistical methods was employed to investigate the effects of diet, phylogeny, and geologic age upon several measures of molar shape, to identify interactions among these factors, and to estimate their relative influence. Classic indices of molar crown shape and cusp relief are highly significantly associated with diet and show no clear phylogenetic or temporal patterning. Correlations with diet are insignificant when phylogenetic effects are controlled, a result which is interpreted as an artifact of the distribution of folivory in the Miocene hominoid radiation. Possession of pronounced molar flare was found to be the primitive condition for Miocene hominoids, but molar flare reduction cannot be considered a crown hominoid synapomorphy. Molar flare is strongly correlated with geologic age but differs significantly among dietary categories when the effects of time are controlled. Among contemporaneous taxa, hard-object feeders consistently show the highest levels of flare. Molar flare reduction is hypothesized to arise from realignment of cusp positions to maximize molar shearing and increase working occlusal surface area, while variation in flare among contemporaneous taxa may be due, at least in part, to enamel thickness variation. The pronounced molar flare of Otavipithecus is interpreted as a primitive retention, although alternative dietary and phylogenetic interpretations cannot be excluded. A dramatic reversal of molar flare reduction in Mio-Pliocene hominins is interpreted as a synapomorphy of the crown hominin clade, thus supporting the hominin status of the Lukeino hominine. The last common ancestor of the Pan-Homo clade is predicted to have possessed relatively non-flaring molars, and implications of this hypothesis for early hominin recognition are discussed.     

Convergent evolution of phenotypic integration and its alignment with morphological diversification in Caribbean Anolis ecomorphs

The adaptive landscape and the G-matrix are keys concepts for understanding how quantitative characters evolve during adaptive radiation. In particular, whether the adaptive landscape can drive convergence of phenotypic integration (i.e., the pattern of phenotypic variation and covariation summarized in the P-matrix) is not well studied. We estimated and compared P for 19 morphological traits in eight species of Caribbean Anolis lizards, finding that similarity in P among species was not correlated with phylogenetic distance. However, greater similarity in P among ecologically similar Anolis species (i.e., the trunk-ground ecomorph) suggests the role of convergent natural selection. Despite this convergence and relatively deep phylogenetic divergence, a large portion of eigenstructure of P is retained among our eight focal species. We also analyzed P as an approximation of G to test for correspondence with the pattern of phenotypic divergence in 21 Caribbean Anolis species. These patterns of covariation were coincident, suggesting that either genetic constraint has influenced the pattern of among-species divergence or, alternatively, that the adaptive landscape has influenced both G and the pattern of phenotypic divergence among species. We provide evidence for convergent evolution of phenotypic integration for one class of Anolis ecomorph, revealing yet another important dimension of evolutionary convergence in this group.     

Morpho morphometrics: Shared ancestry and selection drive the evolution of wing size and shape in Morpho butterflies

Butterfly wings harbor highly diverse phenotypes and are involved in many functions. Wing size and shape result from interactions between adaptive processes, phylogenetic history, and developmental constraints, which are complex to disentangle. Here, we focus on the genus Morpho (Nymphalidae: Satyrinae, 30 species), which presents a high diversity of sizes, shapes, and color patterns. First, we generate a comprehensive molecular phylogeny of these 30 species. Next, using 911 collection specimens, we quantify the variation of wing size and shape across species, to assess the importance of shared ancestry, microhabitat use, and sexual selection in the evolution of the wings. While accounting for phylogenetic and allometric effects, we detect a significant difference in wing shape but not size among microhabitats. Fore and hindwings covary at the individual and species levels, and the covariation differs among microhabitats. However, the microhabitat structure in covariation disappears when phylogenetic relationships are taken into account. Our results demonstrate that microhabitat has driven wing shape evolution, although it has not strongly affected forewing and hindwing integration. We also found that sexual dimorphism of forewing shape and color pattern are coupled, suggesting a common selective force.     

Patterns of morphological integration in the appendicular skeleton of mammalian carnivores

We investigated patterns of evolutionary integration in the appendicular skeleton of mammalian carnivores. The findings are discussed in relation to performance selection in terms of organismal function as a potential mechanism underlying integration. Interspecific shape covariation was quantified by two‐block partial least‐squares (2B‐PLS) analysis of 3D landmark data within a phylogenetic context. Specifically, we compared pairs of anatomically connected bones (within‐limbs) and pairs of both serially homologous and functional equivalent bones (between‐limbs). The statistical results of all the comparisons suggest that the carnivoran appendicular skeleton is highly integrated. Strikingly, the main shape covariation relates to bone robustness in all cases. A bootstrap test was used to compare the degree of integration between specialized cursorial taxa (i.e., those whose forelimbs are primarily involved in locomotion) and noncursorial species (i.e., those whose forelimbs are involved in more functions than their hindlimb) showed that cursors have a more integrated appendicular skeleton than noncursors. The findings demonstrate that natural selection can influence the pattern and degree of morphological integration by increasing the degree of bone shape covariation in parallel to ecological specialization.     

Morphological integration and serial homology: A case study of the cranium and anterior vertebrae in salamanders

Serial homology or the repetition of equivalent developmental units and their derivatives is a phenomenon encountered in a variety of organisms, with the vertebrate axial skeleton as one of the most notable examples. Serially homologous structures can be viewed as an appropriate model system for studying morphological integration and modularity, due to the strong impact of development on their covariation. Here, we explored the pattern of morphological integration of the cranium and the first three serially homologous structures (atlas, first, and second trunk vertebrae) in salamandrid salamanders, using micro-CT scanning and three-dimensional geometric morphometrics. We explored the integration between structures at static and evolutionary levels. Effects of allometry on patterns of modularity were also taken into account. At the static level (within species), we analyzed inter-individual variation in shape to detect functional modules and intra-individual variation to detect developmental modules. Significant integration (based on inter-individual variation) among all structures was detected and allometry is shown to be an important integrating factor. The pattern of intra-individual, asymmetric variation indicates statistically significant developmental integration between the cranium and the atlas and between the first two trunk vertebrae. At the evolutionary level (among species), the cranium, atlas, and trunk vertebrae separate as different modules. Our results show that morphological integration at the evolutionary level coincides with morphological and functional differentiation of the axial skeleton, allowing the more or less independent evolutionary changes of the cranial skeleton and the vertebral column, regardless of the relatively strong integration at the static level. The observed patterns of morphological integration differ across levels, indicating different impacts of developmental and phylogenetic constraints and functional demands.     

Journeys through discrete‐character morphospace: synthesizing phylogeny,tempo, and disparity

Palaeontologists have long employed discrete categorical data to capture morphological variation in fossil species, using the resulting character–taxon matrices to measure evolutionary tempo, infer phylogenies and capture morphological disparity. However, to date these have been seen as separate approaches despite a common goal of understanding morphological evolution over deep time. Here I argue that there are clear advantages to considering these three lines of enquiry in a single space: the phylomorphospace. Conceptually these high‐dimensional spaces capture how a phylogenetic tree explores morphospace and allow us to consider important process questions around evolutionary rates, constraints, convergence and directional trends. Currently the literature contains fundamentally different approaches used to generate such spaces, with no direct comparison between them, despite the differing evolutionary histories they imply. Here I directly compare five different phylomorphospace approaches, three with direct literature equivalents and two that are novel. I use a single empirical case study of coelurosaurian theropod dinosaurs (152 taxa, 853 characters) to show that under many analyses the literature‐derived approaches tend to reflect introduced phylogenetic (rather than the intended morphological) signal. The two novel approaches, which produce limited ancestral state estimates prior to ordination, are able to minimize this phylogenetic signal and thus exhibit more realistic amounts of phylogenetic signal, rate heterogeneity, and convergent evolution.     

Serial homology and the evolution of mammalian limb covariation structure

The tetrapod forelimb and hindlimb are serially homologous structures that share a broad range of developmental pathways responsible for their patterning and outgrowth. Covariation between limbs, which can introduce constraints on the production of variation, is related to the duplication of these developmental factors. Despite this constraint, there is remarkable diversity in limb morphology, with a variety of functional relationships between and within forelimb and hindlimb elements. Here we assess a hierarchical model of limb covariation structure based on shared developmental factors. We also test whether selection for morphologically divergent forelimbs or hindlimbs is associated with reduced covariation between limbs. Our sample includes primates, murines, a carnivoran, and a chiropteran that exhibit varying degrees of forelimb and hindlimb specialization, limb size divergence, and/or phylogenetic relatedness. We analyze the pattern and significance of between-limb morphological covariation with linear distance data collected using standard morphometric techniques and analyzed by matrix correlations, eigenanalysis, and partial correlations. Results support a common limb covariation structure across these taxa and reduced covariation between limbs in nonquadruped species. This result indicates that diversity in limb morphology has evolved without signficant modifications to a common covariation structure but that the higher degree of functional limb divergence in bats and, to some extent, gibbons is associated with weaker integration between limbs. This result supports the hypothesis that limb divergence, particularly selection for increased functional specialization, involves the reduction of developmental factors common to both limbs, thereby reducing covariation.     

Colour, design and reward: phenotypic integration of fleshy fruit displays

The functional or structural linkage among traits [phenotypic integration (PI)] within complex structures can constrain the evolutionary response of individual traits. To analyse whether frugivores with distinct sensory ecology have shaped the patterns of fruit diversification differently, we compared PI values of fleshy fruits that are consumed by birds and mammals. We used phylogenetic comparative analyses of PI among 13 morphological, nutritional and visual fruit traits from 111 Mediterranean plant species. Results showed that morphological traits had higher PI values than nutritional and colour traits. Visual and nutritional traits show positive phylogenetic covariance, while negative covariation occurs between fruits size and nutrients. Importantly, fruits consumed by birds were relatively more integrated than fruits consumed partly or solely by mammals. Hence, we show that major groups of mutualistic frugivores can shape the covariance among some fruit traits differently and thereby influence fruit diversification.     

Morphological innovation and biomechanical diversity in plunge-diving birds

Innovations in foraging behavior can drive morphological diversity by opening up new ways of interacting with the environment, or limit diversity through functional constraints associated with different foraging behaviors. Several classic examples of adaptive radiations in birds show increased variation in ecologically relevant traits. However, these cases primarily focus on geographically narrow adaptive radiations, consider only morphological evolution without a biomechanical approach, or do not investigate tradeoffs with other non-focal traits that might be affected by use of different foraging habitats. Here, we use X-ray microcomputed tomography, biomechanical modeling, and multivariate comparative methods to explore the interplay between foraging behavior and cranial morphology in kingfishers, a global radiation of birds with variable beaks and foraging behaviors, including the archetypal plunge-dive into water. Our results quantify covariation between the shape of the outer keratin covering (rhamphotheca) and the inner skeletal core of the beak, as well as highlight distinct patterns of morphospace occupation for different foraging behaviors and considerable rate variation among these skull regions. We anticipate these findings will have implications for inferring beak shapes in fossil taxa and inform biomimetic design of novel impact-reducing structures.     

Dental topography of platyrrhines and prosimians: Convergence and contrasts

Dental topographic analysis is the quantitative assessment of shape of three‐dimensional models of tooth crowns and component features. Molar topographic curvature, relief, and complexity correlate with aspects of feeding behavior in certain living primates, and have been employed to investigate dietary ecology in extant and extinct primate species. This study investigates whether dental topography correlates with diet among a diverse sample of living platyrrhines, and compares platyrrhine topography with that of prosimians. We sampled 111 lower second molars of 11 platyrrhine genera and 121 of 20 prosimian genera. For each tooth we calculated Dirichlet normal energy (DNE), relief index (RFI), and orientation patch count (OPCR), quantifying surface curvature, relief, and complexity respectively. Shearing ratios and quotients were also measured. Statistical analyses partitioned effects of diet and taxon on topography in platyrrhines alone and relative to prosimians. Discriminant function analyses assessed predictive diet models. Results indicate that platyrrhine dental topography correlates to dietary preference, and platyrrhine‐only predictive models yield high rates of accuracy. The same is true for prosimians. Topographic variance is broadly similar among platyrrhines and prosimians. One exception is that platyrrhines display higher average relief and lower relief variance, possibly related to lower relative molar size and functional links between relief and tooth longevity distinct from curvature or complexity. Explicitly incorporating phylogenetic distance matrices into statistical analyses of the combined platyrrhine‐prosimian sample results in loss of significance of dietary effects for OPCR and SQ, while greatly increasing dietary significance of RFI. Am J Phys Anthropol 153:29–44, 2014. © 2013 Wiley Periodicals, Inc.     

Recent developments in the analysis of comparative data

Comparative methods can be used to test ideas about adaptation by identifying cases of either parallel or convergent evolutionary change across taxa. Phylogenetic relationships must be known or inferred if comparative methods are to separate the cross-taxonomic covariation among traits associated with evolutionary change from that attributable to common ancestry. Only the former can be used to test ideas linking convergent or parallel evolutionary change to some aspect of the environment. The comparative methods that are currently available differ in how they manage the effects brought about by phylogenetic relationships. One method is applicable only to discrete data, and uses cladistic techniques to identify evolutionary events that depart from phylogenetic trends. Techniques for continuous variables attempt to control for phylogenetic effects in a variety of ways. One method examines the taxonomic distribution of variance to identify the taxa within which character variation is small. The method assumes that taxa with small amounts of variation are those in which little evolutionary change has occurred, and thus variation is unlikely to be independent of ancestral trends. Analyses are then concentrated among taxa that show more variation, on the assumption that greater evolutionary change in the character has taken place. Several methods estimate directly the extent to which ancestry can predict the observed variation of a character, and subtract the ancestral effect to reveal variation of phylogeny. Yet another can remove phylogenetic effects if the true phylogeny is known. One class of comparative methods controls for phylogenetic effects by searching for comparative trends within rather than across taxa. With current knowledge of phylogenies, there is a trade-off in the choice of a comparative method: those that control phylogenetic effects with greater certainty are either less applicable to real data, or they make restrictive or untestable assumptions. Those that rely on statistical patterns to infer phylogenetic effects may not control phylogeny as efficiently but are more readily applied to existing data sets.     

Geometric morphometric 3D shape analysis and covariation of human mandibular and maxillary first molars

Dental casts of 160 Greek subjects (80 males, 80 females) were scanned by a structured‐light scanner. The upper and lower right first molar occlusal surface 3D meshes were processed using geometric morphometric methods. A total of 265 and 274 curve and surface sliding semilandmarks were placed on the upper and lower molar surfaces, respectively. Principal component analysis and partial least square analysis were performed to assess shape parameters. Molars tended to vary between an elongated and a more square form. The first two principal components (PCs), comprising almost 1/3 of molar shape variation, were related to mesiodistal–buccolingual ratios and relative cusp position. Distal cusps displayed the greatest shape variability. Molars of males were larger than those of females (2.8 and 3.2% for upper and lower molars respectively), but no shape dimorphism was observed. Upper and lower molar sizes were significantly correlated (r² = 0.689). Allometry was observed for both teeth. Larger lower molars were associated with shorter cusps, expansion of the distal cusp, and constriction of the mesial cusps (predicted variance 3.25%). Upper molars displayed weaker allometry (predicted variance 1.59%). Upper and lower molar shape covariation proved significant (RV = 17.26%, P < 0.0001). The main parameter of molar covariation in partial least square axis 1, contributing to 30% of total covariation, was cusp height, in contrast to the primary variability traits exhibited by PC1 and PC2. The aim of this study was to evaluate shape variation and covariation, including allometry and sexual dimorphism, of maxillary and mandibular first permanent molar occlusal surfaces. Am J Phys Anthropol 152:186–196, 2013. © 2013 Wiley Periodicals, Inc.     

Mapping molar shapes on signaling pathways

A major challenge in evolutionary developmental biology is to understand how genetic mutations underlie phenotypic changes. In principle, selective pressures on the phenotype screen the gene pool of the population. Teeth are an excellent model for understanding evolutionary changes in the genotype-phenotype relationship since they exist throughout vertebrates. Genetically modified mice (mutants) with abnormalities in teeth have been used to explore tooth development. The relationship between signaling pathways and molar shape, however, remains elusive due to the high intrinsic complexity of tooth crowns. This hampers our understanding of the extent to which developmental factors explored in mutants explain developmental and phenotypic variation in natural species that represent the consequence of natural selection. Here we combine a novel morphometric method with two kinds of data mining techniques to extract data sets from the three-dimensional surface models of lower first molars: i) machine learning to maximize classification accuracy of 22 mutants, and ii) phylogenetic signal for 31 Murinae species. Major shape variation among mutants is explained by the number of cusps and cusp distribution on a tooth crown. The distribution of mutant mice in morphospace suggests a nonlinear relationship between the signaling pathways and molar shape variation. Comparative analysis of mutants and wild murines reveals that mutant variation overlaps naturally occurring diversity, including more ancestral and derived morphologies. However, taxa with transverse lophs are not fully covered by mutant variation, suggesting experimentally unexplored developmental factors in the evolutionary radiation of Murines.