Genome-Wide Association Study Reveals Multiple Loci Influencing Normal Human Facial Morphology

Numerous lines of evidence point to a genetic basis for facial morphology in humans, yet little is known about how specific genetic variants relate to the phenotypic expression of many common facial features. We conducted genome-wide association meta-analyses of 20 quantitative facial measurements derived from the 3D surface images of 3118 healthy individuals of European ancestry belonging to two US cohorts. Analyses were performed on just under one million genotyped SNPs (Illumina OmniExpress+Exome v1.2 array) imputed to the 1000 Genomes reference panel (Phase 3). We observed genome-wide significant associations (p < 5 x 10−8) for cranial base width at 14q21.1 and 20q12, intercanthal width at 1p13.3 and Xq13.2, nasal width at 20p11.22, nasal ala length at 14q11.2, and upper facial depth at 11q22.1. Several genes in the associated regions are known to play roles in craniofacial development or in syndromes affecting the face: MAFB, PAX9, MIPOL1, ALX3, HDAC8, and PAX1. We also tested genotype-phenotype associations reported in two previous genome-wide studies and found evidence of replication for nasal ala length and SNPs in CACNA2D3 and PRDM16. These results provide further evidence that common variants in regions harboring genes of known craniofacial function contribute to normal variation in human facial features. Improved understanding of the genes associated with facial morphology in healthy individuals can provide insights into the pathways and mechanisms controlling normal and abnormal facial morphogenesis.

[1]  Wiro J Niessen,et al.  Genetic determination of human facial morphology: links between cleft-lips and normal variation , 2011, European Journal of Human Genetics.

[2]  S. Takada,et al.  The AERO System: A 3D-Like Approach for Recording Gene Expression Patterns in the Whole Mouse Embryo , 2013, PloS one.

[3]  Peter Kraft,et al.  Quality control and quality assurance in genotypic data for genome‐wide association studies , 2010, Genetic epidemiology.

[4]  Seth M. Weinberg,et al.  The 3D Facial Norms Database: Part 1. A Web-Based Craniofacial Anthropometric and Image Repository for the Clinical and Research Community , 2016, The Cleft palate-craniofacial journal : official publication of the American Cleft Palate-Craniofacial Association.

[5]  T. Nakajima,et al.  Association of the growth hormone receptor gene polymorphisms with mandibular height in a Korean population. , 2009, Archives of oral biology.

[6]  C. Klingenberg,et al.  Genetic architecture of mandible shape in mice: effects of quantitative trait loci analyzed by geometric morphometrics. , 2001, Genetics.

[7]  J. Li,et al.  Adjusting multiple testing in multilocus analyses using the eigenvalues of a correlation matrix , 2005, Heredity.

[8]  Jason M. Johnson,et al.  A human homeotic transformation resulting from mutations in PLCB4 and GNAI3 causes auriculocondylar syndrome. , 2012, American journal of human genetics.

[9]  Manuel A. R. Ferreira,et al.  PLINK: a tool set for whole-genome association and population-based linkage analyses. , 2007, American journal of human genetics.

[10]  A. Graham,et al.  A population of caudally migrating cranial neural crest cells: functional and evolutionary implications. , 2001, Developmental biology.

[11]  C. Kokubu,et al.  Pax1 acts as a negative regulator of chondrocyte maturation. , 2013, Experimental cell research.

[12]  B. Papsin,et al.  Novel variants in GNAI3 associated with auriculocondylar syndrome strengthen a common dominant negative effect , 2015, European Journal of Human Genetics.

[13]  I. Krantz,et al.  Gorlin's Syndromes of the Head and Neck , 2010 .

[14]  S. O'Kane,et al.  Temporal and spatial expression of Pax9 and Sonic hedgehog during development of normal mouse palates and cleft palates in TGF-β3 null embryos , 2007 .

[15]  J. Korving,et al.  Severe nasal clefting and abnormal embryonic apoptosis in Alx3/Alx4 double mutant mice. , 2001, Development.

[16]  E. Cuppen,et al.  X-exome sequencing identifies a HDAC8 variant in a large pedigree with X-linked intellectual disability, truncal obesity, gynaecomastia, hypogonadism and unusual face , 2012, Journal of Medical Genetics.

[17]  M. Daly,et al.  Genome-wide association studies for common diseases and complex traits , 2005, Nature Reviews Genetics.

[18]  E. Kobyliansky,et al.  Complex segregation analysis of two principal components derived from horizontal and vertical head size traits , 2006, Annals of human biology.

[19]  M. Pericak-Vance,et al.  A Mutation in the TRPC6 Cation Channel Causes Familial Focal Segmental Glomerulosclerosis , 2005, Science.

[20]  Luisa F. Pallares,et al.  Mapping of Craniofacial Traits in Outbred Mice Identifies Major Developmental Genes Involved in Shape Determination , 2015, PLoS genetics.

[21]  Nicholas J Timpson,et al.  Genome-wide association study of three-dimensional facial morphology identifies a variant in PAX3 associated with nasion position. , 2012, American journal of human genetics.

[22]  D. Jimenez Craniofacial Anthropometry: Practical Measurement of the Head and Face for Clinical, Surgical and Research Use , 1998 .

[23]  S. Weinberg,et al.  Exploratory genotype–phenotype correlations of facial form and asymmetry in unaffected relatives of children with non‐syndromic cleft lip and/or palate , 2014, Journal of anatomy.

[24]  Sinem Bayram,et al.  Relationship between P561T and C422F polymorphisms in growth hormone receptor gene and mandibular prognathism. , 2014, The Angle orthodontist.

[25]  Washington Mio,et al.  Let's Face It—Complex Traits Are Just Not That Simple , 2014, PLoS genetics.

[26]  Michael L. Cunningham,et al.  Quantitative trait loci affecting the 3D skull shape and size in mouse and prioritization of candidate genes in-silico , 2015, Front. Physiol..

[27]  M. Bortolini,et al.  Fibroblast growth factor receptor 1 (FGFR1) variants and craniofacial variation in Amerindians and related populations , 2013, American journal of human biology : the official journal of the Human Biology Council.

[28]  Gabriele Gillessen-Kaesbach,et al.  HDAC8 mutations in Cornelia de Lange Syndrome affect the cohesin acetylation cycle , 2012, Nature.

[29]  E. Olson,et al.  Epigenetic control of skull morphogenesis by histone deacetylase 8. , 2009, Genes & development.

[30]  Ituro Inoue,et al.  Further evidence for an association between mandibular height and the growth hormone receptor gene in a Japanese population. , 2009, American journal of orthodontics and dentofacial orthopedics : official publication of the American Association of Orthodontists, its constituent societies, and the American Board of Orthodontics.

[31]  Family resemblance for components of craniofacial size and shape. , 1985, Journal of craniofacial genetics and developmental biology.

[32]  P. Donnelly,et al.  A Flexible and Accurate Genotype Imputation Method for the Next Generation of Genome-Wide Association Studies , 2009, PLoS genetics.

[33]  G. Morriss-Kay Craniofacial defects in AP-2 null mutant mice. , 1996, BioEssays : news and reviews in molecular, cellular and developmental biology.

[34]  Paul Suetens,et al.  Modeling 3D Facial Shape from DNA , 2014, PLoS genetics.

[35]  E. Zackai,et al.  Rearrangement in the PITX2 and MIPOL1 genes in a patient with a t(4;14) chromosome , 2003, European Journal of Human Genetics.

[36]  J. Marchini,et al.  Genotype Imputation with Thousands of Genomes , 2011, G3: Genes | Genomes | Genetics.

[37]  Marleen de Bruijne,et al.  A Genome-Wide Association Study Identifies Five Loci Influencing Facial Morphology in Europeans , 2012, PLoS genetics.

[38]  H B Wong,et al.  Genetics of cleft lip and cleft palate. , 1987, The Journal of the Singapore Paediatric Society.

[39]  T. Strom,et al.  Expanding the clinical spectrum of the ‘HDAC8‐phenotype’ – implications for molecular diagnostics, counseling and risk prediction , 2016, Clinical genetics.

[40]  Daniel Govier,et al.  The Craniofacial Phenotype of the Crouzon Mouse: Analysis of a Model for Syndromic Craniosynostosis Using Three-Dimensional MicroCT , 2006, The Cleft palate-craniofacial journal : official publication of the American Cleft Palate-Craniofacial Association.

[41]  L. Cooper Treatment of Nonsyndromic Anomalies of Tooth Number , 2015 .

[42]  A. van Daal,et al.  Linkage disequilibrium analysis identifies an FGFR1 haplotype-tag SNP associated with normal variation in craniofacial shape. , 2005, Genomics.

[43]  S. Weinberg,et al.  Face shape of unaffected parents with cleft affected offspring: combining three-dimensional surface imaging and geometric morphometrics. , 2009, Orthodontics & craniofacial research.

[44]  Holger Schwender,et al.  Erratum: A genome-wide association study of cleft lip with and without cleft palate identifies risk variants near MAFB and ABCA4 , 2010, Nature Genetics.

[45]  G. Willems,et al.  A quantitative genetic study of cephalometric variables in twins. , 2001, Clinical orthodontics and research.

[46]  Mauro Santos,et al.  Heritability of human cranial dimensions: comparing the evolvability of different cranial regions , 2009, Journal of anatomy.

[47]  Jay Shendure,et al.  Disruptive CHD8 Mutations Define a Subtype of Autism Early in Development , 2014, Cell.

[48]  Mary L Marazita,et al.  Three‐dimensional morphometric analysis of craniofacial shape in the unaffected relatives of individuals with nonsyndromic orofacial clefts: A possible marker for genetic susceptibility , 2008, American journal of medical genetics. Part A.

[49]  Y Shibasaki,et al.  Growth hormone receptor gene variant and mandibular height in the normal Japanese population. , 2001, American journal of orthodontics and dentofacial orthopedics : official publication of the American Association of Orthodontists, its constituent societies, and the American Board of Orthodontics.

[50]  Hang Zhou,et al.  Detecting Genetic Association of Common Human Facial Morphological Variation Using High Density 3D Image Registration , 2013, PLoS Comput. Biol..

[51]  Maja Šešelj,et al.  Heritability of the Human Craniofacial Complex , 2015, Anatomical record.

[52]  Christopher Lane,et al.  Completing the 3-dimensional picture. , 2008, American journal of orthodontics and dentofacial orthopedics : official publication of the American Association of Orthodontists, its constituent societies, and the American Board of Orthodontics.

[53]  J. Rogers,et al.  A Genomewide Linkage Scan for Quantitative Trait Loci Influencing the Craniofacial Complex in Baboons (Papio hamadryas spp.) , 2008, Genetics.

[54]  J. Helms,et al.  A primary cilia-dependent etiology for midline facial disorders. , 2010, Human molecular genetics.

[55]  L. Farkas,et al.  Anthropometry of the head and face, Second edition , 1994 .

[56]  M. Breen,et al.  The canine homeobox gene MSX2: sequence, chromosome assignment and genetic analysis in dogs of different breeds. , 2001, Animal genetics.

[57]  B. Hallgrímsson,et al.  The brachymorph mouse and the developmental‐genetic basis for canalization and morphological integration , 2006, Evolution & development.

[58]  Harry Hochheiser,et al.  The FaceBase Consortium: a comprehensive program to facilitate craniofacial research. , 2011, Developmental biology.

[59]  P. Scambler,et al.  Frontorhiny, a distinctive presentation of frontonasal dysplasia caused by recessive mutations in the ALX3 homeobox gene. , 2009, American journal of human genetics.

[60]  O. Delaneau,et al.  Supplementary Information for ‘ Improved whole chromosome phasing for disease and population genetic studies ’ , 2012 .

[61]  T. Yoneda,et al.  The transcription factor Znf219 regulates chondrocyte differentiation by assembling a transcription factory with Sox9 , 2010, Journal of Cell Science.

[62]  Kenny Q. Ye,et al.  An integrated map of genetic variation from 1,092 human genomes , 2012, Nature.

[63]  C. Cacou Anthropometry of the head and face , 1995 .

[64]  Holger Schwender,et al.  A genome-wide association study of cleft lip with and without cleft palate identifies risk variants near MAFB and ABCA4 , 2010, Nature Genetics.

[65]  Nathan M. Young,et al.  Quantitative analyses link modulation of sonic hedgehog signaling to continuous variation in facial growth and shape , 2010, Development.

[66]  J. T. Wright,et al.  Craniofacial and Dental Developmental Defects , 2015 .

[67]  Sergey Ermakov,et al.  Family-based study of association between ENPP1 genetic variants and craniofacial morphology , 2010, Annals of human biology.

[68]  Martin S. Taylor,et al.  Heterozygous loss-of-function mutations in YAP1 cause both isolated and syndromic optic fissure closure defects. , 2014, American journal of human genetics.

[69]  H. Hayasaki,et al.  The P561T polymorphism of the growth hormone receptor gene has an inhibitory effect on mandibular growth in young children. , 2009, European journal of orthodontics.

[70]  K. Kratochwil,et al.  Pax9-deficient mice lack pharyngeal pouch derivatives and teeth and exhibit craniofacial and limb abnormalities. , 1998, Genes & development.