Genome-wide study of gene-by-sex interactions identifies risks for cleft palate.

Structural birth defects affect 3-4% of all live births and, depending on the type, tend to manifest in a sex-biased manner. Orofacial clefts (OFCs) are the most common craniofacial structural birth defects and are often divided into cleft lip with or without cleft palate (CL/P) and cleft palate only (CP). Previous studies have found sex-specific risks for CL/P, but these risks have yet to be evaluated in CP. CL/P is more common in males and CP is more frequently observed in females, so we hypothesized there would also be sex-specific differences for CP. Using a trio-based cohort, we performed sex-stratified genome-wide association studies (GWAS) based on proband sex followed by a genome-wide gene-by-sex (G × S) interaction testing. There were 13 loci significant for G × S interactions, with the top finding in LTBP1 (RR = 3.37 [2.04-5.56], p = 1.93 × 10-6). LTBP1 plays a role in regulating TGF-ß bioavailability, and knockdown in both mice and zebrafish lead to craniofacial anomalies. Further, there is evidence for differential expression of LTBP1 between males and females in both mice and humans. Therefore, we tested the association between the imputed genetically regulated gene expression of genes with significant G × S interactions and the CP phenotype. We found significant association for LTBP1 in cell cultured fibroblasts in female probands (p = 0.0013) but not in males. Taken altogether, we show there are sex-specific risks for CP that are otherwise undetectable in a combined sex cohort, and LTBP1 is a candidate risk gene, particularly in females.

Genome scan reveals several loci associated with torus palatinus.

OBJECTIVE: Torus Palatinus (TP) is a common trait with an unclear aetiology. Although prior studies suggest a hereditary component, the genetic factors that influence TP risk remain unknown. The purpose of this study is to identify genetic variants associated with TP. MATERIALS AND METHODS: We assessed the TP status of 829 individuals from various ancestral backgrounds using 3D palate scans. We then carried out a genome-wide association study (GWAS) to identify common variants associated with TP. We also performed gene-based tests across the exome to investigate the role of low-frequency coding variants. RESULTS: Our GWAS did not identify any genome-wide significant signals but identified suggestive associations including hits on chromosomes 2, 5 and 17 with p-values less than 5 × 10-6. Candidate genes at these suggestive loci have been implicated in normal-range craniofacial features, syndromes with facial and oral malformations, and bone density. We did not find evidence that low-frequency coding variants influence TP risk. In addition, we failed to replicate associations identified in prior genetic studies of TP. CONCLUSION: These findings suggest that multiple genes likely influence the development of TP. Independent replication will be required to confirm our suggestive associations.

DNA Methylation Effects on Van der Woude Syndrome Phenotypic Variability.

OBJECTIVE: Van der Woude Syndrome (VWS) presents with combinations of lip pits (LP) and cleft lip and/or cleft palate (CL/P, CPO). VWS phenotypic heterogeneity even amongst relatives, suggests that epigenetic factors may act as modifiers. IRF6, causal for 70% of VWS cases, and TP63 interact in a regulatory loop coordinating epithelial proliferation and differentiation in palatogenesis. We hypothesize that differential DNA methylation within IRF6 and TP63 regulatory regions underlie VWS phenotypic discordance. METHODS: DNA methylation of CpG sites in IRF6 and TP63 promoters and in an IRF6 enhancer element was compared amongst blood or saliva DNA samples of 78 unrelated cases. Analyses were done separately for blood and saliva, within each sex and in combination, and to address cleft type (CL/P ± LP vs. CPO ± LP) and phenotypic severity (any cleft + LP vs. any cleft only). RESULTS: For cleft type, blood samples showed higher IRF6 and TP63 promoter methylation on males with CPO ± LP compared to CL/P ± LP and on individuals with CPO ± LP compared to those with CL/P ± LP, respectively. Saliva samples showed higher IRF6 enhancer methylation on individuals with CPO ± LP compared to CL/P ± LP and contrary to above, lower TP63 promoter methylation on CPO ± LP compared to CL/P ± LP. For phenotypic severity, blood samples showed no differences; however, saliva samples showed higher IRF6 promoter methylation in individuals with any cleft + LP compared to those without lip pits. CONCLUSION: We observed differential methylation in IRF6 and TP63 regulatory regions associated with cleft type and phenotypic severity, indicating that epigenetic changes in IRF6 and TP63 can contribute to phenotypic heterogeneity in VWS.

Genome-wide association study of multiethnic nonsyndromic orofacial cleft families identifies novel loci specific to family and phenotypic subtypes.

Nonsyndromic orofacial clefts (OFCs) are among the most common craniofacial birth defects worldwide, and known to exhibit phenotypic and genetic heterogeneity. Cleft lip plus cleft palate (CLP) and cleft lip only (CL) are commonly combined together as one phenotype (CL/P), separately from cleft palate alone. In comparison, our study analyzes CL and CLP separately. A sample of 2218 CL and CLP cases, 4537 unaffected relatives of cases, and 2673 pure controls with no family history of OFC were selected from the Pittsburgh Orofacial Cleft (Pitt-OFC) multiethnic study.genome-wide association studies were run for seven specific phenotypes created based on the cleft type(s) observed within these families, as well as the combined CL/P phenotype. Five novel genome-wide significant associations, 3q29 (rs62284390), 5p13.2 (rs609659), 7q22.1 (rs6465810), 19p13.3 (rs628271), and 20q13.33 (rs2427238), and nine associations (p ≤ 1.0E-05) within previously confirmed OFC loci-PAX7, IRF6, FAM49A, DCAF4L2, 8q24.21, ARID3B, NTN1, TANC2 and the WNT9B:WNT3 gene cluster-were observed. We also found that single nucleotide polymorphisms within a subset of the associated loci, both previously known and novel, differ substantially in terms of their effects across cleft- or family-specific phenotypes, indicating not only etiologic differences between CL and CLP, but also genetic heterogeneity within each of the two OFC subtypes.

Genome-wide Enrichment of De Novo Coding Mutations in Orofacial Cleft Trios.

Although de novo mutations (DNMs) are known to increase an individual's risk of congenital defects, DNMs have not been fully explored regarding orofacial clefts (OFCs), one of the most common human birth defects. Therefore, whole-genome sequencing of 756 child-parent trios of European, Colombian, and Taiwanese ancestry was performed to determine the contributions of coding DNMs to an individual's OFC risk. Overall, we identified a significant excess of loss-of-function DNMs in genes highly expressed in craniofacial tissues, as well as genes associated with known autosomal dominant OFC syndromes. This analysis also revealed roles for zinc-finger homeobox domain and SOX2-interacting genes in OFC etiology.

Genetic Analyses of Enamel Hypoplasia in Multiethnic Cohorts.

Enamel hypoplasia causes reduction in the thickness of affected enamel and is one of the most common dental anomalies. This defect is caused by environmental and/or genetic factors that interfere with tooth formation, emphasizing the importance of investigating enamel hypoplasia on an epidemiological and genetic level. A genome-wide association of enamel hypoplasia was performed in multiple cohorts, overall comprising 7,159 individuals ranging in age from 7-82 years. Mixed-models were used to test for genetic association while simultaneously accounting for relatedness and genetic population structure. Meta-analysis was then performed. More than 5 million single-nucleotide polymorphisms were tested in individual cohorts. Analyses of the individual cohorts and meta-analysis identified association signals close to genome-wide significance (P < 510-8), and many suggestive association signals (510-8 < P < 510-6) near genes with plausible roles in tooth/enamel development. The strongest association signal (P = 1.5710-9) was observed near BMP2K in one of the individual cohorts. Additional suggestive signals were observed near genes with plausible roles in tooth development in the meta-analysis, such as SLC4A4 which can influence enamel hypoplasia. Additional human genetic studies are needed to replicate these results and functional studies in model systems are needed to validate our findings.

Mutations in the Epithelial Cadherin-p120-Catenin Complex Cause Mendelian Non-Syndromic Cleft Lip with or without Cleft Palate.

Non-syndromic cleft lip with or without cleft palate (NS-CL/P) is one of the most common human birth defects and is generally considered a complex trait. Despite numerous loci identified by genome-wide association studies, the effect sizes of common variants are relatively small, with much of the presumed genetic contribution remaining elusive. We report exome-sequencing results in 209 people from 72 multi-affected families with pedigree structures consistent with autosomal-dominant inheritance and variable penetrance. Herein, pathogenic variants are described in four genes encoding components of the p120-catenin complex (CTNND1, PLEKHA7, PLEKHA5) and an epithelial splicing regulator (ESRP2), in addition to the known CL/P-associated gene, CDH1, which encodes E-cadherin. The findings were also validated in a second cohort of 497 people with NS-CL/P, comprising small families and singletons with pathogenic variants in these genes identified in 14% of multi-affected families and 2% of the replication cohort of smaller families. Enriched expression of each gene/protein in human and mouse embryonic oro-palatal epithelia, demonstration of functional impact of CTNND1 and ESRP2 variants, and recapitulation of the CL/P spectrum in Ctnnd1 knockout mice support a causative role in CL/P pathogenesis. These data show that primary defects in regulators of epithelial cell adhesion are the most significant contributors to NS-CL/P identified to date and that inherited and de novo single gene variants explain a substantial proportion of NS-CL/P.

FoxO6 regulates Hippo signaling and growth of the craniofacial complex.

The mechanisms that regulate post-natal growth of the craniofacial complex and that ultimately determine the size and shape of our faces are not well understood. Hippo signaling is a general mechanism to control tissue growth and organ size, and although it is known that Hippo signaling functions in neural crest specification and patterning during embryogenesis and before birth, its specific role in postnatal craniofacial growth remains elusive. We have identified the transcription factor FoxO6 as an activator of Hippo signaling regulating neonatal growth of the face. During late stages of mouse development, FoxO6 is expressed specifically in craniofacial tissues and FoxO6-/- mice undergo expansion of the face, frontal cortex, olfactory component and skull. Enlargement of the mandible and maxilla and lengthening of the incisors in FoxO6-/- mice are associated with increases in cell proliferation. In vitro and in vivo studies demonstrated that FoxO6 activates Lats1 expression, thereby increasing Yap phosphorylation and activation of Hippo signaling. FoxO6-/- mice have significantly reduced Hippo Signaling caused by a decrease in Lats1 expression and decreases in Shh and Runx2 expression, suggesting that Shh and Runx2 are also linked to Hippo signaling. In vitro, FoxO6 activates Hippo reporter constructs and regulates cell proliferation. Furthermore PITX2, a regulator of Hippo signaling is associated with Axenfeld-Rieger Syndrome causing a flattened midface and we show that PITX2 activates FoxO6 expression. Craniofacial specific expression of FoxO6 postnatally regulates Hippo signaling and cell proliferation. Together, these results identify a FoxO6-Hippo regulatory pathway that controls skull growth, odontogenesis and face morphology.

Genome-Wide Association Study (GWAS) of dental caries in diverse populations.

BACKGROUND: Dental caries is one of the most common chronic diseases and is influenced by a complex interplay of genetic and environmental factors. Most previous genetic studies of caries have focused on identifying genes that contribute to dental caries in specific ethnic groups, usually of European descent. METHODS: The aim of this study is to conduct a genome-wide association study (GWAS) to identify associations affecting susceptibility to caries in a large multiethnic population from Argentina, the Philippines, Guatemala, Hungary, and the USA, originally recruited for studies of orofacial clefts (POFC, N = 3686). Ages of the participants ranged from 2 to 12 years for analysis of the primary dentition, and 18-60 years for analysis of the permanent dentition. For each participant, dental caries was assessed by counts of decayed and filled teeth (dft/DFT) and genetic variants (single nucleotide polymorphisms, SNPs) were genotyped or imputed across the entire genome. Caries was analyzed separately for the primary and permanent dentitions, with age, gender, and presence/absence of any type of OFC treated as covariates. Efficient Mixed-Model Association eXpedited (EMMAX) was used to test genetic association, while simultaneously accounting for relatedness and stratification. RESULTS: We identified several suggestive loci (5 × 10-8 < P < 5 × 10-6) within or near genes with plausible biological roles for dental caries, including a cluster of taste receptor genes (TAS2R38, TAS2R3, TAS2R4, TASR25) on chromosome 7 for the permanent dentition analysis, and DLX3 and DLX4 on chromosome 17 for the primary dentition analysis. Genome-wide significant results were seen with SNPs in the primary dentition only; however, none of the identified genes near these variants have known roles in cariogenesis. CONCLUSION: The results of this study warrant further investigation and may lead to a better understanding of cariogenesis in diverse populations, and help to improve dental caries prediction, prevention, and/or treatment in future.

Mutations in GDF11 and the extracellular antagonist, Follistatin, as a likely cause of Mendelian forms of orofacial clefting in humans.

Cleft lip with or without cleft palate (CL/P) is generally viewed as a complex trait with multiple genetic and environmental contributions. In 70% of cases, CL/P presents as an isolated feature and/or deemed nonsyndromic. In the remaining 30%, CL/P is associated with multisystem phenotypes or clinically recognizable syndromes, many with a monogenic basis. Here we report the identification, via exome sequencing, of likely pathogenic variants in two genes that encode interacting proteins previously only linked to orofacial clefting in mouse models. A variant in GDF11 (encoding growth differentiation factor 11), predicting a p.(Arg298Gln) substitution at the Furin protease cleavage site, was identified in one family that segregated with CL/P and both rib and vertebral hypersegmentation, mirroring that seen in Gdf11 knockout mice. In the second family in which CL/P was the only phenotype, a mutation in FST (encoding the GDF11 antagonist, Follistatin) was identified that is predicted to result in a p.(Cys56Tyr) substitution in the region that binds GDF11. Functional assays demonstrated a significant impact of the specific mutated amino acids on FST and GDF11 function and, together with embryonic expression data, provide strong evidence for the importance of GDF11 and Follistatin in the regulation of human orofacial development.

Genome-wide interaction studies identify sex-specific risk alleles for nonsyndromic orofacial clefts.

Nonsyndromic cleft lip with or without cleft palate (NSCL/P) is the most common craniofacial birth defect in humans and is notable for its apparent sexual dimorphism where approximately twice as many males are affected as females. The sources of this disparity are largely unknown, but interactions between genetic and sex effects are likely contributors. We examined gene-by-sex (G × S) interactions in a worldwide sample of 2,142 NSCL/P cases and 1,700 controls recruited from 13 countries. First, we performed genome-wide joint tests of the genetic (G) and G × S effects genome-wide using logistic regression assuming an additive genetic model and adjusting for 18 principal components of ancestry. We further interrogated loci with suggestive results from the joint test ( p < 1.00 × 10 -5 ) by examining the G × S effects from the same model. Out of the 133 loci with suggestive results ( p < 1.00 × 10 -5 ) for the joint test, we observed one genome-wide significant G × S effect in the 10q21 locus (rs72804706; p = 6.69 × 10 -9 ; OR = 2.62 CI [1.89, 3.62]) and 16 suggestive G × S effects. At the intergenic 10q21 locus, the risk of NSCL/P is estimated to increase with additional copies of the minor allele for females, but the opposite effect for males. Our observation that the impact of genetic variants on NSCL/P risk differs for males and females may further our understanding of the genetic architecture of NSCL/P and the sex differences underlying clefts and other birth defects.

Identification of 16q21 as a modifier of nonsyndromic orofacial cleft phenotypes.

Orofacial clefts (OFCs) are common, complex birth defects with extremely heterogeneous phenotypic presentations. Two common subtypes-cleft lip alone (CL) and CL plus cleft palate (CLP)-are typically grouped into a single phenotype for genetic analysis (i.e., CL with or without cleft palate, CL/P). However, mounting evidence suggests there may be unique underlying pathophysiology and/or genetic modifiers influencing expression of these two phenotypes. To this end, we performed a genome-wide scan for genetic modifiers by directly comparing 450 CL cases with 1,692 CLP cases from 18 recruitment sites across 13 countries from North America, Central or South America, Asia, Europe, and Africa. We identified a region on 16q21 that is strongly associated with different cleft type (P = 5.611 × 10-8 ). We also identified significant evidence of gene-gene interactions between this modifier locus and two recognized CL/P risk loci: 8q21 and 9q22 (FOXE1) (P = 0.012 and 0.023, respectively). Single nucleotide polymorphism (SNPs) in the 16q21 modifier locus demonstrated significant association with CL over CLP. The marker alleles on 16q21 that increased risk for CL were found at highest frequencies among individuals with a family history of CL (P = 0.003). Our results demonstrate the existence of modifiers for which type of OFC develops and suggest plausible elements responsible for phenotypic heterogeneity, further elucidating the complex genetic architecture of OFCs.

Soft tissue nasal asymmetry as an indicator of orofacial cleft predisposition.

The biological relatives of offspring with nonsyndromic orofacial clefts have been shown to exhibit distinctive facial features, including excess asymmetry, which are hypothesized to indicate the presence of genetic risk factors. The significance of excess soft tissue nasal asymmetry in at-risk relatives is unclear and was examined in the present study. Our sample included 164 unaffected parents from families with a history of orofacial clefting and 243 adult controls. Geometric morphometric methods were used to analyze the coordinates of 15 nasal landmarks collected from three-dimensional facial surface images. Following generalized Procrustes analysis, Procrustes ANOVA and MANOVA tests were applied to determine the type and magnitude of nasal asymmetry present in each group. Group differences in mean nasal asymmetry were also assessed via permutation testing. We found that nasal asymmetry in both parents and controls was directional in nature, although the magnitude of the asymmetry was greater in parents. This was confirmed with permutation testing, where the mean nasal asymmetry was significantly different (p < .0001) between parents and controls. The asymmetry was greatest for midline structures and the nostrils. When subsets of parents were subsequently analyzed and compared (parents with bilateral vs. unilateral offspring; parents with left vs. right unilateral offspring), each group showed a similar pattern of asymmetry and could not be distinguished statistically. Thus, the side of the unilateral cleft (right vs. left) in offspring was not associated with the direction of the nasal asymmetry in parents.

Effects of Objective 3-Dimensional Measures of Facial Shape and Symmetry on Perceptions of Facial Attractiveness.

PURPOSE: Meeting patient desires for enhanced facial esthetics requires that providers have standardized and objective methods to measure esthetics. The authors evaluated the effects of objective 3-dimensional (3D) facial shape and asymmetry measurements derived from 3D facial images on perceptions of facial attractiveness. MATERIALS AND METHODS: The 3D facial images of 313 adults in Iowa were digitized with 32 landmarks, and objective 3D facial measurements capturing symmetric and asymmetric components of shape variation, centroid size, and fluctuating asymmetry were obtained from the 3D coordinate data using geo-morphometric analyses. Frontal and profile images of study participants were rated for facial attractiveness by 10 volunteers (5 women and 5 men) on a 5-point Likert scale and a visual analog scale. Multivariate regression was used to identify the effects of the objective 3D facial measurements on attractiveness ratings. RESULTS: Several objective 3D facial measurements had marked effects on attractiveness ratings. Shorter facial heights with protrusive chins, midface retrusion, faces with protrusive noses and thin lips, flat mandibular planes with deep labiomental folds, any cants of the lip commissures and floor of the nose, larger faces overall, and increased fluctuating asymmetry were rated as significantly (P < .001) less attractive. CONCLUSION: Perceptions of facial attractiveness can be explained by specific 3D measurements of facial shapes and fluctuating asymmetry, which have important implications for clinical practice and research.

Candidate gene analyses of 3-dimensional dentoalveolar phenotypes in subjects with malocclusion.

INTRODUCTION: Genetic studies of malocclusion etiology have identified 4 deleterious mutations in genes DUSP6,ARHGAP21, FGF23, and ADAMTS1 in familial Class III cases. Although these variants may have large impacts on Class III phenotypic expression, their low frequency (<1%) makes them unlikely to explain most malocclusions. Thus, much of the genetic variation underlying the dentofacial phenotypic variation associated with malocclusion remains unknown. In this study, we evaluated associations between common genetic variations in craniofacial candidate genes and 3-dimensional dentoalveolar phenotypes in patients with malocclusion. METHODS: Pretreatment dental casts or cone-beam computed tomographic images from 300 healthy subjects were digitized with 48 landmarks. The 3-dimensional coordinate data were submitted to a geometric morphometric approach along with principal component analysis to generate continuous phenotypes including symmetric and asymmetric components of dentoalveolar shape variation, fluctuating asymmetry, and size. The subjects were genotyped for 222 single-nucleotide polymorphisms in 82 genes/loci, and phenotpye-genotype associations were tested via multivariate linear regression. RESULTS: Principal component analysis of symmetric variation identified 4 components that explained 68% of the total variance and depicted anteroposterior, vertical, and transverse dentoalveolar discrepancies. Suggestive associations (P < 0.05) were identified with PITX2, SNAI3, 11q22.2-q22.3, 4p16.1, ISL1, and FGF8. Principal component analysis for asymmetric variations identified 4 components that explained 51% of the total variations and captured left-to-right discrepancies resulting in midline deviations, unilateral crossbites, and ectopic eruptions. Suggestive associations were found with TBX1AJUBA, SNAI3SATB2, TP63, and 1p22.1. Fluctuating asymmetry was associated with BMP3 and LATS1. Associations for SATB2 and BMP3 with asymmetric variations remained significant after the Bonferroni correction (P <0.00022). Suggestive associations were found for centroid size, a proxy for dentoalveolar size variation with 4p16.1 and SNAI1. CONCLUSIONS: Specific genetic pathways associated with 3-dimensional dentoalveolar phenotypic variation in malocclusions were identified.

Maternal alcohol binge-drinking in the first trimester and the risk of orofacial clefts in offspring: a large population-based pooling study.

Using individual participant data from six population-based case-control studies, we conducted pooled analyses to examine maternal alcohol consumption and the risk of clefts among >4600 infants with cleft lip only, cleft lip with cleft palate, or cleft palate only and >10,000 unaffected controls. We examined two first-trimester alcohol measures: average number of drinks/sitting and maximum number of drinks/sitting, with five studies contributing to each analysis. Study-specific odds ratios (ORs) were estimated using logistic regression and pooled to generate adjusted summary ORs. Across studies, 0.9-3.2 % of control mothers reported drinking an average of 5+ drinks/sitting, while 1.4-23.5 % reported drinking a maximum of 5+ drinks/sitting. Compared with non-drinkers, mothers who drank an average of 5+ drinks/sitting were more likely to deliver an infant with cleft lip only (pooled OR 1.48; 95 % confidence intervals 1.01, 2.18). The estimate was higher among women who drank at this level 3+ times (pooled OR 1.95; 1.23, 3.11). Ever drinking a maximum of 5+ drinks/sitting and non-binge drinking were not associated with cleft risk. Repeated heavy maternal alcohol consumption was associated with an increased risk of cleft lip only in offspring. There was little evidence of increased risk for other cleft types or alcohol measures.

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

OBJECTIVES: In humans, there is a large range of variation in the form of the maxillary and mandibular dental arches. This variation can manifest as either prognathism or retrognathism in either or both arches, which can cause malocclusion and lead to abnormal masticatory function. This study aims to identify aspects of variation and morphological integration existing in the dental arches of individuals with different types of malocclusion. METHODS: Coordinate landmark data were collected along the gingival margins of 397 scanned dental casts and then analyzed using geometric morphometric techniques to explore arch form variation and patterns of morphological integration within each malocclusion type. RESULTS: Significant differences were identified between Class II forms (increased projection of upper arch relative to the lower arch) and Class III forms (lower arch projection beyond the upper arch) in symmetrical shape variation, including anteroposterior arch discrepancies and abnormal anterior arch divergence or convergence. Partial least squares analysis demonstrated that Class III dental arches have higher levels of covariance between upper and lower arches (RV = 0.91) compared to the dental arches of Class II (RV = 0.78) and Class I (RV = 0.73). These high levels of covariance, however, are on the lower end of the overall range of possible masticatory blocks, indicating weaker than expected levels of integration. CONCLUSIONS: This study provides evidence for patterns of variation in dental arch shape found in individuals with Class II and Class III malocclusions. Moreover, differences in integration found between malocclusion types have ramifications for how such conditions should be studied and treated. Am. J. Hum. Biol. 28:879-889, 2016. © 2016Wiley Periodicals, Inc.

Exploratory genotype-phenotype correlations of facial form and asymmetry in unaffected relatives of children with non-syndromic cleft lip and/or palate.

Family relatives of children with nonsyndromic cleft lip with or without cleft palate (NSCL/P) who presumably carry a genetic risk yet do not manifest overt oral clefts, often present with distinct facial morphology of unknown genetic etiology. This study investigates distinct facial morphology among unaffected relatives and examines whether candidate genes previously associated with overt NSCL/P and left-right body patterning are correlated with such facial morphology. Cases were unaffected relatives of individuals with NSCL/P (n = 188) and controls (n = 194) were individuals without family history of NSCL/P. Cases and controls were genotyped for 20 SNPs across 13 candidate genes for NSCL/P (PAX7, ABCA4-ARHGAP29, IRF6, MSX1, PITX2, 8q24, FOXE1, TGFB3 and MAFB) and left-right body patterning (LEFTY1, LEFTY2, ISL1 and SNAI1). Facial shape and asymmetry phenotypes were obtained via principal component analyses and Procrustes analysis of variance from 32 coordinate landmarks, digitized on 3D facial images. Case-control comparisons of phenotypes obtained were performed via multivariate regression adjusting for age and gender. Phenotypes that differed significantly (P < 0.05) between cases and controls were regressed on the SNPs one at a time. Cases had significantly (P < 0.05) more profile concavity with upper face retrusion, upturned noses with obtuse nasolabial angles, more protrusive chins, increased lower facial heights, thinner and more retrusive lips and more protrusive foreheads. Furthermore, cases showed significantly more directional asymmetry compared to controls. Several of these phenotypes were significantly associated with genetic variants (P < 0.05). Facial height and width were associated with SNAI1. Midface antero-posterior (AP) projection was associated with LEFTY1. The AP position of the chin was related to SNAI1, IRF6, MSX1 and MAFB. The AP position of the forehead and the width of the mouth were associated with ABCA4-ARHGAP29 and MAFB. Lastly, facial asymmetry was related to LEFTY1, LEFTY2 and SNAI1. This study demonstrates that, genes underlying lip and palate formation and left-right patterning also contribute to facial features characteristic of the NSCL/P spectrum.

Digital imaging analysis to assess scar phenotype.

In order to understand the link between the genetic background of patients and wound clinical outcomes, it is critical to have a reliable method to assess the phenotypic characteristics of healed wounds. In this study, we present a novel imaging method that provides reproducible, sensitive, and unbiased assessments of postsurgical scarring. We used this approach to investigate the possibility that genetic variants in orofacial clefting genes are associated with suboptimal healing. Red-green-blue digital images of postsurgical scars of 68 patients, following unilateral cleft lip repair, were captured using the 3dMD imaging system. Morphometric and colorimetric data of repaired regions of the philtrum and upper lip were acquired using ImageJ software, and the unaffected contralateral regions were used as patient-specific controls. Repeatability of the method was high with intraclass correlation coefficient score > 0.8. This method detected a very significant difference in all three colors, and for all patients, between the scarred and the contralateral unaffected philtrum (p ranging from 1.20(-05) to 1.95(-14) ). Physicians' clinical outcome ratings from the same images showed high interobserver variability (overall Pearson coefficient = 0.49) as well as low correlation with digital image analysis results. Finally, we identified genetic variants in TGFB3 and ARHGAP29 associated with suboptimal healing outcome.

Phenotypic diversity in white adults with moderate to severe Class II malocclusion.

INTRODUCTION: Class II malocclusion affects about 15% of the population in the United States and is characterized by a convex profile and occlusal disharmonies. The specific etiologic mechanisms resulting in the range of Class II dentoskeletal combinations observed are not yet understood. Most studies describing Class II phenotypic diversity have used moderate sample sizes or focused on younger patients who later in life might outgrow their Class II discrepancies; such a focus might also preclude the visualization of adult Class II features. The majority have used simple correlation methods resulting in phenotypes that might not be generalizable to different samples and thus might not be suitable for studies of malocclusion etiology. The purpose of this study was to address these knowledge gaps by capturing the maximum phenotypic variations in a large sample of white Class II subjects selected with strict eligibility criteria and rigorously standardized multivariate reduction analyses. METHODS: Sixty-three lateral cephalometric variables were measured from the pretreatment records of 309 white Class II adults (82 male, 227 female; ages, 16-60 years). Principal component analysis and cluster analysis were used to generate comprehensive phenotypes to identify the most homogeneous groups of subjects, reducing heterogeneity and improving the power of future malocclusion etiology studies. RESULTS: Principal component analysis resulted in 7 principal components that accounted for 81% of the variation. The first 3 components represented variation on mandibular rotation, maxillary incisor angulation, and mandibular length. The cluster analysis identified 5 distinct Class II phenotypes. CONCLUSIONS: A comprehensive spectrum of Class II phenotypic definitions was obtained that can be generalized to other samples to advance our efforts for identifying the etiologic factors underlying Class II malocclusion.