OBJECTIVES: This study aims to explore the association between single nucleotide polymorphisms (SNPs) loci near the haplotype region hg19 chr9:100560865-100660865 of the forkhead box E1 (FOXE1) gene and the occurrence of non-syndromic cleft lip with or without cleft palate (NSCL/P) in western Han Chinese population. METHODS: In the first stage, our study recruited 159 NSCL/P patients and performed targeted region sequencing to screen SNPs loci near the haplotype region of the FOXE1 gene associated with NSCL/P. In the second stage, we selected 21 common SNPs and re-enrolled 1 000 non-syndromic cleft lip only (NSCLO) patients, 1 000 non-syndromic cleft palate only (NSCPO) patients, and 1 000 normal controls to verify the association. PLINK software was used to perform Hardy-Weinberg equilibrium (HWE) test. Association analysis for common variants, gene burden analysis for rare mutations, and function prediction of SNPs with non-synonymous mutations were performed using Mutation Taster and other software programs. RESULTS: In the first stage, 126 variants, including 76 single nucleotide variants and 50 insertion-deletions were identified. All the included SNPs confirmed to HWE, and the results of gene burden analysis and prediction of functional harmfulness for rare variants were not statistically significant. Association analysis showed that rs13292899 of the FOXE1 gene was significantly associated with NSCL/P (P=1.85E-27) and was also correlated with NSCLO (P=6.41E-23) and non-syndromic cleft lip with cleft palate (NSCLP) (P=2.36E-15) subtypes. In the validation phase, rs79268293 (P=0.013, P=0.022), rs10983951 (P=0.009 2, P=0.007 6), rs117227387 (P=0.009 2, P=0.007 6), rs3758250 (P=0.009 2, P=0.007 6), and rs116899397 (P=0.009 2, P=0.007 6) were significantly associated with NSCLO and NSCPO; rs13292899 (P=0.008 5), rs74606599 (P=0.008 3), rs143226042 (P=0.008 3), and rs117236550 (P=0.01) were associated with the occurrence of NSCLO; and rs12343182 (P=0.008 7), rs10119760 (P=0.012), rs10113907 (P=0.012), and rs13299924 (P=0.012) were associated with the occurrence of NSCPO. CONCLUSIONS This study found a new susceptible SNP rs13292899 of the FOXE1 gene that is closely associated with NSCL/P and NSCLO subtype and 13 other SNPs associated with NSCLO or NSCPO.
Female ; Humans ; Male ; China ; Cleft Lip/genetics* ; Cleft Palate/genetics* ; Forkhead Transcription Factors/genetics* ; Haplotypes ; Polymorphism, Single Nucleotide ; East Asian People/genetics*

Muenke syndrome is an autosomal dominant genetic disorder that is typically characterized by unilateral or bilateral coronal synostosis, macrocephaly, midface hypoplasia, and developmental delays. This article reports a case of Muenke syndrome with a soft cleft palate. A heterozygous missense mutation c.749C>G (p.P250A) was identified in the FGFR3 gene through genetic testing. The patient exhibited typical features including coronal synostosis, bilateral hearing loss, right accessory auricle, and developmental delays and underwent surgery to repair the soft cleft palate. Cases of Muenke syndrome with cleft palate in the literature are relatively rare, and common associated symptoms include coronal suture craniosynostosis and hearing impairment. This article reports a differential diagnosis with other craniosynostosis syndromes and provides a reference for clinical diagnosis and treatment.
Humans ; Cleft Palate/surgery* ; Craniosynostoses/diagnosis* ; Mutation, Missense ; Palate, Soft/abnormalities* ; Receptor, Fibroblast Growth Factor, Type 3/genetics*

Here, we report a previously unrecognized syndromic neurodevelopmental disorder associated with biallelic loss-of-function variants in the RBM42 gene. The patient is a 2-year-old female with severe central nervous system (CNS) abnormalities, hypotonia, hearing loss, congenital heart defects, and dysmorphic facial features. Familial whole-exome sequencing (WES) reveals that the patient has two compound heterozygous variants, c.304C>T (p.R102*) and c.1312G>A (p.A438T), in the RBM42 gene which encodes an integral component of splicing complex in the RNA-binding motif protein family. The p.A438T variant is in the RRM domain which impairs RBM42 protein stability in vivo. Additionally, p.A438T disrupts the interaction of RBM42 with hnRNP K, which is the causative gene for Au-Kline syndrome with overlapping disease characteristics seen in the index patient. The human R102* or A438T mutant protein failed to fully rescue the growth defects of RBM42 ortholog knockout ΔFgRbp1 in Fusarium while it was rescued by the wild-type (WT) human RBM42. A mouse model carrying Rbm42 compound heterozygous variants, c.280C>T (p.Q94*) and c.1306_1308delinsACA (p.A436T), demonstrated gross fetal developmental defects and most of the double mutant animals died by E13.5. RNA-seq data confirmed that Rbm42 was involved in neurological and myocardial functions with an essential role in alternative splicing (AS). Overall, we present clinical, genetic, and functional data to demonstrate that defects in RBM42 constitute the underlying etiology of a new neurodevelopmental disease which links the dysregulation of global AS to abnormal embryonic development.
Female ; Animals ; Mice ; Humans ; Child, Preschool ; Intellectual Disability/genetics* ; Heart Defects, Congenital/genetics* ; Facies ; Cleft Palate ; Muscle Hypotonia

OBJECTIVE: To explore the genetic characteristics of a Chinese pedigree affected with van der Woude syndrome (VWS). METHODS: A proband who had visited the Drum Tower Hospital Affiliated to Nanjing University Medical School in May 2020 for "two previous pregnancies with cleft lip and palate" was selected as the study subject. Trio-whole exome sequencing (trio-WES) was carried out for the patient. Candidate variants were verified by Sanger sequencing of her pedigree members (8 individuals from four generations) and bioinformatic analysis. Chromosomal microarray analysis (CMA) was used to rule out copy number variations in the fetuses. RESULTS: Trio-WES revealed that the proband and her father had both harbored a heterozygous c.742G>T (p.G248C) missense variant of the IRF6 gene, for which her mother was of the wild type. The variant was located in a region with important functions and has not been reported previously. Prediction with several software suggested that it is likely to have a significant impact on the protein structure/function and is highly correlated with the specific phenotypes in this pedigree. Sanger sequencing confirmed co-segregation of the genotypes and phenotypes in the pedigree. Based on the guidelines from the American College of Medical Genetics and Genomics (ACMG), this variant was rated as likely pathogenic (PM1+PM2_Supporting+PP1+PP3+PP4). Based on the above results, pre-implantation genetic diagnosis was carried out for the proband, which has led to birth of a healthy offspring with normal results for both site testing and CMA. CONCLUSION The IRF6: c.742G>T (p.G248C) heterozygous variant probably underlay the VWS in this pedigree. Above finding has also enabled reproductive guidance for the proband.
Humans ; Female ; Cleft Lip/genetics* ; Cleft Palate/genetics* ; Pedigree ; DNA Copy Number Variations ; East Asian People ; Interferon Regulatory Factors/genetics* ; Mutation

Non-syndromic oral cleft (NSOC), a common birth defect, remains to be a critical public health problem in China. In the context of adjustment of childbearing policy for two times in China and the increase of pregnancy at older childbearing age, NSOC risk prediction will provide evidence for high-risk population identification and prenatal counseling. Genome-wide association study and second generation sequencing have identified multiple loci associated with NSOC, facilitating the development of genetic risk prediction of NSOC. Despite the marked progress, risk prediction models of NSOC still faces multiple challenges. This paper summarizes the recent progress in research of NSOC risk prediction models based on the results of extensive literature retrieval to provide some insights for the model development regarding research design, variable selection, model-build strategy and evaluation methods.
Humans ; Cleft Palate/genetics* ; Cleft Lip/genetics* ; Genome-Wide Association Study ; Genetic Predisposition to Disease ; Risk Factors ; Polymorphism, Single Nucleotide

The cranial neural crest plays a fundamental role in orofacial development and morphogenesis. Accordingly, mutations with impact on the cranial neural crest and its development lead to orofacial malformations such as cleft lip and palate. As a pluripotent and dynamic cell population, the cranial neural crest undergoes vast transcriptional and epigenomic alterations throughout the formation of facial structures pointing to an essential role of factors regulating chromatin state or transcription levels. Using CRISPR/Cas9-guided genome editing and conditional mutagenesis in the mouse, we here show that inactivation of Kat5 or Ep400 as the two essential enzymatic subunits of the Tip60/Ep400 chromatin remodeling complex severely affects carbohydrate and amino acid metabolism in cranial neural crest cells. The resulting decrease in protein synthesis, proliferation and survival leads to a drastic reduction of cranial neural crest cells early in fetal development and a loss of most facial structures in the absence of either protein. Following heterozygous loss of Kat5 in neural crest cells palatogenesis was impaired. These findings point to a decisive role of the Tip60/Ep400 chromatin remodeling complex in facial morphogenesis and lead us to conclude that the orofacial clefting observed in patients with heterozygous KAT5 missense mutations is at least in part due to disturbances in the cranial neural crest.
Animals ; Mice ; Chromatin Assembly and Disassembly ; Cleft Lip/genetics* ; Cleft Palate/genetics* ; DNA Helicases/metabolism* ; DNA-Binding Proteins ; Neural Crest/metabolism* ; Skull ; Transcription Factors/metabolism*

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.
Asians ; Case-Control Studies ; Cleft Lip/genetics* ; Cleft Palate/genetics* ; Genetic Predisposition to Disease ; Genome-Wide Association Study ; Genotype ; Humans ; Mutation ; Parents ; Polymorphism, Single Nucleotide ; Whole Exome Sequencing

OBJECTIVE: To explore whether WNT signaling pathway genes were associated with non-syndromic oral clefts (NSOC) based on haplotypes analyses among 1 008 Chinese NSOC case-parent trios. METHODS: The genome-wide association study (GWAS) data of 806 Chinese non-syndromic cleft lip with or without cleft palate (NSCL/P) trios and 202 Chinese non-syndromic cleft palate (NSCP) case-parent trios were drawn from the International Consortium to Identify Genes and Interactions Controlling Oral Clefts (ICOCs) study GWAS data set, whose Chinese study population were recruited from four provinces in China, namely Taiwan, Shandong, Hubei, and Sichuan provinces. The process of DNA genotyping was conducted by the Center for Inherited Disease Research in the Johns Hopkins University, using Illumina Human610-Quad v.1_B Bead Chip. The method of sliding windows was used to determine the haplotypes for analyses, including 2 SNPs haplotypes and 3 SNPs haplotypes. Haplotypes with a frequency lower than 1% were excluded for further analyses. To further assess the association between haplotypes and NSOC risks, and the transmission disequilibrium test (TDT) was performed. The Bonferroni method was adopted to correct multiple tests in the study, with which the threshold of statistical significance level was set as P < 0.05 divided by the number of tests, e.g P < 3.47×10^-4 in the current stu-dy. All the statistical analyses were performed by using plink (v1.07). RESULTS: After quality control, a total of 144 single nucleotide polymorphisms (SNPs) mapped in seven genes in WNT signaling pathway were included for the analyses among the 806 Chinese NSCL/P trios and 202 Chinese NSCP trios. A total of 1 042 haplotypes with frequency higher than 1% were included for NSCL/P analyses and another 1 057 haplotypes with frequency higher than 1% were included for NSCP analyses. Results from the TDT analyses showed that a total of 69 haplotypes were nominally associated with the NSCL/P risk among Chinese (P < 0.05). Another 34 haplotypes showed nominal significant association with the NSCP risk among Chinese (P < 0.05). However, none of these haplotypes reached pre-defined statistical significance level after Bonferroni correction (P>3.47×10^-4). CONCLUSION This study failed to observe any statistically significant associations between haplotypes of seven WNT signaling pathway genes and the risk of NSOC among Chinese. Further studies are warranted to replicate the findings here.
Cleft Lip/genetics* ; Cleft Palate/genetics* ; Genetic Predisposition to Disease ; Genome-Wide Association Study ; Genotype ; Haplotypes ; Humans ; Polymorphism, Single Nucleotide ; Wnt Signaling Pathway/genetics*

Genome-wide association studies (GWASs) are the most widely used method to identify genetic risk loci associated with orofacial clefts (OFC). However, despite the increasing size of cohort, GWASs are still insufficient to detect all the heritability, suggesting there are more associations under the current stringent statistical threshold. In this study, we obtained an integrated epigenomic dataset based on the chromatin conformation of a human oral epithelial cell line (HIOEC) using RNA-seq, ATAC-seq, H3K27ac ChIP-seq, and DLO Hi-C. Presumably, this epigenomic dataset could reveal the missing functional variants located in the oral epithelial cell active enhancers/promoters along with their risk target genes, despite relatively less-stringent statistical association with OFC. Taken a non-syndromic cleft palate only (NSCPO) GWAS data of the Chinese Han population as an example, 3664 SNPs that cannot reach the strict significance threshold were subjected to this functional identification pipeline. In total, 254 potential risk SNPs residing in active cis-regulatory elements interacting with 1 718 promoters of oral epithelium-expressed genes were screened. Gapped k-mer machine learning based on enhancers interacting with epithelium-expressed genes along with in vivo and in vitro reporter assays were employed as functional validation. Among all the potential SNPs, we chose and confirmed that the risk alleles of rs560789 and rs174570 reduced the epithelial-specific enhancer activity by preventing the binding of transcription factors related to epithelial development. In summary, we established chromatin conformation datasets of human oral epithelial cells and provided a framework for testing and understanding how regulatory variants impart risk for clefts.
Chromatin ; Cleft Lip/genetics* ; Cleft Palate/genetics* ; Epithelium ; Genome-Wide Association Study ; Humans

Objective: To explore the molecular mechanism of cleft palate in mice induced by 2, 3, 7, 8-tetrachlorodibenzo-p-dioxin (TCDD). Methods: The pregnant mice were randomly divided into TCDD-treated group (n=42) and control group (n=42). TCDD-treated group was given by gavage a single dose of TCDD (64 μg/kg) at 8: 00 AM on gestation day 10 (GD10) and the control group was given by gavage the isopyknic corn oil. At GD13-GD15, the fetal mice palate development was observed by HE staining. The mouse embryonic palatal mesenchymal cell proliferation was detected by 5-bromo-2-deoxyuridine (BrdU) immunofluorescence. The localization and expression of maternally expressed gene3 (MEG3) in mouse embryonic palatal mesenchymal cells was detected by situ hybridization and real-time PCR (RT-PCR). The key protein expressions of transforming growth factor-β (TGF-β)/Smad signaling pathway in mouse embryonic palatal mesenchyme were analyzed by Western blotting. The interaction of MEG3 and TGF-β receptor Ⅰ (TGF-βRⅠ) was examined by RNA binding protein immunoprecipitation (RIP). Results: At GD13 and GD14, compared with the control group, the ratio of BrdU-positive cells in the palatal mesenchyme of TCDD-treated fetuses decreased significantly (GD13, t=6.66, P=0.003; GD14, t=6.56, P=0.003). However, at GD15, the ratio of BrdU-positive cells was significantly increased (t=-5.98, P=0.004). MEG3 was mainly expressed in the nuclei of fetal mouse palatal mesenchymal cells, and the expression of MEG3 in TCDD group was significantly increased at GD13, GD14 and GD15(GD13, t=39.28, P=0.012; GD14, t=18.75, P=0.042; GD15, t=28.36, P=0.045). At GD14, TCDD decreased the levels of p-Smad2 and Smad4 in embryonic palate mesenchymal cells (p-Smad2, t=9.48, P=0.001;Smad4, t=63.10, P=0.001), whereas the expression of Smad7 was significantly increased at GD14 (t=30.77, P<0.001). The results of the RIP experiment showed that the amount of TGF-βRⅠ-bound MEG3 in mouse embryonic palatal mesenchymal cells in the TCDD group (23.940±1.301) was higher than that in the control group (8.537±1.523)(t=24.55, P<0.001). Conclusions: MEG3 is involved in the suppression of mouse embryonic palatal mesenchymal cell proliferation, functioning at least in part via interacting with the TGF-βRⅠ protein and thereby suppressing Smad signaling in the context of TCDD induced cleft palate.
Animals ; Bromodeoxyuridine ; Cleft Palate/genetics* ; Female ; Mice ; Mice, Inbred C57BL ; Palate/metabolism* ; Polychlorinated Dibenzodioxins/toxicity* ; Pregnancy