Genome-wide association study identifies a new susceptibility locus for cleft lip with or without a cleft palate

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

AIMS

To investigate the association between single nucleotide polymorphisms (SNPs) in 3'UTR region of VAX1, SYT14 and PAX7 genes and the risk of non-syndromic cleft palate (NSCLP) in a northwest Chinese population.

MAIN METHODS

A case-control study was conducted in 406 normal controls and 399 NSCLP patients. Using iMLDRTM genotyping technology, eight SNPs of three genes ((rs10787760, rs7086344 at VAX1), (rs1010113, rs851114, and rs485874 at PAX7), and (rs61820397, rs4609425, rs12133399 at SYT14)) were genotyped to investigate the differences in alleles and genotype distribution frequencies between NSCLP patients and healthy controls. RNA Folding Form software was used to predict RNA secondary structure and expression vectors were constructed to explore the function of the relevant SNP. The effect of SNP polymorphism of gene transcription and translation was assessed using qPCR and Western blot analysis.

KEY FINDINGS

Among the eight SNPs of three genes, rs10787760 of VAX1 gene was found to be associated with an increased risk of NSCLP (OR = 1.341, CI = 1.004-1.790) and the GA genotype of rs10787760 increased the risk of cleft lip and/or palate (CL/P) about 1.42 times (p < 0.05), and carrying the A allele might increase the risk of NSCL/P in male (OR = 1.356, 95 % CI = 1.010-1.823). But there was no association observed with cleft palate only (CPO). Cell function experiments revealed that the G to A mutation in rs10787760 up-regulated GFP-VAX1 transcriptional level by 2.39 and 3.13 times in two cell lines respectively, and enhance the protein expression of the VAX1 gene further. RNA secondary structure study showed that the rs10787760 (G > A) had two different secondary structures in 3'UTR region.

SIGNIFICANCE

The rs10787760 variant in the 3'UTR region of VAX1 gene is associated with CL/P in northwest Chinese population. We hypothesize that the machanism of it might be caused by the RNA differenct fold in the 3'UTR region caused by the polymorphism of the gene.

LEVEL OF EVIDENCE

Original Reports.

DeepFace: Deep-learning-based framework to contextualize orofacial-cleft-related variants during human embryonic craniofacial development

Orofacial clefts (OFCs) are among the most common human congenital birth defects. Previous multiethnic studies have identified dozens of associated loci for both cleft lip with or without cleft palate (CL/P) and cleft palate alone (CP). Although several nearby genes have been highlighted, the "casual" variants are largely unknown. Here, we developed DeepFace, a convolutional neural network model, to assess the functional impact of variants by SNP activity difference (SAD) scores. The DeepFace model is trained with 204 epigenomic assays from crucial human embryonic craniofacial developmental stages of post-conception week (pcw) 4 to pcw 10. The Pearson correlation coefficient between the predicted and actual values for 12 epigenetic features achieved a median range of 0.50-0.83. Specifically, our model revealed that SNPs significantly associated with OFCs tended to exhibit higher SAD scores across various variant categories compared to less related groups, indicating a context-specific impact of OFC-related SNPs. Notably, we identified six SNPs with a significant linear relationship to SAD scores throughout developmental progression, suggesting that these SNPs could play a temporal regulatory role. Furthermore, our cell-type specificity analysis pinpointed the trophoblast cell as having the highest enrichment of risk signals associated with OFCs. Overall, DeepFace can harness distal regulatory signals from extensive epigenomic assays, offering new perspectives for prioritizing OFC variants using contextualized functional genomic features. We expect DeepFace to be instrumental in accessing and predicting the regulatory roles of variants associated with OFCs, and the model can be extended to study other complex diseases or traits.

High incidence and geographic distribution of cleft palate cases in Finland are associated with a regulatory variant in IRF6

In Finland the frequency of isolated cleft palate (CP) is higher than that of isolated cleft lip with or without cleft palate (CL/P). This trend contrasts to that in other European countries but its genetic underpinnings are unknown. We performed a genome-wide association study for orofacial clefts, which include CL/P and CP, in the Finnish population. We identified rs570516915, a single nucleotide polymorphism that is highly enriched in Finns and Estonians, as being strongly associated with CP ( P = 5.25 × 10 -34 , OR = 8.65, 95% CI 6.11-12.25), but not with CL/P ( P = 7.2 × 10 -5 ), with genome-wide significance. The risk allele frequency of rs570516915 parallels the regional variation of CP prevalence in Finland, and the association was replicated in independent cohorts of CP cases from Finland ( P = 8.82 × 10 -28 ) and Estonia ( P = 1.25 × 10 -5 ). The risk allele of rs570516915 disrupts a conserved binding site for the transcription factor IRF6 within a previously characterized enhancer upstream of the IRF6 gene. Through reporter assay experiments we found that the risk allele of rs570516915 diminishes the enhancer activity. Oral epithelial cells derived from CRISPR-Cas9 edited induced pluripotent stem cells demonstrate that the CP-associated allele of rs570516915 concomitantly decreases the binding of IRF6 and the expression level of IRF6 , suggesting impaired IRF6 autoregulation as a molecular mechanism underlying the risk for CP.

Shared genetic risk between major orofacial cleft phenotypes in an African population

Nonsyndromic orofacial clefts (NSOFCs) represent a large proportion (70%-80%) of all OFCs. They can be broadly categorized into nonsyndromic cleft lip with or without cleft palate (NSCL/P) and nonsyndromic cleft palate only (NSCPO). Although NSCL/P and NSCPO are considered etiologically distinct, recent evidence suggests the presence of shared genetic risks. Thus, we investigated the genetic overlap between NSCL/P and NSCPO using African genome-wide association study (GWAS) data on NSOFCs. These data consist of 814 NSCL/P, 205 NSCPO cases, and 2159 unrelated controls. We generated common single-nucleotide variants (SNVs) association summary statistics separately for each phenotype (NSCL/P and NSCPO) under an additive genetic model. Subsequently, we employed the pleiotropic analysis under the composite null (PLACO) method to test for genetic overlap. Our analysis identified two loci with genome-wide significance (rs181737795 [p = 2.58E-08] and rs2221169 [p = 4.5E-08]) and one locus with marginal significance (rs187523265 [p = 5.22E-08]). Using mouse transcriptomics data and information from genetic phenotype databases, we identified MDN1, MAP3k7, KMT2A, ARCN1, and VADC2 as top candidate genes for the associated SNVs. These findings enhance our understanding of genetic variants associated with NSOFCs and identify potential candidate genes for further exploration.

Genetic association and functional validation of ZFP36L2 in non-syndromic orofacial cleft subtypes

BACKGROUND

Non-syndromic orofacial cleft (NSOC) is one of the most common craniofacial malformations with complex etiology. This study aimed to explore the role of specific SNPs in ZFP36L2 and its functional relevance in zebrafish models.

METHODS

We analyzed genetic data of the Chinese Han population from two previous GWAS, comprising of 2512 cases and 2255 controls. Based on the Hardy-Weinberg Equilibrium (HWE) and minor allele frequency (MAF), SNPs in the ZFP36L2 were selected for association analysis. In addition, zebrafish models were used to clarify the in-situ expression pattern of zfp36l2 and the impact of its Morpholino-induced knockdown.

RESULTS

Via association analysis, rs7933 in ZFP36L2 was significantly associated with various non-syndromic cleft lip-only subtypes, potentially conferring a protective effect. Zebrafish embryos showed elevated expression of zfp36l2 in the craniofacial region during critical stages of oral cavity formation. Furthermore, Morpholino-induced knockdown of zfp36l2 led to craniofacial abnormalities, including cleft lip, which was partially rescued by the addition of zfp36l2 mRNA.

CONCLUSION

Our findings highlight the significance of ZFP36L2 in the etiology of NSOC, supported by both human genetic association data and functional studies in zebrafish. These results pave the way for further exploration of targeted interventions for craniofacial malformations.

MMP16 as NSCL ± P Susceptible Gene in Western Han Chinese

OBJECTIVE

The role of MMP16 in lip development is unclear. This study aimed to identify nonsyndromic cleft lip with or without palate (NSCL ± P) susceptible loci of MMP16 in western Han Chinese.

DESIGN

We performed targeted sequencing around MMP16 combined with a 2-phase association analysis on common variants. Phase 2 association analysis was performed with NSCL ± P specific subphenotypes (NSCL and NSCLP). Then we used rare variants burden analysis and genotyping, accompanied by motif analysis.

SETTING

This study was completed in a tertiary medical center.

PATIENTS, PARTICIPANTS

Phase 1 targeted sequencing included 159 patients with NSCL ± P and 542 normal controls; phase 2 included 1626 patients with NSCL ± P (1047 NSCL and 579 NSCLP) and 2255 normal controls.

INTERVENTIONS

Venous blood samples were collected from patients and used to extract DNA.

MAIN OUTCOME MEASURES

After Bonferroni correction, phase 1 significant threshold of p-value was 4.28 × 10-5 (0.05/1167 single nucleotide polymorphisms [SNPs]), and phase 2 was .00025 (0.05/200 SNPs). Burden analysis significant threshold p-value was .05.

RESULTS

Common variants phase 1 association analysis identified 11 statistically significant SNPs (lowest p = 1.90 × 10-9, odds ratio (OR) = 0.27, 95% CI: 0.17-0.44), phase 2 replication identified 16 SNPs in NSCL ± P (lowest p = 6.26 × 10-6, OR = 0.77, 95% CI: 0.69-0.86) and 9 in NSCL (lowest p = 8.44 × 10-5, OR = 0.76, 95% CI: 0.66-0.87). Rare variants burden analysis showed no significant results, genotyping results showed they were maternally inherited.

CONCLUSIONS

Our study identified MMP16 susceptible SNPs in NSCL ± P and NSCL, emphasizing its potential role in lip development. Our study also highlighted the importance to perform association analysis with subphenotypes divided.

Ezh2-dependent methylation in oral epithelia promotes secondary palatogenesis

BACKGROUND

In addition to genomic risk variants and environmental influences, increasing evidence suggests epigenetic modifications are important for orofacial development and their alterations can contribute to orofacial clefts. Ezh2 encodes a core catalytic component of the Polycomb repressive complex responsible for addition of methyl marks to Histone H3 as a mechanism of repressing target genes. The role of Ezh2 in orofacial clefts remains unknown.

AIMS

To investigate the epithelial role of Ezh2-dependent methylation in secondary palatogenesis.

METHODS

We used conditional gene-targeting methods to ablate Ezh2 in the surface ectoderm-derived oral epithelium of mouse embryos. We then performed single-cell RNA sequencing combined with immunofluorescence and RT-qPCR to investigate gene expression in conditional mutant palate. We also employed double knockout analyses of Ezh1 and Ezh2 to address if they have synergistic roles in palatogenesis.

RESULTS

We found that conditional inactivation of Ezh2 in oral epithelia results in partially penetrant cleft palate. Double knockout analyses revealed that another family member Ezh1 is dispensable in orofacial development, and it does not have synergistic roles with Ezh2 in palatogenesis. Histochemistry and single-cell RNA-seq analyses revealed dysregulation of cell cycle regulators in the palatal epithelia of Ezh2 mutant mouse embryos disrupts palatogenesis.

CONCLUSION

Ezh2-dependent histone H3K27 methylation represses expression of cell cycle regulator Cdkn1a and promotes proliferation in the epithelium of the developing palatal shelves. Loss of this regulation may perturb movement of the palatal shelves, causing a delay in palate elevation which may result in failure of the secondary palate to close altogether.

Genome-wide meta-analyses identify five new risk loci for nonsyndromic orofacial clefts in the Chinese Han population

BACKGROUND

Nonsyndromic orofacial clefts (NSOFCs) are the most common craniofacial birth malformations in humans and are generally classified as nonsyndromic cleft lip with or without cleft palate (NSCL/P) and nonsyndromic cleft palate only (NSCPO). Genome-wide association studies (GWASs) of NSOFCs have demonstrated multiple risk loci and candidate genes; however, published risk factors are able to explain only a small fraction of the observed NSOFCs heritability.

METHODS

Here, we performed GWASs of 1615 NSCPO cases and 2340 controls, and then conducted genome-wide meta-analyses of NSOFCs, totaling 6812 NSCL/P cases, 2614 NSCPO cases, and 19,165 controls from the Chinese Han population.

RESULTS

We identify 47 risk loci with genome-wide pmeta -value <5.0 × 10-8 , 5 risk loci (1p32.1, 3p14.1, 3p14.3, 3p21.31, and 13q22.1) of which are new. All of the 47 susceptibility loci conjointly account for 44.12% of the NSOFCs' heritability in the Chinese Han population.

CONCLUSION

Our results improve the comprehending of genetic susceptibility to NSOFCs and provide new views into the genetic etiology of craniofacial anomalies.

Translational Research Techniques for the Facial Plastic Surgeon: An Overview

The field of facial plastic and reconstructive surgery (FPRS) is an incredibly diverse, multispecialty field that seeks innovative and novel solutions for the management of physical defects on the head and neck. To aid in the advancement of medical and surgical treatments for these defects, there has been a recent emphasis on the importance of translational research. With recent technological advancements, there are now a myriad of research techniques that are widely accessible for physician and scientist use in translational research. Such techniques include integrated multiomics, advanced cell culture and microfluidic tissue models, established animal models, and emerging computer models generated using bioinformatics. This study discusses these various research techniques and how they have and can be used for research in the context of various important diseases within the field of FPRS.

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

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

Identification of a Novel Variant of PDGFC Associated with Nonsyndromic Cleft Lip and Palate in a Chinese Family

Nonsyndromic cleft lip with or without cleft palate (NSCL/P) accounts for 70% of the total number of patients with cleft lip with or without cleft palate (CL/P) and is the most common type of congenital deformity of the craniomaxillofacial region. In this study, whole exome sequencing (WES) and Sanger sequencing were performed on affected members of a Han Chinese family, and a missense variant in the platelet-derived growth factor C (PDGFC) gene (NM_016205: c.G93T: p.Q31H) was identified to be associated with NSCL/P. Bioinformatic studies demonstrated that the amino acid corresponding to this variation is highly conserved in many mammals and leads to a glutamine-to-histidine substitution in an evolutionarily conserved DNA-binding domain. It was found that the expression of PDGFC was significantly decreased in the dental pulp stem cells (DPSCs) of NSCL/P cases, compared to the controls, and that the variant (NM_016205: c.G93T) reduced the expression of PDGFC. In addition, the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed that Pdgfc deficiency disrupted NSCL/P-related signaling pathways such as the MAPK signaling pathway and cell adhesion molecules. In conclusion, our study identified a missense variant (NM_016205: c.G93T) in exon 1 of PDGFC potentially associated with susceptibility to NSCL/P.

A Variant in the IRF6 Promoter Associated with the Risk for Orofacial Clefting

The single-nucleotide polymorphism (SNP) rs2235371 (IRF6 V274I) is associated with nonsyndromic cleft lip with or without cleft palate (NSCL/P) in Han Chinese and other populations but appears to be without a functional effect. To find the common etiologic variant or variants within the haplotype tagged by rs2235371, we carried out targeted sequencing of an interval containing IRF6 in 159 Han Chinese with NSCL/P. This study revealed that the SNP rs12403599, within the IRF6 promoter, is associated with all phenotypes of NSCL/P, especially nonsyndromic cleft lip (NSCLO) and a subphenotype of it, microform cleft lip (MCL). This association was replicated in 2 additional much larger cohorts of cases and controls from the Han Chinese. Conditional logistic analysis indicated that association of rs2235371 with NSCL/P was lost if rs12403599 was excluded. rs12403599 contributes the most risk to MCL: its G allele is responsible for 38.47% of the genetic contribution to MCL, and the odds ratios of G/C and G/G genotypes were 2.91 and 6.58, respectively, for MCL. To test if rs12403599 is functional, we carried out reporter assays in a fetal oral epithelium cells (GMSM-K). Unexpectedly, the risk allele G yielded higher promoter activity in GMSM-K. Consistent with the reporter studies, expression of IRF6 in lip tissues from NSCLO and MCL patients with the G/G phenotype was higher than in those from patients with the C/C phenotype. These results indicate that rs12403599 is tagging the risk haplotype for NSCL/P better than rs2235371 in Han Chinese and supports investigation of the mechanisms by which the allele of rs12403599 affects IRF6 expression and tests of this association in different populations.

Integrating transcriptomics and genomics to identify fibroblast growth factor/receptor candidate genes for non-syndromic orofacial cleft in Chinese

OBJECTIVES

The objective of this study was to explore the relationship between fibroblast growth factor/receptor (FGF/FGFR) and non-syndromic orofacial cleft (NSOC) in individuals of Han Chinese.

DESIGN

Initially, we performed RNA-Seq between non-syndromic cleft lip only (NSCLO) or non-syndromic cleft palate only (NSCPO) and control groups. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were carried out to evaluate the functions of differentially expressed genes (DEGs) of FGF/FGFR. Then, we selected the most significant DEG FGFR2 and performed an association analysis in Chinese. Linkage disequilibrium (LD) and haplotype analyses were performed with HaploView and PLINK. Additional bioinformatics functional prediction for the notable single nucleotide polymorphisms was performed with HaploReg V4.1 and 3DSNP.

RESULTS

Finally, we identified 32 mRNAs related to FGF/FGFR via RNA-Seq and chose FGFR2 in the subsequent association analysis. Results indicated that the single nucleotide polymorphism (SNP) rs2288336 in FGFR2 contributed significantly to both non-syndromic cleft lip with or without cleft palate (NSCL/P) and NSCLO, with p values of 5.00E-05 (OR = 0.79, 95% CI: 0.70-0.88) and 1.38E-04 (OR = 0.76, 95% CI: 0.65-0.87), respectively. In addition, rs3793893 in FGFR2 was found to be associated with NSCLO, with a p value of 1.02E-04 (OR = 0.67, 95% CI: 0.55-0.82).

CONCLUSIONS

Our research demonstrated that FGFR2 is significantly more involved in NSOC than other FGF/FGFRs in Chinese and further identified rs2288336 and rs3793893 in FGFR2 associated with NSOC subtypes, which provide further evidence for the genetic etiology of NSOC in Han Chinese.

Genetic Variants in Mammalian STE20-like Protein Kinase 2 were associated with risk of NSCL/P

AIM

Mammalian STE20-like protein kinase 2 (MST2) plays an important role in apoptosis and the development of many disorders. Here, we aim to explore if genetic variants in MST2 are associated with the risk of non-syndromic cleft lip with or without palate (NSCL/P).

MATERIALS AND METHODS

The association study was performed in a two-stage study of 1,069 cases and 1,724 controls to evaluate the association between genetic variants in the MST2 and NSCL/P risk. The potential function of the candidate single nucleotide polymorphism (SNP) was predicted using HaploReg, RegulomeDB, and public craniofacial histone chromatin immunoprecipitation sequencing (ChIP-seq) data. Haploview was used to perform the haplotype of risk alleles. The expression quantitative trait loci (eQTL) effect was assessed using the Genotype-Tissue Expression (GTEx) project. Gene expression in mouse embryo tissue was performed using data downloaded from GSE67985. The potential role of candidate gene in the development of NSCL/P was assessed by correlation and enrichment analysis.

RESULTS

Among SNPs in MST2, rs2922070 C allele (P_meta = 2.93E-04) and rs6988087 T allele (P_meta = 1.57E-03) were linked with significantly increased risk of NSCL/P. Rs2922070, rs6988087 and their high linkage disequilibrium (LD) SNPs constituted a risk haplotype of NSCL/P. Individuals carrying 3-4 risk alleles had an elevated risk of NSCL/P compared to those who carried less risk alleles (P = 2.00E-04). The eQTL analysis revealed a significant association between these two variants and MST2 in muscle tissue of the body. The MST2 expressed during mouse craniofacial development and over-expressed in the human orbicularis oris muscle (OOM) of NSCL/P patients compared to controls. MST2 was involved in the development of NSCL/P by regulating the mRNA surveillance pathway, the MAPK signaling pathway, the neurotrophin signaling pathway, the FoxO signaling pathway and the VEGF signaling pathway.

CONCLUSION

MST2 was associated with the development of NSCL/P.

Association between variants around IRF6 and non-syndromic orofacial cleft in Western Han Chinese

OBJECTIVES

Considering limitations of previous studies and differences across populations and subtypes, this study aimed to identify new potential SNPs around IRF6 associated with non-syndromic orofacial cleft (NSOC) in Western Han Chinese.

MATERIALS AND METHODS

We recruited 376 NSOC case-parent trios, including 125 non-syndromic cleft lip only (NSCLO) trios, 151 non-syndromic cleft lip and palate (NSCLP) trios, and 100 non-syndromic cleft palate only (NSCPO) trios. Twenty-two single-nucleotide polymorphisms (SNPs) were genotyped using MassARRAY method. Hardy-Weinberg equilibrium test, allelic transmission disequilibrium test (TDT) analysis, sliding-window haplotype TDT analysis, and tests for parent-of-origin effect were performed using the PLINK software. Pairwise linkage disequilibrium (LD) was computed using the Haploview program.

RESULTS

In TDT analysis, allele A at rs17015217 (p = 0.00011, OR = 0.61 and 95% CI: 0.47-0.78) and allele T at rs12080691 (p = 0.00011, OR = 0.61 and 95% CI: 0.47-0.78) were under-transmitted among NSCLO trios but over-transmitted among NSCPO trios. Haplotypes showing evidence of under-transmission in NSCLO trios were over-transmitted in NSCPO trios. In tests for parent-of-origin effects, T allele at rs12080691 presented paternal under-transmission among NSCLO trios but over-transmission among NSCPO trios.

CONCLUSIONS

Allele A at rs17015217 and allele T at rs12080691 are associated with NSCLO and NSCPO with potential to have opposite effects on two subtypes in this sample from Western Han Chinese.

Genome Analysis Using Whole-Exome Sequencing of Non-Syndromic Cleft Lip and/or Palate from Malagasy Trios Identifies Variants Associated with Cilium-Related Pathways and Asian Genetic Ancestry

Background: Orofacial clefts (OFCs) are common congenital disabilities that can occur as isolated non-syndromic events or as part of Mendelian syndromes. OFC risk factors vary due to differences in regional environmental exposures, genetic variants, and ethnicities. In recent years, significant progress has been made in understanding OFCs, due to advances in sequencing and genotyping technologies. Despite these advances, very little is known about the genetic interplay in the Malagasy population. Methods: Here, we performed high-resolution whole-exome sequencing (WES) on non-syndromic cleft lip with or without palate (nCL/P) trios in the Malagasy population (78 individuals from 26 families (trios)). To integrate the impact of genetic ancestry admixture, we computed both global and local ancestries. Results: Participants demonstrated a high percentage of both African and Asian admixture. We identified damaging variants in primary cilium-mediated pathway genes WNT5B (one family), GPC4 (one family), co-occurrence in MSX1 (five families), WDR11 (one family), and tubulin stabilizer SEPTIN9 (one family). Furthermore, we identified an autosomal homozygous damaging variant in PHGDH (one family) gene that may impact metabiotic activity. Lastly, all variants were predicted to reside on local Asian genetic ancestry admixed alleles. Conclusion: Our results from examining the Malagasy genome provide limited support for the hypothesis that germline variants in primary cilia may be risk factors for nCL/P, and outline the importance of integrating local ancestry components better to understand the multi-ethnic impact on nCL/P.

Identification of putative regulatory single-nucleotide variants in NTN1 gene associated with NSCL/P

Non-syndromic cleft lip with or without cleft palate (NSCL/P) is a common polygenetic disease. Although genome-wide association studies (GWAS) identified NTN1 gene as a high-priority candidate of NSCL/P, the comprehensive genetic architecture of NTN1 weren't yet known. Thus, this study aimed to determine full-scale genetic variants of NTN1 for NSCL/P in Chinese Han people. Initially, targeted sequencing of NTN1 gene was performed on 159 NSCL/P patients to identify susceptible single nucleotide polymorphisms (SNPs) associated with NSCL/P. Then, association analysis and burden analysis were separately used to validate the common variants and rare variants identified among large size of samples (1608 NSCL/P cases and 2255 controls). Additionally, NSCL/P subtype association analysis was applied to elucidate the etiology discrepancy of non-syndromic cleft lip with palate (NSCLP) and non-syndromic cleft lip only (NSCLO). Lastly, bioinformatics analysis was performed to annotate and prioritize candidate variants. We found 15 NSCL/P-associated SNPs including rs4791774 (P = 1.10E-08, OR = 1.467, 95% CI: 1.286~1.673) and rs9788972 (P = 1.28E-07, OR = 1.398, 95% CI : 1.235~1.584) originally detected by previous GWASs in Chinese Han ancestry. Four NSCLO risk-associated SNPs and eight specific NSCLP associated SNPs were found. Three SNPs (rs4791331, rs4791774 and rs9900753) were predicted to locate at regulatory region of NTN1. Our study validated the association between NTN1 gene and pathogenesis of NSCL/P and reinforced the hypothesis that NSCLP have a different etiology from NSCLO. We also identified three putative regulatory SNPs in NTN1 gene.

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.

A Comprehensive Investigation on Potential Risk Factors for NSCL/P in a Rural District of Hebei Province, China

Non syndromic cleft lip with or without palate (NSCL/P), one of the most common birth defects, is closely related to various risk factors. However, information regarding risk factors for NSCL/P in rural districts in China is very limited thus far. The objective of this study was designed to identify the potential risk factors for NSCL/P in rural districts.

A comprehensive retrospective investigation including 435 NSCL/P patients and 402 healthy children was carried out in Hebei Province, China. Multiple logistic regression analysis and transmission disequilibrium test (TDT) were respectively used to identify non-genetic and genetic risk factors for NSCL/P, and then PLINK was used to explore the relationship between non-genetic and genetic risk factors.

The results showed that maternal periconceptional exposure to pesticides and herbicides, as well as low parental education level were involved in the increased risk of NSCL/P, whereas maternal folic acid and multivitamin supplementation use during preconception period were associated with the reduced risk of NSCL/P. TDT analysis identified 2 single nucleotide polymorphisms (SNPs) (rs7078160 and rs4752028) in VAX1 and one SNP (rs17563) in BMP4 as the genetic risk factors for NSCL/P. Further analysis showed that the genetic risk factors were closely related with the negative non-genetic risk factors.

Our study identified the potential risk factors for NSCL/P in rural districts, thus providing a theoretical basis for the prevention of NSCL/P occurrence.

Gene×environment associations in orofacial clefting

This chapter reviews the evidence of gene×environment interactions (G×E) in the etiology of orofacial cleft birth defects (OFCs), specifically cleft lip (CL), cleft palate (CP), and cleft lip with or without cleft palate (CL/P). We summarize the current state of our understanding of the genetic architecture of nonsyndromic OFCs and the evidence that maternal exposures during pregnancy influence risk of OFCs. Further, we present possible candidate gene pathways for these exposures including metabolism of folates, metabolism of retinoids, retinoic acid receptor signaling, aryl hydrocarbon receptor signaling, glucocorticoid receptor signaling, and biotransformation and transport. We review genes in these pathways with prior evidence of association with OFCs, genes with evidence from prior candidate gene G×E studies, and genes identified from genome-wide searches specifically for identifying G×E. Finally, we suggest future directions for G×E research in OFCs.

Chromatin conformation of human oral epithelium can identify orofacial cleft missing functional variants

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

Targeted re-sequencing on 1p22 among non-syndromic orofacial clefts from Han Chinese population

Rs560426 at 1p22 was proved to be associated with NSCL/P (non-syndromic cleft lip with or without the palate) in several populations, including Han Chinese population. Here, we conducted a deep sequencing around rs560426 to locate more susceptibility variants in this region. In total, 2,293 NSCL/P cases and 3,235 normal controls were recruited. After sequencing, association analysis was performed. Western blot, RT-qPCR, HE, immunofluorescence staining, and RNA sequencing were conducted for functional analyses of the selected variants. Association analysis indicated that rs77179923 was the only SNP associated with NSCLP specifically (p = 4.70E-04, OR = 1.84), and rs12071152 was uniquely associated with LCLO (p = 4.00E-04, OR = 1.30, 95%CI: 1.12–1.51). Moreover, de novo harmful rare variant NM_004815.3, NP_004806.3; c.1652G>C, p.R551T in ARHGAP29 resulted in a decreased expression level of ARHGAP29, which in turn affected NSCL/P-related biological processes; however, no overt cleft palate (CP) phenotype was observed. In conclusion, rs12071152 was a new susceptible variant, which is specifically associated with LCLO among the Han Chinese population. Allele A of it could increase the risk of having a cleft baby. Rs77179923 and rare variant NM_004815.3, NP_004806.3; c.1652G>C, p.R551T at 1p22 were both associated with NSCLP among the Han Chinese population. However, this missense variation contributes to no overt CP phenotype due to dosage insufficiency or compensation from other genes.

Whole-genome sequencing reveals de-novo mutations associated with nonsyndromic cleft lip/palate

The majority (85%) of nonsyndromic cleft lip with or without cleft palate (nsCL/P) cases occur sporadically, suggesting a role for de novo mutations (DNMs) in the etiology of nsCL/P. To identify high impact protein-altering DNMs that contribute to the risk of nsCL/P, we conducted whole-genome sequencing (WGS) analyses in 130 African case-parent trios (affected probands and unaffected parents). We identified 162 high confidence protein-altering DNMs some of which are based on available evidence, contribute to the risk of nsCL/P. These include novel protein-truncating DNMs in the ACTL6A, ARHGAP10, MINK1, TMEM5 and TTN genes; as well as missense variants in ACAN, DHRS3, DLX6, EPHB2, FKBP10, KMT2D, RECQL4, SEMA3C, SEMA4D, SHH, TP63, and TULP4. Many of these protein-altering DNMs were predicted to be pathogenic. Analysis using mouse transcriptomics data showed that some of these genes are expressed during the development of primary and secondary palate. Gene-set enrichment analysis of the protein-altering DNMs identified palatal development and neural crest migration among the few processes that were significantly enriched. These processes are directly involved in the etiopathogenesis of clefting. The analysis of the coding sequence in the WGS data provides more evidence of the opportunity for novel findings in the African genome.

PARD3 gene variation as candidate cause of nonsyndromic cleft palate only

Nonsyndromic cleft palate only (NSCP) is a common congenital malformation worldwide. In this study, we report a three‐generation pedigree with NSCP following the autosomal‐dominant pattern. Whole‐exome sequencing and Sanger sequencing revealed that only the frameshift variant c.1012dupG [p. E338Gfs*26] in PARD3 cosegregated with the disease. In zebrafish embryos, ethmoid plate patterning defects were observed with PARD3 ortholog disruption or expression of patient‐derived N‐terminal truncating PARD3 (c.1012dupG), which implicated PARD3 in ethmoid plate morphogenesis. PARD3 plays vital roles in determining cellular polarity. Compared with the apical distribution of wild‐type PARD3, PARD3‐p. E338Gfs*26 mainly localized to the basal membrane in 3D‐cultured MCF‐10A epithelial cells. The interaction between PARD3‐p. E338Gfs*26 and endogenous PARD3 was identified by LC–MS/MS and validated by co‐IP. Immunofluorescence analysis showed that PARD3‐p. E338Gfs*26 substantially altered the localization of endogenous PARD3 to the basement membrane in 3D‐cultured MCF‐10A cells. Furthermore, seven variants, including one nonsense variant and six missense variants, were identified in the coding region of PARD3 in sporadic cases with NSCP. Subsequent analysis showed that PARD3‐p. R133*, like the insertion variant of c.1012dupG, also changed the localization of endogenous full‐length PARD3 and that its expression induced abnormal ethmoid plate morphogenesis in zebrafish. Based on these data, we reveal PARD3 gene variation as a novel candidate cause of nonsyndromic cleft palate only.

Revisiting the embryogenesis of lip and palate development

Clefts of the lip and palate (CLP), the major causes of congenital facial malformation globally, result from failure of fusion of the facial processes during embryogenesis. With a prevalence of 1 in 500-2500 live births, CLP causes major morbidity throughout life as a result of problems with facial appearance, feeding, speaking, obstructive apnoea, hearing and social adjustment and requires complex, multi-disciplinary care at considerable cost to healthcare systems worldwide. Long-term outcomes for affected individuals include increased mortality compared with their unaffected siblings. The frequent occurrence and major healthcare burden imposed by CLP highlight the importance of dissecting the molecular mechanisms driving facial development. Identification of the genetic mutations underlying syndromic forms of CLP, where CLP occurs in association with non-cleft clinical features, allied to developmental studies using appropriate animal models is central to our understanding of the molecular events underlying development of the lip and palate and, ultimately, how these are disturbed in CLP.

Decoding the Human Face: Challenges and Progress in Understanding the Genetics of Craniofacial Morphology

Variations in the form of the human face, which plays a role in our individual identities and societal interactions, have fascinated scientists and artists alike. Here, we review our current understanding of the genetics underlying variation in craniofacial morphology and disease-associated dysmorphology, synthesizing decades of progress on Mendelian syndromes in addition to more recent results from genome-wide association studies of human facial shape and disease risk. We also discuss the various approaches used to phenotype and quantify facial shape, which are of particular importance due to the complex, multipartite nature of the craniofacial form. We close by discussing how experimental studies have contributed and will further contribute to our understanding of human genetic variation and then proposing future directions and applications for the field.

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

OBJECTIVE

A significant genetic association between rs7078160 in VAX1 and NSCL/P has been established through genome-wide association studies (GWAS), and we previously replicated the association in the Chinese population. The critical issue in the post-GWAS era is to identify functional variations that have a real impact on disease in the susceptible regions highlighted by GWAS. This study aimed to elucidate functional variants in VAX1 fully.

MATERIALS AND METHODS

Firstly, target sequencing was performed on 159 NSCL/P patients, followed by association analysis to discover disease-associated single nucleotide polymorphisms (SNPs); we then replicated the findings using a larger sample (1626 cases, 2255 controls) and investigated how candidate SNPs affect disease occurrence using extensive annotation databases. Additionally, we compared the genetic profiles of NSCL/P subtypes.

RESULTS

In this study, 6 SNPs in VAX1 were identified to be associated with NSCL/P in the Western Han Chinese population. Five of them were predicted to influence transcriptional factor biding ability and were expression quantitative trait loci (eQTLs) of nearby genes in multiple tissues.

CONCLUSION

The previously reported association between rs7078160 and NSCL/P was successfully replicated. Moreover, our findings firstly revealed that 5 SNPs in VAX1 are associated with NSCL/P in the Western Han Chinese population.

Association of genetic polymorphisms of VAX1, MAFB, and NTN1 with nonsyndromic cleft lip with or without cleft palate in Chinese population

Nonsyndromic cleft lip with or without palate (NSCL/P) is a common birth defect involving genetic factors. We conducted this case-control study to verify the association of ten single-nucleotide polymorphisms (SNPs) of six genes (VAX1, MAFB, PAX7, ABCA4, NTN1, and NOG) with NSCL/P in the Chinese population. The study included 249 NSCL/P patients, 62 nonsyndromic cleft palate only (NSCPO) patients and 480 controls. Three loci, namely, VAX1 rs7078160, MAFB rs11696257, and NTN1 rs4791774, were associated with NSCL/P (Bonferroni method adjusted p values were 0.020, 0.00031, and 0.030, respectively). We also found that the disease risk of individuals carrying both VAX1 rs7078160 and NTN1 rs4791774 was higher than those carrying only one of them (p = 4.50 × 10^-4 and 6.03 × 10^-3, respectively). SNPs of genes VAX1 rs7078160, MAFB rs11696257, and NTN1 rs4791774 increased NSCL/P risk in the Chinese population.

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

OBJECTIVES

Non-syndromic cleft palate only (NSCPO) is one of the most common craniofacial birth defects with largely undetermined genetic etiology. It has been established that Grainyhead-like 3 (GRHL3) plays an essential role in the pathogenesis of NSCPO. This study aimed to identify and verify the first-reported GRHL3 variant underlying NSCPO among the Chinese cohort.

METHODS

We performed whole-exome sequencing (WES) on a Chinese NSCPO patient and identified a rare variant of GRHL3 (p.Arg391His). A validated deleterious variant p.Arg391Cys was introduced as a positive control. Zebrafish embryos injection, reporter assays, live-cell imaging, and RNA sequencing were conducted to test the pathogenicity of the variants.

RESULTS

Zebrafish embryos microinjection demonstrated that overexpression of the variants could disrupt the normal development of zebrafish embryos. Reporter assays showed that Arg391His disturbed transcriptional activity of GRHL3 and exerted a dominant-negative effect. Interestingly, Arg391His and Arg391Cys displayed distinct nuclear localization patterns from that of wild-type GRHL3 in live-cell imaging. Bulk RNA sequencing suggested that the two variants changed the pattern of gene expression.

CONCLUSIONS

In aggregate, this study identified and characterized a rare GRHL3 variant in NSCPO, revealing the critical role of Arginine 391 in GRHL3. Our findings will help facilitate understanding and genetic counseling of NSCPO.

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

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

Genetic and Epigenetic Studies in Nonsyndromic Oral Clefts

The etiology of non-syndromic oral clefts (NSOFC) is complex with genetics, genomics, epigenetics and stochastics factors playing a role. Several approaches have been applied to understand the etiology of non-syndromic oral clefts. These include linkage, candidate gene association studies, genome-wide association studies, whole genome sequencing, copy number variations and epigenetics. In this review we shared these approaches, genes and loci reported in some studies.

Allele-specific transcription factor binding in a cellular model of orofacial clefting

Non-syndromic cleft lip with/without cleft palate (nsCL/P) is a frequent congenital malformation with multifactorial etiology. While recent genome-wide association studies (GWAS) have identified several nsCL/P risk loci, the functional effects of the associated non-coding variants are largely unknown. Furthermore, additional risk loci remain undetected due to lack of power. As genetic variants might alter binding of transcription factors (TF), we here hypothesized that the integration of data from TF binding sites, expression analyses and nsCL/P GWAS might help to (i) identify functionally relevant variants at GWAS loci, and (ii) highlight novel risk variants that have been previously undetected. Analysing the craniofacial TF TFAP2A in human embryonic palatal mesenchyme (HEPM) cells, we identified 2845 TFAP2A ChIP-seq peaks, several of which were located near nsCL/P candidate genes (e.g. MSX1 and SPRY2). Comparison with independent data suggest that 802 of them might be specific to craniofacial development, and genes near these peaks are enriched in processes relevant to nsCL/P. Integration with nsCL/P GWAS data, however, did not show robust evidence for co-localization of common nsCL/P risk variants with TFAP2A ChIP-seq peaks. This data set represents a new resource for the analyses of craniofacial processes, and similar approaches with additional cell lines and TFs could be applied to generate further insights into nsCL/P etiology.

Poverty and Risk of Cleft Lip and Palate: An Analysis of United States Birth Data

BACKGROUND

The relationship between poverty and incidence of cleft lip and cleft palate remains unclear. The authors investigated the association between socioeconomic status and cleft lip with or without cleft palate and cleft palate only in the United States after controlling for demographic and environmental risk factors.

METHODS

The U.S. 2016 and 2017 natality data were utilized. Proxies for socioeconomic status included maternal education, use of the Special Supplemental Nutrition Program for Women, Infants, and Children, and payment source for delivery. Multiple logistic regression controlled for household demographics, prenatal care, maternal health, and infant characteristics.

RESULTS

Of 6,251,308 live births included, 2984 (0.05 percent) had cleft lip with or without cleft palate and 1180 (0.02 percent) had cleft palate only. Maternal education of bachelor's degree or higher was protective against, and delayed prenatal care associated with, cleft lip with or without cleft palate (adjusted ORs = 0.73 and 1.14 to 1.23, respectively; p < 0.02). Receiving assistance under the Special Supplemental Nutrition Program for Women, Infants, and Children was associated with cleft palate only (adjusted OR = 1.25; p = 0.003). Male sex, first-trimester tobacco smoking, and maternal gestational diabetes were also associated with cleft lip with or without cleft palate (adjusted ORs = 1.60, 1.01, and 1.19, respectively; p < 0.05). Female sex, prepregnancy tobacco smoking, and maternal infections during pregnancy were associated with cleft palate only (adjusted ORs = 0.74, 1.02, and 1.60, respectively; p < 0.05).

CONCLUSIONS

Increased incidence of orofacial clefts was associated with indicators of lower socioeconomic status, with different indicators associated with different cleft phenotypes. Notably, early prenatal care was protective against the development of cleft lip with or without cleft palate.

CLIINCAL QUESTION/LEVEL OF EVIDENCE

Risk, III.

A Novel PAX3 Variant in a Chinese Pedigree with Nonsyndromic Cleft Lip With or Without Palate

Objectives: Nonsyndromic cleft lip with or without palate (NSCL/P) is a common congenital orofacial defect, which is associated with severe disruption of orofacial development. The present study was designed to identify potential underlying gene variants in a Chinese pedigree with NSCL/P, in which the proband and the proband's father were affected. Methods: DNA was extracted from the participants' peripheral venous blood, and whole-exome sequencing was performed on the proband and the proband's parents. Results: After filtering, a paired box gene 3 (PAX3) missense variant (c.92C>G_p.Thr31Ser) was identified, which was verified by Sanger sequencing. This variant, which was not present in 113 unrelated healthy individuals or in a Chinese public database, may affect the transcription inhibition domain of the PAX3 protein. Conservation analysis and in silico predictions suggested that this variant may be evolutionarily conserved and potentially deleterious. In addition, it was reported that mice with PAX3 variants show cleft palates. Thus, the PAX3 missense variant (c.92C>G_p.Thr31Ser) is a candidate causative variant in this family. Conclusions: To the best of our knowledge, the present study is the first to report on a PAX3 variant in a pedigree with NSCL/P. The present study further suggests that PAX3 may be associated with CL/P etiology.

Identification of Novel Genomic Variations in Susceptibility to Nonsyndromic Cleft Lip and Palate Patients

BACKGROUND

Nonsyndromic cleft lip with or without palate (NSCL/P) is a multifactorial and common birth malformation caused by genetic and environmental factors, as well as by teratogens. Genome-wide association studies found genetic variations with modulatory effects of NSCL/P formation in Chinese and Iranian populations. We aimed to identify the susceptibility of single-nucleotide polymorphisms (SNPs) to nonsyndromic cleft lip with or without palate in the Indian population.

MATERIAL AND METHODS

The present study was conducted on NSCL/P cases and controls. Genomic DNA was extracted from peripheral blood and Axiom- Precision Medicine Research Array (PMRA) was performed. The Axiom-PMRA covers 902,527 markers and several thousand novel risk variants. Quality control-passed samples were included for candidate genetic variation identification, gene functional enrichment, and pathway and network analysis.

RESULTS

The genome-wide association study identified fourteen novel candidate gene SNPs that showed the most significant association with the risk of NSCL/P, and eight were predicted to have regulatory sequences.

CONCLUSION

The GWAS study showed novel candidate genetic variations in NSCL/P formations. These findings contribute to the understanding of genetic predisposition to nonsyndromic cleft lip with or without palate.

Advances in research regarding the roles of non-coding RNAs in non-syndromic cleft lip with or without cleft palate: A systematic review

OBJECTIVE

To investigate by means of a literature review, what non-coding RNAs engage in non-syndromic cleft lip with or without cleft palate (NSCL/P) and how they lead to the occurrence of this malformation.

DESIGN

A literature search of online databases (Medline via PubMed, Web of Science, Scopus, and Embase) was performed using appropriate keywords (e.g. non-coding RNA, miRNA, lncRNA, NSCL/P, non-syndromic cleft lip only, and non-syndromic orofacial cleft). The risk of bias in the included studies was then assessed, and a comprehensive review of reported non-coding RNAs associated with NSCL/P was performed.

RESULTS

The initial search retrieved 133 studies reporting non-coding RNAs associated with NSCL/P; after excluding 18 replicates and 77 ineligible studies, 35 remained. Of these, 16 studies fulfilled all the criteria and were included in the systematic review. These studies established the roles of non-coding RNAs in the development of craniofacial structures. The differential expression of these non-coding RNAs could lead to orofacial clefts, indicating their significance in NSCL/P and their profound research value.

CONCLUSION

There is evidence that non-coding RNAs are involved in the formation of NSCL/P. Specifically, they play significant roles in the regulation of genes and signalling pathways related to NSCL/P.

Genetic architecture of orbital telorism

The interocular distance, or orbital telorism, is a distinctive craniofacial trait that also serves as a clinically informative measure. While its extremes, hypo- and hypertelorism, have been linked to monogenic disorders and are often syndromic, little is known about the genetic determinants of interocular distance within the general population. We derived orbital telorism measures from cranial magnetic resonance imaging by calculating the distance between the eyeballs’ centre of gravity, which showed a good reproducibility with an intraclass correlation coefficient of 0.991 (95% confidence interval 0.985–0.994). Heritability estimates were 76% (standard error = 12%) with a family-based method (N = 364) and 39% (standard error = 2.4%) with a single nucleotide polymorphism-based method (N = 34 130) and were unaffected by adjustment for height (model II) and intracranial volume (model III) or head width (model IV). Genome-wide association studies in 34 130 European individuals identified 56 significantly associated genomic loci (P < 5 × 10⁻⁸) across four different models of which 46 were novel for facial morphology, and overall these findings replicated in an independent sample (N = 10 115) with telorism-related horizontal facial distance measures. Genes located nearby these 56 identified genetic loci were 4.9-fold enriched for Mendelian hypotelorism and hypertelorism genes, underlining their biological relevance. This study provides novel insights into the genetic architecture underlying interocular distance in particular, and the face in general, and explores its potential for applications in a clinical setting.

Association Between SNPs in 1q32.2 and NSCL ± P in Han Chinese Population

OBJECTIVES

Non-syndromic cleft lip with or without cleft palate (NSCL ± P) is one of the most common birth malformations. Currently, numerous susceptibility SNPs have been reported by GWA studies, however, the replications of them among NSCL ± P from Han Chinese were very limited.

DESIGN

In this study, we selected 16 SNPs around 1q32.2 based on the published GWA studies and replicated them among 302 trios with NSCL ± P from Han Chinese Population. The genotypic data was analyzed with FBAT, PLINK and R package.

SETTING

The study was conducted in a tertiary medical center.

PATIENTS, PARTICIPANTS

302 patients with CL ± P and their parents.

MAIN OUTCOME MEASURES

To ascertain the genetic variants in 1q32.2 in patients with CL ± P in Han Chinese Population.

INTERVENTIONS

Blood samples were collected.

RESULTS

We found T allele (Z = 4.26, p = 0.00002) and T/T homozygotes (Z = 4.4, p = 0.000011) at rs12063989 was significantly over-transmitted among non-syndromic cleft lip with or without cleft palate (NSCL ± P).

CONCLUSIONS

We found rs12063989 exhibited significant association with the occurrence of NSCL ± P, which would provide new evidence for the future study in the etiology of NSCL ± P.

Genome-wide Gene-by-Sex Interaction Studies Identify Novel Nonsyndromic Orofacial Clefts Risk Locus

Risk loci identified through genome-wide association studies have explained about 25% of the phenotypic variations in nonsyndromic orofacial clefts (nsOFCs) on the liability scale. Despite the notable sex differences in the incidences of the different cleft types, investigation of loci for sex-specific effects has been understudied. To explore the sex-specific effects in genetic etiology of nsOFCs, we conducted a genome-wide gene × sex (GxSex) interaction study in a sub-Saharan African orofacial cleft cohort. The sample included 1,019 nonsyndromic orofacial cleft cases (814 cleft lip with or without cleft palate and 205 cleft palate only) and 2,159 controls recruited from 3 sites (Ethiopia, Ghana, and Nigeria). An additive logistic model was used to examine the joint effects of the genotype and GxSex interaction. Furthermore, we examined loci with suggestive significance (P < 1E-5) in the additive model for the effect of the GxSex interaction only. We identified a novel risk locus on chromosome 8p22 with genome-wide significant joint and GxSex interaction effects (rs2720555, p_2df = 1.16E-08, p_GxSex = 1.49E-09, odds ratio [OR] = 0.44, 95% CI = 0.34 to 0.57). For males, the risk of cleft lip with or without cleft palate at this locus decreases with additional copies of the minor allele (p < 0.0001, OR = 0.60, 95% CI = 0.48 to 0.74), but the effect is reversed for females (p = 0.0004, OR = 1.36, 95% CI = 1.15 to 1.60). We replicated the female-specific effect of this locus in an independent cohort (p = 0.037, OR = 1.30, 95% CI = 1.02 to 1.65), but no significant effect was found for the males (p = 0.29, OR = 0.86, 95% CI = 0.65 to 1.14). This locus is in topologically associating domain with craniofacially expressed and enriched genes during embryonic development. Rare coding mutations of some of these genes were identified in nsOFC cohorts through whole exome sequencing analysis. Our study is additional proof that genome-wide GxSex interaction analysis provides an opportunity for novel findings of loci and genes that contribute to the risk of nsOFCs.

Mycn deficiency underlies the development of orofacial clefts in mice and humans

Non-syndromic cleft lip with or without cleft palate (NSCL/P) is the most common subphenotype of non-syndromic orofacial clefts (NSOFCs) arising from genetic and/or environmental perturbations during embryonic development. We previously identified 2p24.2 as a risk locus associated with NSCL/P in the Chinese Han population, and MYCN is a candidate risk gene in this region. To understand the potential function of MYCN in craniofacial development, we generated Wnt1-Cre;Mycnflox/flox mice that exhibited cleft palate, microglossia, and micrognathia, resembling the Pierre Robin sequence (PRS) in humans. Further analyses indicated that the cleft palate was secondary to the delayed elevation of palatal shelves caused by micrognathia. The micrognathia resulted from impaired chondrogenic differentiation in Merkel's cartilage, which limited tongue development, leading to microglossia. In terms of mechanism, Mycn deficiency in cranial neural crest cells (CNCCs) downregulated Sox9 expression by inhibiting Wnt5a in a CNCC-derived chondrogenic lineage in Merkel's cartilage. To investigate whether MYCN deficiency contributed to NSCL/P, we performed direct sequencing targeting all exons and exon-intron boundaries of MYCN in 104 multiplex families with mendelian NSCL/P and identified a novel pathogenic variant in MYCN. Taken together, our data indicate that ablation of Mycn in mouse CNCCs could resemble PRS by suppressing the Wnt5a-Sox9 signaling pathway in Merkel's cartilage and that mutations in MYCN may be novel potential causes of NSCL/P.

Integrative approaches generate insights into the architecture of non-syndromic cleft lip with or without cleft palate

Non-syndromic cleft lip with or without cleft palate (nsCL/P) is a common congenital facial malformation with a multifactorial etiology. Genome-wide association studies (GWASs) have identified multiple genetic risk loci. However, functional interpretation of these loci is hampered by the underrepresentation in public resources of systematic functional maps representative of human embryonic facial development. To generate novel insights into the etiology of nsCL/P, we leveraged published GWAS data on nsCL/P as well as available chromatin modification and expression data on mid-facial development. Our analyses identified five novel risk loci, prioritized candidate target genes within associated regions, and highlighted distinct pathways. Furthermore, the results suggest the presence of distinct regulatory effects of nsCL/P risk variants throughout mid-facial development and shed light on its regulatory architecture. Our integrated data provide a platform to advance hypothesis-driven molecular investigations of nsCL/P and other human facial defects.

Detecting Gene-Environment Interaction for Maternal Exposures Using Case-Parent Trios Ascertained Through a Case With Non-Syndromic Orofacial Cleft

Two large studies of case-parent trios ascertained through a proband with a non-syndromic orofacial cleft (OFC, which includes cleft lip and palate, cleft lip alone, or cleft palate alone) were used to test for possible gene-environment (G × E) interaction between genome-wide markers (both observed and imputed) and self-reported maternal exposure to smoking, alcohol consumption, and multivitamin supplementation during pregnancy. The parent studies were as follows: GENEVA, which included 1,939 case-parent trios recruited largely through treatment centers in Europe, the United States, and Asia, and 1,443 case-parent trios from the Pittsburgh Orofacial Cleft Study (POFC) also ascertained through a proband with an OFC including three major racial/ethnic groups (European, Asian, and Latin American). Exposure rates to these environmental risk factors (maternal smoking, alcohol consumption, and multivitamin supplementation) varied across studies and among racial/ethnic groups, creating substantial differences in power to detect G × E interaction, but the trio design should minimize spurious results due to population stratification. The GENEVA and POFC studies were analyzed separately, and a meta-analysis was conducted across both studies to test for G × E interaction using the 2 df test of gene and G × E interaction and the 1 df test for G × E interaction alone. The 2 df test confirmed effects for several recognized risk genes, suggesting modest G × E effects. This analysis did reveal suggestive evidence for G × Vitamin interaction for CASP9 on 1p36 located about 3 Mb from PAX7, a recognized risk gene. Several regions gave suggestive evidence of G × E interaction in the 1 df test. For example, for G × Smoking interaction, the 1 df test suggested markers in MUSK on 9q31.3 from meta-analysis. Markers near SLCO3A1 also showed suggestive evidence in the 1 df test for G × Alcohol interaction, and rs41117 near RETREG1 (a.k.a. FAM134B) also gave suggestive significance in the meta-analysis of the 1 df test for G × Vitamin interaction. While it remains quite difficult to obtain definitive evidence for G × E interaction in genome-wide studies, perhaps due to small effect sizes of individual genes combined with low exposure rates, this analysis of two large case-parent trio studies argues for considering possible G × E interaction in any comprehensive study of complex and heterogeneous disorders such as OFC.

Identification of New miRNA-mRNA Networks in the Development of Non-syndromic Cleft Lip With or Without Cleft Palate

Objective: To identify new microRNA (miRNA)-mRNA networks in non-syndromic cleft lip with or without cleft palate (NSCL/P).

Materials and Methods: Overlapping differentially expressed miRNAs (DEMs) were selected from cleft palate patients (GSE47939) and murine embryonic orofacial tissues (GSE20880). Next, the target genes of DEMs were predicted by Targetscan, miRDB, and FUNRICH, and further filtered through differentially expressed genes (DEGs) from NSCL/P patients and controls (GSE42589), MGI, MalaCards, and DECIPHER databases. The results were then confirmed by in vitro experiments. NSCL/P lip tissues were obtained to explore the expression of miRNAs and their target genes.

Results: Let-7c-5p and miR-193a-3p were identified as DEMs, and their overexpression inhibited cell proliferation and promoted cell apoptosis. PIGA and TGFB2 were confirmed as targets of let-7c-5p and miR-193a-3p, respectively, and were involved in craniofacial development in mice. Negative correlation between miRNA and mRNA expression was detected in the NSCL/P lip tissues. They were also associated with the occurrence of NSCL/P based on the MGI, MalaCards, and DECIPHER databases.

Conclusions: Let-7c-5p-PIGA and miR-193a-3p-TGFB2 networks may be involved in the development of NSCL/P.

Orofacial clefts embryology, classification, epidemiology, and genetics

Orofacial clefts (OFCs) rank as the second most common congenital birth defect in the United States after Down syndrome and are the most common head and neck congenital malformations. They are classified as cleft lip with or without cleft palate (CL/P) and cleft palate only (CPO). OFCs have significant psychological and socio-economic impact on patients and their families and require a multidisciplinary approach for management and counseling. A complex interaction between genetic and environmental factors contributes to the incidence and clinical presentation of OFCs. In this comprehensive review, the embryology, classification, epidemiology and etiology of clefts are thoroughly discussed and a "state-of-the-art" snapshot of the recent advances in the genetics of OFCs is presented.

Integrating GWAS and eQTL to predict genes and pathways for non-syndromic cleft lip with or without palate

OBJECTIVE

To explore susceptibility genes and pathways for non-syndromic cleft lip with or without cleft palate (NSCL/P).

MATERIALS AND METHODS

Two genome-wide association studies (GWAS) datasets, including 858 NSCL/P cases and 1,248 controls, were integrated with expression quantitative trait loci (eQTL) dataset identified by Genotype-Tissue Expression (GTEx) project in whole-blood samples. The expression of the candidate genes in mouse orofacial development was inquired from FaceBase. Protein-protein interaction (PPI) network was visualized to identify protein functions. Go and KEGG pathway analyses were performed to explore the underlying risk pathways.

RESULTS

A total of 233 eQTL single-nucleotide polymorphisms (SNPs) in 432 candidate genes were identified to be associated with the risk of NSCL/P. One hundred and eighty-three susceptible genes were expressed in mouse orofacial development according to FaceBase. PPI network analysis highlighted that these genes involved in ubiquitin-mediated proteolysis (KCTD7, ASB1, UBOX5, ANAPC4) and DNA synthesis (XRCC3, RFC3, KAT5, RHNO1) were associated with the risk of NSCL/P. GO and KEGG pathway analyses revealed that the fatty acid metabolism pathway (ACADL, HSD17B12, ACSL5, PPT1, MCAT) played an important role in the development of NSCL/P.

CONCLUSIONS

Our results identified novel susceptibility genes and pathways associated with the development of NSCL/P.