Exceptionally simple, rapidly replaced teeth in sauropod dinosaurs demonstrate a novel evolutionary strategy for herbivory in Late Jurassic ecosystems
BMC Ecology and Evolution volume 21, Article number: 202 (2021)
Dinosaurs dominated terrestrial environments for over 100 million years due in part to innovative feeding strategies. Although a range of dental adaptations was present in Late Jurassic dinosaurs, it is unclear whether dinosaur ecosystems exhibited patterns of tooth disparity and dietary correlation similar to those of modern amniotes, in which carnivores possess simple teeth and herbivores exhibit complex dentitions. To investigate these patterns, we quantified dental shape in Late Jurassic dinosaurs to test relationships between diet and dental complexity.
Here, we show that Late Jurassic dinosaurs exhibited a disparity of dental complexities on par with those of modern saurians. Theropods possess relatively simple teeth, in spite of the range of morphologies tested, and is consistent with their inferred carnivorous habits. Ornithischians, in contrast, have complex dentitions, corresponding to herbivorous habits. The dentitions of macronarian sauropods are similar to some ornithischians and living herbivorous squamates but slightly more complex than other sauropods. In particular, all diplodocoid sauropods investigated possess remarkably simple teeth. The existence of simple teeth in diplodocoids, however, contrasts with the pattern observed in nearly all known herbivores (living or extinct).
Sauropod dinosaurs exhibit a novel approach to herbivory not yet observed in other amniotes. We demonstrate that sauropod tooth complexity is related to tooth replacement rate rather than diet, which contrasts with the results from mammals and saurians. This relationship is unique to the sauropod clade, with ornithischians and theropods displaying the patterns observed in other groups. The decoupling of herbivory and tooth complexity paired with a correlation between complexity and replacement rate demonstrates a novel evolutionary strategy for plant consumption in sauropod dinosaurs.
Extinct dinosaurs exhibit a remarkable variety of dental morphologies ranging from narrow, peg-like teeth to tightly-packed, elaborate dental batteries, which are thought to directly reflect diet . Despite nearly two centuries of study and the discovery of novel tooth forms, our understanding of dinosaur dietary ecologies remains limited compared to other groups, such as mammals [1, 2]. Much of what is known relies on osteological comparisons to extant analogues, frequently noted in early discoveries , microwear analyses [2, 4,5,6,7,8,9], and the rare preservation of stomach contents [10,11,12,13,14]. The few studies that have applied quantitative methods to the study of dinosaur dentitions have yielded tantalizing insights. For example, an investigation of theropod dental disparity linked diet to the survival of neornithine birds after the end-Cretaceous mass extinction . Dental microwear analyses of sauropods strongly suggest that different clades employed differing feeding strategies, which may have permitted these large-bodied animals to live in the same environment [4, 16]. Another common quantitative dental assessment is the reconstruction of tooth replacement patterns and rates (all toothed dinosaurs exhibit polyphyodonty) [2, 17,18,19,20,21,22,23,24,25,26]. Histological investigations demonstrate that dental replacement rates vary substantially across the dinosaur clade, from every 777 days in the theropod Tyrannosaurus to 14 days in the sauropod Nigersaurus [2, 17, 19].
Dental topographic methods have yielded key insights into both modern and past ecosystems. For instance, using a quantitative measure of dental complexity—orientation patch count rotated (OPCR)—Evans et al.  demonstrated that phenotypic tooth complexity is related to diet in living carnivorans and rodents, with carnivores possessing simple teeth and herbivores displaying complex teeth. Using patterns of dental complexity, Wilson et al.  hypothesized that multituberculates ecologically radiated approximately 20 million years prior to the end-Cretaceous mass extinction. To a lesser extent, these methods have also been applied to dentigerous saurians (i.e. lepidosaurs and crocodilians), demonstrating that the relationship between diet and dental complexity originally observed in mammals is found across the amniote clade, in spite of major differences in development and replacement patterns . Despite their popularity and ubiquity, similar quantitative assessments of dental morphology have not been applied to dinosaurs, leaving broad questions unaddressed. For example, were dinosaurian ecosystems characterized by the same range of tooth complexity seen in modern saurian communities, and did dinosaurs exhibit similar correlations between diet and dental complexity as observed in modern amniotes?
The Late Jurassic represents an excellent system for testing hypotheses regarding dinosaur dentition. The paleoenvironment and taxic diversity of both fauna and flora are relatively well known, in part due to an extensive history of worldwide collection, but especially in that of western North America . During this epoch, dinosaurs attained a remarkable diversity, with theropods, sauropods, and ornithischians exhibiting a wide array of taxa and ecologies. In particular, the herbivorous sauropods achieved colossal sizes and an exceptional diversity, with many coexisting diplodocoids and macronarians. In spite of broad similarities in the sauropod bauplan, these two clades are characterized by marked morphological differences in both cranial and postcranial elements. Importantly, these skeletons sometimes include nearly complete cranial material with associated dental elements, which reveal a diverse range of tooth shapes, especially amongst herbivores [19, 31,32,33,34] (Fig. 1). The combination of skeletal and dental microwear differences strongly suggest a partitioning of plant resources. The abundance of dentigerous cranial material has also enabled the estimation of tooth replacement rates across the sauropod and theropod clade [19, 20, 22, 35]. These data, combined with the taxonomic diversity and range of dental morphologies, present an exceptional opportunity to test ecosystem-wide patterns during a relatively narrow time span.
Here, we seek to improve our understanding of dinosaur dietary ecology by testing two hypotheses using quantitative morphological assessments. First, we evaluate the relationship between diet and dental complexity in Late Jurassic dinosaurs using the OPCR method, specifically testing whether carnivores possessed simple teeth and herbivores exhibited complex dentitions. Second, we test the hypothesis that dental complexity values are related to tooth replacement rate. Our results shed light on key evolutionary innovations that allowed multiple herbivorous dinosaur clades to diversify and coexist.
Late Jurassic dinosaurs exhibit a wide range of dental complexities (Fig. 1). Average tooth complexity varies between 9.03 and 33.75 patches per tooth (PPT) when measured with a minimum polygon count threshold of three triangles (3 patch) and a model size of 1000 triangles (Figs. 2 and 3; Additional file 1, Additional file 1: Table S2). At a minimum polygon count threshold of five triangles (5 patch) and a model size of 1000 triangles, dental complexity patterns are consistent, but dental disparity decreases, with a range from 8.60 to 24 patches (Additional file 1: Table S2). This range is similar to those of extant saurians when measured using the same procedures despite differences in gross surface morphology .
Our results demonstrate dental complexities differ between the three major dinosaur clades when measuring the average value of each genus (Kruskal-Wallis: P < 0.01). Pairwise comparisons of 3 patch analyses show no significant difference between mean, median, and distribution of theropod and sauropod dental complexities (Student’s t-test; Mann-Whitney U; Kolmogorov–Smirnov test: P > 0.05, 9999 replicates) (Fig. 4). The mean, median, and distribution of theropod and sauropod tooth complexity are both significantly different from ornithischians (Student’s t-test, Mann-Whitney U, and Kolmogorov–Smirnov: P < 0.01, 9999 replicates). When the individual tooth complexities of macronarian and diplodocoid sauropods are analysed separately there is a significant difference between mean dental complexity of these groups (Student’s t-test: P < 0.01, t = 10.23, 9999 replicates). Macronarians in our sample are characterized by significantly higher complexities (14.28 PPT, 3 patch) than the narrow-crowned teeth of diplodocoids (9.74 PPT, 3 patch), as well as a greater range in complexities. The mean, median, and distribution of both sauropod clades are significantly different from theropods and ornithischians (Student’s t-test, Mann-Whitney U, and Kolmogorov–Smirnov: P < 0.01, 9999 replicates). The dental complexities of saurischians and ornithischians are significantly different for all statistical tests (Student’s t-test, Mann-Whitney U, and Kolmogorov–Smirnov: P < 0.01, 9999 replicates). At 5 patch analyses, results are similar to previous analyses with no significant difference between saurischian clades (Student’s t-test, Mann-Whitney U: P > 0.05). Surface complexity of measured theropods varied between 10.0 and 16.5 PPT (3 patch) when measured at a resolution of 1000 triangles (Figs. 2 and 3). In some cases (e.g. Morphotypes A and K), minor cracks or damaged enamel caused an increase in complexity. In other tooth morphotypes (O, Q, and R) the narrow cutting edge increased complexity by being broken up periodically, resulting in high complexities relative to other carnivorous saurians .
The seven measured sauropod genera display a greater range of dental complexities than those of theropods (Figs. 2 and 4; Additional file 1, Additional file 1: Table S2). The distribution of sauropod dental complexities is somewhat bimodal, with the broad-crowned teeth of macronarians exhibiting higher complexities than the narrow-crowned teeth of diplodocoids. The lowest value measured for a Late Jurassic dinosaur occurs in Tornieria (8 patches), a diplodocid from Tanzania . Other diplodocids (e.g. Apatosaurus, Diplodocus) also possess relatively simple teeth, between 9 and 11 PPT (3 patch) with dental models of 1000 triangles (Figs. 2, 4 and 5). These sauropods are characterized by narrow-crowned, peg-like teeth. The two macronarians included in this study, Brachiosaurus and Camarasaurus, exhibit the highest average dental complexities of our sauropod sample, 14.39 and 18.38 PPT (3 patch), respectively (Figs. 2 and 5). The teeth of these dinosaurs are characterized by a convex labial margin, concave lingual margin, and rugose wrinkled enamel . In both groups, tooth wear, including the development of large facets or the loss of enamel wrinkles, does not significantly impact dental complexity.
Ornithischians possess the most complex teeth of sampled Late Jurassic dinosaurs (Figs. 1, 2 and 4; Additional file 1). Nanosaurus is characterized by labiolingually compressed teeth with up to ten cusps (Additional file 1: Fig. S1C) and exhibits an average complexity of 33.75 PPT (3 patch) at a model size of 1,000 triangles. Fruitadens, one of the smallest known non-avian dinosaurs , possesses an average dental complexity of 23.63 PPT at the same resolution. Unworn teeth of Fruitadens can reach up to 36.5 PPT at the same parameters, rivalling the complexities of herbivorous iguanids . Other ornithischians have lower complexity values. Camptosaurus has an average complexity of 20.15 PPT (3 patch) when measured with a model of 1000 triangles. Gargoyleosaurus, an early ankylosaur, possesses simpler teeth, with a complexity value of 15.68 PPT when measured at the same resolution. Similar to Fruitadens, erupted and worn teeth are simpler than unworn dentitions. In one Camptosaurus specimen (MWC 2), noteworthy ridges and mesiodistal cusps are worn away, simplifying dental morphology and complexity values of this specimen, although wear facets of UMNH VP 16455 do not result in simplified complexities. This simplification contrasts with sauropods, in which tooth wear generally does not dramatically alter dental morphology and complexity.
Across measured dinosaurs, phylogenetic relationships exert less of an influence on OPCR values then predicted under Brownian motion, with K, a measure of phylogenetic signal, less than 1, although phylogenetic effects are still significant (Blomberg’s K = 0.82, P = 0.019) . Similarly, Pagel’s lambda is approximately 1 (λ = 1.04, P < 0.01), suggesting that OPCR values between taxa vary in proportion to their shared phylogenetic history. When measuring sauropods only, there is a significant phylogenetic signal in OPCR values, beyond what would be expected by a Brownian motion model of character evolution (Blomberg’s K = 1.62, P = 0.012). Within Sauropoda there is more variance in dental complexity between clades (i.e. diplodocoids and macronarians) than within them, suggesting that OPCR, at least in sauropods, is associated with phylogenetic relationships. In this case, diplodocoids possess simpler teeth whereas macronarians exhibit more complex dentitions than would be expected if phylogeny were not playing a role. In this reduced dataset, Pagel’s lambda is greater than 1, but this result is marginally non-significant (λ = 1.21, P = 0.059).
To better understand patterns of tooth shape in Late Jurassic dinosaurs, we tested for a correlation between replacement rate (measured in days) and dental complexity [2, 19, 20]. Across all sampled specimens, there is a positive correlation but it is not statistically significant (Pearson’s product-moment correlation: t = 1.99, P = 0.072, r = 0.51; Kendall’s rank correlation τ: z = 1.69, P = 0.09, τ = 0.37). For these data, however, replacement rates have been estimated for ornithischians based on the average of known replacement rates (96 days), which is likely an underestimate [17, 40]. When measuring sauropods only, with both replacement rates of Dicraeosaurus taken into account , there is a significant positive relationship between dental complexity and replacement rate (Pearson’s product-moment correlation: t = 2.97, P = 0.03; r = 0.77; Kendall’s rank correlation τ: z = 2.37, P = 0.02, τ = 0.69).
We then used PGLS and phylogenetic ANOVA to investigate patterns between OPCR and dental replacement rate in a phylogenetic framework. PGLS analyses recover differing results for different datasets. Tests of the complete dinosaur dataset do not recover replacement rate as a significant predictor of OPCR in either 3 or 5 patch datasets using both Brownian motion and Ornstein–Uhlenbeck models (P > 0.05), whereas the Pagel’s “lambda” model does result in a significant predictive relationship (P = 0.01; P = 0.02 for 3 and 5 patch datasets, respectively). A phylogenetic ANOVA found a non-significant relationship in both OPCR datasets (P > 0.05). When investigating only saurischians, in which tooth replacement data are better known, the replacement rate is a significant predictor of OPCR under Brownian motion (P = 0.02 for both 3 and 5 patch datasets), but not Ornstein–Uhlenbeck, Pagel’s “lambda”, or a phylogenetic ANOVA (P > 0.05). This contrasts with the results of sauropods only, in which dental complexity is not a function of replacement rate under Brownian motion, Pagel’s “lambda”, and phylogenetic ANOVA (P > 0.05), but there is a positive, significant relationship with the Ornstein-Uhlenbeck model of evolution (3 patch dataset: P = 0.04). In sauropod-only analyses, if Dicraeosaurus is treated as two separate taxa, because of the two replacement rates for upper and lower dentitions, all PGLS and phylogenetic ANOVA return a significant (P < 0.05) positive relationship between tooth complexity and replacement rate for 3 and 5 patch datasets.
Here, we demonstrate a unique relationship between phenotypic tooth complexity and dietary ecology in Late Jurassic dinosaurs. Because independent evidence of general feeding ecology (e.g. herbivory, carnivory) is known for a variety of dinosaur taxa, diets of these clades can be reconstructed with relative confidence allowing for the testing of relationships between diet and dental complexity. In measured ornithischians and theropods, the previously recovered pattern is maintained; carnivorous theropods have simple teeth, whereas herbivorous ornithischians display the most complex dentitions in sampled dinosaurs (Figs. 1, 2, 4 and 5). Ornithischians, like extant herbivorous saurians , express the greatest range of complexities (Fig. 4), reflecting their heterodont condition. In both living saurians and mammals, there is a positive relationship between dental complexity and the amount of plant matter consumed [27, 29, 41]. Hypercarnivores and carnivores (i.e. animals that primarily consume animal material but also small amounts of plants) possess simple teeth (i.e. low dental complexity). Herbivores exhibit the highest dental complexity values, reflecting additional cusps, crenulations, and grinding surfaces that enlarge the dental surface area and help break down plant material prior to chemical digestion. Herbivorous sauropods, however, complicate this otherwise straightforward pattern by possessing relatively simple dentitions. In fact, with their extremely simple teeth, diplodocoid sauropods appear to exhibit a novel strategy of plant consumption (Figs. 2, 5 and 6). Therefore, our hypothesis that dinosaurs that lived during the Late Jurassic exhibit modern patterns of dental complexity and diet is falsified, at least for the diplodocoids included in this study.
In spite of noteworthy differences in morphology, the dental complexity of diplodocoids closely resembles that of carnivorous theropods, and in some cases is even simpler. The narrow peg-like teeth of this group lack prominent cusps that add complexity to the teeth of other herbivores, both living and extinct [29, 42]. Sauropod herbivory is not in doubt (but see Stokes ), as the group is characterized by a number of features that are related to this diet, including large body size and wide rib cages . The limited role these teeth play in food processing combined with the absence of precise dental occlusion may be expected to result in simpler dentitions in diplodocoids, but these characteristics also are observed in many groups of extant lepidosaurs that are typified by complex dentitions. Together, this suggests that sauropods, especially diplodocoids, are unique in currently sampled taxa, living and extinct, in having low dental complexity values associated with an herbivorous diet.
To investigate this pattern further, we looked at the tooth replacement rates of sauropods [2, 19,20,21, 45]. Replacement rates of dinosaurs are calculated by counting the differences in incremental lines of von Ebner (previously shown to be deposited daily in living crocodilians) between successive tooth generations . There is a marked disparity in dinosaur tooth replacement rates. For example, the narrow peg-like teeth of Diplodocus were replaced at a relatively rapid rate, each tooth every 35 days approximately, whereas the replacement rate is slower in macronarians, which often exhibit a wider and spatulate morphology (e.g. 62 days for Camarasaurus) . In herbivorous taxa with well-documented replacement rates, there is a clear correlation between dental complexity and tooth replacement rate (r = 0.77; τ = 0.69; P < 0.05; Fig. 6). Sauropods with rapid replacement rates tend to have simpler teeth, whereas those that have slower rates possess more complex dentitions (Fig. 6). Analyses that take phylogeny into account demonstrate that this relationship is significant in sauropods. Diplodocoids have the lowest OPCR values of herbivorous saurischians and the fastest rates of tooth replacement. Macronarian sauropods, in contrast, exhibit higher average complexities, longer replacement rates, thicker enamel, and a greater tooth volume . Interestingly, the pattern between complexity and replacement rate appears to extend even to a single taxon. The premaxillary and maxillary teeth of Dicraeosaurus were replaced at a faster rate than those of the dentary, and there is a significant difference in complexity of those elements (10.32 and 12.73 PPT, respectively; P < 0.01). Unfortunately, precise replacement rate data remain unknown for the ornithischians included in this study, but if they are similar to that of sampled relatives (e.g. Edmontonia, Triceratops, Maiasaura) they are longer than sauropods (279–381 days for adult taxa), which would also be consistent with the distinct correlation between higher dental complexity and slower replacement rate .
Narrow, peg-like teeth evolved independently in more than one sauropod lineage [19, 46]. After the disappearance of diplodocids at the beginning of the Early Cretaceous, surviving macronarians developed dentitions that resemble those of diplodocoids (i.e. a similar slenderness index) . In particular, titanosaurs, such as Antarctosaurus, Nemegtosaurus, and Rapetosaurus, developed extremely narrow tooth crowns . Coupled with this change in morphology was a notable increase in the number of replacement teeth (e.g. Museo Provincial Carlos Ameghino, MPCA-79), an indication of high replacement rates, although titanosaur replacement rates remain undocumented histologically [19, 46, 47]. The independent evolution of narrow-crowned teeth strongly implies a simple complexity in titanosaurs. If coupled with high replacement rates (as is hypothesized in MPCA-79 ), this may indicate convergence in herbivory strategies in distantly related sauropod clades. Thus, diplodocoids, and potentially titanosaurs, may have evolved a method of plant consumption and processing not yet documented in either extinct or extant amniotes. Instead of complex teeth that break down plant material, these sauropods utilized extremely simple teeth with narrow bands of enamel that were quickly worn and rapidly replaced. This style of herbivory may have been greatly beneficial because it had a dual advantage for a long-necked plant-eater: enabling resource partitioning with contemporaneous herbivores (other sauropods, ornithischians, and additional amniotes) while also allowing for a lightened skull . Additionally, simple-toothed herbivores independently occur in the Late Jurassic and Late Cretaceous. There was a remarkable shift in plant community composition during this time interval, from gymnosperm-dominated communities to those in which angiosperms were the principal floral component of terrestrial environments . Thus, the presence of this morphology in both environments suggests this approach is not associated with a particular dietary specialization but instead represents an adaptation for a broader feeding strategy.
Previous research interpreted these morphological specializations with low-browsing feeding. The diet of diplodocoids with narrow-crowned teeth has been hypothesized to largely consist of abrasive plants, although the exact food differs [4, 46, 49,50,51]. An analysis of snout shape and microwear of diplodocoids found evidence for multiple feeding strategies, ranging from low to medium browse and both selective and non-selective diets . The plant genus Equisetum (i.e. horsetail) was suggested to be an excellent potential food source in the Morrison Formation environment, in fact it yields energy levels that surpass extant grasses, but a noteworthy downside is that its surfaces are rich in silica [50, 52]. Simple teeth that are rapidly replaced due to extensive wear may be an adaptation to an exceptionally abrasive diet , but this stands in stark contrast to the complex dentitions found in the vast majority of herbivorous mammals, squamates, and extinct crocodyliforms [27, 29, 41, 42]. Complex teeth provide a greater surface area that breaks down plant material prior to ingestion, facilitating the absorption of nutrients. The evolution of a more abrasive diet and the concurrent change in dental complexity has previously been investigated in equids . In equids, as abrasive foods became a greater proportion of diet and body size increased, dental complexity grew through the addition of increasingly fine features . Here, we show the opposite trend in large-bodied taxa that continuously replaced their teeth. Instead of high volume, complex dentitions, browsing sauropods develop simple teeth that were rapidly worn and subsequently replaced. Simple and narrow dentitions, with a lower enamel and dentin volume , that could be rapidly replaced would have permitted the consumption of abrasive plant material that may otherwise have been detrimental to taxa with a more permanent dentition. Simultaneously, smaller teeth would have been advantageous to a retaining a light skull and long neck, which has been hypothesized to be a key factor in sauropod gigantism and success . Interestingly, simple teeth are also known in two species of living herbivorous squamates, the skinks Egernia stokesii and Tiliqua rugosa, but these taxa are closely related to non-herbivorous species [29, 55]. Additionally, replacement rates from skinks remain unknown, hindering our ability to test for the pattern between simple teeth and rapid replacement rates in living herbivorous saurians.
In contrast to the novel patterns observed in sauropods, theropod dental complexity is consistent with that of other carnivorous saurians. This dental simplicity of sampled theropods reinforces previous evidence, including serrations and sharp cusps, that support the view of carnivory across this group during the Late Jurassic. We sampled 12 tooth morphotypes and found that theropods display the lowest disparity of sampled clades in spite of the diverse collection of morphologies exhibited (Fig. 3). The majority of theropod teeth are relatively simple, approximately 10–12 PPT for 3 patch analyses. In part, this is due to fine features, such as serrations, not being detected by OPCR analyses. Although these structures are critical for the consumption of animal material, they are lost during model preparation that permit comparisons between different sized taxa or scanners . Finer morphological elements were also not detected in previous analyses of both extant saurians and extinct crocodyliforms and so this absence is not unique to theropod dinosaurs [29, 42].
Despite the overall low OPCR values of carnivorous theropods, there is a noteworthy range in complexities (i.e. 8.25–16.5 PPT) (Fig. 5). This result demonstrates that dental complexity differences reflect, at least to some degree, the morphological disparity that existed in Late Jurassic theropods. In living saurians, the carnivore dietary category encompasses a broad range of diets . The range in morphology and complexity, especially when coupled with body size variation, likely reflects a similar pattern with theropods, with taxa potentially specializing on different faunal resources (i.e. niche partitioning) . The long-held inference of a uniform diet across theropods of the Late Jurassic should be tested with greater rigor, as the substantial disparity in dental size, shape, and complexity may reflect an otherwise unacknowledged ecological range in this diverse group. In particular, the variation in tooth complexity of theropods present during the narrow temporal range of the Late Jurassic may be a precursor of the later ecological divergence found in Cretaceous theropods. Studies of extant avians show that ecological transitions to omnivory from carnivory are much more common than directly from carnivory to herbivory . Therefore, although Late Jurassic theropods in general possessed simple teeth, the variation in qualitative features, as well as complexity, observed in this clade likely reflect a dietary diversity already in place prior to the more pronounced ecological shifts that occur during the Cretaceous.
The radiation of theropods into novel ecological roles also coincided with numerous shifts in cranial and dental morphology . Although some groups, such as therizinosaurs, evolve heterodonty and even multi-cusped teeth that superficially resemble the morphology of extant iguanids , many others possess simplified, conical-shaped teeth . In some cases, herbivorous habits are associated with tooth loss in theropods . This initial simplification and subsequent loss of teeth, especially in the distal jaw, may represent a similar adaptation towards herbivory as that of diplodocoids. Although this remains untested, it would demonstrate an intriguing pattern in which multiple dinosaur clades independently simplified dentitions for an herbivorous lifestyle, which would stand in stark contrast to the patterns observed in extant herbivores, saurian and mammal alike.
This work demonstrates the remarkable potential of quantitative dental analyses to illuminate dinosaurian palaeoecology. In particular, OPCR analyses have the capacity to address key macroevolutionary questions, measure and describe patterns within and across dinosaur communities around the world, and even allow for direct comparison to modern faunas. Here, we find clear patterns between dinosaur dental complexity and diet. In some cases, carnivorous theropods possess simple teeth and herbivorous ornithischians display more complex tooth shapes, reflecting broad patterns observed in both living and extinct amniotes. Surprisingly, theropods display a wide dental disparity, despite overall low values of tooth complexity. This variation may not only be related to taxonomic diversity but also to dietary differences among Late Jurassic theropods that have remained otherwise undetected. This ecological disparity may provide early clues to the dietary breadth that would characterize the later evolution of this clade during the Cretaceous. In the case of sauropods, however, the widely documented pattern between tooth complexity and diet breaks down, representing a novel strategy of herbivory. Dental complexity is correlated to tooth replacement rate in sauropod dinosaurs, even when phylogeny is taken into account. When combined with other morphological features of the diplodocoid clade, this strategy likely represents an adaptation to massive body sizes and living alongside a diverse fauna of megaherbivores (Additional file 2).
Specimens and sampling
We investigated dental complexities from a broad range of Late Jurassic dinosaurs (Additional file 1: Fig. S1, Table S1). This dataset includes four ornithischians, seven sauropods, two theropods and ten theropod tooth morphotypes, capturing much of the dental disparity of Late Jurassic dinosaurs. We measured a combination of isolated teeth and those preserved within jaws, the latter comprising the majority of the dataset. Although the majority of our dataset comprises specimens from the Morrison Formation of western North America, we include key samples from Europe and Africa that fill in important gaps in species diversity and dental morphologies . Following previous work, we sample teeth along the entire dentigerous element [29, 42].
To model the influence of phylogenetic relationships, we built a time-calibrated phylogenetic tree based on a combination of large, cladistic datasets that analysed broad relationships between major dinosaur groups  and those that focused on specific clades [62,63,64]. While phylogenetic uncertainty and some conflict exist between studies of dinosaur relationships, the phylogenetic positions of taxa sampled in this study are in broad agreement across recent studies. Branch lengths exert considerable influence on character evolution  and unfortunately no definitive time-calibrated phylogeny exists for non-avian dinosaurs. To calibrate our tree, we used estimates based on data presented in Benson et al. , as well as values from Whitlock  and D’Emic  for sauropod relationships. We accounted for differences in the timing of clade origination by manually modifying branch lengths to characterize the effect of different estimates, which did not result in noteworthy differences in statistical analyses. All taxa included in this study are from a relatively narrow temporal range; the ages of the Morrison Formation and the Middle Dinosaur Member of the Tendaguru Formation were estimated to be 152 MA (the Kimmeridgian-Tithonian boundary) based on a narrow range of ages estimated from microfossil and 40Ar/39Ar dating [66, 67]. Theropod tooth morphotypes were not included in phylogenetic analyses. Phylogenies that were used in statistical analyses (see below) were written in Newick format in R and generated using the read.tree function from the ‘ape’ package (version 5.4-1) [68, 69].
We generated 3D models of dinosaur dentition from a combination of microCT and CT scans. CT data of sauropods, theropods, and Fruitadens were taken from previously published sources, whereas Camptosaurus, Gargoyleosaurus, and Nanosaurus were scanned at the University of Southern California Molecular Imaging Center. Individual scanning parameters are available in their original publications and the Additional file (Additional file 1: Table S1). We generated 3D models of both worn and unworn dentitions in Avizo Lite (Version 2020.1; Thermo Fisher Scientific). Following procedures outlined by previous research on mammal and saurian dentitions [36, 70], we standardized surfaces to 10,000 triangles (± 1 triangles) using MeshLab’s  ‘Simplification: Quadratic Edge Collapse Decimation’ tool. The dentitions were then smoothed using the Laplacian smooth function with three steps followed by a second down-sample to 1000 triangles. Dental complexity measurements can be influenced by factors not directly related to surface morphology, including tooth size, digitization method, and number of triangles in each model. Larger teeth produce higher resolution digital models, which, in turn, contain a greater amount of morphological information, a result of more triangles in each model. Thus, larger teeth and models with a greater number of triangles are more likely to result in a higher OPCR value. Similarly, CT scans will likely yield lower OPCR values than microCT scans if scanning the same specimen because the latter captures more information per unit area. Standardizing dental models to a consistent resolution/triangle count largely removes these effects as well as additional minor damage that may have occurred during the fossilization process [36, 70]. In some cases, additional morphological information can be lost during this standardization process, but previous work has demonstrated this did not significantly alter broad complexity patterns . 3D dental models used in this research are available for download at https://www.morphosource.org/projects/000345163.
OPCR and statistical analysis
We measured tooth complexity using the 3D orientation patch count rotated (OPCR) method, which quantitatively assesses dental morphology and allows for teeth with no homologous landmarks to be directly compared [27, 28, 72]. This method assigns a cardinal or ordinal direction to each triangle, groups contiguous triangles with the same orientation together into patches, and then sums the patches. Following this, models are rotated 5.625° and patch counts are recalculated. This latter step is repeated seven times, with the final complexity measurement averaged over a total of eight trials. This method provides a numerical representation of surface shape that is robust to minor differences in model orientation. In sampled living animals, dental complexity tends to be simple in carnivores and complex in herbivores, with dietary generalists frequently falling between the two end-members [27, 29, 41]. We used MorphoTester, a free open-source dental topology software, to measure the dental complexity of each tooth . This program, written originally in the Python programming language, analyses full 3D polygon meshes. This contrasts with a previous application of OPCR, which utilizes Surfer Manipulator and raster-based digital elevation models [27,28,29, 73]. These models differ from 3D polygon meshes in that they associate X and Y values with only a single Z value (i.e. 2.5D) and these two sets of OPCR data are not directly comparable . Surfer Manipulator captures topography of the occlusal margin, the dental surface that interacts with food, whereas MorphoTester measures the entire dental surface. This additional measurement by MorphoTester systematically results in higher complexity measurements relative to Surfer Manipulator [36, 72], however, MorphoTester is freely available and is better able to batch process large numbers of data files. PLY files were uploaded to MorphoTester and evaluated using a minimum patch size of three and five triangles to facilitate comparisons with previous work . Results that utilize minimum patch counts of three triangles capture finer details, such as small cusps, but are also more likely to detect dental damage, whereas minimum patch counts of five triangles are more robust to these issues, but may not identify fine details. We report both sets of results below (Additional file 1: Table S2, Additional file 2). Following previous research on saurians, we averaged the complexity of all teeth for a single taxon (OPCRavg).
We performed statistical analyses in RStudio (version 1.3.1073; R version 4.0.2) and PAST4 [74, 75]. Replacement rates were taken from previous studies on sauropods [19,20,21, 35, 45] and theropods . The dental replacement rate of Suuwassea was estimated based on the average rate of its closest relative Dicraeosaurus, whereas Apatosaurus was estimated based on a similar tooth replacement count between these specimens and Nigersaurus [19, 21]. The significance of the correlation between dental complexity and tooth replacement rate was evaluated using Pearson’s product-moment correlation and Kendall’s rank correlation τ by the cor.test function from the ‘stats’ package (version 4.0). We tested for phylogenetic signal and its significance in OPCR and replacement rate results using the K-statistic and Pagel’s lambda, which were calculated using the phylosig function from the ‘phytools’ package (version 0.7-47) . We used multiple methods implemented in R to assess the relationship between replacement rate and dental complexity in a phylogenetic framework. A phylogenetic ANOVA was preformed using the procD.pgls function (SS.type = I, effect.type = F, iter = 9999) from the package ‘geomorph’ (version 3.3.1), in which average dental complexity of each genus was the response variable and replacement rate in days was the independent variable . We also ran phylogenetic generalized least squares (PGLS) using the ‘ape’ package and its dependencies [69, 78] (version 5.4-1). We first created correlation structures using Pagel’s lambda (corPagel), Brownian motion (corBrownian), and Ornstein–Uhlenbeck (corMartins) models of trait covariance [68, 79]. Following this, we used the gls function (‘nlme’ package, version 3.1–148 ) to perform generalized least squares analysis where dental complexity was a function of replacement rate .
Availability of data and materials
3D models dental models used in this study are available on https://www.morphosource.org/projects/000345163.
Weishampel DB, Dodson P, Osmólska H. The Dinosauria. 2nd ed. Berkeley: University of California Press; 2004.
Sereno PC, Wilson JA, Witmer LM, Whitlock JA, Maga A, Ide O, et al. Structural extremes in a Cretaceous dinosaur. PLoS One. 2007;2:e1230.
Mantell G. Notice on the Iguanodon, a newly discovered fossil reptile, from the sandstone of Tilgate Forest, in Sussex. Philos Trans R Soc London. 1825;115:179–86.
Whitlock JA. Inferences of diplodocoid (Sauropoda: Dinosauria) feeding behavior from snout shape and microwear analyses. PLoS One. 2011;6:e18304.
Williams VS, Barrett PM, Purnell MA. Quantitative analysis of dental microwear in hadrosaurid dinosaurs, and the implications for hypotheses of jaw mechanics and feeding. Proc Natl Acad Sci U S A. 2009;106:11194–11199.
Mallon JC, Anderson JS. The functional and palaeoecological implications of tooth morphology and wear for the megaherbivorous dinosaurs from the Dinosaur Park Formation (Upper Campanian) of Alberta, Canada. PLoS One. 2014;9:e98605.
Varriale FJ. Dental microwear reveals mammal-like chewing in the neoceratopsian dinosaur Leptoceratops gracilis. PeerJ. 2016;4:e2132.
Virág A, Ősi A. Morphometry, microstructure, and wear pattern of neornithischian dinosaur teeth from the Upper Cretaceous Iharkút Locality (Hungary). Anat Rec. 2017;300:1439–1463.
Torices A, Wilkinson R, Arbour VM, Ruiz-Omeñaca JI, Currie PJ. Puncture-and-pull biomechanics in the teeth of predatory coelurosaurian dinosaurs. Curr Biol. 2018;28:1467–1474.
Chen PJ, Dong ZM, Zhen SN. An exceptionally well-preserved theropod dinosaur from the Yixian Formation of China. Nature. 1998;391:147–152.
Varricchio DJ. Gut contents from a Cretaceous tyrannosaurid: implications for theropod dinosaur digestive tracts. J Paleontol. 2001;75:401–406.
Nesbitt SJ, Turner AH, Erickson GM, Norell MA. Prey choice and cannibalistic behaviour in the theropod Coelophysis. Biol Lett. 2006;2:611–614.
Salgado L, Canudo JI, Garrido AC, Moreno-Azanza M, Martínez LCA, Coria RA, et al. A new primitive Neornithischian dinosaur from the Jurassic of Patagonia with gut contents. Sci Rep. 2017;7:42778.
Brown CM, Greenwood DR, Kalyniuk JE, Braman DR, Henderson DM, Greenwood CL, et al. Dietary palaeoecology of an Early Cretaceous armoured dinosaur (Ornithischia; Nodosauridae) based on floral analysis of stomach contents. R Soc Open Sci. 2020;7:200305.
Larson DW, Brown CM, Evans DC. Dental disparity and ecological stability in bird-like dinosaurs prior to the end-Cretaceous mass extinction. Curr Biol. 2016;26:1325–1333.
Fiorillo AR. Dental microwear patterns of the sauropod dinosaurs Camarasaurus and Diplodocus: Evidence for resource partitioning in the Late Jurassic of North America. Hist Biol. 1998;13:1–16.
Erickson GM. Incremental lines of von Ebner in dinosaurs and the assessment of tooth replacement rates using growth line counts. Proc Natl Acad Sci U S A. 1996;93:14623–7.
Wiersma K, Sander PM. The dentition of a well-preserved specimen of Camarasaurus sp.: implications for function, tooth replacement, soft part reconstruction, and food intake. Palaontologische Zeitschrift. 2017;91:145–161.
D’Emic MD, Whitlock JA, Smith KM, Fisher DC, Wilson JA. Evolution of high tooth replacement rates in sauropod dinosaurs. PLoS One. 2013;8:e69235.
Schwarz D, Kosch JCD, Fritsch G, Hildebrandt T. Dentition and tooth replacement of Dicraeosaurus hansemanni (Dinosauria, Sauropoda, Diplodocoidea) from the Tendaguru Formation of Tanzania. J Vertebr Paleontol. 2015;35:e1008134.
McHugh JB. Evidence for niche partitioning among ground-height browsing sauropods from the Upper Jurassic Morrison Formation of North America. Geol Intermt West. 2018;5:95–103.
D’Emic MD, O’Connor PM, Pascucci TR, Gavras JN, Mardakhayava E, Lund EK. Evolution of high tooth replacement rates in theropod dinosaurs. PLoS One. 2019;14:e0224734.
LeBlanc ARH, Reisz RR, Evans DC, Bailleul AM. Ontogeny reveals function and evolution of the hadrosaurid dinosaur dental battery. BMC Evol Biol. 2016;16:152.
Becerra MG, Pol D, Whitlock JA, Porro LB. Tooth replacement in Manidens condorensis: baseline study to address the replacement pattern in dentitions of early ornithischians. Pap Palaeontol. 2021;7:1167–93.
Hanai T, Tsuihiji T. Description of tooth ontogeny and replacement patterns in a juvenile Tarbosaurus bataar (Dinosauria: Theropoda) using CT-scan data. Anat Rec. 2019;302:1210–1225.
Sattler F, Schwarz D. Tooth replacement in a specimen of Tyrannosaurus rex (Dinosauria, Theropoda) from the Hell Creek Formation (Maastrichtian), Montana. Hist Biol. 2019;33:949–972.
Evans AR, Wilson GP, Fortelius M, Jernvall J. High-level similarity of dentitions in carnivorans and rodents. Nature. 2007;445:78–81.
Wilson GP, Evans AR, Corfe IJ, Smits PD, Fortelius M, Jernvall J. Adaptive radiation of multituberculate mammals before the extinction of dinosaurs. Nature. 2012;483:457–60.
Melstrom KM. The relationship between diet and tooth complexity in living dentigerous saurians. J Morphol. 2017;278:500–522.
Foster JR. Jurassic West. Second. Bloomington, Indiana: Indiana University Press; 2020.
Holland WJ. The skull of Diplodocus. Mem Carnegie Museum. 1924;IX:379–403.
Gilmore CW. A nearly complete articulated skeleton of Camarasaurus, a saurischian dinosaur from the Dinosaur National Monument. Ann Carnegie Museum. 1925;10:347–384.
Carpenter K, Miles C, Cloward K. Skull of a Jurassic ankylosaur (Dinosauria). Nature. 1998;393:782–783.
Whitlock JA, Wilson JA, Lamanna MC. Description of a nearly complete juvenile skull of Diplodocus. J Vertebr Paleontol. 2010;30:442–457.
D’Emic MD, Carrano MT. Redescription of brachiosaurid sauropod dinosaur material from the Upper Jurassic Morrison Formation, Colorado, USA. Anat Rec. 2020;303:732–758.
Melstrom KM, Wistort Z. The application of dental complexity metrics on extant saurians. Herpetologica.
Remes K. Revision of the Tendaguru sauropod dinosaur Tornieria africana (Fraas) and its relevance for sauropod paleobiogeography. J Vertebr Paleontol. 2006;26:651–669.
Butler RJ, Galton PM, Porro LB, Chiappe LM, Henderson DM, Erickson GM. Lower limits of ornithischian dinosaur body size inferred from a new upper Jurassic heterodontosaurid from North America. Proc R Soc B Biol Sci. 2010;277:375–381.
Blomberg SP, Garland, Jr. T, Ives AR. Testing for phylogenetic signal in comparative data: Behavioral traits are more labile. Evolution (N Y). 2003;57:717–745.
Sokolskyi TL, Zanno LE, Kosch JCD. Unusual tooth replacement in a new Cenomanian iguanodontian from the Mussentuchit Member of the Cedar Mountain Formation. J Vertebr Paleontol. 2019;2019:195–6.
Pineda-Munoz S, Lazagabaster IA, Alroy J, Evans AR. Inferring diet from dental morphology in terrestrial mammals. Methods Ecol Evol. 2016;8:481–491.
Melstrom KM, Irmis RB. Repeated evolution of herbivorous crocodyliforms during the age of dinosaurs. Curr Biol. 2019;29:2389–2395.
Stokes WL. Fossilized stomach contents of a sauropod dinosaur. Science. 1964;143:576–577.
King GM. Reptiles and herbivory. London: Chapman and Hall; 1996.
Sattler F. Tooth replacement of the sauropod dinosaur Tornieria africana (Fraas) from Tendaguru (Late Jurassic, Tanzania). Bachelor thesis, p. 1–41, 2014.
Chure D, Britt BB, Whitlock JA, Wilson JA. First complete sauropod dinosaur skull from the Cretaceous of the Americas and the evolution of sauropod dentition. Naturwissenschaften. 2010;97:379–391.
Coria RA, Chiappe LM. Tooth replacement in a sauropod premaxilla from the upper Cretaceous of Patagonia, Argentina. Ameghiniana. 2001;38:463–466.
Butler RJ, Barrett PM, Kenrick P, Penn MG. Diversity patterns amongst herbivorous dinosaurs and plants during the Cretaceous: Implications for hypotheses of dinosaur/angiosperm co-evolution. J Evol Biol. 2009;22:446–459.
Upchurch P, Barrett PM. Phylogenetic and taxic perspectives on sauropod diversity. In: Wilson JA, Curry Rogers K, editors. The sauropods evolution: and paleobiology. Berkeley: University of California Press; 2005. p. 104–24.
Gee CT. Dietary options for the sauropod dinosaurs from an integrated botanical and paleobotanical perspective. In: Klein N, Remes K, Gee CT, Sander PM, editors. Biology of the sauropod dinosaurs: Understanding the life of giants. Bloomington: Indiana University Press; 2011. p. 34–56.
Woodruff DC, Carr TD, Storrs GW, Waskow K, Scannella JB, Nordén KK, et al. The smallest diplodocid skull reveals cranial ontogeny and growth-related dietary changes in the largest dinosaurs. Sci Rep. 2018;8:1–12.
Sander PM, Gee CT, Hummel J, Clauss M. Mesozoic plants and dinosaur herbivory. In: Gee CT, editor. Plants in Mesozoic time: morphological innovations, phylogeny, ecosystems. Bloomington: Indiana University Press; 2010. p. 331–359.
Evans AR, Janis CM. The evolution of high dental complexity in the horse lineage. Ann Zool Fennici. 2014;51 April:73–79.
Sander PM, Christian A, Clauss M, Fechner R, Gee CT, Griebeler EM, et al. Biology of the sauropod dinosaurs: The evolution of gigantism. Biol Rev. 2011;86:117–155.
Cooper Jr WE, Vitt LJ. Distribution, extent, and evolution of plant consumption by lizards. J Zool. 2002;257:487–517.
Frederickson JA, Engel MH, Cifelli RL. Niche partitioning in theropod dinosaurs: Diet and habitat preference in predators from the uppermost Cedar Mountain Formation (Utah, U.S.A.). Sci Rep. 2018;8:17872.
Burin G, Kissling WD, Guimarães PR, Şekercioglu CH, Quental TB. Omnivory in birds is a macroevolutionary sink. Nat Commun. 2016;7:11250.
Zanno LE, Makovicky PJ. Herbivorous ecomorphology and specialization patterns in theropod dinosaur evolution. Proc Natl Acad Sci U S A. 2011;108:232–237.
Zanno LE. Osteology of Falcarius utahensis (Dinosauria: Theropoda): Characterizing the anatomy of basal therizinosaurs. Zool J Linn Soc. 2010;158:196–230.
Gerke O, Wings O. Multivariate and cladistic analyses of isolated teeth reveal sympatry of theropod dinosaurs in the Late Jurassic of Northern Germany. PLoS One. 2016;11:e0158334.
Benson RBJ, Campione NE, Carrano MT, Mannion PD, Sullivan C, Upchurch P, et al. Rates of dinosaur body mass evolution indicate 170 Million years of sustained ecological innovation on the avian stem lineage. PLoS Biol. 2014;12:e1001853.
Whitlock JA. A phylogenetic analysis of Diplodocoidea (Saurischia: Sauropoda). Zool J Linn Soc. 2011;161:872–915.
D’Emic MD. The early evolution of titanosauriform sauropod dinosaurs. Zool J Linn Soc. 2012;166:624–71.
Carpenter K, Galton PM. A photo documentation of bipedal ornithischian dinosaurs from the Upper Jurassic Morrison Formation, USA. Geol Intermt West. 2018;5:167–207.
Campione NE, Barrett PM, Evans DC. On the ancestry of feathers in Mesozoic dinosaurs. In: Foth C, RauhutOWM, editors. The evolution of feathers: From their origin to the present. Springer; 2020. p. 213–243.
Bussert R, Heinrich WD, Aberhan M. The Tendaguru Formation (Late Jurassic to Early Cretaceous, southern Tanzania): Definition, palaeoenvironments, and sequence stratigraphy. Foss Rec. 2009;12:141–174.
Trujillo KC, Kowallis BJ. Recalibrated legacy 40Ar/39Ar ages for the Upper Jurassic Morrison Formation, Western Interior, USA. Geol Intermt West. 2015;2:1–8.
Paradis E. Analysis of phylogenetics and evolution with R. 2nd ed. New York: Springer; 2012.
Paradis E, Claude J, Strimmer K. APE: analyses of phylogenetics and evolution in R language. Bioinformatics. 2004;20:289–290.
Berthaume MA, Winchester J, Kupczik K. Effects of cropping, smoothing, triangle count, and mesh resolution on 6 dental topographic metrics. PLoS One. 2019;14:e0216229.
Cignoni P, Callieri M, Corsini M, Dellepiane M, Ganovelli F, Ranzuglia G. MeshLab: An open-source mesh processing tool. Eurographics Ital Chapter Conf. 2008;2008:129–36.
Winchester JM. MorphoTester: An open source application for morphological topographic analysis. PLoS One. 2016;11:e0147649.
Santana SE, Strait S, Dumont ER. The better to eat you with: Functional correlates of tooth structure in bats. Funct Ecol. 2011;25:839–847.
Hammer Ø, Harper DA, Ryan PD. PAST: Paleontological statistics software package for education and data analysis. Paleontol Electron. 2001;4:1–9.
R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. 2020. https://www.R-project.org/.
Revell LJ. phytools: An R package for phylogenetic comparative biology (and other things). Methods Ecol Evol. 2012;3:217–223.
Adams DC, Otárola-Castillo E. geomorph: an R package for the collection and analysis of geometric morphometric shape data. Methods Ecol Evol. 2013;4:393–9.
Felsenstein J. Phylogenies and the comparative method. Am Nat. 1985;125:1–15.
Revell LJ. Phylogenetic signal and linear regression on species data. Methods Ecol Evol. 2010;1:319–329.
Pinheiro J, Bates D, DebRoy S, Sarkar D, Heisterkamp S, Van Willigen B. Package ‘nlme’. 2021. https://0-doi-org.brum.beds.ac.uk/https://cran.r-project.org/web/packages/nlme/nlme.pdf. Accessed 3 May 2021.
We are extraordinarily grateful to those who were willing to share CT data during the Covid-19 pandemic: R. Butler, M. D’Emic, J. McHugh, D. Schwarz, and J. Whitlock. This work would not have been possible without those committed to open science. We are grateful to C. Beno, A. and T. Frisbie, R. Irmis, C. Levitt-Bussian, K. MacKenzie, J. McHugh, E. Melstrom, J. Sertich, and N. Toth for facilitating the borrowing of dinosaur jaw material. We thank H. O’Brien for statistics explanations and T. Skorka for CT assistance. We sincerely appreciate the editing of D. Provete and the reviews of N. Campione, J. Whitlock, and an anonymous reviewer whose excellent advice helped improve the quality of this manuscript.
This work and the Postdoctoral Fellowship of K.M.M. are supported by NHMLAC Dinosaur Institute research funds.
Ethics approval and consent to participate
Consent for publication
The authors have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Dental complexity of the major dinosaur groups, with sauropods broken into their two primary clades included in this study: diplodocoids and macronarians. . Table S1. Information on the Late Jurassic dinosaurs measured. Table S2. Average dental complexity values of Late Jurassic dinosaurs.
OPCR values for each genus and theropod morphotype.
About this article
Cite this article
Melstrom, K.M., Chiappe, L.M. & Smith, N.D. Exceptionally simple, rapidly replaced teeth in sauropod dinosaurs demonstrate a novel evolutionary strategy for herbivory in Late Jurassic ecosystems. BMC Ecol Evo 21, 202 (2021). https://0-doi-org.brum.beds.ac.uk/10.1186/s12862-021-01932-4