SNPs and interaction analyses of IRF6, MSX1 and PAX9 genes in patients with non‑syndromic cleft lip with or without palate

A comprehensive consolidation of data on the relationship between IRF6 polymorphisms and non-syndromic cleft lip/palate susceptibility: From 79 case-control studies

BACKGROUND

Non-syndromic cleft lip with or without cleft palate (NSCL/P) is a prevalent craniofacial birth defect on a global scale. A number of candidate genes have been identified as having an impact on NSCL/P. However, the association between interferon regulatory factor 6 (IRF6) polymorphisms and NSCL/P has yielded inconsistent results, prompting the need for a meta-analysis to obtain more accurate estimates.

METHODS

We conducted a thorough screening of all relevant articles published up until November 15, 2023, in online bibliographic databases. The statistical analysis of the collected data was performed using the Comprehensive Meta-Analysis (Version 4.0) software.

RESULTS

A total of 79 case-control studies, comprising 14,003 cases and 19,905 controls, were included in our analysis. The combined data indicated that the IRF6 rs642961 and rs2235371 polymorphisms were associated with an increased risk of NSCL/P in the overall population. However, no significant association was found between the rs2013162 and rs2235375 polymorphisms and the risk of NSCL/P in the overall population. Furthermore, subgroup analyses revealed significant correlations between the IRF6 rs642961, rs2235371, and rs2235375 polymorphisms and the risk of NSCL/P based on ethnic background and country of origin. Nevertheless, the rs2013162 polymorphism plays a protective role in Caucasians and mixed populations.

CONCLUSIONS

Our collective data indicates a significant association between the rs642961 and rs2235371 polymorphisms and the risk of NSCL/P in the overall population. The rs2235375 polymorphism could influence the susceptibility to NSCL/P based on ethnic background. Meanwhile, the rs2013162 polymorphism provides protective effects in Caucasian, mixed populations, and the Brazilian population.

Association of MSX1 Gene Variants with Nonsyndromic Cleft Lip and/or Palate in the Pakistani Population

OBJECTIVES

This study investigated the association of MSX1 gene variants rs3821949 and rs12532 with nonsyndromic cleft lip and/or palate (NSCL/P) in the Pakistani population.

DESIGN

Comparative cross-sectional study.Setting: Multicenter of CL/P malformation.Patients/Participants: Unrelated Non-Syndromic cleft Lip/Palate patients and healthy controls were enrolled.

METHODS

One hundred (n = 100) subjects with NSCL/P and n = 50 unrelated healthy controls were enrolled in a multicenter comparative cross-sectional study. A tetra amplification refractory mutation system (ARMS) polymerase chain reaction (PCR) was performed to analyze MSXI gene single nucleotide variants (SNVs).

RESULTS

Among 100 NSCL/P subjects, the majority were males (56%; male: female = 1.27: 1). Most of the cases (74%) had cleft lip and palate (CLP) compared to isolated clefts. Genotyping of MSX1 gene variant rs3821949 showed an increased risk for NSCL/P in various genetic models (P < 0.0001), and the A allele exhibited a more than 4-fold increased risk among cases (OR = 4.22: 95% CI = 2.16-8.22; P < 0.0001). Our investigation found no significant difference between the rs12532 variation and NSCL/P.

CONCLUSION

Our study findings suggest that MSX1 gene variants may increase predisposition to NSCL/P in the Pakistani population. Further studies comprising large samples are required to identify the genetic aetiology of NSCL/P among our people.

Genetic Risk Assessment of Nonsyndromic Cleft Lip with or without Cleft Palate by Linking Genetic Networks and Deep Learning Models

Recent deep learning algorithms have further improved risk classification capabilities. However, an appropriate feature selection method is required to overcome dimensionality issues in population-based genetic studies. In this Korean case-control study of nonsyndromic cleft lip with or without cleft palate (NSCL/P), we compared the predictive performance of models that were developed by using the genetic-algorithm-optimized neural networks ensemble (GANNE) technique with those models that were generated by eight conventional risk classification methods, including polygenic risk score (PRS), random forest (RF), support vector machine (SVM), extreme gradient boosting (XGBoost), and deep-learning-based artificial neural network (ANN). GANNE, which is capable of automatic input SNP selection, exhibited the highest predictive power, especially in the 10-SNP model (AUC of 88.2%), thus improving the AUC by 23% and 17% compared to PRS and ANN, respectively. Genes mapped with input SNPs that were selected by using a genetic algorithm (GA) were functionally validated for risks of developing NSCL/P in gene ontology and protein-protein interaction (PPI) network analyses. The IRF6 gene, which is most frequently selected via GA, was also a major hub gene in the PPI network. Genes such as RUNX2, MTHFR, PVRL1, TGFB3, and TBX22 significantly contributed to predicting NSCL/P risk. GANNE is an efficient disease risk classification method using a minimum optimal set of SNPs; however, further validation studies are needed to ensure the clinical utility of the model for predicting NSCL/P risk.

Association of Nucleotide Variants of GRHL3, IRF6, NAT2, SDC2, BCL3, and PVRL1 Genes with Nonsyndromic Cleft Lip With/Without Cleft Palate in Multigenerational Families: A Retrospective Study

Background:

Several genes are associated with the etiology of cleft lip and palate (CLP) in different populations. Many nucleotide variants on genes such as GRHL3, IRF6, NAT2, SDC2, BCL3, and PVRL1 were reported in different populations, but not studied in multigenerational cases in the Indian population.

Aim and Objective:

The aim of this study is to evaluate whether nucleotide variants rs41268753, rs861020, rs1041983, rs1042381, rs2965169, and rs10790332 are involved in the etiology of nonsyndromic CLP (NSCLP) in multigenerational Indian families.

Study Design:

Retrospective genetic study.

Materials and Methods:

Based on inclusion and exclusion criteria, 20 multigenerational families with nonsyndromic cleft lip with or without cleft palate (NSCL/P) were selected. Blood samples from both affected and unaffected participants were collected as a source of genomic DNA. Six nucleotide variants on these genes were genotyped to test for the association with NSCL/P. Genotyping was performed with the MassArray method. Genotype distribution was used to calculate the Hardy–Weinberg equilibrium using PLINK, a whole-genome association analysis toolset. The allelic association was compared among cases and controls using Chi-square test as implemented in PLINK. P ≤ 0.05 indicates statistical differences between groups.

Results:

No significant associations were found between individual single-nucleotide polymorphisms and NSCL/P. The odds ratio was 1.531, 1.198, 0.8082, 1.418, 1, and 0.5929 for polymorphisms rs41268753, rs861020, rs1041983, rs1042381, rs2965169, and rs10790332, respectively.

Conclusion:

Our findings suggest that among the multigenerational families in our population, the high-risk nucleotide variants GRHL3 rs41268753, IRF6 rs861020, NAT2 rs1041983, SDC2 rs1042381, BCL3 rs2965169, and PVRL1 rs10790332 are not associated with increased risk of NSCL/P.

Genetic Factors in Nonsyndromic Orofacial Clefts

Orofacial clefts (OFCs) are the most common congenital birth defects in humans and immediately recognized at birth. The etiology remains complex and poorly understood and seems to result from multiple genetic and environmental factors along with gene-environment interactions. It can be classified into syndromic (30%) and nonsyndromic (70%) clefts. Nonsyndromic OFCs include clefts without any additional physical or cognitive deficits. Recently, various genetic approaches, such as genome-wide association studies (GWAS), candidate gene association studies, and linkage analysis, have identified multiple genes involved in the etiology of OFCs. This article provides an insight into the multiple genes involved in the etiology of OFCs. Identification of specific genetic causes of clefts helps in a better understanding of the molecular pathogenesis of OFC. In the near future, it helps to provide a more accurate diagnosis, genetic counseling, personalized medicine for better clinical care, and prevention of OFCs.

Molecular Diagnostics and In Utero Therapeutics for Orofacial Clefts

Orofacial clefts and their management impose a substantial burden on patients, on their families, and on the health system. Under the current standard of care, affected patients are subjected to a lifelong journey of corrective surgeries and multidisciplinary management to replace bone and soft tissues, as well as restore esthetics and physiologic functions while restoring self-esteem and psychological health. Hence, a better understanding of the dynamic interplay of molecular signaling pathways at critical phases of palate development is necessary to pioneer novel prenatal interventions. Such pathways include transforming growth factor-β (Tgfβ), sonic hedgehog (Shh), wingless-integrated site (Wnt)/β-catenin, bone morphogenetic protein (Bmp), and fibroblast growth factor (Fgf) and its associated receptors, among others. Here, we summarize commonly used surgical methods used to correct cleft defects postnatally. We also review the advances made in prenatal diagnostics of clefts through imaging and genomics and the various in utero surgical corrections that have been attempted thus far. An overview of how key mediators of signaling that drive palatogenesis are emphasized in the context of the framework and rationale for the development and testing of therapeutics in animal model systems and in humans is provided. The pros and cons of in utero therapies that can potentially restore molecular homeostasis needed for the proper growth and fusion of palatal shelves are presented. The theme advanced throughout this review is the need to develop preclinical molecular therapies that could ultimately be translated into human trials that can correct orofacial clefts at earlier stages of development.

DOT: Gene-set analysis by combining decorrelated association statistics

Historically, the majority of statistical association methods have been designed assuming availability of SNP-level information. However, modern genetic and sequencing data present new challenges to access and sharing of genotype-phenotype datasets, including cost of management, difficulties in consolidation of records across research groups, etc. These issues make methods based on SNP-level summary statistics particularly appealing. The most common form of combining statistics is a sum of SNP-level squared scores, possibly weighted, as in burden tests for rare variants. The overall significance of the resulting statistic is evaluated using its distribution under the null hypothesis. Here, we demonstrate that this basic approach can be substantially improved by decorrelating scores prior to their addition, resulting in remarkable power gains in situations that are most commonly encountered in practice; namely, under heterogeneity of effect sizes and diversity between pairwise LD. In these situations, the power of the traditional test, based on the added squared scores, quickly reaches a ceiling, as the number of variants increases. Thus, the traditional approach does not benefit from information potentially contained in any additional SNPs, while our decorrelation by orthogonal transformation (DOT) method yields steady gain in power. We present theoretical and computational analyses of both approaches, and reveal causes behind sometimes dramatic difference in their respective powers. We showcase DOT by analyzing breast cancer and cleft lip data, in which our method strengthened levels of previously reported associations and implied the possibility of multiple new alleles that jointly confer disease risk.

Association between MSX1 rs12532 polymorphism with nonsyndromic unilateral complete cleft lip and palate and tooth agenesis

OBJECTIVES

To investigate the association of MSX1 rs12532 polymorphism with the risk of nonsyndromic unilateral complete cleft lip and palate (NSCLP) and tooth agenesis.

MATERIALS AND METHODS

The study is comprised of 384 individuals divided into 4 groups: group 1, patients with unilateral complete NSCLP and premolar agenesis (n = 57); group 2, patients with unilateral NSCLP without tooth agenesis (n = 117); group 3, patients with premolar agenesis without oral cleft (n = 53) and group 4 (n = 157), a control group with individuals without tooth agenesis and oral cleft. Genotyping of rs12532 was carried out with Taqman chemistry, and associations were investigated using logistic regression analyses.

RESULTS

Overall rs12532 allele and genotype distributions revealed no significant differences between the groups of NSCLP or tooth agenesis.

CONCLUSION

Although our results are consistent with a lack of association of MSX1 rs12532 and the risk of unilateral NSCLP and tooth agenesis, further studies with additional SNPs and a more diverse ethnic cohort are warranted.

Association Between an Interferon Regulatory Factor 6 Gene Polymorphism and Nonsyndromic Cleft Palate Risk

Background: Involvement of interferon regulatory factor 6 (IRF6) gene polymorphisms in nonsyndromic cleft palate (NSCP) risk remains controversial. This investigation was performed to evaluate the relationship between IRF6 gene polymorphisms and NSCP risk. Materials and Methods: Two hundred forty-one patients with NSCP (including 103 complete trio families) were recruited, and 242 unaffected individuals were included as controls. Polymorphisms for the IRF6 rs2235371, rs801619, rs642961, rs44844880, and rs8049367 loci were characterized in both groups. Furthermore, eligible studies were identified from the databases through June 1, 2017, and were included in a meta-analysis to enhance the robustness of our conclusions. Results: The IRF6 rs2235371 A allele and AA genotype in the case group were found at higher frequencies than in the control group (A allele: p < 0.0016; AA genotype: p < 0.0049). The IRF6 rs801619 AA genotype and G allele were associated with NSCP risk (G allele: p < 0.0061; AA genotype: p < 0.0195). At the IRF6 rs642961, rs44844880, and rs8049367 loci genotype and allele frequencies were not statistically different between the NSCP group and normal controls. In the meta-analysis, the IRF6 A/G gene polymorphism (rs2235371) and IRF6 A/G gene polymorphism (rs642961) were associated with NSCP risk in the general population, whereas the IRF6 A/C gene polymorphism (rs2013162) was not. Conclusion: The IRF6 A/G gene polymorphisms at rs2235371 and rs642961, but not the IRF6 A/C gene polymorphism rs2013162, were associated with NSCP risk.

MSX1 gene polymorphisms and non-syndromic cleft lip with or without palate (NSCL/P): A meta-analysis

OBJECTIVE

Non-syndromic cleft lip, with or without cleft palate (NSCL/P), is a common craniofacial birth defect, the risk of which is influenced from multiple genetic loci. Association study outcomes between single nucleotide polymorphisms (SNPs) near the muscle segment homeobox gene 1 (MSX1) and NSCL/P have been inconsistent. This compels a meta-analysis to obtain more precise estimates.

METHODS

From 15 publications, we examined 12 SNPs under six groups (SG), based on linkage disequilibrium. Pooled odds ratios and 95% confidence intervals were calculated under the standard genetic models. The pooled effects were subjected to subgroup, outlier, sensitivity, and funnel plot (publication bias) analyses.

RESULTS

Three of the six SGs showed significant associations. SG1 and SG4 effects indicated reduced risks. SG1 outcomes were attributed to outlier treatment, which the Asian outcomes validated. In contrast, increased risks were observed in SG3. All these significant outcomes were deemed robust by sensitivity analysis with no evidence of publication bias.

CONCLUSIONS

Our study shows eight MSX1 SNPs associated with risk of NSCL/P. SG1 and SG4 carriers are protected (up to 23%), but SG3 carriers are 1.3-fold susceptible. Outlier treatment unmasked the significant associations in SG1. Non-heterogeneity and robustness helped elevate the level of evidence in our significant findings.

A cleft lip and palate gene, Irf6, is involved in osteoblast differentiation of craniofacial bone

BACKGROUND

Interferon regulatory factor 6 (IRF6) plays a critical role in embryonic tissue development, including differentiation of epithelial cells. Besides orofacial clefting due to haploinsufficiency of IRF6, recent human genetic studies indicated that mutations in IRF6 are linked to small mandible and digit abnormalities. The function of IRF6 has been well studied in oral epithelium; however, its role in craniofacial skeletal formation remains unknown. In this study, we investigated the role of Irf6 in craniofacial bone development using comparative analyses between wild-type (WT) and Irf6-null littermate mice.

RESULTS

Immunostaining revealed the expression of IRF6 in hypertrophic chondrocytes, osteocytes, and bone matrix of craniofacial tissues. Histological analysis of Irf6-null mice showed a remarkable reduction in the number of lacunae, embedded osteocytes in matrices, and a reduction in mineralization during bone formation. These abnormalities may explain the decreased craniofacial bone density detected by micro-CT, loss of incisors, and mandibular bone abnormality of Irf6-null mice. To validate the autonomous role of IRF6 in bone, extracted primary osteoblasts from calvarial bone of WT and Irf6-null pups showed no effect on osteoblastic viability and proliferation. However, a reduction in mineralization was detected in Irf6-null cells.

CONCLUSIONS

Altogether, these findings suggest an autonomous role of Irf6 in regulating bone differentiation and mineralization. Developmental Dynamics 248:221-232, 2019. © 2019 Wiley Periodicals, Inc.

The Function and Regulatory Network of Pax9 Gene in Palate Development

Cleft palate, a common congenital deformity, can arise from disruptions in any stage of palatogenesis, including palatal shelf growth, elevation, adhesion, and fusion. Paired box gene 9 (Pax9) is recognized as a vital regulator of palatogenesis with great relevance to cleft palate in humans and mice. Pax9-deficient murine palatal shelves displayed deficient elongation, postponed elevation, failed contact, and fusion. Pax9 is expressed in epithelium and mesenchyme, exhibiting a dynamic expression pattern that changes according to the proceeding of palatogenesis. Recent studies highlighted the Pax9-related genetic interactions and their critical roles during palatogenesis. During palate growth, PAX9 interacts with numerous molecules and members of pathways (e.g., OSR2, FGF10, SHOS2, MSX1, BARX1, TGFβ3, LDB1, BMP, WNT β-catenin dependent, and EDA) in the mesenchyme and functions as a key mediator in epithelial-mesenchymal communications with FGF8, TBX1, and the SHH pathway. During palate elevation, PAX9 is hypothesized to mediate the time point of the elevation event in the anterior and posterior parts of the palatal shelves. The delayed elevation of Pax9 mutant palatal shelves probably results from abnormal expressions of a series of genes ( Osr2 and Bmpr1a) leading to deficient palate growth, abnormal tongue morphology, and altered hyaluronic acid distribution. The interactions between PAX9 and genes encoding the OSR2, TGFβ3, and WNT β-catenin-dependent pathways provide evidence that PAX9 might participate in the regulation of palate fusion. This review summarizes the current understanding of PAX9’s functions and emphasizes the interactions between PAX9 and vital genes during palatogenesis. We hope to provide some clues for further exploration of the function and mechanism of PAX9, especially during palate elevation and fusion events.

Gene-Gene Interactions among SPRYs for Nonsyndromic Cleft Lip/Palate

Nonsyndromic cleft lip with or without cleft palate (NSCL/P) is a common birth defect with a complex genetic architecture. Gene-gene interactions have been increasingly regarded as contributing to the etiology of NSCL/P. A recent genome-wide association study revealed that a novel single-nucleotide polymorphism at SPRY1 in 4q28.1 showed a significant association with NSCL/P. In the current study, we explored the role of 3 SPRY genes in the etiology of NSCL/P by detecting gene-gene interactions: SPRY1, SPRY2, and SPRY4-with SPRY3 excluded due to its special location on the X chromosome. We selected markers in 3 SPRY genes to test for gene-gene interactions using 1,908 case-parent trios recruited from an international consortium established for a genome-wide association study of nonsyndromic oral clefts. As the trios came from populations with different ancestries, subgroup analyses were conducted among Europeans and Asians. Cordell's method based on conditional logistic regression models was applied to test for potential gene-gene interactions via the statistical package TRIO in R software. Gene-gene interaction analyses yielded 10 pairs of SNPs in Europeans and 6 pairs in Asians that achieved significance after Bonferroni correction. The significant interactions were confirmed in the 10,000-permutation tests (empirical P = 0.003 for the most significant interaction). The study identified gene-gene interactions among SPRY genes among 1,908 NSCL/P trios, which revealed the importance of potential gene-gene interactions for understanding the genetic architecture of NSCL/P. The evidence of gene-gene interactions in this study also provided clues for future biological studies to further investigate the mechanism of how SPRY genes participate in the development of NSCL/P.

Association between IRF6 and 8q24 polymorphisms and nonsyndromic cleft lip with or without cleft palate: Systematic review and meta-analysis

BACKGROUND

We conducted a systematic review and meta-analysis of interferon regulatory factor 6 and 8q24 polymorphisms with nonsyndromic cleft lip with/without cleft palate (NSCL/P).

METHODS

Data extraction was independently performed by two reviewers. Genotypic effects of four polymorphisms from 31 studies were pooled separately by ethnicity using a mixed-effect logit model with accounting for heterogeneity.

RESULTS

For rs2235371, AA and GA carried, respectively, 51% (95% confidence interval [CI], 37%-61%) and 42% (95% CI, 32%-50%) lower risks of NSCL/P than GG genotypes in Asians, but these genotypes were not significant in Caucasians. For rs2013162, only AA was significant, that is, carried 0.65 (95% CI, 0.52-0.82) times lower odds than CC in Caucasians but not for Asians. For rs642961, AA and GA genotypes, respectively, carried 2.47 (95% CI, 1.41-4.35) and 1.40 (95% CI, 1.12-1.75) times higher odds in Asian, and 2.03 (95% CI, 1.52-2.71) and 1.58 (95% CI, 1.37-1.82) times higher odds in Caucasians compare with GG genotypes. For rs987525, AA and CA genotypes carried 2.27 (95% CI, 1.43-3.60) and 1.34 (95% CI, 1.02-1.77) times higher odds in Asian, and 5.25 (95% CI, 3.98-6.91) and 2.13 (95% CI-1.82, 2.49) times higher odds in Caucasians, and 1.42 (95% CI, 1.10-1.82) and 1.28 (95% CI, 1.09-1.50) times higher odds in mixed ethnicities compared with CC genotypes. These variant effects remained significant based on applying Bonferroni corrected-thresholds, except in the mixed ethnicity.

CONCLUSION

We show robust variant effects in NSCL/P. Considering them with other genes and risk factors might be useful to improve prediction of NSCL/P occurrence. Birth Defects Research (Part A) 106:773-788, 2016. © 2016 The Authors Birth Defects Research Part A: Clinical and Molecular Teratology Published by Wiley Periodicals, Inc.

WITHDRAWN: The MSX1 gene polymorphisms rs12532 and rs1001179 confer an increased risk of non syndromic orofacial clefts

This article has been withdrawn at the request of the author(s) and/or editor. The Publisher apologizes for any inconvenience this may cause. The full Elsevier Policy on Article Withdrawal can be found at http://www.elsevier.com/locate/withdrawalpolicy.

Association between BMP4 rs17563 polymorphism and NSCL/P risk: a meta-analysis

Objective. To investigate the association between bone morphogenetic protein 4 (BMP4) rs17563 polymorphism and nonsyndromic cleft lip with or without palate (NSCL/P) risk. Methods. Four online databases were researched and the related publications were collected. Odds ratio (OR) with 95% confidence interval (CI) was applied to assess the relationship; publication bias, metaregression, and sensitivity analysis were conducted to guarantee the strength of results. Results. Six published case-control studies were collected. Overall, no significant association between BMP4 rs17563 polymorphism and NSCL/P risk was found. It was notable that significant susceptibility on different ethnicity was observed in the stratified analysis. For Chinese population, the BMP4 rs17563 polymorphism was a significantly increased risk for NSCL/P (C versus T: OR = 1.52, 95% CI = 1.28–1.82, P < 0.01, I² = 0%; CC versus TT: OR = 2.58, 95% CI = 1.74–3.82, P < 0.01, I² = 0%; TC + CC versus TT: OR = 1.45, 95% CI = 1.14–1.84, P < 0.01, I² = 0%; CC versus TT + TC: OR=2.46, 95% CI = 1.46–4.14, P < 0.01, I² = 47.0%). On the contrary, significantly protective effects were found in Brazilian population (C versus T: OR = 0.69, 95% CI = 0.50–0.96, P = 0.03, I² = 68.5%; TC versus TT: OR = 0.52, 95% CI = 0.40–0.68, P < 0.01, I² = 0%; TC + CC versus TT: OR = 0.52, 95% CI = 0.35–0.78, P < 0.010, I² = 54.4%). Conclusion. This meta-analysis indicated that BMP4 rs17563 polymorphism could play a different role during the development of NSCL/P based on ethnicity diversity.