Genetic architecture of mandible shape in mice: effects of quantitative trait loci analyzed by geometric morphometrics

This study introduces a new multivariate approach for analyzing the effects of quantitative trait loci (QTL) on shape and demonstrates this method for the mouse mandible. We quantified size and shape with the methods of geometric morphometrics, based on Procrustes superimposition of five morphological landmarks recorded on each mandible. Interval mapping for F(2) mice originating from an intercross of the LG/J and SM/J inbred strains revealed 12 QTL for size, 25 QTL for shape, and 5 QTL for left-right asymmetry. Multivariate ordination of QTL effects by principal component analysis identified two recurrent features of shape variation, which involved the positions of the coronoid and angular processes relative to each other and to the rest of the mandible. These patterns are reminiscent of the knockout phenotypes of a number of genes involved in mandible development, although only a few of these are possible candidates for QTL in our study. The variation of shape effects among the QTL showed no evidence of clustering into distinct groups, as would be expected from theories of morphological integration. Further, for most QTL, additive and dominance effects on shape were markedly different, implying overdominance for specific features of shape. We conclude that geometric morphometrics offers a promising new approach to address problems at the interface of evolutionary and developmental genetics.     

2,3,7,8-Tetrachlorodibenzo-p-dioxin affects fluctuating asymmetry of molar shape in mice, and an epistatic interaction of two genes for molar size

Fluctuating asymmetry (FA), random variation between left and right sides in a bilaterally symmetrical character, is a commonly used measure of developmental instability that is expected to increase with increasing environmental stress. One potential stressor is 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), a powerful toxicant known to disturb tooth development. In this study, mice in the F(2) generation produced from an intercross between two inbred strains (C57BL/6J and AKR/J) were exposed in utero to TCDD. We hypothesized that TCDD would increase FA in the molars of exposed mice over that of the control mice. In addition, we hypothesized that we would discover genes for molar size, shape or asymmetry whose expression would be affected by TCDD. We detected a very small, but significant, increase in FA of molar shape (but not size) in the TCDD-exposed mice compared to the control mice, although molar size and shape did not differ between these groups. Although we did not uncover any genes that acted differently in the TCDD exposed and control groups, we did identify two genes whose dominance by additive epistatic effect on molar size was affected by TCDD. We concluded that although TCDD may be affecting the expression of some genes governing the development of molars in our population of mice, FA of molar size and shape is not a particularly sensitive indicator of this effect.     

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

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

An epistatic genetic basis for fluctuating asymmetry of tooth size and shape in mice

Although there typically is little additive genetic variation for fluctuating asymmetry (FA), or variation in nondirectional differences between left and right sides of bilateral characters, several investigators have hypothesized that FA may have an epistatic genetic basis. We tested this hypothesis by conducting a whole genome scan of FA of size and shape of the mandibular molars in house mice from an F2 intercross population generated from crossing the Large (LG/J) and Small (SM/J) inbred strains. Although no individual genes (QTLs=quantitative trait loci) on any of the 19 autosomes significantly affected FA for centroid size, and only two affected shape FA, a number of pairwise combinations of QTLs exhibited significant epistasis for FA in both molar size and shape. The QTLs involved in these interactions differed for FA in molar size versus FA in molar shape, but their epistatic contributions to the total variance was nearly the same (about 20%) for FA in both molar characters. It was noted that the genetic architecture of FA in the molar characters, consisting of little or no additive genetic variance but an abundance of epistatic genetic variance, is consistent with that of other typical fitness components such as litter size.     

A search for quantitative trait loci exhibiting imprinting effects on mouse mandible size and shape

Genomic imprinting refers to the pattern of monoallelic parent-of-origin-dependent gene expression where one of the two alleles at a locus is expressed and the other silenced. Although some genes in mice are known to be imprinted, the true scope of imprinting and its impact on the genetic architecture of a wide range of morphometric traits is mostly unknown. We therefore searched for quantitative trait loci (QTL) exhibiting imprinting effects on mandible size and shape traits in a large F(3) population of mice originating from an intercross of the LG/J (Large) and SM/J (Small) inbred strains. We discovered a total of 51 QTL affecting mandible size and shape, 6 of which exhibited differences between reciprocal heterozygotes, the usual signature of imprinting effects. However, our analysis showed that only one of these QTL (affecting mandible size) exhibited a pattern consistent with true imprinting effects, whereas reciprocal heterozygote differences in the other five all were due to maternal genetic effects. We concluded that genomic imprinting has a negligible effect on these specific morphometric traits, and that maternal genetic effects may account for many of the previously reported instances of apparent genomic imprinting.     

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

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

Integration and modularity of quantitative trait locus effects on geometric shape in the mouse mandible

The mouse mandible has long served as a model system for complex morphological structures. Here we use new methodology based on geometric morphometrics to test the hypothesis that the mandible consists of two main modules, the alveolar region and the ascending ramus, and that this modularity is reflected in the effects of quantitative trait loci (QTL). The shape of each mandible was analyzed by the positions of 16 morphological landmarks and these data were analyzed using Procrustes analysis. Interval mapping in the F(2) generation from intercrosses of the LG/J and SM/J strains revealed 33 QTL affecting mandible shape. The QTL effects corresponded to a variety of shape changes, but ordination or a parametric bootstrap test of clustering did not reveal any distinct groups of QTL that would affect primarily one module or the other. The correlations of landmark positions between the two modules tended to be lower than the correlations between arbitrary subsets of landmarks, indicating that the modules were relatively independent of each other and confirming the hypothesized location of the boundary between them. While these results are in agreement with the hypothesis of modularity, they also underscore that modularity is a question of the relative degrees to which QTL contribute to different traits, rather than a question of discrete sets of QTL contributing to discrete sets of traits.     

Lower molar shape and size in prosimian and platyrrhine primates

The goal of this research is to evaluate the relative strength of the influences of diet, size, and phylogenetic signal on dental geometric shape. Accurate comprehension of these factors and their interaction is important for reconstructing diet and deriving characters for a cladistic analysis in fossil primates. Geometric morphometric analysis is used to identify axes of shape variation in the lower second molars of (a) prosimian primates and (b) platyrrhines. Landmarks were placed on µCT‐generated surface renderings. Landmark configurations were aligned using generalized Procrustes analysis. Principal components analysis and phylogenetic principal components analysis (pPCA) were performed on species average landmark co‐ordinates. pPCs were examined with phylogenetic generalized least squares analysis for association with size and with diet. PCs from both phylogenetic and non‐phylogenetic analyses were sufficient to separate species by broad dietary categories, including insectivores and folivores. In neither analysis was pPC1 correlated with tooth size, but some other pPCs were significantly correlated with size. The pattern of association between pPCs and size altered when centroid size and dietary variables were combined in the model; effects of diet factors typically exceeded effects of size. These results indicate a dominant phylogenetic and dietary signal in molar shape but also show some shape change correlated with size in the absence of obvious dietary associations. Geometric morphometric analysis appears to be useful for tracking functional traits in molars, particularly in tracking differences between folivorous and insectivorous species.     

Discrimination of extant Pan species and subspecies using the enamel–dentine junction morphology of lower molars

Previous research has demonstrated that species and subspecies of extant chimpanzees and bonobos can be distinguished on the basis of the shape of their molar crowns. Thus, there is potential for fossil taxa, particularly fossil hominins, to be distinguished at similar taxonomic levels using molar crown morphology. Unfortunately, due to occlusal attrition, the original crown morphology is often absent in fossil teeth, and this has limited the amount of shape information used to discriminate hominin molars. The enamel–dentine junction (EDJ) of molar teeth preserves considerable shape information, particularly in regard to the original shape of the crown, and remains present through the early stages of attrition. In this study, we investigate whether the shape of the EDJ of lower first and second molars can distinguish species and subspecies of extant Pan. Micro‐computed tomography was employed to non‐destructively image the EDJ, and geometric morphometric analytical methods were used to compare EDJ shape among samples of Pan paniscus (N = 17), Pan troglodytes troglodytes (N = 13), and Pan troglodytes verus (N = 18). Discriminant analysis indicates that EDJ morphology distinguishes among extant Pan species and subspecies with a high degree of reliability. The morphological differences in EDJ shape among the taxa are subtle and relate to the relative height and position of the dentine horns, the height of the dentine crown, and the shape of the crown base, but their existence supports the inclusion of EDJ shape (particularly those aspects of shape in the vertical dimension) in the systematic analysis of fossil hominin lower molars. Am J Phys Anthropol, 2009. © 2009 Wiley‐Liss, Inc.     

Genetic analysis of tongue size and taste papillae number and size in recombinant inbred strains of mice

Quantitative trait loci (QTLs) analysis has been used to examine natural variation of phenotypes in the mouse somatosensory cortex, hippocampus, cerebellum, and amygdala. QTL analysis has also been utilized to map and identify genes underlying anatomical features such as muscle, organ, and body weights. However, this methodology has not been previously applied to identification of anatomical structures related to gustatory phenotypes. In this study, we used QTL analysis to map and characterize genes underlying tongue size, papillae number, and papillae area. In a set of 43 BXD recombinant inbred (RI) mice (n = 111) and 2 parental strains (C57BL/6J and DBA/2J; n = 7), we measured tongue length, width, and weight. In a subset of 23 BXD RI mice and the parental mice, we measured filiform and fungiform papillae number and fungiform papillae area. Using QTL linkage analysis (through WebQTL), we detected 2 significant and noninteracting QTLs influencing tongue length on chromosomes 5 and 7. We also found a significant QTL on chromosome 19 underlying fungiform papillae area and a suggestive QTL on chromosome 2 linked to fungiform papillae number. From these QTLs, we identified a number of candidate genes within the QTL intervals that include SRY-box containing gene, nebulin-related anchoring protein, and actin-binding LIM protein 1. This study is an important first step in identifying genetic factors underlying tongue size, papillae size, and papillae number using QTL analysis.     

Sensitivity to the locomotor-stimulant effects of ethanol and allopregnanolone: a quantitative trait locus study of common genetic influence

Previous studies have suggested that common genetic mechanisms influence sensitivity to the locomotor-stimulant effects of ethanol and allopregnanolone. We conducted two quantitative trait locus (QTL) studies to identify chromosomal regions that harbor genes that influence locomotor response to ethanol (2 g/kg) and allopregnanolone (17 mg/kg) using F2 crosses between C57BL/6J and DBA/2J mice. Because our previous data from the BXD recombinant inbred strains had indicated that chromosome 2 contained QTL for sensitivity to the locomotor-stimulant effects of both ethanol and allopregnanolone, we also tested reciprocal chromosome 2 congenic strains for sensitivity to the locomotor-stimulant effects of both drugs. The F2 analysis for ethanol sensitivity identified significant QTL on chromosomes 1 and 2 and suggestive QTL on chromosomes 5 and 9. The analysis of the allopregnanolone F2 study identified suggestive QTL on chromosomes 3, 5 and 12. Suggestive evidence for a female-specific QTL on chromosome 2 was also found. The studies of congenic mouse strains indicated that both the congenic strains captured one or more QTL for sensitivity to the locomotor-stimulant effects of both ethanol (2 g/kg) and allopregnanolone (17 mg/kg). When Fisher's method was used to combine the P values for the RI, F2 and congenic studies of the chromosome 2 QTL, cumulative probability scores of 9.6 x 10(-15) for ethanol and 7.7 x 10(-7) for allopregnanolone were obtained. These results confirm the presence of QTL for ethanol and allopregnanolone sensitivity in a common region of chromosome 2 and suggest possible pleiotropic genetic influence on sensitivity to these drugs.     

Molar development in common chimpanzees (Pan troglodytes)

Numerous studies have reported on enamel and dentine development in hominoid molars, although little is known about intraspecific incremental feature variation. Furthermore, a recent histological study suggested that there is little or no time between age at chimpanzee crown completion and age at molar eruption, which is unlikely given that root growth is necessary for tooth eruption. The study presented here redefines growth standards for chimpanzee molar teeth and examines variation in incremental features. The periodicity of Retzius lines in a relatively large sample was found to be 6 or 7 days. The number of Retzius lines and cuspal enamel thickness both vary within a cusp type, among cusps, and among molars, resulting in marked variation in formation time. Daily secretion rate is consistent within analogous cuspal zones (inner, middle, and outer enamel) within and among cusp types and among molar types. Significantly increasing trends are found from inner to outer cuspal enamel (3 to 5 microns/day). Cuspal initiation and completion sequences also vary, although sequences for mandibular molar cusps are more consistent. Cusp-specific formation time ranges from approximately 2 to 3 years, increasing from M1 to M2, and often decreasing from M2 to M3. These times are intermediate between radiographic studies and a previous histological study, although both formation time within cusps and overlap between molars vary considerably. Cusp-specific (coronal) extension rates range from approximately 4 to 9 microns/day, and root extension rates in the first 5 mm of roots range from 3 to 9 microns/day. These rates are greater in M1 than in M2 or M3, and they are greater in mandibular molars than in respective maxillary molars. This significant enlargement of comparative data on nonhuman primate incremental development demonstrates that developmental variation among cusp and molar types should be considered during interpretations and comparisons of small samples of fossil hominins and hominoids.     

Enamel-dentine junction (EDJ) morphology distinguishes the lower molars of Australopithecus africanus and Paranthropus robustus

Tooth crown morphology plays a central role in hominin systematics, but the removal of the original outer enamel surface by dental attrition often eliminates from consideration the type of detailed crown morphology that has been shown to discriminate among hominin taxa. This reduces the size of samples available for study. The enamel-dentine junction (EDJ) is the developmental precursor and primary contributor to the morphology of the unworn outer enamel surface, and its morphology is only affected after considerable attrition. In this paper, we explore whether the form of the EDJ can be used to distinguish between the mandibular molars of two southern African fossil hominins: Paranthropus (or Australopithecus) robustus and Australopithecus africanus. After micro-computed tomographic scanning the molar sample, we made high-resolution images of the EDJ and used geometric morphometrics to compare EDJ shape differences between species, in addition to documenting metameric variation along the molar row within each species. Landmarks were collected along the marginal ridge that runs between adjacent dentine horns and around the circumference of the cervix. Our results suggest that the morphology of the EDJ can distinguish lower molars of these southern African hominins, and it can discriminate first, second, and third molars within each taxon. These results confirm previous findings that the EDJ preserves taxonomically valuable shape information in worn teeth. Mean differences in EDJ shape, in particular dentine horn height, crown height, and cervix shape, are more marked between adjacent molars within each taxon than for the same molar between the two taxa.     

Molar tooth development in caspase-3 deficient mice

Tooth morphogenesis is accompanied by apoptotic events which show restricted temporospatial patterns suggesting multiple roles in odontogenesis. Dental apoptosis seems to be caspase dependent and caspase-3 has been shown to be activated during dental apoptosis.Caspase-3 mutant mice on different genetic backgrounds were used to investigate alterations in dental apoptosis and molar tooth morphogenesis. Mouse embryos at E15.5 were analyzed to reveal any changes in enamel knots, which are transient structures eliminated by apoptosis. In caspase-3(-/-) mice on the B57BL/6 background, disorganization of the epithelium was found in the original primary enamel knot area and confirmed by altered expression of Shh. Despite this early defect in molar tooth development, these mutants showed correct formation of secondary enamel knots as indicated by Fgf-4 expression. Analyses of adult molar teeth did not reveal any major alterations in tooth shape, enamel structure or pattern when compared to heterozygote littermates. In caspase-3(-/-) mice on the 129X1/SvJ background, no defects in tooth development were found except the position of the upper molars which developed more posteriorly in the oral cavity. This is likely, however, to be a secondary defect caused by a physical squashing of the face by the malformed brain. The results suggest that although caspase-3 becomes activated and may be essential for dental apoptosis, it does not seem fundamental for formation of normal mineralised molar teeth.     

GENE EFFECTS ON A QUANTITATIVE TRAIT: TWO-LOCUS EPISTATIC EFFECTS MEASURED AT MICROSATELLITE MARKERS AND AT ESTIMATED QTL

Most evolutionarily and agriculturally important traits are affected by many genes (quantitative trait loci, or QTL) of relatively small effect. Usually the genetics of these traits are examined by indirect statistical analysis of the covariance among relatives, rather than by direct analyses. We use new analytical and molecular techniques to examine nonadditive interactions of microsatellite markers and estimated QTL that influence adult body weight in mice. Offspring of a cross between a large inbred mouse strain (LG/J) and a small inbred strain (SM/J) were intercrossed to form a segregating F₂ generation. Using 76 microsatellite markers and 19 estimated QTL, we estimate gene-level epistasis and population-level epistasis for body weight at 10 weeks for 534 F₂ mice. Significant epistasis was found for large numbers of the two locus comparisons using both markers and previously detected QTL. There are many genes segregating for adult body weight in this cross and many of these genes appear to interact epistatically. The discovery of potentially extensive epistasis has important implications for evolutionary models.     

Relative warp analysis of cranial and upper molar shape in rock miceApodemus mystacinus sensu lato

We studied phenotypic relationships among 13 samples of two rock mice species:Apodemus mystacinus (Danford and Alston, 1877) from Anatolia (n = 38) andA. epimelas (Nehring, 1902) from the Balkans (n = 71). Cartesian coordinates of landmarks were collected on the skull and on the occlusal projection of the upper molars (18 landmarks). Centroid size (a measure of overall size) suggested that molars vary independently of overall skull size in both species. Discriminant function analysis on relative warp scores classified >80% of specimens into the correct species, with the best results obtained for the ventral aspect of the skull and for molars. Projection of the 1st discriminant function scores against centroid size provided good separation between the two species. Analysis of vector displacements associated with extremes of variation suggested considerable phenetic differences on the ventral side of the skull and in the molar shape of the two species. The great majority of shifts in landmarks were in a longitudinal direction and the rearrangements of molar cusps were more complex than was the case with the cranium. A bivariate plot of the posterior hard palate length against the incisive foramen length separatedA. mystacinus andA. epimelas well.     

Morphology of Liberian Negro deciduous teeth I. Odontometry

An odontometric study of the deciduous teeth of 21 Liberian Negro children is reported. These data were obtained from dental casts alone. Statistical comparisons were made of the mesiodistal diameter alone with two American White populations and the three Japanese populations. Maxillary teeth, with the exception of the cuspids were larger in the Negro sample, while only the Negro mandibular molars were significantly wider. This situation resulted in an alteration of Negro maxillary tooth size order. Comparisons were made with other populations reported in the literature which led us to stress the fact that, as yet, no satisfactory sample of deciduous tooth measurements are available for any population. A discussion of tooth form differences stresses the point that teeth differ in both attributes of form, i.e., in shape as well as in size between the several major racial groups. The processes by which these shape differences are brought about were discussed and it was suggested that differing relative growth rates may well be one cause. An allometric analysis of fetal molar tooth growth showed, again, that a sharp interphase occurs in all teeth, with a subsequent decrease in specific growth rate. This interphase corresponds to the initiation of cuspal calcification.     

A quantitative genetic analysis of male sexual traits distinguishing the sibling species Drosophila simulans and D. sechellia

A quantitative trait locus (QTL) genetic analysis of morphological and reproductive traits distinguishing the sibling species Drosophila simulans and D. sechellia was carried out in a backcross design, using 38 markers with an average spacing of 8.4 cM. The direction of QTL effects for the size of the posterior lobe was consistent across the identified QTL, indicating directional selection for this trait. Directional selection also appears to have acted on testis length, indicating that sexual selection may have influenced many reproductive traits, although other forms of directional selection cannot be ruled out. Sex comb tooth number exhibited high levels of variation both within and among isofemale lines and showed no evidence for directional selection and, therefore, may not have been involved in the early speciation process. A database search for genes associated with significant QTL revealed a set of candidate loci for posterior lobe shape and size, sex comb tooth number, testis length, tibia length, and hybrid male fertility. In particular, decapentaplegic (dpp), a gene known to influence the genital arch, was found to be associated with the largest LOD peak for posterior lobe shape and size.