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 (GxS) interaction testing. There were 13 loci significant for GxS interactions, with the top finding in LTBP1 (RR=3.37 [2.04 - 5.56], p=1.93x10 -6 ). LTBP1 plays a role in regulating TGF-B 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 GxS 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.

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.

Rare Variants Analyses Suggest Novel Cleft Genes in the African Population.

Non-syndromic orofacial clefts (NSOFCs) are common birth defects with a complex etiology. While over 60 common risk loci have been identified, they explain only a small proportion of the heritability for NSOFC. Rare variants have been implicated in the missing heritability. Thus, our study aimed to identify genes enriched with nonsynonymous rare coding variants associated with NSOFCs. Our sample included 814 non-syndromic cleft lip with or without palate (NSCL/P), 205 non-syndromic cleft palate only (NSCPO), and 2150 unrelated control children from Nigeria, Ghana, and Ethiopia. We conducted a gene-based analysis separately for each phenotype using three rare-variants collapsing models: (1) protein-altering (PA), (2) missense variants only (MO); and (3) loss of function variants only (LOFO). Subsequently, we utilized relevant transcriptomics data to evaluate associated gene expression and examined their mutation constraint using the gnomeAD database. In total, 13 genes showed suggestive associations (p = E-04). Among them, eight genes (ABCB1, ALKBH8, CENPF, CSAD, EXPH5, PDZD8, SLC16A9, and TTC28) were consistently expressed in relevant mouse and human craniofacial tissues during the formation of the face, and three genes (ABCB1, TTC28, and PDZD8) showed statistically significant mutation constraint. These findings underscore the role of rare variants in identifying candidate genes for NSOFCs.

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.

Parents and Provider Perspectives on the Return of Genomic Findings for Cleft Families in Africa.

BACKGROUND: Inadequate knowledge among health care providers (HCPs) and parents of affected children limits the understanding and utility of secondary genetic findings (SFs) in under-represented populations in genomics research. SFs arise from deep DNA sequencing done for research or diagnostic purposes and may burden patients and their families despite their potential health importance. This study aims to evaluate the perspective of both groups regarding SFs and their choices in the return of results from genetic testing in the context of orofacial clefts. METHODS: Using an online survey, we evaluated the experiences of 252 HCPs and 197 parents across participating cleft clinics in Ghana and Nigeria toward the return of SFs across several domains. RESULTS: Only 1.6% of the HCPs felt they had an expert understanding of when and how to incorporate genomic medicine into practice, while 50.0% agreed that all SFs should be returned to patients. About 95.4% of parents were willing to receive all the information from genetic testing (including SFs), while the majority cited physicians as their primary information source (64%). CONCLUSIONS: Overall, parents and providers were aware that genetic testing could help in the clinical management of diseases. However, they cited a lack of knowledge about genomic medicine, uncertain clinical utility, and lack of available learning resources as barriers. The knowledge gained from this study will assist with developing guidelines and policies to guide providers on the return of SFs in sub-Saharan Africa and across the continent.

Using machine learning algorithms to investigate factors associated with complete edentulism among older adults in the United States.

AIMS: Edentulism is an incapacitating condition, and its prevalence is unequal among different population groups in the United States (US) despite its declining prevalence. This study aimed to investigate the current prevalence, apply Machine Learning (ML) Algorithms to investigate factors associated with complete tooth loss among older US adults, and compare the performance of the models. METHODS: The cross-sectional 2020 Behavioral Risk Factor Surveillance System (BRFSS) data was used to evaluate the prevalence and factors associated with edentulism. ML models were developed to identify factors associated with edentulism utilizing seven ML algorithms. The performance of these models was compared using the area under the receiver operating characteristic curve (AUC). RESULTS: An overall prevalence of 11.9% was reported. The AdaBoost algorithm (AUC = 84.9%) showed the best performance. Analysis showed that the last dental visit, educational attainment, smoking, difficulty walking, and general health status were among the top factors associated with complete edentulism. CONCLUSION: Findings from our study support the declining prevalence of complete edentulism in older adults in the US and show that it is possible to develop a high-performing ML model to investigate the most important factors associated with edentulism using nationally representative data.

Multidisciplinary approach to community engagement in global public health research.

  Recently, there has been a renewed interest in the role of community engagement in knowledge production and ethical issues such as 'helicopter research', indicating exploitative research activities of some researchers as well as short-term relationships with research communities especially in low- and middle-income countries. This approach is detrimental to both communities and the larger scientific community as this may breed mistrust. Major institutions such as the National Institute of Health and Care Research in the United Kingdom have highlighted the importance of community engagement as a tool to improve the reach, quality, and impact of the research by incorporating the voices and concerns of marginalized communities. Similarly, in its 2022 guidance, the American Society for Human Genetics (ASGH) highlights the need to address underrepresentation in genomics research through community engagement. Establishing ethical and meaningful long-term relationships can be challenging especially for researchers who are not members of the community or those from other countries. This article describes how 'community-engaged research' can address some ethical challenges in global public health in different cultural settings.

DNA methylation effects on Van der Woude Syndrome phenotypic variability.

OBJECTIVE: Van der Woude Syndrome (VWS) classically presents with combinations of lip pits (LP) and orofacial clefts, with marked phenotypic discordance even amongst individuals carrying the same mutation. Such discordance suggests a possible role for epigenetic factors as phenotypic modifiers. Both IRF6 , causal for 70% of VWS cases, and TP63 interact in a regulatory loop to coordinate epithelial proliferation and differentiation for palatogenesis. We hypothesize that differential DNA methylation (DNAm) in CpG sites within regulatory regions of IRF6 and TP63 are associated with VWS phenotypic discordance. METHODS: We measured DNAm levels of CpG sites located in the promoter regions of IRF6 and TP63 and in an IRF6 enhancer element (MCS9.7) in 83 individuals with VWS grouped within 5 phenotypes for primary analysis: 1=CL+/-P+LP, 2=CL+/-P, 3=CP+LP, 4=CP, 5=LP and 2 phenotypes for secondary analysis: 1=any cleft and LP, 2= any cleft without LP. DNA samples were bisulfite converted and pyrosequenced with target-specific primers. Methylation levels were compared amongst phenotypes. RESULTS: CpG sites in the IRF6 promoter showed statistically significant differences in methylation among phenotypic groups in both analyses (P<0.05). Individuals with any form of cleft (Groups 1-4) had significantly higher methylation levels than individuals with lip pits only (Group 5). In the secondary analysis, individuals in Group 1 (cleft+LP) had significantly higher methylation than Group 2 (cleft only). CONCLUSION: Results indicated that hypermethylation of the IRF6 promoter is associated with more severe phenotypes (any cleft +/- lip pits); thus, possibly impacting an already genetically weakened IRF6 protein and leading to a more severe phenotype.

Whole Exome Sequencing Identifies Damaging Variants in Indonesians with Clefts.

OBJECTIVES: The interaction between genomics, genetic and environmental factors have been implicated in non-syndromic orofacial cleft development. In the current study, we investigated the contributions of rare and novel genetic variants in known cleft genes using whole exome sequencing (WES) data of Indonesians with non-syndromic orofacial clefts. DESIGN: WES was conducted on 6 individuals. Variants in their exons were called and annotated. These variants were filtered for novelty and rarity using MAF of 0 and 1%. SETTING: Hospital in Indonesia. PATIENTS/PARTICIPANTS: Indonesians with non-syndromic orofacial clefts. INTERVENTIONS: Deleterious variants were prioritized. Pathogenic amino acid changes effect on protein structure and function were identified using HOPE. MAIN OUTCOME MEASURE(S): Rare and novel variants in known cleft genes were filtered and their deleteriousness were predicted using polyphen, SIFT and CADD. RESULTS: We identified rare (MAF <1%) deleterious variants in 4 craniofacial genes namely MMACHC (rs371937044, MAF = 0.00011). SOS1 (rs190222208, MAF = 0.00045), TULP4 (rs199583035, MAF = 0.067), and MTHFD1L (rs143492706, MAF = 0.0044). MMACHC has a mouse knockout model with facial cleft and failure of palatal fusion. The individual with variant in MMACHC presented with nsCPO. CONCLUSIONS: Our study provides additional evidence for the role of TULP4, SOS1, MTHFD1L and MMACHC genes in nsOFC development. This is the first time MMACHC is implicated in nsOFC development in humans.

Trio-based GWAS identifies novel associations and subtype-specific risk factors for cleft palate.

Cleft palate (CP) is one of the most common craniofacial birth defects; however, there are relatively few established genetic risk factors associated with its occurrence despite high heritability. Historically, CP has been studied as a single phenotype, although it manifests across a spectrum of defects involving the hard and/or soft palate. We performed a genome-wide association study using transmission disequilibrium tests of 435 case-parent trios to evaluate broad risks for any cleft palate (ACP) (n = 435), and subtype-specific risks for any cleft soft palate (CSP), (n = 259) and any cleft hard palate (CHP) (n = 125). We identified a single genome-wide significant locus at 9q33.3 (lead SNP rs7035976, p = 4.24 × 10-8) associated with CHP. One gene at this locus, angiopoietin-like 2 (ANGPTL2), plays a role in osteoblast differentiation. It is expressed both in craniofacial tissue of human embryos and developing mouse palatal shelves. We found 19 additional loci reaching suggestive significance (p < 5 × 10-6), of which only one overlapped between groups (chromosome 17q24.2, ACP and CSP). Odds ratios for the 20 loci were most similar across all 3 groups for SNPs associated with the ACP group, but more distinct when comparing SNPs associated with either subtype. We also found nominal evidence of replication (p < 0.05) for 22 SNPs previously associated with orofacial clefts. Our study to evaluate CP risks in the context of its subtypes and we provide newly reported associations affecting the broad risk for CP as well as evidence of subtype-specific risks.

Clinically actionable secondary findings in 130 triads from sub-Saharan African families with non-syndromic orofacial clefts.

INTRODUCTION: The frequency and implications of secondary findings (SFs) from genomic testing data have been extensively researched. However, little is known about the frequency or reporting of SFs in Africans, who are underrepresented in large-scale population genomic studies. The availability of data from the first whole-genome sequencing for orofacial clefts in an African population motivated this investigation. METHODS: In total, 130 case-parent trios were analyzed for SFs within the ACMG SFv.3.0 list genes. Additionally, we filtered for four more genes (HBB, HSD32B, G6PD and ACADM). RESULTS: We identified 246 unique variants in 55 genes; five variants in four genes were classified as pathogenic or likely pathogenic (P/LP). The P/LP variants were seen in 2.3% (9/390) of the subjects, a frequency higher than ~1% reported for diverse ethnicities. On the ACMG list, pathogenic variants were observed in PRKAG (p. Glu183Lys). Variants in the PALB2 (p. Glu159Ter), RYR1 (p. Arg2163Leu) and LDLR (p. Asn564Ser) genes were predicted to be LP. CONCLUSION: This study provides information on the frequency and pathogenicity of SFs in an African cohort. Early risk detection will help reduce disease burden and contribute to efforts to increase knowledge of the distribution and impact of actionable genomic variants in diverse populations.

Trio-based GWAS identifies novel associations and subtype-specific risk factors for cleft palate.

Orofacial clefts (OFCs) are the most common craniofacial birth defects and are often categorized into two etiologically distinct groups: cleft lip with or without cleft palate (CL/P) and isolated cleft palate (CP). CP is highly heritable, but there are still relatively few established genetic risk factors associated with its occurrence compared to CL/P. Historically, CP has been studied as a single phenotype despite manifesting across a spectrum of defects involving the hard and/or soft palate. We performed GWAS using transmission disequilibrium tests using 435 case-parent trios to evaluate broad risks for any cleft palate (ACP, n=435), as well as subtype-specific risks for any cleft soft palate (CSP, n=259) and any cleft hard palate (CHP, n=125). We identified a single genome-wide significant locus at 9q33.3 (lead SNP rs7035976, p=4.24×10 -8 ) associated with CHP. One gene at this locus, angiopoietin-like 2 ( ANGPTL2 ), plays a role in osteoblast differentiation. It is expressed in craniofacial tissue of human embryos, as well as in the developing mouse palatal shelves. We found 19 additional loci reaching suggestive significance (p<5×10 -6 ), of which only one overlapped between groups (chromosome 17q24.2, ACP and CSP). Odds ratios (ORs) for each of the 20 loci were most similar across all three groups for SNPs associated with the ACP group, but more distinct when comparing SNPs associated with either the CSP or CHP groups. We also found nominal evidence of replication (p<0.05) for 22 SNPs previously associated with cleft palate (including CL/P). Interestingly, most SNPs associated with CL/P cases were found to convey the opposite effect in those replicated in our dataset for CP only. Ours is the first study to evaluate CP risks in the context of its subtypes and we provide newly reported associations affecting the broad risk for CP as well as evidence of subtype-specific risks.

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.

Novel IRF6 variant in orofacial cleft patients from Durban, South Africa.

BACKGROUND: To date, there are over 320 variants identified in the IRF6 gene that cause Van der Woude syndrome or popliteal pterygium syndrome. We sequenced this gene in a South African orofacial cleft cohort to identify the causal IRF6 variants in our population. METHOD: Saliva samples from 100 patients with syndromic and non-syndromic CL ± P were collected. Patients were recruited from the cleft clinics at two public, tertiary hospitals in Durban, South Africa (SA), namely Inkosi Albert Luthuli Central Hospital (IALCH) and KwaZulu-Natal Children's Hospital (KZNCH). We prospectively sequenced the exons of IRF6 in 100 orofacial cleft cases, and where possible, we also sequenced the parents of the individuals to determine the segregation pattern. RESULTS: Two variants were identified; one novel (p.Cys114Tyr) and one known (p.Arg84His) missense variant in IRF6 gene were identified. The patient with the p.Cys114Tyr variant was non-syndromic with no clinical VWS phenotype expected of individuals with IRF6 coding variants, and the patient with the p.Arg84His had phenotypic features of popliteal pterygium syndrome. The p.Arg84His variant segregated in the family, with the father also being affected. CONCLUSIONS: This study provides evidence that IRF6 variants are found in the South African population. Genetic counselling is essential for affected families, particularly in the absence of a known clinical phenotype since it helps with the plans for future pregnancies.

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.

Investigating the relationship between cancer and orofacial clefts using GWAS significant loci for cancers: A case-control and case-triad study.

Background: Several population-based case-control studies have reported concurrent presentation of cancer and congenital malformations. Many associations have been made between oral clefting and cancers, though some of these results are conflicting. Some studies have reported an increased risk of cancer among 1st-degree relatives of cleft cases and vice versa, and also an excess risk of cancers of the breast, lung, and brain among those with oral clefts. This study aimed to determine if the genetic polymorphisms found in some cancers are also associated with orofacial cleft in an African cohort. Methods: The study was a case-control and case-triad study in which cases were 400 individuals clinically diagnosed with non-syndromic cleft lip and/or palate (CL/P), while controls were 450 individuals without CL/P. Samples were obtained from three African countries while DNA extraction, PCR, and genotyping were carried out at the University of Iowa, US. Eleven SNPs in genes coding for SWI/SNF subunits and 13 GWAS significant SNPs for cancers associated with orofacial cleft were selected. Case-control analysis, transmission disequilibrium test (TDT), and DFAM to combine the parent-offspring trio data and unrelated case/control data in a single analysis were carried out using PLINK. Results: For the case-control analyses that included all the clefts and for the CLP subtype, none of the SNPs were statistically significant. Statistically increased risk for the following SNPs rs34775372 (p = 0.02; OR = 1.54, CI:1.07-2.22), rs55658222 (p = 0.009; OR = 2.64, CI:1.28-5.45) and rs72728755 (p = 0.02; OR=2.27, CI:1.17-4.45) was observed with the CL only sub-group. None of these were significant after Bonferoni correction. In the TDT analyses, a significantly reduced risk with rs10941679 (p = 0.003; OR = 0.43, CI:0.24-0.75) was observed and this was significant after Bonferroni correction. The rs10941679 was also significant (p = 0.003) in the DFAM analyses as well even after Bonferroni correction. Conclusion: The results from this study represent an important starting point for understanding the concurrent presentation of some cancers in orofacial clefts, and cancer risks in cleft patients. The associations observed warrant further investigation in a larger cohort and will set the stage for a more mechanistic approach toward understanding the risk for cancers in families with clefts.

Dental anomaly detection using intraoral photos via deep learning.

Children with orofacial clefting (OFC) present with a wide range of dental anomalies. Identifying these anomalies is vital to understand their etiology and to discern the complex phenotypic spectrum of OFC. Such anomalies are currently identified using intra-oral exams by dentists, a costly and time-consuming process. We claim that automating the process of anomaly detection using deep neural networks (DNNs) could increase efficiency and provide reliable anomaly detection while potentially increasing the speed of research discovery. This study characterizes the use of` DNNs to identify dental anomalies by training a DNN model using intraoral photographs from the largest international cohort to date of children with nonsyndromic OFC and controls (OFC1). In this project, the intraoral images were submitted to a Convolutional Neural Network model to perform multi-label multi-class classification of 10 dental anomalies. The network predicts whether an individual exhibits any of the 10 anomalies and can do so significantly faster than a human rater can. For all but three anomalies, F1 scores suggest that our model performs competitively at anomaly detection when compared to a dentist with 8 years of clinical experience. In addition, we use saliency maps to provide a post-hoc interpretation for our model's predictions. This enables dentists to examine and verify our model's predictions.

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.