Non-syndromic cleft lip with or without cleft palate (nsCL/P) ranks among the most common human congenital malformations, and has a multifactorial background in which both exogenous and genetic risk factors act in concert. The present report describes a genome-wide association study (GWAS) involving a total of 285 nsCL/P patients and 1212 controls from the Netherlands and Belgium. Twenty of the 40 previously reported nsC/LP susceptibility loci were replicated, which underlined the validity of this sample. SNV-based analysis of the data identified an as yet unreported suggestive locus at chromosome 16p12.1 (p-value of the lead SNV: 4.17 × 10⁻⁷). This association was replicated in two of three patient/control replication series (Central European and Yemeni). Gene analysis of the GWAS data prioritized SH3PXD2A at chromosome 10q24.33 as a candidate gene for nsCL/P. To date, support for this gene as a cleft gene has been restricted to data from zebrafish and a knockout mouse model. The present GWAS was the first to implicate SH3PXD2A in non-syndromic cleft formation in humans. In summary, although performed in a relatively small sample, the present GWAS generated novel insights into nsCL/P etiology.

CELSR3 codes for a planar cell polarity protein. We describe twelve affected individuals from eleven independent families with bi-allelic variants in CELSR3. Affected individuals presented with an overlapping phenotypic spectrum comprising central nervous system (CNS) anomalies (7/12), combined CNS anomalies and congenital anomalies of the kidneys and urinary tract (CAKUT) (3/12) and CAKUT only (2/12). Computational simulation of the 3D protein structure suggests the position of the identified variants to be implicated in penetrance and phenotype expression. CELSR3 immunolocalization in human embryonic urinary tract and transient suppression and rescue experiments of Celsr3 in fluorescent zebrafish reporter lines further support an embryonic role of CELSR3 in CNS and urinary tract formation.

Orofacial clefting (OFC) is a frequent congenital anomaly and can occur either in the context of underlying syndromes or in isolation (nonsyndromic). The two common OFC phenotypes are cleft lip with/without cleft palate (CL/P) and cleft palate only (CPO). In this study, we searched for penetrant CL/P genes, by evaluating de novo copy number variants (CNV) from an exome sequencing dataset of 50 nonsyndromic patient-parent trios. We detected a heterozygous 86 kb de novo deletion affecting exons 4–11 of ZFHX4, a gene previously associated with OFC. Genetic and phenotypic data from our in-house and the AGORA cohort (710 and 229 individuals with nonsyndromic CL/P) together with literature and database reviews demonstrate that ZFHX4 variants can lead to both nonsyndromic and syndromic forms not only of CL/P but also CPO. Expression analysis in published single-cell RNA-sequencing data (mouse embryo, zebrafish larva) at relevant time-points support an important role of Zfhx4/zfhx4 in craniofacial development. To characterize the role of zfhx4 in zebrafish craniofacial development, we knocked out/down the zebrafish orthologue. Cartilage staining of the zfhx4 CRISPR F0 knockout and morpholino knockdown at 4 days post-fertilization showed an underdeveloped and abnormally shaped ethmoid plate and cartilaginous jaw (resembling micrognathia). While there is evidence for the dominant inheritance of ZFHX4 variants in OFC, we here present a patient with a possible recessive inheritance. In conclusion, ZFHX4 has a highly heterogeneous phenotypic spectrum and variable mode of inheritance. Our data highlight that ZFHX4 should be considered in genetic testing in patients with nonsyndromic clefting.

Background: Nonsyndromic cleft lip with/without cleft palate (nsCL/P) is a congenital malformation of multifactorial etiology. Research has identified >40 genome‐wide significant risk loci, which explain less than 40% of nsCL/P heritability. Studies show that some of the hidden heritability is explained by rare penetrant variants.Methods: To identify new candidate genes, we searched for highly penetrant de novo variants (DNVs) in 50 nsCL/P patient/parent‐trios with a low polygenic risk for the phenotype (discovery).We prioritized DNV‐carrying candidate genes from the discovery for resequencing in independent cohorts of 1010 nsCL/P patients of diverse ethnicities and 1574 population‐matched controls (replication). Segregation analyses and rare variant association in the replication cohort, in combination with additional data (genome‐wide association data, expression, protein–protein‐interactions), were used for final prioritization.Conclusion: In the discovery step, 60 DNVs were identified in 60 genes, including a variant in the established nsCL/P risk gene CDH1. Re‐sequencing of 32 prioritized genes led to the identification of 373 rare, likely pathogenic variants. Finally, MDN1 and PAXIP1 were prioritized as top candidates. Our findings demonstrate that DNV detection, including polygenic risk score analysis, is a powerful tool for identifying nsCL/P candidate genes, which can also be applied to other multifactorial congenital malformations.