Imputation of orofacial clefting data identifies novel risk loci and sheds light on the genetic background of cleft lip ± cleft palate and cleft palate only.

Nonsyndromic cleft lip with or without cleft palate (nsCL/P) is among the most common human birth defects with multifactorial etiology. Here, we present results from a genome-wide imputation study of nsCL/P in which, after adding replication cohort data, four novel risk loci for nsCL/P are identified (at chromosomal regions 2p21, 14q22, 15q24 and 19p13). On a systematic level, we show that the association signals within this high-density dataset are enriched in functionally-relevant genomic regions that are active in both human neural crest cells (hNCC) and mouse embryonic craniofacial tissue. This enrichment is also detectable in hNCC regions primed for later activity. Using GCTA analyses, we suggest that 30% of the estimated variance in risk for nsCL/P in the European population can be attributed to common variants, with 25.5% contributed to by the 24 risk loci known to date. For each of these, we identify credible SNPs using a Bayesian refinement approach, with two loci harbouring only one probable causal variant. Finally, we demonstrate that there is no polygenic component of nsCL/P detectable that is shared with nonsyndromic cleft palate only (nsCPO). Our data suggest that, while common variants are strongly contributing to risk for nsCL/P, they do not seem to be involved in nsCPO which might be more often caused by rare deleterious variants. Our study generates novel insights into both nsCL/P and nsCPO etiology and provides a systematic framework for research into craniofacial development and malformation.

Further evidence for deletions in 7p14.1 contributing to nonsyndromic cleft lip with or without cleft palate.

BACKGROUND: Nonsyndromic cleft with or without cleft palate (nsCL/P) is a common birth defect. Although genome-wide association studies (GWAS) have identified numerous risk variants, a considerable fraction of the genetic heritability remains unknown. The aim of the present study was to replicate a previous finding that de novo deletions in a 62 kb region of chromosome 7p14 are a risk factor for nsCL/P, using an independent cohort. METHODS: Data from a published case-control GWAS cohort of 399 patients and 1318 controls were used. Copy number variant (CNV) detection in the 62 kb candidate region of 7p14 was performed using QuantiSNP. Putative CNVs in probands were verified and validated by quantitative polymerase chain reaction. Segregation analyses were performed in family members for whom DNA was available. RESULTS: Within the 62 kb candidate region, a deletion of 7.4 kb showed association with nsCL/P (13/387 cases, 20/1300 controls, plowest = 0.024, odds ratio = 2.22). In all families with a sporadic case (n = 3), the deletion occurred de novo. In multiplex families, both incomplete segregation and incomplete penetrance were observed. CONCLUSION: The present data support the hypothesis that deletions at 7p14 are a common risk factor for nsCL/P. Genome-wide CNV analyses in nsCL/P cohorts are warranted to explore the functional relevance of these deletions and their contribution to nsCL/P, and to determine exact breakpoints. Birth Defects Research (Part A) 106:767-772, 2016. © 2016 Wiley Periodicals, Inc.

Sequencing the GRHL3 Coding Region Reveals Rare Truncating Mutations and a Common Susceptibility Variant for Nonsyndromic Cleft Palate.

Nonsyndromic cleft lip with/without cleft palate (nsCL/P) and nonsyndromic cleft palate only (nsCPO) are the most frequent subphenotypes of orofacial clefts. A common syndromic form of orofacial clefting is Van der Woude syndrome (VWS) where individuals have CL/P or CPO, often but not always associated with lower lip pits. Recently, ∼5% of VWS-affected individuals were identified with mutations in the grainy head-like 3 gene (GRHL3). To investigate GRHL3 in nonsyndromic clefting, we sequenced its coding region in 576 Europeans with nsCL/P and 96 with nsCPO. Most strikingly, nsCPO-affected individuals had a higher minor allele frequency for rs41268753 (0.099) than control subjects (0.049; p = 1.24 × 10(-2)). This association was replicated in nsCPO/control cohorts from Latvia, Yemen, and the UK (pcombined = 2.63 × 10(-5); ORallelic = 2.46 [95% CI 1.6-3.7]) and reached genome-wide significance in combination with imputed data from a GWAS in nsCPO triads (p = 2.73 × 10(-9)). Notably, rs41268753 is not associated with nsCL/P (p = 0.45). rs41268753 encodes the highly conserved p.Thr454Met (c.1361C>T) (GERP = 5.3), which prediction programs denote as deleterious, has a CADD score of 29.6, and increases protein binding capacity in silico. Sequencing also revealed four novel truncating GRHL3 mutations including two that were de novo in four families, where all nine individuals harboring mutations had nsCPO. This is important for genetic counseling: given that VWS is rare compared to nsCPO, our data suggest that dominant GRHL3 mutations are more likely to cause nonsyndromic than syndromic CPO. Thus, with rare dominant mutations and a common risk variant in the coding region, we have identified an important contribution for GRHL3 in nsCPO.

Meta-analysis Reveals Genome-Wide Significance at 15q13 for Nonsyndromic Clefting of Both the Lip and the Palate, and Functional Analyses Implicate GREM1 As a Plausible Causative Gene.

Nonsyndromic orofacial clefts are common birth defects with multifactorial etiology. The most common type is cleft lip, which occurs with or without cleft palate (nsCLP and nsCLO, respectively). Although genetic components play an important role in nsCLP, the genetic factors that predispose to palate involvement are largely unknown. In this study, we carried out a meta-analysis on genetic and clinical data from three large cohorts and identified strong association between a region on chromosome 15q13 and nsCLP (P = 8.13 × 10(-14) for rs1258763; relative risk (RR): 1.46, 95% confidence interval (CI): 1.32-1.61)) but not nsCLO (P = 0.27; RR: 1.09 (0.94-1.27)). The 5 kb region of strongest association maps downstream of Gremlin-1 (GREM1), which encodes a secreted antagonist of the BMP4 pathway. We show during mouse embryogenesis, Grem1 is expressed in the developing lip and soft palate but not in the hard palate. This is consistent with genotype-phenotype correlations between rs1258763 and a specific nsCLP subphenotype, since a more than two-fold increase in risk was observed in patients displaying clefts of both the lip and soft palate but who had an intact hard palate (RR: 3.76, CI: 1.47-9.61, Pdiff<0.05). While we did not find lip or palate defects in Grem1-deficient mice, wild type embryonic palatal shelves developed divergent shapes when cultured in the presence of ectopic Grem1 protein (P = 0.0014). The present study identified a non-coding region at 15q13 as the second, genome-wide significant locus specific for nsCLP, after 13q31. Moreover, our data suggest that the closely located GREM1 gene contributes to a rare clinical nsCLP entity. This entity specifically involves abnormalities of the lip and soft palate, which develop at different time-points and in separate anatomical regions.

Replication analysis of 15 susceptibility loci for nonsyndromic cleft lip with or without cleft palate in an italian population.

BACKGROUND: Nonsyndromic cleft lip with or without cleft palate (nsCL/P) is one of the most common congenital malformations in humans. Its average global incidence is 1.7 per 1000 live births, with wide variation according to geographical location and ethnicity. Its etiology involves both genetic and environmental factors. The aim of the present study was to confirm genetic association of a selection of 15 candidate nsCL/P loci using an independent sample of the Italian population. METHODS: At least one single-nucleotide polymorphism (SNP) for each locus was genotyped in 380 nuclear trios. RESULTS: Transmission disequilibrium analysis revealed significant associations for three variants at two loci (8q24 and 1p22). Two SNPs at 8q24 showed the strongest level of association, the rs987525 (PTDT = 6.81 × 10(-6) ; homozygous relative risk = 3.60 [95% confidence interval, 2.12-6.13]), and the rs17241253 (PTDT = 1.03 × 10(-5) ; homozygous relative risk = 3.75 [95% confidence interval, 2.10-6.67]). Four additional loci (at 1q32, 3q12, 8q21, and 10q25) achieved nominally significant p-values. Two SNPs at 1p36 achieved p-values of < 0.1. The present data suggest that 6 of the 15 analyzed nsCL/P risk loci contribute significantly to nsCL/P risk in the Italian population. These include the 1p22 locus, which previous research has implicated predominantly in Asian populations. CONCLUSION: Different loci, including 8q24 and 1p22 have been found associated with nsCL/P in multiple populations. Further efforts are needed to identify causative variants and transfer knowledge to clinical application, such as personal genetic risk assessment.