Statistical genetics of molar cusp patterning in pedigreed baboons: implications for primate dental development and evolution

Rules of teeth development align microevolution with macroevolution in extant and extinct primates

Macroevolutionary biologists have classically rejected the notion that higher-level patterns of divergence arise through microevolutionary processes acting within populations. For morphology, this consensus partly derives from the inability of quantitative genetics models to correctly predict the behaviour of evolutionary processes at the scale of millions of years. Developmental studies (evo-devo) have been proposed to reconcile micro- and macroevolution. However, there has been little progress in establishing a formal framework to apply evo-devo models of phenotypic diversification. Here we reframe this issue by asking whether using evo-devo models to quantify biological variation can improve the explanatory power of comparative models, thus helping us bridge the gap between micro- and macroevolution. We test this prediction by evaluating the evolution of primate lower molars in a comprehensive dataset densely sampled across living and extinct taxa. Our results suggest that biologically informed morphospaces alongside quantitative genetics models allow a seamless transition between the micro- and macroscales, whereas biologically uninformed spaces do not. We show that the adaptive landscape for primate teeth is corridor like, with changes in morphology within the corridor being nearly neutral. Overall, our framework provides a basis for integrating evo-devo into the modern synthesis, allowing an operational way to evaluate the ultimate causes of macroevolution.

Microevolutionary variation in molar morphology of Onychomys leucogaster decoupled from genetic structure

In neutral models of quantitative trait evolution, both genetic and phenotypic divergence scale as random walks, producing a correlation between the two measures. However, complexity in the genotype-phenotype map may alter the correlation between genotypic and phenotypic divergence, even when both are evolving neutrally or nearly so. Understanding this correlation between phenotypic and genetic variation is critical for accurately interpreting the fossil record. This study compares the geographic structure and scaling of morphological variation of the shape of the first lower molar of 77 individuals of the northern grasshopper mouse Onychomys leucogaster to genome-wide SNP variation in the same sample. We found strong genetic structure but weak or absent morphological structure indicating that the scaling of each type of variation is decoupled from one another. Low P_ST values relative to F_ST values are consistent with a lack of morphological divergence in contrast to genetic divergence between groups. This lack of phenotypic structure and the presence of notable within-sample phenotypic variance are consistent with uniform selection or constraints on molar shape across a wide geographic and environmental range. Over time, this kind of decoupling may result in patterns of phenotypic stasis masking underlying genetic patterns.

Genetic Correlation, Pleiotropy, and Molar Morphology in a Longitudinal Sample of Australian Twins and Families

This study aims to expand our understanding of the genetic architecture of crown morphology in the human diphyodont dentition. Here, we present bivariate genetic correlation estimates for deciduous and permanent molar traits and evaluate the patterns of pleiotropy within (e.g., m1–m2) and between (e.g., m2–M1) dentitions. Morphology was observed and scored from dental models representing participants of an Australian twin and family study (deciduous n = 290, permanent n = 339). Data collection followed Arizona State University Dental Anthropology System standards. Genetic correlation estimates were generated using maximum likelihood variance components analysis in SOLAR v.8.1.1. Approximately 23% of deciduous variance components models and 30% of permanent variance components models yielded significant genetic correlation estimates. By comparison, over half (56%) of deciduous–permanent homologues (e.g., m2 hypocone–M1 hypocone) were significantly genetically correlated. It is generally assumed that the deciduous and permanent molars represent members of a meristic molar field emerging from the primary dental lamina. However, stronger genetic integration among m2–M1/M2 homologues than among paired deciduous traits suggests the m2 represents the anterior-most member of a “true” molar field. The results indicate genetic factors act at distinct points throughout development to generate homologous molar form, starting with the m2, which is later replaced by a permanent premolariform crown.

The genetic architecture of anterior tooth morphology in a longitudinal sample of Australian twins and families

OBJECTIVE

This study presents a quantitative genetic analysis of human anterior dental morphology in a longitudinal sample of known genealogy. The primary aim of this work is to generate a suite of genetic correlations within and between deciduous and permanent characters to access patterns of integration across the diphyodont dental complex.

DESIGN

Data were recorded from casted tooth crowns representing participants of a long-term Australian twin and family study (deciduous n = 290, permanent n = 339). Morphological trait expression was observed and scored following Arizona State University Dental Anthropology System standards. Bivariate genetic correlations were estimated using maximum likelihood variance decomposition models in SOLAR v.8.1.1.

RESULTS

Genetic correlation estimates indicate high levels of integration between antimeres but low to moderate levels among traits within a tooth row. Only 9% of deciduous model comparisons were significant, while pleiotropy was indicated for one third of permanent trait pairs. Canine characters stood out as strongly integrated, especially in the deciduous dentition. For homologous characters across dentitions (e.g., deciduous i¹ shoveling and permanent I¹ shoveling), ∼70% of model comparisons yielded significant genetic correlations.

CONCLUSIONS

Patterns of genetic correlation suggest a morphological canine module that spans the primary and secondary dentition. Results also point to the existence of a genetic mechanism conserving morphology across the diphyodont dental complex, such that paired deciduous and permanent traits are more strongly integrated than characters within individual tooth rows/teeth.

Heterogeneous frailty and the expression of linear enamel hypoplasia in a genealogical population

OBJECTIVES

Linear enamel hypoplasia (LEH) is a common skeletal marker of physiological stress (e.g., malnutrition or illness) that is studied within and across populations, without reference to familial risk. We examine LEH prevalence in a population with known genealogical relationships to determine the potential influence of genetic heritability and shared environment.

METHODS

LEH data of 239 individuals from a single population were recorded from the Ohio State University Menegaz-Bock collection dental casts. All individuals were of known age, sex, and genealogy. Narrow-sense heritability estimates were obtained for LEH presence and count data from all unworn, fully erupted teeth (excluding third molars) using SOLAR (v.8.1.1). Age, sex, and age-sex interaction were included as covariates. Models were re-run with a household effect variable.

RESULTS

LEH persists across generations in this study population with moderate, significant heritability estimates for presence in four teeth, and count in four teeth (three teeth were significant for both). When a household effect variable was added, no residual heritability remained for LEH count on any tooth. There was no significant household effect for three of the four teeth that had significant heritability estimates for LEH presence. Age was a significant covariate. Further analyses with birth year data revealed a secular trend toward less LEH.

CONCLUSIONS

This study provides evidence for familial risk of LEH (genetic and environmental) that has consequences for the broad use of this skeletal marker of stress. These results have repercussions for archaeological assemblages, or population health studies, where genetic relatives and household groups might be heavily represented.

Gene networks, occlusal clocks, and functional patches: new understanding of pattern and process in the evolution of the dentition

Our understanding of the evolution of the dentition has been transformed by advances in the developmental biology, genetics, and functional morphology of teeth, as well as the methods available for studying tooth form and function. The hierarchical complexity of dental developmental genetics combined with dynamic effects of cells and tissues during development allow for substantial, rapid, and potentially non-linear evolutionary changes. Studies of selection on tooth function in the wild and evolutionary functional comparisons both suggest that tooth function and adaptation to diets are the most important factors guiding the evolution of teeth, yet selection against random changes that produce malocclusions (selectional drift) may be an equally important factor in groups with tribosphenic dentitions. These advances are critically reviewed here.

Baboons as a model to study genetics and epigenetics of human disease

A major challenge for understanding susceptibility to common human diseases is determining genetic and environmental factors that influence mechanisms underlying variation in disease-related traits. The most common diseases afflicting the US population are complex diseases that develop as a result of defects in multiple genetically controlled systems in response to environmental challenges. Unraveling the etiology of these diseases is exceedingly difficult because of the many genetic and environmental factors involved. Studies of complex disease genetics in humans are challenging because it is not possible to control pedigree structure and often not practical to control environmental conditions over an extended period of time. Furthermore, access to tissues relevant to many diseases from healthy individuals is quite limited. The baboon is a well-established research model for the study of a wide array of common complex diseases, including dyslipidemia, hypertension, obesity, and osteoporosis. It is possible to acquire tissues from healthy, genetically characterized baboons that have been exposed to defined environmental stimuli. In this review, we describe the genetic and physiologic similarity of baboons with humans, the ability and usefulness of controlling environment and breeding, and current genetic and genomic resources. We discuss studies on genetics of heart disease, obesity, diabetes, metabolic syndrome, hypertension, osteoporosis, osteoarthritis, and intrauterine growth restriction using the baboon as a model for human disease. We also summarize new studies and resources under development, providing examples of potential translational studies for targeted interventions and therapies for human disease.

Morphometric analysis of molars in a Middle Pleistocene population shows a mosaic of 'modern' and Neanderthal features

Previous studies of upper first molar (M1) crown shape have shown significant differences between Homo sapiens and Homo neanderthalensis that were already present in the European Middle Pleistocene populations, including the large dental sample from Atapuerca-Sima de los Huesos (SH). Analysis of other M1 features such as the total crown base area, cusp proportions, cusp angles and occlusal polygon have confirmed the differences between both lineages, becoming a useful tool for the taxonomic assignment of isolated teeth from Late Pleistocene sites. However, until now the pattern of expression of these variables has not been known for the SH sample. This fossil sample, the largest collection from the European Middle Pleistocene, is generally interpreted as being from the direct ancestors of Neanderthals, and thus is a reference sample for assessing the origin of the Neanderthal morphologies. Surprisingly, our study reveals that SH M(1) s present a unique mosaic of H. neanderthalensis and H. sapiens features. Regarding the cusp angles and the relative occlusal polygon area, SH matches the H. neanderthalensis pattern. However, regarding the total crown base area and relative cusps size, SH M(1) s are similar to H. sapiens, with a small crown area, a strong hypocone reduction and a protocone enlargement, although the protocone expansion in SH is significantly larger than in any other group studied. The SH dental sample calls into question the uniqueness of some so-called modern traits. Our study also sounds a note of caution on the use of M(1) occlusal morphology for the alpha taxonomy of isolated M(1) s.

Modularity in the mammalian dentition: mice and monkeys share a common dental genetic architecture

The concept of modularity provides a useful tool for exploring the relationship between genotype and phenotype. Here, we use quantitative genetics to identify modularity within the mammalian dentition, connecting the genetics of organogenesis to the genetics of population-level variation for a phenotype well represented in the fossil record. We estimated the correlations between dental traits owing to the shared additive effects of genes (pleiotropy) and compared the pleiotropic relationships among homologous traits in two evolutionary distant taxa-mice and baboons. We find that in both mice and baboons, who shared a common ancestor >65 Ma, incisor size variation is genetically independent of molar size variation. Furthermore, baboon premolars show independent genetic variation from incisors, suggesting that a modular genetic architecture separates incisors from these posterior teeth as well. Such genetic independence between modules provides an explanation for the extensive diversity of incisor size variation seen throughout mammalian evolution-variation uncorrelated with equivalent levels of postcanine tooth size variation. The modularity identified here is supported by the odontogenic homeobox code proposed for the patterning of the rodent dentition. The baboon postcanine pattern of incomplete pleiotropy is also consistent with predictions from the morphogenetic field model.

Genetic integration of molar cusp size variation in baboons

Many studies of primate diversity and evolution rely on dental morphology for insight into diet, behavior, and phylogenetic relationships. Consequently, variation in molar cusp size has increasingly become a phenotype of interest. In 2007 we published a quantitative genetic analysis of mandibular molar cusp size variation in baboons. Those results provided more questions than answers, as the pattern of genetic integration did not fit predictions from odontogenesis. To follow up, we expanded our study to include data from the maxillary molar cusps. Here we report on these later analyses, as well as inter-arch comparisons with the mandibular data. We analyzed variation in two-dimensional maxillary molar cusp size using data collected from a captive pedigreed breeding colony of baboons, Papio hamadryas, housed at the Southwest National Primate Research Center. These analyses show that variation in maxillary molar cusp size is heritable and sexually dimorphic. We also estimated additive genetic correlations between cusps on the same crown, homologous cusps along the tooth row, and maxillary and mandibular cusps. The pattern for maxillary molars yields genetic correlations of one between the paracone-metacone and protocone-hypocone. Bivariate analyses of cuspal homologues on adjacent teeth yield correlations that are high or not significantly different from one. Between dental arcades, the nonoccluding cusps consistently yield high genetic correlations, especially the metaconid-paracone and metaconid-metacone. This pattern of genetic correlation does not immediately accord with the pattern of development and/or calcification, however these results do follow predictions that can be made from the evolutionary history of the tribosphenic molar.

Evolution of M1 crown size and cusp proportions in the genus Homo

Previous research into tooth crown dimensions and cusp proportions has proved to be a useful way to identify taxonomic differences in Pliocene and Pleistocene fossil hominins. The present study has identified changes in both M(1) crown size and cusp proportions within the genus Homo, with M(1) overall crown size reduction apparently occurring in two main stages. The first stage (a reduction of ca. 17%) is associated with the emergence of Homo ergaster and Homo erectus sensu stricto. The second stage (a reduction of ca. 10%) occurs in Homo sapiens, but the reduced modern human M(1) tooth crown size was only attained in Upper Paleolithic times. The absolute sizes of the individual cusps are highly positively correlated with overall crown size and dental reduction produces a reduction in the absolute size of each of the cusps. Most of the individual cusps scale isometrically with crown size, but the paracone shows a negative allometric relationship, indicating that the reduction in paracone size is less than in the other M(1) cusps. Thus, the phylogenetically oldest cusp in the upper molars also seems to be the most stable cusp (at least in the M(1)). The most striking change in M(1) cusp proportions is a change in the relative size of the areas of the paracone and metacone. The combination of a small relative paracone and a large relative metacone generally characterizes specimens attributed to early Homo, and the presence of this character state in Australopithecus and Paranthropus suggests it may represent the primitive condition for the later part of the hominin clade. In contrast, nearly all later Homo taxa, with the exception of Homo antecessor, show the opposite condition (i.e. a relatively large paracone and a relatively small metacone). This change in the relationship between the relative sizes of the paracone and metacone is related to an isometric reduction of the absolute size of the metacone. This metacone reduction occurs in the context of relative stability in the paracone as crown size decreases. Among later Homo taxa, both Homo heidelbergensis and Homo neanderthalensis show a further reduction of the metacone and an enlargement of the hypocone. Fossil and contemporary H. sapiens samples show a trend toward increasing the relative size of the protocone and decreasing the relative size of the hypocone. In Europe, modern human M(1) cusp proportions are essentially reached during the Upper Paleolithic. Although some variation was documented among the fossil taxa, we suggest that the relative size of the M(1) paracone and metacone areas may be useful for differentiating the earliest members of our genus from subsequent Homo species.

Quantitative Genetics, Pleiotropy, and Morphological Integration in the Dentition of Papio hamadryas

Variation in the mammalian dentition is highly informative of adaptations and evolutionary relationships, and consequently has been the focus of considerable research. Much of the current research exploring the genetic underpinnings of dental variation can trace its roots to Olson and Miller's 1958 book Morphological Integration. These authors explored patterns of correlation in the post-canine dentitions of the owl monkey and Hyopsodus, an extinct condylarth from the Eocene. Their results were difficult to interpret, as was even noted by the authors, due to a lack of genetic information through which to view the patterns of correlation. Following in the spirit of Olson and Miller's research, we present a quantitative genetic analysis of dental variation in a pedigreed population of baboons. We identify patterns of genetic correlations that provide insight to the genetic architecture of the baboon dentition. This genetic architecture indicates the presence of at least three modules: an incisor module that is genetically independent of the post-canine dentition, and a premolar module that demonstrates incomplete pleiotropy with the molar module. We then compare this matrix of genetic correlations to matrices of phenotypic correlations between the same measurements made on museum specimens of another baboon subspecies and the Southeast Asian colobine Presbytis. We observe moderate significant correlations between the matrices from these three primate taxa. From these observations we infer similarity in modularity and hypothesize a common pattern of genetic integration across the dental arcade in the Cercopithecoidea.

Genetic correlations between mandibular molar cusp areas in baboons

Primate evolutionary studies rely significantly on dental variation given the large role that teeth play in how an organism interacts with its environment (animal and plant) and conspecifics. Variation in cusp size has been shown to vary among primate taxa, although most studies to date focused on extant and extinct hominoids. Here we test the assumed hypothesis that a significant proportion of this variation in baboons is due to the additive effects of genes. We perform quantitative genetic analyses on variation in two-dimensional (2-D) mandibular molar cusp size in a captive pedigreed breeding population of baboons (Papio hamadryas) from the Southwest National Primate Research Center. These analyses show that variation in cusp size is heritable and sexually dimorphic. Additionally, we tested for genetic correlations between cusps on the same crown, between morphological homologues along the tooth row, and between cusp area and crown buccolingual width. We find that four of the six cusp pairs on the first molar have a genetic correlation of one, save for the metaconid-hypoconid and entoconid-hypoconid, which are not statistically different from zero. The second and third molars have lower genetic correlations, although the metaconid-hypoconid correlation is similarly estimated at zero and the entoconid-protoconid correlation is estimated to be one. This cross pattern of genetic and no genetic correlation does not immediately accord with the known pattern of development and/or calcification. We propose two explanative hypotheses.

A multivariate comparison of dental variation in wild and captive populations of baboons (Papio hamadryas)

Phenotypic variation is critical to many aspects of biological research. Use of a captive population to address questions concerning the genetics and evolution of dental variation raises the question of how the pattern of phenotypic variation under study compares with that in a wild population of the same species. Differences in the pattern of variation within wild and captive populations may indicate different genetic and non-genetic factors, and also may have implications for how well the captive group can serve as a model for its wild type relatives. We compared dental size measures from two Papio hamadryas populations, one captive and one wild. Lengths and widths of maxillary and mandibular second molars (M2s) were collected from 630 baboons from a captive pedigreed breeding colony housed at the Southwest National Primate Research Center in San Antonio, Texas, and 125 baboons from a wild population culled from a sisal plantation in Kibweze, Kenya. Although the two populations consistently differed with respect to lengths and widths of the M2s, principal components analyses show that the basic pattern to variation in these molar crown traits is remarkably similar in both populations; and linear functions based on these measures cannot reliably discriminate between the two groups. This similarity in the pattern of variation among these dental crown measures in these two groups suggests that analyses to dissect their genetic architecture in captive populations is likely to be highly relevant to dental variation in wild baboons as well.

Associations between Carabelli trait and cusp areas in human permanent maxillary first molars

Few dental anthropological studies have investigated the associations between tooth crown size and crown traits in humans using quantitative methods. We tested several hypotheses about overall crown size, individual cusp areas, and expression of Carabelli cusps in human permanent first molars by obtaining data from standardized occlusal photographs of 308 Australians of European descent (171 males and 137 females). Specifically, we aimed to calculate the areas of the four main molar cusps, and also Carabelli cusp, and to compare the relative variability of cusp areas in relation to timing of development. We also aimed to compare cusp areas between males and females and to describe how Carabelli cusp interacted with other molar cusps. Measurements included maximum crown diameters (mesiodistal and buccolingual crown diameters), the areas of the four main cusps, and the area of Carabelli cusp. The pattern of relative variability in absolute areas of molar cusps corresponded with their order of formation, the first-forming paracone displaying the least variation, and the last-forming Carabelli cusp showing the greatest. Overall crown size and areas of individual cusps all showed sexual dimorphism, with values in males exceeding those in females. Sexual dimorphism was smallest for paracone area and greatest for Carabelli cusp area. Overall crown size and cusp areas were larger in individuals displaying a Carabelli cusp, especially the hypocone area. Although the combined area of the protocone and a Carabelli cusp was greater in cuspal forms than noncuspal forms, protocone area alone was significantly smaller in the former. Our findings lead us to propose that, in individuals with the genotype for Carabelli trait expression, larger molar crowns are more likely to display Carabelli cusps, whereas molars with smaller crowns are more likely to display reduced forms of expression of the trait. We suggest that the pattern of folding of the internal enamel epithelium in developing molar crowns, particularly in the protocone region, can be modified by a developing Carabelli cusp.

New approach to quantifying developmental variation in the dentition using serial microtomographic imaging

Dental morphogenesis and cellular differentiation are expressed in the fully formed tooth by the topography of the dentin-enamel junction and outer enamel surface. These boundaries can be differentiated using a variety of imaging systems. In this study, we used serial microCT imaging to provide accurate 3D reconstructions of developing lower human second deciduous molars. These were used to quantify the volume of enamel and dentin of individual cusps in relation to basal crown height. As growth and differentiation proceed apically, the spatial orientation of cusp tips and their bases were used to estimate their order of initiation and coalescence. We found that the order of coalescence differed from the order of initiation. We also found that dentin cusp height and volume as well as rate and quantity of enamel apposition varied along mesio-distal and bucco-lingual axes, and were independent of order of initiation and duration of growth. These results demonstrate that the potential for variation in crown size and form is maintained throughout development. We propose that the microCT model developed in this study constitutes a new approach for the investigation of developmental variation and its contribution to phylogenetic variation expressed in crown form and size.

Development and phenotypic correlations: the evolution of tooth shape in Sorex araneus

Do morphogenetic processes cause common patterns of phenotypic covariation, and do those patterns evolve over microevolutionary timescales? Evolution of molar shape variance-covariance (P) matrixes was studied in five populations of the common shrew, Sorex araneus. P matrix evolution was assessed using matrix correlation, matrix disparity, and common principal component analysis (CPCA). Significant changes in covariance structure were found among the populations, but the differences were small. A computer model was used to estimate the theoretical covariance introduced into the phenotype by developmental interactions. Molar developmental processes explained some of the covariance in the shrew samples, especially as measured by matrix correlation, but the proportion was relatively small. Developmental principal components (PCs) were only infrequently associable with common principal components. The results suggest that molar shape P matrixes can evolve quickly in a manner only loosely constrained by development, and that their shared covariance is probably dominated by factors more proximate than development. Rarefaction showed that sample size severely affected P comparisons when n < 15 for matrix correlation and disparity, and when n < 30 for CPCA. Among CPCA evaluation criteria, Akaike Information Criterion performed better than jump-up at n < 30, but worse at n > 30.