Non-Syndromic Cleft Lip with or without Cleft Palate: Genome-Wide Association Study in Europeans Identifies a Suggestive Risk Locus at 16p12.1 and Supports SH3PXD2A as a Clefting Susceptibility Gene

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.

Gene sequencing applications to combat oral-cavity related disorders: a systematic review with meta-analysis

Gene sequencing (GS) has numerous applications in combatting oral-cavity related disorders, including identifying genetic risk factors for diseases, developing targeted therapies, and improving diagnostic methods. It can help identify specific genetic mutations or variations that increase the risk of developing oral-cavity related disorders, such as oral cancer, periodontal disease, and cleft lip and palate. By the means of the following investigation, our primary objective was to assess the impact of GS technique in diagnosing and potentially treating diseases of the oral cavity by the means of a systematic review and meta-analysis. We commenced by defining the terms "gene sequencing," "oral cavity," and "disorders" as the important elements in our investigation's subject. Next, relevant databases like PubMed, Scopus, Embase, Web of Science, and Google Scholar were searched using keywords and synonyms for each concept, such as "genomic sequencing," "DNA sequencing," "oral health," "oral diseases," "dental caries," "periodontal disease," "oral cancer," and "salivary gland disorders." We combined several search terms, such as "gene sequencing AND oral disorders AND periodontal disease" or "oral cancer OR genomic sequencing," to further hone your search results using Boolean operators like "AND" and "OR." The oral cavity analysis obtained by CS in the selected articles revealed that most of the disorders were, in fact, a direct causal event influenced by the oral microbiome. Moreover, each sampled oral cavity evidenced a different microbial community, which predicted the precipitation of benign as well as malignant conditions, though not on a definitive basis. In the last ten years, genomic sequencing had advanced remarkably as majority of our selected studies observed, making it possible to diagnose and treat a variety of oral and maxillofacial disorders, including cancer. It was also used to ascertain a person's genetic make-up as well as to spot numerous genetic abnormalities that can predispose individuals to diseases. Understanding the different sequencing techniques and the resulting genetic anomalies may help with their clinical application and lead to an improvement in illness diagnosis and prognosis as a whole in the field of dentistry.

Unveiling dysregulated lncRNAs and networks in non-syndromic cleft lip with or without cleft palate pathogenesis

Non-syndromic cleft lip with or without cleft palate (NSCL/P) is a common congenital facial malformation with a complex, incompletely understood origin. Long noncoding RNAs (lncRNAs) have emerged as pivotal regulators of gene expression, potentially shedding light on NSCL/P's etiology. This study aimed to identify critical lncRNAs and construct regulatory networks to unveil NSCL/P's underlying molecular mechanisms. Integrating gene expression profiles from the Gene Expression Omnibus (GEO) database, we pinpointed 30 dysregulated NSCL/P-associated lncRNAs. Subsequent analyses enabled the creation of competing endogenous RNA (ceRNA) networks, lncRNA-RNA binding protein (RBP) interaction networks, and lncRNA cis and trans regulation networks. RT-qPCR was used to examine the regulatory networks of lncRNA in vivo and in vitro. Furthermore, protein levels of lncRNA target genes were validated in human NSCL/P tissue samples and murine palatal shelves. Consequently, two lncRNAs and three mRNAs: FENDRR (log2FC = - 0.671, P = 0.040), TPT1-AS1 (log2FC = 0.854, P = 0.003), EIF3H (log2FC = - 1.081, P = 0.041), RBBP6 (log2FC = 0.914, P = 0.037), and SRSF1 (log2FC = 0.763, P = 0.026) emerged as potential contributors to NSCL/P pathogenesis. Functional enrichment analyses illuminated the biological functions and pathways associated with these lncRNA-related networks in NSCL/P. In summary, this study comprehensively delineates the dysregulated transcriptional landscape, identifies associated lncRNAs, and reveals pivotal sub-networks relevant to NSCL/P development, aiding our understanding of its molecular progression and setting the stage for further exploration of lncRNA and mRNA regulation in NSCL/P.

Ezh2-dependent methylation in oral epithelia promotes secondary palatogenesis

BACKGROUND

In addition to genomic risk variants and environmental influences, increasing evidence suggests epigenetic modifications are important for orofacial development and their alterations can contribute to orofacial clefts. Ezh2 encodes a core catalytic component of the Polycomb repressive complex responsible for addition of methyl marks to Histone H3 as a mechanism of repressing target genes. The role of Ezh2 in orofacial clefts remains unknown.

AIMS

To investigate the epithelial role of Ezh2-dependent methylation in secondary palatogenesis.

METHODS

We used conditional gene-targeting methods to ablate Ezh2 in the surface ectoderm-derived oral epithelium of mouse embryos. We then performed single-cell RNA sequencing combined with immunofluorescence and RT-qPCR to investigate gene expression in conditional mutant palate. We also employed double knockout analyses of Ezh1 and Ezh2 to address if they have synergistic roles in palatogenesis.

RESULTS

We found that conditional inactivation of Ezh2 in oral epithelia results in partially penetrant cleft palate. Double knockout analyses revealed that another family member Ezh1 is dispensable in orofacial development, and it does not have synergistic roles with Ezh2 in palatogenesis. Histochemistry and single-cell RNA-seq analyses revealed dysregulation of cell cycle regulators in the palatal epithelia of Ezh2 mutant mouse embryos disrupts palatogenesis.

CONCLUSION

Ezh2-dependent histone H3K27 methylation represses expression of cell cycle regulator Cdkn1a and promotes proliferation in the epithelium of the developing palatal shelves. Loss of this regulation may perturb movement of the palatal shelves, causing a delay in palate elevation which may result in failure of the secondary palate to close altogether.

Polygenic risk impacts PDGFRA mutation penetrance in non-syndromic cleft lip and palate

Non-syndromic cleft lip with or without cleft palate (NSCL/P) is a common, severe craniofacial malformation that imposes significant medical, psychosocial and financial burdens. NSCL/P is a multifactorial disorder with genetic and environmental factors playing etiologic roles. Currently, only 25% of the genetic variation underlying NSCL/P has been identified by linkage, candidate gene and genome-wide association studies. In this study, whole-genome sequencing and genome-wide genotyping followed by polygenic risk score (PRS) and linkage analyses were used to identify the genetic etiology of NSCL/P in a large three-generation family. We identified a rare missense variant in PDGFRA (c.C2740T; p.R914W) as potentially etiologic in a gene-based association test using pVAAST (P = 1.78 × 10-4) and showed decreased penetrance. PRS analysis suggested that variant penetrance was likely modified by common NSCL/P risk variants, with lower scores found among unaffected carriers. Linkage analysis provided additional support for PRS-modified penetrance, with a 7.4-fold increase in likelihood after conditioning on PRS. Functional characterization experiments showed that the putatively causal variant was null for signaling activity in vitro; further, perturbation of pdgfra in zebrafish embryos resulted in unilateral orofacial clefting. Our findings show that a rare PDGFRA variant, modified by additional common NSCL/P risk variants, have a profound effect on NSCL/P risk. These data provide compelling evidence for multifactorial inheritance long postulated to underlie NSCL/P and may explain some unusual familial patterns.

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.

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.

[Exploring the association between de novo mutations and non-syndromic cleft lip with or without palate based on whole exome sequencing of case-parent trios]

OBJECTIVE

To explore the association between de novo mutations (DNM) and non-syndromic cleft lip with or without palate (NSCL/P) using case-parent trio design.

METHODS

Whole-exome sequencing was conducted for twenty-two NSCL/P trios and Genome Analysis ToolKit (GATK) was used to identify DNM by comparing the alleles of the cases and their parents. Information of predictable functions was annotated to the locus with SnpEff. Enrichment analysis for DNM was conducted to test the difference between the actual number and the expected number of DNM, and to explore whether there were genes with more DNM than expected. NSCL/P-related genes indicated by previous studies with solid evidence were selected by literature reviewing. Protein-protein interactions analysis was conducted among the genes with protein-altering DNM and NSCL/P-related genes. R package "denovolyzeR" was used for the enrichment analysis (Bonferroni correction: P=0.05/n, n is the number of genes in the whole genome range). Protein-protein interactions among genes with DNM and genes with solid evidence on the risk factors of NSCL/P were predicted depending on the information provided by STRING database.

RESULTS

A total of 339 908 SNPs were qualified for the subsequent analysis after quality control. The number of high confident DNM identified by GATK was 345. Among those DNM, forty-four DNM were missense mutations, one DNM was nonsense mutation, two DNM were splicing site mutations, twenty DNM were synonymous mutations and others were located in intron or intergenic regions. The results of enrichment analysis showed that the number of protein-altering DNM on the exome regions was larger than expected (P < 0.05), and five genes (KRTCAP2, HMCN2, ANKRD36C, ADGRL2 and DIPK2A) had more DNM than expected (P < 0.05/(2×19 618)). Protein-protein interaction analysis was conducted among forty-six genes with protein-altering DNM and thirteen genes associated with NSCL/P selected by literature reviewing. Six pairs of interactions occurred between the genes with DNM and known NSCL/P-related genes. The score measuring the confidence level of the predicted interaction between RGPD4 and SUMO1 was 0.868, which was higher than the scores for other pairs of genes.

CONCLUSION

Our study provided novel insights into the development of NSCL/P and demonstrated that functional analyses of genes carrying DNM were warranted to understand the genetic architecture of complex diseases.

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.

Uncovering the Pathogenesis of Orofacial Clefts Using Bioinformatics Analysis

OBJECTIVE

Many genes have been found to be associated with the occurrence of the orofacial clefts (OFC). The links between these pathogenic genes are rarely studied. In this study, bioinformatics analysis were performed in order to find associations between OFC-related genes and provide new ideas for etiology study of OFCs.

METHODS

Orofacial clefts-related genes were searched and identified from the Online Mendelian Inheritance of Man (OMIM.org). These genes were then analyzed by bioinformatics methods, including protein-protein interaction network, functional enrichment analysis, module analysis, and hub genes analysis.

RESULTS

After searching the database of OMIM.org and removing duplicate results, 279 genes were finally obtained. These genes were involved to 369 pathways in biological process, 56 in cell component, 64 in molecular function, and 45 in the Kyoto Encyclopedia of Genes and Genomes. Most identified genes were significantly enriched in embryonic appendage morphogenesis (29.17%), embryonic limb morphogenesis (6.06%), and limb development (4.33%) for biological process (Fig. 5A); ciliary tip (42.86%), MKS complex (28.57%), ciliary basal body (14.29%), and ciliary membrane (14.29%) for cell component. The top 10 hub genes were identified, including SHH, GLI2, PTCH1, SMAD4, FGFR1, BMP4, SOX9, SOX2, RUNX2, and CDH1.

CONCLUSIONS

Bioinformatics methods were used to analyze OFC-related genes in this study, including hub gene identifying and analysis, protein-protein interaction network construction, and functional enrichment analysis. Several potential mechanisms related to occurrence of OFCs were also discussed. These results may be helpful for further studies of the etiology of OFC.

Genome-wide association study in patients with posterior urethral valves

Congenital lower urinary tract obstructions (LUTO) are most often caused by posterior urethral valves (PUV), a male limited anatomical obstruction of the urethra affecting 1 in 4,000 male live births. Little is known about the genetic background of PUV. Here, we report the largest genome-wide association study (GWAS) for PUV in 4 cohorts of patients and controls. The final meta-analysis included 756 patients and 4,823 ethnicity matched controls and comprised 5,754,208 variants that were genotyped or imputed and passed quality control in all 4 cohorts. No genome-wide significant locus was identified, but 33 variants showed suggestive significance (P < 1 × 10^-5). When considering only loci with multiple variants residing within < 10 kB of each other showing suggestive significance and with the same effect direction in all 4 cohorts, 3 loci comprising a total of 9 variants remained. These loci resided on chromosomes 13, 16, and 20. The present GWAS and meta-analysis is the largest genetic study on PUV performed to date. The fact that no genome-wide significant locus was identified, can be explained by lack of power or may indicate that common variants do not play a major role in the etiology of PUV. Nevertheless, future studies are warranted to replicate and validate the 3 loci that yielded suggestive associations.

Association between PTCH1 and RAD54B Single-Nucleotide Polymorphisms and Non-syndromic Orofacial Clefts in the Northeast Population of Iran

Background

Non-Syndromic Cleft Lip with or without cleft Palate (NSCL/P) is a common developmental disorder of the head and neck with a multifactorial etiology. The current study aimed to evaluate the potential association of PTCH1 (rs10512248) and RAD54B (rs12681366) polymorphisms with NSCL/P in the Northeast Iranian population.

Methods

In the present study, blood samples were taken from 122 subjects with NSCL/P and 161 healthy controls. Polymerase Chain Reaction (PCR) followed by Restriction Fragment Length Polymorphism (RFLP) were used to conduct genotyping of single-nucleotide polymorphisms.

Results

Although differences were observed between cases and controls in rs10512248 and rs12681366, our data did not support a significant association of these polymorphisms with NSCL/P in our population.

Conclusion

Our findings suggest that polymorphisms of rs10512248 and rs12681366 may not be potential risk factors for NSCL/P in the Northeast Iranian population due to the multifactorial and multiethnicity characteristics of some genes.

Genetic architecture of orbital telorism

The interocular distance, or orbital telorism, is a distinctive craniofacial trait that also serves as a clinically informative measure. While its extremes, hypo- and hypertelorism, have been linked to monogenic disorders and are often syndromic, little is known about the genetic determinants of interocular distance within the general population. We derived orbital telorism measures from cranial magnetic resonance imaging by calculating the distance between the eyeballs’ centre of gravity, which showed a good reproducibility with an intraclass correlation coefficient of 0.991 (95% confidence interval 0.985–0.994). Heritability estimates were 76% (standard error = 12%) with a family-based method (N = 364) and 39% (standard error = 2.4%) with a single nucleotide polymorphism-based method (N = 34 130) and were unaffected by adjustment for height (model II) and intracranial volume (model III) or head width (model IV). Genome-wide association studies in 34 130 European individuals identified 56 significantly associated genomic loci (P < 5 × 10⁻⁸) across four different models of which 46 were novel for facial morphology, and overall these findings replicated in an independent sample (N = 10 115) with telorism-related horizontal facial distance measures. Genes located nearby these 56 identified genetic loci were 4.9-fold enriched for Mendelian hypotelorism and hypertelorism genes, underlining their biological relevance. This study provides novel insights into the genetic architecture underlying interocular distance in particular, and the face in general, and explores its potential for applications in a clinical setting.

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.

Genome-Wide Association Study of Non-syndromic Orofacial Clefts in a Multiethnic Sample of Families and Controls Identifies Novel Regions

Orofacial clefts (OFCs) are among the most prevalent craniofacial birth defects worldwide and create a significant public health burden. The majority of OFCs are non-syndromic and vary in prevalence by ethnicity. Africans have the lowest prevalence of OFCs (~ 1/2,500), Asians have the highest prevalence (~1/500), Europeans and Latin Americans lie somewhere in the middle (~1/800 and 1/900, respectively). Thus, ethnicity appears to be a major determinant of the risk of developing OFC. The Pittsburgh Orofacial Clefts Multiethnic study was designed to explore this ethnic variance, comprising a large number of families and individuals (~12,000 individuals) from multiple populations worldwide: US and Europe, Asians, mixed Native American/Caucasians, and Africans. In this current study, we analyzed 2,915 OFC cases, 6,044 unaffected individuals related to the OFC cases, and 2,685 controls with no personal or family history of OFC. Participants were grouped by their ancestry into African, Asian, European, and Central and South American subsets, and genome-wide association run on the combined sample as well as the four ancestry-based groups. We observed 22 associations to cleft lip with or without cleft palate at 18 distinct loci with p-values < 1e-06, including 10 with genome-wide significance (<5e-08), in the combined sample and within ancestry groups. Three loci - 2p12 (rs62164740, p = 6.27e-07), 10q22.2 (rs150952246, p = 3.14e-07), and 10q24.32 (rs118107597, p = 8.21e-07) are novel. Nine were in or near known OFC loci - PAX7, IRF6, FAM49A, DCAF4L2, 8q24.21, NTN1, WNT3-WNT9B, TANC2, and RHPN2. The majority of the associations were observed only in the combined sample, European, and Central and South American groups. We investigated whether the observed differences in association strength were (a) purely due to sample sizes, (b) due to systematic allele frequency difference at the population level, or (c) due to the fact certain OFC-causing variants confer different amounts of risk depending on ancestral origin, by comparing effect sizes to observed allele frequencies of the effect allele in our ancestry-based groups. While some of the associations differ due to systematic differences in allele frequencies between groups, others show variation in effect size despite similar frequencies across ancestry groups.

Co-occurrence of orofacial clefts and clubfoot phenotypes in a sub-Saharan African cohort: Whole-exome sequencing implicates multiple syndromes and genes

Background

Orofacial clefts (OFCs) are congenital malformations of the face and palate, with an incidence of 1 per 700 live births. Clubfoot or congenital talipes equinovarus (CTEV) is a three‐dimensional abnormality of the leg, ankle, and feet that leads to the anomalous positioning of foot and ankle joints and has an incidence of 1 per 1000 live births. OFCs and CTEV may occur together or separately in certain genetic syndromes in addition to other congenital abnormalities. Here, we sought to decipher the genetic etiology of OFC and CTEV that occurred together in six probands.

Methods

At the time of recruitment, the most clinically obvious congenital anomalies in these individuals were the OFC and CTEV. We carried out whole‐exome sequencing (WES) on DNA samples from probands and available parents employing the Agilent SureSelect XT kit and Illumina HiSeq2500 platform, followed by bioinformatics analyses. WES variants were validated by clinical Sanger Sequencing.

Results

Of the six probands, we observed probable pathogenic genetic variants in four. In three probands with probable pathogenic genetic variants, each individual had variants in three different genes, whereas one proband had probable pathogenic variant in just one gene. In one proband, we observed variants in DIS3L2, a gene associated with Perlman syndrome. A second proband had variants in EPG5 (associated with Vici Syndrome), BARX1 and MKI67, while another proband had potentially etiologic variants in FRAS1 (associated with Fraser Syndrome 1), TCOF1 (associated with Treacher Collins Syndrome 1) and MKI67. The last proband had variants in FRAS1, PRDM16 (associated with Cardiomyopathy, dilated, 1LL/Left ventricular noncompaction 8) and CHD7 (associated with CHARGE syndrome/Hypogonadotropic hypogonadism 5 with or without anosmia).

Conclusion

Our results suggest that clubfoot and OFCs are two congenital abnormalities that can co‐occur in certain individuals with varying genetic causes and expressivity, warranting the need for deep phenotyping.

Orofacial clefts embryology, classification, epidemiology, and genetics

Orofacial clefts (OFCs) rank as the second most common congenital birth defect in the United States after Down syndrome and are the most common head and neck congenital malformations. They are classified as cleft lip with or without cleft palate (CL/P) and cleft palate only (CPO). OFCs have significant psychological and socio-economic impact on patients and their families and require a multidisciplinary approach for management and counseling. A complex interaction between genetic and environmental factors contributes to the incidence and clinical presentation of OFCs. In this comprehensive review, the embryology, classification, epidemiology and etiology of clefts are thoroughly discussed and a "state-of-the-art" snapshot of the recent advances in the genetics of OFCs is presented.

Genetics and signaling mechanisms of orofacial clefts

Craniofacial development involves several complex tissue movements including several fusion processes to form the frontonasal and maxillary structures, including the upper lip and palate. Each of these movements are controlled by many different factors that are tightly regulated by several integral morphogenetic signaling pathways. Subject to both genetic and environmental influences, interruption at nearly any stage can disrupt lip, nasal, or palate fusion and result in a cleft. Here, we discuss many of the genetic risk factors that may contribute to the presentation of orofacial clefts in patients, and several of the key signaling pathways and underlying cellular mechanisms that control lip and palate formation, as identified primarily through investigating equivalent processes in animal models, are examined.