Interaction between IRF6 and TGFA genes contribute to the risk of nonsyndromic cleft lip/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.

Identification of copy neutral loss of heterozygosity on chromosomes 1p, 1q, and 6p among nonsyndromic cleft lip and/or without cleft palate with hypodontia

BACKGROUND

Nonsyndromic cleft lip and/or without cleft palate (NSCL/P) with or without hypodontia is a common developmental aberration in humans and animals. This study aimed to identify the loss of heterozygosity (LOH) involved in hypodontia and NSCL/P pathogenesis.

METHODS

This is a cross-sectional study that conducted genome-wide copy number analysis using CytoScan 750K array on salivary samples from Malay subjects with NSCL/P with or without hypodontia aged 7-13 years. To confirm the significant results, simple logistic regression was employed to conduct statistical data analysis using SPSS software.

RESULTS

The results indicated the most common recurrent copy neutral LOH (cnLOH) observed at 1p33-1p32.3, 1q32.2-1q42.13 and 6p12.1-6p11.1 loci in 8 (13%), 4 (7%), and 3 (5%) of the NSCL/P subjects, respectively. The cnLOHs at 1p33-1p32.3 (D1S197), 1q32.2-1q42.13 (D1S160), and 6p12.1-6p11.1 (D1S1661) were identified observed in NSCL/P and noncleft children using microsatellite analysis markers as a validation analysis. The regions affected by the cnLOHs at 1p33-1p32.3, 1q32.2-1q42.13, and 6p12.1-6p11.1 loci contained selected genes, namely FAF1, WNT3A and BMP5, respectively. There was a significant association between the D1S197 (1p33-32.3) markers containing the FAF1 gene among NSCL/P subjects with or without hypodontia compared with the noncleft subjects (p-value = 0.023).

CONCLUSION

The results supported the finding that the genetic aberration on 1p33-32.3 significantly contributed to the development of NSCL/P with or without hypodontia. These results have an exciting prospect in the promising field of individualized preventive oral health care.

Orofacial Clefts: Genetics of Cleft Lip and Palate

Orofacial clefting is considered one of the commonest birth defects worldwide. It presents as cleft lip only, isolated cleft palate or cleft lip and palate. The condition has a diverse genetic background influenced by gene-gene and gene-environment interaction, resulting in two main types, syndromic and nonsyndromic orofacial clefts. Orofacial clefts lead to significant physiological difficulties that affect feeding, speech and language development and other developmental aspects, which results in an increased social and financial burden on the affected individuals and their families. The management of cleft lip and palate is solely based on following a multidisciplinary team approach. In this narrative review article, we briefly summarize the different genetic causes of orofacial clefts and discuss some of the common syndromes and the approach to the management of orofacial clefts.

The IRF6 AP-2α binding site polymorphism relate to the severity of non-syndromic orofacial cleft of Indonesian patients

BACKGROUND

IRF6 AP-2α binding site polymorphism is known as IRF6 rs642961. It has been associated with a nonsyndromic orofacial cleft (NS OFC). This study aimed to determine the IRF6 rs642961 as a risk factor associated with NS OFC and its phenotypes.

METHODS

The case-control design used for 264 subjects consists of 158 NS CLP subjects (42 CU CLP, 34 CB CLP, 33 CLO, 49 CPOs) and 106 healthy controls. The DNA is extracted from venous blood. The segment of IRF6 rs642961 amplified by polymerase chain reaction (PCR) followed by restriction fragment length of polymorphisms (RFLPs) used the MspI digestion enzyme. The qPCR method to identify the mRNA expression levels of the IRF6 gene rs642961 was analyzed by the Livak method.

RESULTS

The study results show that in NS CB CLP phenotype as the most severe phenotype of NS OFC, the Odds Ratio (OR) of A mutant allele was 5.094 (CI=1.456-17.820; P=0.011) and the OR of AA homozygous mutant genotype was 13.481 (CI=2.648-68.635; P=0.001). There are different levels of mRNA expression changes from NS OFC and its phenotypes. It is substantial among the 2^-ΔΔCt and the group of AA, GA, and GG genotypes (P<0.05); in the NS CPO phenotype, it shows IRF6 mRNA under-expression in GA, AA genotypes while in other phenotypes it shows IRF6 mRNA overexpression.

CONCLUSIONS

The IRF6 AP-2α binding site polymorphism is strongly associated with the severity of NS OFC, and this polymorphism has a functional role in affecting IRF6 mRNA expression that is variable in each phenotype.

Dental anomalies in non-syndromic orofacial clefts: a clinical approach

Orofacial clefts (OC) represent some of the most common congenital malformations and affect multiple structures in the craniofacial region. There are a wide range of morphological OC types within the spectrum of both non-syndromic OC (NSOC) and syndromic OC, including cleft lip (CL); cleft lip and palate (CLP) and cleft palate (CP). Here, we describe the most frequent dental anomalies seen in the permanent dentition of individuals with NSOC, comparing them between the three main cleft types (CL, CLP, CP). We present these findings from the perspective of prevalence relating to each anomaly, as well as the clinical characteristics and potential impact in the rehabilitation process. The investigation of dental anomalies associated with NSOC is important, helping to expand the phenotypic characterization of NSOC, improve initial diagnostics and establish appropriate rehabilitation and counseling.

Monogenic Versus Multifactorial Inheritance in the Development of Isolated Cleft Palate: A Whole Genome Sequencing Study

Craniofacial morphogenesis is highly complex, as is the anatomical region involved. Errors during this process, resulting in orofacial clefts, occur in more than 400 genetic syndromes. Some cases of cleft lip and/or palate (CLP) are caused by mutations in single genes; however, complex interactions between genetic and environmental factors are considered to be responsible for the majority of non-syndromic CLP development. The aim of the current study was to identify genetic risk factors in patients with isolated cleft palate (CP) by whole genome sequencing. Patients with isolated CP (n = 30) recruited from the Riga Cleft Lip and Palate Centre, Institute of Stomatology, Riga, were analyzed by whole genome sequencing. Pathogenic or likely pathogenic variants were discovered in genes associated with CP (TBX22, COL2A1, FBN1, PCGF2, and KMT2D) in five patients; hence, rare disease variants were identified in 17% of patients with non-syndromic isolated CP. Our results were relevant to routine genetic counselling practice and genetic testing recommendations. Based on our data, we propose that all newborns with orofacial clefts should be offered genetic testing, at least for a panel of known CLP genes. Only if the results are negative and there is no suggestive family history or additional clinical symptoms (which would support additional exome or genome-wide investigation), should multifactorial empiric recurrence risk prediction tools be applied for families.

Dental caries and periodontitis risk factors in cleft lip and palate patients

Cleft lip and palate (CLP) is the most common congenital facial malformation and has a significant developmental, physical, and psychological impact on those with the deformity and their families. Risk factors contributing to CLP may conclude as genetic factors and environmental factors. The anatomical and morphological abnormalities related to CLP are favorable for dental plaque accumulation on the tooth surface. Therefore, patients with CLP undergo poorer oral hygiene and higher susceptibility to dental caries and periodontitis. In this review, we aim to conclude and update probable causes underlying the association between CLP and poor oral health and provide novel ideas of targeted early prevention for such oral diseases.

Identification of Copy Number Variation Among Nonsyndromic Cleft Lip and or Without Cleft Palate With Hypodontia: A Genome-Wide Association Study

Nonsyndromic cleft lip and or without cleft palate (NSCL/P) with the hypodontia is a common developmental abnormality in humans and animals. This study identified the genetic aberration involved in both NSCL/P and hypodontia pathogenesis. A cross-sectional study using genome-wide study copy number variation-targeted CytoScan 750K array carried out on salivary samples from 61 NSCL/P and 20 noncleft with and without hypodontia Malay subjects aged 7-13 years old. Copy number variations (CNVs) of SKI and fragile histidine triad (FHIT) were identified in NSCL/P and noncleft children using quantitative polymerase chain reaction (qPCR) as a validation analysis. Copy number calculated (CNC) for each gene determined with Applied Biosystems CopyCaller Software v2.0. The six significant CNVs included gains (12q14.3, 15q26.3, 1p36.32, and 1p36.33) and losses (3p14.2 and 4q13.2) in NSCL/P with hypodontia patients compared with the NSCL/P only. The genes located in these regions encoded LEMD3, IGF1R, TP73, SKI, FHIT, and UGT2β15. There were a significant gain and loss of both SKI and FHIT copy number in NSCL/P with hypodontia compared with the noncleft group (p < 0.05). The results supported that CNVs significantly furnish to the development of NSCL/P with hypodontia.

Scandcleft randomized trials of primary surgery for unilateral cleft lip and palate: impact of maxillary dental agenesis on craniofacial growth and dental arch relationship in 8 year olds

BACKGROUND

It is suggested that dental agenesis affects maxillary protrusion and dental arch relationship in children with unilateral cleft lip and palate (UCLP). In addition, an association between the need for orthognathic surgery and dental agenesis is reported.

AIM

The aim was to study the impact of maxillary dental agenesis on craniofacial growth and dental arch relationship in 8-year-old children with UCLP.

SUBJECTS AND METHODS

The sample consisted of individuals with UCLP from Scandcleft randomized trials. The participants had available data from diagnosis of maxillary dental agenesis as well as cephalometric measurements (n = 399) and GOSLON assessment (n = 408) at 8 years of age.

RESULTS

A statistically significant difference was found for ANB between individuals with agenesis of two or more maxillary teeth (mean 1.52°) in comparison with those with no or only one missing maxillary tooth (mean 3.30° and 2.70°, respectively). Mean NSL/NL was lower among individuals with agenesis of two or more maxillary teeth (mean 9.90°), in comparison with individuals with no or one missing maxillary tooth (mean 11.46° and 11.45°, respectively). The number of individuals with GOSLON score 4-5 was 47.2% in the group with two or more missing maxillary teeth and 26.1% respectively 26.3% in the groups with no or one missing maxillary tooth. No statistically significant difference was found in the comparison between individuals with no agenesis or with agenesis solely of the cleft-side lateral.

CONCLUSION

Maxillary dental agenesis impacts on craniofacial growth as well as dental arch relationship and should be considered in orthodontic treatment planning.

Tooth agenesis: What do we know and is there a connection to cancer?

Like all developmental processes, odontogenesis is highly complex and dynamically regulated, with hundreds of genes co-expressed in reciprocal networks. Tooth agenesis (missing one or more/all teeth) is a common human craniofacial anomaly and may be caused by genetic variations and/or environmental factors. Variants in PAX9, MSX1, AXIN2, EDA, EDAR, and WNT10A genes are associated with tooth agenesis. Currently, variants in ATF1, DUSP10, CASC8, IRF6, KDF1, GREM2, LTBP3, and components and regulators of WNT signaling WNT10B, LRP6, DKK, and KREMEN1 are at the forefront of interest. Due to the interconnectedness of the signaling pathways of carcinogenesis and odontogenesis, tooth agenesis could be a suitable marker for early detection of cancer predisposition. Variants in genes associated with tooth agenesis could serve as prognostic or therapeutic targets in cancer. This review aims to summarize existing knowledge of development and clinical genetics of teeth. Concurrently, the review proposes possible approaches for future research in this area, with particular attention to roles in monitoring, early diagnosis and therapy of tumors associated with defective tooth development.

Gene-environment interaction in molar-incisor hypomineralization

Molar incisor hypomineralization (MIH) is an enamel condition characterized by lesions ranging in color from white to brown which present rapid caries progression, and mainly affects permanent first molars and incisors. These enamel defects usually occur when there are disturbances during the mineralization or maturation stage of amelogenesis. Both genetic and environmental factors have been suggested to play roles in MIH’s development, but no conclusive risk factors have shown the source of the disease. During head and neck development, the interferon regulatory factor 6 (IRF6) gene is involved in the structure formation of the oral and maxillofacial regions, and the transforming growth factor alpha (TGFA) is an essential cell regulator, acting during proliferation, differentiation, migration and apoptosis. In this present study, it was hypothesized that these genes interact and contribute to predisposition of MIH. Environmental factors affecting children that were 3 years of age or older were also hypothesized to play a role in the disease etiology. Those factors included respiratory issues, malnutrition, food intolerance, infection of any sort and medication intake. A total of 1,065 salivary samples from four different cohorts were obtained, and DNA was extracted from each sample and genotyped for nine different single nucleotide polymorphisms. Association tests and logistic regression implemented in PLINK were used for analyses. A potential interaction between TGFA rs930655 with all markers tested in the cohort from Turkey was identified. These interactions were not identified in the remaining cohorts. Associations (p<0.05) between the use of medication after three years of age and MIH were also found, suggesting that conditions acquired at the age children start to socialize might contribute to the development of MIH.

Genetic variation of IRF6 and TGFA genes in an HIV-exposed newborn with non-syndromic cleft lip palate

We present here the first reported case of a non-syndromic cleft lip and palate (NSCLP) in an HIV-exposed newborn of a mother on antiretroviral therapy (ART) in Indonesia. Genetic testing was performed to confirm a suspected genetic condition. Genomic DNA was extracted from the blood, and genetic variations of the interferon regulatory factor 6 (IRF6) rs642961 (Mspl) (G>A) and transforming growth factor alpha (TGFA) BamHI (rs11466297, A>C) and RsaI (rs3732248, C>T) were performed by PCR-RFLP and IRF6 gene analysis by PCR sequencing. Genotyping of DNA sequence variants in the IRF6 gene showed both parents had genotype GA, while the child had genotype GG (genotype wild type). There was no difference observed in the TGFA BamHI gene variant between the child and her mother and father that were wild-type polymorphisms (normal), while the Rsa1 polymorphisms of them were heterozygotes. A genetic variant of IRF6 might be a protective factor for NSCLP, while Rsa1 gene variant (A) allele can be considered to be the risk factor associated with NSCLP development. This case report also highlights the possible etiologic role of ART in NSCLP; therefore, early control of adverse effects of ART might be an important factor in decreasing the incidence of the congenital anomalies in HIV-infected children.

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.

[Evaluating the effect of WNT pathway genes considering interactions on the risk of non-syndromic oral clefts among Chinese populations]

OBJECTIVE

In this study, we used genome-wide association study (GWAS) data to explore whether WNT pathway genes were associated with non-syndromic oral clefts (NSOC) considering gene-gene interaction and gene-environment interaction.

METHODS

We conducted the analysis using 806 non-syndromic cleft lip with or without cleft palate (NSCL/P) case-parent trios and 202 non-syndromic cleft palate (NSCP) case-parent trios among Chinese populations selected from an international consortium established for a GWAS of non-syndromic oral clefts. Genotype data and maternal environmental exposures were collected through DNA samples and questionnaires. Conditional Logistic regression models were adopted to explore gene-gene interaction and gene-environment in teraction using trio package in R software. The threshold of significance level was set as 3.47×10^-4 using Bonferroni correction.

RESULTS

A total of 144 single nucleotide polymorphisms (SNPs) in seven genes passed the quality control process in NSCL/P trios and NSCP trios, respectively. Totally six pairs of SNPs interactions showed statistically significant SNP-SNP interaction (P < 3.47×10^-4) after Bonferroni correction, which were rs7618735 (WNT5A) and rs10848543 (WNT5B), rs631948 (WNT11) and rs556874 (WNT5A), and rs631948 (WNT11) and rs472631 (WNT5A) among NSCL/P trios; rs589149 (WNT11) and rs4765834 (WNT5B), rs1402704 (WNT11) and rs358792 (WNT5A), and rs1402704 (WNT11) and rs358793 (WNT5A) among NSCP trios, respectively. In addition, no significant result was found for gene-environment interaction analysis in both of the NSCL/P trios and NSCP trios.

CONCLUSION

Though this study failed to detect significant association based on gene-environment interactions of seven WNT pathway genes and the risk of NSOC, WNT pathway genes may influence the risk of NSOC through potential gene-gene interaction.

Genetic variants in tooth agenesis-related genes might be also involved in tooth size variations

OBJECTIVE

The present study aimed to evaluate if genetic variants in PAX9, MSX1, TGFα, FGF3, FGF10, FGF13, GLI2 and GLI3 are involved in TS of permanent teeth.

MATERIALS AND METHODS

Pretreatment dental records from orthodontic patients were assessed prior to recruitment. Patients with tooth agenesis and congenital anomalies (including oral cleft) and/or syndromes were excluded. Dental casts were used to measure the maximum crown dimensions of all fully erupted permanent teeth except second and third molars in mesiodistal direction. Teeth with caries, occlusal wear, mesiodistal restorations, and obvious deformities were not evaluated. Genomic DNA samples were used for genotyping. The allelic discrimination of 13 genetic variants was performed. The associations between TS and genotype were analyzed by linear regression, adjusted by gender at a significance level of p ≤ 0.05.

RESULTS

Genetic polymorphisms in the tooth agenesis-related genes studied here were associated with increased and decreased TS, in both maxilla and mandible (p < 0.05).

CONCLUSION

This study reported associations of novel tooth agenesis-related gene variants with permanent tooth size variations.

CLINICAL RELEVANCE

The presence of some genetic variants could allow the prediction of permanent tooth size.

Evidence for craniofacial enhancer variation underlying nonsyndromic cleft lip and palate

Nonsyndromic cleft lip with or without cleft palate (NSCLP) is a common birth defect for which only ~ 20% of the underlying genetic variation has been identified. Variants in noncoding regions have been increasingly suggested to contribute to the missing heritability. In this study, we investigated whether variation in craniofacial enhancers contributes to NSCLP. Candidate enhancers were identified using VISTA Enhancer Browser and previous publications. Prioritization was based on patterning defects in knockout mice, deletion/duplication of craniofacial genes in animal models and results of whole exome/whole genome sequencing studies. This resulted in 20 craniofacial enhancers to be investigated. Custom amplicon-based sequencing probes were designed and used for sequencing 380 NSCLP probands (from multiplex and simplex families of non-Hispanic white (NHW) and Hispanic ethnicities) using Illumina MiSeq. The frequencies of identified variants were compared to ethnically matched European (CEU) and Los Angeles Mexican (MXL) control genomes and used for association analyses. Variants in mm427/MSX1 and hs1582/SPRY1 showed genome-wide significant association with NSCLP (p ≤ 6.4 × 10^-11). In silico analysis showed that these enhancer variants may disrupt important transcription factor binding sites. Haplotypes involving these enhancers and also mm435/ABCA4 were significantly associated with NSCLP, especially in NHW (p ≤ 6.3 × 10^-7). Importantly, groupwise burden analysis showed several enhancer combinations significantly over-represented in NSCLP individuals, revealing novel NSCLP pathways and supporting a polygenic inheritance model. Our findings support the role of craniofacial enhancer sequence variation in the etiology of NSCLP.

Scandcleft randomized trials of primary surgery for unilateral cleft lip and palate: dental anomalies in 8-year olds

BACKGROUND

Children born with unilateral cleft lip and palate (UCLP) are reported to display several dental anomalies including agenesis, supernumeraries, as well as variations in dental size, shape, and path of eruption. The extensive sample of individuals with UCLP included in the Scandcleft randomized control trials offers the opportunity to study more rare conditions, which is seldom possible with limited samples.

OBJECTIVES

The aim was to study dental anomalies at 8 years of age in children born with UCLP included in the Scandcleft randomized control trials.

METHODS

Panoramic and intraoral radiographs from 425 individuals (279 males and 146 females) with a mean age of 8.1 years were assessed by four orthodontists regarding dental anomalies.

RESULTS

Agenesis was found in 52.5 per cent and supernumerary teeth in 16.9 per cent of the participants. The cleft lateral was missing in 43.8 per cent and was found peg shaped in 44.7 per cent. The distribution of ectopic eruption was 14.6 per cent, mainly affecting maxillary first molars, while transposition was found in 3.4 per cent of the individuals. In addition, infraocclusion of one or several primary molars was registered in 7.2 per cent of the participants.

CONCLUSION

We conclude that 8-year-old children born with UCLP display multiple dental anomalies. The Scandcleft sample allowed rarely studied conditions such as infraocclusion of primary molars and transposition to be studied in children born with UCLP.

TRIAL REGISTRATION

ISRCTN29932826.

PBX-WNT-P63-IRF6 pathway in nonsyndromic cleft lip and palate

Nonsyndromic cleft lip and palate (NSCLP) is one of the most common craniofacial anomalies in humans, affecting more than 135,000 newborns worldwide. NSCLP has a multifactorial etiology with more than 50 genes postulated to play an etiologic role. The genetic pathway comprised of Pbx-Wnt-p63-Irf6 genes was shown to control facial morphogenesis in mice and proposed as a regulatory pathway for NSCLP. Based on these findings, we investigated whether variation in PBX1, PBX2, and TP63, and their proposed interactions were associated with NSCLP. Fourteen single nucleotide variants (SNVs) in/nearby PBX1, PBX2, and TP63 were genotyped in 780 NSCLP families of nonHispanic white (NHW) and Hispanic ethnicities. Family-based association tests were performed for individual SNVs stratified by ethnicity and family history of NSCLP. Gene-gene interactions were also tested. A significant association was found for PBX2 rs3131300 and NSCLP in combined Hispanic families (p = .003) while nominal association was found for TP63 rs9332461 in multiplex Hispanic families (p = .005). Significant haplotype associations were observed for PBX2 in NHW (p = .0002) and Hispanic families (p = .003), and for TP63 in multiplex Hispanic families (.003). An independent case-control group was used to validate findings, and significant associations were found with PBX1 rs6426870 (p = .007) and TP63 rs9332461 (p = .03). Gene-gene interactions were detected between PBX1/PBX2/TP63 with IRF6 in NHW families, and between PBX1 with WNT9B in both NHW and Hispanic families (p < .0018). This study provides the first evidence for a role of PBX1 and PBX2, additional evidence for the role of TP63, and support for the proposed PBX-WNT-TP63-IRF6 regulatory pathway in the etiology of NSCLP.

IRF6 polymorphisms in Brazilian patients with non-syndromic cleft lip with or without palate

INTRODUCTION

Non-syndromic orofacial clefts have a complex etiology due to the contribution from both genetic and environmental risk factors, as well as the interaction between them. Among the more than 15 susceptibility loci for non-syndromic orofacial clefts with considerable statistical and biological support, the IRF6 is the most validated gene by the majority of studies. Nonetheless, in genetically heterogeneous populations such as Brazilian, the confirmation of association between non-syndromic orofacial clefts and IRF6 common variants is not a consolidated fact and unrecognized IRF6 variants are poorly investigated.

OBJECTIVE

The aim of this study was to investigate the association of IRF6 polymorphisms with non-syndromic orofacial clefts development in a population from northeast Brazil.

METHODS

Blood samples of 186 non-syndromic orofacial clefts patients and 182 controls from Rio Grande do Norte, Brazil, were obtained to analyze IRF6 polymorphisms (rs2235371, rs642961, rs2236907, rs861019, and rs1044516) by real-time polymerase chain reaction. Non-syndromic orofacial clefts patients were classified in cleft lip and palate, cleft palate only and cleft lip only groups.

RESULTS

The genotype and allele frequencies of single nucleotide polymorphism rs2235371 in IRF6 showed significant differences in patients with cleft palate when compared to the controls, whereas no association was shown between rs642961, rs2236907, rs861019, and rs1044516 and non-syndromic orofacial clefts.

CONCLUSION

The association found between rs2235371 and isolated cleft palate should be interpreted with caution due to the low number of individuals investigated, and more studies with larger sample size are needed to confirm these association. In addition, there is a lack of association of the rs642961, rs2236907 and rs861019 polymorphisms with non-syndromic orofacial clefts susceptibility.

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.

Exploring the interaction between FGF Genes and T-box genes among chinese nonsyndromic cleft lip with or without cleft palate case-parent trios

Nonsyndromic cleft lip with or without cleft palate (NSCL/P) is a common birth defect. Genetic variants causing syndromic orofacial clefts can also contribute to the etiology of NSCL/P. The purpose of the present study was to explore gene-gene (G × G) interaction using common single nucleotide polymorphic (SNP) markers in fibroblast growth factor (FGF) family and its receptors and T-box genes, which were associated with syndromic orofacial clefts. Our study was conducted in 806 Chinese NSCL/P case-parent trios drawn from an international consortium. A total of 252 SNPs in FGF8, FGF10, FGFR1, FGFR2, and TBX5 passed the quality control criteria and were included in the analysis. The interactions between SNPs in different genes were assessed using Cordell's method, which fitted a conditional logistic regression model. The analysis was performed using the R-package trio (Version 3.8.0). Bonferroni correction was used to adjust for multiple comparisons, and the overall significance threshold was set as P = 1.98 × 10^-4 (0.05/252). Conditional logistic regression revealed the most significant interaction between rs2330542 in FGF10 and rs1946295 in TBX5, which remained significant (P = 9.63 × 10^-6 ) after Bonferroni correction. The relative risk of allele C in rs2330542 (FGF10) was 1.02 (95%CI 0.81-1.28), while the relative risk was 1.42 (95%CI 1.03-1.97) when the exposure was a combination of allele C in rs2330542 and allele A in rs1946295 (TBX5). Our findings confirmed the importance of considering G × G interaction when exploring the genetic risk factors of NSCL/P. Further investigations are warranted to validate the potential interaction and reveal the biological function of FGF10/FGFR2/TBX5. Environ. Mol. Mutagen. 2019. © 2019 Wiley Periodicals, Inc.

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.

Evaluating the effect of nicotinic cholinergic receptor genes on the risk of nonsyndromic cleft lip with or without cleft palate

OBJECTIVE

Multiple studies have suggested nonsyndromic cleft lip with or without cleft palate (NSCL/P), and lung cancer may have common genetic etiology. Previous studies have showed genetic variants in nicotinic cholinergic receptor genes (CHRNs) may influence risk of lung cancer. We aimed to explore the effect of CHRNs on risk of NSCL/P considering gene-gene (GxG) interaction for these genes.

SUBJECTS AND METHODS

We selected 120 markers in 14 CHRNs to test for GxG interaction using 806 Chinese case-parent trios recruited from an international consortium established for a GWAS of oral clefts.

RESULTS

Totally, two pairs of SNPs yielded significant GxG interactions after Bonferroni correction (rs935865 and rs2337980 with p = 4.04 × 10^-5 , rs2741335 and rs3743077 with p = 4.80 × 10^-4 ), and these pairwise interactions were confirmed in permutation tests. In addition, the relative risk (RR) of the putative interaction between rs935865 and rs2337980 was 1.10 (95% CI: 0.92~1.31).

CONCLUSIONS

While the single SNP association and the gene-environment interaction analysis of 14 CHRN genes yielded no signal, this study did demonstrate the importance of considering potential GxG interaction for exploring etiology of NSCL/P. This study suggests an important role for particular combinations of SNPs in CHRN genes in influencing risk to NSCL/P, which needs further study.

Evidence of interaction between genes in the folate/homocysteine metabolic pathway in controlling risk of non-syndromic oral cleft

OBJECTIVE

Little consistent evidence is available for the association between the risk of non-syndromic cleft lip with or without cleft palate (NSCL/P) and any of the individual genes in the folate/homocysteine metabolic pathway. We investigated the genes in the folate pathway to further clarify its potential influence on the risk of NSCL/P considering gene-gene (G×G) interaction.

SUBJECTS AND METHODS

We selected markers in 18 genes from the pathway and applied Cordell's method to test for G×G interaction using 1,908 NSCL/P case-parent trios ascertained in an international consortium where a genomewide association study (GWAS) of oral clefts was conducted.

RESULTS

We found intriguing signals among Asian and European ancestry groups for G×G interaction between markers in betaine-homocysteine methyltransferase gene (BHMT/BHMT2) and dimethylglycine dehydrogenase gene (DMGDH) attaining genomewide significance. In the pooled data, the top significant interaction was found between rs13158309 (BHMT) and rs10514154 (DMGDH, p = 1.45 × 10^-12 ).

CONCLUSIONS

Our study illustrated the importance of taking into account potential G×G interaction for genetic association analysis in NSCL/P, and this study suggested both BHMT/BHMT2 and DMGDH should be considered as candidate genes for NSCL/P in future studies.

Gene-gene interaction between MSX1 and TP63 in Asian case-parent trios with nonsyndromic cleft lip with or without cleft palate

BACKGROUND

Small ubiquitin-like modification, also known as sumoylation, is a crucial post-translational regulatory mechanisms involved in development of the lip and palate. Recent studies reported two sumoylation target genes, MSX1 and TP63, to have achieved genome-wide level significance in tests of association with nonsyndromic clefts. Here, we performed a candidate gene analysis considering gene-gene and gene-environment interaction for SUMO1, MSX1, and TP63 to further explore the etiology of nonsyndromic cleft lip with or without cleft palate (NSCL/P).

METHODS

A total of 130 single-nucleotide polymorphisms (SNPs) in or near SUMO1, MSX1, and TP63 was analyzed among 1,038 Asian NSCL/P trios ascertained through an international consortium. Conditional logistic regression models were used to explore gene-gene (G × G) and gene-environment (G × E) interaction involving maternal environmental tobacco smoke and multivitamin supplementation. Bonferroni correction was used for G × E analysis and permutation tests were used for G × G analysis.

RESULTS

While transmission disequilibrium tests and gene-environment interaction analysis showed no significant results, we did find signals of gene-gene interaction between SNPs near MSX1 and TP63. Three pairwise interactions yielded significant p values in permutation tests (rs884690 and rs9290890 with p = 9.34 × 10^-5 and empirical p = 1.00 × 10^-4 , rs1022136 and rs4687098 with p = 2.41 × 10^-4 and empirical p = 2.95 × 10^-4 , rs6819546 and rs9681004 with p = 5.15 × 10^-4 and empirical p = 3.02 × 10^-4 ).

CONCLUSION

Gene-gene interaction between MSX1 and TP63 may influence the risk of NSCL/P in Asian populations. Our study provided additional understanding of the genetic etiology of NSCL/P and underlined the importance of considering gene-gene interaction in the etiology of this common craniofacial malformation.

Interferon Regulatory Factor 6 Is Necessary for Salivary Glands and Pancreas Development

Interferon regulatory factor 6 ( IRF6) acts as a tumor suppressor and controls cell differentiation in ectodermal and craniofacial tissues by regulating expression of target genes. Haploinsufficiency of IRF6 causes Van der Woude and popliteal pterygium syndrome, 2 syndromic forms of cleft lip and palate. Around 85% of patients with Van der Woude express pits on the lower lip that continuously or intermittently drain saliva, and patients with the common cleft lip and palate have a higher prevalence of dental caries and gingivitis. This study aims to identify the role of IRF6 in development of exocrine glands, specifically the major salivary glands. Our transgenic mouse model that expresses LacZ reporter under the control of the human IRF6 enhancer element showed high expression of IRF6 in major and minor salivary glands and ducts. Immunostaining data also confirmed the endogenous expression of IRF6 in the developing ductal, serous, and mucous acinar cells of salivary glands. As such, we hypothesized that Irf6 is important for proper development of salivary glands and potentially other exocrine glands. Loss of Irf6 in mice causes an increase in the proliferation level of salivary cells, disorganized branching morphogenesis, and a lack of differentiated mucous acinar cells in submandibular and sublingual glands. Expression and localization of the acinar differentiation marker MIST1 were altered in Irf6-null salivary gland and pancreas. The RNA-Seq analysis demonstrated that 168 genes are differentially expressed and confer functions associated with transmembrane transporter activity, spliceosome, and transcriptional regulation. Furthermore, expression of genes involved in the EGF pathway-that is, Ereg, Ltbp4, Matn1, Matn3, and Tpo-was decreased at embryonic day 14.5, while levels of apoptotic proteins were elevated at postnatal day 0. In conclusion, our data report a novel role of Irf6 in exocrine gland development and support a rationale for performing exocrine functional tests for patients with IRF6-damaging mutations.

Gene-gene interaction of single nucleotide polymorphisms in 16p13.3 may contribute to the risk of non-syndromic cleft lip with or without cleft palate in Chinese case-parent trios

Non-syndromic cleft lip with or without cleft palate (NSCL/P) is a common birth defect with a complex and heterogeneous etiology. A recent genome-wide association study (GWAS) among Chinese populations has identified a new region at 16p13.3 as being associated with NSCL/P, which requires further replication. Here, we attempted to replicate and further clarify the genetic association between this region and NSCL/P, as well as testing for potential gene-gene (G × G) and gene-environment (G × E) interactions. We conducted transmission disequilibrium tests on 69 single nucleotide polymorphisms (SNPs) mapping to 16p13.3 among 806 Chinese case-parent trios ascertained through an international consortium where a GWAS of oral clefts was conducted. G × G, as well as G × E interactions involving maternal environmental tobacco smoke (ETS) and multivitamin supplementation, were explored using conditional logistic regression model. We applied Cordell's method as implemented in the R package TRIO to test for possible interactions. While no SNPs showed evidence of linkage and association with NSCL/P after Bonferroni correction, we found signals of G × G interactions between SNPs in 16p13.3. Nine pairs of SNP-SNP interactions attained significance after Bonferroni correction, among which the most significant interaction was found between rs2072346 (ADCY9) and rs11646137 (intergenic region, P = 7.2 × 10^-5 ). Linkage disequilibrium (LD) analysis revealed only low level of LD between these SNPs. This study failed to confirm the significant association between SNPs within 16p13.3 and the risk of NSCL/P, but underlined the importance of taking into account potential G × G interactions for the genetic association analysis of NSCL/P.

Tooth agenesis and orofacial clefting: genetic brothers in arms?

Tooth agenesis and orofacial clefts represent the most common developmental anomalies and their co-occurrence is often reported in patients as well in animal models. The aim of the present systematic review is to thoroughly investigate the current literature (PubMed, EMBASE) to identify the genes and genomic loci contributing to syndromic or non-syndromic co-occurrence of tooth agenesis and orofacial clefts, to gain insight into the molecular mechanisms underlying their dual involvement in the development of teeth and facial primordia. Altogether, 84 articles including phenotype and genotype description provided 9 genomic loci and 26 gene candidates underlying the co-occurrence of the two congenital defects: MSX1, PAX9, IRF6, TP63, KMT2D, KDM6A, SATB2, TBX22, TGFα, TGFβ3, TGFβR1, TGFβR2, FGF8, FGFR1, KISS1R, WNT3, WNT5A, CDH1, CHD7, AXIN2, TWIST1, BCOR, OFD1, PTCH1, PITX2, and PVRL1. The molecular pathways, cellular functions, tissue-specific expression and disease association were investigated using publicly accessible databases (EntrezGene, UniProt, OMIM). The Gene Ontology terms of the biological processes mediated by the candidate genes were used to cluster them using the GOTermMapper (Lewis-Sigler Institute, Princeton University), speculating on six super-clusters: (a) anatomical development, (b) cell division, growth and motility, (c) cell metabolism and catabolism, (d) cell transport, (e) cell structure organization and (f) organ/system-specific processes. This review aims to increase the knowledge on the mechanisms underlying the co-occurrence of tooth agenesis and orofacial clefts, to pave the way for improving targeted (prenatal) molecular diagnosis and finally to reflect on therapeutic or ultimately preventive strategies for these disabling conditions in the future.

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.

Brazilian multicenter study of association between polymorphisms in CRISPLD2 and JARID2 and non-syndromic oral clefts

BACKGROUND

Variants in the cysteine-rich secretory protein LCCL domain containing 2 gene (CRISPLD2) and in the jumonji, AT-rich interaction domain 2 gene (JARID2) were previously shown to influence non-syndromic oral cleft susceptibility. Herein, we performed a case-control study to examine the potential association of single-nucleotide polymorphisms (SNPs) in CRISPLD2 and JARID2 with non-syndromic cleft lip and/or palate (NSCL/P) in the Brazilian population. Given the ethnicity-dependent genetic predisposition to NSCL/P, we performed a structured analysis taking into account the genomic ancestry variation of each individual.

METHODS

Four SNPs in CRISPLD2 (rs1546124, rs8061351, rs2326398, and rs4783099) and four in JARID2 (rs915344, rs2299043, rs2237138, and rs2076056), that were previously reported to be associated with NSCL/P, were genotyped in 785 Brazilian patients with NSCL/P (549 with cleft lip with or without cleft palate-NSCL ± P, and 236 with cleft palate only-NSCPO) and 693 unaffected Brazilian controls. Genomic ancestry was assessed with a set of 40 biallelic short insertion/deletion variants previously validated as ancestry informative markers of the Brazilian population.

RESULTS

After adjustment of ancestry variations, allelic analysis revealed marginal associations between the CRISPLD2 rs4783099 T allele and increased risk for NSCPO (OR: 1.31, 95% CI: 1.05-1.62, P = 0.01) and between JARID2 rs2237138 and decreased NSCL ± P risk (OR: 0.80, 95% CI: 0.67-0.97, P = 0.02). Haplotype analysis indicated a lack of association between JARID2 haplotypes and non-syndromic oral cleft risk.

CONCLUSIONS

Our results suggest that CRISPLD2 rs4783099 may represent a risk factor for NSCPO while JARID2 rs2237138 shows a protective effect against NSCL ± P in the Brazilian population.

Gene Interactions Provide Evidence for Signaling Pathways Involved in Cleft Lip/Palate in Humans

Nonsyndromic cleft lip with or without cleft palate (NSCL/P) is a common craniofacial birth defect that has a complex etiology. Genome-wide association studies have recently identified new loci associated with NSCL/P, but these loci have not been analyzed in a Mexican Mestizo population. A complex etiology implies the presence of genetic interactions, but there is little available information regarding this in NSCL/P, and no signaling pathway has been clearly implicated in humans. Here, we analyzed the associations of 24 single nucleotide polymorphisms (SNPs) with NSCL/P in a Mexican Mestizo population (133 cases, 263 controls). The multifactorial dimensionality reduction method was used to examine gene-gene and gene-folic acid consumption interactions for the 24 SNPs analyzed in this study and for 2 additional SNPs that had previously been genotyped in the same study population. Six SNPs located in paired box 7, ventral anterior homeobox 1, sprouty RTK signaling antagonist 2, bone morphogenetic protein 4, and tropomyosin 1 genes were associated with higher risks of NSCL/P (P = 0.0001 to 0.04); 2 SNPs, 1 each in netrin 1 and V-maf avian musculoaponeurotic fibrosarcoma oncogene homolog B, were associated with a lower risk of NSCL/P (P = 0.013 to 0.03); and 2 SNPs, 1 each in ATP binding cassette subfamily A member 4 (ABCA4) and noggin, showed associations with NSCL/P that approached the threshold of significance (P = 0.056 to 0.07). In addition, 6 gene-gene interactions (P = 0.0001 to 0.001) and an ABCA4-folic acid consumption interaction (P < 0.0001) were identified. On the basis of these results, combined with those of previous association studies in the literature and biological characterizations of murine models, we propose an interaction network in which interferon regulatory factor 6 plays a central role in the etiology of NSCL/P.

Spectrum of Dental Phenotypes in Nonsyndromic Orofacial Clefting

Children with oral clefts show a wide range of dental anomalies, adding complexity to understanding the phenotypic spectrum of orofacial clefting. The evidence is mixed, however, on whether the prevalence of dental anomalies is elevated in unaffected relatives and is mostly based on small samples. In the largest international cohort to date of children with nonsyndromic clefts, their relatives, and controls, this study characterizes the spectrum of cleft-related dental anomalies and evaluates whether families with clefting have a significantly higher risk for such anomalies compared with the general population. A total of 3,811 individuals were included: 660 cases with clefts, 1,922 unaffected relatives, and 1,229 controls. Dental anomalies were identified from in-person dental exams or intraoral photographs, and case-control differences were tested using χ(2) statistics. Cases had higher rates of dental anomalies in the maxillary arch than did controls for primary (21% vs. 4%, P = 3 × 10(-8)) and permanent dentitions (51% vs. 8%, P = 4 × 10(-62)) but not in the mandible. Dental anomalies were more prevalent in cleft lip with cleft palate than other cleft types. More anomalies were seen in the ipsilateral side of the cleft. Agenesis and tooth displacements were the most common dental anomalies found in case probands for primary and permanent dentitions. Compared with controls, unaffected siblings (10% vs. 2%, P = 0.003) and parents (13% vs. 7%, P = 0.001) showed a trend for increased anomalies of the maxillary permanent dentition. Yet, these differences were nonsignificant after multiple-testing correction, suggesting genetic heterogeneity in some families carrying susceptibility to both overt clefts and dental anomalies. Collectively, the findings suggest that most affected families do not have higher genetic risk for dental anomalies than the general population and that the higher prevalence of anomalies in cases is primarily a physical consequence of the cleft and surgical interventions.

Parental risk factors for oral clefts among Central Africans, Southeast Asians, and Central Americans

Background

Several lifestyle and environmental exposures have been suspected as risk factors for oral clefts, although few have been convincingly demonstrated. Studies across global diverse populations could offer additional insight given varying types and levels of exposures.

Methods

We performed an international case–control study in the Democratic Republic of the Congo (133 cases, 301 controls), Vietnam (75 cases, 158 controls), the Philippines (102 cases, 152 controls), and Honduras (120 cases, 143 controls). Mothers were recruited from hospitals and their exposures were collected from interviewer‐administered questionnaires. We used logistic regression modeling to estimate odds ratios (OR) and 95% confidence intervals (CI).

Results

Family history of clefts was strongly associated with increased risk (maternal: OR = 4.7; 95% CI, 3.0–7.2; paternal: OR = 10.5; 95% CI, 5.9–18.8; siblings: OR = 5.3; 95% CI, 1.4–19.9). Advanced maternal age (5 year OR = 1.2; 95% CI, 1.0–1.3), pregestational hypertension (OR = 2.6; 95% CI, 1.3–5.1), and gestational seizures (OR = 2.9; 95% CI, 1.1–7.4) were statistically significant risk factors. Lower maternal (secondary school OR = 1.6; 95% CI, 1.2–2.2; primary school OR = 2.4, 95% CI, 1.6–2.8) and paternal education (OR = 1.9; 95% CI, 1.4–2.5; and OR = 1.8; 95% CI, 1.1–2.9, respectively) and paternal tobacco smoking (OR = 1.5, 95% CI, 1.1–1.9) were associated with an increased risk. No other significant associations between maternal and paternal factors were found; some environmental factors including rural residency, indoor cooking with wood, chemicals and water source appeared to be associated with an increased risk in adjusted models.

Conclusion

Our study represents one of the first international studies investigating risk factors for clefts among multiethnic underserved populations. Our findings suggest a multifactorial etiology including both maternal and paternal factors. Birth Defects Research (Part A) 103:863–879, 2015. © 2015 The Authors Birth Defects Research Part A: Clinical and Molecular Teratology Published by Wiley Periodicals, Inc.

IRF6 Sequencing in Interrupted Clefting

In a retrospective review of patients seen at the University of Alabama at Birmingham Cleft and Craniofacial Center, four patients with rare interrupted clefting were identified who had undergone genetic testing. Each of these patients had a typical cleft lip, with intact hard palate and cleft of the soft palate. Given this picture of mixed clefting, IRF6 sequencing was done and was negative for mutations in all four patients. As genetic testing for single-gene mutations and exome sequencing become clinically available, it may be possible to identify novel mutations responsible for this previously unreported type of interrupted clefting.

Accelerated wound closure in vitro by fibroblasts from a subgroup of cleft lip/palate patients: role of transforming growth factor-α

In a fraction of patients surgically treated for cleft lip/palate, excessive scarring disturbs maxillary growth and dento-alveolar development. Since certain genes are involved in craniofacial morphogenesis as well as tissue repair, a primary defect causing cleft lip/palate could lead to altered wound healing. We performed in vitro wound healing assays with primary lip fibroblasts from 16 cleft lip/palate patients. Nine foreskin fibroblast strains were included for comparison. Cells were grown to confluency and scratch wounds were applied; wound closure was monitored morphometrically over time. Wound closure rate showed highly significant differences between fibroblast strains. Statistically, fibroblast strains from the 25 individuals could be divided into three migratory groups, namely “fast”, “intermediate”, and “slow”. Most cleft lip/palate fibroblasts were distributed between the “fast” (5 strains) and the “intermediate” group (10 strains). These phenotypes were stable over different cell passages from the same individual. Expression of genes involved in cleft lip/palate and wound repair was determined by quantitative PCR. Transforming growth factor-α mRNA was significantly up-regulated in the “fast” group. 5 ng/ml transforming growth factor-α added to the culture medium increased the wound closure rate of cleft lip/palate strains from the “intermediate” migratory group to the level of the “fast”, but had no effect on the latter group. Conversely, antibody to transforming growth factor-α or a specific inhibitor of its receptor most effectively reduced the wound closure rate of “fast” cleft lip/palate strains. Thus, fibroblasts from a distinct subgroup of cleft lip/palate patients exhibit an increased migration rate into wounds in vitro, which is linked to higher transforming growth factor-α expression and attenuated by interfering with its signaling.

Association of Transforming Growth Factor Alpha and Methylenetetrahydrofolate reductase gene variants with nonsyndromic cleft lip and palate in the Indian population

Objectives:

The aim was to evaluate the relationship of the K-primer variant of the transforming growth factor-alpha (TGF-α) gene and C677T variant of the methylenetetrahydrofolate reductase (MTHFR) gene with nonsyndromic cleft lip and palate (CL/P) in the Indian population.

Setting and Sample Population:

The study group consisted of DNA samples of 25 subjects with nonsyndromic CL with or without cleft palate and 25 unrelated controls, already existing in the Department of Orthodontics, D.A.P.M.R.V. Dental College, Bengaluru, Karnataka, India.

Materials and Methods:

The DNA samples were divided into two categories: Group A which included the 25 subjects with nonsyndromic CL/P; and Group B, which consisted of the 25 unrelated controls. The polymerase chain reaction (PCR) test was done for amplification of the region of interest from the DNA samples. Restriction digestion was then performed on the amplified product using the restriction enzyme HinfI, separately for each of the variants. The digested PCR products were separated into channels on a 1.5% agarose gel containing ethidium bromide in an electrophoretic chamber. A U.V. transilluminator was used to see the specific bands of base pairs of the digested PCR products.

Results:

In Group A, the TGF-α gene variant was present in 16 subjects (P = 0.001) and MTHFR gene variant was present in 8 subjects (P = 0.185). A combination of both gene variants were present in seven subjects, which was an interesting finding. In Group B, four subjects tested positive for the TGF-α and MTHFR gene variants.

Conclusions:

The TGF-α gene variant and a combination of TGF-α + MTHFR gene variants significantly contribute to the development of nonsyndromic CL/P and can be considered as genetic markers for Indian population. The MTHFR gene variant, though a minor risk factor, cannot be considered as a genetic marker.

Molecular and bioenergetic differences between cells with African versus European inherited mitochondrial DNA haplogroups: implications for population susceptibility to diseases

The geographic origins of populations can be identified by their maternally inherited mitochondrial DNA (mtDNA) haplogroups. This study compared human cybrids (cytoplasmic hybrids), which are cell lines with identical nuclei but mitochondria from different individuals with mtDNA from either the H haplogroup or L haplogroup backgrounds. The most common European haplogroup is H while individuals of maternal African origin are of the L haplogroup. Despite lower mtDNA copy numbers, L cybrids had higher expression levels for nine mtDNA-encoded respiratory complex genes, decreased ATP (adenosine triphosphate) turnover rates and lower levels of reactive oxygen species production, parameters which are consistent with more efficient oxidative phosphorylation. Surprisingly, GeneChip arrays showed that the L and H cybrids had major differences in expression of genes of the canonical complement system (5 genes), dermatan/chondroitin sulfate biosynthesis (5 genes) and CCR3 (chemokine, CC motif, receptor 3) signaling (9 genes). Quantitative nuclear gene expression studies confirmed that L cybrids had (a) lower expression levels of complement pathway and innate immunity genes and (b) increased levels of inflammation-related signaling genes, which are critical in human diseases. Our data support the hypothesis that mtDNA haplogroups representing populations from different geographic origins may play a role in differential susceptibilities to diseases.

Contribution of transforming growth factor α polymorphisms to nonsyndromic orofacial clefts: a HuGE review and meta-analysis

We performed a meta-analysis of the association of transforming growth factor α gene (TGFA) polymorphisms with the risk of cleft lip with or without cleft palate (CL/P) or cleft palate (CP). In total, data from 29 studies were pooled for the following 3 polymorphisms: TGFA/TaqI, TGFA/BamHI, and TGFA/RasI in the TGFA gene. A fixed-effects or random-effects model was used to calculate the pooled odds ratios based on the results from the heterogeneity tests. A significantly increased CL/P or CP risk was observed in persons carrying a C2 allele at the TaqI polymorphism (odds ratio (OR) = 1.70, 95% confidence interval (CI): 1.41, 2.05) compared with those with a C1 allele (OR = 1.57, 95% CI: 1.23, 2.01). For the TGFA/BamHI polymorphism, carriers of the minor A1 allele had an estimated relative decrease in CL/P risk (OR = 0.44, 95% CI: 0.30, 0.64). These associations remained significant when only high-quality studies were included. However, no significant association was observed between the TGFA/RasI variant and CL/P risk. In summary, this meta-analysis provided a robust estimate of the positive association of the TGFA/TaqI polymorphism with both CL/P and CP and suggests that persons with an A1 allele may have a markedly decreased risk of CL/P.

Search for genetic modifiers of IRF6 and genotype-phenotype correlations in Van der Woude and popliteal pterygium syndromes

Van der Woude syndrome is the most common form of syndromic orofacial clefting, accounting for 1-2% of all patients with cleft lip and/or cleft palate. Van der Woude and popliteal pterygium syndromes are caused by mutations in IRF6, but phenotypic variability within and among families with either syndrome suggests that other genetic factors contribute to the phenotypes. The aim of this study was to identify common variants acting as genetic modifiers of IRF6 as well as genotype-phenotype correlations based on mutation type and location. We identified an association between mutations in the DNA-binding domain of IRF6 and limb defects (including pterygia). Although we did not detect formally significant associations with the genes tested, borderline associations suggest several genes that could modify the VWS phenotype, including FOXE1, TGFB3, and TFAP2A. Some of these genes are hypothesized to be part of the IRF6 gene regulatory network and may suggest additional genes for future study when larger sample sizes are also available. We also show that families with the Van de Woude phenotype but in whom no mutations have been identified have a lower frequency of cleft lip, suggesting there may be locus and/or mutation class differences in Van de Woude syndrome.

TGFA and IRF6 contribute to the risk of nonsyndromic cleft lip with or without cleft palate in northeast China

Nonsyndromic cleft lip with or without cleft palate (NSCL/P) are common birth defects with a complex etiology. Multiple interacting loci and possible environmental factors influence the risk of NSCL/P. 12 single nucleotide polymorphisms (SNPs) in 7 candidate genes were tested using an allele-specific primer extension for case-control and case-parent analyses in northeast China in 236 unrelated patients, 185 mothers and 154 fathers, including 128 complete trios, and 400 control individuals. TGFA and IRF6 genes showed a significant association with NSCL/P. In IRF6, statistical evidence of an association between rs2235371 (p = 0.003), rs2013162 (p<0.0001) and NSCL/P was observed in case-control analyses. Family based association tests (FBATs) showed over-transmission of the C allele at the rs2235371 polymorphism (p = 0.007). In TGFA, associations between rs3771494, rs3771523 (G3822A), rs11466285 (T3851C) and NSCL/P were observed in case-control and FBAT analyses. Associations between other genes (BCL3, TGFB3, MTHFR, PVRL1 and SUMO1) and NSCL/P were not detected.

Smad4-Irf6 genetic interaction and TGFβ-mediated IRF6 signaling cascade are crucial for palatal fusion in mice

Cleft palate is one of the most common human birth defects and is associated with multiple genetic and environmental risk factors. Although mutations in the genes encoding transforming growth factor beta (TGFβ) signaling molecules and interferon regulatory factor 6 (Irf6) have been identified as genetic risk factors for cleft palate, little is known about the relationship between TGFβ signaling and IRF6 activity during palate formation. Here, we show that TGFβ signaling regulates expression of Irf6 and the fate of the medial edge epithelium (MEE) during palatal fusion in mice. Haploinsufficiency of Irf6 in mice with basal epithelial-specific deletion of the TGFβ signaling mediator Smad4 (Smad4(fl/fl);K14-Cre;Irf6(+/R84C)) results in compromised p21 expression and MEE persistence, similar to observations in Tgfbr2(fl/fl);K14-Cre mice, although the secondary palate of Irf6(+/R84C) and Smad4(fl/fl);K14-Cre mice form normally. Furthermore, Smad4(fl/fl);K14-Cre;Irf6(+/R84C) mice show extra digits that are consistent with abnormal toe and nail phenotypes in individuals with Van der Woude and popliteal pterygium syndromes, suggesting that the TGFβ/SMAD4/IRF6 signaling cascade might be a well-conserved mechanism in regulating multiple organogenesis. Strikingly, overexpression of Irf6 rescued p21 expression and MEE degeneration in Tgfbr2(fl/fl);K14-Cre mice. Thus, IRF6 and SMAD4 synergistically regulate the fate of the MEE, and TGFβ-mediated Irf6 activity is responsible for MEE degeneration during palatal fusion in mice.