Association between FOXE1 and non-syndromic orofacial clefts in a northeastern Chinese population

Association of selected gene variants with nonsyndromic orofacial clefts in Kuwait

INTRODUCTION AND OBJECTIVES

Non-syndromic orofacial clefts (NSOFCs) are complex congenital abnormalities involving both environmental and genetic factors involved in orofacial development. This study aimed to investigate the genetic association of specific genetic variants at different CYRIA gene loci with the development of NSOFCs in Kuwait.

METHODS

Four genetic variants (rs7552, rs3758249, rs3821949, and rs3917201) at four selected gene loci (CYRIA, FOXE1, MSX1, and TGFB3) were genotyped in a total of 240 DNA samples (patients (n = 114) and random controls (n = 126)) employing TaqMan® allele discrimination assay. For each variant and its genotype, the frequencies were determined and tested for Hardy-Weinberg Equilibrium. Genotype frequencies was compared between patients and controls using Pearson's test. Logistic regression analyses were employed to test for the associations of the four selected variants with the occurrence of NSOFCSs.

RESULTS

Significant differences in the distribution of genotypes between cases and controls, rs7552, rs3821949, and rs3917201 were found to have a positive association with NSOFCs. After adjusting for gender, the GG genotype of the rs7552 variant, the AG genotype of the rs3821949 variant, and the CC genotype of the rs3917201 variant showed nearly a two-fold increased risk of NSOFC (p < 0.05).

CONCLUSION

This study reports significant findings on the contribution and modest effect of CYRIA rs7552, MSX1 rs3821949, and TGFB3 rs3917201 in the development of NSOFCs. Our findings provide further evidence on the molecular mechanism and the role of the selected genes in NSOFCs.

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.

Immunohistochemical Evaluation of BARX1, DLX4, FOXE1, HOXB3, and MSX2 in Nonsyndromic Cleft Affected Tissue

Background

Nonsyndromic craniofacial clefts are relatively common congenital malformations which could create a significant negative effect on the health status and life quality of affected individuals within the pediatric population. Multiple cleft candidate genes and their coded proteins have been described with their possible involvement during cleft formation. Some of these proteins like Homeobox Protein BarH-like 1 (BARX1), Distal-Less Homeobox 4 (DLX4), Forkhead Box E1 (FOXE1), Homeobox Protein Hox-B3 (HOXB3), and Muscle Segment Homeobox 2 (MSX2) have been associated with the formation of craniofacial clefts. Understanding the pathogenetic mechanisms of nonsyndromic craniofacial cleft formation could provide a better knowledge in cleft management and could be a possible basis for development and improvement of cleft treatment options. This study investigates the presence of BARX1, DLX4, FOXE1, HOXB3, and MSX2 positive cells by using immunohistochemistry in different types of cleft-affected tissue while determining their possible connection with cleft pathogenesis process.

Materials and Methods

Craniofacial cleft tissue material was obtained during cleft-correcting surgery from patients with nonsyndromic craniofacial cleft diagnosis. Tissue material was gathered from patients who had unilateral cleft lip (n=36), bilateral cleft lip (n=13), and cleft palate (n=26). Control group (n=7) tissue material was received from individuals without any craniofacial clefts. The number of factor positive cells in the control group and patient group tissue was evaluated by using the semiquantitative counting method. Data was evaluated with the use of nonparametric statistical methods.

Results

Statistically significant differences were identified between the number of BARX1, FOXE1, HOXB3, and MSX2-containing cells in controls and cleft patient groups but no statistically significant difference was found for DLX4. Statistically significant correlations between the evaluated factors were also notified in cleft patient groups.

Conclusions

HOXB3 could be more associated with morphopathogenesis of unilateral cleft lip during postnatal course of the disorder. FOXE1 and BARX1 could be involved with both unilateral and bilateral cleft lip morphopathogenesis. The persistence of MSX2 in all evaluated cleft types could indicate its possible interaction within multiple cleft types. DLX4 most likely is not involved with postnatal cleft morphopathogenesis process.

Non-syndromic Cleft Palate: An Overview on Human Genetic and Environmental Risk Factors

The epithelial and mesenchymal cells involved in early embryonic facial development are guided by complex regulatory mechanisms. Any factor perturbing the growth, approach and fusion of the frontonasal and maxillary processes could result in orofacial clefts that represent the most common craniofacial malformations in humans. The rarest and, probably for this reason, the least studied form of cleft involves only the secondary palate, which is posterior to the incisive foramen. The etiology of cleft palate only is multifactorial and involves both genetic and environmental risk factors. The intention of this review is to give the reader an overview of the efforts made by researchers to shed light on the underlying causes of this birth defect. Most of the scientific papers suggesting potential environmental and genetic causes of non-syndromic cleft palate are summarized in this review, including genome-wide association and gene–environment interaction studies.

Association between forkhead box E1 polymorphisms and risk of non-syndromic cleft lip with or without cleft palate: A meta-analysis

OBJECTIVE

The purpose of the present work was to investigate the association between forkhead box E1 (FOXE1) and the risk of non-syndromic cleft lip with or without cleft palate (NSCL/P).

MATERIALS AND METHODS

Relevant studies were searched in several professional databases up to 31 July 2019. The pooled odds ratios (ORs) and 95% confidence intervals (95% CIs) were calculated using a fixed-effect model or a random-effect model to analyse the relationship between FOXE1 polymorphisms and NSCL/P.

RESULTS

A total of four single nucleotide polymorphisms (SNPs), including rs3758249, rs4460498, rs1443434 and rs10217225, were analysed. The overall findings showed that FOXE1 rs4460498 was statistically associated with NSCL/P (including cleft lip with or without cleft palate (CL/P) and cleft palate only (CPO)). Genotypes CC and CT of rs4460498 were significantly more closely correlated with NSCL/P (including CL/P and CPO) than genotype TT (NSCL/P: TT vs CC, OR = 0.630, P = .000; TT vs TC + CC, OR = 0.775, P = .020; CL/P: TT vs CC, OR = 0.664, P = .000; TT vs TC + CC, OR = 0.738, P = .006. CPO: TT vs CC, OR = 0.761, P = .027; TT vs TC + CC, OR = 0.792, P = .045). For rs10217225, only the TT genotype might have contributed to the elevated risk of CL/P (TT vs CC OR = 2.236, P = .000). The other FOXE1 polymorphisms were not associated with NSCLP, CL/P or CPO.

CONCLUSION

The meta-analysis provided confirmation that the polymorphism of FOXE1 rs10217225 was correlated with an increased risk of CL/P, and the polymorphism of FOXE1 rs4460498 was a protective factor for NSCL/P, including CLP and CPO.

A systematic review and meta-analysis on protective role of forkhead box E1 (FOXE1) polymorphisms in susceptibility to non-syndromic cleft lip/palate

INTRODUCTION

Several environmental and genetic factors have a role in the aetiology of non-syndromic cleft lip/palate (NSCL/P). This meta-analysis evaluated the association of rs3758249 and rs4460498 forkhead box E1 (FOXE1) polymorphisms with the NSCL/P risk.

MATERIALS AND METHODS

The Scopus, Cochrane Library, Web of Science, and PubMed databases were searched for articles published until March 2019. The analyses were performed by Review Manager 5.3 using the crude odds ratio (OR) and 95% confidence interval (CI) for a strong association between FOXE1 polymorphisms and the risk of NSCL/P.

RESULTS

Out of 161 articles retrieved from the databases, four case-control articles were involved in the meta-analysis. The pooled ORs of rs4460498 polymorphism based on allelic, homozygous, heterozygous, dominant, and recessive models were 0.74 (95% CI: 0.69, 0.80; P<0.00001), 0.43 (95% CI: 0.30, 0.61; P<0.00001), 0.66 (95% CI: 0.55, 0.80; P<0.0001), 0.66 (95% CI: 0.59, 0.73; P<0.00001), and 0.70 (95% CI: 0.60, 0.82; P<0.0001), respectively; whereas, the pooled OR of rs3758249 polymorphism were 0.86 (95% CI: 0.71, 1.04; P=0.12), 0.68 (95% CI: 0.57, 0.82; P<0.0001), 0.79 (95% CI: 0.57, 1.09; P=0.15), 0.79 (95% CI: 0.58, 1.08; P=0.14), and 0.80 (95% CI: 0.68, 0.95; P=0.010) for the afore-mentioned models, respectively.

CONCLUSIONS

The results showed that the T allele, TT, and CT genotypes of rs4460498 polymorphism were significantly associated with a decreased risk of NSCL/P; whereas, for rs3758249 polymorphism, only the AA genotype had a significant protective role in NSCL/P. Thus, FOXE1 is strongly associated with NSCL/P in the populations.

Association of SNP rs1867277 in FOXE1 Gene and Cleft Lip with or without Cleft Palate in a Han Chinese Population

BACKGROUND

The genetic factors causing cleft lip with or without cleft palate (CL ± P) are still unclear. The SNPs in FOXE1 gene were associated with CL ± P. However, the results have been inconsistent.

OBJECTIVE

We explored the associations of four SNPs in FOXE1 gene and CL ± P by a family based study.

MATERIALS AND METHODS

128 children with CL ± P and their parents were recruited. rs3758249 and rs1867277 were genotyped by high-resolution melting curve (HRM) method, whereas rs1443434 and rs907577 were genotyped by Sequenom MassARRAY® method. The software PLINK, FBAT and FAMHAP were used for analyzing data.

RESULTS

rs1867277 was associated with CL ± P (P_m = 0.0395). The patients were divided into two subgroups, individuals with cleft lip only and persons with cleft lip and palate. There were no associations in subgroup analyses.

CONCLUSION

We confirmed the association of FOXE1 gene and CL ± P by a family based study. For the first time, rs1867277 was significantly associated with CL ± P.

Prevalence and Parental Risk Factors for Speech Disability Associated with Cleft Palate in Chinese Children-A National Survey

Although the prevalence of oral clefts in China is among the highest worldwide, little is known about the prevalence of speech disability associated with cleft palate in Chinese children. The data for this study were collected from the Second China National Sample Survey on Disability, and identification of speech disability associated with cleft palate was based on consensus manuals. Logistic regression was used to estimate odds ratios (ORs) and 95% confidence intervals (CIs). A weighted number of 112,070 disabled children affected by cleft palate were identified, yielding a prevalence of 3.45 per 10,000 children (95% CI: 3.19–3.71). A history of speech disability in the mother (OR = 20.266, 95% CI 5.788–70.959, p < 0.0001), older paternal child-bearing age (OR = 1.061, 95% CI 1.017–1.108, p = 0.0065, per year increase in age), and lower parental education (maternal: OR = 3.424, 95% CI 1.082–10.837, p = 0.0363; paternal: OR = 2.923, 95% CI 1.245–6.866, p = 0.0138) were strongly associated with risk of speech disability associated with cleft palate in the offspring. Our results showed that maternal speech disability, older paternal child-bearing age, and lower levels of parental education were independent risk factors for speech disability associated with cleft palate for children in China. These findings may have important implications for health disparities and prevention.

Genetic variation of FOXE1 and risk for orofacial clefts in a California population

We investigated whether orofacial clefts are associated with polymorphic variation within and around FOXE1. This California population-based case control study focused on white Hispanic and white nonHispanic infants among which there were 262 infants with cleft lip with or without cleft palate (CL/P), 103 with cleft palate only (CPO), and 382 unaffected controls. These cases and controls were genotyped for 13 SNPs across 220 Kb at the FOXE1 Locus. We observed associations with multiple FOXE1 SNPs for CL/P and for CPO, especially for the Hispanic study population. Increased risks were associated with the more common allele for all SNPs tested. Our results implicate FOXE1 as an important locus whose polymorphic variation increases risks for all types of isolated clefts, and opens a new biological pathway to investigate in efforts to understand genetic factors underlying human clefting. © 2016 Wiley Periodicals, Inc.