Orofacial clefts (OFCs), which include non-syndromic cleft lip with or without cleft palate (CL/P), are among the most common birth defects in humans, affecting approximately 1 in 700 newborns. CL/P is phenotypically heterogeneous and has a complex etiology caused by genetic and environmental factors. Previous genome-wide association studies (GWASs) have identified at least 15 risk loci for CL/P. As these loci do not account for all of the genetic variance of CL/P, we hypothesized the existence of additional risk loci. We conducted a multiethnic GWAS in 6480 participants (823 unrelated cases, 1700 unrelated controls and 1319 case-parent trios) with European, Asian, African and Central and South American ancestry. Our GWAS revealed novel associations on 2p24 near FAM49A, a gene of unknown function (P = 4.22 × 10), and 19q13 near RHPN2, a gene involved in organizing the actin cytoskeleton (P = 4.17 × 10). Other regions reaching genome-wide significance were 1p36 (PAX7), 1p22 (ARHGAP29), 1q32 (IRF6), 8q24 and 17p13 (NTN1), all reported in previous GWASs. Stratification by ancestry group revealed a novel association with a region on 17q23 (P = 2.92 × 10) among individuals with European ancestry. This region included several promising candidates including TANC2, an oncogene required for development, and DCAF7, a scaffolding protein required for craniofacial development. In the Central and South American ancestry group, significant associations with loci previously identified in Asian or European ancestry groups reflected their admixed ancestry. In summary, we have identified novel CL/P risk loci and suggest new genes involved in craniofacial development, confirming the highly heterogeneous etiology of OFCs.

Nonsyndromic orofacial clefts (OFCs) are a heterogeneous group of common craniofacial birth defects with complex etiologies that include genetic and environmental risk factors. OFCs are commonly categorized as cleft lip with or without cleft palate (CL/P) and cleft palate alone (CP), which have historically been analyzed as distinct entities. Genes for both CL/P and CP have been identified via multiple genome-wide linkage and association studies (GWAS), however, altogether, known variants account for a minority of the estimated heritability in risk to these craniofacial birth defects. We performed genome-wide meta-analyses of CL/P, CP, and all OFCs across two large, multiethnic studies. We then performed population specific meta-analyses in sub-samples of Asian and European ancestry. In addition to observing associations with known variants, we identified a novel genome-wide significant association between SNPs located in an intronic TP63 enhancer and CL/P (p = 1.16 × 10−8). Several novel loci with compelling candidate genes approached genome-wide significance on 4q21.1 (SHROOM3), 12q13.13 (KRT18), and 8p21 (NRG1). In the analysis of all OFCs combined, SNPs near FOXE1 reached genome-wide significance (p = 1.33 × 10−9). Our results support the highly heterogeneous nature of OFCs and illustrate the utility of meta-analysis for discovering new genetic risk factors.

Cleft palate (CP) is a common birth defect occurring in 1 in 2,500 live births. Approximately half of infants with CP have a syndromic form, exhibiting other physical and cognitive disabilities. The other half have nonsyndromic CP, and to date, few genes associated with risk for nonsyndromic CP have been characterized. To identify such risk factors, we performed a genome-wide association study of this disorder. We discovered a genome-wide significant association with a missense variant in GRHL3 (p.Thr454Met [c.1361C>T]; rs41268753; p = 4.08 × 10(-9)) and replicated the result in an independent sample of case and control subjects. In both the discovery and replication samples, rs41268753 conferred increased risk for CP (OR = 8.3, 95% CI 4.1-16.8; OR = 2.16, 95% CI 1.43-3.27, respectively). In luciferase transactivation assays, p.Thr454Met had about one-third of the activity of wild-type GRHL3, and in zebrafish embryos, perturbed periderm development. We conclude that this mutation is an etiologic variant for nonsyndromic CP and is one of few functional variants identified to date for nonsyndromic orofacial clefting. This finding advances our understanding of the genetic basis of craniofacial development and might ultimately lead to improvements in recurrence risk prediction, treatment, and prognosis.

Phenotypic heterogeneity is a hallmark of complex traits, and genetic studies of such traits may focus on them as a single diagnostic entity or by analyzing specific components. For example, in orofacial clefting (OFC), three subtypescleft lip (CL), cleft lip and palate (CLP), and cleft palate (CP) have been studied separately and in combination. To further dissect the genetic architecture of OFCs and how a given associated locus may be contributing to distinct subtypes of a trait we developed a framework for quantifying and interpreting evidence of subtype-specific or shared genetic effects in complex traits. We applied this technique to create a "cleft map" of the association of 30 genetic loci with three [TB]

Orofacial clefts (OFCs) are congenital dysmorphologies of the human face and oral cavity, with a global incidence of 1 per 700 live births. These anomalies exhibit a multifactorial pattern of inheritance, with genetic and environmental factors both playing crucial roles. Many loci have been implicated in the etiology of nonsyndromic cleft lip with or without cleft palate (NSCL/P) in populations of Asian and European ancestries, through genome-wide association studies and candidate gene studies. However, few populations of African descent have been studied to date. Here, the authors show evidence of an association of some loci with NSCL/P and nonsyndromic cleft palate only (NSCPO) in cohorts from Africa (Ghana, Ethiopia, and Nigeria). The authors genotyped 48 single-nucleotide polymorphisms that were selected from previous genome-wide association studies and candidate gene studies. These markers were successfully genotyped on 701 NSCL/P and 163 NSCPO cases, 1,070 unaffected relatives, and 1,078 unrelated controls. The authors also directly sequenced 7 genes in 184 nonsyndromic OFC (NSOFC) cases and 96 controls from Ghana. Population-specific associations were observed in the casecontrol analyses of the subpopulations, with West African subpopulations (Ghana and Nigeria) showing a similar pattern of associations. In meta-analyses of the case-control cohort, PAX7 (rs742071, P = 5.10 × 10 -3 ), 8q24 (rs987525, P = 1.22 × 10 -3 ), and VAX1 (rs7078160, P = 0.04) were nominally associated with NSCL/P, and MSX1 (rs115200552, P = 0.01), TULP4 (rs651333, P = 0.04), CRISPLD2 (rs4783099, P = 0.02), and NOG1 (rs17760296, P = 0.04) were nominally associated with NSCPO. Moreover, 7 loci exhibited evidence of threshold overtransmission in NSOFC cases through the transmission disequilibrium test and through analyses of the family-based association for disease traits. Through DNA sequencing, the authors also identified 2 novel, rare, potentially pathogenic variants (p.Asn323Asp and p.Lys426IlefsTer6) in ARHGAP29. In conclusion, the authors have shown evidence for the association of many loci with NSCL/P and NSCPO. To the best of this knowledge, this study is the first to demonstrate any of these association signals in any African population.