The Inhibitory Cascade Model is Not a Good Predictor of Molar Size Covariation

Sex-biased sampling may influence Homo naledi tooth size variation

A frequent source of debate in paleoanthropology concerns the taxonomic unity of fossil assemblages, with many hominin samples exhibiting elevated levels of variation that can be interpreted as indicating the presence of multiple species. By contrast, the large assemblage of hominin fossils from the Rising Star cave system, assigned to Homo naledi, exhibits a remarkably low degree of variation for most skeletal elements. Many factors can contribute to low sample variation, including genetic drift, strong natural selection, biased sex ratios, and sampling of closely related individuals. In this study, we tested for potential sex-biased sampling in the Rising Star dental sample. We compared coefficients of variation for the H. naledi teeth to those for eight extant hominoid samples. We used a resampling procedure that generated samples from the extant taxa that matched the sample size of the fossil sample for each possible Rising Star dental sex ratio. We found that variation at four H. naledi tooth positions-I2, M1, P4, M1-is so low that the possibility that one sex is represented by few or no individuals in the sample cannot be excluded. Additional evidence is needed to corroborate this inference, such as ancient DNA or enamel proteome data, and our study design does not address other potential factors that would account for low sample variation. Nevertheless, our results highlight the importance of considering the taphonomic history of a hominin assemblage and suggest that sex-biased sampling is a plausible explanation for the low level of phenotypic variation found in some aspects of the current H. naledi assemblage.

Hypoconulid loss in cercopithecins: Functional and developmental considerations

Cercopithecins differ from papionins in lacking a M3 hypoconulid. Although this loss may be related to dietary differences, the functional and developmental ramifications of hypoconulid loss are currently unclear. The following makes use of dental topographic analysis to quantify shape variation in a sample of cercopithecin M3s, as well as in a sample of Macaca, which has a hypoconulid. To help understand the consequences of hypoconulid loss, Macaca M3s were virtually cropped to remove the hypoconulid and were also subjected to dental topographic analysis. The patterning cascade model and the inhibitory cascade model attempt to explain variation in cusp pattern and molar proportions, respectively. These models have both previously been used to explain patterns of variation in cercopithecines, but have not been examined in the context of hypoconulid loss. For example, previous work suggests that earlier developing cusps impact the development of later developing cusps (i.e., the hypoconulid) and that cercopithecines do not conform to the predictions of the inhibitory cascade model in that the size of the molars is not linear moving distally. Results of the current study suggest that the loss of the hypoconulid is associated with a reduction in dental topography among cercopithecins, which is potentially related to diet, although the connection to diet is not necessarily clear. Results also suggest that the loss of the hypoconulid can be explained by the patterning cascade model, and that hypoconulid loss explains the apparent lack of support for the inhibitory cascade model among cercopithecines. These findings highlight the importance of a holistic approach to studying variation in molar proportions and developmental models.

A tooth crown morphology framework for interpreting the diversity of primate dentitions

Variation in tooth crown morphology plays a crucial role in species diagnoses, phylogenetic inference, and the reconstruction of the evolutionary history of the primate clade. While a growing number of studies have identified developmental mechanisms linked to tooth size and cusp patterning in mammalian crown morphology, it is unclear (1) to what degree these are applicable across primates and (2) which additional developmental mechanisms should be recognized as playing important roles in odontogenesis. From detailed observations of lower molar enamel-dentine junction morphology from taxa representing the major primate clades, we outline multiple phylogenetic and developmental components responsible for crown patterning, and formulate a tooth crown morphology framework for the holistic interpretation of primate crown morphology. We suggest that adopting this framework is crucial for the characterization of tooth morphology in studies of dental development, discrete trait analysis, and systematics.

Bat teeth illuminate the diversification of mammalian tooth classes

Tooth classes are an innovation that has contributed to the evolutionary success of mammals. However, our understanding of the mechanisms by which tooth classes diversified remain limited. We use the evolutionary radiation of noctilionoid bats to show how the tooth developmental program evolved during the adaptation to new diet types. Combining morphological, developmental and mathematical modeling approaches, we demonstrate that tooth classes develop through independent developmental cascades that deviate from classical models. We show that the diversification of tooth number and size is driven by jaw growth rate modulation, explaining the rapid gain/loss of teeth in this clade. Finally, we mathematically model the successive appearance of tooth buds, supporting the hypothesis that growth acts as a key driver of the evolution of tooth number and size. Our work reveal how growth, by tinkering with reaction/diffusion processes, drives the diversification of tooth classes and other repeated structure during adaptive radiations.

Keeping 21st Century Paleontology Grounded: Quantitative Genetic Analyses and Ancestral State Reconstruction Re-Emphasize the Essentiality of Fossils

Advances in genetics and developmental biology are revealing the relationship between genotype and dental phenotype (G:P), providing new approaches for how paleontologists assess dental variation in the fossil record. Our aim was to understand how the method of trait definition influences the ability to reconstruct phylogenetic relationships and evolutionary history in the Cercopithecidae, the Linnaean Family of monkeys currently living in Africa and Asia. We compared the two-dimensional assessment of molar size (calculated as the mesiodistal length of the crown multiplied by the buccolingual breadth) to a trait that reflects developmental influences on molar development (the inhibitory cascade, IC) and two traits that reflect the genetic architecture of postcanine tooth size variation (defined through quantitative genetic analyses: MMC and PMM). All traits were significantly influenced by the additive effects of genes and had similarly high heritability estimates. The proportion of covariate effects was greater for two-dimensional size compared to the G:P-defined traits. IC and MMC both showed evidence of selection, suggesting that they result from the same genetic architecture. When compared to the fossil record, Ancestral State Reconstruction using extant taxa consistently underestimated MMC and PMM values, highlighting the necessity of fossil data for understanding evolutionary patterns in these traits. Given that G:P-defined dental traits may provide insight to biological mechanisms that reach far beyond the dentition, this new approach to fossil morphology has the potential to open an entirely new window onto extinct paleobiologies. Without the fossil record, we would not be able to grasp the full range of variation in those biological mechanisms that have existed throughout evolution.

Testing the inhibitory cascade model in a recent human sample

The Inhibitory Cascade Model was proposed by Kavanagh and colleagues (Nature, 449, 427-433 [2007]) after their experimental studies on the dental development of murine rodent species. These authors described an activator-inhibitor mechanism that has been employed to predict evolutionary size patterns of mammalian teeth, including hominins. In the present study, we measured the crown area of the three lower permanent molars (M1, M2, and M3) of a large recent modern human sample of male and female individuals from a collection preserved at the Institute of Anthropology of the University of Coimbra (Portugal). The main aim of the present study is to test if the size molar patterns observed in this human sample fits the Inhibitory Cascade Model. For this purpose, we first measured the crown area in those individuals preserving the complete molar series. Measurements were taken in photographs, using a planimeter and following well-tested techniques used in previous works. We then plot the M3 /M1 and M2 /M1  size ratios. Our results show that the premise of the Inhibitory Cascade Model, according to which the average of the crown area of M2 is approximately one-third of the sum of the crown area of the three molars, is fulfilled. However, our results also show that the individual values of a significant number of males and females are out of the 95% confidence interval predicted by the Inhibitory Cascade Model in rodents. As a result, the present analyses suggest that neither the sample of males, nor that of females, nor the pooled sample fits the Inhibitory Cascade Model. It is important to notice that, although this model has been successfully tested in a large number of current human populations, to the best of our knowledge this is the first study in which individual data have been obtained in a recent human population rather than using the average of the sample. Our results evince that, at the individual level, some factors not yet known could interfere with this model masking the modulation of the size on the molar series in modern humans. We suggest that the considerable delay in the onset of M3 formation in modern humans could be related to a weakening of the possible activation/inhibition process for this tooth. Finally, and in support of our conclusions, we have checked that the absolute and relative size of M1 and M2 is not related to the M3 agenesis in our sample. In line with other studies in primates, our results do not support the Inhibitory Cascade Model in a recent human sample. Further research is needed to better understand the genetic basis of this mechanism and its relationship to the phenotype. In this way, we may be able to find out which evolutionary changes may be responsible for the deviations observed in many species, including Homo sapiens.

Unexpected variation of human molar size patterns

A tenet of mammalian, including primate dental evolution, is the Inhibitory Cascade Model, where first molar (M1) size predicts in a linear cline the size and onset time of the second (M2) and third (M3) molars: a larger M1 portends a progressively smaller and later-developing M2 and M3. In contemporary modern Homo sapiens, later-developing M3s are less likely to erupt properly. The Inhibitory Cascade Model is also used to predict molar sizes of extinct taxa, including fossil Homo. The extent to which Inhibitory Cascade Model predictions hold in contemporary H. sapiens molars is unclear, including whether this tenet informs about molar initiation, development, and eruption. We tested these questions here. In our radiographic sample of 323 oral quadrants and molar rows from contemporary humans based on mesiodistal crown lengths, we observed the distribution of molar proportions with a central tendency around parity (M1 = M2 = M3) that parsed into 13 distinct molar size ratio patterns. These patterns presented at different frequencies (e.g., M1 > M2 > M3 in about one-third of cases) that reflected whether the molar row was located in the maxilla or mandible and included both linear (e.g., M1 < M2 < M3) and nonlinear molar size ratio progressions (e.g., M1 > M2 < M3). Up to four patterns were found in the same subject's mouth. Lastly, M1 size alone does not predict M3 size, developmental timing, or eruption; rather, M2 size is integral to predicting M3 size. Our study indicates that human molar size is genetically 'softwired' and sensitive to factors local to the human upper jaw vs. lower jaw. The lack of a single stereotypical molar size ratio for contemporary H. sapiens suggests that predictions of fossil H. sapiens molar sizes using the Inhibitory Cascade Model must be made with caution.

Segmental series and size: clade-wide investigation of molar proportions reveals a major evolutionary allometry in the dentition of placental mammals

Iterative segments such as teeth or limbs are a widespread characteristic of living organisms. While their proportions may be governed by similar developmental rules in vertebrates, there is no emerging pattern as regards their relation to size. Placental mammals span eight orders of magnitude in body size and show a wide spectrum of dietary habits associated with size and reflected in their dentitions, especially molars. Although variation in size constitutes an important determinant for variation in biological traits, few major allometric trends have been documented on placental molars so far. Molar proportions have been intensively explored in placentals in relation to developmental models, but often at a small phylogenetic scale. Here, we analyzed the diversity of upper molar proportions in relation to absolute size in a large sample of placental species (n = 286) encompassing most of the group's dental diversity. Our phylogenetically informed analyses revealed a twofold pattern of evolutionary integration among upper molars: while molars covary in size with each other, their proportions covary with the absolute size of the entire molar field. With increasing absolute size, posterior molars increase in size relative to anterior ones, meaning that large-sized species have relatively large rear molars while the opposite is true for small-sized species. The directionality of proportional increase in the molar row exhibits a previously unsuspected allometric patterning among placentals, showing how large-scale variations in size may have influenced variation in dental morphology. This finding provides new evidence that processes regulating the size of individual molars are integrated with overall patterns of growth and calls for further testing of allometric variation in the dentition and in other segmental series of the vertebrate body.

Mammalian molar complexity follows simple, predictable patterns

Identifying developmental explanations for the evolution of complex structures like mammalian molars is fundamental to studying phenotypic variation. Previous study showed that a "morphogenetic gradient" of molar proportions was explained by a balance between inhibiting/activating activity from earlier developing molars, termed the inhibitory cascade model (ICM). Although this model provides an explanation for variation in molar proportions, what remains poorly understood is if molar shape, or specifically complexity (i.e., the number of cusps, crests), can be explained by the same developmental model. Here, we show that molar complexity conforms to the ICM, following a linear, morphogenetic gradient along the molar row. Moreover, differing levels of inhibiting/activating activity produce contrasting patterns of molar complexity depending on diet. This study corroborates a model for the evolution of molar complexity that is developmentally simple, where only small-scale developmental changes need to occur to produce change across the entire molar row, with this process being mediated by an animal's ecology. The ICM therefore provides a developmental framework for explaining variation in molar complexity and a means for testing developmental hypotheses in the broader context of mammalian evolution.

Genetic correlations in the rhesus macaque dentition

Quantitative genetic analyses can indicate how complex traits respond to natural selection by demonstrating the genetic relationships between features that constrain their evolution. Genetic correlations between dental measurements have been estimated previously in baboons, humans, and tamarins and indicate variable patterns of modularity by tooth type across these taxa. Here, heritabilities of, and genetic correlations between, linear dental measurements were estimated from the Cayo Santiago rhesus macaques (Macaca mulatta). Relationships between the genetic correlation matrix and matrices designed to test hypotheses of modularity by tooth type, region, function, and development were assessed using a random skewers approach. Dental measurements were found to be moderately to highly heritable, with 24 of 28 heritability estimates differing significantly (p < 0.05) from zero. Almost all genetic correlations between dental dimensions were positive. The genetic correlation matrix was most similar to a regionally modular matrix, with distinct anterior and postcanine tooth modules. This pattern is consistent with previous quantitative genetic analyses of baboons and previous phenotypic analyses of cercopithecoid primates. The existence of a genetic module for the canines and honing premolar was not supported. Ongoing selection pressures, rather than strong genetic constraints, are likely necessary to preserve functional relationships between the canines and honing premolar based on these findings. The genetic correlation matrix of the Cayo Santiago rhesus macaques mirrors patterns of phenotypic correlations observed for cercopithecoid primates broadly and demonstrates that genetic contributions to these patterns may be fairly stable over the course of cercopithecoid evolution. The quantitative genetic study of additional taxa will be necessary to determine whether the regional modularity of baboons and macaques, or the integrated pattern of humans and tamarins, is shared more broadly across primates.

A genotype:phenotype approach to testing taxonomic hypotheses in hominids

Paleontology has long relied on assumptions about the genetic and developmental influences on skeletal variation. The last few decades of developmental genetics have elucidated the genetic pathways involved in making teeth and patterning the dentition. Quantitative genetic analyses have refined this genotype:phenotype map even more, especially for primates. We now have the ability to define dental traits with a fair degree of fidelity to the underlying genetic architecture; for example, the molar module component (MMC) and the premolar-molar module (PMM) that have been defined through quantitative genetic analyses. We leverage an extensive dataset of extant and extinct hominoid dental variation to explore how these two genetically patterned phenotypes have evolved through time. We assess MMC and PMM to test the hypothesis that these two traits reveal a more biologically informed taxonomy at the genus and species levels than do more traditional measurements. Our results indicate that MMC values for hominids fall into two categories and that Homo is derived compared with earlier taxa. We find a more variable, species-level pattern for PMM. These results, in combination with previous research, demonstrate that MMC reflects the phenotypic output of a more evolutionarily stable, or phylogenetically congruent, genetic mechanism, and PMM is a reflection of a more evolutionarily labile mechanism. These results suggest that the human lineage since the split with chimpanzees may not represent as much genus-level variation as has been inferred from traits whose etiologies are not understood.

Mammal Molar Size Ratios and the Inhibitory Cascade at the Intraspecific Scale

Mammalian molar crowns form a module in which measurements of size for individual teeth within a tooth row covary with one another. Molar crown size covariation is proposed to fit the inhibitory cascade model (ICM) or its variant the molar module component (MMC) model, but the inability of the former model to fit across biological scales is a concern in the few cases where it has been tested in Primates. The ICM has thus far failed to explain patterns of intraspecific variation, an intermediate biological scale, even though it explains patterns at both smaller organ-level and larger between-species biological scales. Studies of this topic in a much broader range of taxa are needed, but the properties of a sample appropriate for testing the ICM at the intraspecific level are unclear. Here, we assess intraspecific variation in relative molar sizes of the cotton mouse, Peromyscus gossypinus, to further test the ICM and to develop recommendations for appropriate sampling protocols in future intraspecific studies of molar size variation across Mammalia. To develop these recommendations, we model the sensitivity of estimates of molar ratios to sample size and simulate the use of composite molar rows when complete ones are unavailable. Similar to past studies on primates, our results show that intraspecific variance structure of molar ratios within the rodent P. gossypinus does not meet predictions of the ICM or MMC. When we extend these analyses to include the MMC, one model does not fit observed patterns of variation better than the other. Standing variation in molar size ratios is relatively constant across mammalian samples containing all three molars. In future studies, analyzing average ratio values will require relatively small minimum sample sizes of two or more complete molar rows. Even composite-based estimates from four or more specimens per tooth position can accurately estimate mean molar ratios. Analyzing variance structure will require relatively large sample sizes of at least 40-50 complete specimens, and composite molar rows cannot accurately reconstruct variance structure of ratios in a sample. Based on these results, we propose guidelines for intraspecific studies of molar size covariation. In particular, we note that the suitability of composite specimens for averaging mean molar ratios is promising for the inclusion of isolated molars and incomplete molar rows from the fossil record in future studies of the evolution of molar modules, as long as variance structure is not a key component of such studies.