Genome-wide analysis of copy-number variation in humans with cleft lip and/or cleft palate identifies COBLL1, RIC1, and ARHGEF38 as clefting genes.

Cleft lip with or without cleft palate (CL/P) is a common birth defect with a complex, heterogeneous etiology. It is well established that common and rare sequence variants contribute to the formation of CL/P, but the contribution of copy-number variants (CNVs) to cleft formation remains relatively understudied. To fill this knowledge gap, we conducted a large-scale comparative analysis of genome-wide CNV profiles of 869 individuals from the Philippines and 233 individuals of European ancestry with CL/P with three primary goals: first, to evaluate whether differences in CNV number, amount of genomic content, or amount of coding genomic content existed within clefting subtypes; second, to assess whether CNVs in our cohort overlapped with known Mendelian clefting loci; and third, to identify unestablished Mendelian clefting genes. Significant differences in CNVs across cleft types or in individuals with non-syndromic versus syndromic clefts were not observed; however, several CNVs in our cohort overlapped with known syndromic and non-syndromic Mendelian clefting loci. Moreover, employing a filtering strategy relying on population genetics data that rare variants are on the whole more deleterious than common variants, we identify several CNV-associated gene losses likely driving non-syndromic clefting phenotypes. By prioritizing genes deleted at a rare frequency across multiple individuals with clefts yet enriched in our cohort of individuals with clefts compared to control subjects, we identify COBLL1, RIC1, and ARHGEF38 as clefting genes. CRISPR-Cas9 mutagenesis of these genes in Xenopus laevis and Danio rerio yielded craniofacial dysmorphologies, including clefts analogous to those seen in human clefting disorders.

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.

Genome-wide study of gene-by-sex interactions identifies risks for cleft palate.

Structural birth defects affect 3-4% of all live births and, depending on the type, tend to manifest in a sex-biased manner. Orofacial clefts (OFCs) are the most common craniofacial structural birth defects and are often divided into cleft lip with or without cleft palate (CL/P) and cleft palate only (CP). Previous studies have found sex-specific risks for CL/P, but these risks have yet to be evaluated in CP. CL/P is more common in males and CP is more frequently observed in females, so we hypothesized there would also be sex-specific differences for CP. Using a trio-based cohort, we performed sex-stratified genome-wide association studies (GWAS) based on proband sex followed by a genome-wide gene-by-sex (G × S) interaction testing. There were 13 loci significant for G × S interactions, with the top finding in LTBP1 (RR = 3.37 [2.04-5.56], p = 1.93 × 10-6). LTBP1 plays a role in regulating TGF-ß bioavailability, and knockdown in both mice and zebrafish lead to craniofacial anomalies. Further, there is evidence for differential expression of LTBP1 between males and females in both mice and humans. Therefore, we tested the association between the imputed genetically regulated gene expression of genes with significant G × S interactions and the CP phenotype. We found significant association for LTBP1 in cell cultured fibroblasts in female probands (p = 0.0013) but not in males. Taken altogether, we show there are sex-specific risks for CP that are otherwise undetectable in a combined sex cohort, and LTBP1 is a candidate risk gene, particularly in females.

Axenfeld-Rieger syndrome: more than meets the eye.

BACKGROUND: Axenfeld-Rieger syndrome (ARS) is characterised by typical anterior segment anomalies, with or without systemic features. The discovery of causative genes identified ARS subtypes with distinct phenotypes, but our understanding is incomplete, complicated by the rarity of the condition. METHODS: Genetic and phenotypic characterisation of the largest reported ARS cohort through comprehensive genetic and clinical data analyses. RESULTS: 128 individuals with causative variants in PITX2 or FOXC1, including 81 new cases, were investigated. Ocular anomalies showed significant overlap but with broader variability and earlier onset of glaucoma for FOXC1-related ARS. Systemic anomalies were seen in all individuals with PITX2-related ARS and the majority of those with FOXC1-related ARS. PITX2-related ARS demonstrated typical umbilical anomalies and dental microdontia/hypodontia/oligodontia, along with a novel high rate of Meckel diverticulum. FOXC1-related ARS exhibited characteristic hearing loss and congenital heart defects as well as previously unrecognised phenotypes of dental enamel hypoplasia and/or crowding, a range of skeletal and joint anomalies, hypotonia/early delay and feeding disorders with structural oesophageal anomalies in some. Brain imaging revealed highly penetrant white matter hyperintensities, colpocephaly/ventriculomegaly and frequent arachnoid cysts. The expanded phenotype of FOXC1-related ARS identified here was found to fully overlap features of De Hauwere syndrome. The results were used to generate gene-specific management plans for the two types of ARS. CONCLUSION: Since clinical features of ARS vary significantly based on the affected gene, it is critical that families are provided with a gene-specific diagnosis, PITX2-related ARS or FOXC1-related ARS. De Hauwere syndrome is proposed to be a FOXC1opathy.

Perceptions and beliefs of community gatekeepers about genomic risk information in African cleft research.

BACKGROUND: A fundamental ethical issue in African genomics research is how socio-cultural factors impact perspectives, acceptance, and utility of genomic information, especially in stigmatizing conditions like orofacial clefts (OFCs). Previous research has shown that gatekeepers (e.g., religious, political, family or community leaders) wield considerable influence on the decision-making capabilities of their members, including health issues. Thus, their perspectives can inform the design of engagement strategies and increase exposure to the benefits of genomics testing/research. This is especially important for Africans underrepresented in genomic research. Our study aims to investigate the perspectives of gatekeepers concerning genomic risk information (GRI) in the presence of OFCs in a sub-Saharan African cohort. METHODS: Twenty-five focus group discussions (FGDs) consisting of 214 gatekeepers (religious, community, ethnic leaders, and traditional birth attendants) in Lagos, Nigeria, explored the opinions of participants on genomic risk information (GRI), OFC experience, and the possibility of involvement in collaborative decision-making in Lagos, Nigeria. Transcripts generated from audio recordings were coded and analyzed in NVivo using thematic analysis. RESULTS: Three main themes-knowledge, beliefs, and willingness to act-emerged from exploring the perspective of gatekeepers about GRI in this group. We observed mixed opinions regarding the acceptance of GRI. Many participants believed their role is to guide and support members when they receive results; this is based on the level of trust their members have in them. However, participants felt they would need to be trained by medical experts to do this. Also, religious and cultural beliefs were crucial to determining participants' understanding of OFCs and the acceptance and utilization of GRI. CONCLUSIONS: Incorporating cultural sensitivity into public engagement could help develop appropriate strategies to manage conflicting ideologies surrounding genomic information in African communities. This will allow for more widespread access to the advances in genomics research in underrepresented populations. We also recommend a synergistic relationship between community health specialists/scientists, and community leaders, including spiritual providers to better understand and utilize GRI.

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.

Genome-wide association study of multiethnic nonsyndromic orofacial cleft families identifies novel loci specific to family and phenotypic subtypes.

Nonsyndromic orofacial clefts (OFCs) are among the most common craniofacial birth defects worldwide, and known to exhibit phenotypic and genetic heterogeneity. Cleft lip plus cleft palate (CLP) and cleft lip only (CL) are commonly combined together as one phenotype (CL/P), separately from cleft palate alone. In comparison, our study analyzes CL and CLP separately. A sample of 2218 CL and CLP cases, 4537 unaffected relatives of cases, and 2673 pure controls with no family history of OFC were selected from the Pittsburgh Orofacial Cleft (Pitt-OFC) multiethnic study.genome-wide association studies were run for seven specific phenotypes created based on the cleft type(s) observed within these families, as well as the combined CL/P phenotype. Five novel genome-wide significant associations, 3q29 (rs62284390), 5p13.2 (rs609659), 7q22.1 (rs6465810), 19p13.3 (rs628271), and 20q13.33 (rs2427238), and nine associations (p ≤ 1.0E-05) within previously confirmed OFC loci-PAX7, IRF6, FAM49A, DCAF4L2, 8q24.21, ARID3B, NTN1, TANC2 and the WNT9B:WNT3 gene cluster-were observed. We also found that single nucleotide polymorphisms within a subset of the associated loci, both previously known and novel, differ substantially in terms of their effects across cleft- or family-specific phenotypes, indicating not only etiologic differences between CL and CLP, but also genetic heterogeneity within each of the two OFC subtypes.

Rare variant modifier analysis identifies variants in SEC24D associated with orofacial cleft subtypes.

As one of the most common structural birth defects, orofacial clefts (OFCs) have been studied for decades, and recent studies have demonstrated that there are genetic differences between the different phenotypic presentations of OFCs. However, the contribution of rare genetic variation genome-wide to different subtypes of OFCs has been understudied, with most studies focusing on common genetic variation or rare variation within targeted regions of the genome. Therefore, we used whole-genome sequencing data from the Gabriella Miller Kids First Pediatric Research Program to conduct a gene-based burden analysis to test for genetic modifiers of cleft lip (CL) vs cleft lip and palate (CLP). We found that there was a significantly increased burden of rare variants in SEC24D in CL cases compared to CLP cases (p = 6.86 [Formula: see text] 10-7). Of the 15 variants within SEC24D, 53.3% were synonymous, but overlapped a known craniofacial enhancer. We then tested whether these variants could alter predicted transcription factor binding sites (TFBS), and found that the rare alleles destroyed binding sites for 9 transcription factors (TFs), including Pax1 (p = 0.0009), and created binding sites for 23 TFs, including Pax6 (p = 6.12 [Formula: see text] 10-5) and Pax9 (p = 0.0001), which are known to be involved in normal craniofacial development, suggesting a potential mechanism by which these synonymous variants could have a functional impact. Overall, this study indicates that rare genetic variation may contribute to the phenotypic heterogeneity of OFCs and suggests that regulatory variation may also contribute and warrant further investigation in future studies of genetic variants controlling risk to OFC.

Genome-wide Enrichment of De Novo Coding Mutations in Orofacial Cleft Trios.

Although de novo mutations (DNMs) are known to increase an individual's risk of congenital defects, DNMs have not been fully explored regarding orofacial clefts (OFCs), one of the most common human birth defects. Therefore, whole-genome sequencing of 756 child-parent trios of European, Colombian, and Taiwanese ancestry was performed to determine the contributions of coding DNMs to an individual's OFC risk. Overall, we identified a significant excess of loss-of-function DNMs in genes highly expressed in craniofacial tissues, as well as genes associated with known autosomal dominant OFC syndromes. This analysis also revealed roles for zinc-finger homeobox domain and SOX2-interacting genes in OFC etiology.

Rare variants found in clinical gene panels illuminate the genetic and allelic architecture of orofacial clefting.

PURPOSE: Orofacial clefts (OFCs) are common birth defects including cleft lip, cleft lip and palate, and cleft palate. OFCs have heterogeneous etiologies, complicating clinical diagnostics because it is not always apparent if the cause is Mendelian, environmental, or multifactorial. Sequencing is not currently performed for isolated or sporadic OFCs; therefore, we estimated the diagnostic yield for 418 genes in 841 cases and 294 controls. METHODS: We evaluated 418 genes using genome sequencing and curated variants to assess their pathogenicity using American College of Medical Genetics criteria. RESULTS: 9.04% of cases and 1.02% of controls had "likely pathogenic" variants (P < .0001), which was almost exclusively driven by heterozygous variants in autosomal genes. Cleft palate (17.6%) and cleft lip and palate (9.09%) cases had the highest yield, whereas cleft lip cases had a 2.80% yield. Out of 39 genes with likely pathogenic variants, 9 genes, including CTNND1 and IRF6, accounted for more than half of the yield (4.64% of cases). Most variants (61.8%) were "variants of uncertain significance", occurring more frequently in cases (P = .004), but no individual gene showed a significant excess of variants of uncertain significance. CONCLUSION: These results underscore the etiological heterogeneity of OFCs and suggest sequencing could reduce the diagnostic gap in OFCs.

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.

A mutational hotspot in AMOTL1 defines a new syndrome of orofacial clefting, cardiac anomalies, and tall stature.

AMOTL1 encodes angiomotin-like protein 1, an actin-binding protein that regulates cell polarity, adhesion, and migration. The role of AMOTL1 in human disease is equivocal. We report a large cohort of individuals harboring heterozygous AMOTL1 variants and define a core phenotype of orofacial clefting, congenital heart disease, tall stature, auricular anomalies, and gastrointestinal manifestations in individuals with variants in AMOTL1 affecting amino acids 157-161, a functionally undefined but highly conserved region. Three individuals with AMOTL1 variants outside this region are also described who had variable presentations with orofacial clefting and multi-organ disease. Our case cohort suggests that heterozygous missense variants in AMOTL1, most commonly affecting amino acid residues 157-161, define a new orofacial clefting syndrome, and indicates an important functional role for this undefined region.

FAT4 identified as a potential modifier of orofacial cleft laterality.

Orofacial clefts (OFCs) are common (1 in 700 births) congenital malformations that include a cleft lip (CL) and cleft lip and palate (CLP). These OFC subtypes are also heterogeneous themselves, with the CL occurring on the left, right, or both sides of the upper lip. Unilateral CL and CLP have a 2:1 bias towards left-sided clefts, suggesting a nonrandom process. Here, we performed a study of left- and right-sided clefts within the CL and CLP subtypes to better understand the genetic factors controlling cleft laterality. We conducted genome-wide modifier analyses by comparing cases that had right unilateral CL (RCL; N = 130), left unilateral CL (LCL; N = 216), right unilateral CLP (RCLP; N = 416), or left unilateral CLP (LCLP; N = 638), and identified a candidate region on 4q28, 400 kb downstream from FAT4, that approached genome-wide significance for LCL versus RCL (p = 8.4 × 10-8 ). Consistent with its potential involvement as a genetic modifier of CL, we found that Fat4 exhibits a specific domain of expression in the mesenchyme of the medial nasal processes that form the median upper lip. Overall, these results suggest that the epidemiological similarities in left- to right-sided clefts in CL and CLP are not reflected in the genetic association results.

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.

Interactions between PDA-associated polymorphisms and genetic ancestry alter ductus arteriosus gene expression.

BACKGROUND: DNA polymorphisms in PTGIS and TFAP2B have been identified as risk factors for patent ductus arteriosus (PDA) in a population composed of preterm infants with European genetic ancestry but not in more genetically diverse populations. GOAL: To determine if the effects of TFAP2B and PTGIS polymorphisms on ductus arteriosus (DA) gene expression differ based on genetic ancestry. METHODS: DA from 273 human second trimester fetuses were genotyped for TFAP2B and PTGIS polymorphisms and for polymorphisms distributing along genetic ancestry lines. RT-PCR was used to measure the RNA expression of 49 candidate genes involved with DA closure. RESULTS: Seventeen percent of the DA analyzed were of European ancestry. In multivariable regression analyses we found consistent associations between four PDA-related TFAP2B polymorphisms (rs2817399(A), rs987237(G), rs760900(C), and rs2817416(C)) and expression of the following genes: EPAS1, CACNB2, ECE1, KCNA2, ATP2A3, EDNRA, EDNRB, BMP9, and BMP10, and between the PTGIS haplotype rs493694(G)/rs693649(A) and PTGIS and NOS3. These changes only occurred in DA with European ancestry. No consistent positive or negative associations were found among DA samples unless an interaction between the polymorphisms and genetic ancestry was taken into account. CONCLUSION: PTGIS and TFAP2B polymorphisms were associated with consistent changes in DA gene expression when present in fetuses with European ancestry. IMPACT: DNA polymorphisms in PTGIS and TFAP2B have been identified as risk factors for patent ductus arteriosus (PDA) in a population composed primarily of preterm infants with European genetic ancestry but not in more genetically diverse populations. The same PTGIS and TFAP2B polymorphisms are associated with changes in ductus gene expression when present in ductus from fetuses with European genetic ancestry. No consistent associations with gene expression can be found unless an interaction between the polymorphisms and genetic ancestry is taken into account.

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.

Genome-Wide Scan for Parent-of-Origin Effects in a sub-Saharan African Cohort With Nonsyndromic Cleft Lip and/or Cleft Palate (CL/P).

OBJECTIVE: Nonsyndromic cleft lip and/or cleft palate (NSCL/P) have multifactorial etiology where genetic factors, gene-environment interactions, stochastic factors, gene-gene interactions, and parent-of-origin effects (POEs) play cardinal roles. POEs arise when the parental origin of alleles differentially impacts the phenotype of the offspring. The aim of this study was to identify POEs that can increase risk for NSCL/P in humans using a genome-wide dataset. METHODS: The samples (174 case-parent trios from Ghana, Ethiopia, and Nigeria) included in this study were from the African only genome wide association studies (GWAS) that was published in 2019. Genotyping of individual DNA using over 2 million multiethnic and African ancestry-specific single-nucleotide polymorphisms from the Illumina Multi-Ethnic Genotyping Array v2 15070954 A2 (genome build GRCh37/hg19) was done at the Center for Inherited Diseases Research. After quality control checks, PLINK was employed to carry out POE analysis employing the pooled subphenotypes of NSCL/P. RESULTS: We observed possible hints of POEs at a cluster of genes at a 1 mega base pair window at the major histocompatibility complex class 1 locus on chromosome 6, as well as at other loci encompassing candidate genes such as ASB18, ANKEF1, AGAP1, GABRD, HHAT, CCT7, DNMT3A, EPHA7, FOXO3, lncRNAs, microRNA, antisense RNAs, ZNRD1, ZFAT, and ZBTB16. CONCLUSION: Findings from our study suggest that some loci may increase the risk for NSCL/P through POEs. Additional studies are required to confirm these suggestive loci in NSCL/P etiology.

Damaging Mutations in AFDN 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.

Using an aquatic model, Xenopus laevis, to uncover the role of chromodomain 1 in craniofacial disorders.

The chromodomain family member chromodomain 1 (CHD1) has been shown to have numerous critical molecular functions including transcriptional regulation, splicing, and DNA repair. Complete loss of function of this gene is not compatible with life. On the other hand, missense and copy number variants of CHD1 can result in intellectual disabilities and craniofacial malformations in human patients including cleft palate and Pilarowski-Bjornsson Syndrome. We have used the aquatic developmental model organism Xenopus laevis, to determine a specific role for Chd1 in such cranioafcial disorders. Protein and gene knockdown techniques in Xenopus, including antisense oligos and mosaic Crispr/Cas9-mediated mutagenesis, recapitulated the craniofacial defects observed in humans. Further analysis indicated that embryos deficient in Chd1 had defects in cranial neural crest development and jaw cartilage morphology. Additionally, flow cytometry and immunohistochemistry revealed that decreased Chd1 resulted in increased in apoptosis in the developing head. Together, these experiments demonstrate that Chd1 is critical for fundamental processes and cell survival in craniofacial development. We also presented evidence that Chd1 is regulated by retinoic acid signaling during craniofacial development. Expression levels of chd1 mRNA, specifically in the head, were increased by RAR agonist exposure and decreased upon antagonist treatment. Subphenotypic levels of an RAR antagonist and Chd1 morpholinos synergized to result in orofacial defects. Further, RAR DNA binding sequences (RAREs) were detected in chd1 regulatory regions by bioinformatic analysis. In summary, by combining human genetics and experiments in an aquatic model we now have a better understanding of the role of CHD1 in craniofacial disorders.

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.