Abstract
Orofacial clefts (OFCs) are one of the most common structural birth defects, with the prevalence of OFC varying across populations, and studies on the causes of OFCs in diverse populations are necessary, but still limited. We analyzed whole genome sequencing data on 419 parent-child trios from the Philippines, a population with a particularly high rate of OFC. To identify novel genes for OFCs, we studied both common variation and de novo variants (DNVs). We identified a significant enrichment in both loss-of-function (N = 62; P = 8.34 × 10-5) and protein-altering DNVs (N = 394; P = 1.49 × 10-7) among OFC probands. Among the genes individually enriched for DNVs was GRHL2 (P = 6.60 × 10-6), where there were two DNVs, a stop-gain and a frameshift deletion. We then queried OFC trios from other cohorts in the Gabriella Miller Kids First program (total N = 1254) and GeneMatcher and identified an 89 kb de novo deletion in GRHL2 and a de novo 8q22.3 microdeletion with one breakpoint in GRHL2. Additionally, within the common variant analyses we found significant gene x gene interactions with GRHL2. GRHL2 is a conserved transcription factor involved in embryonic development, with truncating mutations causing autosomal dominant progressive hearing loss and missense variants causing autosomal recessive ectodermal dysplasia. Heterozygous variation in its homolog, GRHL3, causes Van der Woude syndrome and isolated cleft palate. Additionally, mice deficient for either Grhl2 or Grhl3 have craniofacial anomalies, including facial and palatal clefts, strongly supporting GRHL2 as a risk locus for OFCs.