Correlation between FOXN3-SIN3A complex expression in peripheral blood and non-syndromic cleft lip and palate in Xinjiang.

OBJECTIVES: This work aimed to study the correlation between FOXN3-SIN3A complex expression and non-syndromic oral clefts (NSOC) in Xinjiang. METHODS: In this study, 60 patients with NSOC attending the People's Hospital of Xinjiang Uygur Autonomous Region were recruited into the case group, including 30 cleft lip with or without cleft palate (NSCL/P), 30 cleft palate only (CPO), and 30 healthy children in the control group. The expression levels of FOXN3, SIN3A, and NEAT1 in peripheral blood of each group were detected by high-throughput second-generation sequencing technology and quantitative reverse transcription polymerase chain reaction (RT-qPCR). Receiver operating characteristic (ROC) curve and area under the curve (AUC) were used to analyze the diagnostic efficiency of NSOC. RESULTS: The comparison of the NSOC and control groups showed that FOXN3, SIN3A, and NEAT1 genes increased compared with the control group. The differences were all statistically significant (P<0.05). The AUCs of FOXN3, SIN3A, and NEAT1 in the NSCL/P group were 0.933 [95%CI=(0.864, 1.000)], 0.822 [(95%CI=(0.713, 0.932)], and 1.000[95%CI= (1.000, 1.000)], respectively. The AUCs of FOX-N3, SIN3A, and NEAT1 in the CPO group were 0.891 [95%CI=(0.806, 0.976)], 0.688 [95%CI=(0.552, 0.824)], and 1.000 [95%CI=(1.000, 1.000)], respectively. CONCLUSIONS: The results showed a correlation between the rising gene expression of FOXN3, SIN3A, and NEAT1 in peripheral blood and the occurrence of NSOC in Xinjiang. This work provides a theoretical basis for further study of the FOXN3-SIN3A complex as biomarkers to facilitate the early screening, disease prediction, and early prevention of NSOC.

[Gene-gene/gene-environment interaction of transforming growth factor-ß signaling pathway and the risk of non-syndromic oral clefts].

OBJECTIVE: To explore the association between polymorphisms of transforming growth factor-ß (TGF-ß) signaling pathway and non-syndromic cleft lip with or without cleft palate (NSCL/P) among Asian populations, while considering gene-gene interaction and gene-environment interaction. METHODS: A total of 1 038 Asian NSCL/P case-parent trios were ascertained from an international consortium, which conducted a genome-wide association study using a case-parent trio design to investigate the genes affec-ting risk to NSCL/P. After stringent quality control measures, 343 single nucleotide polymorphism (SNP) spanning across 10 pivotal genes in the TGF-ß signaling pathway were selected from the original genome-wide association study(GWAS) dataset for further analysis. The transmission disequilibrium test (TDT) was used to test for SNP effects. The conditional Logistic regression models were used to test for gene-gene interaction and gene-environment interaction. Environmental factors collected for the study included smoking during pregnancy, passive smoking during pregnancy, alcohol intake during pregnancy, and vitamin use during pregnancy. Due to the low rates of exposure to smoking during pregnancy and alcohol consumption during pregnancy (<3%), only the interaction between maternal smoking during pregnancy and multivitamin supplementation during pregnancy was analyzed. The threshold for statistical significance was rigorously set at P =1.46×10-4, applying Bonferroni correction to account for multiple testing. RESULTS: A total of 23 SNPs in 4 genes yielded nominal association with NSCL/P (P<0.05), but none of these associations was statistically significant after Bonferroni' s multiple test correction. However, there were 6 pairs of SNPs rs4939874 (SMAD2) and rs1864615 (TGFBR2), rs2796813 (TGFB2) and rs2132298 (TGFBR2), rs4147358 (SMAD3) and rs1346907 (TGFBR2), rs4939874 (SMAD2) and rs1019855 (TGFBR2), rs4939874 (SMAD2) and rs12490466 (TGFBR2), rs2009112 (TGFB2) and rs4075748 (TGFBR2) showed statistically significant SNP-SNP interaction (P<1.46×10-4). In contrast, the analysis of gene-environment interactions did not yield any significant results after being corrected by multiple testing. CONCLUSION: The comprehensive evaluation of SNP associations and interactions within the TGF-ß signaling pathway did not yield any direct associations with NSCL/P risk in Asian populations. However, the significant gene-gene interactions identified suggest that the genetic architecture influencing NSCL/P risk may involve interactions between genes within the TGF-ß signaling pathway. These findings underscore the necessity for further investigations to unravel these results and further explore the underlying biological mechanisms.

The potential up-regulation risk of 3' UTR SNP (rs10787760 G > A) for the VAX1 gene is associated with NSCLP in the northwest Chinese population.

AIMS: To investigate the association between single nucleotide polymorphisms (SNPs) in 3'UTR region of VAX1, SYT14 and PAX7 genes and the risk of non-syndromic cleft palate (NSCLP) in a northwest Chinese population. MAIN METHODS: A case-control study was conducted in 406 normal controls and 399 NSCLP patients. Using iMLDRTM genotyping technology, eight SNPs of three genes ((rs10787760, rs7086344 at VAX1), (rs1010113, rs851114, and rs485874 at PAX7), and (rs61820397, rs4609425, rs12133399 at SYT14)) were genotyped to investigate the differences in alleles and genotype distribution frequencies between NSCLP patients and healthy controls. RNA Folding Form software was used to predict RNA secondary structure and expression vectors were constructed to explore the function of the relevant SNP. The effect of SNP polymorphism of gene transcription and translation was assessed using qPCR and Western blot analysis. KEY FINDINGS: Among the eight SNPs of three genes, rs10787760 of VAX1 gene was found to be associated with an increased risk of NSCLP (OR = 1.341, CI = 1.004-1.790) and the GA genotype of rs10787760 increased the risk of cleft lip and/or palate (CL/P) about 1.42 times (p < 0.05), and carrying the A allele might increase the risk of NSCL/P in male (OR = 1.356, 95 % CI = 1.010-1.823). But there was no association observed with cleft palate only (CPO). Cell function experiments revealed that the G to A mutation in rs10787760 up-regulated GFP-VAX1 transcriptional level by 2.39 and 3.13 times in two cell lines respectively, and enhance the protein expression of the VAX1 gene further. RNA secondary structure study showed that the rs10787760 (G > A) had two different secondary structures in 3'UTR region. SIGNIFICANCE: The rs10787760 variant in the 3'UTR region of VAX1 gene is associated with CL/P in northwest Chinese population. We hypothesize that the machanism of it might be caused by the RNA differenct fold in the 3'UTR region caused by the polymorphism of the gene. LEVEL OF EVIDENCE: Original Reports.

Integrative Multi-omics Analysis Identifies Genetic Variants Contributing to Non-syndromic Cleft Lip with or without Cleft Palate.

OBJECTIVE: To provide novel insights into the aetiology of non-syndromic cleft lip with or without cleft palate (NSCL/P) by integrating multi-omics data and exploring susceptibility genes associated with NSCL/P. METHODS: A two-stage genome-wide association study (GWAS) of NSCL/P was performed, involving a total of 1,069 cases and 1,724 controls. Using promoter capture Hi-C (pCHi-C) datasets in human embryonic stem cells (hESC) and chromatin immunoprecipitation sequencing (ChIP-seq) in craniofacial tissues, we filtered out single nucleotide polymorphisms (SNPs) with active cis-regulation and their target genes. Additionally, we employed expression quantitative trait loci (eQTL) analysis to identify candidate genes. RESULTS: Thirteen SNPs were identified as cis-regulation units associated with the risk of NSCL/P. Five of these were proven to be active in chromatin states in early human craniofacial development (rs7218002: odds ratio [OR] 1.50, P = 8.14E-08; rs835367: OR 0.78, P = 3.48E- 05; rs77022994: OR 0.55, P = 1.05E-04; rs961470: OR 0.73, P = 1.38E-04; rs17314727: OR 0.73, P = 1.85E-04). Additionally, pCHi-C and eQTL analysis prioritised three candidate genes (rs7218002: NTN1, rs835367: FGGY, LINC01135). NTN1 and FGGY were expressed in mouse orofacial development. Deficiencies in NTN1, FGGY and LINC01135 were associated with cleft palate and cleft lip, abnormal facial shape and bifid uvula, and abnormality of the face, respectively. CONCLUSION: Our study identified five SNPs (rs7218002, rs835367, rs77022994, rs961470 and rs17314727) and three susceptibility genes (NTN1, FGGY and LINC01135) associated with NSCL/P. These findings contribute to a better understanding of the genetic factors involved.

Progress of epigenetic modification of SATB2 gene in the pathogenesis of non-syndromic cleft lip and palate.

Non-syndromic Cleft Lip and Palate (NSCLP) is one of the most common congenital craniofacial malformations. However, there is no enough knowledge about its mechanism, even through many relevant studies verify that cleft lip and palate is caused by interactions between environmental and genetic factors. SATB2 gene is one of the most common candidate genes of NSCLP, and the development of epigenetics provides a new direction on pathogenesis of cleft lip and palate. This review summarizes SATB2 gene in the pathogenesis of non-syndromic cleft lip and palate, expecting to provide strategies to prevent and treat cleft and palate in the future.

Genetic variants in BCL-2 family genes influence the risk of non-syndromic cleft lip with or without cleft palate.

BACKGROUND: The BCL-2 family is crucial for cell death regulation and is involved in development, tissue homeostasis, and immunity. This study aimed to investigate the association between genetic variants in BCL-2 family genes and non-syndromic cleft lip with or without cleft palate (NSCL/P) risk. METHODS: A two-stage case-control study was conducted in this association study. Gene-based analysis using Multi-marker Analysis of GenoMic Annotation was performed in the first stage cohort, which included 565 cases and 1269 controls. A logistic regression model was employed to assess the effect of single nucleotide polymorphisms (SNPs) on susceptibility to NSCL/P. Candidate SNPs were replicated by extra dbGaP case-parent trios. Haploreg, RegulomeDB, and UCSC Genome Browser were used to identify enhancer effects of promising SNPs. Bulk RNA sequencing data obtained from the Gene Expression Omnibus was used to identify co-expressed genes. Single-cell RNA sequencing dataset was used to infer the cell population of the candidate gene. The "Monocle" package was used to analyze the pseudotime cell trajectories. RESULTS: Rs3943258 located in the enhancer region was associated with the risk of NSCL/P (Pmeta = 5.66 × 10-04 ) and exhibited an eQTL effect for BCL2 (P = 3.96 × 10-02 ). Co-expression and pathway enrichment analysis revealed that genes related to Bcl2 were significantly enriched in the PI3K-Akt signaling pathway, MAPK signaling pathway, and Wnt signaling pathway. Five cell clusters were identified in single-cell RNA sequencing, and Bcl2 was mainly located in the mesenchyme. CONCLUSION: The rs3943258 located within BCL2 was probably related to NSCL/P susceptibility.

Genetic Inheritance Models of Non-Syndromic Cleft Lip with or without Palate: From Monogenic to Polygenic.

Non-syndromic cleft lip with or without palate (NSCL/P) is a prevalent birth defect that affects 1/500-1/1400 live births globally. The genetic basis of NSCL/P is intricate and involves both genetic and environmental factors. In the past few years, various genetic inheritance models have been proposed to elucidate the underlying mechanisms of NSCL/P. These models range from simple monogenic inheritance to more complex polygenic inheritance. Here, we present a comprehensive overview of the genetic inheritance model of NSCL/P exemplified by representative genes and regions from both monogenic and polygenic perspectives. We also summarize existing association studies and corresponding loci of NSCL/P within the Chinese population and highlight the potential of utilizing polygenic risk scores for risk stratification of NSCL/P. The potential application of polygenic models offers promising avenues for improved risk assessment and personalized approaches in the prevention and management of NSCL/P individuals.

Development of Artificial Neural Network-Based Prediction Model for Evaluation of Maxillary Arch Growth in Children with Complete Unilateral Cleft Lip and Palate.

INTRODUCTION: Cleft lip and palate (CLP) are the most common congenital craniofacial deformities that can cause a variety of dental abnormalities in children. The purpose of this study was to predict the maxillary arch growth and to develop a neural network logistic regression model for both UCLP and non-UCLP individuals. METHODS: This study utilizes a novel method incorporating many approaches, such as the bootstrap method, a multi-layer feed-forward neural network, and ordinal logistic regression. A dataset was created based on the following factors: socio-demographic characteristics such as age and gender, as well as cleft type and category of malocclusion associated with the cleft. Training data were used to create a model, whereas testing data were used to validate it. The study is separated into two phases: phase one involves the use of a multilayer neural network and phase two involves the use of an ordinal logistic regression model to analyze the underlying association between cleft and the factors chosen. RESULTS: The findings of the hybrid technique using ordinal logistic regression are discussed, where category acts as both a dependent variable and as the study's output. The ordinal logistic regression was used to classify the dependent variables into three categories. The suggested technique performs exceptionally well, as evidenced by a Predicted Mean Square Error (PMSE) of 2.03%. CONCLUSION: The outcome of the study suggests that there is a strong association between gender, age, and cleft. The difference in width and length of the maxillary arch in UCLP is mainly related to the severity of the cleft and facial growth pattern.

Rs28446116 in PTCH1 is associated with non-syndromic cleft lip with or without palate in the Ningxia population, China.

OBJECTIVES: To investigate the association between PTCH1 single nucleotide polymorphism(SNP) and non-syndromic cleft lip with or without palate (NSCL/P) in the Ningxia Hui Autonomous region and predict the function of single nucleotide polymorphism through bioinformatics analysis. DESIGN: A case-control analysis of 31 single nucleotide polymorphism locus alleles on PTCH1 gene (504 cases and 455 controls) was performed to explore the association between PTCH1 gene polymorphisms and non-syndromic cleft lip with or without palate in Ningxia region. Transcription factors, 3D single nucleotide polymorphism and other related information of single nucleotide polymorphism loci with statistical significance were screened by the case-control experiments, and then analyzed the corresponding transcription factors through the NCBI database. RESULTS: The case-control study showed that 5 of the 31 single nucleotide polymorphism loci rs357564 (P = 0.0233), rs1805155 (P = 0.0371), rs28446116 (P = 0.0408), rs2282041 (P = 0.0439), rs56119276 (P = 0.0256) had statistically significant differences in allele frequencies between the case and control groups. Bioinformatics analysis revealed that EP300 and RUNX3, among the transcription factors associated with rs28446116, may be associated with the development of non-syndromic cleft lip with or without palate. CONCLUSION: PTCH1 gene may be associated with the occurrence of non-syndromic cleft lip with or without palate in the Ningxia region, which may be related to the role of EP300 and RUNX3 in the development of cleft lip and palate.

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

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

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

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

Mechanical dissection and culture of mouse cranial neural crest cells.

Owing to the contribution of cranial neural crest cells (CNCCs) to the majority of craniofacial structures, they have been studied extensively for the pathogenesis of craniofacial diseases. To investigate and summarize how to isolate and culture the CNCCs from wild-type mice, a literature search was performed in online databases (PubMed and Web of Science) using optimized keywords "mouse," "cranial neural crest cell" and "culture." The literature was checked by two investigators according to the screening and exclusion criteria. Initially, 197 studies were retrieved from PubMed and 169 from Web of Science, and after excluding replicate studies, 293 articles were considered. Finally, 17 studies met all the criteria and were included in this review. The results showed that obtaining purified stem cells and balancing the need to promote cell growth and prevent unwanted early cell differentiation were the two key points in the isolation and culture of CNCCs. However, no standard criteria are available for answering these questions. Thus, it is important to emphasize the necessity for standardization of CNCC isolation, culture, and identification in research on craniofacial diseases.

BMP4 Regulates EMT to be Involved in non-Syndromic Cleft lip With or Without Palate.

OBJECTIVE: In the previous study, we identified bone morphogenetic protein 4 (BMP4) responsible for non-syndromic cleft lip with or without cleft palate (NSCL/P). We aimed to elucidate the effects and mechanisms of BMP4 on epithelial-mesenchymal transition (EMT) through Smad1 signaling pathway to be involved in NSCL/P. METHODS: The human oral epidermoid carcinoma cells (KBs) were transfected with plasmids or small interfering RNA (siRNA) to build the models. The migration of the cells was evaluated by transwell assay. Western blotting and quantitative real-time reverse transcription-polymerase chain reaction (qRT-PCR) were used to detect the expressions of BMP4, E-cadherin, N-cadherin, EMT-related transcription factors snal1 and snal2, matrix metalloproteinase 2 (MMP2), MMP9, Smad1, and phosphorylated Smad1. RESULTS: In the overexpression group, the migration number of cells was increased significantly. The protein expression of E-cadherin was decreased significantly, while the protein expression level of the N-cadherin was increased significantly. The protein and mRNA expressions of MMP2, MMP9, snal1, and snal2 were significantly higher. The expression level of Smad1 was not significantly changed, while the phosphorylation of Smad1 was significantly increased. In the BMP4-siRNA group, the migrating number cells was significantly decreased. The protein expression of E-cadherin was increased significantly, while the expression of N-cadherin was significantly decreased. The protein and mRNA expressions of MMP2, MMP9, snal1, and snal2 were significantly lower than that of the control group. The expressions of Smad1 and phosphorylation of Smad1 were not significantly changed. CONCLUSIONS: BMP4 enhances cell migration and promotes cell EMT through Smad1 signaling pathway. Abnormal BMP4 mediates migration and EMT through other relevant signaling pathways resulting in NSCL/P. The study provides new insight into the mechanisms of NSCL/P associated with BMP4.n.

Two rare variants reveal the significance of Grainyhead-like 3 Arginine 391 underlying non-syndromic cleft palate only.

OBJECTIVES: Non-syndromic cleft palate only (NSCPO) is one of the most common craniofacial birth defects with largely undetermined genetic etiology. It has been established that Grainyhead-like 3 (GRHL3) plays an essential role in the pathogenesis of NSCPO. This study aimed to identify and verify the first-reported GRHL3 variant underlying NSCPO among the Chinese cohort. METHODS: We performed whole-exome sequencing (WES) on a Chinese NSCPO patient and identified a rare variant of GRHL3 (p.Arg391His). A validated deleterious variant p.Arg391Cys was introduced as a positive control. Zebrafish embryos injection, reporter assays, live-cell imaging, and RNA sequencing were conducted to test the pathogenicity of the variants. RESULTS: Zebrafish embryos microinjection demonstrated that overexpression of the variants could disrupt the normal development of zebrafish embryos. Reporter assays showed that Arg391His disturbed transcriptional activity of GRHL3 and exerted a dominant-negative effect. Interestingly, Arg391His and Arg391Cys displayed distinct nuclear localization patterns from that of wild-type GRHL3 in live-cell imaging. Bulk RNA sequencing suggested that the two variants changed the pattern of gene expression. CONCLUSIONS: In aggregate, this study identified and characterized a rare GRHL3 variant in NSCPO, revealing the critical role of Arginine 391 in GRHL3. Our findings will help facilitate understanding and genetic counseling of NSCPO.

Target sequencing reveals the association between variants in VAX1 and NSCL/P in Chinese population.

OBJECTIVE: A significant genetic association between rs7078160 in VAX1 and NSCL/P has been established through genome-wide association studies (GWAS), and we previously replicated the association in the Chinese population. The critical issue in the post-GWAS era is to identify functional variations that have a real impact on disease in the susceptible regions highlighted by GWAS. This study aimed to elucidate functional variants in VAX1 fully. MATERIALS AND METHODS: Firstly, target sequencing was performed on 159 NSCL/P patients, followed by association analysis to discover disease-associated single-nucleotide polymorphisms (SNPs); we then replicated the findings using a larger sample (1626 cases, 2255 controls) and investigated how candidate SNPs affect disease occurrence using extensive annotation databases. Additionally, we compared the genetic profiles of NSCL/P subtypes. RESULTS: In this study, 6 SNPs in VAX1 were identified to be associated with NSCL/P in the Western Han Chinese population. Five of them were predicted to influence transcriptional factor-biding ability and were expression quantitative trait loci (eQTLs) of nearby genes in multiple tissues. CONCLUSION: The previously reported association between rs7078160 and NSCL/P was successfully replicated. Moreover, our findings firstly revealed that 5 SNPs in VAX1 are associated with NSCL/P in the Western Han Chinese population. LEVEL OF EVIDENCE: Original Reports.

Reading-related Brain Function Restored to Normal After Articulation Training in Patients with Cleft Lip and Palate: An fMRI Study.

Cleft lip and/or palate (CLP) are the most common craniofacial malformations in humans. Speech problems often persist even after cleft repair, such that follow-up articulation training is usually required. However, the neural mechanism behind effective articulation training remains largely unknown. We used fMRI to investigate the differences in brain activation, functional connectivity, and effective connectivity across CLP patients with and without articulation training and matched normal participants. We found that training promoted task-related brain activation among the articulation-related brain networks, as well as the global attributes and nodal efficiency in the functional-connectivity-based graph of the network. Our results reveal the neural correlates of effective articulation training in CLP patients, and this could contribute to the future improvement of the post-repair articulation training program.

Rs9891446 in NTN1 is associated with right-side cleft lip in Han Chinese population.

OBJECTIVES: This study aimed to reveal the association between single-nucleotide polymorphism of NTN1 and subtypes of non-syndromic cleft lip with or without cleft palate (NSCL/P) in Han Chinese Population. DESIGN: Initially, we selected three single-nucleotide polymorphisms (SNP) (rs4791331, rs4791774 and rs9891446) in NTN1 from previous genetic studies. Then we recruited two Han Chinese cohorts (2004 cases and 1823 controls) and divided cases into subgroups: non-syndromic right-side cleft lip, non-syndromic left-side cleft lip, non-syndromic bilateral cleft lip and non-syndromic cleft lip with palate to further evaluate the associations between the subtypes of NSCL/P and SNPs in NTN1. PLINK and Haploview program were utilized to analyze the data. RESULTS: In the association analysis under additive model, we found that G allele at rs9891446 could specifically increase the risk of right-side cleft lip (P = 0.0073, OR = 1.44, 95%CI: 1.1-1.88), which was consistent with the results of association analysis under genotypic model. CONCLUSION: This study showed that rs9891446 of NTN1 was specifically associated with right-side cleft lip in Han Chinese, which indicates that different subtypes of non-syndromic cleft lip have distinct genetic background.

[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.

Association of rs2013162 and rs2235375 Polymorphisms in IRF6 Gene with Susceptibility to Non-Syndromic Cleft Lip and Palate.

Background: Non-syndromic cleft lip occurs by the interaction of environmental and genetic factors. The purpose of the current study was to analyze the association of Single Nucleotide Polymorphisms (SNPs) in IRF6 and NSCL/P in an Iranian population. Methods: A group of 105 children with NSCL/P and 185 normal controls were included in the current study. Genotyping of IRF6 rs2013162 and rs2235375 was performed by Polymerase Chain Reaction-Restriction Fragment Length Polymorphism (PCR-RFLP) method. Results: A substantial association of AA and CA genotypes in rs2013162 with the risk of NSCL/P (AA vs. CC; OR=2.36; 95%CI [1.05-5.29], p=0.004; and CA vs. CC; OR=0.47; 95%CI [0.28-0.79], p=0.018) was found. However, there were no important associations between A allele and risk of NSCL/P (p=0.980). According to logistic regression analysis results, subjects with GG genotype and G allele in rs2235375 polymorphism had increased risk of NSCL/P. Conclusion: The IRF6 polymorphisms are associated with the susceptibility to NSCL/P in Iranian population.

Targeted Re-Sequencing of the 2p21 Locus Identifies Non-Syndromic Cleft Lip Only Novel Susceptibility Gene ZFP36L2.

rs7590268 present on the 2p21 locus was identified to be associated with non-syndromic cleft lip with or without cleft palate (NSCL/P) in several populations, including the Chinese Han population, indicating that 2p21 was a susceptibility locus for NSCL/P. However, previous studies have only identified common single-nucleotide polymorphism (SNP) within the THADA gene, neglecting the rare variants and other genes in 2p21; thus, this study was designed to investigate additional variants and novel susceptibility genes in 2p21. A total of 159 NSCL/P patients and 542 controls were recruited in the discovery phase, whereas 1830 NSCL/P patients and 2,436 controls were recruited in the replication phase. After targeted region sequencing, we performed association and burden analyses for the common and rare variants, respectively. Furthermore, RNA-seq, proliferation assay and cell cycle analysis were performed to clarify the possible function of the candidate gene ZFP36L2. Association analysis showed that four SNPs were specifically associated with non-syndromic cleft lip only (NSCLO) and two SNPs were associated with both NSCLO and NSCL/P. Burden analysis indicated that ZFP36L2 was associated with NSCLO (p = .0489, OR = 2.41, 95% CI: 0.98-5.90). Moreover, SNPs in the ZFP36L2 targeted gene JUP were also associated with NSCLO. ZFP36L2 also inhibited cell proliferation and induced G2 phase arrest in the GMSM-K cell line. Therefore, we proposed that ZFP36L2 is a novel susceptibility gene of NSCLO in the 2p21 locus, which could lead to NSCLO by modulating cell proliferation and cycle.