Vol 1, No 1 (2013) ISSN 2167-8677 (online) DOI 10.5195/d3000.2013.14 http://dentistry3000.pitt.edu This work is licensed under a Creative Commons Attribution 3.0 United States License. This site is published by the University Library System, University of Pittsburgh as part of its D-Scribe Digital Publishing Program and is cosponsored by the University of Pittsburgh Press. Heritability of face shape in twins: a preliminary study using 3D ster- ophotogrammetry and geometric morphometrics Seth M. Weinberg 1 , Trish E. Parsons 1 , Mary L. Marazita 1 , Brion S. Maher 2 1Center for Craniofacial and Dental Genetics, Department of Oral Biology, School of Dental Medicine, University of Pittsburgh, Pittsburgh, PA, USA 2Department of Mental Health, Bloomberg School of Public Health, Johns Hopkins University, Baltimore, MD, USA Abstract Introduction: Previous research suggests that aspects of facial surface morphology are her- itable. Traditionally, heritability studies have used a limited set of linear distances to quan- tify facial morphology and often employ statistical methods poorly designed to deal with biological shape. In this preliminary report, we use a combination of 3D photogrammetry and landmark-based morphometrics to explore which aspects of face shape show the strongest evidence of heritability in a sample of twins. Methods: 3D surface images were obtained from 21 twin pairs (10 monozygotic, 11 same- sex dizygotic). Thirteen 3D landmarks were collected from each facial surface and their co- ordinates subjected to geometric morphometric analysis. This involved superimposing the individual landmark configurations and then subjecting the resulting shape coordinates to a principal components analysis. The resulting PC scores were then used to calculate rough narrow-sense heritability estimates. Results: Three principal components displayed evidence of moderate to high heritability and were associated with variation in the breadth of orbital and nasal structures, upper lip height and projection, and the vertical and forward projection of the root of the nose due to variation in the position of nasion. Conclusions: Aspects of facial shape, primarily related to variation in length and breadth of central midfacial structures, were shown to demonstrate evidence of strong heritability. An improved understanding of which facial features are under strong genetic control is an im- portant step in the identification of specific genes that underlie normal facial variation. Citation: Weinberg SM, Parsons TE, Marazita ML, Ma- her BS. (2013) Heritability of face shape in twins: a pre- liminary study using 3D sterophotogrammerty and ge- ometric morphometrics. Dentistry 3000 1:a004 doi: 10.5195/d3000.2013.14 Received: August 19, 2013 Accepted: November 15, 2013 Published: November 25, 2013 Copyright: ©2013 Weinberg et al. This is an open- access article licensed under a Creative Commons Attribution Work 3.0 United States License. Email: smwst46@pitt.edu Introduction Variation in human facial morphology is strongly influenced by genetics. This con- tention is supported by decades of accumu- lated anthropometric research on twin con- cordance and parent-offspring resemblance [1-16]. To summarize the major findings of the previous research: 1) estimates of herit- ability for individual or composite quantita- tive facial features generally range from moderate (0.40) to high (0.80); 2) additive genetic effects account for much of the ob- served phenotypic variance in facial traits; 3) measured variables in each of the princi- pal anatomical directions (horizontal, verti- cal and projective) display strong evidence of heritability; 4) vertically oriented measures (e.g., facial height measures) tend to be among the most heritable facial traits. Determining the genetic contribution for specific facial traits remains a relevant area of concern in a number of fields that deal with human morphology, including paleoan- thropology and orthodontics. While the majority of craniofacial heritabil- ity studies are limited to measures of skele- tal structure, most often derived from radi- ographs, a smaller number of studies have focused solely on external soft-tissue fea- tures [17-35]. These studies typically use traditional anthropometric methods (i.e., direct measurement with calipers) to quan- tify the head and face, reporting similar heritability estimates for these surface traits as for the underlying skull. 3D surface im- aging represents an alternative, and far less invasive, method for capturing quantitative data on the facial soft tissues [36,37]. Re- cent advances in capture speed, accuracy, and affordability have rendered traditional methods of direct facial measurement large- ly obsolete. Moreover, by combining 3D surface imaging technology with powerful methods of landmark-based shape analysis, it is possible to quantify subtle aspects of facial morphology and variation often missed with conventional anthropometrics. Nevertheless, only two published studies to date have utilized 3D surface imaging to evaluate the genetic basis of facial features [38,39]. Although both studies found a very high degree of similarity among monozygot- ic twins compared with same-sex dizygotic pairs, with the central midfacial region ex- hibiting the greatest level of concordance, the quantitative analysis was limited to a small set of simple linear distances, which are often inadequate for describing complex 3D features. http://dentistry3000.pitt.edu/ http://creativecommons.org/licenses/by-nc-nd/3.0/us/ http://www.library.pitt.edu/ http://www.pitt.edu/ http://www.library.pitt.edu/articles/digpubtype/index.html http://upress.pitt.edu/ http://creativecommons.org/licenses/by-nc-nd/3.0/us/ Heritability of face shape in twins: a preliminary study using 3D sterophotogrammetry and geometric morphometrics Vol 1, No 1 (2013) DOI 10.5195/d3000.2013.14 http://dentistry3000.pitt.edu 2 In this preliminary investigation, we utilize 3D digital stereophotogrammetry in concert with geometric morphometrics to quantify midfacial structures in a sample of monozy- gotic (MZ) and same-sex dizygotic (DZ) twins. The main objective is to identify as- pects of midfacial shape that show the strongest evidence of heritability. Materials and Methods A set of 10 monozygotic twin pairs (20 indi- viduals) and 11 same-sex dizygotic twin pairs (22 individuals) were ascertained at the 2005 annual Twins Days Festival in Twinsburg, Ohio. The twin sample was comprised of Caucasian boys and girls be- tween the ages of 5 and 12 years (mean age = 9.3 years 2.2). Following informed con- sent, a 3D facial photo was obtained from each child as part of a broader assessment of phenotypic traits. Medical history was also collected from all participants to rule out any congenital conditions or trauma that might affect facial morphology. Rou- tine orthodontic interventions were not considered as a basis for exclusion. Parents of same-sex twins completed a standard 15- item questionnaire regarding the twins’ similarities and differences to determine zygosity [40]. The algorithm for determin- ing zygosity has been described elsewhere [41]. It has been previously demonstrated that this short questionnaire has an accura- cy of 94% for zygosity determination when compared with molecular genetic results [41,42]. All aspects of this study were ap- proved by an Institutional Review Board (IRB). Facial surface images were obtained via 3D digital stereophotogrammetry (FaceCam 250, Genex Technologies, Kensington, MD). Briefly, the imaging process involves pro- jecting a structured light pattern onto the facial surface, which is then recorded by multiple CCD cameras with overlapping fields of view. By extracting positional in- formation from the deformation of the pro- jected light pattern by the natural contours of the target surface, the location of points (> 300,000) distributed over the entire sur- face is determined. From this collection of points, an accurate 3D model of the facial surface is then generated (Figure 1). Color and texture information is captured simul- taneously and is mapped onto the underly- ing surface geometry, providing lifelike ren- dering. The entire image acquisition pro- cess takes place in less than 0.5 seconds. The accuracy and precision of measures derived these 3D surface models has been explored previously by the authors and found to very high [43,44]. Thirteen surface landmarks were collected from each 3D facial model (Figure 1) and their XYZ coordinate locations saved. These landmarks correspond primarily to central- ly located midfacial structures and are de- fined as in traditional surface anthropome- try [45]. There were several factors that limited the number and kinds of landmarks we were able to collect. For many subjects, the quality of the scans was compromised below the labial fissure, eliminating the opportunity to capture landmarks on the lower lips or chin. Furthermore, a number of subjects failed to maintain open eyes during the capture, making the identifica- tion of exocanthion impossible. A geometric morphometric approach was utilized to perform shape analysis of the landmark coordinate data [46,47]. First, the individu- al landmark configurations for all 42 partic- ipants were scaled and aligned through a least-squares Procrustes superimposition [48]. This procedure results in a new set of 3D coordinates (Procrustes coordinates), preserving only shape information. To con- trol for any remaining allometry effects in the data, centroid size was regressed against the Procrustes shape coordinates, and the resulting residuals were saved for further analysis. Principal components analysis (PCA) was then applied to the al- lometry-adjusted Procrustes coordinates in order to describe the shape variation pre- sent in the dataset. In the context of shape analysis, PCA reduces the dimensionality of the Procrustes coordinate data into a more manageable number of uncorrelated sum- mary variables (components), each captur- ing a distinct mode of shape variation [49]. Because the geometric morphometric ap- proach preserves the intrinsic geometry present in landmark coordinate data, shape variation associated with a given principal component (PC) can be visualized intuitive- ly as a displacement of points in 3D space [50,51]. To further enhance visualization, facial surface warps can be generated based on the PCA results. Figure 1. 3D facial surface showing the midfacial landmarks used in the present study. Midline points: n=nasion, prn=pronasale, sn=subnasale. Bilateral points: en=endocanthion, al=alare, sbal=subalare, cph=crista philtri, ch=chelion. http://dentistry3000.pitt.edu/ Heritability of face shape in twins: a preliminary study using 3D sterophotogrammetry and geometric morphometrics Vol 1, No 1 (2013) DOI 10.5195/d3000.2013.14 http://dentistry3000.pitt.edu 3 Similar to traditional PCA or factor analysis, every subject in the dataset receives a score on each extracted PC. These PC scores can then be used to compute narrow-sense her- itability (G) estimates with the following formula: G = 2(ICCMZ – ICCDZ) where ICC is the intraclass correlation coef- ficient. This is a very crude method for es- timating heritability, often used when am- ple sizes preclude the use of more refined methods, such as structural equation mod- eling. All geometric morphometric analyses (Procrustes analysis, allometry regression, and PCA) were carried out using the pro- gram MorphoJ [52]. Facial surface warps based on the PCA results were generated using the program Landmark v3.6 (http://www.idav.ucdavis.edu/research/Ev oMorph). Results A total of 17 PCs were extracted; ICCs were only calculated from scores on the first 9 PCs, which accounted for approximately 90% of the total shape variance. Three of the derived shape PCs displayed evidence of moderate to high heritability (PC4, PC5, and PC7); the remaining six PCs were unreveal- ing. The heritability estimates for PC4, PC5, and PC7 are provided in Table 1; although the values reported here approach 1.0, it should be noted these are still very crude approximations. We attempted to use boot- strap resampling techniques to obtain em- pirical confidence intervals for these esti- mates. These resulting confidence intervals approached or included zero and were like- ly affected by the small sample size from which the permutations were based. Indi- vidually, these PCs accounted for a relative- ly small amount of the total shape variance (PC4 = 6.8%, PC5 = 5.7%, PC7 = 4.0%). PC4 was associated with a complex suite of shape variations. In the upper face, varia- tion in the lateral position of the left/right endocanthion points (relating to the degree of hyper/hypotelorism) was observed. Var- iation in nasal breadth, height, and projec- tion was due primarily to shifts in the verti- cal and anterior-posterior position of nasi- on, the lateral position of the alar cartilages (left/right alare), and the anterior-posterior position of the nasal tip (pronasale). In the orolabial region, PC4 was associated with variation in the length of the philtrum and vertical height of the upper lip; this was coupled with relative changes in the vertical position of the left and right chelion. PC5 was associated mainly with the vertical and anterior-posterior position of nasion, the anterior-posterior position of the orbits (left/right endocanthion), and the breadth of the nasal floor and philtrum. Variation in shape associated with PC4 and PC5 are de- picted as surface warps in Figure 2. PC7 (not shown in Figure 2) was related primar- ily to the vertical position of landmarks defining the interorbital septum (nasion and the left/right endocanthion) and varia- tion in the nasolabial angle due to the ante- rior position of the upper lip relative to the nose. The first few PCs, accounting for the majority of shape variation in the sample, did not demonstrate strong evidence of heritability. Instead, shape variation along these PCs was related more to sex and age/maturity related factors. The source of the variation for the remaining PCs was not immediately apparent. Discussion The results of the present study revealed evidence of heritable shape variation in central midfacial structures, including the interorbital region, nose, and upper lip. Principal shape components associated with the degree of horizontal separation between the eyes, the length, breadth and projection of the nose, and the height and projection of the upper lip all demonstrated high herita- bility. These findings are largely supported by more traditional anthropometric studies that utilize standard linear distances to quantify facial surface structures [18-20,22- 35]. There is also substantial overlap be- tween the results of the present study and those of prior 3D facial surface studies on twin pairs. Burke [38] employed an early form of stereophotogrammetry to assess facial similarity in a sample of 18 like-sexed twin pairs. Of the 13 facial measures exam- ined, the degree of intrapair difference was significantly greater in DZ compared to MZ twins for eight variables, most of which involved central midfacial structures includ- ing the nose, interorbital region, and upper Table 1. Intraclass correlation coefficients (ICC) and G estimates for principal components showing evi- dence of high heritability. Shape Component Variation Explained ICCMZ ICCDZ Ga PC4 6.8% 0.783* 0.013 1.0b PC5 5.7% 0.771* 0.314 0.9 PC7 4.0% 0.675* 0.284 1.0b * p < 0.05 for test of within group non-zero ICC a G = 2(ICCMZ – ICCDZ) b Heritability estimate capped at 1.0 Figure 2. Facial surface warps modeling the shape variation associated with PC4 (along the horizontal axis) and PC5 (along the vertical axis). The centrally located face represents the average or consensus facial surface for the sample. Movement along the principal axes of shape variation in either direction away from the consensus face is accompanied by the 3D displacement of landmarks and modeled as a deformation of the consensus 3D surface. This deformation reflects the major modes of shape variation associated with PC4 and PC5, respectively. Figure 1. 3D facial surface showing the midfacial landmarks used in the present study. Midline points: n = nasion, prn = pronasale, sn = subnasale. Bilateral points: en = endocanthion, al = alare, sbal = subalare, cph = crista philtri, ch = chelion http://dentistry3000.pitt.edu/ http://www.idav.ucdavis.edu/research/EvoMorph http://www.idav.ucdavis.edu/research/EvoMorph Heritability of face shape in twins: a preliminary study using 3D sterophotogrammetry and geometric morphometrics Vol 1, No 1 (2013) DOI 10.5195/d3000.2013.14 http://dentistry3000.pitt.edu 4 lip. Using a similar study design, Naini and Moss [39] reported comparable results in a sample of 20 like-sexed twin pairs using 3D laser scanning technology. Kau et al. [53] utilized 3D laser surface scanning to evalu- ate facial resemblance in a single pair of monozygotic twins; within-pair differences at points spread over the facial surface were generally at the sub-millimeter level (the chin region showed the greatest disparity at 1.5-1.9mm). Several prior studies report the highest heritability estimates for vertically oriented facial measures [6,9-11,24,54], although there is considerable disagreement on this point [7,14,33-35]. The current findings present a slightly more complex picture. PC4 was associated with simultaneous changes in the relative length and breadth of the midface. This pattern of correlated shifts in length and breadth contributed to an overall shift in midface shape, from nar- row and long to short and broad (see Figure 2). The remaining two PCs showing evi- dence of heritability, however, were related primarily to the vertical and anterior posi- tion of the midline landmark nasion. Thus, there is only partial support for the conten- tion that vertical aspects of the face show higher heritability relative to variation along other principal anatomical directions. Determining the heritability of quantitative facial traits represents an important step in the effort to identify specific genes that un- derlie typical variation in face shape. To date, limited work in this area has been carried out on animal models. Quantitative trait loci (QTL) have been identified in mice that relate to variation in mandibular [55- 57] and maxillary shape [58,59]. In particu- lar, QTLs on several chromosomes have been associated with variation in midfacial breadth and length (projection) in the mouse skull [58,59]. In dogs, polymor- phisms in specific candidate genes (e.g., Tcof1, Fgf8, Runx2) have been shown to influence the extreme variation in facial form that characterizes certain breeds [60- 63]. In humans, a handful of recent studies have also begun to uncover the genetic basis of facial shape [64,65]. Interestingly, both of the aforementioned human studies have independently identified associations be- tween SNP variants in the PAX3 gene and variation in nasal root morphology, specifi- cally relating to the relative position of na- sion and left/right endocanthion. Variation in the spatial position of these same points was associated with all three relevant PCs in the present study, suggesting that our pre- liminary results, although quite limited, may have some external validity. It is important to keep in mind that alt- hough several derived shape PCs showed seemingly very high levels of heritability, these estimates are based on a simplistic approach to modeling heritability in a small sample and ultimately could not be substan- tiated using resampling methods. Conse- quently, the heritability values for PC4, PC5, and PC7 are likely to be overestimates of the true additive genetic effects on pheno- typic similarity. Rather than focusing on the values in an absolute sense, a more prudent interpretation might be that PC4, PC5, and PC7 demonstrated the most compelling evidence of heritability, relative to the other components of shape variation derived here. In future studies, larger samples will hopefully permit the use of more sophisti- cated methods for estimating heritability, including those based on structural equa- tion modeling. Conclusion This study represents a preliminary attempt to estimate the heritability of facial shape using a combination of 3D stereophoto- grammetry and landmark-based geometric morphometrics in a sample of twin pairs. Our results indicate that variation in both the length and breadth of central midfacial structures showed statistical evidence of moderate to high heritability. These results are in line with previous reports, including recent studies that have identified specific genetic variants underlying facial shape in humans. Nevertheless, due to the small size of the sample of twin pairs available for analysis, the results here should be viewed as preliminary and, as such, interpreted with caution. Acknowledgements The authors would like to thank Carla Bran- don and Kathy Bardi for their help recruit- ing subjects and Nicole Scott for helping to collect 3D images. We also thank the direc- tors of the Twins Days Festival for granting permission to conduct this research. This work was supported by the following grant from the National Institute of Dental and Craniofacial Research: R01-DE016148. Conflict of interest: There are no conflicts of interest to declare. References 1. Genetic determinants of cranio-facial morphol- ogy: a twin study. Nakata M, Yu PL, Davis B, Nance WE. Ann Hum Genet 1974;37:431-43. 2. Multivariate analysis of craniofacial measure- ments in twin and family data. Nakata M, Yu PL, Nance WE. Am J Phys Anthropol 1974;41:423-9. 3. A family study of craniofacial dimensions in the Burlington Growth Center sample. Saunders SR, Popovich F, Thompson GW. Am J Orthod 1980;78:394-403. 4. Sibling correlations for cranial measurements from serial radiographs. Byard PJ, Lewis AB, Ohtsuki F, Siervogel RM, Roche AF. J Craniofac Genet Dev Biol 1984;4:265-9. 5. Craniofacial morphology and occlusal variation in monozygotic and dizygotic twins. Lobb WK. Angle Orthodont 1987;57:219-33. 6. Comparison of some cephalometric distances and corresponding facial proportions with re- gard to heritability. Lundstrom A, McWilliam J. Eur J Orthod 1988;10:27-9. 7. Heritability of craniometric and occlusal varia- bles: a longitudinal sib analysis. Harris EF, John- son MG. Am J Orthod Dentofacial Orthop 1991;99:258-68. 8. Heritability of cephalometric and occlusal vari- ables as assessed from siblings with overt mal- occlusions. King L, Harris EF, Tolley EA. Am J Orthod Dentofacial Orthop 1993;104:121-31. 9. Heritability of 39 orthodontic cephalometric parameters on MZ, DZ twins and MN-paired singletons. Manfredi C, Martina R, Grossi GB, Giuliani M. Am J Orthod Dentofacial Orthop 1997;111:44-51. 10. A genetic study of anteroposterior and vertical facial proportions using model-fitting. Savoye I, Loos R, Carels C, Derom C, Vlietinck R. Angle Or- thodont 1998;68:467-70. 11. A quantitative genetic study of cephalometric variables in twins. Carels C, Van Cauwenberghe N, Savoye I, et al. Orthod Craniofacial Res 2001;4:130-40. 12. Heritability of craniofacial characteristics be- tween parents and offspring from lateral cepha- lograms. Johannsdottir B, Thorarinsson F, Thordarson A, Magnusson TE. Am J Orthod Dentofacial Orthop 2005;127:200-7. 13. Maximum likelihood estimation of human cra- niometric heritabilities. Carson EA. Am J Phys Anthropol 2006;131:169-80. 14. Heritability of facial proportions and soft-tissue profile characteristics in Turkish Anatolian sib- lings. Baydaş B, Erdem A, Yavuz I, Ceylan I. Am J Orthod Dentofacial Orthop 2007;131:504-9. 15. Quantitative genetics of modern human cranial variation. Sherwood RJ, Duren DL, Demerath EW, Czerwinski SA, Siervogel RM, Towne B. J Hum Evol 2008;54:909-14. 16. Heritability of human cranial dimensions: com- paring the evolvability of different cranial re- gions. Martínez-Abadías N, Esparza M, Sjøvold T, González-José R, Santos M, Hernández M. J Anat 2009;214:19-35. 17. Correlation of brothers in factor scores. Howells WW. Am J Phys Anthropol 1953;11:121-40. 18. The heritability of certain anthropometric char- acters as ascertained from measurements of twins. Clark PJ. Am J Hum Genet 1956;8:49-54. 19. Genetic Basis of Morphological Variation. Os- borne RH, De George FV. Cambridge: Harvard University Press; 1959. 20. A comparison of identical and fraternal twins on some anthropometric measures.Vandenberg SG, Strandskov HH. Hum Biol 1964;36:45-52. http://dentistry3000.pitt.edu/ Heritability of face shape in twins: a preliminary study using 3D sterophotogrammetry and geometric morphometrics Vol 1, No 1 (2013) DOI 10.5195/d3000.2013.14 http://dentistry3000.pitt.edu 5 21. Genetic and environmental influences on mor- phological characteristics. Susanne C. Ann Hum Biol 1975;2:279-87. 22. Heritability of anthropometric characters. Su- sanne C. Hum Biol 1977;49:573-80. 23. A family study of anthropometric traits in a Punjabi community: II. An investigation of fa- milial transmission. Byard PJ, Sharma K, Russell JM, Rao DC. Am J Phys Anthropol 1984;64:97- 104. 24. Family resemblance for components of cranio- facial size and shape. Byard PJ, Poosha DVR, Satyanarayana M, Rao DC. J Craniofac Genet Dev Biol 1985;5:229-38. 25. Family resemblance for cranio-facial measure- ments in Velanti Brahmins from Andhra Pra- desh, India. Poosha DV, Byard PJ, Satyanaraya- na M, Rice JP, Rao DC. Am J Phys Anthropol 1984;65:15-22. 26. A family study of anthropometric traits in a Punjabi community: I. Introduction and familial correlations. Sharma K, Byard PJ, Russell JM, Rao DC. Am J Phys Anthropol 1984;63:389-95. 27. Testing for the presence of genetic variance in factors of face measurements of Belgian twins. Hauspie RC, Susanne C, Defrise-Gussenhoven E. Ann Hum Biol 1985;12:429-40. 28. Transmissible and nontransmissible compo- nents of anthropometric variation in the Alex- anderwohl Mennonites: I. Description and fa- milial correlations. Devor EJ, McGue M, Craw- ford MH, Lin PM. Am J Phys Anthropol 1986;69:71-82. 29. Transmissible and nontransmissible compo- nents of anthropometric variation in the Alex- anderwohl Mennonites: II. Resolution by path analysis. Devor EJ, McGue M, Crawford MH, Lin PM. Am J Phys Anthropol 1986;69:83-92. 30. Transmission of human craniofacial dimensions. Devor EJ. J Craniofac Genet Dev Biol 1987;7:95- 106. 31. Genetics of human body size and shape: com- plex segregation analysis. Livshits G, Otremski I, Kobyliansky E. Ann Hum Biol 1995;22:13-27. 32. Similarities in anthropometrical traits of chil- dren and their parents in a Bulgarian popula- tion. Nikolova M. Ann Hum Genet 1996;60:517- 25. 33. Heritability of anthropometric phenotypes in caste populations of Visakhapatnam, India. Ar- ya R, Duggirala R, Comuzzie AG, Puppala S, Mo- dem S, Busi BR, Crawford MH. Hum Biol 2002;74:325-44. 34. Quantitative genetic study of head size related phenotypes in ethnically homogeneous Chu- vasha pedigrees. Ermakov S, Kobyliansky E, Livshits G. Ann Hum Biol 2005;32:585-98. 35. Genetic determination of head-size-related anthropometric traits in an ethnically homoge- neous sample of 373 Indian pedigrees of West Bengal. Karmakar B, Ermakov S, Yakovenko K, Kobyliansky E. Hum Biol 2007;79:501-14. 36. Three-dimensional surface imaging: limitations and considerations from the anthropometric perspective. Weinberg SM, Kolar JC. J Craniofac Surg 2005;16:847-51. 37. Measuring adult facial morphology in three dimensions. Kau CH, Richmond S, Savio C, Mal- lorie C. Angle Orthodont 2006;76:771-6. 38. Intrapair facial differences in twins. Burke PH. Acta Genet Med Gemellol (Roma) 1989;38:37- 47. 39. Three-dimensional assessment of the relative contribution of genetics and environment to various facial parameters with the twin meth- od. Naini FB, Moss JP. Am J Orthod Dentofacial Orthop 2004;126:655-65. 40. Diagnosis of twin zygosity. Nichols RC, Bilbro WC. Acta Genet Stat Med 1966;16:265-75. 41. Pittsburgh Registry of Infant Multiplets (PRIM). Strassberg M, Peters K, Marazita M, et al. Twin Res 2002;5:499-501. 42. Environmental and genetics influence on di- mensions of perceived parenting: a twin study. Rowe DC. Dev Psychol 1981;17:203-8. 43. Digital three-dimensional photogrammetry: evaluation of anthropometric precision and ac- curacy using a Genex 3D camera system. Wein- berg SM, Scott NM, Neiswanger K, Brandon CA, Marazita ML. Cleft Palate Craniofac J 2004;41:507-18. 44. Anthropometric precision and accuracy of digital three-dimensional photogrammetry: comparing the Genex and 3dMD imaging sys- tems to one another and to direct anthropome- try. Weinberg SM, Naidoo S, Govier DP, Martin RA, Kane AA, Marazita ML. J Craniofac Surg 2006;17:477-83. 45. Anthropometry of the Head and Face. Farkas LG. New York: Raven Press; 1994. 46. Statistical Shape Analysis. Dryden IL, Mardia KV. Chichester: Wiley; 1998. 47. In: Slice DE, ed. Modern Morphometrics in Physical Anthropology. Slice DE. Modern mor- phometrics. New York: Kluwer Academ- ic/Plenum Publishers; 2005:1-45. 48. Extensions of the Procrustes method for the optimal superimposition of landmarks. Rohlf F, Slice D. Syst Zool 1990;39:40-59. 49. Geometric Morphometrics for Biologists: A Primer. Zelditch ML, Swiderski DL, Sheets HD, Fink WL. Amsterdam: Elsevier Academic Press; 2004. 50. Facial growth in Cercocebus torquantus: an application of three-dimensional geometric morphometric techniques to the study of mor- phological variation. O'Higgins P, Jones N. J An- at 1998;193:251-72. 51. Geometric morphometric study of the regional variation of modern human craniofacial form. Hennessy RJ, Stringer CB. Am J Phys Anthropol 2002;117:37-48. 52. MorphoJ: an integrated software package for geometric morphometrics. Klingenberg CP. Mol Ecol Resour 2011;11:353-7. 53. The investigation of the changing facial appear- ance of identical twins employing a three- dimensional laser imaging system. Kau CH, Zhurov A, Bibb R, Hunter L, Richmond S. Orthod Craniofacial Res 2005;8:85-90. 54. A study of the inheritance of craniofacial char- acteristics as seen in lateral cephalograms of 72 like-sexed twins. Hunter WS. Rep Congr Eur Or- thod Soc 1965;41:59-70. 55. Pleiotropic effects of individual gene loci on mandibular morphology. Cheverud JM, Rout- man EJ, Irschick DJ. Evolution 1997;51:2006-16. 56. Quantitative trait loci for directional but not fluctuating asymmetry of mandibular charac- ters in mice. Leamy LJ, Pomp D, Eisen EJ, Cheve- rud JM. Genet Res 2000;76:27-40. 57. Genetic architecture of mandible shape in mice: effects of quantitative trait loci analyzed by ge- ometric morphometrics. Klingenberg CP, Leamy LJ, Routman EJ, Cheverud JM. Genetics 2001;157:785-802. 58. ENU large-scale mutagenesis and quantitative trait linkage (QTL) analysis in mice: novel tech- nologies for searching polygenetic determi- nants of craniofacial abnormalities. Nishimura I, Drake TA, Lusis AJ, Lyons KM, Nadeau JH, Zernik J. Crit Rev Oral Biol Med 2003;14:320-30. 59. A genome segment on mouse chromosome 12 determines maxillary growth. Oh J, Wang CJ, Poole M, et al. J Dent Res 2007;86:1203-6. 60. The canine homeobox gene MSX2: sequence, chromosome assignment and genetic analysis in dogs of different breeds. Haworth K, Breen M, Binns M, Hopkinson DA, Edwards YH. Anim Genet 2001;32:32-6. 61. Canine TCOF1; cloning, chromosome assign- ment and genetic analysis in dogs with different head types. Haworth KE, Islam I, Breen M, et al. Mamm Genome 2001;12:622-9. 62. Characterization of the genomic canine Fgf8 locus and screen for genetic variants in 4 dogs with different face types. Haworth KE, Healy C, McGonnell IM, Binns M, Sharpe PT. DNA Seq 2007;18:209-19. 63. Molecular origins of rapid and continuous mor- phological evolution. Fondon JW, Garner HR. Proc Natl Acad Sci USA 2004;101:18058-63. 64. Genome-wide association study of three- dimensional facial morphology identifies a vari- ant in PAX3 associated with nasion position. Paternoster L, Zhurov AI, Toma AM, Kemp JP, St Pourcain B, Timpson NJ, McMahon G, McArdle W, Ring SM, Smith GD, Richmond S, Evans DM. Am J Hum Genet 2012;90:478-85. 65. A genome-wide association study identifies five loci influencing facial morphology in Europe- ans.Liu F, van der Lijn F, Schurmann C, Zhu G, Chakravarty MM, Hysi PG, Wollstein A, Lao O, de Bruijne M, Ikram MA, van der Lugt A, Rivad- eneira F, Uitterlinden AG, Hofman A, Niessen WJ, Homuth G, de Zubicaray G, McMahon KL, Thompson PM, Daboul A, Puls R, Hegenscheid K, Bevan L, Pausova Z, Medland SE, Montgom- ery GW, Wright MJ, Wicking C, Boehringer S, Spector TD, Paus T, Martin NG, Biffar R, Kayser M. PLoS Genet 2012;8:e1002932. http://dentistry3000.pitt.edu/