Craniofacial variability and morphological integration in mice susceptible to cleft lip and palate

Using data-driven phenotyping to investigate the impact of sex on 3D human facial surface morphology

The effects of sex on human facial morphology have been widely documented. Because sexual dimorphism is relevant to a variety of scientific and applied disciplines, it is imperative to have a complete and accurate account of how and where male and female faces differ. We apply a comprehensive facial phenotyping strategy to a large set of existing 3D facial surface images. We investigate facial sexual dimorphism in terms of size, shape, and shape variance. We also assess the ability to correctly assign sex based on shape, both for the whole face and for subregions. We applied a predefined data-driven segmentation to partition the 3D facial surfaces of 2446 adults into 63 hierarchically linked regions, ranging from global (whole face) to highly localized subparts. Each facial region was then analyzed with spatially dense geometric morphometrics. To describe the major modes of shape variation, principal components analysis was applied to the Procrustes aligned 3D points comprising each of the 63 facial regions. Both nonparametric and permutation-based statistics were then used to quantify the facial size and shape differences and visualizations were generated. Males were significantly larger than females for all 63 facial regions. Statistically significant sex differences in shape were also seen in all regions and the effects tended to be more pronounced for the upper lip and forehead, with more subtle changes emerging as the facial regions became more granular. Males also showed greater levels of shape variance, with the largest effect observed for the central forehead. Classification accuracy was highest for the full face (97%), while most facial regions showed an accuracy of 75% or greater. In summary, sex differences in both size and shape were present across every part of the face. By breaking the face into subparts, some shape differences emerged that were not apparent when analyzing the face as a whole. The increase in facial shape variance suggests possible evolutionary origins and may offer insights for understanding congenital facial malformations. Our classification results indicate that a high degree of accuracy is possible with only parts of the face, which may have implications for biometrics applications.

What’s Shape Got to Do With It? Examining the Relationship Between Facial Shape and Orofacial Clefting

Nonsyndromic orofacial clefts belong to a class of congenital malformations characterized by a complex and multifactorial etiology. During early facial development, multiple factors can disrupt fusion leading to a cleft; this includes the shape of the embryonic face. The face shape hypothesis (FSH) of orofacial clefting emerged in the 1960s, influenced by morphological differences observed within affected families, comparative studies of mouse models, and advances in modeling genetic liability for complex traits in populations. For the past five decades, studies have documented changes in the shape or spatial arrangement of facial prominences in embryonic mice and altered post-natal facial shape in individuals at elevated risk for orofacial clefting due to their family history. Moreover, recent studies showing how genes that impact facial shape in humans and mice are providing clues about the genetic basis of orofacial clefting. In this review, I discuss the origins of the FSH, provide an overview of the supporting evidence, and discuss ways in which the FSH can inform our understanding of orofacial clefting.

Canalization and developmental stability of the yellow-necked mouse (Apodemus flavicollis) mandible and cranium related to age and nematode parasitism

BACKGROUND

Mammalian mandible and cranium are well-established model systems for studying canalization and developmental stability (DS) as two elements of developmental homeostasis. Nematode infections are usually acquired in early life and increase in intensity with age, while canalization and DS of rodent skulls could vary through late postnatal ontogeny. We aimed to estimate magnitudes and describe patterns of mandibular and cranial canalization and DS related to age and parasite intensity (diversity) in adult yellow-necked mice (Apodemus flavicollis).

RESULTS

We found the absence of age-related changes in the levels of canalization for mandibular and cranial size and DS for mandibular size. However, individual measures of mandibular and cranial shape variance increased, while individual measures of mandibular shape fluctuating asymmetry (FA) decreased with age. We detected mandibular and cranial shape changes during postnatal ontogeny, but revealed no age-related dynamics of their covariance structure among and within individuals. Categories regarding parasitism differed in the level of canalization for cranial size and the level of DS for cranial shape. We observed differences in age-related dynamics of the level of canalization between non-parasitized and parasitized animals, as well as between yellow-necked mice parasitized by different number of nematode species. Likewise, individual measures of mandibular and cranial shape FA decreased with age for the mandible in the less parasitized category and increased for the cranium in the most parasitized category.

CONCLUSIONS

Our age-related results partly agree with previous findings. However, no rodent study so far has explored age-related changes in the magnitude of FA for mandibular size or mandibular and cranial FA covariance structure. This is the first study dealing with the nematode parasitism-related canalization and DS in rodents. We showed that nematode parasitism does not affect mandibular and cranial shape variation and covariance structure among and within individuals. However, parasite intensity (diversity) is related to ontogenetic dynamics of the levels of canalization and DS. Overall, additional studies on animals from natural populations are required before drawing some general conclusions.

Phenotypic integration in feliform carnivores: Covariation patterns and disparity in hypercarnivores versus generalists

The skeleton is a complex arrangement of anatomical structures that covary to various degrees depending on both intrinsic and extrinsic factors. Among the Feliformia, many species are characterized by predator lifestyles providing a unique opportunity to investigate the impact of highly specialized hypercarnivorous diet on phenotypic integration and shape diversity. To do so, we compared the shape of the skull, mandible, humerus, and femur of species in relation to their feeding strategies (hypercarnivorous vs. generalist species) and prey preference (predators of small vs. large prey) using three-dimensional geometric morphometric techniques. Our results highlight different degrees of morphological integration in the Feliformia depending on the functional implication of the anatomical structure, with an overall higher covariation of structures in hypercarnivorous species. The skull and the forelimb are not integrated in generalist species, whereas they are integrated in hypercarnivores. These results can potentially be explained by the different feeding strategies of these species. Contrary to our expectations, hypercarnivores display a higher disparity for the skull than generalist species. This is probably due to the fact that a specialization toward high-meat diet could be achieved through various phenotypes. Finally, humeri and femora display shape variations depending on relative prey size preference. Large species feeding on large prey tend to have robust long bones due to higher biomechanical constraints.

Morphological association between the muscles and bones in the craniofacial region

The strains of inbred laboratory mice are isogenic and homogeneous for over 98.6% of their genomes. However, geometric morphometric studies have demonstrated clear differences among the skull shapes of various mice strains. The question now arises: why are skull shapes different among the mice strains? Epigenetic processes, such as morphological interaction between the muscles and bones, may cause differences in the skull shapes among various mice strains. To test these predictions, the objective of this study is to examine the morphological association between a specific part of the skull and its adjacent muscle. We examined C57BL6J, BALB/cA, and ICR mice on embryonic days (E) 12.5 and 16.5 as well as on postnatal days (P) 0, 10, and 90. As a result, we found morphological differences between C57BL6J and BALB/cA mice with respect to the inferior spine of the hypophyseal cartilage or basisphenoid (SP) and the tensor veli palatini muscle (TVP) during the prenatal and postnatal periods. There was a morphological correlation between the SP and the TVP in the C57BL6J, BALB/cA, and ICR mice during E15 and P0. However, there were not correlation between the TVP and the SP during P10. After discectomy, bone deformation was associated with a change in the shape of the adjacent muscle. Therefore, epigenetic modifications linked to the interaction between the muscles and bones might occur easily during the prenatal period, and inflammation seems to allow epigenetic modifications between the two to occur.

Methods for Studying Morphological Integration and Modularity

Morphological integration and modularity are closely related concepts about how different traits of an organism are correlated. Integration is the overall pattern of intercorrelation; modularity is the partitioning of integration into evolutionarily or developmentally independent blocks of traits. Modularity and integration are usually studied using quantitative phenotypic data, which can be obtained either from extant or fossil organisms. Many methods are now available to study integration and modularity, all of which involve the analysis of patterns found in trait correlation or covariance matrices. We review matrix correlation, random skewers, fluctuating asymmetry, cluster analysis, Euclidean distance matrix analysis (EDMA), graphical modelling, two-block partial least squares, RV coefficients, and theoretical matrix modelling and discuss their similarities and differences. We also review different coefficients that are used to measure correlations. We apply all the methods to cranial landmark data from and ontogenetic series of Japanese macaques,Macaca fuscatato illustrate the methods and their individual strengths and weaknesses. We conclude that the exploratory approaches (cluster analyses of various sorts) were less informative and less consistent with one another than were the results of model testing or comparative approaches. Nevertheless, we found that competing models of modularity and integration are often similar enough that they are not statistically distinguishable; we expect, therefore, that several models will often be significantly correlated with observed data.

Development Shapes a Consistent Inbreeding Effect in Mouse Crania of Different Line Crosses

Development translates genetic variation into a multivariate pattern of phenotypic variation, distributing it among traits in a nonuniform manner. As developmental processes are largely shared within species, this suggests that heritable phenotypic variation will be patterned similarly, in spite of the different segregating alleles. To investigate developmental effect on the variational pattern in the shape of the mouse skull across genetically differentiated lines, we employed the full set of reciprocal crosses (a.k.a. diallel) between eight inbred mouse strains of the Collaborative Cross Project. We used geometric morphometrics and multivariate analysis to capture cranial size and shape changes in 8 parentals and their 54 F1 crosses. The high heterozygosity generated in the F1 crosses allowed us to compare the multivariate deviations of the F1 phenotypes from the expected midparental phenotypes in different haplotype combinations. In contrast to body weight, we found a high degree of nonadditive deviation in craniofacial shape. Whereas the phenotypic and genetic divergence of parental strains manifested in high dimensionality of additive effects, the nonadditive deviations exhibited lesser dimensionality and in particular a strikingly coherent direction in shape space. We interpret this finding as evidence for a strong structuring effect of a relatively small set of developmental processes on the mapping of genetic to phenotypic variation.

Effects of cranial integration on hominid endocranial shape

Because brains do not fossilize, the internal surface of the braincase (endocast) serves as an important source of information about brain growth, development, and evolution. Recent studies of endocranial morphology and development in great apes, fossil hominins, and modern humans have revealed taxon-specific differences. However, it remains to be investigated to which extent differences in endocranial morphology reflect differences in actual brain morphology and development, and to which extent they reflect different interactions of the brain and its case with the cranial base and face. Here we address this question by analyzing the effects of cranial integration on endocranial morphology. We test the 'spatial packing' and 'facial orientation' hypotheses, which propose that size and orientation of the neurocranium relative to the viscerocranium influence endocranial shape. Results show that a substantial proportion of endocranial shape variation along and across ontogenetic trajectories is due to cranial integration. Specifically, the uniquely globular shape of the human endocast mainly results from the combination of an exceptionally large brain with a comparatively small face. Overall, thus, cranial integration has pervasive effects on endocranial morphology, and only a comparatively small proportion of inter- and intra-taxon variation can directly be associated with variation in brain morphology.

Patterns of morphological integration in the dental arches of individuals with malocclusion

OBJECTIVES

In humans, there is a large range of variation in the form of the maxillary and mandibular dental arches. This variation can manifest as either prognathism or retrognathism in either or both arches, which can cause malocclusion and lead to abnormal masticatory function. This study aims to identify aspects of variation and morphological integration existing in the dental arches of individuals with different types of malocclusion.

METHODS

Coordinate landmark data were collected along the gingival margins of 397 scanned dental casts and then analyzed using geometric morphometric techniques to explore arch form variation and patterns of morphological integration within each malocclusion type.

RESULTS

Significant differences were identified between Class II forms (increased projection of upper arch relative to the lower arch) and Class III forms (lower arch projection beyond the upper arch) in symmetrical shape variation, including anteroposterior arch discrepancies and abnormal anterior arch divergence or convergence. Partial least squares analysis demonstrated that Class III dental arches have higher levels of covariance between upper and lower arches (RV = 0.91) compared to the dental arches of Class II (RV = 0.78) and Class I (RV = 0.73). These high levels of covariance, however, are on the lower end of the overall range of possible masticatory blocks, indicating weaker than expected levels of integration.

CONCLUSIONS

This study provides evidence for patterns of variation in dental arch shape found in individuals with Class II and Class III malocclusions. Moreover, differences in integration found between malocclusion types have ramifications for how such conditions should be studied and treated. Am. J. Hum. Biol. 28:879-889, 2016. © 2016Wiley Periodicals, Inc.

A Multivariate Analysis of Unilateral Cleft Lip and Palate Facial Skeletal Morphology

Unilateral cleft lip and palate (UCLP) occurs when the maxillary and nasal facial prominences fail to fuse correctly during development, resulting in a palatal cleft and clefted soft and hard tissues of the dentoalveolus. The UCLP deformity may compromise an individual's ability to eat, chew, and speak. In this retrospective cross-sectional study, cone beam computed tomography (CBCT) images of 7-17-year-old individuals born with UCLP (n = 24) and age- and sex-matched controls (n = 24) were assessed. Coordinate values of three-dimensional anatomical landmarks (n = 32) were recorded from each CBCT image. Data were evaluated using principal coordinates analysis (PCOORD) and Euclidean distance matrix analysis (EDMA). Approximately 40% of morphometric variation is captured by PCOORD axes 1-3, and the negative and positive ends of each axis are associated with specific patterns of morphological differences. Approximately 36% of facial skeletal measures significantly differ by confidence interval testing (α = 0.10) between samples. Although significant form differences occur across the facial skeleton, strong patterns of morphological differences were localized to the lateral and superioinferior aspects of the nasal aperture, particularly on the clefted side of the face. The UCLP deformity strongly influences facial skeletal morphology of the midface and oronasal facial regions, and to a lesser extent the upper and lower facial skeletons. The pattern of strong morphological differences in the oronasal region combined with differences across the facial complex suggests that craniofacial bones are integrated and covary, despite influences from the congenital cleft.

Morphologic variability of nonsyndromic operated patients affected by cleft lip and palate: a geometric morphometric study

INTRODUCTION

In this study, we compared patterns of morphologic variations of the craniofacial skeleton between patients affected by clefts who were operated on and unaffected subjects, aiming to discuss possible morpho-functional consequences of treatment in craniofacial development.

METHODS

The lateral cephalograms of 76 subjects, comprising patients with operated unilateral cleft lip and palate (OpC) and a group matched for sex and age without cleft, were used. Thirteen landmarks were used as variables in geometric morphometric tests quantifying and describing overall shape variation, differences between group means, allometry, and upper-lower face covariation.

RESULTS

The OpC group showed broader shape variations including noncleft group characteristics, but mainly a retrognathic maxilla, a vertically elongated face, a more open mandibular angle, and a more closed basicranial angle. Group means differed mainly in the maxillomandibular relationships. Allometry differed between groups, with the smallest OpC patients showing the most altered morphology. Upper and lower face covariation was stronger in the OpC group, showing mainly vertical changes in the anterior face.

CONCLUSIONS

Operated patients affected by clefts achieve a broad range of morphologies; the most altered were found in those with skeletal Class III and small size. Furthermore, their strongest upper and lower face shape covariation suggests that a harmonic dental occlusion could be a key factor in achieving "normal" craniofacial morphology.

The morphology of the mouse masticatory musculature

The mouse has been the dominant model organism in studies on the development, genetics and evolution of the mammalian skull and associated soft-tissue for decades. There is the potential to take advantage of this well studied model and the range of mutant, knockin and knockout organisms with diverse craniofacial phenotypes to investigate the functional significance of variation and the role of mechanical forces on the development of the integrated craniofacial skeleton and musculature by using computational mechanical modelling methods (e.g. finite element and multibody dynamic modelling). Currently, there are no detailed published data of the mouse masticatory musculature available. Here, using a combination of micro-dissection and non-invasive segmentation of iodine-enhanced micro-computed tomography, we document the anatomy, architecture and proportions of the mouse masticatory muscles. We report on the superficial masseter (muscle, tendon and pars reflecta), deep masseter, zygomaticomandibularis (anterior, posterior, infraorbital and tendinous parts), temporalis (lateral and medial parts), external and internal pterygoid muscles. Additionally, we report a lateral expansion of the attachment of the temporalis onto the zygomatic arch, which may play a role in stabilising this bone during downwards loading. The data presented in this paper now provide a detailed reference for phenotypic comparison in mouse models and allow the mouse to be used as a model organism in biomechanical and functional modelling and simulation studies of the craniofacial skeleton and particularly the masticatory system.

The Developmental Basis of Quantitative Craniofacial Variation in Humans and Mice

The human skull is a complex and highly integrated structure that has long held the fascination of anthropologists and evolutionary biologists. Recent studies of the genetics of craniofacial variation reveal a very complex and multifactorial picture. These findings contrast with older ideas that posit much simpler developmental bases for variation in cranial morphology such as the growth of the brain or the growth of the chondrocranium relative to the dermatocranium. Such processes have been shown to have major effects on cranial morphology in mice. It is not known, however, whether they are relevant to explaining normal phenotypic variation in humans. To answer this question, we obtained vectors of shape change from mutant mouse models in which the developmental basis for the craniofacial phenotype is known to varying degrees, and compared these to a homologous dataset constructed from human crania obtained from a single population with a known genealogy. Our results show that the shape vectors associated with perturbations to chondrocranial growth, brain growth, and body size in mice do largely correspond to axes of covariation in humans. This finding supports the view that the developmental basis for craniofacial variation funnels down to a relatively small number of key developmental processes that are similar across mice and humans. Understanding these processes and how they influence craniofacial shape provides fundamental insights into the developmental basis for evolutionary change in the human skull as well as the developmental-genetic basis for normal phenotypic variation in craniofacial form.

Orofacial Clefts at Bugando Medical Centre: Associated Factors and Postsurgical Complications

Objective

To determine factors associated with orofacial clefts and postsurgical complications of cleft lip and palate repair surgeries in northwestern Tanzania.

Methods

This was a cohort study involving patients with orofacial clefts. Associated factors (family history of orofacial clefts, maternal use of alcohol and cigarette smoking during pregnancy) were obtained through interviews with accompanying parents. Antenatal cards were used to obtain maternal age at birth and birth weight.

Results

Ninety-four patients with different orofacial clefts were seen. Among them, 46.8% (44/94), 13.8% (13/94), and 39.4% (37/94) had cleft lip, cleft palate, and cleft lip and palate, respectively. About 15% of orofacial cleft cases had a positive family history of orofacial clefts. Among these, 7.4% had an affected relative on the maternal side, 4.3% had an affected relative on the paternal side, and 3.2% had an affected sibling. This difference was statistically significant (chi-square = 27.7, p < .001). Orofacial cleft was significantly associated with order of birth (chi-square = 21.0, p < .001). Postoperative complications observed included palatal fistula and philtrum dehiscence.

Conclusion

Family history of orofacial clefts and order of birth were significantly associated with orofacial clefts in northwestern Tanzania. These factors have been associated with risk of orofacial clefts elsewhere and suggest a hereditary role in the etiology of orofacial cleft. Palatal fistula and philtrum dehiscence were postsurgical complications observed in orofacial clefts patients who had primary surgery past the recommended age. These complications could have resulted from delayed surgery and absence of presurgical procedures.

If the skull fits: magnetic resonance imaging and microcomputed tomography for combined analysis of brain and skull phenotypes in the mouse

The mammalian brain and skull develop concurrently in a coordinated manner, consistently producing a brain and skull that fit tightly together. It is common that abnormalities in one are associated with related abnormalities in the other. However, this is not always the case. A complete characterization of the relationship between brain and skull phenotypes is necessary to understand the mechanisms that cause them to be coordinated or divergent and to provide perspective on the potential diagnostic or prognostic significance of brain and skull phenotypes. We demonstrate the combined use of magnetic resonance imaging and microcomputed tomography for analysis of brain and skull phenotypes in the mouse. Co-registration of brain and skull images allows comparison of the relationship between phenotypes in the brain and those in the skull. We observe a close fit between the brain and skull of two genetic mouse models that both show abnormal brain and skull phenotypes. Application of these three-dimensional image analyses in a broader range of mouse mutants will provide a map of the relationships between brain and skull phenotypes generally and allow characterization of patterns of similarities and differences.

Vitamin B-complex application promotes secondary palate development in a palate organ model of the A/WySnJ mouse

PURPOSE

This study analyzed the direct influence of vitamin B-complex supplements (Polybion N, Merck Pharma GmbH, Germany) in medium on secondary palatal development in palatal organ cultures of A/WySnJ mice. Because of positive clinical experiences with prophylactic vitamin B substitution in mothers of cleft-related families, the direct influence of the vitamin B-complex on palatal tissue was analyzed.

MATERIALS AND METHODS

The inbred A/WySnJ mouse strain shows a highly spontaneous, genetically determined clefting rate of 20% to 44%. One hundred seventy-seven A/WySnJ fetuses were microdissected on gestational day 14.3 before the occurrence of palatal fusion. Palatal organ cultures were prepared and incubated in chemically defined serum-free medium with different concentrations (0.1% and 1.0%) of the vitamin B-complex Polybion N for 72 hours. Palatal development was analyzed microscopically according to the 6-step visual scale that describes the approximation of palatal shelves during development.

RESULTS

At the beginning of the experiment (gestational day 14.3), the palatal development of all specimens used for in vitro organ culture showed a clear approach of the palatal shelves at stage II (2.25±0.78). Seventy-two hours after in vitro cultivation, the palatal shelves of the organ cultures supplemented with the vitamin B-complex showed significant growth (0.1%, P=.00017; 1.0%, P=.00078), whereas the untreated control group remained at initial developmental stage II (P=.291).

CONCLUSIONS

The results of this in vitro study suggest a significant positive influence of vitamin B supplementation on palatal shelf development in organ culture. Further studies will focus on the vitamin B concentration in the amniotic fluid of dams with or without cleft in their offspring.

Missense mutation in the second RNA binding domain reveals a role for Prkra (PACT/RAX) during skull development

Random chemical mutagenesis of the mouse genome can causally connect genes to specific phenotypes. Using this approach, reduced pinna (rep) or microtia, a defect in ear development, was mapped to a small region of mouse chromosome 2. Sequencing of this region established co-segregation of the phenotype (rep) with a mutation in the Prkra gene, which encodes the protein PACT/RAX. Mice homozygous for the mutant Prkra allele had defects not only in ear development but also growth, craniofacial development and ovarian structure. The rep mutation was identified as a missense mutation (Serine 130 to Proline) that did not affect mRNA expression, however the steady state level of RAX protein was significantly lower in the brains of rep mice. The mutant protein, while normal in most biochemical functions, was unable to bind dsRNA. In addition, rep mice displayed altered morphology of the skull that was consistent with a targeted deletion of Prkra showing a contribution of the gene to craniofacial development. These observations identified a specific mutation that reduces steady-state levels of RAX protein and disrupts the dsRNA binding function of the protein, demonstrating the importance of the Prkra gene in various aspects of mouse development.

Evolution of adaptive phenotypic variation patterns by direct selection for evolvability

A basic assumption of the Darwinian theory of evolution is that heritable variation arises randomly. In this context, randomness means that mutations arise irrespective of the current adaptive needs imposed by the environment. It is broadly accepted, however, that phenotypic variation is not uniformly distributed among phenotypic traits, some traits tend to covary, while others vary independently, and again others barely vary at all. Furthermore, it is well established that patterns of trait variation differ among species. Specifically, traits that serve different functions tend to be less correlated, as for instance forelimbs and hind limbs in bats and humans, compared with the limbs of quadrupedal mammals. Recently, a novel class of genetic elements has been identified in mouse gene-mapping studies that modify correlations among quantitative traits. These loci are called relationship loci, or relationship Quantitative Trait Loci (rQTL), and affect trait correlations by changing the expression of the existing genetic variation through gene interaction. Here, we present a population genetic model of how natural selection acts on rQTL. Contrary to the usual neo-Darwinian theory, in this model, new heritable phenotypic variation is produced along the selected dimension in response to directional selection. The results predict that selection on rQTL leads to higher correlations among traits that are simultaneously under directional selection. On the other hand, traits that are not simultaneously under directional selection are predicted to evolve lower correlations. These results and the previously demonstrated existence of rQTL variation, show a mechanism by which natural selection can directly enhance the evolvability of complex organisms along lines of adaptive change.

In vitro investigation of the secondary palate development in two strains of mice

The pathogenesis of cleft lip and palate (CL/P) is studied in animal experiments. This study revealed significant differences in foetal secondary palate development in two strains of mice (NMRI, A/WySnJ) using a palatal organ model. Palatal shelves of 114 NMRI embryos, resistant to cleft occurrence, and 93 A/WySnJ embryos, a strain with a high spontaneous CL/P rate, were micro-dissected at 14.25 GD (gestational day), before palatal fusion takes place. After cultivation in serum-free medium, palatal development was investigated microscopically and scored in a six-step system. At death (14.25 GD) the palatal shelves of the NMRI embryos (mean 3.5) were significant more developed than those of A/WySnJ (mean 2.7; p=0.05). After incubation, 53% (60/114) NMRI and 14% (13/93) A/WySnJ cultures had over two-thirds fusion to stage V-VI, therefore in 17% NMRI (19/114) and 1% A/WySnJ cultures (1/93) fusion was macroscopically complete. 62% of the A/WySnJ cultures showed no significant development in vitro (mean 2.84; p=0.094). There is a significant palatal development difference between normally developed NMRI (mean 4.45, p=0.05) and CL/P appearance in A/WySnJ mice (mean 2.84). Palatal development of both strains was significantly delayed in organ culture (p=0.05). The A/WySnJ strain was more susceptible to manipulation and vulnerable.

An assessment of orofacial clefts in Tanzania

BACKGROUND

Clefts of the lip (CL), the palate (CP), or both (CLP) are the most common orofacial congenital malformations found among live births, accounting for 65% of all head and neck anomalies. The frequency and pattern of orofacial clefts in different parts of the world and among different human groups varies widely. Generally, populations of Asian or Native American origin have the highest prevalence, while Caucasian populations show intermediate prevalence and African populations the lowest. To date, little is known regarding the epidemiology and pattern of orofacial clefts in Tanzania.

METHODS

A retrospective descriptive study was conducted at Bugando Medical Centre to identify all children with orofacial clefts that attended or were treated during a period of five years. Cleft lip and/or palate records were obtained from patient files in the Hospital's Departments of Surgery, Paediatrics and medical records. Age at presentation, sex, region of origin, type and laterality of the cleft were recorded. In addition, presence of associated congenital anomalies or syndromes was recorded.

RESULTS

A total of 240 orofacial cleft cases were seen during this period. Isolated cleft lip was the most common cleft type followed closely by cleft lip and palate (CLP). This is a departure from the pattern of clefting reported for Caucasian and Asian populations, where CLP or isolated cleft palate is the most common type. The distribution of clefts by side showed a statistically significant preponderance of the left side (43.7%) (χ2 = 92.4, p < 0.001), followed by the right (28.8%) and bilateral sides (18.3%). Patients with isolated cleft palate presented at very early age (mean age 1.00 years, SE 0.56). Associated congenital anomalies were observed in 2.8% of all patients with orofacial clefts, and included neural tube defects, Talipes and persistent ductus arteriosus.

CONCLUSIONS

Unilateral orofacial clefts were significantly more common than bilateral clefts; with the left side being the most common affected side. Most of the other findings did not show marked differences with orofacial cleft distributions in other African populations.

Three-dimensional morphometric analysis of brain shape in nonsyndromic orofacial clefting

Previous studies report structural brain differences in individuals with nonsyndromic orofacial clefts (NSOFC) compared with healthy controls. These changes involve non-uniform shifts in tissue volume within the cerebral cortex and cerebellum, suggesting that the shape of the brain may be altered in cleft-affected individuals. To test this hypothesis, a landmark-based morphometric approach was utilized to quantify and compare brain shape in a sample of 31 adult males with cleft lip with or without cleft palate (CL/P), 14 adult males with cleft palate only (CPO) and 41 matched healthy controls. Fifteen midline and surface landmarks were collected from MRI brain scans and the resulting 3D coordinates were subjected to statistical shape analysis. First, a geometric morphometric analysis was performed in three steps: Procrustes superimposition of raw landmark coordinates, omnibus testing for group difference in shape, followed by canonical variates analysis (CVA) of shape coordinates. Secondly, Euclidean distance matrix analysis (EDMA) was carried out on scaled inter-landmark distances to identify localized shape differences throughout the brain. The geometric morphometric analysis revealed significant differences in brain shape among all three groups (P < 0.001). From CVA, the major brain shape changes associated with clefting included selective enlargement of the anterior cerebrum coupled with a relative reduction in posterior and/or inferior cerebral portions, changes in the medio-lateral position of the cerebral poles, posterior displacement of the corpus callosum, and reorientation of the cerebellum. EDMA revealed largely similar brain shape changes. Thus, compared with controls, major brain shape differences were present in adult males with CL/P and CPO. These results both confirm and expand previous findings from traditional volumetric studies of the brain in clefting and provide further evidence that the neuroanatomical phenotype in individuals with NSOFC is a primary manifestation of the defect and not a secondarily acquired characteristic.

Face shape of unaffected parents with cleft affected offspring: combining three-dimensional surface imaging and geometric morphometrics

OBJECTIVE

Various lines of evidence suggest that face shape may be a predisposing factor for non-syndromic cleft lip with or without cleft palate (CL/P). In the present study, 3D surface imaging and statistical shape analysis were used to evaluate face shape differences between the unaffected (non-cleft) parents of individuals with CL / P and unrelated controls.

METHODS

Sixteen facial landmarks were collected from 3D captures of 80 unaffected parents and 80 matched controls. Prior to analysis, each unaffected parent was assigned to a subgroup on the basis of prior family history (positive or negative). A geometric morphometric approach was utilized to scale and superimpose the landmark coordinate data (Procrustes analysis), test for omnibus group differences in face shape, and uncover specific modes of shape variation capable of discriminating unaffected parents from controls.

RESULTS

Significant disparity in face shape was observed between unaffected parents and controls (p < 0.01). Notably, these changes were specific to parents with a positive family history of CL/P. Shape changes associated with CL/P predisposition included marked flattening of the facial profile (midface retrusion), reduced upper facial height, increased lower facial height, and excess interorbital width. Additionally, a sex-specific pattern of parent-control difference was evident in the transverse dimensions of the nasolabial complex.

CONCLUSIONS

The faces of unaffected parents from multiplex cleft families displayed meaningful shape differences compared with the general population. Quantitative assessment of the facial phenotype in cleft families may enhance efforts to discover the root causes of CL/P.

Deciphering the Palimpsest: Studying the Relationship Between Morphological Integration and Phenotypic Covariation

Organisms represent a complex arrangement of anatomical structures and individuated parts that must maintain functional associations through development. This integration of variation between functionally related body parts and the modular organization of development are fundamental determinants of their evolvability. This is because integration results in the expression of coordinated variation that can create preferred directions for evolutionary change, while modularity enables variation in a group of traits or regions to accumulate without deleterious effects on other aspects of the organism. Using our own work on both model systems (e.g., lab mice, avians) and natural populations of rodents and primates, we explore in this paper the relationship between patterns of phenotypic covariation and the developmental determinants of integration that those patterns are assumed to reflect. We show that integration cannot be reliably studied through phenotypic covariance patterns alone and argue that the relationship between phenotypic covariation and integration is obscured in two ways. One is the superimposition of multiple determinants of covariance in complex systems and the other is the dependence of covariation structure on variances in covariance-generating processes. As a consequence, we argue that the direct study of the developmental determinants of integration in model systems is necessary to fully interpret patterns of covariation in natural populations, to link covariation patterns to the processes that generate them, and to understand their significance for evolutionary explanation.

A novel 3-D image-based morphological method for phenotypic analysis

A new approach for the study of geometric morphometrics is presented based on well-established image processing techniques in a novel combination to support high-throughput analysis necessary for large-scale determination of genotype-phenotype relationships. The method retains full 3-D data, and avoids manual landmark selection. Micro-computed tomography images are superimposed into a common orientation by rigid image registration with an isotropic scale factor. An average sample shape is determined by averaging the intensities of corresponding voxels of the registered images, and shape variation is determined by calculating the image gradient of the average shape. Localized shape differences between mean images or between an individual and a group mean are identified and quantified by surface-to-surface distance measures of superimposed images. Validation was performed using geometric shapes of known dimensions as well as biological samples of C57 BL/6 J and A/WySnJ mouse skulls, and shape variation of the mouse skulls was consistent with previously published results. Although the image gradient is sensitive to both image registration and filtration of the average image, the effect can be minimized by consistent use of image analysis parameters. While the proposed approach deviates from well-established landmark-based geometric morphometric tools, it is not intended to replace these current methods. Rather, it will be an important contribution to provide high-throughput screening in large-scale studies focused on understanding genotype-phenotype relationships so that subsequent morphometric approaches using established techniques can be better focused.

Short-faced mice and developmental interactions between the brain and the face

The length of the face represents an important axis of variation in mammals and especially in primates. Mice with mutations that produce variation along this axis present an opportunity to study the developmental factors that may underlie evolutionary change in facial length. The Crf4 mutant, obtained from the C57BL/6J (wt/wt) background by chemical mutagenesis by the Baylor Mouse Mutagenesis Resource, is reported to have a short-faced phenotype. As an initial step towards developing this model, we performed 3D geometric morphometric comparisons of Crf4 mice to C57BL/6J wild-type mice focusing on three stages of face development and morphology--embryonic (GD 9.5-12), neonatal, and adult. Morphometric analysis of adult Crf4 mutants revealed that in addition to a shortened face, these mice exhibit a significant reduction in brain size and basicranial length. These same features also differ at the neonatal stage. During embryonic face formation, only dimensions related to brain growth were smaller, whereas the Crf4 face actually appeared advanced relative to the wild-type at the same somite stage. These results show that aspects of the Crf4 phenotype are evident as early as embryonic face formation. Based on our anatomical findings we hypothesize that the reduction in facial growth in Crf4 mice is a secondary consequence of reduction in the growth of the brain. If correct, the Crf4 mutant will be a useful model for studying the role of epigenetic interactions between the brain and face in the evolutionary developmental biology of the mammalian craniofacial complex as well as human craniofacial dysmorphology.

Spatial packing, cranial base angulation, and craniofacial shape variation in the mammalian skull: testing a new model using mice

The hypothesis that variation in craniofacial shape within and among species is influenced by spatial packing has a long history in comparative anatomy, particularly in terms of primates. This study develops and tests three alternative models of spatial packing to address how and to what extent the cranial base angle is influenced by variation in brain and facial size. The models are tested using mouse strains with different mutations affecting craniofacial growth. Although mice have distinctive crania with small brains, long faces, and retroflexed cranial bases, the results of the study indicate that the mouse cranial base flexes to accommodate larger brain size relative to cranial base length. In addition, the mouse cranial base also extends, but to a lesser degree, to accommodate larger face size relative to cranial base length. In addition, interactions between brain size, face size, and the widths and lengths of the components of the cranial base account for a large percentage of variation in cranial base angle. The results illustrate the degree to which the cranial base is centrally embedded within the covariation structure of the craniofacial complex as a whole.

Quantification and variation in experimental studies of morphogenesis

The application of quantitative methods has long been the norm in comparative and evolutionary studies of morphology, but within the field of experimental embryology mathematical descriptions of anatomical form are seldom calculated, and morphological variation within treatment groups is rarely taken into account. Here we argue that many of the analytical techniques that are commonly applied in other areas of morphological research are also well suited for experimental studies of anatomical development. The application of these methodologies shows promise for augmenting such endeavors by enhancing researchers' ability to detect morphological patterns, account for developmental variation, and employ statistical methods. We review selected studies of experimental morphogenesis that underscore the potential of quantitative methods to reveal important aspects of anatomical development and growth. These examples demonstrate the benefits of quantifying ontogenetic data and accounting for developmental variation, and we suggest that the adoption of such practices by researchers performing experimental studies of morphogenesis will enhance our understanding of the processes by which genetic changes affect anatomical formation.

Three-dimensional morphometric analysis of craniofacial shape in the unaffected relatives of individuals with nonsyndromic orofacial clefts: a possible marker for genetic susceptibility

Numerous studies have described altered patterns of craniofacial form in the unaffected relatives of individuals with nonsyndromic clefts. Unfortunately, results from such studies have been highly variable and have failed to provide a reliable method for differentiating "at-risk" relatives from controls. In the present study, we compared craniofacial shape between a sample of unaffected relatives (33 females; 14 males) from cleft multiplex families and an equal number of age/sex/ethnicity-matched controls. Sixteen x,y,z facial landmark coordinates derived from 3D photogrammetry were analyzed via Euclidean Distance Matrix Analysis, while 14 additional linear distances were analyzed via t tests. A subset of variables was then entered into a discriminant function analysis (DFA). Compared to controls, female unaffected relatives demonstrated increased upper facial width, midface reduction and lateral displacement of the alar cartilage. DFA correctly classified 70% of female unaffected relatives and 73% of female controls. Male unaffected relatives demonstrated increased upper facial and cranial base width, increased lower facial height and decreased upper facial height compared with controls. DFA correctly classified 86% of male unaffected relatives and 93% of male controls. In both sexes, upper facial width contributed most to group discrimination. Following DFA, unaffected relatives were assigned to risk/liability classes based on the degree of phenotypic divergence from controls. Results indicate that craniofacial shape differences characterizing unaffected relatives are partly sex-specific and are in broad agreement with previous reports. These findings further suggest that a quantitative assessment of the craniofacial phenotype may allow for the identification of susceptible individuals within nonsyndromic cleft families.

Cranial bone morphometric study among mouse strains

BACKGROUND

Little is known about the molecular mechanism which regulates how the whole cranium is shaped. Mouse models currently available for genetic research include several hundreds of unique inbred strains and genetically engineered mutants. By cross comparing their genomic structures, we can elucidate the cause of any differences in the phenotype between two strains. The craniometry of subspecies, or closely related species, of mice provide a good systemic model to study the relationship between genetic variance and cranial shape evolution. The lack of a quantified framework for comparing and analyzing mouse cranial shape has been a problem. For this reason, we performed quantitative analysis of cranial shape morphology between several mouse strains.

RESULTS

This article reports on a craniometric assay of seven mouse strains: four inbred strains (C57BL/6J, BALB/cA, C3H/HeJ, and CBA/JNCr) from Mus musculus domesticus (M. m. domesticus); one closed colony strain (ICR) from M. m. domesticus; one inbred strain (MSM/Ms) from Mus musculus molossinus; and, Mus spretus as a strain from a species other than M. m. domesticus. We performed linear measurements and geometric morphometrics. Geometric morphometrics revealed that the cranial characteristics of each strains were clearly distinguishable. We obtained mean scores for each species using the tpsRelw Program and plotted them.

CONCLUSION

Geometric morphometrics proved to be useful for identifying and classifying variations in form, and it revealed that M. spretus has a slender cranium when compared with our other strains. The mean cranial shape of C3H or CBA was more similar to MSM/Ms, which is derived from M. m. molossinus, than to either C57BL/6J, BALB, or ICR which are derived from M. m. domesticus. Future work in this field will aid in elucidating the mechanism of whole cranial shape regulation.

Phenotypic variability and craniofacial dysmorphology: increased shape variance in a mouse model for cleft lip

Cleft lip and palate (CL/P), as is true of many craniofacial malformations in humans, is etiologically complex and highly variable in expression. A/WySn mice are an intriguing model for human CL/P because they develop this dysmorphology with a variable expression pattern, incomplete penetrance and frequent unilateral expression on a homogeneous genetic background. The developmental basis for this variation in expression is unknown, but of great significance for understanding such expression patterns in humans. As a step towards this goal, this study used three-dimensional geometric morphometric and novel high throughput morphometric techniques based on three-dimensional computed microtomography of mouse embryos to analyze craniofacial shape variation during primary palate formation. Our analysis confirmed previous findings based on two-dimensional analyses that the midface in A/WySn embryos, and the maxillary prominence in particular, is relatively reduced in size and appears to be developmentally delayed. In addition, we find that shape variance is increased in A/WySn embryos during primary palate formation compared to both C57BL/6J mice and the F1 crosses between these strains. If the reduction in midfacial growth caused by the Wnt9b hypomorphic mutation pushes A/WySn mice closer on average to the threshold for cleft lip formation, the elevated shape variance may explain why some, but not all, embryos develop the dysmorphology in a genetically homogeneous inbred line of mice.

Genetics of microenvironmental canalization in Arabidopsis thaliana

Canalization is a fundamental feature of many developmental systems, yet the genetic basis for this property remains elusive. We examine the genetic basis of microenvironmental canalization in the model plant Arabidopsis thaliana, focusing on differential developmental stability between genotypes in one fitness and four quantitative morphological traits. We measured developmental stability in genetically identical replicates of two populations of recombinant inbred (RI) lines and one population of geographically widespread accessions of A. thaliana grown in two different photoperiod-controlled environments. We were able to map quantitative trait loci associated with developmental stability. We also identified a candidate gene, ERECTA, that may contribute to microenvironmental canalization in rosette leaf number under long-day photoperiods, and analysis of mutant lines indicates that the er-105 allele results in increased canalization for this trait. ERECTA, which encodes a signaling protein, appears to act as an ecological amplifier by transducing developmental noise (e.g., microenvironmental variation) into phenotypic differentiation. We also measured genotypic selection on four plant architecture traits and find evidence for selection for both increased and decreased canalization at various traits.

Comparative morphometrics of embryonic facial morphogenesis: implications for cleft-lip etiology

Cleft lip (CL) with or without cleft palate (CL[P]) has a complex etiology but is thought to be due to either genetic or environmentally induced disruptions of developmental processes affecting the shape and size of the facial prominences (medial nasal, lateral nasal, and maxilla). Recent advances in landmark-based morphometrics enable a rigorous reanalysis of phenotypic shape variation associated with facial clefting. Here we use geometric morphometric (GM) tools to characterize embryonic shape variation in the midface and head of six strains of mice that are both cleft-liable (A, A/WySn, CL/Fr) and normal (BALB/cBy, C57BL, CD1). Data were comprised of two-dimensional landmarks taken from frontal and lateral photographs of embryos spanning the time period in which the facial prominences fuse (GD10-12). Results indicate that A/- strain mice, and particularly A/WySn, have overall smaller midfaces compared to other strains. The A/WySn strain also has significant differences in facial shape related to retarded development. Overall, CL/Fr strain mice are normal-sized, but tend to have undersized maxillary prominences that do not project anteriorly and have a small nasal contact area. These results suggest that the etiology of clefting differs in A/WySn and CL/Fr strains, with the former strain suffering disruptions to developmental processes affecting overall size (e.g., neural crest migration deficiencies and lower mitotic activity), while the latter strain has defects restricted to the shape and size of the maxilla. A combination of molecular experimentation and phenotypic analysis of shape is required to test these hypotheses further.

Epigenetic interactions and the structure of phenotypic variation in the cranium

Understanding the developmental and genetic basis for evolutionarily significant morphological variation in complex phenotypes such as the mammalian skull is a challenge because of the sheer complexity of the factors involved. We hypothesize that even in this complex system, the expression of phenotypic variation is structured by the interaction of a few key developmental processes. To test this hypothesis, we created a highly variable sample of crania using four mouse mutants and their wild-type controls from similar genetic backgrounds with developmental perturbations to particular cranial regions. Using geometric morphometric methods we compared patterns of size, shape, and integration in the sample within and between the basicranium, neurocranium, and face. The results highlight regular and predictable patterns of covariation among regions of the skull that presumably reflect the epigenetic influences of the genetic perturbations in the sample. Covariation between relative widths of adjoining regions is the most dominant factor, but there are other significant axes of covariation such as the relationship between neurocranial size and basicranial flexion. Although there are other sources of variation related to developmental perturbations not analyzed in this study, the patterns of covariation created by the epigenetic interactions evident in this sample may underlie larger scale evolutionary patterns in mammalian craniofacial form.

Epidermal growth factor receptor and transforming growth factor-beta signaling contributes to variation for wing shape in Drosophila melanogaster

Wing development in Drosophila is a common model system for the dissection of genetic networks and their roles during development. In particular, the RTK and TGF-beta regulatory networks appear to be involved with numerous aspects of wing development, including patterning, cell determination, growth, proliferation, and survival in the developing imaginal wing disc. However, little is known as to how subtle changes in the function of these genes may contribute to quantitative variation for wing shape, per se. In this study 50 insertional mutations, representing 43 loci in the RTK, Hedgehog, TGF-beta pathways, and their genetically interacting factors were used to study the role of these networks on wing shape. To concurrently examine how genetic background modulates the effects of the mutation, each insertion was introgressed into two wild-type genetic backgrounds. Using geometric morphometric methods, it is shown that the majority of these mutations have profound effects on shape but not size of the wing when measured as heterozygotes. To examine the relationships between how each mutation affects wing shape hierarchical clustering was used. Unlike previous observations of environmental canalization, these mutations did not generally increase within-line variation relative to their wild-type counterparts. These results provide an entry point into the genetics of wing shape and are discussed within the framework of the dissection of complex phenotypes.

Canalization and developmental stability in the Brachyrrhine mouse

The semi-dominant Br mutation affects presphenoid growth, producing the facial retrognathism and globular neurocranial vault that characterize heterozygotes. We analysed the impact of this mutation on skull shape, comparing heterozygotes to wildtype mice, to determine if the effects are skull-wide or confined to the sphenoid region targeted by the mutation. In addition, we examined patterns of variability of shape for the skull as a whole and for three regions (basicranium, face and neurocranium). We found that the Br mice differed significantly from wildtype mice in skull shape in all three regions as well as in the shape of the skull as a whole. However, the significant increases in variance and fluctuating asymmetry were found only in the basicranium of mutant mice. These results suggest that the mutation has a significant effect on the underlying developmental architecture of the skull, which produces an increase in phenotypic variability that is localized to the anatomical region in which the mean phenotype is most dramatically affected. These results suggest that the same developmental mechanisms that produce the change in phenotypic mean also produce the change in variance.

The brachymorph mouse and the developmental-genetic basis for canalization and morphological integration

Although it is well known that many mutations influence phenotypic variability as well as the mean, the underlying mechanisms for variability effects are very poorly understood. The brachymorph (bm) phenotype results from an autosomal recessive mutation in the phosphoadenosine-phosphosulfate synthetase 2 gene (Papps2). A major cranial manifestation is a dramatic reduction in the growth of the chondrocranium which results from undersulfation of glycosaminoglycans (GAGs) in the cartilage matrix. We found that this reduction in the growth of the chondrocranium is associated with an altered pattern of craniofacial shape variation, a significant increase in phenotypic variance and a dramatic increase in morphological integration for craniofacial shape. Both effects are largest in the basicranium. The altered variation pattern indicates that the mutation produces developmental influences on shape that are not present in the wildtype. As the mutation dramatically reduces sulfation of GAGs, we infer that this influence is variation among individuals in the degree of sulfation, or variable expressivity of the mutation. This variation may be because of genetic variation at other loci that influence sulfation, environmental effects, or intrinsic effects. We infer that chondrocranial development exhibits greater sensitivity to variation in the sulfation of chondroitin sulfate when the degree of sulfation is low. At normal levels, sulfation probably contributes minimally to phenotypic variation. This case illustrates canalization in a particular developmental-genetic context.