The developmental basis for scaling of mammalian tooth size.

When evolution leads to differences in body size, organs generally scale along. A well-known example of the tight relationship between organ and body size is the scaling of mammalian molar teeth. To investigate how teeth scale during development and evolution, we compared molar development from initiation through final size in the mouse and the rat. Whereas the linear dimensions of the rat molars are twice that of the mouse molars, their shapes are largely the same. Here, we focus on the first lower molars that are considered the most reliable dental proxy for size-related patterns due to their low within-species variability. We found that scaling of the molars starts early, and that the rat molar is patterned equally as fast but in a larger size than the mouse molar. Using transcriptomics, we discovered that a known regulator of body size, insulin-like growth factor 1 (Igf1), is more highly expressed in the rat molars compared to the mouse molars. Ex vivo and in vivo mouse models demonstrated that modulation of the IGF pathway reproduces several aspects of the observed scaling process. Furthermore, analysis of IGF1-treated mouse molars and computational modeling indicate that IGF signaling scales teeth by simultaneously enhancing growth and by inhibiting the cusp-patterning program, thereby providing a relatively simple mechanism for scaling teeth during development and evolution. Finally, comparative data from shrews to elephants suggest that this scaling mechanism regulates the minimum tooth size possible, as well as the patterning potential of large teeth.

Chromosomal neighbourhoods allow identification of organ specific changes in gene expression.

Although most genes share their chromosomal neighbourhood with other genes, distribution of genes has not been explored in the context of individual organ development; the common focus of developmental biology studies. Because developmental processes are often associated with initially subtle changes in gene expression, here we explored whether neighbouring genes are informative in the identification of differentially expressed genes. First, we quantified the chromosomal neighbourhood patterns of genes having related functional roles in the mammalian genome. Although the majority of protein coding genes have at least five neighbours within 1 Mb window around each gene, very few of these neighbours regulate development of the same organ. Analyses of transcriptomes of developing mouse molar teeth revealed that whereas expression of genes regulating tooth development changes, their neighbouring genes show no marked changes, irrespective of their level of expression. Finally, we test whether inclusion of gene neighbourhood in the analyses of differential expression could provide additional benefits. For the analyses, we developed an algorithm, called DELocal that identifies differentially expressed genes by comparing their expression changes to changes in adjacent genes in their chromosomal regions. Our results show that DELocal removes detection bias towards large changes in expression, thereby allowing identification of even subtle changes in development. Future studies, including the detection of differential expression, may benefit from, and further characterize the significance of gene-gene neighbour relationships.

A simple rule governs the evolution and development of hominin tooth size.

The variation in molar tooth size in humans and our closest relatives (hominins) has strongly influenced our view of human evolution. The reduction in overall size and disproportionate decrease in third molar size have been noted for over a century, and have been attributed to reduced selection for large dentitions owing to changes in diet or the acquisition of cooking. The systematic pattern of size variation along the tooth row has been described as a 'morphogenetic gradient' in mammal, and more specifically hominin, teeth since Butler and Dahlberg. However, the underlying controls of tooth size have not been well understood, with hypotheses ranging from morphogenetic fields to the clone theory. In this study we address the following question: are there rules that govern how hominin tooth size evolves? Here we propose that the inhibitory cascade, an activator-inhibitor mechanism that affects relative tooth size in mammals, produces the default pattern of tooth sizes for all lower primary postcanine teeth (deciduous premolars and permanent molars) in hominins. This configuration is also equivalent to a morphogenetic gradient, finally pointing to a mechanism that can generate this gradient. The pattern of tooth size remains constant with absolute size in australopiths (including Ardipithecus, Australopithecus and Paranthropus). However, in species of Homo, including modern humans, there is a tight link between tooth proportions and absolute size such that a single developmental parameter can explain both the relative and absolute sizes of primary postcanine teeth. On the basis of the relationship of inhibitory cascade patterning with size, we can use the size at one tooth position to predict the sizes of the remaining four primary postcanine teeth in the row for hominins. Our study provides a development-based expectation to examine the evolution of the unique proportions of human teeth.

System-level analyses of keystone genes required for mammalian tooth development.

When a null mutation of a gene causes a complete developmental arrest, the gene is typically considered essential for life. Yet, in most cases, null mutations have more subtle effects on the phenotype. Here we used the phenotypic severity of mutations as a tool to examine system-level dynamics of gene expression. We classify genes required for the normal development of the mouse molar into different categories that range from essential to subtle modification of the phenotype. Collectively, we call these the developmental keystone genes. Transcriptome profiling using microarray and RNAseq analyses of patterning stage mouse molars show highly elevated expression levels for genes essential for the progression of tooth development, a result reminiscent of essential genes in single-cell organisms. Elevated expression levels of progression genes were also detected in developing rat molars, suggesting evolutionary conservation of this system-level dynamics. Single-cell RNAseq analyses of developing mouse molars reveal that even though the size of the expression domain, measured in the number of cells, is the main driver of organ-level expression, progression genes show high cell-level transcript abundances. Progression genes are also upregulated within their pathways, which themselves are highly expressed. In contrast, a high proportion of the genes required for normal tooth patterning are secreted ligands that are expressed in fewer cells than their receptors and intracellular components. Overall, even though expression patterns of individual genes can be highly different, conserved system-level principles of gene expression can be detected using phenotypically defined gene categories.

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.

Modeling enamel matrix secretion in mammalian teeth.

The most mineralized tissue of the mammalian body is tooth enamel. Especially in species with thick enamel, three-dimensional (3D) tomography data has shown that the distribution of enamel varies across the occlusal surface of the tooth crown. Differences in enamel thickness among species and within the tooth crown have been used to examine taxonomic affiliations, life history, and functional properties of teeth. Before becoming fully mineralized, enamel matrix is secreted on the top of a dentine template, and it remains to be explored how matrix thickness is spatially regulated. To provide a predictive framework to examine enamel distribution, we introduce a computational model of enamel matrix secretion that maps the dentine topography to the enamel surface topography. Starting from empirical enamel-dentine junctions, enamel matrix deposition is modeled as a diffusion-limited free boundary problem. Using laboratory microCT and synchrotron tomographic data of pig molars that have markedly different dentine and enamel surface topographies, we show how diffusion-limited matrix deposition accounts for both the process of matrix secretion and the final enamel distribution. Simulations reveal how concave and convex dentine features have distinct effects on enamel surface, thereby explaining why the enamel surface is not a straightforward extrapolation of the dentine template. Human and orangutan molar simulations show that even subtle variation in dentine topography can be mapped to the enamel surface features. Mechanistic models of extracellular matrix deposition can be used to predict occlusal morphologies of teeth.

Replaying evolutionary transitions from the dental fossil record.

The evolutionary relationships of extinct species are ascertained primarily through the analysis of morphological characters. Character inter-dependencies can have a substantial effect on evolutionary interpretations, but the developmental underpinnings of character inter-dependence remain obscure because experiments frequently do not provide detailed resolution of morphological characters. Here we show experimentally and computationally how gradual modification of development differentially affects characters in the mouse dentition. We found that intermediate phenotypes could be produced by gradually adding ectodysplasin A (EDA) protein in culture to tooth explants carrying a null mutation in the tooth-patterning gene Eda. By identifying development-based character inter-dependencies, we show how to predict morphological patterns of teeth among mammalian species. Finally, in vivo inhibition of sonic hedgehog signalling in Eda null teeth enabled us to reproduce characters deep in the rodent ancestry. Taken together, evolutionarily informative transitions can be experimentally reproduced, thereby providing development-based expectations for character-state transitions used in evolutionary studies.

Mechanical constraint from growing jaw facilitates mammalian dental diversity.

Much of the basic information about individual organ development comes from studies using model species. Whereas conservation of gene regulatory networks across higher taxa supports generalizations made from a limited number of species, generality of mechanistic inferences remains to be tested in tissue culture systems. Here, using mammalian tooth explants cultured in isolation, we investigate self-regulation of patterning by comparing developing molars of the mouse, the model species of mammalian research, and the bank vole. A distinct patterning difference between the vole and the mouse molars is the alternate cusp offset present in the vole. Analyses of both species using 3D reconstructions of developing molars and jaws, computational modeling of cusp patterning, and tooth explants cultured with small braces show that correct cusp offset requires constraints on the lateral expansion of the developing tooth. Vole molars cultured without the braces lose their cusp offset, and mouse molars cultured with the braces develop a cusp offset. Our results suggest that cusp offset, which changes frequently in mammalian evolution, is more dependent on the 3D support of the developing jaw than other aspects of tooth shape. This jaw-tooth integration of a specific aspect of the tooth phenotype indicates that organs may outsource specific aspects of their morphology to be regulated by adjacent body parts or organs. Comparative studies of morphologically different species are needed to infer the principles of organogenesis.

Differential tissue growth and cell adhesion alone drive early tooth morphogenesis: An ex vivo and in silico study.

From gastrulation to late organogenesis animal development involves many genetic and bio-mechanical interactions between epithelial and mesenchymal tissues. Ectodermal organs, such as hairs, feathers and teeth are well studied examples of organs whose development is based on epithelial-mesenchymal interactions. These develop from a similar primordium through an epithelial folding and its interaction with the mesenchyme. Despite extensive knowledge on the molecular pathways involved, little is known about the role of bio-mechanical processes in the morphogenesis of these organs. We propose a simple computational model for the biomechanics of one such organ, the tooth, and contrast its predictions against cell-tracking experiments, mechanical relaxation experiments and the observed tooth shape changes over developmental time. We found that two biomechanical processes, differential tissue growth and differential cell adhesion, were enough, in the model, for the development of the 3D morphology of the early tooth germ. This was largely determined by the length and direction of growth of the cervical loops, lateral folds of the enamel epithelium. The formation of these cervical loops was found to require accelerated epithelial growth relative to other tissues and their direction of growth depended on specific differential adhesion between the three tooth tissues. These two processes and geometrical constraints in early tooth bud also explained the shape asymmetry between the lateral cervical loops and those forming in the anterior and posterior of the tooth. By performing mechanical perturbations ex vivo and in silico we inferred the distribution and direction of tensile stresses in the mesenchyme that restricted cervical loop lateral growth and forced them to grow downwards. Overall our study suggests detailed quantitative explanations for how bio-mechanical processes lead to specific morphological 3D changes over developmental time.

On the difficulty of increasing dental complexity.

One of the fascinating aspects of the history of life is the apparent increase in morphological complexity through time, a well known example being mammalian cheek tooth evolution. In contrast, experimental studies of development more readily show a decrease in complexity, again well exemplified by mammalian teeth, in which tooth crown features called cusps are frequently lost in mutant and transgenic mice. Here we report that mouse tooth complexity can be increased substantially by adjusting multiple signalling pathways simultaneously. We cultured teeth in vitro and adjusted ectodysplasin (EDA), activin A and sonic hedgehog (SHH) pathways, all of which are individually required for normal tooth development. We quantified tooth complexity using the number of cusps and a topographic measure of surface complexity. The results show that whereas activation of EDA and activin A signalling, and inhibition of SHH signalling, individually cause subtle to moderate increases in complexity, cusp number is doubled when all three pathways are adjusted in unison. Furthermore, the increase in cusp number does not result from an increase in tooth size, but from an altered primary patterning phase of development. The combination of a lack of complex mutants, the paucity of natural variants with complex phenotypes, and our results of greatly increased dental complexity using multiple pathways, suggests that an increase may be inherently different from a decrease in phenotypic complexity.

Adaptive radiation of multituberculate mammals before the extinction of dinosaurs.

The Cretaceous-Paleogene mass extinction approximately 66 million years ago is conventionally thought to have been a turning point in mammalian evolution. Prior to that event and for the first two-thirds of their evolutionary history, mammals were mostly confined to roles as generalized, small-bodied, nocturnal insectivores, presumably under selection pressures from dinosaurs. Release from these pressures, by extinction of non-avian dinosaurs at the Cretaceous-Paleogene boundary, triggered ecological diversification of mammals. Although recent individual fossil discoveries have shown that some mammalian lineages diversified ecologically during the Mesozoic era, comprehensive ecological analyses of mammalian groups crossing the Cretaceous-Paleogene boundary are lacking. Such analyses are needed because diversification analyses of living taxa allow only indirect inferences of past ecosystems. Here we show that in arguably the most evolutionarily successful clade of Mesozoic mammals, the Multituberculata, an adaptive radiation began at least 20 million years before the extinction of non-avian dinosaurs and continued across the Cretaceous-Paleogene boundary. Disparity in dental complexity, which relates to the range of diets, rose sharply in step with generic richness and disparity in body size. Moreover, maximum dental complexity and body size demonstrate an adaptive shift towards increased herbivory. This dietary expansion tracked the ecological rise of angiosperms and suggests that the resources that were available to multituberculates were relatively unaffected by the Cretaceous-Paleogene mass extinction. Taken together, our results indicate that mammals were able to take advantage of new ecological opportunities in the Mesozoic and that at least some of these opportunities persisted through the Cretaceous-Paleogene mass extinction. Similar broad-scale ecomorphological inventories of other radiations may help to constrain the possible causes of mass extinctions.

A computational model of teeth and the developmental origins of morphological variation.

The relationship between the genotype and the phenotype, or the genotype-phenotype map, is generally approached with the tools of multivariate quantitative genetics and morphometrics. Whereas studies of development and mathematical models of development may offer new insights into the genotype-phenotype map, the challenge is to make them useful at the level of microevolution. Here we report a computational model of mammalian tooth development that combines parameters of genetic and cellular interactions to produce a three-dimensional tooth from a simple tooth primordia. We systematically tinkered with each of the model parameters to generate phenotypic variation and used geometric morphometric analyses to identify, or developmentally ordinate, parameters best explaining population-level variation of real teeth. To model the full range of developmentally possible morphologies, we used a population sample of ringed seals (Phoca hispida ladogensis). Seal dentitions show a high degree of variation, typically linked to the lack of exact occlusion. Our model suggests that despite the complexity of development and teeth, there may be a simple basis for dental variation. Changes in single parameters regulating signalling during cusp development may explain shape variation among individuals, whereas a parameter regulating epithelial growth may explain serial, tooth-to-tooth variation along the jaw. Our study provides a step towards integrating the genotype, development and the phenotype.

Tooth shape formation and tooth renewal: evolving with the same signals.

Teeth are found in almost all vertebrates, and they therefore provide a general paradigm for the study of epithelial organ development and evolution. Here, we review the developmental mechanisms underlying changes in tooth complexity and tooth renewal during evolution, focusing on recent studies of fish, reptiles and mammals. Mammals differ from other living vertebrates in that they have the most complex teeth with restricted capacity for tooth renewal. As we discuss, however, limited tooth replacement in mammals has been compensated for in some taxa by the evolution of continuously growing teeth, the development of which appears to reuse the regulatory pathways of tooth replacement.

Ectodysplasin regulates activator-inhibitor balance in murine tooth development through Fgf20 signaling.

Uncovering the origin and nature of phenotypic variation within species is the first step in understanding variation between species. Mouse models with altered activities of crucial signal pathways have highlighted many important genes and signal networks regulating the morphogenesis of complex structures, such as teeth. The detailed analyses of these models have indicated that the balanced actions of a few pathways regulating cell behavior modulate the shape and number of teeth. Currently, however, most mouse models studied have had gross alteration of morphology, whereas analyses of more subtle modification of morphology are required to link developmental studies to evolutionary change. Here, we have analyzed a signaling network involving ectodysplasin (Eda) and fibroblast growth factor 20 (Fgf20) that subtly affects tooth morphogenesis. We found that Fgf20 is a major downstream effector of Eda and affects Eda-regulated characteristics of tooth morphogenesis, including the number, size and shape of teeth. Fgf20 function is compensated for by other Fgfs, in particular Fgf9 and Fgf4, and is part of an Fgf signaling loop between epithelium and mesenchyme. We showed that removal of Fgf20 in an Eda gain-of-function mouse model results in an Eda loss-of-function phenotype in terms of reduced tooth complexity and third molar appearance. However, the extra anterior molar, a structure lost during rodent evolution 50 million years ago, was stabilized in these mice.

Genetic comparisons yield insight into the evolution of enamel thickness during human evolution.

Enamel thickness varies substantially among extant hominoids and is a key trait with significance for interpreting dietary adaptation, life history trajectory, and phylogenetic relationships. There is a strong link in humans between enamel formation and mutations in the exons of the four genes that code for the enamel matrix proteins and the associated protease. The evolution of thick enamel in humans may have included changes in the regulation of these genes during tooth development. The cis-regulatory region in the 5' flank (upstream non-coding region) of MMP20, which codes for enamelysin, the predominant protease active during enamel secretion, has previously been shown to be under strong positive selection in the lineages leading to both humans and chimpanzees. Here we examine evidence for positive selection in the 5' flank and 3' flank of AMELX, AMBN, ENAM, and MMP20. We contrast the human sequence changes with other hominoids (chimpanzees, gorillas, orangutans, gibbons) and rhesus macaques (outgroup), a sample comprising a range of enamel thickness. We find no evidence for positive selection in the protein-coding regions of any of these genes. In contrast, we find strong evidence for positive selection in the 5' flank region of MMP20 and ENAM along the lineage leading to humans, and in both the 5' flank and 3' flank regions of MMP20 along the lineage leading to chimpanzees. We also identify putative transcription factor binding sites overlapping some of the species-specific nucleotide sites and we refine which sections of the up- and downstream putative regulatory regions are most likely to harbor important changes. These non-coding changes and their potential for differential regulation by transcription factors known to regulate tooth development may offer insight into the mechanisms that allow for rapid evolutionary changes in enamel thickness across closely-related species, and contribute to our understanding of the enamel phenotype in hominoids.

Algorithms to automatically quantify the geometric similarity of anatomical surfaces.

We describe approaches for distances between pairs of two-dimensional surfaces (embedded in three-dimensional space) that use local structures and global information contained in interstructure geometric relationships. We present algorithms to automatically determine these distances as well as geometric correspondences. This approach is motivated by the aspiration of students of natural science to understand the continuity of form that unites the diversity of life. At present, scientists using physical traits to study evolutionary relationships among living and extinct animals analyze data extracted from carefully defined anatomical correspondence points (landmarks). Identifying and recording these landmarks is time consuming and can be done accurately only by trained morphologists. This necessity renders these studies inaccessible to nonmorphologists and causes phenomics to lag behind genomics in elucidating evolutionary patterns. Unlike other algorithms presented for morphological correspondences, our approach does not require any preliminary marking of special features or landmarks by the user. It also differs from other seminal work in computational geometry in that our algorithms are polynomial in nature and thus faster, making pairwise comparisons feasible for significantly larger numbers of digitized surfaces. We illustrate our approach using three datasets representing teeth and different bones of primates and humans, and show that it leads to highly accurate results.

Wnt target gene Ascl4 is dispensable for skin appendage development.

The development of skin appendages, including hair follicles, teeth and mammary glands is initiated through the formation of the placode, a local thickening of the epithelium. The Wnt/ß-catenin signaling cascade is an evolutionary conserved pathway with an essential role in placode morphogenesis, but its downstream targets and their exact functions remain ill defined. In this study, we identify Achaete-scute complex-like 4 (Ascl4) as a novel target of the Wnt/ß-catenin pathway and demonstrate its expression pattern in the signaling centers of developing hair follicles and teeth. Ascl transcription factors belong to the superfamily of basic helix-loop-helix transcriptional regulators involved in cell fate determination in many tissues. However, their specific role in the developing skin remains largely unknown. We report that Ascl4 null mice have no overt phenotype. Absence of Ascl4 did not impair hair follicle morphogenesis or hair shaft formation suggesting that it is non-essential for hair follicle development. No tooth or mammary gland abnormalities were detected either. We suggest that other transcription factors may functionally compensate for the absence of Ascl4, but further research is warranted to assess this possibility.

Vole genomics links determinate and indeterminate growth of teeth.

Continuously growing teeth are an important innovation in mammalian evolution, yet genetic regulation of continuous growth by stem cells remains incompletely understood. Dental stem cells responsible for tooth crown growth are lost at the onset of tooth root formation. Genetic signaling that initiates this loss is difficult to study with the ever-growing incisor and rooted molars of mice, the most common mammalian dental model species, because signals for root formation overlap with signals that pattern tooth size and shape (i.e., cusp patterns). Different species of voles (Cricetidae, Rodentia, Glires) have evolved rooted and unrooted molars that have similar size and shape, providing alternative models for studying roots. We assembled a de novo genome of Myodes glareolus, a vole with high-crowned, rooted molars, and performed genomic and transcriptomic analyses in a broad phylogenetic context of Glires (rodents and lagomorphs) to assess differential selection and evolution in tooth forming genes. We identified 15 dental genes with changing synteny relationships and six dental genes undergoing positive selection across Glires, two of which were undergoing positive selection in species with unrooted molars, Dspp and Aqp1. Decreased expression of both genes in prairie voles with unrooted molars compared to bank voles supports the presence of positive selection and may underlie differences in root formation. Bulk transcriptomics analyses of embryonic molar development in bank voles also demonstrated conserved patterns of dental gene expression compared to mice, with species-specific variation likely related to developmental timing and morphological differences between mouse and vole molars. Our results support ongoing evolution of dental genes across Glires, revealing the complex evolutionary background of convergent evolution for ever-growing molars.

Predicting evolutionary patterns of mammalian teeth from development.

One motivation in the study of development is the discovery of mechanisms that may guide evolutionary change. Here we report how development governs relative size and number of cheek teeth, or molars, in the mouse. We constructed an inhibitory cascade model by experimentally uncovering the activator-inhibitor logic of sequential tooth development. The inhibitory cascade acts as a ratchet that determines molar size differences along the jaw, one effect being that the second molar always makes up one-third of total molar area. By using a macroevolutionary test, we demonstrate the success of the model in predicting dentition patterns found among murine rodent species with various diets, thereby providing an example of ecologically driven evolution along a developmentally favoured trajectory. In general, our work demonstrates how to construct and test developmental rules with evolutionary predictability in natural systems.

High-level similarity of dentitions in carnivorans and rodents.

The study of mammalian evolution depends greatly on understanding the evolution of teeth and the relationship of tooth shape to diet. Links between gross tooth shape, function and diet have been proposed since antiquity, stretching from Aristotle to Cuvier, Owen and Osborn. So far, however, the possibilities for exhaustive, quantitative comparisons between greatly different tooth shapes have been limited. Cat teeth and mouse teeth, for example, are fundamentally distinct in shape and structure as a result of independent evolutionary change over tens of millions of years. There is difficulty in establishing homology between their tooth components or in summarizing their tooth shapes, yet both carnivorans and rodents possess a comparable spectrum of dietary specializations from animals to plants. Here we introduce homology-free techniques to measure the phenotypic complexity of the three-dimensional shape of tooth crowns. In our geographic information systems (GIS) analysis of 441 teeth from 81 species of carnivorans and rodents, we show that the surface complexity of tooth crowns directly reflects the foods they consume. Moreover, the absolute values of dental complexity for individual dietary classes correspond between carnivorans and rodents, illustrating a high-level similarity between overall tooth shapes despite a lack of low-level similarity of specific tooth components. These results suggest that scale-independent forces have determined the high-level dental shape in lineages that are widely divergent in size, ecology and life history. This link between diet and phenotype will be useful for inferring the ecology of extinct species and illustrates the potential of fast-throughput, high-level analysis of the phenotype.