Exome sequencing provides additional evidence for the involvement of ARHGAP29 in Mendelian orofacial clefting and extends the phenotypic spectrum to isolated cleft palate

The heterogeneous genetic architectures of orofacial clefts

Orofacial clefts (OFCs) are common, affecting 1:1000 live births. OFCs occur across a phenotypic spectrum - including cleft lip (CL), cleft lip and palate (CLP), or cleft palate (CP) - and can be further subdivided based on laterality, severity, or specific structures affected. Herein we review what is known about the genetic architecture underlying each of these subtypes, considering both shared and subtype-specific risks. While there are more known genetic similarities between CL and CLP than CP, recent research supports both shared and subtype-specific genetic risk factors within and between phenotypic classifications of OFCs. Larger sample sizes and deeper phenotyping data will be of increasing importance for the discovery of novel genetic risk factors for OFCs and various subtypes going forward.

Shared genetic risk between major orofacial cleft phenotypes in an African population

Nonsyndromic orofacial clefts (NSOFCs) represent a large proportion (70%-80%) of all OFCs. They can be broadly categorized into nonsyndromic cleft lip with or without cleft palate (NSCL/P) and nonsyndromic cleft palate only (NSCPO). Although NSCL/P and NSCPO are considered etiologically distinct, recent evidence suggests the presence of shared genetic risks. Thus, we investigated the genetic overlap between NSCL/P and NSCPO using African genome-wide association study (GWAS) data on NSOFCs. These data consist of 814 NSCL/P, 205 NSCPO cases, and 2159 unrelated controls. We generated common single-nucleotide variants (SNVs) association summary statistics separately for each phenotype (NSCL/P and NSCPO) under an additive genetic model. Subsequently, we employed the pleiotropic analysis under the composite null (PLACO) method to test for genetic overlap. Our analysis identified two loci with genome-wide significance (rs181737795 [p = 2.58E-08] and rs2221169 [p = 4.5E-08]) and one locus with marginal significance (rs187523265 [p = 5.22E-08]). Using mouse transcriptomics data and information from genetic phenotype databases, we identified MDN1, MAP3k7, KMT2A, ARCN1, and VADC2 as top candidate genes for the associated SNVs. These findings enhance our understanding of genetic variants associated with NSOFCs and identify potential candidate genes for further exploration.

Damaging Mutations in Contribute to Risk of Nonsyndromic Cleft Lip With or Without Cleft Palate

Novel or rare damaging mutations have been implicated in the developmental pathogenesis of nonsyndromic cleft lip with or without cleft palate (nsCL ± P). Thus, we investigated the human genome for high-impact mutations that could explain the risk of nsCL ± P in our cohorts.

We conducted next-generation sequencing (NGS) analysis of 130 nsCL ± P case-parent African trios to identify pathogenic variants that contribute to the risk of clefting. We replicated this analysis using whole-exome sequence data from a Brazilian nsCL ± P cohort. Computational analyses were then used to predict the mechanism by which these variants could result in increased risks for nsCL ± P.

We discovered damaging mutations within the AFDN gene, a cell adhesion molecule (CAMs) that was previously shown to contribute to cleft palate in mice. These mutations include p.Met1164Ile, p.Thr453Asn, p.Pro1638Ala, p.Arg669Gln, p.Ala1717Val, and p.Arg1596His. We also discovered a novel splicing p.Leu1588Leu mutation in this protein. Computational analysis suggests that these amino acid changes affect the interactions with other cleft-associated genes including nectins (PVRL1, PVRL2, PVRL3, and PVRL4) CDH1, CTNNA1, and CTNND1.

This is the first report on the contribution of AFDN to the risk for nsCL ± P in humans. AFDN encodes AFADIN, an important CAM that forms calcium-independent complexes with nectins 1 and 4 (encoded by the genes PVRL1 and PVRL4). This discovery shows the power of NGS analysis of multiethnic cleft samples in combination with a computational approach in the understanding of the pathogenesis of nsCL ± P.

CACNA1S mutation-associated dental anomalies: A calcium channelopathy

OBJECTIVES

To identify the molecular etiology of distinct dental anomalies found in eight Thai patients and explore the mutational effects on cellular functions.

MATERIALS AND METHODS

Clinical and radiographic examinations were performed for eight patients. Whole exome sequencing, mutant protein modelling, qPCR, western blot analysis, scratch assays, immunofluorescence, confocal analysis, in situ hybridization, and scanning electron micrography of teeth were done.

RESULTS

All patients had molars with multiple supernumerary cusps, single-cusped premolars, and a reduction in root number. Mutation analysis highlighted a heterozygous c.865A>G; p.Ile289Val mutation in CACNA1S in the patients. CACNA1S is a component of the slowly inactivating L-type voltage-dependent calcium channel. Mutant protein modeling suggested that the mutation might allow leakage of Ca2+ or other cations, or a tightening, to restrict calcium flow. Immunohistochemistry analysis showed expression of Cacna1s in the developing murine tooth epithelium during stages of crown and root morphogenesis. In cell culture, the mutation resulted in abnormal cell migration of transfected CHO cells compared to wildtype CACNA1S, with changes to the cytoskeleton and markers of focal adhesion.

CONCLUSIONS

The malformations observed in our patients suggest a role for calcium signaling in organization of both cusps and roots, affecting cell dynamics within the dental epithelium.

Rare variants found in multiplex families with orofacial clefts: Does expanding the phenotype make a difference?

Exome sequencing (ES) is now a relatively straightforward process to identify causal variants in Mendelian disorders. However, the same is not true for ES in families where the inheritance patterns are less clear, and a complex etiology is suspected. Orofacial clefts (OFCs) are highly heritable birth defects with both Mendelian and complex etiologies. The phenotypic spectrum of OFCs may include overt clefts and several subclinical phenotypes, such as discontinuities in the orbicularis oris muscle (OOM) in the upper lip, velopharyngeal insufficiency (VPI), microform clefts or bifid uvulas. We hypothesize that expanding the OFC phenotype to include these phenotypes can clarify inheritance patterns in multiplex families, making them appear more Mendelian. We performed exome sequencing to find rare, likely causal genetic variants in 31 multiplex OFC families, which included families with multiple individuals with OFCs and individuals with subclinical phenotypes. We identified likely causal variants in COL11A2, IRF6, SHROOM3, SMC3, TBX3, and TP63 in six families. Although we did not find clear evidence supporting the subclinical phenotype hypothesis, our findings support a role for rare variants in the etiology of OFCs.

Multimodal spatiotemporal transcriptomic resolution of embryonic palate osteogenesis

The terminal differentiation of osteoblasts and subsequent formation of bone marks an important phase in palate development that leads to the separation of the oral and nasal cavities. While the morphogenetic events preceding palatal osteogenesis are well explored, major gaps remain in our understanding of the molecular mechanisms driving the formation of this bony union of the fusing palate. Through bulk, single-nucleus, and spatially resolved RNA-sequencing analyses of the developing secondary palate, we identify a shift in transcriptional programming between embryonic days 14.5 and 15.5 pinpointing the onset of osteogenesis. We define spatially restricted expression patterns of key osteogenic marker genes that are differentially expressed between these developmental timepoints. Finally, we identify genes in the palate highly expressed by palate nasal epithelial cells, also enriched within palatal osteogenic mesenchymal cells. This investigation provides a relevant framework to advance palate-specific diagnostic and therapeutic biomarker discovery.

Genetics and orofacial clefts: a clinical perspective

Orofacial clefts (OFCs) are the most common congenital craniofacial anomaly seen in humans. Most OFCs are sporadic and isolated - these are thought to be multifactorial in origin. Chromosomal and monogenic variants account for the syndromic forms and for some of the non-syndromic inherited forms. This review discusses the importance of genetic testing and the current clinical strategy to deliver a genomics service that is of direct benefit to patients and their families.

Integrated spatiotemporal transcriptomic resolution of embryonic palate osteogenesis

The differentiation of osteoblasts and the subsequent formation of bone marks an important terminal phase in palate formation that leads to the separation of the oral and nasal cavities. While the developmental events that precede palatal osteogenesis are well explored, major gaps remain in our understanding of the molecular mechanisms that lead to the bony union of fusing palatal shelves. Herein, the timeline of osteogenic transcriptional programming is unveiled in the embryonic palate by way of integrated bulk, single-cell, and spatially resolved RNA-seq analyses. We define spatially restricted expression patterns of key marker genes, both regulatory and structural, that are differentially expressed during palatal fusion, including the identification of several novel genes ( Deup1, Dynlrb2, Lrrc23 ) spatially restricted in expression to the palate, providing a relevant framework for future studies that identify new candidate genes for cleft palate anomalies in humans as well as the timing of mammalian embryonic palatal osteogenesis.

Identification of Novel Risk Variants of Non-Syndromic Cleft Palate by Targeted Gene Panel Sequencing

Non-syndromic cleft palate (ns-CP) has a genetically heterogeneous aetiology. Numerous studies have suggested a crucial role of rare coding variants in characterizing the unrevealed component of genetic variation in ns-CP called the "missing heritability". Therefore, this study aimed to detect low-frequency variants that are implicated in ns-CP aetiology in the Polish population. For this purpose, coding regions of 423 genes associated with orofacial cleft anomalies and/or involved with facial development were screened in 38 ns-CP patients using the next-generation sequencing technology. After multistage selection and prioritisation, eight novel and four known rare variants that may influence an individual's risk of ns-CP were identified. Among detected alternations, seven were located in novel candidate genes for ns-CP, including COL17A1 (c.2435-1G>A), DLG1 (c.1586G>C, p.Glu562Asp), NHS (c.568G>C, p.Val190Leu-de novo variant), NOTCH2 (c.1997A>G, p.Tyr666Cys), TBX18 (c.647A>T, p.His225Leu), VAX1 (c.400G>A, p.Ala134Thr) and WNT5B (c.716G>T, p.Arg239Leu). The remaining risk variants were identified within genes previously linked to ns-CP, confirming their contribution to this anomaly. This list included ARHGAP29 (c.1706G>A, p.Arg569Gln), FLNB (c.3605A>G, Tyr1202Cys), IRF6 (224A>G, p.Asp75Gly-de novo variant), LRP6 (c.481C>A, p.Pro161Thr) and TP63 (c.353A>T, p.Asn118Ile). In summary, this study provides further insights into the genetic components contributing to ns-CP aetiology and identifies novel susceptibility genes for this craniofacial anomaly.

A Comprehensive Genetic Analysis of Slovenian Families with Multiple Cases of Orofacial Clefts Reveals Novel Variants in the Genes IRF6, GRHL3, and TBX22

Although the aetiology of non-syndromic orofacial clefts (nsOFCs) is usually multifactorial, syndromic OFCs (syOFCs) are often caused by single mutations in known genes. Some syndromes, e.g., Van der Woude syndrome (VWS1; VWS2) and X-linked cleft palate with or without ankyloglossia (CPX), show only minor clinical signs in addition to OFC and are sometimes difficult to differentiate from nsOFCs. We recruited 34 Slovenian multi-case families with apparent nsOFCs (isolated OFCs or OFCs with minor additional facial signs). First, we examined IRF6, GRHL3, and TBX22 by Sanger or whole exome sequencing to identify VWS and CPX families. Next, we examined 72 additional nsOFC genes in the remaining families. Variant validation and co-segregation analysis were performed for each identified variant using Sanger sequencing, real-time quantitative PCR and microarray-based comparative genomic hybridization. We identified six disease-causing variants (three novel) in IRF6, GRHL3, and TBX22 in 21% of families with apparent nsOFCs, suggesting that our sequencing approach is useful for distinguishing syOFCs from nsOFCs. The novel variants, a frameshift variant in exon 7 of IRF6, a splice-altering variant in GRHL3, and a deletion of the coding exons of TBX22, indicate VWS1, VWS2, and CPX, respectively. We also identified five rare variants in nsOFC genes in families without VWS or CPX, but they could not be conclusively linked to nsOFC.

Rare variants found in multiplex families with orofacial clefts: Does expanding the phenotype make a difference?

Whole-exome sequencing (WES) is now a relatively straightforward process to identify causal variants in Mendelian disorders. However, the same is not true for WES in families where the inheritance patterns are less clear, and a complex etiology is suspected. Orofacial clefts (OFCs) are highly heritable birth defects with both Mendelian and complex etiologies. The phenotypic spectrum of OFCs may include overt clefts and several subclinical phenotypes, such as discontinuities in the orbicularis oris muscle (OOM) in the upper lip, velopharyngeal insufficiency (VPI), microform clefts or bifid uvulas. We hypothesize that expanding the OFC phenotype to include these phenotypes can clarify inheritance patterns in multiplex families, making them appear more Mendelian. We performed whole-exome sequencing to find rare, likely causal genetic variants in 31 multiplex OFC families, which included families with multiple individuals with OFCs and individuals with subclinical phenotypes. We identified likely causal variants in COL11A2, IRF6, KLF4, SHROOM3, SMC3, TP63 , and TBX3 in seven families. Although we did not find clear evidence supporting the subclinical phenotype hypothesis, our findings support a role for rare variants in the etiology of OFCs.

ARHGAP29 is required for keratinocyte proliferation and migration

BACKGROUND

RhoA GTPase plays critical roles in actin cytoskeletal remodeling required for controlling a diverse range of cellular functions including cell proliferation, cell adhesions, migration and changes in cell shape. RhoA cycles between an active GTP-bound and an inactive GDP-bound form, a process that is regulated by guanine nucleotide exchange factors (GEFs), and GTPase-activating proteins (GAPs). ARHGAP29 is a GAP expressed in keratinocytes of the skin and is decreased in the absence of Interferon Regulator Factor 6, a critical regulator of cell proliferation and migration. However, the role for ARHGAP29 in keratinocyte biology is unknown.

RESULTS

Novel ARHGAP29 knockdown keratinocyte cell lines were generated using both CRISPR/Cas9 and shRNA technologies. Knockdown cells exhibited significant reduction of ARHGAP29 protein (50-80%) and displayed increased filamentous actin (stress fibers), phospho-myosin light chain (contractility), cell area and population doubling time. Furthermore, we found that ARHGAP29 knockdown keratinocytes displayed significant delays in scratch wound closure in both single cell and collective cell migration conditions. Particularly, our results show a reduction in path lengths, speed, directionality and persistence in keratinocytes with reduced ARHGAP29. The delay in scratch closure was rescued by both adding back ARHGAP29 or adding a ROCK inhibitor to ARHGAP29 knockdown cells.

CONCLUSIONS

These data demonstrate that ARHGAP29 is required for keratinocyte morphology, proliferation and migration mediated through the RhoA pathway.

Ultra-Rare Variants Identify Biological Pathways and Candidate Genes in the Pathobiology of Non-Syndromic Cleft Palate Only

In recent decades, many efforts have been made to elucidate the genetic causes of non-syndromic cleft palate (nsCPO), a complex congenital disease caused by the interaction of several genetic and environmental factors. Since genome-wide association studies have evidenced a minor contribution of common polymorphisms in nsCPO inheritance, we used whole exome sequencing data to explore the role of ultra-rare variants in this study. In a cohort of 35 nsCPO cases and 38 controls, we performed a gene set enrichment analysis (GSEA) and a hypergeometric test for assessing significant overlap between genes implicated in nsCPO pathobiology and genes enriched in ultra-rare variants in our cohort. GSEA highlighted an enrichment of ultra-rare variants in genes principally belonging to cytoskeletal protein binding pathway (Probability Density Function corrected p-value = 1.57 × 10^-4); protein-containing complex binding pathway (p-value = 1.06 × 10^-2); cell adhesion molecule binding pathway (p-value = 1.24 × 10^-2); ECM-receptor interaction pathway (p-value = 1.69 × 10^-2); and in the Integrin signaling pathway (p-value = 1.28 × 10^-2). Two genes implicated in nsCPO pathobiology, namely COL2A1 and GLI3, ranked among the genes (n = 34) with nominal enrichment in the ultra-rare variant collapsing analysis (Fisher's exact test p-value < 0.05). These genes were also part of an independent list of genes highly relevant to nsCPO biology (n = 25). Significant overlap between the two sets of genes (hypergeometric test p-value = 5.86 × 10^-3) indicated that enriched genes are likely to be implicated in physiological palate development and/or the pathological processes of oral clefting. In conclusion, ultra-rare variants collectively impinge on biological pathways crucial to nsCPO pathobiology and point to candidate genes that may contribute to the individual risk of disease. Sequencing can be an effective approach to identify candidate genes and pathways for nsCPO.

Variant analyses of candidate genes in orofacial clefts in multi-ethnic populations

OBJECTIVES

Cleft lip with/without cleft palate and cleft palate only is congenital birth defects where the upper lip and/or palate fail to fuse properly during embryonic facial development. Affecting ~1.2/1000 live births worldwide, these orofacial clefts impose significant social and financial burdens on affected individuals and their families. Orofacial clefts have a complex etiology resulting from genetic variants combined with environmental covariates. Recent genome-wide association studies and whole-exome sequencing for orofacial clefts identified significant genetic associations and variants in several genes. Of these, we investigated the role of common/rare variants in SHH, RORA, MRPL53, ACVR1, and GDF11.

MATERIALS AND METHODS

We sequenced these five genes in 1255 multi-ethnic cleft lip with/without palate and cleft palate only samples in order to find variants that may provide potential explanations for the missing heritability of orofacial clefts. Rare and novel variants were further analyzed using in silico predictive tools.

RESULTS

Ninteen total variants of interest were found, with variant types including stop-gain, missense, synonymous, intronic, and splice-site variants. Of these, 3 novel missense variants were found, one in SHH, one in RORA, and one in GDF11.

CONCLUSION

This study provides evidence that variants in SHH, RORA, MRPL53, ACVR1, and GDF11 may contribute to risk of orofacial clefts in various populations.

Targeted re-sequencing on 1p22 among non-syndromic orofacial clefts from Han Chinese population

Rs560426 at 1p22 was proved to be associated with NSCL/P (non-syndromic cleft lip with or without the palate) in several populations, including Han Chinese population. Here, we conducted a deep sequencing around rs560426 to locate more susceptibility variants in this region. In total, 2,293 NSCL/P cases and 3,235 normal controls were recruited. After sequencing, association analysis was performed. Western blot, RT-qPCR, HE, immunofluorescence staining, and RNA sequencing were conducted for functional analyses of the selected variants. Association analysis indicated that rs77179923 was the only SNP associated with NSCLP specifically (p = 4.70E-04, OR = 1.84), and rs12071152 was uniquely associated with LCLO (p = 4.00E-04, OR = 1.30, 95%CI: 1.12–1.51). Moreover, de novo harmful rare variant NM_004815.3, NP_004806.3; c.1652G>C, p.R551T in ARHGAP29 resulted in a decreased expression level of ARHGAP29, which in turn affected NSCL/P-related biological processes; however, no overt cleft palate (CP) phenotype was observed. In conclusion, rs12071152 was a new susceptible variant, which is specifically associated with LCLO among the Han Chinese population. Allele A of it could increase the risk of having a cleft baby. Rs77179923 and rare variant NM_004815.3, NP_004806.3; c.1652G>C, p.R551T at 1p22 were both associated with NSCLP among the Han Chinese population. However, this missense variation contributes to no overt CP phenotype due to dosage insufficiency or compensation from other genes.

Whole-genome sequencing reveals de-novo mutations associated with nonsyndromic cleft lip/palate

The majority (85%) of nonsyndromic cleft lip with or without cleft palate (nsCL/P) cases occur sporadically, suggesting a role for de novo mutations (DNMs) in the etiology of nsCL/P. To identify high impact protein-altering DNMs that contribute to the risk of nsCL/P, we conducted whole-genome sequencing (WGS) analyses in 130 African case-parent trios (affected probands and unaffected parents). We identified 162 high confidence protein-altering DNMs some of which are based on available evidence, contribute to the risk of nsCL/P. These include novel protein-truncating DNMs in the ACTL6A, ARHGAP10, MINK1, TMEM5 and TTN genes; as well as missense variants in ACAN, DHRS3, DLX6, EPHB2, FKBP10, KMT2D, RECQL4, SEMA3C, SEMA4D, SHH, TP63, and TULP4. Many of these protein-altering DNMs were predicted to be pathogenic. Analysis using mouse transcriptomics data showed that some of these genes are expressed during the development of primary and secondary palate. Gene-set enrichment analysis of the protein-altering DNMs identified palatal development and neural crest migration among the few processes that were significantly enriched. These processes are directly involved in the etiopathogenesis of clefting. The analysis of the coding sequence in the WGS data provides more evidence of the opportunity for novel findings in the African genome.

PARD3 gene variation as candidate cause of nonsyndromic cleft palate only

Nonsyndromic cleft palate only (NSCP) is a common congenital malformation worldwide. In this study, we report a three‐generation pedigree with NSCP following the autosomal‐dominant pattern. Whole‐exome sequencing and Sanger sequencing revealed that only the frameshift variant c.1012dupG [p. E338Gfs*26] in PARD3 cosegregated with the disease. In zebrafish embryos, ethmoid plate patterning defects were observed with PARD3 ortholog disruption or expression of patient‐derived N‐terminal truncating PARD3 (c.1012dupG), which implicated PARD3 in ethmoid plate morphogenesis. PARD3 plays vital roles in determining cellular polarity. Compared with the apical distribution of wild‐type PARD3, PARD3‐p. E338Gfs*26 mainly localized to the basal membrane in 3D‐cultured MCF‐10A epithelial cells. The interaction between PARD3‐p. E338Gfs*26 and endogenous PARD3 was identified by LC–MS/MS and validated by co‐IP. Immunofluorescence analysis showed that PARD3‐p. E338Gfs*26 substantially altered the localization of endogenous PARD3 to the basement membrane in 3D‐cultured MCF‐10A cells. Furthermore, seven variants, including one nonsense variant and six missense variants, were identified in the coding region of PARD3 in sporadic cases with NSCP. Subsequent analysis showed that PARD3‐p. R133*, like the insertion variant of c.1012dupG, also changed the localization of endogenous full‐length PARD3 and that its expression induced abnormal ethmoid plate morphogenesis in zebrafish. Based on these data, we reveal PARD3 gene variation as a novel candidate cause of nonsyndromic cleft palate only.

Genetic and Epigenetic Studies in Nonsyndromic Oral Clefts

The etiology of non-syndromic oral clefts (NSOFC) is complex with genetics, genomics, epigenetics and stochastics factors playing a role. Several approaches have been applied to understand the etiology of non-syndromic oral clefts. These include linkage, candidate gene association studies, genome-wide association studies, whole genome sequencing, copy number variations and epigenetics. In this review we shared these approaches, genes and loci reported in some studies.

To Stick or Not to Stick: Adhesions in Orofacial Clefts

Morphogenesis requires a tight coordination between mechanical forces and biochemical signals to inform individual cellular behavior. For these developmental processes to happen correctly the organism requires precise spatial and temporal coordination of the adhesion, migration, growth, differentiation, and apoptosis of cells originating from the three key embryonic layers, namely the ectoderm, mesoderm, and endoderm. The cytoskeleton and its remodeling are essential to organize and amplify many of the signaling pathways required for proper morphogenesis. In particular, the interaction of the cell junctions with the cytoskeleton functions to amplify the behavior of individual cells into collective events that are critical for development. In this review we summarize the key morphogenic events that occur during the formation of the face and the palate, as well as the protein complexes required for cell-to-cell adhesions. We then integrate the current knowledge into a comprehensive review of how mutations in cell-to-cell adhesion genes lead to abnormal craniofacial development, with a particular focus on cleft lip with or without cleft palate.

Generating Zebrafish RNA-Less Mutant Alleles by Deleting Gene Promoters with CRISPR/Cas9

Danio rerio (zebrafish), traditionally used in forward genetic screens, has in the last decade become a popular model for reverse genetic studies with the introduction of TALENS, zinc finger nucleases, and CRISPR/Cas9. Unexpectedly, homozygous frameshift mutations generated by these tools frequently result in phenotypes that are less penetrant than those seen in embryos injected with antisense morpholino oligonucleotides targeting the same gene. One explanation for the difference is that some frameshift mutations result in nonsense-mediated decay of the gene transcript, a process which can induce expression of homologous genes. This form of genetic compensation, called transcriptional adaptation, does not occur when the mutant allele results in no RNA transcripts being produced from the targeted gene. Such RNA-less mutants can be generated by deleting a gene's promoter using a pair of guide RNAs and Cas9 protein. Here, we present a protocol and use it to generate alleles of arhgap29b and slc41a1 that lack detectable zygotic transcription. In the case of the arhgap29b mutant, an emerging phenotype did not segregate with the promoter deletion mutation, highlighting the potential for off-target mutagenesis with these tools. In summary, this chapter describes a method to generate zebrafish mutants that avoid a form of genetic compensation that occurs in many frameshift mutants.

Genetic models and approaches to study orofacial clefts

INTRODUCTION

Orofacial clefts (OFCs) are common craniofacial birth defects with heterogeneous phenotype and etiology. Geneticists have applied nearly every available method and technology to further our understanding of the genetic architectures of OFCs.

OBJECTIVE

This review describes the evidence for a genetic etiology in OFCs, statistical genetic approaches employed to identify genetic causes, and how the results have shaped our current understanding of the genetic architectures of syndromic and nonsyndromic OFCs.

CONCLUSION

There has been rapid progress towards elucidating the genetic architectures of OFCs due to the availability of large collections of DNA samples from cases, controls, and families with OFCs and the consistent adoption of new methodologies and novel statistical approaches as they are developed. Genetic studies have identified rare and common variants influencing risk of OFCs in both Mendelian and complex forms of OFCs, blurring the distinctions traditional categories used in genetic studies and clinical medicine.

Genome-wide Gene-by-Sex Interaction Studies Identify Novel Nonsyndromic Orofacial Clefts Risk Locus

Risk loci identified through genome-wide association studies have explained about 25% of the phenotypic variations in nonsyndromic orofacial clefts (nsOFCs) on the liability scale. Despite the notable sex differences in the incidences of the different cleft types, investigation of loci for sex-specific effects has been understudied. To explore the sex-specific effects in genetic etiology of nsOFCs, we conducted a genome-wide gene × sex (GxSex) interaction study in a sub-Saharan African orofacial cleft cohort. The sample included 1,019 nonsyndromic orofacial cleft cases (814 cleft lip with or without cleft palate and 205 cleft palate only) and 2,159 controls recruited from 3 sites (Ethiopia, Ghana, and Nigeria). An additive logistic model was used to examine the joint effects of the genotype and GxSex interaction. Furthermore, we examined loci with suggestive significance (P < 1E-5) in the additive model for the effect of the GxSex interaction only. We identified a novel risk locus on chromosome 8p22 with genome-wide significant joint and GxSex interaction effects (rs2720555, p_2df = 1.16E-08, p_GxSex = 1.49E-09, odds ratio [OR] = 0.44, 95% CI = 0.34 to 0.57). For males, the risk of cleft lip with or without cleft palate at this locus decreases with additional copies of the minor allele (p < 0.0001, OR = 0.60, 95% CI = 0.48 to 0.74), but the effect is reversed for females (p = 0.0004, OR = 1.36, 95% CI = 1.15 to 1.60). We replicated the female-specific effect of this locus in an independent cohort (p = 0.037, OR = 1.30, 95% CI = 1.02 to 1.65), but no significant effect was found for the males (p = 0.29, OR = 0.86, 95% CI = 0.65 to 1.14). This locus is in topologically associating domain with craniofacially expressed and enriched genes during embryonic development. Rare coding mutations of some of these genes were identified in nsOFC cohorts through whole exome sequencing analysis. Our study is additional proof that genome-wide GxSex interaction analysis provides an opportunity for novel findings of loci and genes that contribute to the risk of nsOFCs.

Pleiotropy method reveals genetic overlap between orofacial clefts at multiple novel loci from GWAS of multi-ethnic trios

Based on epidemiologic and embryologic patterns, nonsyndromic orofacial clefts- the most common craniofacial birth defects in humans- are commonly categorized into cleft lip with or without cleft palate (CL/P) and cleft palate alone (CP), which are traditionally considered to be etiologically distinct. However, some evidence of shared genetic risk in IRF6, GRHL3 and ARHGAP29 regions exists; only FOXE1 has been recognized as significantly associated with both CL/P and CP in genome-wide association studies (GWAS). We used a new statistical approach, PLACO (pleiotropic analysis under composite null), on a combined multi-ethnic GWAS of 2,771 CL/P and 611 CP case-parent trios. At the genome-wide significance threshold of 5 × 10-8, PLACO identified 1 locus in 1q32.2 (IRF6) that appears to increase risk for one OFC subgroup but decrease risk for the other. At a suggestive significance threshold of 10-6, we found 5 more loci with compelling candidate genes having opposite effects on CL/P and CP: 1p36.13 (PAX7), 3q29 (DLG1), 4p13 (LIMCH1), 4q21.1 (SHROOM3) and 17q22 (NOG). Additionally, we replicated the recognized shared locus 9q22.33 (FOXE1), and identified 2 loci in 19p13.12 (RAB8A) and 20q12 (MAFB) that appear to influence risk of both CL/P and CP in the same direction. We found locus-specific effects may vary by racial/ethnic group at these regions of genetic overlap, and failed to find evidence of sex-specific differences. We confirmed shared etiology of the two OFC subtypes comprising CL/P, and additionally found suggestive evidence of differences in their pathogenesis at 2 loci of genetic overlap. Our novel findings include 6 new loci of genetic overlap between CL/P and CP; 3 new loci between pairwise OFC subtypes; and 4 loci not previously implicated in OFCs. Our in-silico validation showed PLACO is robust to subtype-specific effects, and can achieve massive power gains over existing approaches for identifying genetic overlap between disease subtypes. In summary, we found suggestive evidence for new genetic regions and confirmed some recognized OFC genes either exerting shared risk or with opposite effects on risk to OFC subtypes.

Integrative approaches generate insights into the architecture of non-syndromic cleft lip with or without cleft palate

Non-syndromic cleft lip with or without cleft palate (nsCL/P) is a common congenital facial malformation with a multifactorial etiology. Genome-wide association studies (GWASs) have identified multiple genetic risk loci. However, functional interpretation of these loci is hampered by the underrepresentation in public resources of systematic functional maps representative of human embryonic facial development. To generate novel insights into the etiology of nsCL/P, we leveraged published GWAS data on nsCL/P as well as available chromatin modification and expression data on mid-facial development. Our analyses identified five novel risk loci, prioritized candidate target genes within associated regions, and highlighted distinct pathways. Furthermore, the results suggest the presence of distinct regulatory effects of nsCL/P risk variants throughout mid-facial development and shed light on its regulatory architecture. Our integrated data provide a platform to advance hypothesis-driven molecular investigations of nsCL/P and other human facial defects.

Identification of Novel Variants in Cleft Palate-Associated Genes in Brazilian Patients With Non-syndromic Cleft Palate Only

The identification of genetic risk factors for non-syndromic oral clefts is of great importance for better understanding the biological processes related to this heterogeneous and complex group of diseases. Herein we applied whole-exome sequencing to identify potential variants related to non-syndromic cleft palate only (NSCPO) in the multiethnic Brazilian population. Thirty NSCPO samples and 30 sex- and genetic ancestry-matched healthy controls were pooled (3 pools with 10 samples for each group) and subjected to whole-exome sequencing. After filtering, the functional affects, individually and through interactions, of the selected variants and genes were assessed by bioinformatic analyses. As a group, 399 variants in 216 genes related to palatogenesis/cleft palate, corresponding to 6.43%, were exclusively identified in the NSCPO pools. Among those genes are 99 associated with syndromes displaying cleft palate in their clinical spectrum and 92 previously related to cleft lip palate. The most significantly biological processes and pathways overrepresented in the NSCPO-identified genes were associated with the folic acid metabolism, highlighting the interaction between LDL receptor-related protein 6 (LRP6) and 5-methyltetrahydrofolate-homocysteine methyltransferase (MTR) that interconnect two large networks. This study yields novel data on characterization of specific variants and complex processes and pathways related to NSCPO, including many variants in genes of the folate/homocysteine pathway, and confirms that variants in genes related to syndromic cleft palate and cleft lip-palate may cause NSCPO.

Clinically Responsive Genomic Analysis Pipelines: Elements to Improve Detection Rate and Efficiency

Massively parallel sequencing has markedly improved mendelian diagnostic rates. This study assessed the effects of custom alterations to a diagnostic genomic bioinformatic pipeline in response to clinical need and derived practice recommendations relative to diagnostic rates and efficiency. The Genomic Annotation and Interpretation Application (GAIA) bioinformatics pipeline was designed to detect panel, exome, and genome sample integrity and prioritize gene variants in mendelian disorders. Reanalysis of selected negative cases was performed after improvements to the pipeline. GAIA improvements and their effect on sensitivity are described, including addition of a PubMed search for gene-disease associations not in the Online Mendelian Inheritance of Man database, inclusion of a process for calling low-quality variants (known as QPatch), and gene symbol nomenclature consistency checking. The new pipeline increased the diagnostic rate and reduced staff costs, resulting in a saving of US$844.34 per additional diagnosis. Recommendations for genomic analysis pipeline requirements are summarized. Clinically responsive bioinformatics pipeline improvements increase diagnostic sensitivity and increase cost-effectiveness.

Detecting Gene-Environment Interaction for Maternal Exposures Using Case-Parent Trios Ascertained Through a Case With Non-Syndromic Orofacial Cleft

Two large studies of case-parent trios ascertained through a proband with a non-syndromic orofacial cleft (OFC, which includes cleft lip and palate, cleft lip alone, or cleft palate alone) were used to test for possible gene-environment (G × E) interaction between genome-wide markers (both observed and imputed) and self-reported maternal exposure to smoking, alcohol consumption, and multivitamin supplementation during pregnancy. The parent studies were as follows: GENEVA, which included 1,939 case-parent trios recruited largely through treatment centers in Europe, the United States, and Asia, and 1,443 case-parent trios from the Pittsburgh Orofacial Cleft Study (POFC) also ascertained through a proband with an OFC including three major racial/ethnic groups (European, Asian, and Latin American). Exposure rates to these environmental risk factors (maternal smoking, alcohol consumption, and multivitamin supplementation) varied across studies and among racial/ethnic groups, creating substantial differences in power to detect G × E interaction, but the trio design should minimize spurious results due to population stratification. The GENEVA and POFC studies were analyzed separately, and a meta-analysis was conducted across both studies to test for G × E interaction using the 2 df test of gene and G × E interaction and the 1 df test for G × E interaction alone. The 2 df test confirmed effects for several recognized risk genes, suggesting modest G × E effects. This analysis did reveal suggestive evidence for G × Vitamin interaction for CASP9 on 1p36 located about 3 Mb from PAX7, a recognized risk gene. Several regions gave suggestive evidence of G × E interaction in the 1 df test. For example, for G × Smoking interaction, the 1 df test suggested markers in MUSK on 9q31.3 from meta-analysis. Markers near SLCO3A1 also showed suggestive evidence in the 1 df test for G × Alcohol interaction, and rs41117 near RETREG1 (a.k.a. FAM134B) also gave suggestive significance in the meta-analysis of the 1 df test for G × Vitamin interaction. While it remains quite difficult to obtain definitive evidence for G × E interaction in genome-wide studies, perhaps due to small effect sizes of individual genes combined with low exposure rates, this analysis of two large case-parent trio studies argues for considering possible G × E interaction in any comprehensive study of complex and heterogeneous disorders such as OFC.

Identification of a Novel Variant of ARHGAP29 in a Chinese Family with Nonsyndromic Cleft Lip and Palate

Background

Cleft lip with or without cleft palate (CL/P) is the most common facial birth defect, with a worldwide incidence of 1 in 700-1000 live births. CL/P can be divided into syndromic CL/P (SCL/P) and nonsyndromic CL/P (NSCL/P). Genetic factors are an important component to the etiology of NSCL/P. ARHGAP29, one of the NSCL/P disease-causing genes, mediates the cyclical regulation of small GTP binding proteins such as RhoA and plays an essential role in cellular shape, proliferation, and craniofacial development.

Methods

The present study investigated a Chinese family with NSCL/P and explored potential pathogenic variants using whole-exome sequencing (WES). Variants were screened and filtered through bioinformatic analysis and prediction of variant pathogenicity. Cosegregation was subsequently conducted.

Results

We identified a novel heterozygous missense variant of ARHGAP29 (c.2615C > T, p.A872V) in a Chinese pedigree with NSCL/P.

Conclusion

We detected the disease-causing variant in this NSCL/P family. Our identification expands the genetic spectrum of ARHGAP29 and contributes to novel approaches to the genetic diagnosis and counseling of CL/P families.

Non-syndromic Cleft Palate: An Overview on Human Genetic and Environmental Risk Factors

The epithelial and mesenchymal cells involved in early embryonic facial development are guided by complex regulatory mechanisms. Any factor perturbing the growth, approach and fusion of the frontonasal and maxillary processes could result in orofacial clefts that represent the most common craniofacial malformations in humans. The rarest and, probably for this reason, the least studied form of cleft involves only the secondary palate, which is posterior to the incisive foramen. The etiology of cleft palate only is multifactorial and involves both genetic and environmental risk factors. The intention of this review is to give the reader an overview of the efforts made by researchers to shed light on the underlying causes of this birth defect. Most of the scientific papers suggesting potential environmental and genetic causes of non-syndromic cleft palate are summarized in this review, including genome-wide association and gene–environment interaction studies.

A novel splicing mutation of ARHGAP29 is associated with nonsyndromic cleft lip with or without cleft palate

BACKGROUND

Nonsyndromic cleft lip with or without cleft palate (NSCL/P) is one of the most common birth defects, and occurs in approximately 1/700 live births worldwide. The correlation between the ABCA4-ARHGAP29 region and NSCL/P was first identified by genome-wide association studies (GWAS), but few reports have examined NSCL/P caused by ARHGAP29 mutations in the Chinese population.

METHODS

We performed chromosome microarray analysis (CMA) for two consecutive abnormal fetuses and whole exome sequencing (WES) for the family, including 3 patients and 2 normal family members, Sanger sequencing and RT-PCR were used to confirm the mutation.

RESULTS

We identified a novel splice donor mutation (ARHGAP29 c.1920 + 1G > A) in two consecutive NSCL/P fetuses, and the variant was inherited from the mother and grandfather. The mutation caused abnormal skipping of exon 17, and the mRNA level of ARHGAP29 was significantly decreased compared to the wild type.

CONCLUSIONS

In this study, we successfully diagnosed the genetic cause of NSCL/P in a family and first report that the c.1920 + 1G > A mutation in ARHGAP29 is associated with NSCL/P. Our study enriches the genetic landscape of NSCL/P, extends the mutation spectrum of ARHGAP29, and provides a new direction for the diagnosis of NSCL/P in patients and its prenatal diagnosis in fetuses.

SPECC1L regulates palate development downstream of IRF6

SPECC1L mutations have been identified in patients with rare atypical orofacial clefts and with syndromic cleft lip and/or palate (CL/P). These mutations cluster in the second coiled-coil and calponin homology domains of SPECC1L and severely affect the ability of SPECC1L to associate with microtubules. We previously showed that gene-trap knockout of Specc1l in mouse results in early embryonic lethality. We now present a truncation mutant mouse allele, Specc1lΔC510, that results in perinatal lethality. Specc1lΔC510/ΔC510 homozygotes showed abnormal palate rugae but did not show cleft palate. However, when crossed with a gene-trap allele, Specc1lcGT/ΔC510 compound heterozygotes showed a palate elevation delay with incompletely penetrant cleft palate. Specc1lcGT/ΔC510 embryos exhibit transient oral epithelial adhesions at E13.5, which may delay shelf elevation. Consistent with oral adhesions, we show periderm layer abnormalities, including ectopic apical expression of adherens junction markers, similar to Irf6 hypomorphic mutants and Arhgap29 heterozygotes. Indeed, SPECC1L expression is drastically reduced in Irf6 mutant palatal shelves. Finally, we wanted to determine if SPECC1L deficiency also contributed to non-syndromic (ns) CL/P. We sequenced 62 Caucasian, 89 Filipino, 90 Ethiopian, 90 Nigerian and 95 Japanese patients with nsCL/P and identified three rare coding variants (p.Ala86Thr, p.Met91Iso and p.Arg546Gln) in six individuals. These variants reside outside of SPECC1L coiled-coil domains and result in milder functional defects than variants associated with syndromic clefting. Together, our data indicate that palate elevation is sensitive to deficiency of SPECC1L dosage and function and that SPECC1L cytoskeletal protein functions downstream of IRF6 in palatogenesis.

Likely Pathogenic Variants in One Third of Non-Syndromic Discontinuous Cleft Lip and Palate Patients

Oral clefts are composed of cleft of the lip, cleft of the lip and palate, or cleft of the palate, and they are associated with a wide range of expression and severity. When cleft of the palate is associated with cleft of the lip with preservation of the primary palate, it defines an atypical phenotype called discontinuous cleft. Although this phenotype may represent 5% of clefts of the lip and/or palate (CLP), it is rarely specifically referred to and its pathophysiology is unknown. We conducted whole exome sequencing (WES) and apply a candidate gene approach to non-syndromic discontinuous CLP individuals in order to identify genes and deleterious variants that could underlie this phenotype. We discovered loss-of-function variants in two out of the seven individuals, implicating FGFR1 and DLG1 genes, which represents almost one third of this cohort. Whole exome sequencing of clinically well-defined subgroups of CLP, such as discontinuous cleft, is a relevant approach to study CLP etiopathogenesis. It could facilitate more accurate clinical, epidemiological and fundamental research, ultimately resulting in better diagnosis and care of CLP patients. Non-syndromic discontinuous cleft lip and palate seems to have a strong genetic basis.

A systematic genetic analysis and visualization of phenotypic heterogeneity among orofacial cleft GWAS signals

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 subtypes-cleft 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 OFC subtypes. In addition to new associations, we found loci with subtype-specific effects (e.g., GRHL3 [CP], WNT5A [CLP]), as well as loci associated with two or all three subtypes. We cross-referenced these results with mouse craniofacial gene expression datasets, which identified additional promising candidate genes. However, we found no strong correlation between OFC subtypes and expression patterns. In aggregate, the cleft map revealed that neither subtype-specific nor shared genetic effects operate in isolation in OFC architecture. Our approach can be easily applied to any complex trait with distinct phenotypic subgroups.

iSyTE 2.0: a database for expression-based gene discovery in the eye

Although successful in identifying new cataract-linked genes, the previous version of the database iSyTE (integrated Systems Tool for Eye gene discovery) was based on expression information on just three mouse lens stages and was functionally limited to visualization by only UCSC-Genome Browser tracks. To increase its efficacy, here we provide an enhanced iSyTE version 2.0 (URL: http://research.bioinformatics.udel.edu/iSyTE) based on well-curated, comprehensive genome-level lens expression data as a one-stop portal for the effective visualization and analysis of candidate genes in lens development and disease. iSyTE 2.0 includes all publicly available lens Affymetrix and Illumina microarray datasets representing a broad range of embryonic and postnatal stages from wild-type and specific gene-perturbation mouse mutants with eye defects. Further, we developed a new user-friendly web interface for direct access and cogent visualization of the curated expression data, which supports convenient searches and a range of downstream analyses. The utility of these new iSyTE 2.0 features is illustrated through examples of established genes associated with lens development and pathobiology, which serve as tutorials for its application by the end-user. iSyTE 2.0 will facilitate the prioritization of eye development and disease-linked candidate genes in studies involving transcriptomics or next-generation sequencing data, linkage analysis and GWAS approaches.

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.

Hypoxia suppresses myofibroblast differentiation by changing RhoA activity

Fibroblasts show a high range of phenotypic plasticity, including transdifferentiation into myofibroblasts. Myofibroblasts are responsible for generation of the contraction forces that are important for wound healing and scar formation. Overactive myofibroblasts, by contrast, are involved in abnormal scarring. Cell stretching and extracellular signals such as transforming growth factor β can induce the myofibroblastic program, whereas microenvironmental conditions such as reduced tissue oxygenation have an inhibitory effect. We investigated the effects of hypoxia on myofibroblastic properties and linked this to RhoA activity. Hypoxia reversed the myofibroblastic phenotype of primary fibroblasts. This was accompanied by decreased αSMA (ACTA2) expression, alterations in cell contractility, actin reorganization and RhoA activity. We identified a hypoxia-inducible induction of ARHGAP29, which is critically involved in myocardin-related transcription factor-A (MRTF-A) signaling, the differentiation state of myofibroblasts and modulates RhoA activity. This novel link between hypoxia and MRTF-A signaling is likely to be important for ischemia-induced tissue remodeling and the fibrotic response.This article has an associated First Person interview with the first author of the paper.

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.

Zebrafish models of orofacial clefts

Zebrafish is a model organism that affords experimental advantages toward investigating the normal function of genes associated with congenital birth defects. Here we summarize zebrafish studies of genes implicated in orofacial cleft (OFC). The most common use of zebrafish in this context has been to explore the normal function an OFC-associated gene product in craniofacial morphogenesis by inhibiting expression of its zebrafish ortholog. The most frequently deployed method has been to inject embryos with antisense morpholino oligonucleotides targeting the desired transcript. However, improvements in targeted mutagenesis strategies have led to widespread adoption of CRISPR/Cas9 technology. A second application of zebrafish has been for functional assays of gene variants found in OFC patients; such in vivo assays are valuable because the success of in silico methods for testing allele severity has been mixed. Finally, zebrafish have been used to test the tissue specificity of enhancers that harbor single nucleotide polymorphisms associated with risk for OFC. We review examples of each of these approaches in the context of genes that are implicated in syndromic and non-syndromic OFC. Developmental Dynamics 246:897-914, 2017. © 2017 Wiley Periodicals, Inc.

Candidate Genes for Nonsyndromic Cleft Palate Detected by Exome Sequencing

Nonsyndromic cleft palate only (nsCPO) is a facial malformation that has a livebirth prevalence of 1 in 2,500. Research suggests that the etiology of nsCPO is multifactorial, with a clear genetic component. To date, genome-wide association studies have identified only 1 conclusive common variant for nsCPO, that is, a missense variant in the gene grainyhead-like-3 ( GRHL3). Thus, the underlying genetic causes of nsCPO remain largely unknown. The present study aimed at identifying rare variants that might contribute to nsCPO risk, via whole-exome sequencing (WES), in multiply affected Central European nsCPO pedigrees. WES was performed in 2 affected first-degree relatives from each family. Variants shared between both individuals were analyzed for their potential deleterious nature and a low frequency in the general population. Genes carrying promising variants were annotated for 1) reported associations with facial development, 2) multiple occurrence of variants, and 3) expression in mouse embryonic palatal shelves. This strategy resulted in the identification of a set of 26 candidate genes that were resequenced in 132 independent nsCPO cases and 623 independent controls of 2 different ethnicities, using molecular inversion probes. No rare loss-of-function mutation was identified in either WES or resequencing step. However, we identified 2 or more missense variants predicted to be deleterious in each of 3 genes ( ACACB, PTPRS, MIB1) in individuals from independent families. In addition, the analyses identified a novel variant in GRHL3 in 1 patient and a variant in CREBBP in 2 siblings. Both genes underlie different syndromic forms of CPO. A plausible hypothesis is that the apparently nonsyndromic clefts in these 3 patients might represent hypomorphic forms of the respective syndromes. In summary, the present study identified rare variants that might contribute to nsCPO risk and suggests candidate genes for further investigation.

Identification of common non-coding variants at 1p22 that are functional for non-syndromic orofacial clefting

Genome-wide association studies (GWAS) do not distinguish between single nucleotide polymorphisms (SNPs) that are causal and those that are merely in linkage-disequilibrium with causal mutations. Here we describe a versatile, functional pipeline and apply it to SNPs at 1p22, a locus identified in several GWAS for non-syndromic cleft lip with or without cleft palate (NS CL/P). First we amplified DNA elements containing the ten most-highly risk-associated SNPs and tested their enhancer activity in vitro, identifying three SNPs with allele-dependent effects on such activity. We then used in vivo reporter assays to test the tissue-specificity of these enhancers, chromatin configuration capture to test enhancer-promoter interactions, and genome editing in vitro to show allele-specific effects on ARHGAP29 expression and cell migration. Our results further indicate that two SNPs affect binding of CL/P-associated transcription factors, and one affects chromatin configuration. These results translate risk into potential mechanisms of pathogenesis.

Genetics and genomics etiology of nonsyndromic orofacial clefts

Orofacial clefts (OFC) are complex birth defects. Studies using contemporary genomic techniques, bioinformatics, and statistical analyses have led to appreciable advances in identifying the causes of syndromic forms of clefts. This commentary gives an overview of the important cleft gene discoveries found using various genomic methods and tools.