Facial phenotypes in subgroups of prepubertal boys with autism spectrum disorders are correlated with clinical phenotypes
© Aldridge et al; licensee BioMed Central Ltd. 2011
Received: 27 May 2011
Accepted: 14 October 2011
Published: 14 October 2011
The brain develops in concert and in coordination with the developing facial tissues, with each influencing the development of the other and sharing genetic signaling pathways. Autism spectrum disorders (ASDs) result from alterations in the embryological brain, suggesting that the development of the faces of children with ASD may result in subtle facial differences compared to typically developing children. In this study, we tested two hypotheses. First, we asked whether children with ASD display a subtle but distinct facial phenotype compared to typically developing children. Second, we sought to determine whether there are subgroups of facial phenotypes within the population of children with ASD that denote biologically discrete subgroups.
The 3dMD cranial System was used to acquire three-dimensional stereophotogrammetric images for our study sample of 8- to 12-year-old boys diagnosed with essential ASD (n = 65) and typically developing boys (n = 41) following approved Institutional Review Board protocols. Three-dimensional coordinates were recorded for 17 facial anthropometric landmarks using the 3dMD Patient software. Statistical comparisons of facial phenotypes were completed using Euclidean Distance Matrix Analysis and Principal Coordinates Analysis. Data representing clinical and behavioral traits were statistically compared among groups by using χ2 tests, Fisher's exact tests, Kolmogorov-Smirnov tests and Student's t-tests where appropriate.
First, we found that there are significant differences in facial morphology in boys with ASD compared to typically developing boys. Second, we also found two subgroups of boys with ASD with facial morphology that differed from the majority of the boys with ASD and the typically developing boys. Furthermore, membership in each of these distinct subgroups was correlated with particular clinical and behavioral traits.
Boys with ASD display a facial phenotype distinct from that of typically developing boys, which may reflect alterations in the prenatal development of the brain. Subgroups of boys with ASD defined by distinct facial morphologies correlated with clinical and behavioral traits, suggesting potentially different etiologies and genetic differences compared to the larger group of boys with ASD. Further investigations into genes involved in neurodevelopment and craniofacial development of these subgroups will help to elucidate the causes and significance of these subtle facial differences.
Autism is a spectrum of disorders (ASDs) united by a common Diagnostic and Statistical Manual of Mental Disorders, Fourth Edition (DSM-IV)-defined  behavioral phenotype. Research into this disorder is increasingly focused on both genetic causes and neuroanatomical bases for the behavioral phenotypes. Thus far attempts to discover major autism susceptibility genes have been largely unsuccessful, with approximately only 15% to 20% of cases of autism linked to specific genes, chromosomal aneuploidy or recognized syndromes [2–4]. The rest remain idiopathic. A primary reason for the lack of progress in understanding the etiology and genetic underpinnings of ASD is undoubtedly the significant heterogeneity within both behavioral and clinical phenotypes.
The neurodevelopmental model of ASD  suggests that changes in embryonic developmental patterns result in the spectrum of ASD phenotypes and that these changes may result from permutations of genetics, the environment, or the interaction of the two. This model posits that the brain is altered during embryonic development, a time when the brain is intimately tied to developing facial tissues via genetic signaling, biomechanical and biochemical mechanisms [6–13]. The face and brain form a vast but segmented population of cells whose differentiation and identity are established through intricate signaling mechanisms. It has been stated that the brain is the foundation on which the various parts of the developing face grow . The phrase, "The face predicts the brain" , has been employed frequently to explain developmental disorders such as holoprosencephaly. Thus changes to the developing brain may be reflected in the face [15–17].
The face develops from populations of neural crest cells migrating from the neural tube into developing embryonic facial prominences. These neural crest cells interact with the developing brain via both physical contact and genetic signaling. Previous research has shown that the expression of Sonic hedgehog (SHH), fibroblast growth factor 8 (FGF8) and bone morphogenetic proteins (BMPs outline developmental interactions between the face and the brain in discrete temporal and regional patterns [18–22]. The connection of brain and facial phenotypes to their underlying genetic bases involves hierarchies of complex regulatory cascades, nested epigenetic networks and ever-changing patterns of cross-talk between molecules, cells and tissues throughout development [6, 7, 12, 23]. Localization of facial phenotypic variations to specific areas of the face may reveal potential candidate genes and/or pathways targeted in the development of the brain in autism.
Given the clear evidence that the embryological face and anterior brain emerge and develop in exquisite intimacy [16, 17], facial phenotypes can serve as accessible and informative indices of brain phenotypes in neurodevelopmental disorders. We propose that altered expression of genes involved in the development of neural tube structures and overlying facial prominences may result in distinct facial and neural phenotypes in ASD. We hypothesize that there are common autism-causing genes that affect early brain development and simultaneously the facial phenotype. Defining one or more common facial phenotypes within ASD will provide a new physical biomarker that can be used to improve ASD diagnoses, with all the associated benefits related to prognosis, recurrence counseling, choice of subgroup-appropriate therapies and the possibility of developing a screening tool to assist in early diagnosis.
Precise measures of phenotypes and innovative methods of analysis are integral to discovering the nature of developmental contributions to phenotypic variation. Using state-of-the-art three-dimensional photography to obtain facial images to precisely measure facial phenotypes, we tested two hypotheses. First, children with ASD display a subtle but distinct facial phenotype compared to typically developing (TD) children. Second, there are subgroups of facial phenotypes within the population of children with ASD that denote biologically discrete subgroups.
Previous work has suggested that there are autism facial phenotypes with a developmental basis. In an epidemiological study of facial photographs of children with autism and developmental disabilities, Rodier and colleagues [24, 25] reported a facial phenotype common in autism consisting of decreased interpupillary distance (although intercanthic distance was not decreased), ptosis, strabismus, lop ears and hypotonia of the lower face. They postulated that these minor anomalies arise as the face is closing and the cranial nerves are invading mesenchyme that will develop into the muscle, skeletal and dermal tissues of the head. Rodier and colleagues [26, 27] also suggested that facial phenotype might allow researchers to pick out the children whose autism is due to mutations in the homeobox genes which control the development of both the brainstem and face. Though this group of genes is important for embryological development, we now know that the face and the rostral brain are not patterned by genes in the homeobox family . Additionally, Hammond et al.  studied a group of boys, ages 2 to 18 years from families with at least two affected family members, finding minor shape differences in comparisons of the mean facial phenotype of these boys to a control mean facial phenotype, with their major findings emphasizing significant facial asymmetry in boys with ASD and their family members and suggesting a shared developmental basis for these phenotypes.
Study sample and recruitment
Study sample and age range
Age in years (mean ± SD)
Boys with ASD
10.14 ± 1.318
Typically developing boys
10.19 ± 1.261
Participants with ASD (n = 64) were recruited through the Thompson Center for Autism and Neurodevelopmental Disorders. All participants were screened before inclusion and met the following criteria: individuals were male; of Caucasian ethnicity; had not worn dental braces; were prepubertal (by parent report); were able to sit relatively still for picture-taking; had been diagnosed with Autistic Disorder, Asperger syndrome, or pervasive developmental disorder-not otherwise specified (PDD-NOS) according to the DSM-IV criteria prior to the day of the study; and had no additional syndrome diagnoses. Boys with fragile X syndrome and/or chromosomal disorders, including copy number variants (CNV), generalized dysmorphology or gestational age less than 35 weeks were excluded.
Of the 64 boys with ASD, 36 had completed the Simons Simplex Collection (SSC) protocol, which includes the Autism Diagnostic Interview, Revised (ADI-R)  and the Autism Diagnostic Observation Schedule (ADOS) , which were used in conjunction with the clinical judgment of one of the authors (JHM) to make the diagnosis of ASD. The 28 boys recruited through the Autism Medical Clinic were diagnosed on the basis of the DSM-IV criteria using a center-specific protocol based on the ADI-R together with the clinical judgment of the same author (JHM). The boys were assessed for generalized dysmorphology using the Autism Dysmorphology Measure .
TD boys (n = 41) were recruited from the Columbia, MO, USA, community via a notice published in the University of Missouri online information email and by word of mouth. Participants were screened using the same criteria described in the preceding paragraph, with the exception of a diagnosis of ASD. We chose TD boys as the control group, as our hypothesis was that development in ASD deviates from normal development. Samples were not matched for IQ, since this would not have allowed us to interpret how ASD deviates from the normal developmental trajectory. Recruitment and data collection procedures were carried out in accordance with approved Institutional Review Board protocols.
Three-dimensional stereophotogrammetric imaging
Prior to each session, all cameras were calibrated and tested to ensure that the images collected were consistent and usable. Individuals participating in the study were brought into the camera area and asked to sit as still as possible, look directly at one camera marked with a sticker and maintain a neutral expression (verbal instructions included closed mouth, no smile, no visible teeth and no raised eyebrows). Once the participant was comfortable and able to sit still, collection of the images began. Multiple pictures of each child were taken to ensure that the image used for analysis adequately captured all of the facial areas needed for landmarking.
Anthropometric landmark data collection
Three-dimensional anthropometric landmarks acquired from 3dMD images
To determine the reliability of data collection, an error study was performed. Four trials of landmark coordinate data were collected from two 3dMD images by both raters. Coordinate data were converted to all possible linear distances among the landmarks (all linear distances between 17 landmarks, resulting in 136 linear distances). Means, standard deviations and values of standard deviations as percentages of the linear distances were calculated for each linear distance in each observer's trials. The results derived by both raters are presented.
The ranges of the standard deviations expressed as percentages of linear distance were 0.19% to 8.11% for rater 1 and 0.19% to 14.3% for rater 2. Of the 136 total linear distances evaluated, linear distances with standard deviations greater than 5% of the mean totaled seven for rater 1 and eight for rater 2, leaving 129 and 128 linear distances, respectively, with less than a 5% error for each rater, respectively. Of the seven and eight linear distances with greater degrees of error, three of them were shared by the two raters.
The results of this error study indicate that the landmark coordinate data and the linear distances calculated from them can be collected with a very low degree of error. Therefore, the data collected by the two raters in this study are highly precise and repeatable.
Morphometric data analysis
The landmark coordinate data collected using 3dMDpatient software were analyzed using Euclidean Distance Matrix Analysis (EDMA) , which is a linear distance-based morphometric method that does not rely on registration or fitting criteria [43, 44]. Linear distances calculated between all possible pairs of landmarks were compared across samples as ratios. EDMA represents the form of each individual as a form matrix (FM), which is the set of all possible linear distances between the facial landmarks. Average FMs for each sample, that is, ASD and controls, are compared as ratios of like linear distances. This set of ratios of corresponding linear distances is called a "form difference matrix" (FDM). If a ratio in a FDM is equal to 1, then the faces being compared do not differ for that discrete linear distance. If the ratio is above 1, the linear distance is greater in the face used as the numerator. Likewise, if the ratio is less than 1, the linear distance is greater in the face used as the denominator. We used a nonparametric bootstrapping algorithm to calculate confidence intervals for each discrete linear distance to test for the significance of localized form differences . The null hypothesis is that each discrete linear distance is similar for the two samples. Individual linear distances were considered significantly different if the calculated two-tailed 90% confidence interval did not include 1.0. Evaluation of confidence intervals for differences in specific linear distances enables localization of differences to specific facial regions.
This test of empirical differences in shape between samples is based on marginal confidence intervals of the bootstrap estimates of the linear distances between unique pairs of landmarks. Bonferroni-type corrections are not needed for these marginal confidence intervals, because in this approach multiple tests of linear distance differences using the same data are not conducted. Instead, with each bootstrapping step, all measures are estimated for an individual and tested in a high-dimensional space where each dimension represents a unique linear distance. The low-dimensional projection of these results for each linear distance is reported (see [43, 45] for details).
This method has been used in numerous previous studies to compare facial morphologies (for example, see [42–49]) as well as the morphologies of other anatomic regions. A validation of this method for the data set in this study was performed to ensure that differences found in comparing the boys with ASD to TD boys using EDMA were not spurious. The group of TD boys was split into two randomly assigned age-equivalent groups. These two groups were then compared, and confidence intervals were calculated for each linear distance. On the basis of the results of these analyses, we determined that there were three significantly different linear distances among the total of 136 compared (2.2% of 136). These results show that fewer were statistically different than would be expected by chance (that is, 5%), demonstrating that this method is both sensitive and specific.
Principal Coordinates Analysis (PCOORD) application of EDMA was then performed on the scaled data for all participants in both groups [43–45, 50]. This procedure is a form of clustering analysis that detects groups of forms with similar shapes and identifies linear distances that are influential in forming the defined clusters. In this procedure, the distribution of participants in multidimensional morphological space is examined. Unlike the form difference analyses described above, PCOORD compares individuals rather than samples. Axes are fitted through the shape space of this analysis such that the first axis accounts for the majority of the variation, the second axis accounts for the second-largest amount of variation and so on. These axes are referred to as the "first principal axis," the "second principal axis," and so forth. The position of participants along these axes is defined in terms of the linear distances between landmarks most highly correlated with these axes. Therefore, participants who cluster along a particular axis are similar in terms of the linear distances correlated with that axis. This analysis was performed to determine whether participants clustered on the basis of facial morphology, to identify the metrics that contributed to determination of the clusters and to explore the nature of the clusters to formulate hypotheses about the pattern of differences in the development of the face in children with ASD. The PCOORD analyses were performed on data that were scaled for differences in size. To do this, the FM for each individual was scaled such that each linear distance was divided by the geometric mean of all linear distances within that individual's FM. Thus each participant's data were scaled using a unique scaling factor. The geometric mean was chosen as a surrogate for size [51–53].
We analyzed all of the participants in the study to determine (1) whether there are aspects of facial morphology that distinguish the facial phenotypes of boys with ASD compared to TD boys and (2) whether there are facial subgroups within the ASD cohort that differ in their associated clinical and behavioral parameters.
Behavioral and medical data
Each of the boys was evaluated for characteristics of their ASD diagnosis (social function, verbal function, repetitive behavior and language level), behavioral problems (aggression, attention deficits and self-injurious behaviors), outcome measures (IQ, communication, daily living skills, socialization and Vineland Adaptive Behavior Scale composite scores), the clinical course of their disorder (age at onset and presence of regression at onset), medical and neurological variables (seizures, electroencephalogram results, hypotonia, hypertonia, clumsiness, vision or hearing problems, tics, enuresis, handedness, feeding difficulties in infancy and allergies), physical morphology (head circumference, height, weight and dysmorphology) and family history of autism and related neuropsychiatric disorders among first-degree relatives.
The tests administered to all or the majority of participants included the ADI-R , ADOS , Social Communication Questionnaire (SCQ) , Vineland Adaptive Behavior Scale II , Peabody Picture Vocabulary Test (PPVT) , Child Behavior Checklist (CBCL) , an age- and development-appropriate IQ test (Full Scale IQ (FSIQ), Verbal IQ (VIQ), Nonverbal IQ (NVIQ)) and the Autism Dysmorphology Measure . Not all measures of IQ were available for a small number of boys. NVIQ was available for the entire sample, FSIQ was available for all but one boy and VIQ was missing for four of the boys. Comprehensive prenatal, perinatal, teratogen exposure, development, general health, neurological and family histories (including income and education), were obtained using either the SSC Medical History or the Thompson Center Medical History, which record similar information. All participants received complete medical and neurological examinations, including assessment of growth and dysmorphology.
Statistical comparisons of facial phenotypes with clinical and behavioral phenotypes
We compared clinical and behavioral traits to determine whether there were significant correlations between subgroup membership and the variables described in the preceding subsection. Continuous random variables were summarized by their mean, standard deviation and range. For categorical random variables, univariate comparisons of subgroup 1, subgroup 2 and the remainder were made using the χ2 test or Fisher's exact test. For continuous variables, comparisons were made using Student's t-test. Because the IQ score is skewed, comparisons of IQ scores were made using the Kolmogorov-Smirnov test and the t-test after log-transforming IQ.
Euclidean Distance Matrix Analysis
We did not find a statistically significant difference when we compared the face sizes of ASD and TD boys (t-test; P = 0.301). This suggests that the morphological differences found in the EDMA comparisons were not due to differences in size.
Principal Coordinates Analysis
The second principal axis accounted for 11.37% of the variance in the study population. Seven of the individuals with ASD (7.8% of the ASD study group) clustered at the strongly negative end of this axis. The linear distances that were reduced in this subgroup relative to the remainder of the ASD sample population spanned the area from the inferior aspect of the nasal region to the philtrum and lateral mouth. The linear distances that were increased in this subgroup relative to the remainder of the ASD sample population (1) spanned the breadth of the upper face and (2) connected the lateral aspect of the upper face with the inferior aspect of the nasal region (Figure 6B).
Correlations between axis score and age (r2 = 0.125 for axis 1 and r2 = 0.068 for axis 2) and between axis score and head size as measured by the geometric mean of all possible linear distances for each individual (r2 = 0.285 for axis 1 and r2 = 0.039 for axis 2) were low. These results suggest that group membership is not due to differences in age or head size.
Comparisons of clinical and behavioral phenotypes among facial morphology subgroups
Clinical and behavioral characteristics of autism spectrum disorder facial subgroups
Subgroup 1 n= 12 (18.8%)
Subgroup 2 n= 5 (7.8%)
Remainder n= 47 (73.4%)
SCQ lifetime, mean (SD)
FSIQ, mean (SD)
36 to 113
38 to 127
31 to 130
VIQ, mean (SD)
13 to 122
23 to 121
28 to 138
NVIQ, mean (SD)
34 to 119
53 to 121
33 to 129
Mean age at first word month (SD)
Mean age at first phrase month (SD)
CBCL total problems t-score ≥65
Infant feeding problems
SRS mother, mean (SD)
No differences were found on the language measures, including the PPVT scores, Vineland Adaptive Behavior Scale II mean communication scores, ADI-R B language scores and the ADOS modules used. For all of the language measures, however, subgroup 2 subjects scored slightly better than subgroup 1. The age at first words and phrases illustrates the trend toward earlier language development in subgroup 2. No significant differences were found in the Vineland Adaptive Behavior Scale II scores, including Communication, Daily living skills, Socialization and Composite. No differences were seen for demographic or socioeconomic variables, including parent education or income. Since all subjects were male and 56% were also enrolled in the SSC, genetic indicators could be assessed. To be eligible for the SSC, there could be no history of autism among first- or second-degree relatives and no close relatives with major neuropsychiatric disorders. The only significant difference in the subgroups was a significantly higher Social Responsiveness Scale (SRS) mean score for mothers in subgroup 1.
We proposed that ASD might be associated with a distinctive facial phenotype and hypothesized that we might be able to identify subgroups within ASD based on facial phenotypes. Using sophisticated facial phenotyping based on three-dimensional stereophotogrammetric imaging and advanced statistical analyses, we developed a research methodology that allowed us to test and garner support for both hypotheses. In addition, we maximized our chances of success by selecting a relatively homogeneous ASD population that included only boys of self-reported Caucasian ethnicity, with essential autism, defined as having no discernible dysmorphology or microcephaly and within a limited age range.
We found that essential autism in boys is associated with a distinctive facial phenotype characterized by an increased breadth of the mouth, orbits and upper face, combined with a flattened nasal bridge and reduced height of the philtrum and maxillary region. This facial phenotype is similar to the one we recognized clinically (JHM) and may be the "beautiful face" mentioned by Kanner . Embryologically, this facial phenotype is indicative of a perturbation of the FNP. It can be explained by reduction in the superoinferior dimension of the midline structures derived from the FNP, an increase in the subnasal portions of the FNP and a concomitant increase in the breadth of the upper face.
The common facial phenotype described by Rodier and colleagues  includes a reduced interpupillary distance with no difference in intercanthic distance. In contrast, in our present study, we found a narrowing of the intercanthic distance, or mild hypotelorism. These findings are complementary in that we found an overall decrease in intercanthic distance, which potentially translates to decreased interpupillary distance, although we did not directly measure that distance. Our findings are also in line with those of Hammond et al. . However, the sample included in the Hammond et al.  study consisted of a group of boys, ages 2 to 18 years from families with at least two affected family members. Our study extends the findings of both Rodier and colleagues  and Hammond et al.  by quantifying precisely localized differences and variations in facial phenotypes in a homogeneous group of boys with essential autism.
Our findings further extend previous work in that we have discovered two subgroups of boys with ASD who displayed unique facial phenotypes, which correspond to distinct clinical phenotypes, compared to both the majority of boys with ASD and TD boys. Subgroup 1 displays decreased height of the facial midline and increased breadth of the mouth as well as the length and height of the chin. These regions of the face develop primarily from the FNP and midline portions of the MAND prominences of the embryonic face. Subgroup 2 displays increased breadth of the upper face in combination with decreased height of the philtrum. Both of these regions develop from the embryonic FNP.
The results of the tests of our hypotheses indicate that boys with ASD have an altered developmental pattern of the structures derived from the embryonic FNP and the MAND. It is well-documented that the developing FNP is derived from localized, specific cell populations under patterned genetic control. A number of developmental genes have been implicated in patterning the outgrowth of FNP, including FGF8, SHH and BMP2.
Neural crest cells that ultimately make up the FNP migrate over the forebrain to become encased within the neural ectoderm of the forebrain and the facial ectoderm . FGF8 plays a chemoattractive role in neural crest cell migration . Signals from the forebrain neuroectoderm are essential for the survival of the neural crest cells of the FNP , including SHH-dependent signaling from forebrain . Other studies have shown that SHH provides a key signal in regulating facial neural crest cell survival and patterning  and regulates BMP2 expression in the middle and upper face . Furthermore, it has been shown in an animal model that decreased SHH signaling leads to narrowing of the FNP and hypotelorism . The MAXs and LNPs are responsive to WNT signaling . WNT signaling leads to elevated cell proliferation and migration of neural crest cells . No midline structures are responsive to WNT.
Correct patterning and development of the forebrain also requires a balance of these genetic signaling factors , including SHH and the FGF, BMP and WNT families [20, 60]. SHH is expressed in the facial ectoderm and neuroectoderm at various developmental stages and acts synergistically with FGF8 in both face and brain development . SHH expression in facial ectoderm affects the expression of FGF8 in the brain . Likewise, expression of FGF8 by the forebrain is stimulated by the presence of neural crest cells .
It is clear that development of the face and brain is an interactive process, both anatomically and genetically. Gene expression studies have shown that facial and neural tube development are intimately interrelated. Altered gene expression patterns are associated with alterations in face and brain development. The reverse is also true: Altered face and brain development is associated with alterations in gene expression. We know that the brain is altered in people with ASD (reviewed in ), and the results of our present study show that the face is also affected. However, the sequence of events leading to these differences is unclear and may differ among the various subgroups described herein.
There is evidence derived from genetic studies implicating the developmental genes that control the patterning of the FNP and forebrain which causes autism. The SHH gene, though not identified by autism linkage or association studies, is functionally related to the Patched gene (PTCHD1), which is a strong autism candidate gene. Investigators who conducted CNV studies  first identified PTCHD1 gene microdeletions and missense mutations in males with ASD . SHH signaling repression is relieved when SHH binds to PTCHD1. BMP is part of the functional face and brain patterning network interacting with SHH and FGF8 to maintain brain and face patterning. Bakrania et al.  evaluated gene expression in embryos and demonstrated cotemporal and cospatial expression of BMP4 and SHH signaling genes. It is expected that sophisticated brain and/or face functional transcriptome studies may be useful in further linking genes involved in simultaneous face and brain development.
The clinical and behavioral differences that we identified between subgroup 1 and subgroup 2 boys support our hypothesis that the subgroups are biologically and etiologically distinctive. Subgroup 1 appears to be more severely autistic, with only 8% diagnosed with Asperger syndrome compared to 60% in subgroup 2, and has higher SCQ lifetime scores than subgroup 2 (25.3 vs 17.7). The percentage of IQ scores less than 70 were higher in subgroup 1 than in the other groups: 50% on the FSIQ and 45% on the NVIQ, compared to only 20% on both the FSIQ and NVIQ (one of five) in subgroup 2. and 15% FSIQ and 12% NVIQ in the ASD remainder group. Verbal IQ scores were also lower in subgroup 1 but did not reach statistical significance. Subgroup 1 also displayed several features predictive of poor outcome, including a higher risk for seizures and increased incidence of language regression at ASD onset. The observation of significantly higher SRS scores in mothers of boys in subgroup 1 is interesting but unexplained at this time.
By contrast, subgroup 2 appears to be aligned more with an Asperger syndrome diagnosis, which was made in 60% of subgroup 2 participants compared to only 8% in subgroup 1 and 32% in the remainder. Lifetime SCQ scores in subgroup 2 were also lower (17.7) than those in subgroup 1 (25.9) and the remainder (20.6). Consistent with the Asperger syndrome diagnosis, boys in subgroup 2 spoke their first words significantly earlier than boys in subgroup 1 and the remainder, and they were significantly more likely to be macrocephalic (80% vs 17%) compared to subgroup 1. Though IQ score differences did not reach significance within this small subgroup, boys in subgroup 2 had consistently higher FSIQ, NVIQ and VIQ scores than boys in subgroup 1.
Though we maximized our study outcomes by using a relatively homogeneous group of boys of Caucasian ancestry with essential autism within a narrow age range, the subject group was imperfect in a number of ways. Though 56% of the boys had participated in the SSC, not every participant completed the entire test battery, which slightly decreased the number of subjects who could be analyzed statistically for some comparisons. However, all of the four boys without VIQ scores and the one boy without a FSIQ score clustered morphologically within the main group and not within either subgroup. Thus our findings are highly unlikely to be affected by the small number of missing data points. In addition, the SSC population is biased toward a higher-functioning group of boys, which tended to shift the cognitive and outcome curves. The higher-functioning population recruited from the SSC may, however, have helped delineate subgroup 2 by increasing the proportion of subjects with Asperger syndrome. In addition, the number of boys in each subgroup was small. One additional observation is that in subgroup 2, which contained only five boys, one of the five was an outlier with significantly lower scores on IQ, language and outcome measures. This suggests that in future studies of larger numbers of subjects, we may find that subgroup 2 will be dissected into subgroups 2A and 2B. Finally, although one of the strengths of our study is that it comprised a homogeneous group of boys, that is, narrow age range, single sex, limited ethnic diversity and diagnosis made by a single clinician, it remains to be seen whether our findings will be consistent in a more heterogeneous population.
Differences in facial morphology may reflect alterations in embryologic brain development in children with ASD. Our results suggest potential differences in etiologies for the various subgroups of children. Further investigations into brain morphology will help to elucidate the causes and significance of these subtle differences. Verification of the role of a number of neurodevelopmental candidate genes may also be expedited by restricting analyses to the more homogeneous autism subgroups described herein. Likewise, based on our understanding of facial and neural development, identification of specific neurodevelopmental genes responsible for autism suggests which regions of the embryonic brain are most apt to be affected, providing potential target structures for future investigation.
Written informed consent was obtained from the parents of the patients for publication of this research article and any accompanying images. A copy of the written consent is available for review by the Editor-in-Chief of this journal.
This work was supported in part by Autism Spectrum Disorder Research Program grant W91ZSQ8006N604, the Simons Foundation and the Thompson Center Research Scholar Funds. The Departments of Pathology and Anatomical Sciences, Child Health, Radiology and Computer Engineering at the University of Missouri also provided financial support. Members of the Thompson Center staff obtained some of the 3dMD images. Written consent for publication was obtained from the patients or their relatives. Finally, we are extremely grateful to the families for their participation.
- American Psychiatric Association: Diagnostic and Statistical Manual of Mental Disorders, Fourth Edition (DSM-IV). 1994, Washington, DC: American Psychiatric AssociationGoogle Scholar
- Arndt TL, Stodgell CJ, Rodier PM: The teratology of autism. Int J Dev Neurosci. 2005, 23: 189-199. 10.1016/j.ijdevneu.2004.11.001.View ArticlePubMedGoogle Scholar
- O'Roak BJ, State MW: Autism genetics: strategies, challenges, and opportunities. Autism Res. 2008, 1: 4-17. 10.1002/aur.3.View ArticlePubMedGoogle Scholar
- Schaefer G, Mendelsohn N: Genetics evaluation for the etiologic diagnosis of autism spectrum disorders. Genet Med. 2008, 10: 4-12. 10.1097/GIM.0b013e31815efdd7.View ArticlePubMedGoogle Scholar
- Bauman M, Kemper T: Neuroanatomic observations of the brain in autism: a review and future directions. Int J Dev Neurosci. 2005, 23: 183-187. 10.1016/j.ijdevneu.2004.09.006.View ArticlePubMedGoogle Scholar
- Carroll SB, Grenier JK, Weatherbee SD: From DNA to Diversity: Molecular Genetics and the Evolution of Animal Design. 2001, Oxford: Blackwell PublishingGoogle Scholar
- Davidson EH: Genomic Regulatory Systems: Development and Evolution. 2001, San Diego: Academic PressGoogle Scholar
- Francis-West PH, Robson L, Evans DJR: Craniofacial Development: The Tissue and Molecular Interactions That Control Development of the Head. 2003, Berlin: SpringerView ArticleGoogle Scholar
- Jabs EW, Li X, Scott AF, Meyers G, Chen W, Eccles M, Mao JI, Charnas LR, Jackson CE, Jaye M: Jackson-Weiss and Crouzon syndromes are allelic with mutations in fibroblast growth factor receptor 2. Nat Genet. 1994, 8: 275-279. 10.1038/ng1194-275.View ArticlePubMedGoogle Scholar
- Klaauw C: Size and position of the functional components of the skull. Arch Neerl Zool. 1948-1952, 7: 16-37.View ArticleGoogle Scholar
- Mao JJ, Wang X, Mooney MP, Kopher RA, Nudera JA: Strain induced osteogenesis of the craniofacial suture upon controlled delivery of low-frequency cyclic forces. Front Biosci. 2003, 8: a10-a17. 10.2741/917.View ArticlePubMedGoogle Scholar
- Weiss KM: The phenogenetic logic of life. Nat Rev Genet. 2005, 6: 36-45. 10.1038/nrg1502.View ArticlePubMedGoogle Scholar
- Yu JC, Lucas JH, Fryberg K, Borke JL: Extrinsic tension results in FGF-2 release, membrane permeability change, and intracellular Ca++ increase in immature cranial sutures. J Craniofac Surg. 2001, 12: 391-398. 10.1097/00001665-200107000-00018.View ArticlePubMedGoogle Scholar
- Boughner JC, Wat S, Diewert VM, Young NM, Browder LW, Hallgrímsson B: Short-faced mice and developmental interactions between the brain and the face. J Anat. 2008, 213: 646-662. 10.1111/j.1469-7580.2008.00999.x.PubMed CentralView ArticlePubMedGoogle Scholar
- deMyer W, Zeman W, Palmer CG: The face predicts the brain: diagnostic significance of median facial anomalies for holoprosencephaly (arhinencephaly). Pediatrics. 1964, 34: 256-263.PubMedGoogle Scholar
- Diewert VM, Lozanoff S, Choy V: Computer reconstructions of human embryonic craniofacial morphology showing changes in relations between the face and brain during primary palate formation. J Craniofac Genet Dev Biol. 1993, 13: 193-201.PubMedGoogle Scholar
- Hennessy RJ, McLearie S, Kinsella A, Waddington JL: Facial shape and asymmetry by 3D laser surface scanning covary with cognition in a sexually dimorphic manner. J Neuropsychiatr Clin Neurosci. 2006, 18: 73-80. 10.1176/appi.neuropsych.18.1.73.View ArticleGoogle Scholar
- Foppiano S, Hu D, Marcucio RS: Signaling by bone morphogenetic proteins directs formation of an ectodermal signaling center that regulates craniofacial development. Dev Biol. 2007, 312: 103-114. 10.1016/j.ydbio.2007.09.016.PubMed CentralView ArticlePubMedGoogle Scholar
- Hu D, Marcucio RS, Helms JA: A zone of frontonasal ectoderm regulates patterning and growth in the face. Development. 2003, 130: 1749-1758. 10.1242/dev.00397.View ArticlePubMedGoogle Scholar
- Le Douarin NM, Brito JM, Creuzet S: Role of the neural crest in face and brain development. Brain Res Rev. 2007, 55: 237-247. 10.1016/j.brainresrev.2007.06.023.View ArticlePubMedGoogle Scholar
- Marcucio RS, Cordero DR, Hu D, Helms JA: Molecular interactions coordinating the development of the forebrain and face. Dev Biol. 2005, 284: 48-61. 10.1016/j.ydbio.2005.04.030.View ArticlePubMedGoogle Scholar
- Szabo-Rogers HL, Smithers LE, Yakob W, Liu KJ: New directions in craniofacial morphogenesis. Dev Biol. 2010, 341: 84-94. 10.1016/j.ydbio.2009.11.021.View ArticlePubMedGoogle Scholar
- Weiss KM, Buchanan AV: Genetics and the Logic of Evolution. 2004, Hoboken, NJ: Wiley-LissView ArticleGoogle Scholar
- Rodier PM: Environmental causes of central nervous system maldevelopment. Pediatrics. 2004, 113 (4 Suppl): 1076-1083.PubMedGoogle Scholar
- Rodier PM, Bryson SE, Welch JP: Minor malformations and physical measurements in autism: data from Nova Scotia. Teratology. 1997, 55: 319-325. 10.1002/(SICI)1096-9926(199705)55:5<319::AID-TERA4>3.0.CO;2-U.View ArticlePubMedGoogle Scholar
- Rodier PM: Converging evidence for brain stem injury in autism. Dev Psychopathol. 2002, 14: 537-557.View ArticlePubMedGoogle Scholar
- Rodier PM, Ingram JL, Tisdale B, Nelson S, Romano J: Embryological origin for autism: developmental anomalies of the cranial nerve motor nuclei. J Comp Neurol. 1996, 370: 247-261. 10.1002/(SICI)1096-9861(19960624)370:2<247::AID-CNE8>3.0.CO;2-2.View ArticlePubMedGoogle Scholar
- Kiecker C, Lumsden A: Compartments and their boundaries in vertebrate brain development. Nat Rev Neurosci. 2005, 6: 553-564. 10.1038/nrn1702.View ArticlePubMedGoogle Scholar
- Hammond P, Forster-Gibson C, Chudley AE, Allanson JE, Hutton TJ, Farrell SA, McKenzie J, Holden JJ, Lewis ME: Face-brain asymmetry in autism spectrum disorders. Mol Psychiatry. 2008, 13: 614-623. 10.1038/mp.2008.18.View ArticlePubMedGoogle Scholar
- Abzhanov A, Cordero DR, Sen J, Tabin CJ, Helms JA: Cross-regulatory interactions between Fgf8 and Shh in the avian frontonasal prominence. Congenit Anom (Kyoto). 2007, 47: 136-148. 10.1111/j.1741-4520.2007.00162.x.View ArticleGoogle Scholar
- Tapadia MD, Cordero DR, Helms JA: It's all in your head: new insights into craniofacial development and deformation. J Anat. 2005, 207: 461-477. 10.1111/j.1469-7580.2005.00484.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Brugmann SA, Goodnough LH, Gregorieff A, Leucht P, ten Berge D, Fuerer C, Clevers H, Nusse R, Helms JA: Wnt signaling mediates regional specification in the vertebrate face. Development. 2007, 134: 3283-3295. 10.1242/dev.005132.View ArticlePubMedGoogle Scholar
- Farkas LG, Posnick JC: Growth and development of regional units in the head and face based on anthropometric measurements. Cleft Palate Craniofac J. 1992, 29: 301-302. 10.1597/1545-1569(1992)029<0301:GADORU>2.3.CO;2.View ArticlePubMedGoogle Scholar
- Nellhaus G: Head circumference from birth to eighteen years: practical composite international and interracial graphs. Pediatrics. 1968, 41: 106-114.PubMedGoogle Scholar
- Dekaban AS: Changes in brain weights during the span of human life: relation of brain weights to body heights and body weights. Ann Neurol. 1978, 4: 345-356. 10.1002/ana.410040410.View ArticlePubMedGoogle Scholar
- Miroue M, Rosenberg L: The human facial sutures: a morphologic and histologic study of changes from 20 to 95 years. MSD thesis. 1975, University of WashingtonGoogle Scholar
- Lord C, Rutter M, Le Couteur A: Autism Diagnostic Interview-Revised: a revised version of a diagnostic interview for caregivers of individuals with possible pervasive developmental disorders. J Autism Dev Disord. 1994, 24: 659-685. 10.1007/BF02172145.View ArticlePubMedGoogle Scholar
- Lord C, Rutter M, Goode S, Heemsbergen J, Jordan H, Mawhood L, Schopler E: Autism diagnostic observation schedule: a standardized observation of communicative and social behavior. J Autism Dev Disord. 1989, 19: 185-212. 10.1007/BF02211841.View ArticlePubMedGoogle Scholar
- Miles JH, Takahashi TN, Hong J, Munden N, Flournoy N, Braddock SR, Martin RA, Bocian ME, Spence MA, Hillman RE, Farmer JE: Development and validation of a measure of dysmorphology: useful for autism subgroup classification. Am J Med Genet A. 2008, 146A: 1101-1116. 10.1002/ajmg.a.32244.View ArticlePubMedGoogle Scholar
- Farkas LG: Anthropometry of the head and face in clinical practice. Anthropometry of the Head and Face. Edited by: Farkas LG. 1994, New York: Raven Press, 71-77. 2Google Scholar
- Aldridge K, Boyadjiev SA, Capone GT, DeLeon VB, Richtsmeier JT: Precision and error of three-dimensional phenotypic measures acquired from 3dMD photogrammetric images. Am J Med Genet A. 2005, 138A: 247-253. 10.1002/ajmg.a.30959.View ArticlePubMedGoogle Scholar
- Weinberg SM, Scott NM, Neiswanger K, Brandon CA, Marazita ML: Digital three-dimensional photogrammetry: evaluation of anthropometric precision and accuracy using a Genex 3D camera system. Cleft Palate Craniofac J. 2004, 41: 507-518. 10.1597/03-066.1.View ArticlePubMedGoogle Scholar
- Lele S, Richtsmeier JT: An Invariant Approach to the Statistical Analysis of Shapes. 2001, London: Chapman & Hall/CRC PressView ArticleGoogle Scholar
- Lele S: Euclidean Distance Matrix Analysis (EDMA): Estimation of mean form and mean form difference. Math Geol. 1993, 25: 573-602. 10.1007/BF00890247.View ArticleGoogle Scholar
- Lele S, Richtsmeier JT: Euclidean Distance Matrix Analysis: confidence intervals for form and growth differences. Am J Phys Anthropol. 1995, 98: 73-86. 10.1002/ajpa.1330980107.View ArticlePubMedGoogle Scholar
- McIntyre GT, Mossey PA: Asymmetry of the craniofacial skeleton in the parents of children with a cleft lip, with or without a cleft palate, or an isolated cleft palate. Eur J Orthod. 2010, 32: 177-185. 10.1093/ejo/cjp067.View ArticlePubMedGoogle Scholar
- O'Grady KF, Antonyshyn OM: Facial asymmetry: three-dimensional analysis using laser surface scanning. Plast Reconstr Surg. 1999, 104: 928-937. 10.1097/00006534-199909040-00006.View ArticlePubMedGoogle Scholar
- Weinberg SM, Naidoo SD, Bardi KM, Brandon CA, Neiswanger K, Resick JM, Martin RA, Marazita ML: Face shape of unaffected parents with cleft affected offspring: combining three-dimensional surface imaging and geometric morphometrics. Orthod Craniofac Res. 2009, 12: 271-281. 10.1111/j.1601-6343.2009.01462.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Weinberg SM, Neiswanger K, Richtsmeier JT, Maher BS, Mooney MP, Siegel MI, Marazita ML: Three-dimensional morphometric analysis of craniofacial shape in the unaffected relatives of individuals with nonsyndromic orofacial clefts: a possible marker for genetic susceptibility. Am J Med Genet A. 2008, 146A: 409-420. 10.1002/ajmg.a.32177.View ArticlePubMedGoogle Scholar
- Richtsmeier JT, Cole TM, Krovitz G, Valeri CJ, Lele S: Preoperative morphology and development in sagittal synostosis. J Craniofac Genet Dev Biol. 1998, 18: 64-78.PubMedGoogle Scholar
- Darroch JN, Mosimann JE: Canonical and principal components of shape. Biometrika. 1985, 72: 241-252. 10.1093/biomet/72.2.241.View ArticleGoogle Scholar
- Falsetti AB, Jungers WL, Colle TM: Morphometrics of the callitrichid forelimb: a case study in size and shape. Int J Primat. 1993, 14: 551-572. 10.1007/BF02215447.View ArticleGoogle Scholar
- Jungers WL, Falsetti AB, Wall CE: Shape, relative size, and size-adjustments in morphometrics. Yearb Phys Anthropol. 1995, 38: 137-161. 10.1002/ajpa.1330380608.View ArticleGoogle Scholar
- Rutter M, Bailey A, Lord C: Social Communication Questionnaire. 2003, Los Angeles: Western Psychological ServicesGoogle Scholar
- Sparrow S, Balla D, Cicchetti D: Vineland Adaptive Behavior Scales. 1984, Circle Pines, MN: American Guidance ServiceGoogle Scholar
- Dunn L, Dunn L: Peabody Picture Vocabulary Test-Revised Manual. 1981, Circle Pines, MN: American Guidance ServiceGoogle Scholar
- Achenbach TM, Edelbrock C: Manual for the Child Behavior Checklist and Revised Behavior Profile. 1983, Burlington, VT: Department of Psychiatry, University of VermontGoogle Scholar
- Kanner L: Autistic disturbances of affective contact. Nerv Child. 1943, 2: 217-250.Google Scholar
- Young NM, Chong HJ, Hu D, Hallgrímsson B, Marcucio RS: Quantitative analyses link modulation of sonic hedgehog signaling to continuous variation in facial growth and shape. Development. 2010, 137: 3405-3409. 10.1242/dev.052340.PubMed CentralView ArticlePubMedGoogle Scholar
- Liu B, Rooker SM, Helms JA: Molecular control of facial morphology. Semin Cell Dev Biol. 2010, 21: 309-313. 10.1016/j.semcdb.2009.09.002.PubMed CentralView ArticlePubMedGoogle Scholar
- Creuzet S, Schuler B, Couly G, Le Douarin NM: Reciprocal relationships between Fgf8 and neural crest cells in facial and forebrain development. Proc Natl Acad Sci USA. 2004, 101: 4843-4847. 10.1073/pnas.0400869101.PubMed CentralView ArticlePubMedGoogle Scholar
- Verhoeven JS, De Cock P, Lagae L, Sunaert S: Neuroimaging of autism. Neuroradiology. 2010, 52: 3-14. 10.1007/s00234-009-0583-y.View ArticlePubMedGoogle Scholar
- Marshall CR, Noor A, Vincent JB, Lionel AC, Feuk L, Skaug J, Shago M, Moessner R, Pinto D, Ren Y, Thiruvahindrapduram B, Fiebig A, Schreiber S, Friedman J, Ketelaars CEJ, Vos YJ, Ficicioglu C, Kirkpatrick S, Nicolson R, Sloman L, Summers A, Gibbons CA, Teebi A, Chitayat D, Weksberg R, Thompson A, Vardy C, Crosbie V, Luscombe S, Baatjes R, Zwaigenbaum L, Roberts W, Fernandez B, Szatmari P, Scherer SW: Structural variation of chromosomes in autism spectrum disorder. Am J Hum Genet. 2008, 82: 477-488. 10.1016/j.ajhg.2007.12.009.PubMed CentralView ArticlePubMedGoogle Scholar
- Noor A, Whibley A, Marshall CR, Gianakopoulos PJ, Piton A, Carson AR, Orlic-Milacic M, Lionel AC, Sato D, Pinto D, Drmic I, Noakes C, Senman L, Zhang X, Mo R, Gauthier J, Crosbie J, Pagnamenta AT, Munson J, Estes AM, Fiebig A, Franke A, Schreiber S, Stewart AF, Roberts R, McPherson R, Guter SJ, Cook EH, Dawson G, Schellenberg GD, Battaglia A, Maestrini E, Autism Genome Project Consortium, Jeng L, Hutchison T, Rajcan-Separovic E, Chudley AE, Lewis SM, Liu X, Holden JJ, Fernandez B, Zwaigenbaum L, Bryson SE, Roberts W, Szatmari P, Gallagher L, Stratton MR, Gecz J, Brady AF, Schwartz CE, Schachar RJ, Monaco AP, Rouleau GA, Hui CC, Raymond FL, Scherer SW, Vincent JB: Disruption at the PTCHD1 locus on Xp22.11 in autism spectrum disorder and intellectual disability. Sci Transl Med. 2010, 2: 49ra68-10.1126/scitranslmed.3001267.PubMed CentralView ArticlePubMedGoogle Scholar
- Bakrania P, Efthymiou M, Klein JC, Salt A, Bunyan DJ, Wyatt A, Ponting CP, Martin A, Williams S, Lindley V, Gilmore J, Restori M, Robson AG, Neveu MM, Holder GE, Collin JR, Robinson DO, Farndon P, Johansen-Berg H, Gerrelli D, Ragge NK: Mutations in BMP4 cause eye, brain, and digit developmental anomalies: overlap between the BMP4 and hedgehog signaling pathways. Am J Hum Genet. 2008, 82: 304-319. 10.1016/j.ajhg.2007.09.023.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.